当液滴撞上"微锥森林"——锥阵列表面如何操控液滴撞击动力学

学术论文 Apr 1, 2025

当液滴撞上"微锥森林"——微结构表面如何操控液滴的命运

一滴水从高处落下,撞击桌面后摊开、回缩、再弹起——这个日常生活中毫不起眼的瞬间,背后却隐藏着丰富的物理图景。液滴撞击固体表面的行为,从飞溅(splashing)到沉积(deposition)再到回弹(rebound),深刻影响着喷墨打印、防冰涂层、热交换器等众多工业应用。那么,如果我们在表面"种"上一片微小的锥形阵列,液滴的命运又将如何改变?

研究背景:微结构表面的"魔法"

在自然界中,荷叶表面的微纳结构让水珠轻松滚落,蝴蝶翅膀的纹理使其保持洁净——这些现象启示科学家:表面微结构能够显著改变液滴与固体表面的相互作用。以往的研究主要关注超疏水表面(接触角大于150°),但对于弱疏水性(接触角约120°)的微结构表面,液滴的动力学行为仍不清晰。更重要的是,实验手段难以窥探液滴与微结构间隙中气体的复杂交互,而传统的数值方法又面临精度或稳定性的瓶颈。

本研究采用改进的颜色梯度格子玻尔兹曼方法(colour-gradient lattice Boltzmann method, LBM),系统模拟了液滴撞击平表面和锥阵列表面的动力学过程。该方法具有三大优势:对微小液滴的溶解极低、严格保证每种流体的质量守恒、在各种粘度比下保持数值稳定性。

液滴撞击锥阵列表面示意图
图1 液滴撞击锥阵列表面的计算模型示意图:(a) 三维计算域,(b) x-z截面,(c) y-z截面,(d) 锥体几何与网格配置。灰色代表固体表面,绿色代表液滴。

雷诺数的双面效应:从"助推器"到"绊脚石"

雷诺数(Reynolds number, Re)表征惯性力与粘滞力的比值。研究发现,锥阵列对液滴回弹的影响随雷诺数的改变而"翻转":

低雷诺数(Re = 10):锥阵列促进液滴回弹。锥体的存在增大了液滴变形程度,使更多动能转化为界面能;同时减小了液滴与底面的接触面积,降低了吸附力。两者协同,让液滴弹得更高。

高雷诺数(Re = 350):锥阵列反而抑制液滴回弹。在高惯性条件下,平表面上的液滴变形持续增大,界面能的积累更充分;而锥阵列"托住"了液滴,限制了其进一步变形,导致动能向界面能的转化效率降低。

不同雷诺数下液滴撞击平表面与锥阵列表面的对比
图2 液滴撞击平表面(上行)与锥阵列表面(下行)的对比:(a) Re=10,锥阵列促进回弹;(b) Re=350,锥阵列抑制回弹。

韦伯数的调控:弹跳高度的"钟形曲线"

韦伯数(Weber number, We)反映惯性力与表面张力的比值。固定Re=350,研究We的影响发现:

对于平表面,最大回弹高度随We先增后减,在We≈20处达到峰值。We较小时,增大We使液滴变形加剧、界面能积累增多,有利于回弹;We过大时,液滴上界面直接触碰底面,吸附力急剧增大,回弹反而减弱。

对于锥阵列表面,回弹高度随We单调递减。锥纤维始终限制液滴变形,而表面张力随We增大而减小,导致存储的界面能持续降低。

不同韦伯数下液滴撞击的瞬时形貌
图3 不同韦伯数下液滴撞击平表面(上行)与锥阵列表面(下行)的瞬时形貌:(a) We=5,(b) We=20,(c) We=35。注意高We时锥阵列表面出现的气泡捕获现象。

相图:液滴命运的"地图"

通过系统性参数扫描,研究绘制了液滴回弹模式的相图(phase diagram)——这就像一张液滴命运的"地图",清晰标出了沉积区与回弹区的边界。

对于平表面,液滴在Re≤150时始终沉积于表面;当Re≥175后,回弹区域随Re增大而显著扩展。对于锥阵列表面,临界韦伯数在Re=10~300范围内保持恒定,但当Re进一步增大时临界值下降——惯性力主导了液滴的脱附行为。特别值得注意的是,只要We超过10,锥阵列就能有效防止液滴粘附于底面。

平表面回弹模式相图
图4a 平表面回弹模式相图:Re≤150时液滴始终沉积
锥阵列表面回弹模式相图
图4b 锥阵列表面回弹模式相图:We>10时即可有效防止液滴粘附

展望:从模拟到设计

这项研究不仅揭示了锥阵列微结构对液滴动力学的双重调控机制,更为纹理表面的优化设计提供了定量指导。当需要液滴快速脱离表面(如防冰、自清洁)时,低Re条件下锥阵列表现优异;而需要液滴滞留表面(如微流控、涂层固化)时,高Re下锥阵列的抑制效应恰好可用。

未来,研究团队计划采用局部网格细化技术以更精确地捕捉锥体几何细节,并探索倾斜锥阵列——让液滴不仅在垂直方向回弹,还能在水平方向定向输运。这将为微流控芯片上的液滴操控开辟新的可能。

————————————————————————

论文信息

Zhang J, Zhao W, Liu H, Wang D, Cui H, Chen L. Numerical study of droplets impacting on flat and cone-arrayed surfaces. International Communications in Heat and Mass Transfer, 2025, 163: 108729.

DOI: 10.1016/j.icheatmasstransfer.2025.108729

Tags