液滴撞击加热壁面雾化弹起模式及驻留时间

荣松; 沈世全; 王天友; 车志钊

doi:10.7498/aps.68.20190097

摘要

Abstract

作者及机构信息

Authors and contacts

文章全文

1. 引　言
2. 实验方法
3. 实验结果与讨论
3.1 液滴撞击形貌及驻留时间
3.2 液滴驻留时间减少机理及两种弹起模式的过渡
4. 结　论

参考文献

施引文献

摘要

液滴撞击壁面是一种常见的现象, 减少撞击过程中液滴在壁面的驻留时间在很多应用中非常重要. 本文在对液滴撞击加热壁面的动力学过程的研究中, 发现了一种能够大幅减少液滴驻留时间的弹起模式, 即雾化弹起模式. 与液滴撞击壁面过程中的回缩弹起模式相比, 雾化弹起模式的液滴驻留时间能减少约40%. 与其他试图减少液滴驻留时间的手段相比, 本文中驻留时间的减少仅需加热壁面, 因此该方法容易实施而且非常可靠. 最后, 本文对雾化弹起模式下液滴驻留时间减少的机理进行了讨论, 提出了一个针对膜态沸腾区域雾化弹起模式发生条件的简化理论模型, 并得出了膜态沸腾区域液滴从雾化弹起模式到回缩弹起模式过渡边界的幂次关系, 该幂次关系与实验数据吻合良好. 此外, 该理论模型也能解释过渡沸腾区域两种弹起模式的过渡边界的变化规律.

关键词:

Abstract

The impact of droplets on surfaces is a ubiquitous phenomenon, and reducing the droplet residence time is the aim of many studies because of the potential applications in self-cleaning, anti-icing, corrosion resistance, etc. This study identifies a mode of droplet bouncing (bouncing-with-spray) that can reduce the residence time significantly. And compared with the way of using complex microstructures on the substrate employed in previous studies, simply heating the substrate to reduce the residence time is novel and simple. The dimensionless residence time decreases down to about 40% compared with that from the traditional retraction-bouncing mode. The reduction in the residence time is due to the burst of vapor bubbles in the liquid film, which results in holes forming in the liquid film and consequently the liquid film recoiling from the holes. The reduction in the recoiling distance leads to the reduction in the recoiling time. Then a simplified theoretical model with considering the energy balance and the critical condition of the bubble burst is proposed. According to this theoretical model, a scaling law is proposed for the transition boundary between the retraction-bouncing mode and the bouncing-with-spray mode in the film boiling regime, and it accords well with our experimental data. This model can also explain the transition boundary between these two modes in the transition boiling regime.

Keywords:

PACS:

47.55.D-	(Drops and bubbles)
47.55.dp	(Cavitation and boiling)
47.55.Ca	(Gas/liquid flows)
47.55.N-	(Interfacial flows)

作者及机构信息

天津大学, 内燃机燃烧学国家重点实验室, 天津　300072

通信作者: 车志钊, chezhizhao@tju.edu.cn

基金项目: 国家自然科学基金(批准号: 51676137)、天津市自然科学基金(批准号: 16JCYBJC41100)和国家杰出青年科学基金(批准号: 51525603)资助的课题.

Authors and contacts

State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China

Corresponding author: Che Zhi-Zhao, chezhizhao@tju.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant No. 51676137), the Natural Science Foundation of Tianjin City, China (Grant No. 16JCYBJC41100), and the National Science Fund for Distinguished Young Scholars, China (Grant No. 51525603).

文章全文

1. 引　言

液滴撞击壁面现象广泛见于自然界和工业应用中, 包括但不限于喷雾冷却、喷墨打印及内燃机缸内燃油喷雾撞壁^[1-3]. 液滴撞击不同温度的壁面的过程会极大地受到传热和相变的影响. 当壁面温度足够高时, 液滴和壁面之间会形成一层蒸汽膜, 即Leidenfrost效应^[4]. 一般来说, 对于不同的壁面温度, 液滴沸腾会经历和池沸腾一致的沸腾形式, 即核态沸腾、过渡沸腾和膜态沸腾^[5-9]. Tran等^[10]以液滴底部是否与壁面接触(即是否存在子液滴喷射的现象)作为判别液滴处于膜态沸腾和其他沸腾区域的依据. 很多文献也依照这种方法判别液滴是否处于膜态沸腾区域, 并进一步研究了液滴的动态Leidenfrost温度会受到表面结构^[11,12]以及撞击韦伯数^[7]等因素的影响. Shirota等^[13]和van Limbeek等^[14]则通过全内反射技术, 对3个沸腾区域的液滴底部与壁面的接触情况进行了直接测量. 他们观察到在核态沸腾区域, 液滴底部一直完全与壁面接触; 在过渡沸腾区域, 由于液滴底部形成的蒸汽膜不稳定, 液滴底部中央区域会和壁面部分接触; 在膜态沸腾区域, 液滴底部不会接触壁面.

在液滴撞击壁面过程中, 液滴与壁面的接触时间(contact time)是一个重要的参数. 当一个无黏液滴撞击在超疏水表面上时, 其与壁面的接触时间受Rayleigh时间尺度的限制^[15]. 由于液滴接触时间会影响到液滴与壁面之间的物质、动量及能量的交换, 在很多应用中, 减少液滴接触时间是很重要的, 例如防结冰、自清洁、防腐蚀等. 因此, 大量研究者采用了多种方法去克服液滴接触时间的理论极限^[15], 例如, 使用带有微脊的表面^[16]、使用带有亚毫米级柱子的超疏水表面^[17]、弯曲表面上的不对称撞击^[18]及使用鸡蛋形状的液滴^[19]等.

在液滴撞击加热壁面过程中, 由于在膜态沸腾区域液滴看似与壁面接触, 而事实上, 液滴和壁面之间会形成一层蒸汽膜, 因此采用驻留时间(residence time)来表征液滴在撞击壁面过程中的停留时间而没有采用接触时间. 本文发现了一种由蒸汽泡破碎而引起的液滴弹起模式, 即雾化弹起模式, 该模式能极大地减少液滴的驻留时间. 与其他减少液滴驻留时间的手段相比, 本文中驻留时间的减少仅需加热壁面, 避免了复杂的表面微结构加工或者表面改性的工艺、超疏水表面性能老化等问题, 因此该方法容易实施而且非常可靠. 本文不仅给出了雾化弹起模式下液滴驻留时间减少的程度, 而且对驻留时间减少的机理进行了讨论, 提出了一个针对膜态沸腾区域雾化弹起模式发生条件的简化理论模型, 并得出了膜态沸腾区域液滴从雾化弹起模式到回缩弹起模式过渡边界的幂次关系. 此外, 该理论模型也能解释过渡沸腾区域两种弹起模式的过渡边界的变化规律.

2. 实验方法

