全文HTML
--> --> -->为深入研究液滴的振动特性, 大多数文章通过实验研究液滴的振动特性. Brunet等[14], Noblin等[15]及Dong等[16]对在倾斜表面、垂直及水平方向上振动的液滴运动特征进行了研究. 周建臣等[17]实验研究了液滴在不同频率下的垂直振动特性, 研究发现液滴对外界驱动的不同响应与接触线的振荡行为、变形程度密切相关. 近些年诸多****通过实验研究自由液滴表面出现重力波的振动频率范围[18-21]. Noblin等[18]研究了接触角和滞后作用对液滴振动模态的影响, 实验通过改变垂直振动的频率和振幅观察到接触线固着-移动两种类型的振动模态. Shin和Lim[19], Kim和Lim[20]以及Park等[21]针对疏水玻璃片上垂直振动液滴内部的流动规律进行了研究, 实验研究得到液滴垂直振动的2, 4, 6, 8四种振动模态, 在模态2和模态4下液滴内部呈现“Y”型流动, 而在模态6和模态8频率下呈现对称的涡流动, 此外通过粒子图像测速法测量(particle image velocimetry, PIV)发现模态4的流速大于模态2, 模态6和8的流速几乎相同. Ramos[22]针对超疏水表面微液滴在不同振幅下的振动特性, 得出超疏水表面微液滴可看作自由液滴的结论. 只有少数文章通过数值模拟研究固体表面上液滴的振动问题[23-25]. James等[25]实现了固体表面上单个液滴在低频和小振幅垂直振动下的数值模拟, 研究发现振动导致自由液滴表面形成重力波, 当振幅高于临界值时会引发液滴破裂.
基于以上实验研究发现, 液滴内部流动具有清晰的三维特征, 并且接触角变化是关键参数之一. 目前关于液滴振动的研究主要是通过实验测定振动频率、振幅和周期等因素对液滴内部流动产生的影响[19-21]. 而在实验研究中, 由于透镜效应, 导致无法准确测量到液滴自由界面附近边缘层的流动特征. 因此, 本文在考虑动态接触角的情况下, 建立了三维液滴振动模型, 实现了液滴在疏水表面上垂直振动的三维数值模拟, 获得了液滴在不同频率下垂直振动的四种模态, 预测了不同共振模态下液滴的内部流场结构, 有助于弥补实验测量时由于存在透镜效应而导致液滴自由界面附近边缘层流场不准确的局限性.
2.1.物理问题及假定
当基底受到强迫振动时, 液滴随基底垂直上下振动, 当驱动加速度超过某一临界值时, 液滴的内表面和外表面由于振动的影响产生压力差, 在压力差的作用下形成表面波. 随着振动的进行, 表面波叠加形成叶瓣, 导致液滴发生形变. 液滴在上下振动时, 其内部的流动受到几何形状的限制, 在重力和黏性力的共同作用下, 形成涡流动或回流. 故在不同频率下的垂直振动, 液滴内部形成不同的流场结构.根据已有的实验研究[19-21]可知液滴在特定的频率下会呈现共振模态. 当液滴处于共振模态时, 表面波叠加形成叶瓣, 叶瓣与叶瓣之间会出现节点, 液滴的振动模态阶数等于节点数. 在共振模态时液滴表面形成的叶瓣是以中心线为轴对称分布, 叶瓣大小和数量相同, 所以液滴内部的微流动是轴对称的. 本文主要研究液滴在垂直振动下的2, 4, 6和8四种共振模态及其内部微流动, 故对三维液滴的振动过程做出如下假设:
1)液滴内部的流动是不可压缩的, 并且在垂直振动下是轴对称的;
2)由于液滴的振动周期远远小于同质量同体积的液滴蒸发的时间, 所以液滴振动过程不考虑自然蒸发;
3)由于本次模拟的液滴半径约为1.14 mm, 小于毛细长度





