Y型微通道内双重乳液流动破裂机理*

俞炜 邓梓龙† 吴苏晨 于程 王超

1) (东南大学能源与环境学院, 能源热转换及其过程测控教育部重点实验室, 南京 210096)

2) (扬州大学水利与能源动力工程学院, 扬州 225127)

基于体积分数法建立了Y型微通道中双重乳液流动非稳态理论模型, 数值模拟研究了Y型微通道内双重乳液破裂情况, 详细分析了双重乳液流经Y型微通道时的流场信息以及双重乳液形变参数演化特性, 定量地给出了双重乳液流动破裂的驱动以及阻碍作用, 揭示了双重乳液破裂流型的内在机理. 研究结果表明: 流经Y型微通道时, 双重乳液受上游压力驱动产生形变, 形变过程中乳液两端界面张力差阻碍双重乳液形变破裂, 两者正相关; 隧道的出现将减缓双重乳液外液滴颈部收缩速率以及沿流向拉伸的速率, 并减缓了内液滴沿流向拉伸的速率, 其对于内液滴颈部收缩速率影响不大; 隧道破裂和不破裂工况临界线可以采用幂律关系式 进行预测, 隧道破裂和阻塞破裂工况临界线可以采用线性关系 描述; 与单乳液运动相图相比, 双重乳液运动相图各工况的分界线关系式系数 和 均相应增大.

1 引 言

双重乳液是一种分散相液滴中包裹着更小液滴的结构化流体, 其独特的壳-核结构使之具有更高的灵活性和可控性, 在化学反应[1,2]、功能材料合成[3-5]、药物输运[6,7]等领域有着广泛的应用前景.自20年前问世以来, 微流控技术获得了长足发展,可以实现对多相流体在微尺度上的精确操控, 并且具有装置高度集成、实验试剂消耗少等优势. 因此,微流控技术已经成为制备双重乳液的一种有效手段, 并引起了众多学者的广泛关注[8-10]. 然而, 为了保证对多相流体的精确和稳定操控, 微流控装置内各相工质通常处于流速较低的稳定层流状态, 这在一定程度上限制了微流控乳液制备的产率. 为了实现乳液的高通量制备, 研究者们多利用微流控分裂结构对制备生成的乳液进行分裂倍增, 从而实现乳液产率的提升[11,12]. 因此, 深入认识微流控分裂结构中乳液分裂过程的流体动力学机理十分必要, 这不仅对完善乳液多相流体动力学理论具有重要的科学意义, 并且可为工程中微流控乳液制备产率的提升提供关键的技术支撑.

为实现乳液尺寸的精确操控和乳液产量的倍增, 研究者们设计出多种微流控分裂流道结构, 如十字通道[13]、T 型通道[14-23]、Y 型通道[24-29]、带阻碍结构通道[14,30]等. Link等[14]首次将T型微通道用于液滴被动破裂, 并将一系列T型通道相连形成树状网络, 可在不改变高分散相体积分数的情况下实现小液滴的大量制备. 自Link以后, 众多学者继续致力于通过实验研究T型微通道内液滴破裂特性和机理. 研究结果表明: 液滴在T型微通道内的被动破裂存在两种流型, 分别为隧道破裂流型和阻塞破裂流型[16]; T型微通道中液滴破裂与否, 主要取决于液滴初始长度和连续相流体的毛细数, 两者之间存在指数关系, 以此表征破裂与不破裂之间的流型临界线[14,16,21]. 考虑到已有的微流控装置制备双重乳液存在产率低的问题, 并且T型微通道不易实现双重乳液分裂, Abate和Weitz[26]将Link[14]等的工作进行了拓展, 将Y型通道作为树状网络的单元结构, 实现了液滴和双重乳液的高通量制备. 其研究结果表明, 主通道中连续相流速增大,双重乳液形变加快, 有助于双乳液分裂. 最近, 不同于传统的2D光刻技术, 新型的3D对称毛细管被用于液滴和双重乳液的分裂, 可实现高质量的乳液可控倍增[31]. 然而, 针对双重乳液在特征几何结构微通道内的流动破裂特性实验研究还未见到. 由于双重乳液相比单液滴具有更复杂的核壳结构[32],因此有必要对双重乳液变形破裂过程中内外液滴界面的相互影响开展深入讨论与分析.