实验装置示意图见图1. 注射器中的液体在注射泵(Harvard Apparatus, Pump 11 elite Pico plus)的推力下以约4 μL·s^–1流速经过一根注射管后到达平口针头顶端, 当液滴所受重力大于表面张力时, 液滴从针头顶端脱离. 随后液滴撞击到一片被铜板加热的抛光硅表面上(表面粗糙度小于0.5 nm). 实验中将两根K型热电偶分别置于距撞击点1 cm的两侧, 并通过PID控制器实现对壁面温度的控制(精度 $\pm \;1\;{\rm ^\circ C}$ ). 实验通过高速相机(Photron Fastcam SA1.1)捕捉水平视角和倾斜视角下的液滴撞壁过程. 背景光源为一个发光二极管(LED)灯. 在进行多次重复实验的条件下, 通过一个Matlab程序进行图像处理, 获得液滴的尺寸(水和甘油的混合物液滴直径为 ${D_0} =$ 2.8 mm, 水滴为2.2—3.4 mm, 乙醇液滴为1.6 mm, 标准差变化范围为0.020—0.022 mm)、液滴的撞击速度(V = 0.9—2.1 m·s^–1, 标准差变化范围为0.006—0.009 m·s^–1)及撞击过程中的液滴轮廓变化情况. 其中液滴的撞击速度, 通过在图像处理中计算撞击前两帧图片中液滴质心位移和时间间隔的比值得到. 本文通过改变甘油和水的混合物中甘油的质量分数来实现液滴黏度的改变, 实验中采用液体的物性参数见表1. 韦伯数用来定义液滴的动能和表面能的比值, 即

$图 1 液滴撞击加热壁面的实验装置示意图　(a)水平视角; (b)倾斜视角\r\nFig. 1. Schematic diagram of the experimental setup for droplet impact on a heated surface: (a) Side-view; (b) aerial-view.$

图 1 液滴撞击加热壁面的实验装置示意图　(a)水平视角; (b)倾斜视角

Fig. 1. Schematic diagram of the experimental setup for droplet impact on a heated surface: (a) Side-view; (b) aerial-view.

下载: 全尺寸图片幻灯片

表 1 本实验中不同液体物性参数表

Table 1. Properties of the liquids used in this study.

液体种类	密度 $\rho {\rm /kg} \cdot {{\rm m}^{{\rm - 3}}}$	表面张力 $\sigma /{\rm mN} \cdot {{\rm m}^{ - 1}}$	动力黏度 $\mu /{\rm mPa} \cdot {\rm s}$	汽化潜热 $\ell /{\rm kJ} \cdot {\rm k}{{\rm g}^{ - 1}}$
水	998	72	1.00	2257
水 + 甘油(28 wt%)	1068	70	2.35
水 + 甘油(34 wt%)	1083	70	3.00
水 + 甘油(39 wt%)	1097	70	3.63
水 + 甘油(56 wt%)	1143	68	9.00
水 + 甘油(75 wt%)	1195	66	42.47
乙醇	790	22	1.07	853

下载: 导出CSV

| 显示表格

$We = \frac{{\rho {D_0}{V^2}}}{\sigma },$

(1)

其中, ρ为液体密度, σ为液体表面张力. 实验中不同种类液体的韦伯数变化范围为10—180.

3. 实验结果与讨论

3.1 液滴撞击形貌及驻留时间

本文首先讨论不同沸腾区域的液滴撞击形貌. 通过全内反射技术^[13,14]可以观察到在核态沸腾和过渡沸腾区域, 液滴底部和壁面的接触情况存在明显的区别: 在核态沸腾区域, 液滴底部会一直附着在壁面上; 在过渡沸腾区域, 液滴底部的中央区域会和壁面存在部分接触. 接触情况的不同会影响子液滴的喷射和液滴在铺展后的动力学行为: 对于核态沸腾区域, 在整个液膜上会产生大量的子液滴喷射, 并且液滴在铺展后, 底部仍然会附着在壁面上; 对于过渡沸腾区域, 子液滴喷射会集中在液膜中央部分, 并且液滴在铺展后, 会发生回缩并弹起离开壁面. 由于采用的表面为光学不透明的硅表面, 因此, 本文将子液滴喷射的不同现象及液滴铺展后是否弹起, 作为区分三种沸腾区域的依据.

对于三种不同的沸腾区域, 通过改变壁面温度及液滴的撞击韦伯数可以获得不同的液滴撞击形貌, 而不同的液滴撞击形貌对于液滴在壁面上的驻留时间有很大的影响. 一般来讲, 液滴在壁面上的驻留时间定义为视觉上从液滴与壁面接触到液滴离开壁面时的时间间隔.

在核态沸腾区域, 由于液滴底部会和壁面直接接触, 当液滴撞击壁面时, 液滴底部的液体会迅速汽化从而形成蒸汽泡^[8], 随后蒸汽泡会在浮力的作用下上升至液滴的自由表面破碎, 并产生大量的子液滴喷射. 如图2(a)所示, 可以看到在液滴铺展过程中存在大量子液滴喷射的现象, 这种现象可以解释为: 当上升到液滴自由表面的蒸汽泡破碎后, 破碎点处会在表面波的作用下产生射流喷射, 进一步喷射的射流会断裂形成子液滴. 在液滴铺展完成后(例如图2(a)中5.1 ms时), 液滴的底部液体仍然附着在壁面上, 直到蒸发完毕.

$图 2 在3个沸腾区域水滴撞击加热壁面的形貌图, 子液滴的喷射如红色圆圈及箭头所示　(a)核态沸腾区域, $We = 37$, $T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (b)过渡沸腾区域回缩弹起模式, $We = 37$, $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (c)过渡沸腾区域雾化弹起模式, 其中(c1)为$We = 69$, $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$, (c2)为$We = 105$, $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (d)膜态沸腾区域回缩弹起模式, 其中(d1)为$We = 25$, $T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$, (d2)为$We = 25$, $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (e)膜态沸腾区域雾化弹起模式, 其中(e1)为$We = 60$, $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$, (e2)为$We = 129$, $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$\r\nFig. 2. Morphologies of water droplet impacting on heated surfaces in three boiling regimes, the red circles and arrows highlight the spray-like ejection of tiny droplets: (a) Nucleate boiling, $We = 37$ and $T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (b) retraction-bouncing mode in transition boiling regime, $We = 37$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$; (c) bouncing-with-spray mode in transition boiling regime, $We = 69$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (c1), and $We = 105$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (c2); (d) retraction-bouncing mode in film boiling regime, $We = 25$ and $T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (d1), and $We = 25$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (d2); (e) bouncing-with-spray mode in film boiling regime, $We = 60$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (e1), and $We = 129$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (e2).$

图 2 在3个沸腾区域水滴撞击加热壁面的形貌图, 子液滴的喷射如红色圆圈及箭头所示　(a)核态沸腾区域,

$We = 37$

$T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (b)过渡沸腾区域回缩弹起模式,

$We = 37$

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (c)过渡沸腾区域雾化弹起模式, 其中(c1)为

$We = 69$

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

, (c2)为

$We = 105$

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (d)膜态沸腾区域回缩弹起模式, 其中(d1)为

$We = 25$

$T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

, (d2)为

$We = 25$

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (e)膜态沸腾区域雾化弹起模式, 其中(e1)为

$We = 60$

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

, (e2)为

$We = 129$

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

Fig. 2. Morphologies of water droplet impacting on heated surfaces in three boiling regimes, the red circles and arrows highlight the spray-like ejection of tiny droplets: (a) Nucleate boiling,

$We = 37$

and

$T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (b) retraction-bouncing mode in transition boiling regime,

$We = 37$

and

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

; (c) bouncing-with-spray mode in transition boiling regime,

$We = 69$

and

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (c1), and

$We = 105$

and

$T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (c2); (d) retraction-bouncing mode in film boiling regime,

$We = 25$

and

$T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (d1), and

$We = 25$

and

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (d2); (e) bouncing-with-spray mode in film boiling regime,

$We = 60$

and

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (e1), and

$We = 129$

and

$T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

for (e2).

下载: 全尺寸图片幻灯片

