液滴撞击疏水球面的高速可视化实验研究

2023-01-10王程遥

上海电力大学学报 2022年6期

王程遥, 张魏, 彭程

(上海电力大学能源与机械工程学院, 上海 200090)

液滴撞击固体表面的现象非常常见,广泛存在于工业生产中,例如喷墨打印、喷雾冷却、燃烧室燃烧等[1-3]。液滴撞击固体表面受到液滴性质[4]、撞击速度[5]、撞击表面[6]和周边环境[7]等多种影响因素的共同作用,可能出现铺展、部分反弹、完全反弹和飞溅等行为特征[8-10]。液滴撞击壁面是一个复杂的多相流过程,涉及了液滴的瞬态行为、内在动力学机理等关键科学问题,引起了国内外学者的广泛关注。在喷雾冷却、化学涂层等领域,液滴在弯曲表面上的撞击过程通常出现在较小的局部区域,而球形表面是一种典型的、规则的曲面,液滴在球面上的碰撞动力学研究对理解液滴在曲面上的碰撞机理具有重要意义[11]。

近年来,国内外科研工作者对液滴撞击固体表面进行了大量的实验研究和数值模拟研究。CHEN M J等人[12]研究了液滴撞击后的润湿过程,对液滴撞击平面和倾斜表面时的三相接触线和动态接触角进行了观测,分析了液滴尺寸、表面润湿性、倾斜角度等因素对液滴撞击行为的影响。TABBARA H等人[13]利用数值模拟方法,模拟了金属液滴在固态平面上的撞击、铺展和凝固的瞬间过程,研究了增加撞击速度后金属凝固过程的变化情况,分析了金属液滴撞击时可能发生的断裂现象。QIN M X等人[14]通过实验研究了不同液滴撞击受热不锈钢平面后的铺展和弹跳过程,分析了韦伯数和莱顿弗罗斯特温度点对液滴撞击行为的影响。ZHU Y等人[15]用数值模拟方法研究了液滴撞击球体的问题,发现液滴铺展阶段惯性力处于主导地位,并对球面上液膜的厚度进行了预测,且预测结果与实验结果吻合度较高。HAO J G[16]研究了表面粗糙度对液滴撞击产生飞溅的影响,实验结果表明,粗糙度略微增大时会显著增强prompt飞溅和corona飞溅,而粗糙度明显增大时会促进prompt飞溅并抑制corona飞溅,而粗糙度略微增大会增强飞溅的原因是润湿面积随着表面粗糙度的增大而减小。IKALO等人[17-19]系统性地研究了水和甘油液滴以不同角度撞击不同材质固体表面后的行为特征,定性分析了各因素对液滴撞击行为的影响。

目前,大部分的研究仍着眼于液滴撞击平面,而实际过程中撞击面多为带有一定曲率的球面,考虑到实际应用,本文搭建了可视化实验平台,对液滴撞击球面的动力学行为开展实验研究,观测液滴以不同撞击速度(0.46 m/s,0.79 m/s,1.08 m/s)撞击不同直径(8 mm,14 mm,20 mm)的疏水球面,定量表征液滴撞击过程的动态行为,并分析液滴撞击速度和液滴-球体直径比对液滴撞击行为特征的影响。

1 实验系统与方法

液滴撞击球面实验系统由计算机、高速摄像机、光源、注射泵、注射器、可调节高度铁架台和不同直径疏水表面球体组成,具体如图1所示。

图1 实验系统

实验时由注射泵推动注射器产生液滴,液滴自由落体撞击目标球体,高速摄影机拍摄画面后传输到计算机后处理图片,抓取数据。液滴撞击速度由调节掉落高度来控制。不同直径疏水球体是将玻璃球经过表面处理后得到,制备得到的球面静态接触角为105°±3°,属于疏水表面。实验过程中的参数表征如图2所示。

图2 实验各参数

选取液滴下降过程中最临近撞击时刻的连续相邻2张图片,如图2(a)所示。根据帧率计算连续相邻2张图片的时间间隔,记为Δt,液滴中心点在Δt内下降的高度记为Δx,液滴速度v的计算公式为

(1)

取相同工况下多次实验测得速度的平均值,可认为是液滴撞击时的真实速度。

在液滴撞击固体表面的过程中,液滴的铺展因子是表征液滴铺展行为的重要参数,其计算公式为

(2)

(3)

式中:D*——液滴无量纲铺展因子;

d——液滴在疏水球面上的铺展直径,即液滴在球面上所铺展的弧长;

d0——液滴当量直径[20];

DV,DH——液滴竖直直径和水平直径。

DV和DH如图2(b)所示。经过多次测量取平均值,实验中液滴当量直径为2.04 mm。液滴无量纲中心高度h*是另一个表征液滴铺展行为的重要参数,其计算公式为

(4)

式中:h——液滴中心高度。

h如图2(c)所示。另外,无量纲时间定义t*为

(5)

式中:t——液滴撞击时间,液滴与表面接触时t=0。

液滴-球体直径比ε为

(6)

式中:D——撞击球体直径。

2 实验结果与分析

2.1 液滴撞击过程动态行为研究

为研究液滴撞击疏水球面后的动态行为特征,本文通过将液滴以不同速度撞击球体和相同直径液滴撞击不同直径球面开展了实验研究。图3给出了液滴撞击疏水球面的动态演化过程(v=0.79 m/s,ε=0.30)。图4是此工况下液滴撞击疏水球面行为特征的变化趋势。