与实验研究相比, 计算流体动力学方法为深入探究双重乳液多相流动特性提供了另一个研究途径, 此方法可以详细给出实验研究中较难获得的速度分布、压力梯度以及界面形变的实时情况[33]. 而充分了解这些流场和压力场的演化规律将有助于优化微流控装置的设计. 目前, 相场方法已被用于模拟研究T型分叉处的液滴流动破裂特性[15,23]. 研究表明, 液滴是否破裂与连续相和离散相黏度比、毛细数和液滴长度有着重要关联[15,23], 在大黏度比和大毛细数下液滴更易于破裂, 而液滴尺寸减小使得液滴不易破裂[23]. 此外, Carlson 等[25]运用该方法研究Y型分叉微通道中液滴破裂特性,其工况分为破裂工况和不破裂工况, 两种工况的临界线与液滴长度以及毛细数相关, 该研究认为Rayleigh-Plateau不稳定性是液滴破裂的驱动机制. 梁宏等[27]也采用此方法模拟了液滴在Y型分叉微通道中的迁移过程. 结果表明: 当毛细数足够小时, 液滴则滞留在分叉口处, 不发生破裂. 除了相场方法, 已有基于体积分数法(VOF)方法研究对称T型微通道中液滴破裂和不破裂特性的工作, 与实验研究结果类似, 发现了三种不同破裂机制, 分别为隧道破裂、不连续阻塞破裂和阻塞破裂[18,19,22], 液滴破裂特性与毛细数以及液滴长度密切相关[17-19,22], 毛细数增大会加快液滴破裂[20,22].随后, 该方法还被用于研究不同分叉角度微通道中的液滴破裂特性, 结果表明, 两个子通道之间的夹角越小, 液滴分裂得更快[29]. 此外, 已有文献基于数值方法研究双重乳液在不同制备结构中的生成过程[34−39]、在剪切拉伸流中的变形破裂过程[32,40,41]、在剪切拉伸流中双重乳液的碰撞过程[42,43]等, 但尚缺乏采用数值方法全面深入了解双重乳液在Y型分叉微通道内流动破裂特性的研究.

综上所述, 具有Y型分叉结构的微通道为精确操控乳液尺寸和高通量制备微小乳液提供了有效途径. 与T型微通道相比, Y型微通道更利于单乳液分裂, 并已被成功应用于双重乳液分裂实验.因此, 充分了解Y型微通道内双重乳液破裂这一过程有着重要的科学意义和工程应用价值. 然而, 现有研究主要集中于单乳液在Y型分叉通道中流动破裂特性的研究. 此外, 与单乳液相比, 双重乳液中内液滴的存在增加了内外液滴间的相互作用[32,40,41]. 目前, 关于内液滴对双重乳液在Y型微通道中流动破裂特性影响的研究还有所不足, 其内在机理有待进一步阐明. 为此, 本文基于VOF方法模拟研究了双重乳液在Y型微通道中破裂的机理, 详细分析了双重乳液流经Y型微通道时的流场信息以及双重乳液形变参数的变化, 并给出了Y型微通道内双重乳液运动相图.

2 数学模型

本文建立了三相流体流动的二维数学模型以研究双重乳液在Y型分叉微通道内的动力学行为特性. 计算区域如图 1 所示, 通道尺寸为 w0= 200,w1= 0.75w0, l0= 12w0, l1= 10w0, l2= 8w0,=30°. 内相流体 (密度, 黏度)、中间相流体 (密度, 黏度)和外相流体 (密度, 黏度)由主通道入口处注入并且在表面张力的作用下形成双重乳液, 在外相流体的带动下向下游运动. 内液滴和外液滴的尺寸由内相和中间相流体的注入时间控制.

