搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

液-液两相液层间传质过程的Rayleigh-Bénard-Marangoni对流特性

陈俊 沈超群 王贺 张程宾

液-液两相液层间传质过程的Rayleigh-Bénard-Marangoni对流特性

陈俊, 沈超群, 王贺, 张程宾
PDF
HTML
导出引用
导出核心图
  • 传质引发的Rayleigh-Bénard-Marangoni对流(RBM对流)对化工传递过程有着显著影响. 但是, 已有的相关研究多集中于气-液体系, 并且有限的针对液-液体系的相关研究尚缺乏对RBM对流演化及其引发的界面扰动行为的深入分析. 因此, 本文基于阴影法设计搭建了竖直狭缝内液-液两相液层间传质过程的RBM对流特性可视化实验平台, 并实验观测了水-甲苯-丙酮三元体系中丙酮组分扩散传质时出现的RBM对流结构以及其向下层水相主体的发展演变过程, 探讨了水相丙酮初始浓度、甲苯相丙酮初始浓度以及甲苯层厚度对RBM对流特性和液-液界面形貌的影响. 研究表明: 在Rayleigh-Taylor不稳定性作用下, 水相上层密度(重力)分层“界面”下凸沉降形成波浪形丘状“界面”, 并随着“界面”处密度与压力失调的加剧而演变成羽状流; 因羽流区“界面”不同浓度梯度引起的传质特性差异, 羽状流又可以演变成弱羽状流和强羽状流两种形态; 当丙酮浓度梯度增大到一定程度后, 近界面处短时间内产生大量RBM对流结构, 且结构间相互影响增强而聚并成对流团, 并随着传质过程的进行, 逐渐演变成独立的强羽状流; RBM对流强度与上下液层丙酮浓度梯度大小呈正相关关系, 且液-液界面粗糙度及其非稳态波动随着丙酮浓度梯度的增加而增大.
      通信作者: 张程宾, cbzhang@seu.edu.cn
    • 基金项目: 国家自然科学基金委员会-中国工程物理研究院联合基金(批准号: U1530260)、国家自然科学基金(批准号: 51706193)和江苏省高校自然科学研究项目(批准号: 17KJB470014)资助的课题.
    [1]

    Yao F, Chen Y P, Peterson G P 2013 Int. J. Heat Mass Transfer 64 418

    [2]

    Liu X D, Chen Y P, Shi M H 2013 Int. J. Therm. Sci. 65 224

    [3]

    Chen Y P, Cheng P 2005 Int. J. Heat Mass Transfer 32 931

    [4]

    Bodenschatz E, Pesch W, Ahlers G 2000 Annu. Rev. Fluid Mech. 32 709

    [5]

    王飞, 彭岚, 张全壮, 刘佳 2015 物理学报 64 140202

    Wang F, Peng L, Zhang Q Z, Liu J 2015 Acta Phys. Sin. 64 140202

    [6]

    翟薇, 王楠, 魏炳波 2007 物理学报 56 2353

    Zhai W, Wang N, Wei B B 2007 Acta Phys. Sin. 56 2353

    [7]

    Schwabe D 1999 Adv. Space Res. 24 1347

    [8]

    Touazi O, Chénier E, Doumenc F, Guerrier B 2010 Int. J. Heat Mass Transfer 53 656

    [9]

    Chen J, Yang C, Mao Z S 2015 Eur. Phys. J. Spec. Top. 224 389

    [10]

    Bo Z, Mao S, Han Z J, Cen K F, Chen J H, Kostya O 2015 Chem. Soc. Rev. 44 2018

    [11]

    张婷, 施保昌, 柴振华 2015 物理学报 64 254701

    Zhang T, Shi B C, Chai Z H 2015 Acta Phys. Sin. 64 254701

    [12]

    郑连存, 盛晓艳, 张欣欣 2006 物理学报 55 5298

    Zheng L C, Sheng X Y, Zhang X X 2006 Acta Phys. Sin. 55 5298

    [13]

    Kline J L, Hager J D 2016 Matter Radiat. Extremes 2 16

    [14]

    Dong S X, Han W, Liu M F, Zhang Z W, Li B, Ge L Q 2016 Colloids Surf. A 509 32

    [15]

    Bai L, Zhao S F, Fu Y H, Cheng Y 2016 Biochem. Eng. J. 298 281

    [16]

    Sun Z F 2012 Chem. Eng. Sci. 68 579

    [17]

    Liu C, Zeng A, Yuan X, Yu G 2008 Chem. Eng. Res. Des. 86 201

    [18]

    Alvarez-Herrera C, Moreno-Hernández D, Barrientos-García B, Guerrero-Viramontes J A 2005 Opt. Laser Technol. 41 233

    [19]

    Piekarska W, Kubiak M 2013 Appl. Math. Modell. 37 2051

    [20]

    Szymczyk J A 1991 Can. J. Chem. Eng. 69 1233

    [21]

    Okhotsimskii A, Hozawa M 1998 Chem. Eng. Sci. 53 2547

    [22]

    王勇, 张泽廷 2002 北京化工大学学报 29 11

    Wang Y, Zhang Z T 2002 J. Beijing Univ. Chem. Technol. 29 11

    [23]

    Sun Z F, Yu K T, S Y W, Miao Y Z 2002 Ind. Eng. Chem. Res. 41 1905

    [24]

    沙勇, 李樟云, 林芬芬, 吐芬, 肖宗源, 叶李艺 2010 化工学报 61 844

    Sha Y, Li Z Y, Lin F F, Tu F, Xiao Z Y, Ye L Y 2010 J. Chem. Ind. Eng. 61 844

    [25]

    Orell A, Westwater J W 1961 AlChE J. 8 350

    [26]

    Zhang S H, Wang Z M, Su Y F 1990 Chem. Eng. Res. Des. 68 84

    [27]

    Guzun-Stoica A, Kurzeluk M, Floarea O 2000 Chem. Eng. Sci. 55 3813

    [28]

    Kostarev K G, Shmyrov A V, Zuev A L, Viviani A 2011 Exp. Fluids 51 457

    [29]

    Chen Y, Cheng P 2005 Int. Commun. Heat Mass Transfer 32 175

    [30]

    Agble D, Mendes-Tatsis M A 2000 Int. J. Heat Mass Transfer 43 1025

    [31]

    Shi Y, Kerstin E 2007 Chin. J. Chem. Eng. 15 748

    [32]

    Chen Y P, Liu X D, Shi M H 2013 Appl. Phys. Lett. 102 051609

    [33]

    Chen Y P, Liu X D, Zhao Y J 2015 Appl. Phys. Lett. 106 141601

    [34]

    Chen Y P, Wu L Y, Zhang L 2015 Int. J. Heat Mass Transfer 82 42

    [35]

    Sharp D H 1984 Physica D 12 3

    [36]

    Roberts M S, Jacobs J W 2015 J. Fluid Mech. 787 50

    [37]

    胡楠, 张会书, 傅强, 李陆星, 袁希钢 2016 化工学报 68 584

    Hu L, Zhang H S, Fu Q, Li L X, Yuan X G 2016 J. Chem. Ind. Eng. 68 584

    [38]

    Puthenveettil B A, Arakeri J H 2005 J. Fluid Mech. 542 217

    [39]

    Yang C, Tartaglino U, Persson B N 2006 J. Phys. Rev. Lett. 97 11

  • 图 1  RBM对流的实验系统图 (a) 阴影法实验系统图; (b) 示踪粒子法实验系统图; (c) 玻璃狭缝尺寸图

    Fig. 1.  Schematic diagram of the experimental system for RBM convection: (a) Schematic diagram of the experimental system based on shadowgraph method; (b) schematic diagram of the experimental system based on particle tracer method; (c) size of the glass slit.

    图 2  T = 20 ℃时不同水相丙酮浓度下水-甲苯两相间的界面张力系数

    Fig. 2.  Interfacial tension coefficient between water and toluene phases under different acetone concentrations of aqueous phase with T = 20 ℃.

    图 3  T = 20 ℃时不同丙酮浓度的水相溶液的密度

    Fig. 3.  Density of aqueous solution with different acetone concentrations with T = 20 ℃.

    图 4  传质过程引起的密度分层示意图与实验结果 (a) 密度分层示意图; (b) 实验图像

    Fig. 4.  Schematic diagram and experimental result of density stratification caused by mass transfer: (a) Schematic diagram of density stratification; (b) experimental image.

    图 6  丘状“界面”的形成过程(水相丙酮初始体积浓度${\varphi _0} = 5\% $, 甲苯相丙酮初始体积浓度${\varphi _1} = 0\% $) (a) t = 13 s; (b) t = 30 s; (c) t = 36 s; (d) t = 42 s

    Fig. 6.  The forming process of the mound “interface” (the initial volume concentration of acetone in aqueous phase ${\varphi _0} = 5\% $, the initial volume concentration of acetone in the toluene phase ${\varphi _1} = 0\% $): (a) t = 13 s; (b) t = 30 s; (c) t = 36 s; (d) t = 42 s.

    图 5  密度分层引起的Rayleigh-Taylor不稳定性示意图 $\omega $是涡流, P是压力, $\rho $是密度, u是速度, g是重力加速度; 粗的环形箭头表示涡旋产生的速度场

    Fig. 5.  Schematic diagram of Rayleigh-Taylor instability caused by density stratification. $\omega $ is vorticity, P is pressure, $\rho $ is density, u is velocity and g is acceleration of gravity; the thick circular arrows represent the velocity field created by the vortex.

    图 7  羽状流的演变过程 (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 0 s; (b) t = 18 s; (c) t = 20 s; (d) t = 22 s

    Fig. 7.  The evolution of the plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 0 s; (b) t = 18 s; (c) t = 20 s; (d) t = 22 s.

    图 8  弱羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 348 s; (b) t = 356 s; (c) t = 363 s; (d) t = 370 s

    Fig. 8.  The evolution of the weak plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 348 s; (b) t = 356 s; (c) t = 363 s; (d) t = 370 s.

    图 10  强羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 22 s; (b) t = 23 s; (c) t = 26 s

    Fig. 10.  The evolution of the strong plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 22 s; (b) t = 23 s; (c) t = 26 s.

    图 9  羽状流的速度矢量及涡量云图(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $)

    Fig. 9.  The velocity vector and vorticity contours of the plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $).

    图 11  强羽状流向弱羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 228 s; (b) t = 234 s; (c) t = 238 s; (d) t = 245 s

    Fig. 11.  The evolution of the strong plume flow to the weak plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 228 s; (b) t = 234 s; (c) t = 238 s; (d) t = 245 s.

    图 12  传质初期对流结构的聚并过程(${\varphi _0} \!=\! 30\% $, ${\varphi _1} \!=\! 0\% $)  (a) t = 15 s; (b) t = 20 s; (c) t = 24 s; (d) t = 30 s

    Fig. 12.  Convergence process of convective structure at the beginning of the mass transfer (${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $): (a) t = 15 s; (b) t = 20 s; (c) t = 24 s; (d) t = 30 s.

    图 13  对流团的消失以及强羽状流的出现(${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $)  (a) t = 101 s; (b) t = 196 s; (c) t = 259 s; (d) t = 350 s

    Fig. 13.  The disappearance of convective cloud and the appearance of the strong plume (${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $): (a) t = 101 s; (b) t = 196 s; (c) t = 259 s; (d) t = 350 s.

    图 14  水相丙酮初始浓度对第一个RBM对流结构向下延伸速度的影响

    Fig. 14.  The influence of initial concentration of acetone in aqueous phase on the elongation velocity of the first RBM convective structure.

    图 15  水相丙酮初始浓度对羽状流数量的影响

    Fig. 15.  The influence of initial concentration of acetone in aqueous phase on the number of the plumes.

    图 16  水相丙酮初始浓度对水-甲苯界面形貌的影响 (a)界面粗糙度; (b)界面波动程度

    Fig. 16.  The influence of initial concentration of acetone in aqueous phase on water-toluene interface morphology: (a) Interfacial roughness; (b) the degree of interface fluctuation.

    图 17  不同甲苯相丙酮初始浓度下的投影图像(t = 35 s, ${\varphi _0} = 20\% $) (a) ${\varphi _1} = 5\% $; (b) ${\varphi _1} = 7.5\% $; (c) ${\varphi _1} = 10\% $; (d) ${\varphi _1} = 15\% $

    Fig. 17.  Schlieren images under different initial concentrations of acetone in the toluene phase (t = 35 s, ${\varphi _0} = 20\% $): (a) ${\varphi _1} = 5\% $ (b) ${\varphi _1} = 7.5\% $; (c) ${\varphi _1} = 10\% $; (d) ${\varphi _1} = 15\% $.

    图 18  甲苯相丙酮初始浓度对第一个RBM对流结构向下延伸速度的影响

    Fig. 18.  The influence of initial concentration of acetone in the toluene phase on the elongation velocity of the first RBM convective structure.

    图 19  甲苯相丙酮初始浓度对羽状流数量的影响

    Fig. 19.  The influence of initial concentration of acetone in the toluene phase on the number of the plumes.

    图 20  甲苯相丙酮初始浓度对水-甲苯界面形貌的影响 (a)界面粗糙度; (b)界面波动程度

    Fig. 20.  The influence of initial concentration of acetone in the toluene phase on water-toluene interface morphology: (a) Interfacial roughness; (b) the degree of interface fluctuation.

    图 21  甲苯层厚度对羽状流数量的影响

    Fig. 21.  The influence of thickness of toluene layer on the number of the plumes.

    图 22  甲苯层厚度对水-甲苯界面形貌的影响

    Fig. 22.  The influence of thickness of toluene layer on water - toluene interface morphology.

    表 1  实验试剂的物性参数(T = 20 ℃, P = 0.1 MPa)

    Table 1.  The physical parameters of experimental reagents (T = 20 ℃, P = 0.1 MPa).

    实验试剂ρ/kg·m–3µ/10–4Pa·s
    99810.04
    甲苯8675.86
    丙酮7903.26
    下载: 导出CSV
  • [1]

    Yao F, Chen Y P, Peterson G P 2013 Int. J. Heat Mass Transfer 64 418

    [2]

    Liu X D, Chen Y P, Shi M H 2013 Int. J. Therm. Sci. 65 224

    [3]

    Chen Y P, Cheng P 2005 Int. J. Heat Mass Transfer 32 931

    [4]

    Bodenschatz E, Pesch W, Ahlers G 2000 Annu. Rev. Fluid Mech. 32 709

    [5]

    王飞, 彭岚, 张全壮, 刘佳 2015 物理学报 64 140202

    Wang F, Peng L, Zhang Q Z, Liu J 2015 Acta Phys. Sin. 64 140202

    [6]

    翟薇, 王楠, 魏炳波 2007 物理学报 56 2353

    Zhai W, Wang N, Wei B B 2007 Acta Phys. Sin. 56 2353

    [7]

    Schwabe D 1999 Adv. Space Res. 24 1347

    [8]

    Touazi O, Chénier E, Doumenc F, Guerrier B 2010 Int. J. Heat Mass Transfer 53 656

    [9]

    Chen J, Yang C, Mao Z S 2015 Eur. Phys. J. Spec. Top. 224 389

    [10]

    Bo Z, Mao S, Han Z J, Cen K F, Chen J H, Kostya O 2015 Chem. Soc. Rev. 44 2018

    [11]

    张婷, 施保昌, 柴振华 2015 物理学报 64 254701

    Zhang T, Shi B C, Chai Z H 2015 Acta Phys. Sin. 64 254701

    [12]

    郑连存, 盛晓艳, 张欣欣 2006 物理学报 55 5298

    Zheng L C, Sheng X Y, Zhang X X 2006 Acta Phys. Sin. 55 5298

    [13]

    Kline J L, Hager J D 2016 Matter Radiat. Extremes 2 16

    [14]

    Dong S X, Han W, Liu M F, Zhang Z W, Li B, Ge L Q 2016 Colloids Surf. A 509 32

    [15]

    Bai L, Zhao S F, Fu Y H, Cheng Y 2016 Biochem. Eng. J. 298 281

    [16]

    Sun Z F 2012 Chem. Eng. Sci. 68 579

    [17]

    Liu C, Zeng A, Yuan X, Yu G 2008 Chem. Eng. Res. Des. 86 201

    [18]

    Alvarez-Herrera C, Moreno-Hernández D, Barrientos-García B, Guerrero-Viramontes J A 2005 Opt. Laser Technol. 41 233

    [19]

    Piekarska W, Kubiak M 2013 Appl. Math. Modell. 37 2051

    [20]

    Szymczyk J A 1991 Can. J. Chem. Eng. 69 1233

    [21]

    Okhotsimskii A, Hozawa M 1998 Chem. Eng. Sci. 53 2547

    [22]

    王勇, 张泽廷 2002 北京化工大学学报 29 11

    Wang Y, Zhang Z T 2002 J. Beijing Univ. Chem. Technol. 29 11

    [23]

    Sun Z F, Yu K T, S Y W, Miao Y Z 2002 Ind. Eng. Chem. Res. 41 1905

    [24]

    沙勇, 李樟云, 林芬芬, 吐芬, 肖宗源, 叶李艺 2010 化工学报 61 844

    Sha Y, Li Z Y, Lin F F, Tu F, Xiao Z Y, Ye L Y 2010 J. Chem. Ind. Eng. 61 844

    [25]

    Orell A, Westwater J W 1961 AlChE J. 8 350

    [26]

    Zhang S H, Wang Z M, Su Y F 1990 Chem. Eng. Res. Des. 68 84

    [27]

    Guzun-Stoica A, Kurzeluk M, Floarea O 2000 Chem. Eng. Sci. 55 3813

    [28]

    Kostarev K G, Shmyrov A V, Zuev A L, Viviani A 2011 Exp. Fluids 51 457

    [29]

    Chen Y, Cheng P 2005 Int. Commun. Heat Mass Transfer 32 175

    [30]

    Agble D, Mendes-Tatsis M A 2000 Int. J. Heat Mass Transfer 43 1025

    [31]

    Shi Y, Kerstin E 2007 Chin. J. Chem. Eng. 15 748

    [32]

    Chen Y P, Liu X D, Shi M H 2013 Appl. Phys. Lett. 102 051609

    [33]

    Chen Y P, Liu X D, Zhao Y J 2015 Appl. Phys. Lett. 106 141601

    [34]

    Chen Y P, Wu L Y, Zhang L 2015 Int. J. Heat Mass Transfer 82 42

    [35]

    Sharp D H 1984 Physica D 12 3

    [36]

    Roberts M S, Jacobs J W 2015 J. Fluid Mech. 787 50

    [37]

    胡楠, 张会书, 傅强, 李陆星, 袁希钢 2016 化工学报 68 584

    Hu L, Zhang H S, Fu Q, Li L X, Yuan X G 2016 J. Chem. Ind. Eng. 68 584

    [38]

    Puthenveettil B A, Arakeri J H 2005 J. Fluid Mech. 542 217

    [39]

    Yang C, Tartaglino U, Persson B N 2006 J. Phys. Rev. Lett. 97 11

  • [1] 陈书赢, 王海斗, 徐滨士, 康嘉杰. 基于分形理论的超音速等离子喷涂层界面结合行为研究. 物理学报, 2014, 63(15): 156801. doi: 10.7498/aps.63.156801
    [2] 蔡继兴, 郭明, 渠旭, 李贺, 金光勇. 激光诱导等离子体的气体动力学和燃烧波扩展速度研究. 物理学报, 2017, 66(9): 094202. doi: 10.7498/aps.66.094202
    [3] 张敏梁, 田煜, 蒋继乐, 孟永钢, 温诗铸. 极板形貌修饰对电流变液/极板界面滑移抑制实验研究. 物理学报, 2009, 58(12): 8394-8399. doi: 10.7498/aps.58.8394
    [4] 王友琴, 陈式刚. Rayleigh-Bénard对流中的波长增长现象与最大熵产生判据. 物理学报, 1983, 32(2): 209-215. doi: 10.7498/aps.32.209
    [5] 检测超声组. 液面法超声全息成象技术. 物理学报, 1975, 141(1): 12-20. doi: 10.7498/aps.24.12
    [6] 王飞, 彭岚, 张全壮, 刘佳. 水平温差对环形浅液池内Marangoni-热毛细对流的影响. 物理学报, 2015, 64(14): 140202. doi: 10.7498/aps.64.140202
    [7] 郑小青, 杨洋, 孙得彦. 模型二元有序合金固液界面结构的分子动力学研究. 物理学报, 2013, 62(1): 017101. doi: 10.7498/aps.62.017101
    [8] 小松启, 金蔚青. 固液界面温度的一种测量方法. 物理学报, 1985, 34(9): 1166-1172. doi: 10.7498/aps.34.1166
    [9] 周化光, 林鑫, 王猛, 黄卫东. Cu固液界面能的分子动力学计算. 物理学报, 2013, 62(5): 056803. doi: 10.7498/aps.62.056803
    [10] 尹灵康, 徐顺, Seongmin Jeong, Yongseok Jho, 王健君, 周昕. 广义等温等压系综-分子动力学模拟全原子水的气液共存形貌. 物理学报, 2017, 66(13): 136102. doi: 10.7498/aps.66.136102
    [11] 李志宏, 丁召, 汤佳伟, 王一, 罗子江, 马明明, 黄延彬, 张振东, 郭祥. Ga液滴沉积速率对GaAs/GaAs (001)量子双环形貌的影响. 物理学报, 2019, 68(18): 183601. doi: 10.7498/aps.68.20190615
    [12] 涂相征. 稳定自然对流下的温度梯度液相外延. 物理学报, 1982, 31(1): 78-89. doi: 10.7498/aps.31.78
    [13] 潘宵, 鞠焕鑫, 冯雪飞, 范其瑭, 王嘉兴, 杨耀文, 朱俊发. F8BT薄膜表面形貌及与Al形成界面的电子结构和反应. 物理学报, 2015, 64(7): 077304. doi: 10.7498/aps.64.077304
    [14] 叶贞成, 蔡 钧, 张书令, 刘洪来, 胡 英. 方阱链流体在固液界面分布的密度泛函理论研究. 物理学报, 2005, 54(9): 4044-4052. doi: 10.7498/aps.54.4044
    [15] 王晓冬, 董 鹏, 陈胜利, 仪桂云. 亚微米聚苯乙烯微球在气-液界面组装的机理研究. 物理学报, 2007, 56(3): 1831-1836. doi: 10.7498/aps.56.1831
    [16] 王晓冬, 董 鹏, 陈胜利, 仪桂云. 亚微米聚苯乙烯微球在气-液界面组装的机理研究. 物理学报, 2007, 56(5): 3017-3021. doi: 10.7498/aps.56.3017
    [17] 张碧星, 王文龙. 凹面相控阵聚焦声场在液固界面上的反射和折射. 物理学报, 2008, 57(6): 3613-3619. doi: 10.7498/aps.57.3613
    [18] 张旋, 张天赐, 葛际江, 蒋平, 张贵才. 表面活性剂对气-液界面纳米颗粒吸附规律的影响. 物理学报, 2020, 69(2): 026801. doi: 10.7498/aps.69.20190756
    [19] 周广刚, 卢贵武, 矫玉秋, 李英峰, 王坤, 于养信. KDP晶体固-液界面吸附行为的分子模拟研究. 物理学报, 2012, 61(1): 010204. doi: 10.7498/aps.61.010204
    [20] 葛宋, 陈民. 接触角与液固界面热阻关系的分子动力学模拟. 物理学报, 2013, 62(11): 110204. doi: 10.7498/aps.62.110204
  • 引用本文:
    Citation:
计量
  • 文章访问数:  184
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-07-05
  • 修回日期:  2019-01-19
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-04-01

液-液两相液层间传质过程的Rayleigh-Bénard-Marangoni对流特性

  • 1. 东南大学能源与环境学院, 能源热转换及其过程测控教育部重点实验室, 南京 210096
  • 2. 扬州大学水利与能源动力工程学院, 扬州 225127
  • 通信作者: 张程宾, cbzhang@seu.edu.cn
    基金项目: 国家自然科学基金委员会-中国工程物理研究院联合基金(批准号: U1530260)、国家自然科学基金(批准号: 51706193)和江苏省高校自然科学研究项目(批准号: 17KJB470014)资助的课题.

摘要: 传质引发的Rayleigh-Bénard-Marangoni对流(RBM对流)对化工传递过程有着显著影响. 但是, 已有的相关研究多集中于气-液体系, 并且有限的针对液-液体系的相关研究尚缺乏对RBM对流演化及其引发的界面扰动行为的深入分析. 因此, 本文基于阴影法设计搭建了竖直狭缝内液-液两相液层间传质过程的RBM对流特性可视化实验平台, 并实验观测了水-甲苯-丙酮三元体系中丙酮组分扩散传质时出现的RBM对流结构以及其向下层水相主体的发展演变过程, 探讨了水相丙酮初始浓度、甲苯相丙酮初始浓度以及甲苯层厚度对RBM对流特性和液-液界面形貌的影响. 研究表明: 在Rayleigh-Taylor不稳定性作用下, 水相上层密度(重力)分层“界面”下凸沉降形成波浪形丘状“界面”, 并随着“界面”处密度与压力失调的加剧而演变成羽状流; 因羽流区“界面”不同浓度梯度引起的传质特性差异, 羽状流又可以演变成弱羽状流和强羽状流两种形态; 当丙酮浓度梯度增大到一定程度后, 近界面处短时间内产生大量RBM对流结构, 且结构间相互影响增强而聚并成对流团, 并随着传质过程的进行, 逐渐演变成独立的强羽状流; RBM对流强度与上下液层丙酮浓度梯度大小呈正相关关系, 且液-液界面粗糙度及其非稳态波动随着丙酮浓度梯度的增加而增大.

English Abstract

    • 物质传递存在于很多化工生产过程中, 如蒸馏、萃取、吸收和解吸等[1-3]. 物质传递会改变两相界面处的流体物化性质, 从而形成密度差和界面张力梯度, 并诱使界面附近产生Rayleigh-Bénard效应[4]和Marangoni效应[5,6], 进而促使界面处的流体产生对流. 这两种效应通常相互耦合、相互影响, 被称之为Rayleigh-Bénard-Marangoni (RBM)效应[7,8]. RBM效应能改变流体流动状态和界面对流演变行为, 导致传质界面形貌变化, 从而影响界面处表面张力分布. 而不均匀的表面张力能够引起界面处流体的湍动, 促进流体界面的更新[9,10], 进而增强化工过程的传质效率[11,12]. 因此深入研究RBM效应对化工非稳态传质过程有重要的实际指导意义[13,14]. 此外RBM效应是强非线性传质流动耦合现象, 是非线性理论良好的研究对象[15,16], 所以深入研究RBM效应还能进一步完善传质与界面科学基础理论, 有着重要的科学意义.

      直观的实验观测是研究RBM效应的一种重要方法, 其中纹影法[17,18]、阴影法[19]以及干涉法[20]为非介入式可视化实验观测方法, 对实验对象的流动及传质过程没有影响, 因而被广泛应用于RBM效应的实验研究. Okhotsimskii和Hozawa[21]运用纹影技术在多种有机溶剂吸收CO2的过程中观测到向液相主体发展的羽状对流; 王勇和张泽廷[22]与Sun等[23]利用纹影法在平行于界面方向上观察了有机溶剂在水平层流运动过程中吸收CO2的传质过程, 得到了滚筒状和多边形的Rayleigh对流结构; 沙勇等[24]利用阴影法研究丙酮-水溶液解吸CO2传质过程, 在垂直于界面方向上观察到了羽状对流结构. 物质传递过程中RBM效应的研究大多采用气液传质体系, 相较之下, 液-液传质体系的研究还很欠缺. Orell和Westwater[25]利用纹影法研究了乙二醇-乙酸-乙酸乙酯3组分的液-液萃取过程, 观察到了有序的多边形和细胞状对流结构; Zhang等[26]通过纹影法利用水-乙酸-甲苯液-液系统研究了双组分扩散引发的界面湍动现象, 发现只有乙酸浓度高于某一值时才会引发界面湍动; Guzun-Stoica等[27]提出了一种实时动态全息干涉法, 能够直接观察液-液过程中的界面湍动现象, 在研究水-乙酸-乙酸丁酯体系时发现了细胞状对流结构; Kostarev等[28]利用干涉法研究异丙醇水溶液中的氯苯液滴, 观察到细胞状对流结构. 可是这些液-液体系的RBM效应研究都局限于定性地观测界面处的湍动现象, 缺乏对界面湍动的定量化研究, 并且所观察到的现象都是界面对流充分发展后的形态[29,30], 而对于RBM对流由界面向主流发展过程的研究十分匮乏, 对其演化规律也缺乏深入认识.

      为此, 本文以水-甲苯-丙酮(溶质)三元体系为研究对象, 基于阴影法建立竖直狭缝内液-液两相液层间传质过程的RBM对流可视化实验平台, 实验观测RBM对流向下层主流区发展过程中出现的典型对流结构, 探索RBM对流向主流区演变发展的规律, 并通过量化主流区流场速度分布、投影图像特征及液-液界面粗糙度, 研究水相、甲苯相丙酮初始浓度以及甲苯液层厚度对RBM对流强度和液-液界面形貌的影响规律.

    • 本文利用阴影法观测竖直玻璃狭缝内液-液两相液层间传质过程中的RBM对流, 实验系统如图1(a)所示, 主要包括点光源、玻璃狭缝、投影屏幕以及高速显微成像系统四部分. 传质过程的不均匀将引起狭缝内工质浓度分布的差异, 进而造成非均匀折射率分布. 因此, 光透过玻璃狭缝内的流体介质后, 会发生不同角度的偏转, 从而在屏幕上形成明暗不均的投影图像. 实验使用显微镜(型号: AT-X M100 PRO)观测投影图像的演变过程, 并采用与显微镜配套安装的高速CCD(型号: FASTCAM SA1.1)实时记录及保存. 本文实验采用激光点光源发出的发散光成像, 在屏幕上得到的是实验观测区域的放大图像, 能够清晰地观察和记录实验出现的微小流动结构.

      图  1  RBM对流的实验系统图 (a) 阴影法实验系统图; (b) 示踪粒子法实验系统图; (c) 玻璃狭缝尺寸图

      Figure 1.  Schematic diagram of the experimental system for RBM convection: (a) Schematic diagram of the experimental system based on shadowgraph method; (b) schematic diagram of the experimental system based on particle tracer method; (c) size of the glass slit.

      此外, 实验中采用粒子示踪法[31]获取详细的流场信息, 实验系统如图1(b)所示. 为此, 下层丙酮水溶液中预先均匀分散质量浓度为0.01%的PSP微球(直径${d_{\rm{p}}} = 5\;{\text{μ}}{\rm{m}}$, 密度${\rho _{\rm{p}}} = 1.04\;{\rm{g}} \cdot {\rm{c}}{{\rm{m}}^{ - 3}}$, 不溶于丙酮、甲苯等有机溶剂), 当其被激光照射时, 借助CCD高速显微系统(显微镜型号: NAVITAR 1-50487A), 可实时观测并记录粒子运动图像, 再结合PIV技术便可定量分析该片区域的流场信息.

      玻璃狭缝由两块光学玻璃构成, 玻璃两侧和底部用U型硅胶垫片密封, 具体尺寸如图1(c)所示. 玻璃狭缝被竖直安装在三维可移动平台上, 置于点光源与屏幕之间, 调节三者相对距离在屏幕上形成清晰放大的投影图像. 实验时先用注射器向玻璃狭缝中注入一定量的丙酮水溶液, 静置一段时间, 待其界面平稳后, 再用两台LSP01-1BH型注射泵从狭缝两侧同时注入甲苯, 并实时拍摄记录投影屏上的投影图像. 拍摄记录时, CCD拍摄速率设为50帧. 除溶质初始浓度外, 保证其他实验参数一致. 实验环境温度T = 20 ℃, 丙酮水溶液注入量为1.6 mL, 甲苯注液速率均设为0.54 mL/min. 实验前将甲苯与去离子水进行预饱和, 以排除实验过程中这两者之间的传质对溶质扩散传质的影响. 实验中所使用的玻璃狭缝为上部开口结构, 从而使扩散至甲苯中的丙酮可以继续向外界大空间环境中不断挥发, 以保证本实验观测的时间范围内丙酮持续由水向甲苯中扩散传质. 本文将注液开始时刻作为起始时刻. 实验所用分析纯(纯度 ≥ 99.7%)甲苯、分析纯(纯度 ≥ 99.7%)丙酮以及去离子水(纯度 ≥ 99.99%)的具体物性参数见表1. 水-甲苯两相界面张力系数$\gamma $和水相溶液密度$\rho $随丙酮体积浓度${\varphi _0}$的变化关系如图2图3所示.

      实验试剂ρ/kg·m–3µ/10–4Pa·s
      99810.04
      甲苯8675.86
      丙酮7903.26

      表 1  实验试剂的物性参数(T = 20 ℃, P = 0.1 MPa)

      Table 1.  The physical parameters of experimental reagents (T = 20 ℃, P = 0.1 MPa).

      图  2  T = 20 ℃时不同水相丙酮浓度下水-甲苯两相间的界面张力系数

      Figure 2.  Interfacial tension coefficient between water and toluene phases under different acetone concentrations of aqueous phase with T = 20 ℃.

      图  3  T = 20 ℃时不同丙酮浓度的水相溶液的密度

      Figure 3.  Density of aqueous solution with different acetone concentrations with T = 20 ℃.

    • 由上述实验原理以及图2图3可知, 由于下层水相中的丙酮向上层甲苯相的不断扩散传质[32]以及该过程存在着不可避免的非均匀性[33,34], 下部水相液层和界面处分别会产生密度和界面张力差异, 进而分别诱发Rayleigh-Bénard效应和Marangoni效应. 这两种效应相互耦合、相互影响, 促使水相内部产生RBM效应, 而由RBM效应产生的水相液层内的RBM对流过程会使其内部各种组分浓度重新分布, 从而使其在光源照射下产生折射率差异, 并反映在投影屏上明暗交替的图像之中. 因此, 通过考察投影图像特征和发展规律即可对RBM对流特性加以研究.

    • 考虑到实验中注液带来的扰动以及液层厚度的不断变化将不利于上部液层RBM对流结构的实验观测与定量化分析. 因此, 本文集中对丙酮组分跨界面传质时下层水相溶液中的RBM对流特性进行研究.

      图4所示, 随着丙酮向上层甲苯液层扩散传质, 靠近界面处的水相丙酮浓度减小, 并且, 由于扩散过程缓慢, 该处近界面处水相中因扩散传质而损失的丙酮组分得不到来自于下部水相的及时补充. 因此, 传质过程初期水相溶液上部出现了浓度分层的现象, 并且由于丙酮浓度不同造成的折射率差异使得投影图像上近界面处水相上下两层之间产生了一条水相内部“界面”. 而由图3可知, 水相中丙酮浓度的分层必然使其内部产生低密度流体支撑高密度流体的(密度)重力分层现象. 此外, 界面处丙酮传质存在着不可避免的非均匀性, 导致水相内部“界面”处产生不规则的扰动, 如图4(b)所示.

      图  4  传质过程引起的密度分层示意图与实验结果 (a) 密度分层示意图; (b) 实验图像

      Figure 4.  Schematic diagram and experimental result of density stratification caused by mass transfer: (a) Schematic diagram of density stratification; (b) experimental image.

      在这样的条件下, 水溶液内部“界面”上便诱发产生了Rayleigh-Taylor不稳定性[35], “界面”上的扰动单元会造成界面上的密度(重力)梯度与静压梯度失调(界面弯曲造成密度梯度和压力梯度不作用在同一水平线上), 导致上层液体克服流体黏性阻力不均匀地侵入下层液体, 并在界面各处产生斜压转矩而产生涡流[36], 如图5所示. 涡流会使扰动“界面”的凹凸程度不断扩大, 从而使扰动“界面”在下降的过程中呈现出向下凸出的非均匀波浪形丘状“界面”, 如图6所示.

      图  6  丘状“界面”的形成过程(水相丙酮初始体积浓度${\varphi _0} = 5\% $, 甲苯相丙酮初始体积浓度${\varphi _1} = 0\% $) (a) t = 13 s; (b) t = 30 s; (c) t = 36 s; (d) t = 42 s

      Figure 6.  The forming process of the mound “interface” (the initial volume concentration of acetone in aqueous phase ${\varphi _0} = 5\% $, the initial volume concentration of acetone in the toluene phase ${\varphi _1} = 0\% $): (a) t = 13 s; (b) t = 30 s; (c) t = 36 s; (d) t = 42 s.

      图  5  密度分层引起的Rayleigh-Taylor不稳定性示意图 $\omega $是涡流, P是压力, $\rho $是密度, u是速度, g是重力加速度; 粗的环形箭头表示涡旋产生的速度场

      Figure 5.  Schematic diagram of Rayleigh-Taylor instability caused by density stratification. $\omega $ is vorticity, P is pressure, $\rho $ is density, u is velocity and g is acceleration of gravity; the thick circular arrows represent the velocity field created by the vortex.

      随着丘状“界面”结构的不断下移, “界面”处密度梯度与压力梯度失调进一步加剧, “界面”各处涡流强度不断提升, 进而使扰动界面不稳定性进一步加剧, 部分丘状 “界面”下端的前凸幅度进一步增大. 随着时间的推移, 当丘状 “界面”下端前凸处从山丘状演变成柱状时, 羽状流形成(如图7所示).

      图  7  羽状流的演变过程 (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 0 s; (b) t = 18 s; (c) t = 20 s; (d) t = 22 s

      Figure 7.  The evolution of the plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 0 s; (b) t = 18 s; (c) t = 20 s; (d) t = 22 s.

      实验中, 在不同初始丙酮浓度下, 由于密度梯度不同引起的Rayleigh-Taylor不稳定性强度不同, 使得羽状流呈现出弱羽状流和强羽状流两种典型形态. 弱羽状流出现在水相溶液中丙酮浓度比较低的情况(如图8(a), t = 348 s时). 由于此时水相与甲苯相内的丙酮浓度梯度较小, 传质较弱, 其产生的水相液层内竖直密度差较小, 造成Rayleigh-Taylor不稳定性强度较弱, 羽状流向下延伸较慢. 另外, 起初羽状流内外流体的丙酮浓度差异也较小. 因此, 在羽状流缓慢的下降行程中, 羽流区下端流体可与其四周的大空间流体实现充分的传质, 使羽状流下端流体浓度逐渐增大, 直至与其周围主流区流体的丙酮浓度一致. 由于浓度差异消失, 这样由浓度差异导致的液体折射率差异也消失, 这部分羽状流结构会从屏幕上逐渐变浅并最终消失, 从而呈现出颜色由根部向下端逐渐变浅的羽毛状结构, 即弱羽状流, 如图8(d)所示.

      图  8  弱羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 348 s; (b) t = 356 s; (c) t = 363 s; (d) t = 370 s

      Figure 8.  The evolution of the weak plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 348 s; (b) t = 356 s; (c) t = 363 s; (d) t = 370 s.

      为进一步认识羽状流的产生机理, 采用粒子示踪法获取羽状流区域的详细流场信息, 如图9所示. 羽状流区域的流场可划分为羽流区和涡流区两部分. 其中, 涡流区分布在羽流区两侧, 羽流区内流体自上而下运动, 两个涡流区内的流体做涡旋运动且方向相反, 漩涡区内涡量绝对值高于其他区域. 与羽状流的下降速度相比, 涡流区流体流速明显较小, 导致涡流区产生高静压, 羽流区产生低静压, 两者间的压差产生挤压力, 将羽状流颈部压细, 并将流体压向末端. 受到下层相对静止流体的阻碍作用, 快速流向下端的流体流速明显减缓. 这样, 在颈部两侧压差、下层流体阻碍的耦合作用下, 羽状流末端出现体积增大和扁平化的趋势. 特别是(如图10所示), 当水相中丙酮浓度较高时, 羽状流下降快, 导致上述两种效应的耦合作用增强, 使羽状流末端发展成宽扁的弧形帽状, 从而使羽流区整体投影图像呈现出倒蘑菇状(如图10(c), t = 26 s时), 此时, 强羽状流形成.

      图  10  强羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 22 s; (b) t = 23 s; (c) t = 26 s

      Figure 10.  The evolution of the strong plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 22 s; (b) t = 23 s; (c) t = 26 s.

      图  9  羽状流的速度矢量及涡量云图(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $)

      Figure 9.  The velocity vector and vorticity contours of the plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $).

      强羽状流在向下发展过程中, 周边水相内的丙酮组分不断向羽流区进行传质, 导致其内部丙酮组分浓度升高. 因此, 在强羽状流产生的后期, 由于靠近帽状头部的下半部分沉降时间较长, 上述传质过程进行得比较充分, 头部丙酮浓度率先与周边水相内的丙酮浓度趋于一致(折射率趋于一致), 在投影上逐渐消失. 此时, 对流结构从蘑菇云状 变成了羽毛状, 即强羽状流演变为弱羽状流, 如图11 所示.

      图  11  强羽状流向弱羽状流的演变过程(${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $) (a) t = 228 s; (b) t = 234 s; (c) t = 238 s; (d) t = 245 s

      Figure 11.  The evolution of the strong plume flow to the weak plume flow (${\varphi _0} = 15\% $, ${\varphi _1} = 0\% $): (a) t = 228 s; (b) t = 234 s; (c) t = 238 s; (d) t = 245 s.

      当丙酮初始浓度梯度较大(如${\varphi _0} = 30\% $, ${\varphi _1} = $0%)时, 丙酮向甲苯传质得很剧烈, 传质导致的密度差以及初始扰动的频率很大, 其产生的Rayleigh-Taylor不稳定性很强, 造成液-液界面处产生大量微小的细胞状对流结构. 由于这些对流结构数量大、彼此间距很小, 且其两侧的夹带流相互作用强, 导致这些对流结构在向下延伸时, 快速聚并成大对流团(如图12红色方框所示), 而不能充分发展成独立的羽状流形态.

      图  12  传质初期对流结构的聚并过程(${\varphi _0} \!=\! 30\% $, ${\varphi _1} \!=\! 0\% $)  (a) t = 15 s; (b) t = 20 s; (c) t = 24 s; (d) t = 30 s

      Figure 12.  Convergence process of convective structure at the beginning of the mass transfer (${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $): (a) t = 15 s; (b) t = 20 s; (c) t = 24 s; (d) t = 30 s.

      随着传质过程的进行, 水相与甲苯相丙酮浓度梯度逐渐降低, 传质过程减弱, Rayleigh-Taylor不稳定性减弱, 导致界面处产生的对流结构数量减少、间距增大、相互叠加现象逐渐减弱. 这时, 如图13所示, 卷吸对流团的混乱程度降低, 且其体积逐渐变小直到消失, 最终形成独立的羽状流结构(如图13(d), t = 350 s时), 即强羽状流.

      图  13  对流团的消失以及强羽状流的出现(${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $)  (a) t = 101 s; (b) t = 196 s; (c) t = 259 s; (d) t = 350 s

      Figure 13.  The disappearance of convective cloud and the appearance of the strong plume (${\varphi _0} = 30\% $, ${\varphi _1} = 0\% $): (a) t = 101 s; (b) t = 196 s; (c) t = 259 s; (d) t = 350 s.

      由此可知, 上述各种RBM对流结构主要是由界面处传质导致的界面张力和密度差异引起的, 同时可随着界面处传质特性的时空变化而呈现出不同的形态并产生相互演化.

    • 为了探讨溶质初始浓度对RBM对流强度的影响, 本研究借助图像处理手段, 采用对流结构延伸速度v以及羽状流个数n对RBM对流特征进行定量表征. 其中, v定义为第一个RBM对流结构向下延伸的速度, 表征了其向下层主流区发展的快慢; n定义为羽状流的羽流区“枝杈”个数.

      图14给出了水相丙酮初始浓度对对流结构延伸速度的影响. 如图14所示, 随着水相丙酮初始浓度的增大, 水相与甲苯相丙酮浓度梯度增大, 界面处丙酮传质强度增强, 造成如图4所示的水相近界面处的密度分层现象加剧, 且高、低密度层间“界面”处的扰动增强, 从而诱发了更强的Rayleigh-Taylor不稳定性, 导致水相近界面处上层高密度液体更快地侵入下层低密度液体, 因此, 对流结构向下延伸的速度加快.

      图  14  水相丙酮初始浓度对第一个RBM对流结构向下延伸速度的影响

      Figure 14.  The influence of initial concentration of acetone in aqueous phase on the elongation velocity of the first RBM convective structure.

      图15给出了水相丙酮初始浓度对羽状流数量的影响. 由图可知, 随着水相丙酮初始浓度的增加, 初始阶段界面处的传质更为剧烈, 造成水相内近界面处的扰动频率升高[37], 使得羽状流产生的“种子”数量增加, 羽状流数目随之升高; 随着时间的推移, 不同强度的羽状流先后消失-再生-消失, 如此循环便造成羽流区“枝杈”个数出现波动, 并且, 随着丙酮初始浓度的增加, 羽状流消失-再生-消失的演变周期缩短, 使其羽流区“枝杈”数量波动不断增大. 另一方面, 随着羽流区“枝杈”数量的增多, 其相对间距减小, 造成彼此间聚并增强[38]. 因此, 羽流区“枝杈”数量的增幅会随着初始浓度的增大而减小(${\varphi _0} > 10\% $). 特别是, 当羽状流聚并速度大于其生成速度时, 羽流区“枝杈”总数甚至会变少(如${\varphi _0} = 20\% $${\varphi _0} = 15\% $相比).

      图  15  水相丙酮初始浓度对羽状流数量的影响

      Figure 15.  The influence of initial concentration of acetone in aqueous phase on the number of the plumes.

      此外, RBM对流结构的形成和演变都会造成近界面流场的变化, 从而导致液-液界面形貌产生变化. 因此, 采用界面粗糙度RC [39]定量表征液-液界面形貌

      ${RC} = \frac{1}{{{N_1}}}\sum\limits_{i = 1}^{{N_1}} {\left| {{y_i} - \bar y} \right|} ,$

      式中, N1为界面投影下边界的像素点数, yi$\bar y$分别是各像素点的纵坐标值以及纵坐标平均值. 此外, 采用传质开始后某一时间段内不同时刻相界面粗糙度的标准差$\sigma $(RC )表征该时间段内粗糙度的波动程度:

      $\sigma (RC) = \sqrt{ \displaystyle\frac{1}{N_2 - 1} \sum_{i = 1}^{N_2} (RC_i - \overline{RC})^2}, $

      其中, N2为某一时间段内取样时刻数(本文取样时间段为0—400 s, 取样时刻数为21, 取样时刻间的时间步长为20 s), RCi$\overline {{RC}} $分别是各个时刻下粗糙度值以及取样时间段内所有时刻粗糙度的算数平均值.

      由上述讨论可知, 水相丙酮初始浓度增大时, 由于对流结构数量变多且发展速度加快, 近界面处流场扰动增强, 从而导致不同时刻下的液-液界面变得更加粗糙, 如图16(a)所示. 此外, 同一工况下, 随着时间推移, 对流结构类型和强度不断变化, 导致界面粗糙度也在不断地变化, 从而引起界面波动. 而随着丙酮初始浓度增大, 同一时间段内(如0—400 s)RBM对流结构的种类愈加丰富, 且不同种类之间相互演化行为更为频繁, 造成对流结构形态更为复杂且数量波动变大, 导致界面粗糙度的数值波动会变大, 如图16(b)所示.

      图  16  水相丙酮初始浓度对水-甲苯界面形貌的影响 (a)界面粗糙度; (b)界面波动程度

      Figure 16.  The influence of initial concentration of acetone in aqueous phase on water-toluene interface morphology: (a) Interfacial roughness; (b) the degree of interface fluctuation.

    • 甲苯相丙酮初始浓度${\varphi _1}$的变化可以改变水-甲苯界面两侧的丙酮浓度梯度, 进而影响丙酮的传质强度, 故而成为影响RBM对流特性的重要因素. 为此, 本文研究了甲苯相丙酮初始浓度对RBM对流特性的影响, 结果如图17所示. 由图可知, 在保证水相丙酮初始浓度不变的前提下, 随着甲苯相初始丙酮浓度的提高, 界面两侧丙酮浓度梯度减小, 界面传质强度降低, 由其引发的RBM对流现象逐渐减弱. 特别当甲苯相初始丙酮浓度达到${\varphi _1} = $15%时(见图17(d)), 下部水相液层内部的RBM对流现象几乎消失. 这是由于此时界面两侧丙酮初始浓度梯度较小, 跨界面传质较弱, 在水相溶液内部产生的竖直浓度梯度较小. 再者, 下层水相溶液内部的不同浓度层间也一直存在着溶质扩散, 从而进一步减小了水相溶液内部竖直方向上的浓度梯度. 这些因素共同使得界面丙酮扩散传质而产生的界面张力和密度(重力)梯度无法克服黏性阻力而诱发RBM对流的产生.

      图  17  不同甲苯相丙酮初始浓度下的投影图像(t = 35 s, ${\varphi _0} = 20\% $) (a) ${\varphi _1} = 5\% $; (b) ${\varphi _1} = 7.5\% $; (c) ${\varphi _1} = 10\% $; (d) ${\varphi _1} = 15\% $

      Figure 17.  Schlieren images under different initial concentrations of acetone in the toluene phase (t = 35 s, ${\varphi _0} = 20\% $): (a) ${\varphi _1} = 5\% $ (b) ${\varphi _1} = 7.5\% $; (c) ${\varphi _1} = 10\% $; (d) ${\varphi _1} = 15\% $.

      为进一步展示甲苯相丙酮初始浓度对RBM对流特性的影响, 图18图20分别给出了${\varphi _0} = 20\% $时不同甲苯相丙酮初始浓度下RBM对流结构延伸速度、羽状流数量和水-甲苯界面形貌的变化. 如上所述, 甲苯相丙酮初始浓度提高引发界面处丙酮传质强度降低, 造成水相近界面处的密度分层现象减弱, 高、低密度层间“界面”处的扰动频率减小. 因此, RBM对流结构延伸速度降低(见图18), 羽状流数目减少(见图19). 再者, 界面处丙酮传质强度的降低造成RBM对流紊乱程度变小, 近界面处流场扰动减弱, 从而使得液-液界面粗糙度降低且粗糙度数值的波动减小, 如图20所示.

      图  18  甲苯相丙酮初始浓度对第一个RBM对流结构向下延伸速度的影响

      Figure 18.  The influence of initial concentration of acetone in the toluene phase on the elongation velocity of the first RBM convective structure.

      图  19  甲苯相丙酮初始浓度对羽状流数量的影响

      Figure 19.  The influence of initial concentration of acetone in the toluene phase on the number of the plumes.

      图  20  甲苯相丙酮初始浓度对水-甲苯界面形貌的影响 (a)界面粗糙度; (b)界面波动程度

      Figure 20.  The influence of initial concentration of acetone in the toluene phase on water-toluene interface morphology: (a) Interfacial roughness; (b) the degree of interface fluctuation.

    • 值得注意的是, 由于本实验中上部甲苯液层的体积有限, 因此丙酮向甲苯相的扩散传质将使得甲苯相中的丙酮浓度(分压力)升高, 从而降低水相与甲苯相间的丙酮浓度梯度, 进而影响到RBM对流特性. 为此, 图21图22分别给出了${\varphi _0} = 20\% $时甲苯液层厚度d对RBM对流特性的影响.

      图  21  甲苯层厚度对羽状流数量的影响

      Figure 21.  The influence of thickness of toluene layer on the number of the plumes.

      图  22  甲苯层厚度对水-甲苯界面形貌的影响

      Figure 22.  The influence of thickness of toluene layer on water - toluene interface morphology.

      由图可知, 在同样的丙酮跨界面扩散传质强度下, 液层厚度大(体积大)的甲苯相中丙酮浓度的上升幅度要小于液层厚度小(体积小)时的上升幅度. 因此, 较厚的甲苯液层能够维持界面内外较大的丙酮浓度差异, 丙酮向上传质的强度更大, 其引发的Rayleigh-Taylor不稳定性和RBM对流现象更加明显. 因此, 当甲苯液膜较厚时, RBM对流羽状流数目较多(见图21), 水-甲苯界面粗糙度较大且其数值波动更为明显(见图22).

    • 本文基于阴影法设计搭建了竖直狭缝内液-液两相液层间传质过程的RBM对流特性可视化实验平台, 着重研究了水-甲苯-丙酮三元体系中丙酮组分扩散传质时出现的RBM对流结构及其向下层水相主体的发展演变过程, 通过改变水相、甲苯相中丙酮初始浓度和甲苯液层厚度, 探讨了其对RBM对流特性以及液-液界面形貌的影响, 获得以下主要结论.

      1) 丙酮跨液-液界面传质引起近界面处水相流体内的密度(重力)分层和扰动, 造成密度(重力)分层“界面”上产生Rayleigh-Taylor不稳定性, 进而促使密度(重力)分层“界面”下凸沉降形成丘状“界面”, 随着丘状“界面”处密度与压力失调的加剧而演变成羽状流;

      2) 因中心羽流区与两侧涡流区所产生的压差挤压作用和下层相对静止流体的阻碍作用, 羽流区由柱状发展成倒蘑菇状, 从而形成强羽状流, 并且随着羽流区与主流区传质过程的进行, 而进一步演化成羽毛状, 形成弱羽状流; 当丙酮浓度梯度增大到一定程度后, 近界面处短时间内产生大量RBM对流结构, 且结构间相互影响增强而聚并成对流团, 并随着传质过程的进行, 逐渐演变成独立的强羽状流;

      3) RBM对流强度与上下液层丙酮浓度梯度大小呈正相关关系: 随着界面处丙酮浓度梯度的增大, RMB对流结构延伸速度逐渐加快; 羽状流枝杈数量增多且传质过程中其数值波动变大, 由此近界面处流场扰动增强, 导致不同时刻液-液界面粗糙度相对增大, 且界面粗糙度的数值波动变大.

      衷心感谢东南大学能源与环境学院邓梓龙博士在实验系统设计、搭建方面的指导和有益建议.

参考文献 (39)

目录

    /

    返回文章
    返回