基于以上假设, 本文采用VOF-CSF模型[26,27]对三维液滴的振动过程进行模拟.
2
2.2.VOF模型
本文模拟的液滴振动涉及气液两相和相界面的形变, 需要得到气液两相的交界面, 因此可以采用VOF模型[26]. 设定空气为主相, 液态水为次相. 气液两相共用一套动量方程, 通过引进相体积分数这一变量, 实现对每一个计算单元相界面的追踪, 气液界面采用几何重建法处理. 模型控制方程如下.连续性方程:









2
2.3.连续表面张力(CSF)模型
动量方程(3)式中的表面张力












2
2.4.接触角模型
(9)式取决于接触角




2
2.5.网格划分与边界条件
如图1所示, 计算区域为正方体空间, 尺寸为 3 mm × 3 mm × 3 mm, 整个计算空间可满足液滴在振动条件下的振动幅度范围和计算要求. 体积为5


Figure1. 3D computational domain (3 mm × 3 mm × 3 mm).
底面采用无滑移边界条件, 即液滴和底面接触处的流体速度等于底面的速度. 计算工质为空气和水, 空气为主相, 液态水为次相, 通过改变振动加速度函数的共振频率, 研究不同共振频率对液滴形态变化和内部流动规律的影响. 在基板上施加一个振动加速度, 使液滴做垂直振动:





Lamb[10]总结出自由微液滴共振模态下的不同共振频率计算公式, 称为Rayleigh共振频率, 表示为









Mode n | Rayleigh方程 计算值 f/Hz | 共振频率 实验值 f/Hz |
2 | 96 | 85 |
4 | 288 | 226 |
6 | 526 | 469 |
8 | 804 | 635 |
表1共振频率的理论值和实验值对比
Table1.Comparisons of theoretical and experimental results for resonance frequency of a 5

3.1.振动模态
图2是体积为5

Figure2. Comparisons of the resonance modes (2, 4, 6, and 8) of a vibrating water drop on a hydrophobic surface: (a) Experimental results[20]; (b) simulation results.
2
3.2.液滴振动模态的动态演化过程
图3是模拟的液滴分别在2, 4, 6和8模态一个周期内的形状变化与实验结果[21]的对比图. 由图3可知, 模拟所得到的液滴形态的变化与实验结果符合较好. 随着共振模态阶数的增加, 液滴的振动周期越来越短, 约为12.4, 4.0, 2.0和1.2 ms. 这一现象与共振频率以及振动幅度大小有关, 模态阶数越高, 共振频率越大, 振动周期越短. 前一个周期的最后相位与后一周期的第一相位重合. 如图3所示, 液滴在模态6和8中形貌变化更丰富, 在模态2和4中则变化的相对简单. 因为在高阶模态中, 共振频率较大, 液滴表面形成的波不断叠加形成的叶瓣数更多, 导致模态6和8中的液滴呈现更丰富的形貌变化.

Figure3. Comparisons between the experimental[21] and simulation results of the droplet shape evolution during one cycle of the vibration at resonance modes 2, 4, 6 and 8.
2
3.3.液滴振动过程的内部流场结构
目前关于液滴垂直振动的研究主要是通过实验测定液滴振动周期、形状、内部流动方式, 由于实验测量时存在透镜效应, 导致无法准确测量液滴自由界面附近边缘层的流动特征. 为了进一步了解液滴垂直振动时的运动规律, 模拟三维液滴在共振模态2, 4, 6和8下的运动情况, 获得液滴内部的流场结构.图4是液滴在模态2, 4, 6和8下内部流场的三维视图, 可以看出在每种模态下, 液滴内部的微流动都呈现轴对称的涡流动. 液滴在模态2和4下, 液滴中心流动向上呈现“Y”型, 周围伴随着涡流动; 而在模态6和8时, 液滴中心流动向上到达液滴顶部, 周围同样伴随着涡流动. 为了进一步了解液滴内部的流场结构, 在图5中对液滴内部的二维剖面进行分析并与实验结果对比. 图5是2, 4, 6和8模态下液滴内部的速度场模拟结果与实验结果[20]对比, 左侧是液滴内部流动的可视化实验结果, 右侧是模拟的液滴内部速度流线图. 随着模态阶数的增加, 多样化的流动模式得到了体现.