图3 液滴撞击疏水球面的动态演化过程

图4 液滴撞击疏水球面行为特征的变化趋势

液滴撞击疏水球面后的行为大致可分为铺展、回缩和振荡-稳定3个阶段。铺展阶段从液滴下端接触球面开始,液滴达到最大铺展结束。达到最大铺展同时回缩阶段开始,当液滴铺展因子进行最大幅度缩小后且尚未重新增大时,回缩阶段结束。振荡-稳定阶段液滴持续小幅振荡并最终趋于静止。当t*=0时,液滴与球面发生撞击,以撞击点为圆心向四周铺展,在铺展早期液滴上端依然保持球状,在较短时间内球面出现表面波且呈阶梯状。随着铺展的进行,液滴中心高度逐渐降低,从凸起最终变为凹陷。最终液滴铺展因子达到最大值,D*=1.81,此时液滴中心高度低、边缘高,液滴中心出现空气柱,变为环状液滴,如图3(t*=1.60时)所示。此时,图4中液滴无量纲铺展因子D*最大,液滴无量纲中心高度h*最小。这是由于边缘液滴的动能已经耗尽,而中心液滴在动能的驱使下依然向四周扩散,导致中心液量减少从而凹陷。液滴铺展因子达到最大值后,液滴的铺展阶段结束,回缩阶段开始。在液滴表面能的作用下液滴回流,中心部位重新升高,h*由0.29增大到0.58,从凹陷转为凸起,直到液滴中心高度达到最大值。此时液滴铺展因子在回缩状态下可达到最小值,D*=1.09,液滴回缩阶段结束。之后液滴会经历多次的铺展-回缩的振荡周期,并最终趋于平衡稳定状态。这就是振荡-稳定阶段。在此阶段,液滴铺展因子小于铺展阶段结束时的最大值,大于回缩状态结束时的最小值,波动幅度逐渐减小趋于稳定,回缩时的液滴中心高度小于最大值且随时间逐渐降低。由图4可知,经过多次振荡后,液滴固液气三相接触线不再移动,液滴铺展因子稳定时,D*=1.22,液滴上端振动幅度逐渐减小但振荡时间长于三相接触线,液滴中心高度迟于液滴铺展因子达到稳定,稳定时液滴中心高度为h*=0.40。最终液滴静止在球面上呈椭球状,达到平衡稳定状态。

2.2 撞击速度对液滴撞击行为特征的影响

撞击速度对液滴撞击疏水球面后的行为特征有着明显的影响。图5给出了液滴以不同速度撞击相同直径疏水球面的动态演化过程(ε=0.11)。图6是不同速度下液滴撞击疏水球面行为特征的变化趋势。

图5 液滴以不同速度撞击疏水球面的动态演化过程

图6 液滴以不同速度撞击疏水球面特征的变化趋势

图5中,第1～3张为液滴的铺展阶段,第4～5张为回缩阶段,第6张为稳定状态。由图6(a)可知,液滴中心高度在撞击后经过多次振荡,振荡幅度随撞击时间的增大而减小。由图6(b)可知:液滴铺展因子在液滴撞击球面后迅速增大,达到最大值后减小并最终趋于稳定;液滴铺展因子最大值与速度成正比,撞击速度越大,铺展因子最大值就越大,撞击速度分别为0.46 m/s,0.79 m/s,1.08 m/s时,铺展因子最大值分别为1.43,1.85,2.16;而不同撞击速度的液滴最终达到稳定状态时铺展因子较为接近,D*的波动范围在1.13～1.18,表面液滴动能对液滴的最终状态没有明显影响。对比图6(a)和图6(b)可以发现:固液气三相接触线的稳定要早于上层液滴的稳定;不同撞击速度时,液滴铺展因子达到最大值的无量纲时间较为接近,因此铺展速度随撞击速度的增大而增大,且速度变化对铺展时间影响较小;液滴撞击速度增加时,液滴惯性力增大,带来的动能增大,液滴铺展因子最大值增大,需要克服液滴黏性耗散所需的能量也增大;动能转变而来的表面能也增大,回缩阶段的液滴中心高度随之增加。

2.3 液滴球体直径比对液滴撞击行为的影响

图7为撞击速度v=0.79 m/s时,液滴撞击不同直径疏水球面的动态演化过程。图8为其动态演化过程中液滴撞击不同直径的疏水球面行为特征的变化趋势。

图7 液滴撞击不同直径疏水球面的动态演化过程

图8 液滴撞击不同直径疏水球面行为特征的变化趋势

由图8(a)可知,液滴中心高度经历多次振荡后趋于稳定,撞击速度v相同、液滴-球体直径比ε不同时,振荡周期和幅度都接近。由图8(b)可知:撞击速度v相同、液滴-球体直径比ε不同时,液滴撞击球面铺展直径最大值与液滴-球体直径比ε呈反比,ε增大时,液滴铺展直径最大值小幅减小,ε为0.11,0.17,0.30时,D*分别为1.81,1.85,1.86;随着液滴-球体直径比ε的增大,液滴各阶段持续时间减少,振荡稳定阶段振荡幅度增大;液滴-球体直径比ε不同时,液膜在球面上最终稳定的铺展因子接近,说明液滴撞击球面后,最终稳定附着在球面的润湿面积差异较小。因此,液滴-球体直径比ε增大对液滴撞击过程起到了促进的作用,各阶段持续时间小幅度缩短,除此之外,并无明显影响。