锥形纤维上自驱动液滴的动力学行为:格子玻尔兹曼模拟研究

学术论文 Aug 22, 2023

一滴水的"自助旅行":锥形纤维上的自驱动液滴

你是否注意过,仙人掌的刺尖上,清晨的雾水会神奇地沿着刺身向根部滑动?蜘蛛丝上凝结的水珠,也会沿着丝线从细端滑向粗端。这不是重力使然——在微观世界里,一滴水可以凭借"自身的力量"逆流而上,从纤维的尖端一路攀升到粗端。这种看似违反常识的现象,背后隐藏着流体力学中一个精妙的物理机制。

锥形纤维上自驱动液滴示意图
图1 锥形纤维上液滴自驱动示意图。前移接触角(advancing contact angle, θa)和后退接触角(receding contact angle, θr)的不对称,驱动液滴从锥尖向锥底自发迁移。

挑战:两种观点,谁是对的?

关于液滴为什么能在锥形纤维上自发移动,学术界长期存在两种不同的解释。

第一种是拉普拉斯压力梯度模型(Laplace pressure gradient model)。这个观点认为,液滴内部存在从锥尖到锥底的压力梯度,就像一条倾斜管道中水流总是从高处流向低处一样,压力梯度"推动"液滴向锥底迁移。然而,这个模型的推导过程缺乏严格的数学证明——它仅仅是把圆柱纤维上的压力公式,简单地"替换"了几个变量后推广到锥形纤维上。

第二种是移动接触线模型(moving contact line model)。这个模型着眼于液滴前后两端与纤维表面的"接触线"。在锥形纤维上,纤维半径沿轴向逐渐增大,打破了液滴前后形状的对称性:液滴前端的表观接触角大于平衡接触角(equilibrium contact angle, θeq),驱动前接触线向锥底移动;而液滴后端的表观接触角小于平衡接触角,同样驱动后接触线向锥底移动。两者"合力"之下,液滴整体向锥底迁移。

这两种观点哪个更准确?对于大锥角、大变形液滴的情况,之前的数值和理论研究几乎没有涉及。这正是本文要回答的核心问题。

方法:格子玻尔兹曼——在格子棋盘上模拟流体

作者采用了改进的三维格子玻尔兹曼方法(lattice Boltzmann method, LBM)中的颜色梯度模型(color-gradient model)来模拟这一过程。与传统的流体计算方法不同,LBM不在网格上直接求解纳维-斯托克斯方程(Navier-Stokes equations),而是在一个个"格子节点"上,让代表流体的"粒子"进行碰撞和迁移,通过统计大量粒子的集体行为来还原宏观流场。这种方式天然适合捕捉复杂的界面运动和多相流动。

这项工作的关键创新在于:将润湿边界条件(wetting boundary condition)融入了改进的颜色梯度模型中,使其能够处理任意复杂几何表面上移动的接触线,同时还能正确模拟密度比高达100的两相流动。此前的方法要么只能处理静态接触角,要么无法应对曲面上的接触线移动。

液滴在锥形纤维上自发迁移的时间序列
图6 液滴在锥形纤维上自发迁移的仿真时间序列(0-500 μs)。上排为Chen等人(2022)的结果,下排为本文改进LBM方法的模拟结果,两者高度吻合。

发现一:速度先升后降,还有小波动

无论改变哪个参数——邦德数(Bond number, Bo,表征重力与毛细力的相对强弱)、表面润湿性(contact angle, θ)、锥形半角(cone hemi-angle, α)还是液滴体积——液滴的攀升速度都呈现出同一个规律:先加速、后减速,并且在减速过程中出现速度波动。

自驱动液滴在不同时刻的形态变化
图11 Bo = 0.0225条件下自驱动液滴在不同时刻的形态变化。展示了液滴从启动到平衡的完整动力学过程。

为什么会波动?答案是:前移和后退接触角并非同时达到平衡值。想象一个跑步者——他的左右脚步幅不同步时,跑动节奏就会出现抖动。液滴也是如此:前接触角先趋向平衡,后接触角还在追赶,这种"步调不一致"导致了速度的微小振荡。

发现二:重力越强,走得越近

当邦德数增大(即重力相对于毛细力的作用增强)时,液滴攀升的最大高度和润湿面积都单调减小。这不难理解——重力像一根拉住液滴的绳索,越重,液滴走得越近。

不同邦德数下液滴的平衡形态
图12 不同邦德数(Bo = 0.135, 0.09, 0.045, 0.009)下液滴的平衡形态。邦德数越小(毛细力越强),液滴向锥底迁移得越远。

发现三:润湿越强,走得越远

表面润湿性的影响恰好相反:接触角越小(即表面越亲液),液滴攀升越高,润湿面积也越大。毛细力可以简化为 σl·cosθ——与接触角的余弦值成正比。θ越小,cosθ越大,驱动力越强。这暗示了一个实用技巧:少量表面活性剂(surfactant)可以显著降低接触角,从而增强液滴在锥形纤维上的输运能力。

发现四:存在最优锥角

最引人注目的发现是:液滴攀升高度随锥形半角的增大先升后降,在 α = 2.5° 时达到最大值

最优锥角时的液滴平衡形态
图17b 锥形半角 α = 2.5° 时液滴的平衡形态。这是液滴攀升高度最大的最优锥角。

为什么不是越尖越好?当锥角极小时,纤维接近圆柱形,前后对称性几乎消失,驱动力微弱;而当锥角增大时,液滴与纤维表面的接触面积增加,黏性阻力也随之增大。2.5°恰好是驱动力增长和阻力增长之间"甜蜜的平衡点"。

这个最优锥角的发现对实际应用意义重大——在设计仿生集水结构或微流控芯片时,选择恰当的锥角能让液滴自动传输最远的距离,实现最高效率的无泵输运。

展望:从锥形纤维到阶梯状结构

这项研究不仅澄清了自驱动液滴的物理机制——移动接触线模型在大锥角和大变形条件下依然适用——还系统地揭示了各参数对液滴动力学的影响规律。未来,作者计划将研究拓展到锥形结构阵列阶梯状锥形纤维(ladderlike conical fiber),这些设计有望在质量输运、快速冷却、被动集水和微化学合成等领域发挥重要作用。

回到开头那个日常场景:仙人掌刺尖上的水珠,正是凭借毛细力的不对称驱动,一步步从刺尖滑向刺根,最终汇聚到仙人掌的表面,滋养这株沙漠中的生命。理解这个微观世界的"自助旅行",我们就能更好地模仿自然、服务人类。


论文信息

Zhang J, Shen H, Cui H, Chen L, Chen L. The dynamic behavior of a self-propelled droplet on a conical fiber: A lattice Boltzmann study. Phys. Fluids, 2023, 35: 082119.

DOI: 10.1063/5.0164908

Tags