2.1 控制方程

不相混且不可压缩的三相流体流动可用VOF表征. 在流场中定义一个体积率函数. 计算网格中= 0, 表示网格内不存在 j相流体; 网格中= 1, 表示网格内充满 j相流体; 0 << 1表示网格内为相界面. 此方法通过研究网格中每一相的体积率函数来确定相界面, 能捕捉各种不规则形状的自由表面, 可较好地模拟双重乳液在流场中的生成[34-39]、变形[32,40,41]、碰撞[42,43]等复杂的变化. 在每个网格中, 各相流体的体积率函数总和为1, 即. 在整个计算区域内满足以下控制方程:

图1 计算区域示意图Fig. 1. Geometrical description of the numerical domain.

体积率函数方程

式中t为时间(s), U为速度(m/s), p为压力(Pa),为密度 (kg/m3),为黏度 (Pa·s), g 为重力加速度 (m/s2),为表面张力 (N/m3). 因为通道尺度为微米级, 所以重力加速度的影响可忽略不计.和可由下式计算:

表1 数值模拟中各相流体的物性参数Table 1. The properties of the fluids used for numerical simulation.

2.2 边界条件

因为数值模拟中雷诺数小于60, 所以本文采用层流模型. 主通道入口处采用周期性边界条件以控制各相流体注入通道的体积. 如图2所示, 主通道入口处流体速度为匀速u0. 双重乳液形成的具体步骤如下: 步骤 (a), 当 0 < t < tm, 中间相流体流入; 步骤 (b), 当 tm< t < tm+i, 内相和中间相流体同时流入 (win= 0.8w0); 步骤 (c), 当 tm+i< t

通道出口处为一个大气压 (1 atm = 1.013 ×105Pa). 通道的内壁面采用无滑移边界条件

2.3 数值求解方法

图2 一个生成周期内双重乳液生成过程Fig. 2. Schematic of the generation process of double emulsion in one generation cycle.

采用有限体积差分法同时离散求解VOF模型中的控制方程. 求解时, 采用二阶迎风差分格式离散动量方程, 压力场与速度场之间的耦合采用SIMPLE算法, 压力场的离散采用体积力分数法,相界面采用分段线性界面重构法(PLIC)[45]进行几何重构. 控制方程离散后, 采用Gauss-Seidel方法进行迭代求解. 为了保证计算的收敛性和收敛速度 , 各 项 松 弛 因 子 设 为 : 0.2(压 力 ), 0.3(密 度),0.3(体积力)和0.2(动量). 采用相对准则判断收敛,即当每个计算步长内各计算单元中所有变量的相对残差之和小于初始值0.1%时判定计算收敛. 计算中的时间步长依据全局库朗数不大于0.2的规则进行自动调整.

本文采用四边形网格对图1所示的计算区域进行网格划分. 对近壁面附近的网格进行局部加密以确保能捕捉到双重乳液与壁面之间的液膜. 计算前采用四套不同的网格数进行网格独立性检验. 如图3所示, 当计算单元总数大于315610时, 随着网格数的增加, 双重乳液形貌基本重合. 考虑计算成本, 本文采用的网格数为315610.

图3 当 t* = 3 时, 不同网格数下双重乳液的形貌 (Ca =0.01, Voi = 1.3, l* = 1.6)Fig. 3. Grid independence test results at t* = 3 (Ca = 0.01,Voi = 1.3, l* = 1.6).

2.4 模型验证