在过渡沸腾区域, 液滴底部和壁面之间会形成一层不稳定的蒸汽层, 液滴中央区域与壁面之间存在局部接触. 在撞击速度较小时, 例如在 $We = 37$ , $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}{\kern 1pt}$ 时, 如图2(b)所示, 可以看到液滴在铺展后会发生回缩, 随后在撞击时刻12.0 ms时弹起离开壁面, 并且这个过程在液膜中央区域会伴随着子液滴喷射的现象(红色圆圈所示). 这种液滴在达到最大铺展后回缩弹起离开壁面的弹起模式被称为回缩弹起模式. 当撞击速度较大时, 例如在 $We = 69$ , $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}{\kern 1pt}$ 时, 如图2(c1)所示, 可以看到液滴会在没有明显回缩的情况下, 在6.5 ms时弹起离开壁面, 并且同样在液膜中央区域伴随着子液滴喷射的现象(红色箭头所示). 我们称这种弹起模式为雾化弹起, 并且这种弹起模式下的液滴驻留时间要远远小于回缩弹起模式下的驻留时间(6.5 ms vs. 12.0 ms). 当撞击速度进一步增大时, 例如在 $We = 105$ , $T = 380\; ^\circ {\rm C}$ 时, 如图2(c2)所示, 可以看到液滴外围指迹结构的形成以及由指迹结构破碎而产生的小液滴, 这种现象可以解释为液滴在较大的惯性下发生溅射^[20], 但是液滴的中央部分仍然在没有明显回缩的情况下弹起离开壁面, 其驻留时间也比回缩弹起模式下的驻留时间小很多(6.3 ms vs. 12 ms). 因此该模式也被称为雾化弹起模式.

在膜态沸腾区域, 液滴底部和壁面之间会形成一层稳定的蒸汽膜, 即Leidenfrost效应^[4], 并且壁面传递给液滴的热流密度会随之减少. 一般来说, 没有子液滴喷射现象产生时的壁面温度定义为液滴的动态Leidenfrost温度^[7,10-12]. 在这个沸腾区域同样可以观察到上述的两种弹起模式. 对于回缩弹起模式, 相比图2(d2), 图2(d1)所呈现的现象发生于一个相对较低的壁面温度( $400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ vs. $440\; ^\circ {\rm C}$ ), 但有同样的撞击韦伯数( $We = 25$ ). 可以看到在图2(d2)中, 子液滴喷射的现象再次发生, 并且与过渡沸腾区域相比, 子液滴喷射的现象会在整个液膜上发生, 而与核态沸腾相比, 子液滴的粒径更小且分布更均匀. 这可以解释为: 在膜态沸腾区域虽然液滴底部存在一层完整的蒸汽膜, 但是随着壁面温度的升高, 液滴内部的蒸汽泡增长速率也随之增大, 从而使得蒸汽泡脱离蒸汽膜然后上升到自由表面破碎这一过程变得更容易^[10]. 当撞击速度进一步增加时, 例如在 $We = 60$ , $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}{\kern 1pt}$ 时, 如图2(e1)所示, 可以看到雾化弹起模式. 当 $We \!=\! 129$ , $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}{\kern 1pt}$ 时, 则可以看到液滴溅射现象及雾化弹起的弹起模式.

需要指出的是, Tran等^[10]对水滴撞击加热壁面的沸腾区域的分类为: 当壁面温度较低时, 依据是否存在子液滴喷射的现象, 将其分为接触沸腾和膜态沸腾; 当壁面温度较高时, 在膜态沸腾区域由于水滴会再次发生子液滴喷射的现象, 故进一步将膜态沸腾区域分为了稳定膜态沸腾(gentle film boiling)和雾化膜态沸腾(spray film boiling). 本文为了更清楚地表征上述的两种弹起模式及驻留时间的变化, 通过考虑液滴底部与壁面的接触情形对子液滴喷射现象及液滴铺展后是否弹起的影响, 将文献[10]中的接触沸腾区域分类为核态沸腾区域和过渡沸腾区域, 但并没有对膜态沸腾区域进行文献[10]中的分类, 而是依据两种弹起模式对过渡沸腾区域和膜态沸腾区域进行了分类.

上述造成液滴驻留时间减少的雾化弹起模式在本研究采用的所有液体中都可以观察到. 为了确定雾化弹起模式下液滴驻留时间的影响参数, 实验采用了不同比例的甘油和水的混合物来改变液滴的黏度, 并固定液滴的直径 ${D_0}$ 和撞击速度 $V$ , 本文发现该弹起模式下的液滴驻留时间不受液滴黏度的影响, 如图3所示.

$图 3 液体黏度对雾化弹起模式下驻留时间的影响(液滴直径为2.8 mm, 撞击速度为$1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$)\r\nFig. 3. Effect of liquid viscosity on the residence time of the bouncing-with-spray mode. The droplet diameter is 2.8 mm, and the impact velocity is $1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$.$

图 3 液体黏度对雾化弹起模式下驻留时间的影响(液滴直径为2.8 mm, 撞击速度为

$1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$

)

Fig. 3. Effect of liquid viscosity on the residence time of the bouncing-with-spray mode. The droplet diameter is 2.8 mm, and the impact velocity is

$1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$

下载: 全尺寸图片幻灯片

液滴撞击超疏水表面或者在膜态沸腾区域液滴撞击加热壁面的回缩弹起模式的驻留时间受液滴尺寸的影响, 且与液滴自由振荡时间成比例^[21], 即

$\tau = \sqrt {\frac{{\rho R_0^3}}{\sigma }} .$

(2)

通过利用该时间尺度对液滴的驻留时间进行无量纲化, 研究者们得出的无量纲驻留时间为2.2—3.2^[12,22,23]. 在本实验中, 回缩弹起模式下的无量纲驻留时间为2.3, 为了更清楚地表示液滴驻留时间的减少程度, 将雾化弹起模式下的驻留时间进行相同的无量纲化处理, 得出其无量纲驻留时间为1.4, 如图4所示. 因此, 相比于回缩弹起模式, 雾化弹起模式下液滴的驻留时间大约减少40%.

$图 4 回缩弹起模式(空心圆)和雾化弹起模式(实心圆)下的无量纲驻留时间比较\r\nFig. 4. Comparison of the dimensionless residence time in the retraction-bouncing mode (hollow circles) and in the bouncing-with-spray mode (solid dots).$

图 4 回缩弹起模式(空心圆)和雾化弹起模式(实心圆)下的无量纲驻留时间比较

Fig. 4. Comparison of the dimensionless residence time in the retraction-bouncing mode (hollow circles) and in the bouncing-with-spray mode (solid dots).

下载: 全尺寸图片幻灯片

3.2 液滴驻留时间减少机理及两种弹起模式的过渡

