表面活性剂如何影响微通道中液滴自发断裂?——从T形通道到不规则通道的数值模拟研究

学术论文 Jan 23, 2026

想象你在厨房里用水龙头往两个方向的岔管里注水——当水流足够快时,水柱会在岔口处一分为二。如果把管道缩小到头发丝的粗细,换成两种互不相溶的液体(比如油和水),再往油滴表面撒上一层洗洁精般的"表面活性剂"(surfactant),液滴还会乖乖分裂吗?分裂的临界条件又是什么?这正是张景刚、刘海湖、王东、崔海航和陈力五位研究者在Journal of Fluid Mechanics上发表的最新成果所回答的核心问题。

微流控中的液滴分裂难题

微流控芯片(microfluidic chip)是当今生物医学和化学分析的"迷你工厂",其核心操作之一就是让液滴在T形微通道(T-junction microchannel)中被动分裂为两个子液滴。前人已经发现,液滴颈部在收缩过程中会经历一个关键转变:从惯性主导阶段过渡到界面张力主导阶段,一旦越过临界点,颈部就会因Rayleigh-Plateau不稳定性而自发坍缩,瞬间断裂。

然而,关于这个临界转变的前提条件,学界存在分歧。有人认为颈部曲率超过其他界面区域即可触发断裂,有人则主张颈部的毛细压力必须超过内部Laplace压力差。更关键的是,表面活性剂如何影响这一临界条件,此前从未被系统研究过。而实际应用中,表面活性剂几乎是标配——它既能降低界面张力、调控液滴大小,又可能因Marangoni效应带来意想不到的后果。

Figure 6: 液滴界面在临界时刻前后的演化过程
图1 液滴界面在临界时刻前后的演化:左图为x-z截面,右图为y-z截面。当流动恰好在临界时刻停止时,颈部迅速坍缩至断裂。

格子Boltzmann方法:从微观粒子到宏观液滴

研究团队采用了他们近期发展的格子Boltzmann方法(Lattice Boltzmann Method, LBM)来模拟这一过程。简单来说,LBM不直接求解Navier-Stokes方程,而是让大量"虚拟粒子"在规则网格上碰撞和迁移,它们的集体行为自然涌现出流体的宏观运动。这种方法特别擅长处理多相流和复杂几何边界,而且天然适合并行计算。

为了研究液滴自发断裂的临界条件,研究者巧妙地设计了一套"停流模拟"(stop-flow simulation)方案:先让液滴在T形通道中正常流动,在某一时刻突然将全场速度归零、出口恢复大气压,然后观察液滴颈部是回缩恢复还是继续坍缩断裂。通过二分法逐步逼近,就能精确锁定临界颈部厚度δcr*。

关键发现一:触发区域——毛细压力突破Laplace压力的"突破口"

对于规则矩形T形通道中的无表面活性剂液滴,研究揭示了一个清晰的自发断裂前提:在液滴颈部的特定区域("触发区域"),局部毛细压力必须超过颈部的Laplace压力差。毛细压力驱动颈部收缩,Laplace压力差则抵抗收缩。当触发区域的毛细压力"翻越"Laplace压力这道防线时,颈部就会不可逆地坍缩。

这一判据比此前提出的"曲率主导"和"压力平衡"判据更加精确和普适。研究还发现,临界颈部厚度δcr*随通道高宽比χ的增大而单调递增,且与解析公式δcr* = χ/(χ+1)吻合良好。

Figure 11: 不同高宽比下临界时刻的液滴形态
图2 不同高宽比(χ)下,液滴在临界时刻的形态。随着χ增大,颈部平均半径增大,临界颈部厚度相应增加。

关键发现二:表面活性剂降低临界厚度,且效应随高宽比增强

加入表面活性剂后,界面张力降低,毛细压力整体减弱。要达到自发断裂所需的毛细压力,就需要更大的界面曲率——也就是更细的颈部。因此,表面活性剂使临界颈部厚度δcr*减小

更有趣的是,虽然不同高宽比下液滴界面吸附的表面活性剂浓度几乎相同,但表面活性剂的效应却随χ的增大而增强。研究者解释道:随着χ增大,临界颈部厚度增大,曲率对毛细压力的贡献减弱,界面张力的贡献相对增强。既然表面活性剂主要通过降低界面张力起作用,它在界面张力"话语权"更大的高χ通道中自然影响更显著。

Figure 13: 有无表面活性剂体系下临界颈部厚度随高宽比的变化
图3 无表面活性剂(黑色虚线)与有表面活性剂(红色实线)体系下,临界颈部厚度δ*cr随高宽比χ的变化。表面活性剂使δ*cr降低,且差异随χ增大而增大。

关键发现三:不规则通道中的线性和指数衰减

研究进一步在T形通道的上壁引入收缩段,构建不规则微通道。收缩深度d*和长度l*是两个关键几何参数。

对于收缩深度的影响,研究发现:在无表面活性剂体系中,δcr*随d*线性递减;而在有表面活性剂体系中,δcr*随d*呈指数衰减。这一差异源于表面活性剂效应随d*增大而逐渐减弱——当曲率贡献因收缩段加深而增强时,界面张力贡献的权重下降,表面活性剂"降界面张力"的作用空间被压缩。

Figure 19: 不同收缩深度下液滴在临界时刻的轮廓
图4 不同收缩深度(d*)下,液滴在临界时刻的x-z截面轮廓。黑色虚线为无表面活性剂,红色实线为有表面活性剂。收缩越深,颈部上方曲率越大。

对于收缩长度的影响,两个体系中δcr*均随l*线性递增,但有表面活性剂体系的斜率更陡,说明表面活性剂效应随l*增大而逐渐减弱。

基于这些规律,研究者提出了四个经验公式,分别预测两种体系下δcr*随d*和l*的变化,拟合优度R²均超过0.84。

Figure 22: 有无表面活性剂体系下临界颈部厚度随收缩深度的对比
图5 无表面活性剂与有表面活性剂体系下,δ*cr随收缩深度d*的对比。两者差异先增后减,在d*≈0.2处达到最大。

展望:从经验公式到统一理论

这项工作首次系统揭示了表面活性剂对T形微通道中液滴自发断裂临界条件的影响机制,提出的"触发区域毛细压力判据"为理解和预测微流控液滴分裂提供了更精确的理论工具。四条经验公式也为工程设计提供了实用的预测手段。

当然,经验公式终究是经验性的。研究者也坦诚指出,当前还缺乏一个能够统一描述复杂几何和表面活性剂条件下颈部自发断裂动力学的理论框架。未来,他们计划推导一个广义物理模型,而本文揭示的经验规律将成为验证该模型不可或缺的基准。

对于微流控芯片的设计者而言,这项研究传递了一个实用信息:表面活性剂的用量并非越多越好——当浓度超过临界胶束浓度(CMC)后,继续添加表面活性剂对临界颈部厚度几乎没有影响。精准控制表面活性剂浓度,才能在降低界面张力和避免浪费之间找到最优平衡。


论文信息

Zhang J, Liu H, Wang D, Cui H, Chen L. Modelling surfactant effects on the onset of spontaneous droplet breakup in both regular and irregular T-junction microchannels. J. Fluid Mech., 2026, 1027: A45.

DOI: 10.1017/jfm.2025.11098

Tags