为验证本文所建立的数学模型的正确性, 依据文献[41]中的实验, 基于上文的二维数学模型数值模拟了如图4所示的剪切流场中双重乳液形变, 并且与实验结果进行了比对. 模拟中内、中、外三相流体分别为去离子水 (20 ℃ 下= 0.001 Pa·s,998 kg/m3)、Ucon 润滑油(20 ℃ 下= 0.125 Pa·s,= 991 kg/m3)和不同分子量的硅油混合物(20 ℃ 下= 5.2 Pa·s,= 989 kg/m3). 表面张力系数= 0.024 N/m,= 0.003 N/m. 计算区域的几何尺寸为 W × H = 20Ro× 8Ro, Ro为外液滴半径. 如图5所示, 对于双重乳液的形貌和内、外液滴的 Taylor[46]稳态形变参数 D(Di= (Li– Bi)/(Li+ Bi), Do= (Lo– Bo)/(Lo+ Bo)), 数值模拟结果与实验结果吻合较好, 这表明本文的数学模型能有效地预测Y型微通道内双重乳液的形变.

图4 剪切流场下双重乳液形变研究示意图 (a) 计算区域示意图; (b) 双重乳液形变参数示意图Fig. 4. Schematic of deformed double emulsion in steady shear flow: (a) Schematic of computational domain; (b) schematics of deformation parameters of the inner and outer droplets, respectively.

图5 模拟结果与实验结果[41]对比 (a) 双重乳液的形变参数 D 随 Ca 的变化; (b) 双重乳液形貌对比Fig. 5. Comparison of steady deformation of double emulsion between simulation and experiment[41]: (a) Steady deformation of double emulsion in the function of Ca; (b) comparison of droplet morphology reconstructed from numerical simulation with experimental snapshots.

3 结果分析

基于上述多相流模型, 本文开展了双重乳液在Y型分叉微通道内流动特性的研究. 在文献[16,22, 23]中, 液滴在流经T型分叉微通道时会呈现三种典型流型, 分别为阻塞破裂 (obstructed breakup)、隧道破裂 (tunnel breakup)和不破裂(non-breakup)流型. 与液滴情况类似, 双重乳液流经Y型分叉微通道时, 同样存在这三种典型的流型. 本文针对这三种典型流型的动力学信息(界面形貌、压力分布、速度分布等)瞬时演化规律进行详细分析, 以探究双重乳液破裂过程中的驱动和阻碍作用变化情况, 揭示双重乳液破裂的机理.

3.1 阻塞破裂

当双重乳液流经Y型分叉处时, 上游压力驱动双重乳液形变, 此时双重乳液前端以及尾部界面形貌发生变化(图6(a)), 而双重乳液界面张力将阻碍该过程的发生, 该过程中速度场演化如图6(b)所示. 为阐明上游驱动压力以及双重乳液界面张力之间相互关系, 图6—8分别给出了阻塞破裂工况中流场整体压力分布情况、双重乳液界面张力以及上游压力瞬时演化曲线. 如图6(a)所示, 可将双重乳液流经Y型分叉处分为三个阶段: entering阶段、squeezing阶段以及post-breakup阶段.

图6 阻塞破裂工况 Y 型微通道中压力场与相界面演化 (Ca = 0.01, Voi = 1.3, l* = 2.1)Fig. 6. Evolution of the interface profile and pressure field during obstructed breakup in a Y-junction (Ca = 0.01, Voi = 1.3, l* = 2.1).

图7 阻塞破裂工况乳液前端及尾部界面张力演化情况 (Ca = 0.01, Voi = 1.3, l* = 2.1) (a) 乳液前端界面张力; (b) 乳液尾部界面张力; (c) 乳液前端与尾部界面张力之差; (d) 特征时刻乳液前端与尾部界面张力的示意图Fig. 7. Evolution of the pressure for obstructed breakup (Ca = 0.01, Voi = 1.3, l* = 2.1): (a) The Laplace pressure of the forefront droplet interface; (b) the Laplace pressure of the rear droplet interface; (c) the Laplace pressure difference between the forefront and rear droplet interfaces; (d) schematics of and at different times.