图5为倾斜视角下液滴雾化弹起模式的图片序列, 可以看到, 在液滴撞击过程中, 液膜上会产生一个接一个的孔洞, 并伴随着子液滴喷射的现象. 这些孔洞在表面张力的作用下迅速扩大, 整个液膜会从这些孔洞回缩成几个部分(如图5中红色箭头所示), 随后弹起离开壁面. 由于液滴的回缩距离减少, 因此回缩时间也相应减少, 进而液滴在壁面上的驻留时间得到减少. 这种液滴驻留时间减少的机理与液滴撞击带有微脊的超疏水表面的驻留时间的减少机理^[16]相似, 即壁面上的微脊会造成液滴在回缩过程中破碎成几个部分, 从而减少了液滴的回缩距离, 最终导致了驻留时间的减少. 对于雾化弹起模式, 液膜的破碎是由于液膜中蒸汽泡破碎造成的, 所以这个过程会伴随着子液滴喷射的现象. 实验中在更大的放大倍率下, 能够看到在液膜出现孔洞时, 即在撞击时刻为4.07 ms时存在子液滴喷射现象, 如图6箭头所示. 另外, 文献[16]中通过改变液滴撞击在壁面上微脊的位置, 可以实现对液膜上破碎地点的控制, 而本文由蒸汽泡破碎所产生的孔洞位置和数目具有随机性, 难以通过一个简化的理论模型对液滴驻留时间的减少程度进行直接预测. 下面主要针对膜态沸腾区域驻留时间减少的雾化弹起模式的发生条件进行分析.

$图 5 倾斜视角下雾化弹起模式的序列图片($We = 82$, $T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$)\r\nFig. 5. Image series of the bouncing-with-spray mode in aerial view ($We = 82$ and $T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$).$

图 5 倾斜视角下雾化弹起模式的序列图片(

$We = 82$

$T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

)