Figure4. The three-dimensional flow field inside the droplet at modes 2, 4, 6, and 8.

Figure5. Comparisons between the experimental[20]and simulation results of the internal flow structure of the droplet during the vibration at modes 2, 4, 6 and 8: (a) Mode 2; (b) mode 4; (c) mode 6; (d) mode 2.
图5(a)是液滴在模态2振动时内部的速度场模拟结果. 液滴在振动开始时, 底端的流体受力向上流动, 顶端流体则因为重力影响向下流动, 并在液滴上部形成对称的回流. 随着振动时间的增加, 液滴的主体流动由中间向上, 经过液滴上侧叶瓣叠加的节点处, 再沿着液滴表面回到三相接触线上, 同时在重力和黏性力的共同作用下, 在液滴中心轴两侧形成对称的涡流动. 图5(b)显示了模态4的内部流场, 与模态2流动方式类似, 内部流动皆是自底部中心向上的循环流动, 最后返回液滴的中心. 模态4中的流型更长, 且产生的循环流动位置和大小与模态2不一样, 在液滴上部会形成一个涡环. 故液滴在模态2和模态4中内部流动都呈现“Y”型流动.
图5(c),(d)分别是模态6和模态8的内部流动情况. 如图所示, 与模态2和模态4的流型不同, 以液滴中心线为轴形成两个大尺寸的涡流动, 并且在模态6时, 液滴边缘处出现涡流动, 而在模态8时, 液滴边缘处则没有出现涡流动. 出现这种情况的原因是液滴在模态8时表面形成的叶瓣较多且较小, 几何空间不足以诱导形成回流. 因此, 在模态6和模态8中液滴的内部流动循环变大形成对称的涡流动而不是“Y”型流动.
为了进一步探究液滴的内部流动规律, 对液滴中心底部的垂直速度进行了计算. 在实验测量时, 由于存在透镜效应, 通常是选取某一特定区域测量速度值, 无法准确测量液滴内部完整区域的速度值[20,21]. 为了与实验测量速度值对比, 选取了同样区域计算了液滴的平均垂直速度, 表2为液滴在四种共振模态下周期内的垂直速度平均值的实验值[20]与模拟计算值. 从表2中可以看出, 模拟值与实验值呈现同样的趋势, 在模态4时液滴速度平均值约为模态2速度值的两倍; 而模态6和模态8的速度差值不大. 模拟结果与实验结果数值上存在一定的差距. 从速度大小和量级上来看, 模拟的液滴速度应该小于或接近基板速度. 通过估算基板速度发现模拟得到的液滴局部区域平均速度值和基板速度值接近, 因此理论上模拟结果是准确的. 实验和模拟数据之间的差异可能主要是实验条件和测量误差导致.
模态 | 模拟速度值/(mm·s–1) | 实验速度值/(mm·s–1) |
2 | 1.76 | 0.36 |
4 | 3.37 | 0.79 |
6 | 10.43 | 3.87 |
8 | 11.05 | 3.61 |
表2液滴在模态2, 4, 6和8下中心底部的平均垂直速度的实验[20]与模拟对比
Table2.Comparisons between the experimental[20] and simulation results of averaged vertical velocity at the central bottom region of the droplet in modes 2, 4, 6, and 8.
图6为液滴完整区域在2, 4, 6和8模态下的平均速度随时间的变化曲线图. 从图6中可以看出, 四种模态下的液滴平均速度皆呈现正弦波动变化趋势, 且平均速度的峰值按模态阶数的高低排序, 在模态8时的液滴平均速度峰值达到最大. 在模态2时由于振动频率较低, 液滴内部速度变化趋于平缓, 而其他高频模态下液滴速度波动较大.