当双重乳液完全破裂, 进入post-breakup阶段, 双重乳液尾部回缩, 界面曲面半径减小, 尾部界面张力增加, 从而流阻减小, 上游压力pinlet相应降低. 该预测与图8中pinlet变化规律相一致.

图 8 阻 塞 破 裂 工 况 入 口 与 出 口 压 力 演 化 情 况 (Ca =0.01, Voi = 1.3, l* = 2.1)Fig. 8. Evolution of the inlet pressure and outlet pressure of the Y-junction for obstructed breakup (Ca = 0.01, Voi =1.3, l* = 2.1).

为深入理解双重乳液破裂的机理, 采用无量纲特征参数定量描述squeezing阶段中双重乳液形貌(包括内液滴以及外液滴)的演化过程(见图9),分别为外液滴颈部厚度、外液滴前端运动距离、内液滴颈部厚度和内液滴前端运动距离.以及分别代表外液滴和内液滴前端接触分叉尖角的时刻, 而和分别为所对应时刻的外液滴和内液滴颈部厚度. 如图9(a)所示,外液滴颈部厚度随时间减小直至双重乳液破裂, 插图中给出了和的对数坐标图.根据指数标度律可判定squeezing阶段包含两个子过程. 子过程 I 中和呈线性关系, 而在子过程 II中两者呈指数关系~()4/5, 这与T型微通道中液滴阻塞破裂的特性类似[23]. 两个子过程中颈部厚度随时间的变化规律不同, 这主要是由于外液滴尾部界面张力发生了改变. 与颈部厚度变化不同的是, 外液滴前端运动 距 离 随 时 间 线 性 增 加 , 有≈(见图 9(b)). 此时, 外液滴前端以 2u0/3 恒定速度运动, 这主要由于该工况中子通道完全被外液滴阻塞所导致. 图9(c)和图9(d)分别给出了内液滴颈部厚度和内液滴前端运动距离随时间的演化曲线, 其演化特性与外液滴一致. 内液滴颈部厚度随时间缩小, 也分为两个子过程,

图9 双重乳液无量纲特征参数在 squeezing 阶段内演化情况 (Ca = 0.01, Voi = 1.3, l* = 2.1) (a)外液滴颈部厚度 , 插图中给出了 与 的对数坐标图; (b) 外液滴前端运动距离 , 插图中给出了 与 的对数坐标图;(c) 内液滴颈部厚度 , 插图中给出了 与 的对数坐标图; (d) 内液滴前端运动距离 , 插图中给出了 与的对数坐标图Fig. 9. Evolution of the dimensionless characteristic parameters in the squeezing stage for obstructed breakup (Ca = 0.01, Voi = 1.3,l* = 2.1): (a)The neck thickness of outer droplet , inset is the same data as log( ) versus log( ); (b) the distance travelled by the tip of outer droplet , the same data as log( ) versus log( ); (c) the neck thickness of inner droplet , inset is the same data as log( ) versus log( ); (b) the distance travelled by the tip of inner droplet , the same data as log( ) versus log( ).

3.2 隧道破裂

从流动行为角度分析, 隧道破裂流型在entering和post-breakup阶段与阻塞破裂流型非常类似,但在squeezing阶段两者有明显区别(见图10). 在squeezing阶段时刻, 隧道破裂流型中乳液与子通道壁面间开始出现两条隧道. 从压力演化角度分析, 隧道破裂流型在 entering和 squeezing阶段均展现出一定的不同(见图11和图12).

图10 隧道破裂工况 Y 型微通道中压力场与相界面演化 (Ca = 0.01, Voi = 1.3, l* = 1.3)Fig. 10. Evolution of the interface profile and pressure field during tunnel breakup in a Y-junction (Ca = 0.01, Voi = 1.3, l* = 1.3).

与阻塞破裂类似, 进入 post-breakup 阶段, 双重乳液将完全破裂, 其尾部回缩, 界面曲面半径减小, 尾部界面张力增加, 流阻减小, 上游压力pinlet相应降低(图11(b)中—时间段).