Fig. 5. Image series of the bouncing-with-spray mode in aerial view (

$We = 82$

and

$T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

下载: 全尺寸图片幻灯片

$图 6 更大的放大倍率下液膜孔洞形成的图片序列(We =122, $T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$)\r\nFig. 6. Image series of hole formation in the liquid film using greater magnification ($We = 122$ and $T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$).$

图 6 更大的放大倍率下液膜孔洞形成的图片序列(We =122,

$T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

)

Fig. 6. Image series of hole formation in the liquid film using greater magnification (

$We = 122$

and

$T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

下载: 全尺寸图片幻灯片

孔洞形成的临界条件即为回缩弹起模式到雾化弹起模式下的过渡. 在液滴撞击壁面过程中, 随着液滴在壁面的铺展, 蒸汽泡的尺寸随之增大, 而液膜随之变薄. 对于两种弹起模式的过渡, 可以认为当蒸汽泡的半径和液膜在最大铺展时的最小厚度相当时, 由于蒸汽泡的破碎, 会在液膜上留下孔洞, 如图7(a)所示. 当蒸汽泡的半径小于液膜的最小厚度时, 液膜上不会产生孔洞; 而当蒸汽泡的半径大于液膜的厚度时, 蒸汽泡则会在液滴达到最大铺展前破碎, 并在液膜上留下孔洞.

$图 7 (a)蒸汽泡破碎理论模型的示意图; (b)采用液滴自由振荡周期无量纲化后的液滴最大铺展时间; (c)对数坐标下膜态沸腾区域水滴的无量纲最大铺展直径和撞击韦伯数的关系图, ${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$; (d)对数坐标下水滴在膜态沸腾区域由回缩弹起到雾化弹起模式过渡的$\Delta T \sim {V^{ - 1}}$边界线, 液滴直径为2.2 mm; (e)水滴撞击加热硅片的相图及从回缩弹起到雾化弹起模式下的过渡, 图中红色虚线为示意, 液滴直径为2.2 mm, “过渡”为过渡腾区域, “膜态”指膜态沸腾区域\r\nFig. 7. (a) Schematic diagram for the theoretical model for the burst of vapor bubbles; (b) the time of maximum spreading diameter of water droplets normalized by droplet free oscillation period; (c) log-log plot for the normalized maximal spreading diameter of water droplets versus We in the film boiling regime, showing ${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$; (d) log-log plot for $\Delta T \sim {V^{ - 1}}$ boundary line between the retraction-bouncing mode and the bouncing-with-spray mode in the film boiling regime for water droplets, the droplet diameter is 2.2 mm; (e) regime map for a water droplet impacting on a heated silicon wafer and the transition between the retraction-bouncing mode and the bouncing-with-spray mode, the red dashed line is only for eye guidance, the droplet diameter is 2.2 mm.$

图 7 (a)蒸汽泡破碎理论模型的示意图; (b)采用液滴自由振荡周期无量纲化后的液滴最大铺展时间; (c)对数坐标下膜态沸腾区域水滴的无量纲最大铺展直径和撞击韦伯数的关系图,

${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$

; (d)对数坐标下水滴在膜态沸腾区域由回缩弹起到雾化弹起模式过渡的

$\Delta T \sim {V^{ - 1}}$

边界线, 液滴直径为2.2 mm; (e)水滴撞击加热硅片的相图及从回缩弹起到雾化弹起模式下的过渡, 图中红色虚线为示意, 液滴直径为2.2 mm, “过渡”为过渡腾区域, “膜态”指膜态沸腾区域

${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$

; (d) log-log plot for

$\Delta T \sim {V^{ - 1}}$

boundary line between the retraction-bouncing mode and the bouncing-with-spray mode in the film boiling regime for water droplets, the droplet diameter is 2.2 mm; (e) regime map for a water droplet impacting on a heated silicon wafer and the transition between the retraction-bouncing mode and the bouncing-with-spray mode, the red dashed line is only for eye guidance, the droplet diameter is 2.2 mm.

下载: 全尺寸图片幻灯片

蒸汽泡的形成和增大是由于壁面对液滴的热传递造成的. 考虑到能量守恒, 壁面传递的热量应该等于蒸汽泡的汽化潜热, 即

$h\Delta T\Delta \tau A \sim \frac{2}{3}{\text{π}}r_{\rm{v}}^3\ell {\rho _{\rm{v}}},$

(3)

其中, h为平均传热系数, $\Delta T = T - {T_{\rm s}}$ 为壁面的过热度, T和 ${T_{\rm s}}$ 分别为壁面的温度及液体的饱和温度, A为蒸汽泡的热传递面积, $\ell$ 为液滴的汽化潜热, 下标v表示“蒸汽”. 热传递的时间间隔 $\Delta \tau$ 可由液滴达到最大铺展时的时间 ${\tau _0}$ 给出, 如图7(b)所示, 其并不随着撞击速度的改变而改变, 并且和液滴的自由振荡周期成比例^[24]. 在临界条件下, 蒸汽泡的半径 ${r_{\rm{v}}}$ 和液膜的最小厚度 ${b_{\min }}$ 相当, 其中, 液膜的最小厚度可以表示为^[25,26]

${b_{\min }} \sim R{e^{{{ - 2} / 5}}}{D_0}.$

(4)

为了获得(3)式中的传热系数h, 考虑在膜态沸腾下, 蒸汽层内部的热传递主要以导热形式进行^[27], 假设液滴底部蒸汽层具有均匀厚度 $\delta$ , 则通过蒸汽层单位面积的热流密度可以表示为 ${{{\lambda_{\rm{v}}}\Delta T} / \delta }$ , 其中 ${\lambda_{\rm{v}}}$ 为蒸汽的热导率, 因此液滴流体的汽化质量速率可以表示为 ${{{\lambda_{\rm{v}}}\Delta T} / {\delta \ell }}$ . 考虑在蒸汽层半径r的区域内, 向外流动的蒸汽和液滴产生的蒸汽体积平衡, 有

${\text{π}}{r^2}\frac{{{\lambda_{\rm{v}}}\Delta T}}{{\delta \ell }} = 2{\text{π}}ru\left( r \right)\delta {\rho_{\rm{v}}}.$

(5)

故蒸汽的径向速度可以表示为

$u\left( r \right) = \frac{{{\lambda_{\rm{v}}}\Delta T}}{{2{\rho_{\rm{v}}}\ell }}\frac{r}{{{\delta ^2}}}.$

(6)

由于蒸汽层的厚度很小, 故可以将其内部的流动视为泊肃叶流动. 因此考虑蒸汽的径向流动为黏性流动, 则蒸汽层内的平均径向速度可以表示为

$u\left( r \right) = - \frac{{{\delta ^2}}}{{12{\mu_{\rm{v}}}}}\frac{{{\rm{d}}p}}{{{\rm{d}}r}}.$

(7)

把(6)式代入到(7)式中, 并对其进行积分, 可以获得蒸汽层内压力分布为

$p(r) = {p_a} + \frac{{3{\lambda_{\rm{v}}}{\mu_{\rm{v}}}\Delta T}}{{{\rho_{\rm{v}}}\ell }}\frac{{{R^2} - {r^2}}}{{{\delta ^4}}},$

(8)

其中, ${p_a}$ 为大气压力, R为液膜的半径. 通过对压力梯度 $\left( {p\left( r \right) - {p_a}} \right)$ 进行积分, 可以获得使液膜上升的力, 并且这个力与液滴的重力和浮力之差 $\left[ {{\text{π}}{R^2}b\left( {{\rho _{\rm{l}}} - {\rho_{\rm{v}}}} \right)g} \right]$ 平衡, 如图7(a)所示. 则蒸汽层的厚度可以表示为

$\delta = {\left[ {\frac{{3{\lambda_{\rm{v}}}{\mu_{\rm{v}}}\Delta T{R^2}}}{{2{\rho_{\rm{v}}}\left( {{\rho _{\rm{l}}} - {\rho_{\rm{v}}}} \right)b\ell g}}} \right]^{\frac{1}{4}}}.$

(9)

因此, 平均传热系数可以表示为

$h = \frac{{{\lambda_{\rm{v}}}}}{\delta } = {\left[ {\frac{{2\lambda_{\rm{v}}^3{\rho_{\rm{v}}}\left( {{\rho _{\rm{l}}} - {\rho_{\rm{v}}}} \right)b\ell g}}{{3{\mu_{\rm{v}}}\Delta T{R^2}}}} \right]^{\frac{1}{4}}}.$

(10)

在从回缩弹起模式到雾化弹起模式过渡的临界条件下, (10)式中的液膜半径R可以近似为液滴的最大铺展半径, 即

${{{R_{\max }}} / {{R_0}}} = {{{D_{\max }}} / {{D_0}}} \sim W{e^n}.$

(11)

如图7(c)所示, 其中 $n = 0.35$ . 事实上, 关于液滴撞击加热壁面情况下液滴铺展尺度的研究可以在大量的文献中找到^{[10,12,22,23,28]}, 并且指数n会因壁面的结构不同而不同, 本实验所求得的指数n与Tran等^[10]的研究接近, 即 $n = 0.39$ .

通过将(10)式代入(3)式, 并分别代入(4)和(11)式所求得临界条件下的b和R, 可以得到

${\left[\!{\frac{{\lambda_{\rm{v}}^3\left( {{\rho _l} \!-\! {\rho_{\rm{v}}}}\!\right)g\Delta {T^3}}}{{{\mu_{\rm{v}}}{\ell ^3}\rho_{\rm{v}}^3}}} \right]^{\frac{1}{4}}}{\tau _0}A \sim W{e^{\frac{n}{2}}}R{e^{ - 1.1}}R_0^{\frac{{13}}{4}}.$

(12)

引入一个修正的格拉晓夫数 $G{r^*} =$ ${{g{\rho_{\rm{v}}}\left( {{\rho _{\rm{l}}} - {\rho_{\rm{v}}}} \right)R_0^3} / {\mu_{\rm{v}}^2}}$ , 蒸汽的普朗特数 $P{r_{\rm{v}}} =$ ${{{c_{p{\rm{v}}}}{\mu_{\rm{v}}}} / {{\lambda_{\rm{v}}}}}$ , 雅各比数 $Ja = {{{c_{p{\rm{l}}}}\Delta T} / \ell }$ , 以及一个修正的傅里叶数 $F{o^*} = {{{\lambda_{\rm{v}}}{\tau _0}A} / {({\rho_{\rm{v}}}}}{c_{p{\rm{v}}}}R_0^4)$ , 可以得出如下无量纲关系式:

${\left( {{{{c_{p{\rm{v}}}}} / {{c_{p{\rm{l}}}}}}} \right)^{\frac{3}{4}}}{\left( {G{r^*}P{r_{\rm{v}}}} \right)^{\frac{1}{4}}}F{o^*}J{a^{\frac{3}{4}}} \sim W{e^{\frac{n}{2}}}R{e^{ - 1.1}}.$

(13)

对于不同的液体, 假定蒸汽泡的传热面积A为定值. 因此, $G{r^*}$ , $P{r_{\rm{v}}}$ 和 $F{o^*}$ 主要受液体的物性及液滴初始尺寸的影响, 而对于相同液体及相同的初始液滴尺寸, 液滴的物性参数及液滴达到最大铺展尺寸的时间 ${\tau _0}$ 可视为常数. 代入图7(c)中的 $n = 0.35$ , 则可以得到

$\Delta T \sim {V^{ - 1}}.$

(14)

为了验证该幂次关系的有效性, 膜态沸腾下的实验数据在对数坐标下给出, 如图7(d)所示, 可以看到该幂次关系与实验数据吻合良好. 另外, 该理论模型同样可以解释在过渡沸腾区域, 从回缩弹起到雾化弹起模式过渡的壁面温度随着撞击韦伯数的增加而增加, 如图7(e)所示的相图. 在过渡沸腾区域, 壁面传递给液滴的热流密度会随着壁面温度的上升而下降, 故在一定时间间隔 ${\tau _0}$ 内, 蒸汽泡的最大尺寸会随着壁面温度的升高而减小, 因此液膜需要变得更薄, 即更大的韦伯数, 以满足临界条件下 ${r_{\rm v}} \sim {b_{\min }}$ . 所以在过渡沸腾区域, 从回缩弹起到雾化弹起模式过渡的撞击韦伯数会随着壁面温度的上升而增加.

4. 结　论

本文通过高速摄影技术对不同液体种类的液滴撞击加热壁面进行了研究, 主要结论如下:

1)发现了一种由于液膜破碎而导致的弹起模式, 即雾化弹起模式, 该模式能够极大地减少液滴在壁面上的驻留时间, 并且与复杂的表面微结构加工或者表面改性的方式相比, 本文中仅依赖加热壁面的方式以减少液滴驻留时间的方法更容易实现和可靠;

2)雾化弹起模式下液滴在壁面上的驻留时间不受液滴黏度的影响, 和液滴自由振荡周期成比例, 并且和回缩弹起模式相比, 雾化弹起模式下液滴的无量纲驻留时间能减少约40%;

3)揭示了雾化弹起模式下液滴驻留时间减少的机理, 即由于液膜内部蒸汽泡的破碎导致孔洞的形成, 进而液膜从孔洞处回缩, 由于液滴的回缩距离减少, 因此其回缩时间减少, 最终导致了液滴在壁面上的驻留时间减少;

4)依据能量守恒及蒸汽泡破碎的临界条件 ${r_{\rm v}} \sim {b_{\min }}$ 提出了一个简化理论模型, 依据该模型得出了膜态沸腾区域两种弹起模式的过渡边界的幂次关系, 即 $\Delta T \sim {V^{ - 1}}$ , 并且与实验数据吻合良好;

5)在液体种类和初始液滴尺寸相同的条件下, 依据该理论模型对过渡沸腾区域的两种弹起模式的过渡边界进行了讨论.

参考文献

[1]	毕菲菲, 郭亚丽, 沈胜强, 陈觉先, 李熠桥 2012 物理学报 61 184702 Google Scholar Bi F F, Guo Y L, Shen S Q, Chen J X, Li Y Q 2012 Acta Phys. Sin. 61 184702 Google Scholar
[2]	马理强, 常建忠, 刘汉涛, 刘谋斌 2012 物理学报 61 054701 Google Scholar Ma L Q, Chang J Z, Liu H T, Liu M B 2012 Acta Phys. Sin. 61 054701 Google Scholar
[3]	沈胜强, 张洁珊, 梁刚涛 2015 物理学报 64 134704 Google Scholar Shen S Q, Zhang J S, Liang G T 2015 Acta Phys. Sin. 64 134704 Google Scholar
[4]	Quéré D 2013 Annu. Rev. Fluid Mech. 45 197 Google Scholar
[5]	Bernardin J D, Stebbins C J, Mudawar I 1997 Int. J. Heat Mass Transfer 40 247 Google Scholar
[6]	Castanet G, Lienart T, Lemoine F 2009 Int. J. Heat Mass Transfer 52 670 Google Scholar
[7]	Bertola V 2015 Int. J. Heat Mass Transfer 85 430 Google Scholar
[8]	Breitenbach J, Roisman I V, Tropea C 2017 Phy. Rev. Fluids 2 074301 Google Scholar
[9]	Liang G T, Mudawar I 2017 Int. J. Heat Mass Transfer 106 103 Google Scholar
[10]	Tran T, Staat H J J, Prospeeretti A, Sun C, Lohse D 2012 Phys. Rev. Lett. 108 036101 Google Scholar
[11]	Nair H, Staat H J J, Tran T, van Houselt A, Prosperetti A, Lohse D, Sun C 2014 Soft Matter 10 2102 Google Scholar
[12]	Tran T, Staat H J J, Susarrey-Arce A, Foertsch T C, van Houselt A, Gardeniers H J G E, Prosperetti A, Lohse D, Sun C 2013 Soft Matter 9 3272 Google Scholar
[13]	Shirota M, van Limbeek M A J, Sun C, Prosperetti A, Lohse D 2016 Phys. Rev. Lett. 116 064501 Google Scholar
[14]	van Limbeek M A J, Hoefnagels P B J, Sun C, Lohse D 2017 Soft Matter 13 7514 Google Scholar
[15]	Richard D, Clanet C, Quéré D 2002 Nature 417 811
[16]	Bird J C, Dhiman R, Kwon H M, Varanasi K K 2013 Nature 503 385 Google Scholar
[17]	Liu Y, Moevius L, Xu X, Qian T, Yeomans J M, Wang Z 2014 Nat. Phys. 10 515 Google Scholar
[18]	Liu Y, Andrew M, Li J, Yeomans J M, Wang Z 2015 Nat. Commun. 6 10034 Google Scholar
[19]	Yun S 2018 Soft Matter 14 4946 Google Scholar
[20]	Roisman I V, Horvat K, Tropea C 2006 Phys. Fluids 18 102104 Google Scholar
[21]	Rayleigh L 1879 Proc. Roy. Soc. London 29 71 Google Scholar
[22]	Biance A L, Chevy F, Clanet C, Lagubeau G, Quéré D 2006 J. Fluid Mech. 554 47 Google Scholar
[23]	Liang G T, Shen S Q, Guo Y L, Zhang J L 2016 Int. J. Heat Mass Transfer 100 48 Google Scholar
[24]	Biance A L, Pirat C, Ybert C 2011 Phys. Fluids 23 022104 Google Scholar
[25]	Eggers J, Fontelos M A, Josserand C, Zaleski S 2010 Phys. Fluids 22 062101 Google Scholar
[26]	Fedorchenko A I, Wang A B, Wang Y H 2005 Phys. Fluids 17 093104 Google Scholar
[27]	Labeish V G 1994 Exp. Therm. Fluid Sci. 8 181 Google Scholar
[28]	Antonimi C, Bernagozzi I, Jung S, Poilikakos D, Marengo M 2013 Phys. Rev. Lett. 111 014501 Google Scholar

施引文献

期刊类型引用(9)

1.	李聪，王佳利，王志敏，李源慧，陶振翔，吴建松，杨溢. 液滴撞击加热壁面的实验平台设计与教学应用. 实验技术与管理. 2025(01): 238-245 . 百度学术
2.	冯山青，龚路远，权生林，郭亚丽，沈胜强. 纳米液滴撞击高温平板壁的分子动力学模拟. 物理学报. 2024(10): 113-122 . 百度学术
3.	鲁杰，李亚磊，徐龙，郝继光. 液滴撞击倾斜表面铺展研究. 实验流体力学. 2023(06): 42-50 . 百度学术
4.	王开珉，张玉杰，康培森，刘宏升，刘晓华. 液滴撞击加热亲水管壁后的反弹和中心射流. 浙江大学学报(工学版). 2022(06): 1191-1198 . 百度学术
5.	李林豪，杨秀峰，SOHAG Md M.A.，王晓亮，刘青泉. 双液滴连续撞击高温壁面的SPH-ASR模拟研究. 中国科学:物理学力学天文学. 2022(10): 37-49 . 百度学术
6.	周鑫，马小晶，胡丽娜，刘佳，张博文. 不同壁面条件下液滴撞击铺展特性的模拟研究. 计算力学学报. 2022(06): 761-767 . 百度学术
7.	春江，王瑾萱，徐晨，温荣福，兰忠，马学虎. 液滴撞击超亲水表面的最大铺展直径预测模型. 物理学报. 2021(10): 248-258 . 百度学术
8.	沈学峰，曹宇，王军锋，刘海龙. 剪切变稀液滴撞击不同浸润性壁面的数值模拟研究. 物理学报. 2020(06): 158-167 . 百度学术
9.	王凯，冮铁强. 当量滴径对于液滴驻留液面时间尺度影响的研究. 科技创新导报. 2020(18): 91-92+94 . 百度学术

其他类型引用(6)

图 1 液滴撞击加热壁面的实验装置示意图　(a)水平视角; (b)倾斜视角

Fig. 1. Schematic diagram of the experimental setup for droplet impact on a heated surface: (a) Side-view; (b) aerial-view.

下载: 全尺寸图片幻灯片

图 2 在3个沸腾区域水滴撞击加热壁面的形貌图, 子液滴的喷射如红色圆圈及箭头所示　(a)核态沸腾区域, $We = 37$ , $T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (b)过渡沸腾区域回缩弹起模式, $We = 37$ , $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (c)过渡沸腾区域雾化弹起模式, 其中(c1)为 $We = 69$ , $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ , (c2)为 $We = 105$ , $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (d)膜态沸腾区域回缩弹起模式, 其中(d1)为 $We = 25$ , $T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ , (d2)为 $We = 25$ , $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (e)膜态沸腾区域雾化弹起模式, 其中(e1)为 $We = 60$ , $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ , (e2)为 $We = 129$ , $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$

Fig. 2. Morphologies of water droplet impacting on heated surfaces in three boiling regimes, the red circles and arrows highlight the spray-like ejection of tiny droplets: (a) Nucleate boiling, $We = 37$ and $T = 260{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (b) retraction-bouncing mode in transition boiling regime, $We = 37$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ; (c) bouncing-with-spray mode in transition boiling regime, $We = 69$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (c1), and $We = 105$ and $T = 380{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (c2); (d) retraction-bouncing mode in film boiling regime, $We = 25$ and $T = 400{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (d1), and $We = 25$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (d2); (e) bouncing-with-spray mode in film boiling regime, $We = 60$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (e1), and $We = 129$ and $T = 440{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ for (e2).

下载: 全尺寸图片幻灯片

图 3 液体黏度对雾化弹起模式下驻留时间的影响(液滴直径为2.8 mm, 撞击速度为 $1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$ )

Fig. 3. Effect of liquid viscosity on the residence time of the bouncing-with-spray mode. The droplet diameter is 2.8 mm, and the impact velocity is $1.2{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm m} \cdot {{\rm s}^{{\rm - 1}}}$ .

下载: 全尺寸图片幻灯片

图 4 回缩弹起模式(空心圆)和雾化弹起模式(实心圆)下的无量纲驻留时间比较

Fig. 4. Comparison of the dimensionless residence time in the retraction-bouncing mode (hollow circles) and in the bouncing-with-spray mode (solid dots).

下载: 全尺寸图片幻灯片

图 5 倾斜视角下雾化弹起模式的序列图片( $We = 82$ , $T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ )

Fig. 5. Image series of the bouncing-with-spray mode in aerial view ( $We = 82$ and $T = 450{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ).

下载: 全尺寸图片幻灯片

图 6 更大的放大倍率下液膜孔洞形成的图片序列(We =122, $T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ )

Fig. 6. Image series of hole formation in the liquid film using greater magnification ( $We = 122$ and $T = 480{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ^\circ {\rm C}$ ).

下载: 全尺寸图片幻灯片

图 7 (a)蒸汽泡破碎理论模型的示意图; (b)采用液滴自由振荡周期无量纲化后的液滴最大铺展时间; (c)对数坐标下膜态沸腾区域水滴的无量纲最大铺展直径和撞击韦伯数的关系图, ${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$ ; (d)对数坐标下水滴在膜态沸腾区域由回缩弹起到雾化弹起模式过渡的 $\Delta T \sim {V^{ - 1}}$ 边界线, 液滴直径为2.2 mm; (e)水滴撞击加热硅片的相图及从回缩弹起到雾化弹起模式下的过渡, 图中红色虚线为示意, 液滴直径为2.2 mm, “过渡”为过渡腾区域, “膜态”指膜态沸腾区域

Fig. 7. (a) Schematic diagram for the theoretical model for the burst of vapor bubbles; (b) the time of maximum spreading diameter of water droplets normalized by droplet free oscillation period; (c) log-log plot for the normalized maximal spreading diameter of water droplets versus We in the film boiling regime, showing ${{{D_{\max }}} / {{D_0}}} \sim W{e^n}$ ; (d) log-log plot for $\Delta T \sim {V^{ - 1}}$ boundary line between the retraction-bouncing mode and the bouncing-with-spray mode in the film boiling regime for water droplets, the droplet diameter is 2.2 mm; (e) regime map for a water droplet impacting on a heated silicon wafer and the transition between the retraction-bouncing mode and the bouncing-with-spray mode, the red dashed line is only for eye guidance, the droplet diameter is 2.2 mm.

下载: 全尺寸图片幻灯片

表 1 本实验中不同液体物性参数表

Table 1. Properties of the liquids used in this study.

液体种类	密度 $\rho {\rm /kg} \cdot {{\rm m}^{{\rm - 3}}}$	表面张力 $\sigma /{\rm mN} \cdot {{\rm m}^{ - 1}}$	动力黏度 $\mu /{\rm mPa} \cdot {\rm s}$	汽化潜热 $\ell /{\rm kJ} \cdot {\rm k}{{\rm g}^{ - 1}}$
水	998	72	1.00	2257
水 + 甘油(28 wt%)	1068	70	2.35
水 + 甘油(34 wt%)	1083	70	3.00
水 + 甘油(39 wt%)	1097	70	3.63
水 + 甘油(56 wt%)	1143	68	9.00
水 + 甘油(75 wt%)	1195	66	42.47
乙醇	790	22	1.07	853

下载: 导出CSV

[1]	毕菲菲, 郭亚丽, 沈胜强, 陈觉先, 李熠桥 2012 物理学报 61 184702 Google Scholar Bi F F, Guo Y L, Shen S Q, Chen J X, Li Y Q 2012 Acta Phys. Sin. 61 184702 Google Scholar
[2]	马理强, 常建忠, 刘汉涛, 刘谋斌 2012 物理学报 61 054701 Google Scholar Ma L Q, Chang J Z, Liu H T, Liu M B 2012 Acta Phys. Sin. 61 054701 Google Scholar
[3]	沈胜强, 张洁珊, 梁刚涛 2015 物理学报 64 134704 Google Scholar Shen S Q, Zhang J S, Liang G T 2015 Acta Phys. Sin. 64 134704 Google Scholar
[4]	Quéré D 2013 Annu. Rev. Fluid Mech. 45 197 Google Scholar
[5]	Bernardin J D, Stebbins C J, Mudawar I 1997 Int. J. Heat Mass Transfer 40 247 Google Scholar
[6]	Castanet G, Lienart T, Lemoine F 2009 Int. J. Heat Mass Transfer 52 670 Google Scholar
[7]	Bertola V 2015 Int. J. Heat Mass Transfer 85 430 Google Scholar
[8]	Breitenbach J, Roisman I V, Tropea C 2017 Phy. Rev. Fluids 2 074301 Google Scholar
[9]	Liang G T, Mudawar I 2017 Int. J. Heat Mass Transfer 106 103 Google Scholar
[10]	Tran T, Staat H J J, Prospeeretti A, Sun C, Lohse D 2012 Phys. Rev. Lett. 108 036101 Google Scholar
[11]	Nair H, Staat H J J, Tran T, van Houselt A, Prosperetti A, Lohse D, Sun C 2014 Soft Matter 10 2102 Google Scholar
[12]	Tran T, Staat H J J, Susarrey-Arce A, Foertsch T C, van Houselt A, Gardeniers H J G E, Prosperetti A, Lohse D, Sun C 2013 Soft Matter 9 3272 Google Scholar
[13]	Shirota M, van Limbeek M A J, Sun C, Prosperetti A, Lohse D 2016 Phys. Rev. Lett. 116 064501 Google Scholar
[14]	van Limbeek M A J, Hoefnagels P B J, Sun C, Lohse D 2017 Soft Matter 13 7514 Google Scholar
[15]	Richard D, Clanet C, Quéré D 2002 Nature 417 811
[16]	Bird J C, Dhiman R, Kwon H M, Varanasi K K 2013 Nature 503 385 Google Scholar
[17]	Liu Y, Moevius L, Xu X, Qian T, Yeomans J M, Wang Z 2014 Nat. Phys. 10 515 Google Scholar
[18]	Liu Y, Andrew M, Li J, Yeomans J M, Wang Z 2015 Nat. Commun. 6 10034 Google Scholar
[19]	Yun S 2018 Soft Matter 14 4946 Google Scholar
[20]	Roisman I V, Horvat K, Tropea C 2006 Phys. Fluids 18 102104 Google Scholar
[21]	Rayleigh L 1879 Proc. Roy. Soc. London 29 71 Google Scholar
[22]	Biance A L, Chevy F, Clanet C, Lagubeau G, Quéré D 2006 J. Fluid Mech. 554 47 Google Scholar
[23]	Liang G T, Shen S Q, Guo Y L, Zhang J L 2016 Int. J. Heat Mass Transfer 100 48 Google Scholar
[24]	Biance A L, Pirat C, Ybert C 2011 Phys. Fluids 23 022104 Google Scholar
[25]	Eggers J, Fontelos M A, Josserand C, Zaleski S 2010 Phys. Fluids 22 062101 Google Scholar
[26]	Fedorchenko A I, Wang A B, Wang Y H 2005 Phys. Fluids 17 093104 Google Scholar
[27]	Labeish V G 1994 Exp. Therm. Fluid Sci. 8 181 Google Scholar
[28]	Antonimi C, Bernagozzi I, Jung S, Poilikakos D, Marengo M 2013 Phys. Rev. Lett. 111 014501 Google Scholar

[1]	王浩, 徐进良. 油面上相邻Leidenfrost液滴的相互作用及运动机制. 物理学报, 2023, 72(5): 054401. doi: 10.7498/aps.72.20221822
[2]	唐修行, 陈泓樾, 王婧婧, 王志军, 臧渡洋. 表面活性剂液滴过渡沸腾的Marangoni效应与二次液滴形成. 物理学报, 2023, 72(19): 196801. doi: 10.7498/aps.72.20230919
[3]	李逢超, 付宇, 李超, 杨建刚, 胡春波. 铝液滴撞击曲面的流动特性分析. 物理学报, 2022, 71(18): 184701. doi: 10.7498/aps.71.20220442
[4]	潘伶, 张昊, 林国斌. 纳米液滴撞击柱状固体表面动态行为的分子动力学模拟. 物理学报, 2021, 70(13): 134704. doi: 10.7498/aps.70.20210094
[5]	刘联胜, 刘轩臣, 贾文琪, 田亮, 杨华, 段润泽. 小液滴撞击壁面传热特性数值分析. 物理学报, 2021, 70(7): 074402. doi: 10.7498/aps.70.20201354
[6]	沈学峰, 曹宇, 王军锋, 刘海龙. 剪切变稀液滴撞击不同浸润性壁面的数值模拟研究. 物理学报, 2020, 69(6): 064702. doi: 10.7498/aps.69.20191682
[7]	杨亚晶, 梅晨曦, 章旭东, 魏衍举, 刘圣华. 液滴撞击液膜的穿越模式及运动特性. 物理学报, 2019, 68(15): 156101. doi: 10.7498/aps.68.20190604
[8]	赵可, 佘阳梓, 蒋彦龙, 秦静, 张振豪. 液氮滴撞击壁面相变行为的数值研究. 物理学报, 2019, 68(24): 244401. doi: 10.7498/aps.68.20190945
[9]	李玉杰, 黄军杰, 肖旭斌. 液滴撞击圆柱内表面的数值研究. 物理学报, 2018, 67(18): 184701. doi: 10.7498/aps.67.20180364
[10]	黄虎, 洪宁, 梁宏, 施保昌, 柴振华. 液滴撞击液膜过程的格子Boltzmann方法模拟. 物理学报, 2016, 65(8): 084702. doi: 10.7498/aps.65.084702
[11]	白玲, 李大鸣, 李彦卿, 王志超, 李杨杨. 基于范德瓦尔斯表面张力模式液滴撞击疏水壁面过程的研究. 物理学报, 2015, 64(11): 114701. doi: 10.7498/aps.64.114701
[12]	戴剑锋, 樊学萍, 蒙波, 刘骥飞. 单液滴撞击倾斜液膜飞溅过程的耦合Level Set-VOF模拟. 物理学报, 2015, 64(9): 094704. doi: 10.7498/aps.64.094704
[13]	沈胜强, 张洁珊, 梁刚涛. 液滴撞击加热壁面传热实验研究. 物理学报, 2015, 64(13): 134704. doi: 10.7498/aps.64.134704
[14]	郭亚丽, 魏兰, 沈胜强, 陈桂影. 双液滴撞击平面液膜的流动与传热特性. 物理学报, 2014, 63(9): 094702. doi: 10.7498/aps.63.094702
[15]	梁刚涛, 沈胜强, 郭亚丽, 陈觉先, 于欢, 李熠桥. 实验观测液滴撞击倾斜表面液膜的特殊现象. 物理学报, 2013, 62(8): 084707. doi: 10.7498/aps.62.084707
[16]	梁刚涛, 郭亚丽, 沈胜强. 液滴撞击液膜的射流与水花形成机理分析. 物理学报, 2013, 62(2): 024705. doi: 10.7498/aps.62.024705
[17]	梁刚涛, 郭亚丽, 沈胜强. 液滴低速撞击润湿球面现象观测分析. 物理学报, 2013, 62(18): 184703. doi: 10.7498/aps.62.184703
[18]	毕菲菲, 郭亚丽, 沈胜强, 陈觉先, 李熠桥. 液滴撞击固体表面铺展特性的实验研究. 物理学报, 2012, 61(18): 184702. doi: 10.7498/aps.61.184702
[19]	张林, 孔红艳, 杨国健. 约束阱中受激发原子的集体反弹效应所导致的自组织行为. 物理学报, 2006, 55(10): 5122-5128. doi: 10.7498/aps.55.5122
[20]	刘辽, 李鉴曾. 迹反常与时间对称性反弹宇宙模型. 物理学报, 1982, 31(12): 96-99. doi: 10.7498/aps.31.96

期刊类型引用(9)

1.	李聪，王佳利，王志敏，李源慧，陶振翔，吴建松，杨溢. 液滴撞击加热壁面的实验平台设计与教学应用. 实验技术与管理. 2025(01): 238-245 . 百度学术
2.	冯山青，龚路远，权生林，郭亚丽，沈胜强. 纳米液滴撞击高温平板壁的分子动力学模拟. 物理学报. 2024(10): 113-122 . 百度学术
3.	鲁杰，李亚磊，徐龙，郝继光. 液滴撞击倾斜表面铺展研究. 实验流体力学. 2023(06): 42-50 . 百度学术
4.	王开珉，张玉杰，康培森，刘宏升，刘晓华. 液滴撞击加热亲水管壁后的反弹和中心射流. 浙江大学学报(工学版). 2022(06): 1191-1198 . 百度学术
5.	李林豪，杨秀峰，SOHAG Md M.A.，王晓亮，刘青泉. 双液滴连续撞击高温壁面的SPH-ASR模拟研究. 中国科学:物理学力学天文学. 2022(10): 37-49 . 百度学术
6.	周鑫，马小晶，胡丽娜，刘佳，张博文. 不同壁面条件下液滴撞击铺展特性的模拟研究. 计算力学学报. 2022(06): 761-767 . 百度学术
7.	春江，王瑾萱，徐晨，温荣福，兰忠，马学虎. 液滴撞击超亲水表面的最大铺展直径预测模型. 物理学报. 2021(10): 248-258 . 百度学术
8.	沈学峰，曹宇，王军锋，刘海龙. 剪切变稀液滴撞击不同浸润性壁面的数值模拟研究. 物理学报. 2020(06): 158-167 . 百度学术
9.	王凯，冮铁强. 当量滴径对于液滴驻留液面时间尺度影响的研究. 科技创新导报. 2020(18): 91-92+94 . 百度学术

其他类型引用(6)

计量

文章访问数: 9834
PDF下载量: 142
被引次数: 15

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板