Figure6. The variations of the velocity with time at modes 2, 4, 6, and 8.
2
3.4.液滴振动过程中动态接触角的变化
液滴在随基板垂直振动过程中, 由于接触角的变化量不为零, 因此有必要考虑接触角变化对液滴振动形态演变的影响. 如图7所示, 考虑动态接触角模拟得到的模态2液滴振动幅值随时间变化曲线(Num.1)和考虑静态接触角模拟得到的模态2液滴振动幅值曲线(Num.2)具有相似的趋势, 但在峰值上有一定的差异, 曲线(Num.1)和实验测得的幅值变化曲线(Exp.)更为接近. 不考虑动态接触角就会导致模拟的液滴高度偏低, 从而影响准确模拟液滴形态演化. 所以, 模拟中考虑了动态接触角的影响.
Figure7. The variations of droplet vibration amplitude with time at mode 2. Num.1: dynamic wetting angle; Num.2: static contact angle; Exp.: experiment[20].
图8给出了在4种模态下液滴的动态接触角随时间的变化规律. 从图8中可以看到, 在四种不同模态下, 液滴的动态接触角都是以先增大到波峰再减小到波谷这样的趋势变化. 当施加振动加速度后, 液滴的接触线受振动作用, 液滴从接触线附近开始加速运动, 动态接触角开始增大, 此时动态接触角大于静态接触角为前进接触角. 随着振动持续液滴表面形成波瓣导致液滴变形, 液滴边缘的压力逐渐增大, 并且受到黏性阻力和毛细力的阻碍作用, 液滴的运动开始减速, 前进接触角开始减小, 造成接触线停止前进运动, 此时动态接触角接近于静态接触角115°. 由于液滴边缘积累了一定的压力, 在这个压力的作用下, 液滴开始做回复运动, 并造成接触线的加速后退, 动态接触角变小, 此时动态接触角小于静态接触角115°为后退接触角; 随着回复运动的进行, 液滴边缘压力变小, 接触线的后退开始减速, 后退接触角增大直至接近静态接触角. 随着液滴振动的进行, 接触线的前进与后退运动重复之前的变化规律, 动态接触角也随之变化.

Figure8. The variations of the dynamic contact angle with time at modes 2, 4, 6, and 8.
从图8中还可以看到, 在前进接触角阶段, 接触线的运动速度比较大, 因此动态接触角偏离静态接触角比较多. 在本文的计算条件下, 最大接触角偏离静态接触角接近6°. 综合图7和图8结果, 考虑动态接触角对于准确模拟液滴形态演化尤为重要.
2
3.5.液滴振动过程中润湿面积的变化
液滴在振动过程中会经历铺展和收缩两个过程, 图9为振动周期内液滴润湿面积在2, 4, 6和8共振模态下的变化曲线, 润湿面积皆呈现正弦波动变化趋势, 液滴初始润湿面积约为3.4 mm2. 从图9中可以看出, 振动开始时, 液滴的润湿面积逐渐增大, 在大约


Figure9. The variations of the wetting area with time at modes 2, 4, 6, and 8.
1) 液滴在特定的频率下发生共振并且呈现特殊的模态形状, 较高的共振模态阶数有更多的节点叶瓣. 随着模态数的增加, 液滴表面形成的叶瓣数目增加, 尺寸逐渐减小. 在四种共振模态下, 模态阶数越高, 共振频率越大, 周期越短.
2) 在所有模态中, 液滴的内部流动方式呈现轴对称的. 在模态2和模态4时, 内部流动从液滴底部中心开始, 经过液滴上侧叶瓣叠加的节点处, 呈现“Y”型流动方式; 在模态6和模态8中, 液滴内部流动呈现对称的涡流动方式. 且共振模态阶数越高, 液滴内部速度平均值越大.
3) 动态接触角对液滴振动幅值变化有一定的影响, 且动态接触角明显偏离静态接触角, 表明考虑动态接触角对取得正确数值模拟结果的必要性.