为深入探究“隧道效应”对双重乳液破裂行为的影响, 图13给出了隧道破裂工况中典型几何参数(包括外液滴颈部厚度、外液滴前端运动距离、内液滴颈部厚度和内液滴前端运动距离)在squeezing阶段随时间演化的情况. 与阻塞破裂类似,随时间演化分为两个子过程, 子过程 I 中和呈线性关系, 而在子过程II中两者呈指数关系(见图13(a)插图). 不同之处在于, 子过程 II的指数约为3.75/5, 小于阻塞破裂的4/5. 这表明隧道的出现使得双重乳液外液滴颈部收缩速率降低, 这也间接说明隧道内润滑流动引起的剪切作用并没促进外液滴的挤压形变. 图13(b)给出了外液滴前端运动距离随时间的演化曲线, 如图所示, 在隧道破裂流型中,演化趋势前期为线性变化, 随后变为指数变化, 指数为4/5. 这表明隧道的出现也降低了外液滴沿流向的拉伸速率.

图11 隧道破裂工况乳液前端及尾部界面张力演化情况(Ca = 0.01, Voi = 1.3, l* = 1.3) (a) 乳液前 端 界面张力 ;(b) 乳液尾部界面张力; (c) 乳液前端与尾部界面张力之差Fig. 11. Evolution of the pressure for tunnel breakup (Ca =0.01, Voi = 1.3, l* = 1.3): (a) The Laplace pressure of the forefront droplet interface; (b) the Laplace pressure of the rear droplet interface; (c) the Laplace pressure difference between the forefront and rear droplet interfaces.

图12 隧道破裂工况入口与出口压力演化情况 (Ca =0.01, Voi = 1.3, l* = 1.3)Fig. 12. Evolution of the inlet pressure and outlet pressure of the Y-junction for tunnel breakup (Ca = 0.01, Voi = 1.3,l* = 1.3).

3.3 不破裂

当双重乳液尺寸较小时, Y型分叉微通道中出现不破裂流型, 此时乳液将经历三个阶段, 依次为:entering阶段、sliding阶段、和non-breakup阶段.其中, entering阶段与隧道破裂类似, 后两个阶段是其特有的(见图14和图15). 在entering阶段,双重乳液进入渐扩过渡通道, 由于此时双重乳液长度小于通道宽度, 外液滴界面没有显著形变, 因此上游压力pinlet维持不变. 进入sliding阶段, 当双重乳液前端前进撞击分叉尖角, 双重乳液发生形变, 由圆形被挤压为椭圆形. 在—阶段, 双重乳液内部压力增加, 同时上游压力增加, 在时刻到达最大值. 在—阶段, 由于不对称,双重乳液滑向左侧子通道, 此时双重乳液内存在不对称涡流, 双重乳液内部压力不再增加, 上游压力也开始回落. 值得注意的是, 由于双重乳液滑向单侧子通道, 两个子通道出口压力不再完全相同, 出现略微差异. 当双重乳液完全越过主通道中心线, 双重乳液内部压力以及上游压力也回落到初始水平.

3.4 相 图

图13 隧道破裂工况双重乳液无量纲特征参数在 squeezing 阶段内演化情况 (Ca = 0.01, Voi = 1.3, l* = 1.3) (a) 外液滴颈部厚度 , 插图中给出了 与 的对数坐标图; (b) 外液滴前端运动距离 , 插图中给出了 与 的对数坐标图; (c) 内液滴颈部厚度 , 插图中给出了 与 的对数坐标图; (d) 内液滴前端运动距离 , 插图中给出了 与 的对数坐标图Fig. 13. Evolution of the dimensionless characteristic parameters in the squeezing stage for tunnel breakup (Ca = 0.01, Voi = 1.3, l*= 1.3): (a) The neck thickness of outer droplet , inset is the same data as log( ) versus log( ); (b) the distance travelled by the tip of outer droplet , the same data as log( ) versus log( ); (c) the neck thickness of inner droplet , inset is the same data as log( ) versus log( ); (b) the distance travelled by the tip of inner droplet , the same data as log( ) versus log( ).

图14 不破裂工况 Y 型微通道中压力场、流场及相界面演化 (Ca = 0.01, Voi = 1.3, l* = 0.4)Fig. 14. Evolution of the interface profile, pressure field and flow field during non-breakup in a Y-junction (Ca = 0.01, Voi = 1.3, l*= 0.4).

图15 不破裂工况入口与出口压力演化情况Fig. 15. Evolution of the inlet pressure and outlet pressure of the Y-junction for non-breakup.

图16 乳液在 Y 型微通道中流动相图 (a) 双重乳液;(b) 单乳液Fig. 16. Phase diagrams of droplet behaviours through the Y-junction: (a) Double emulsion; (b) single phase emulsion.

图16 给出了Y型微通道内双重乳液和单乳液运动相图. 图中x轴坐标为主通道中乳液无量纲长度, y轴坐标为主通道毛细数, 所得相图结构与前人研究结果类似. 图中采用幂函数描述隧道破裂和不破裂工况临界线, 其中双重乳液和单乳液运动相图临界线分别为 l*= 0.09Ca–0.49和l*= 0.07Ca–0.49; 采用线性关系描述隧道破裂和阻塞破裂工况临界线, 其中双重乳液和单乳液运动相图临界线分别为 l*= 1.66 和 l*= 1.53. 由图 16 可知, 当保持 Ca 不变, 随着 l*的减小, 流型由隧道破裂逐渐转变为不破裂. 这是由于随着乳液尺寸减小, 乳液与壁面间易出现隧道, 从而会降低乳液形变速率, 最终使得乳液不破裂. 而当保持l*不变, 随着Ca减小, 流型由隧道破裂逐渐转变为不破裂. 这是由于在相同表面张力情况下毛细数越大即外流体流速越大, 从而对应更大的上游压力,而后者正是Y型微通道中乳液产生形变的主要驱动力, 因此毛细数的增加利于乳液破裂发生. 此外,对比两个相图可知, 当单乳液变为双重乳液, 各工况的分界线均向右移动, 不破裂工况与隧道破裂工况分界线位置右移表明双重乳液比单乳液更不易破裂; 隧道破裂与阻塞破裂分界线右移则是由于相同工况下双重乳液上游压力水平高于单乳液工况,从而更容易出现隧道.

4 结 论

本文基于VOF液/液相界面追踪方法, 建立Y型微通道中双重乳液流动破裂非稳态理论模型,数值模拟研究了Y型微通道内双重乳液破裂情况.为揭示双重乳液破裂流型的内在机理, 详细分析了双重乳液流经Y型微通道时的流场信息以及双重乳液形变参数演化, 并给出了Y型微通道内双重乳液运动相图. 研究结果表明:

1) 在Y型微通道中, 双重乳液运动有三种流型, 即阻塞破裂流型、隧道破裂流型和不破裂流型;双重乳液阻塞破裂或隧道破裂过程均可以分为entering, squeezing 和 post-breakup 三个阶段; 而双重乳液不破裂过程经历三个阶段为entering阶段、sliding阶段和non-breakup阶段;

2) 在阻塞破裂和隧道破裂工况的entering阶段以及squeezing阶段, 双重乳液均受上游压力驱动产生形变, 形变过程中乳液前端与尾部界面张力之差阻碍双重乳液形变破裂, 且两者正相关; 与阻塞破裂工况相比, 隧道破裂工况中隧道的出现将减缓双重乳液外液滴颈部收缩速率以及沿流向拉伸速率, 并减缓了内液滴沿流向拉伸速率, 对于内液滴颈部收缩速率影响不大; 当减小毛细数或乳液初始长度, 隧道破裂工况将过渡到不破裂工况, 此时乳液不再破裂;