EXPERIMENTAL STUDY ON VISCOUS NEWTONIAN DROPLET IMPACTS ON DRY OR PRE-WETTED MESHES
-
摘要: 液滴撞击网面现象广泛存在于自然界和一系列应用中, 液滴撞网后会穿透破碎产生二次液滴或不破碎全部附着在网面上, 两种情况下都会残留液体在网面而形成预湿, 影响后续撞击结果, 但前人研究集中于低黏性液滴撞击干燥网面, 黏性牛顿流体液滴撞击干燥或预湿网面的演化与机理仍有待探索. 文章采用高速阴影成像技术, 研究了黏性液滴(甘油水溶液)撞击干燥和预湿网面形成液指和破碎的演化规律, 考虑了网面结构尺寸、液滴黏性及撞击前网面上预湿液膜厚度对撞击结果的影响. 实验结果表明, 液滴撞击干燥网面后形成液指的最大长度随网孔宽度降低、液滴黏性增加而减小; 液滴黏性增加、网孔宽度减小均会抑制液滴对干燥网面的完全穿透; 预湿液膜高度的增加抑制液滴对网面的完全穿透, 并使不完全穿透时形成液指的最大长度减小. 建立了考虑液滴黏性、网孔宽度和网面预湿的液滴撞击网面后不完全穿透时形成液指的最大长度预测模型, 以及出现完全穿透时的临界参数理论预测模型, 模型预测结果均与实验结果吻合良好.Abstract: Droplet impacts on meshes are ubiquitous both in nature and in a variety of applications. The impact may lead to liquid penetration through the mesh and formation of secondary droplets underneath the mesh, or spreading on the mesh and no occurrence of penetration. In either situation, liquid remains on the meshes and forms pre-wetted meshes after impacts, leading to different following impact outcomes compared with impacts on dry meshes. However, previous studies focused on impacts of low-viscosity droplets on dry meshes. The evolution and mechanism of viscous Newtonian droplets impacts on dry or pre-wetted meshes remain to be explored. In this paper, the liquid fingers and the fragmentation occurred underneath the mesh following a viscous droplet (aqueous glycerol solution) impacting a dry or pre-wetted mesh are investigated using high-speed shadow imaging technology, with special attention paid on the influence of mesh size, droplet viscosity and pre-wetted liquid film thickness on impact outcomes. It was observed that both a decrease of mesh size and an increase of droplet viscosity resulted in a decrease of maximum length of liquid finger and suppressed a complete penetration through a dry mesh, an increase of pre-wetted liquid film thickness suppressed a complete penetration and resulted in a decrease of maximum length of liquid finger. Considering the influence of mesh size, droplet viscosity, and pre-wetted liquid film thickness, theoretical models for predicting the maximum length of liquid finger when incomplete penetration occurs and for predicting the threshold parameters for complete penetration is proposed and is validated by comparisons with experimental results.
-
Keywords:
- viscous droplets /
- pre-wetted meshes /
- liquid fingers /
- secondary droplets
-
引 言
液滴撞击网面现象广泛存在于自然界和一系列工农业应用中, 如喷涂[1-6]、防护衣物开发[7]、传热[8]、印刷[9-11]、两相分离[12-17]、口罩的开发与使用[18-21]和农药喷洒[3]等. 液滴撞击网面经历复杂的动态过程, 除了在网面上产生液滴撞击固体表面时的飞溅、铺展和反弹[22-33]现象外, 还会在网面下产生穿透和破碎, 形成二次液滴. 液滴撞网结果通常分为以下3种: 韦伯数较低时, 撞击液滴在网面上扩散, 不穿透网面, 定义为不穿透; 在韦伯数较高时, 液滴撞击网面后在网面下方形成液指, 液指不断裂并最终回缩到网面上表面, 定义为不完全穿透; 进一步提高韦伯数到大于某一临界值Wep, 液指会破裂并在网面下产生二次液滴, 这种现象被定义为完全穿透[34-38].
液滴撞击网面产生的二次液滴在有些场合是需要的, 在另一些场合又会产生负面影响, 因而工程中存在认清其规律和机理, 从而实现调控的迫切需求. 如, 对于防护衣物开发和口罩的使用, 二次液滴的出现会显著降低防护效果; 而对于喷涂和印刷, 二次液滴是有助于用更少液体实现更大覆盖面积的目标的. 临界韦伯数Wep通常被用来界定从不完全穿透到完全穿透的转变, 当撞击韦伯数大于Wep时, 完全穿透现象出现, 撞击会形成二次液滴; 反之, 撞击不会形成二次液滴[34]. 对于特定的液滴, 也可用临界速度Up来界定两种现象的转变[35,38].
随着高速摄影技术的发展, 直接观测液滴撞网现象变得可能, 对该现象的研究在过去的20年中引起了学者们的广泛关注. Lorenceau等[38]对液滴撞击单孔进行研究, 首次定义了Up, 并提出了一个基于动压和毛细压力平衡的穿透预测模型, 还根据实验结果拟合了模型中的系数. 此后, 类似模型被广泛应用于分析网面材料[37,39-43]、网面润湿性[35-36,43-48]、网面形态[49-52]、网面倾角[53]和网面预湿[19,34]等因素对临界韦伯数Wep或速度Up的影响, 建立了适应这些特定因素的穿透预测模型[34-36,52].
前述研究都是基于低黏性液滴(大部分是水)撞击开展的, 关注点仍是网面属性对撞击结果的影响. 液滴属性[54-58]对撞击结果的影响最近也引起了学者们的注意, 如Blackwell等[54]、Mehrizi等[55]研究了非牛顿流体液滴撞击网面的现象和规律; Wang等[56]和Vontas等[57]采用数值模拟方法研究了液滴黏性的影响; Abouelsoud等[58]实验研究了不同黏性的低表面张力硅油滴撞击亲油网的现象和规律, 预测了最大液指长度. 然而, 前述研究都没有关注网面预湿的影响, 实际应用中的网面只要经过一次液滴撞击, 就会被预湿[19,34]; 针对预湿的研究[19,34]虽有开展, 但是都没有研究液滴黏性的影响, 尤其是大表面张力液滴黏性的影响.
本文利用高速阴影成像技术, 研究了大表面张力液滴(甘油水溶液)黏性、网面尺寸和网面预湿对液滴撞网结果的影响, 获得了系列条件下的最大液指长度和形成完全穿透所需的临界参数, 分析了关键因素的影响规律和机制, 建立了考虑液滴黏性、网面尺寸和网面预湿的最大液指长度和形成完成穿透所需临界参数的预测模型, 预测结果与实验结果吻合良好, 为工程应用提供了实验依据和理论工具.
1. 实验设置
实验装置如图1所示. 利用注射泵缓慢驱动注射器在毛细管口产生甘油水溶液液滴, 液滴在重力作用下加速并撞击网面. 实验中用到了10种具有不同甘油质量浓度ψ的甘油水溶液, 分别使用不同质量比的分析纯甘油和去离子水混合形成, 其物理属性[59]如表1所示. 液滴由相同尺寸的毛细管生成, 直径D是由图像处理获得液滴接触网面前的水平宽度Dx和垂直高度Dy, 再由式(1)计算得出(D = 2.5 ± 0.15 mm). 液滴撞击速度U是由图像处理获得液滴与网面接触前10帧内液滴上下沿移动距离的平均值, 再除以10帧对应的飞行时间计算得出, 速度通过升降平台改变针头到网面的距离X来调控(U = 0.96 ~ 3.82 m/s)
表 1 甘油水液滴的物理属性(25 °C)Table 1. The physical properties of aqueous glycerol solutions at 25 °Cψ/wt% μ/(mPa·s) ρ/(kg·m−3) σ/(mN·m−1) 50 5.01 1113.41 68.05 64 11.54 1150.67 66.95 70 18.07 1167.41 66.85 75 27.73 1181.74 66.54 78 36.94 1190.50 66.16 80 45.37 1196.42 66.77 81.5 53.35 1200.90 65.43 83 63.16 1205.41 65.52 84 71.06 1208.43 66.20 85 80.16 1211.48 66.13 $$ D = \sqrt[3]{{D_x^2{D_y}}} $$ (1) 本文用到以下无量纲参数: 雷诺数Re = ρUω/μ, 韦伯数We = ρU2ω/σ, 毛细数Ca = μU/σ. 式中, ρ为液滴密度, U为液滴撞网速度, ω为网孔宽度(如图1(a)所示), μ为液体黏性系数, σ为液体表面张力系数, 由Dataphysics OCA20接触角测量仪测定.
被碰撞网面为商用紫铜网, 网面一端使用夹具固定, 另一端利用砝码张紧网面(如图1(b)所示), 使网面保持平整并减少液滴撞网过程中网面振动对液滴破碎的影响. 实验中用到了5种目数Nm的网面, 它们的网孔宽度ω和网线直径d如表2所示. 为考量网面润湿性对实验结果的影响[60-63], 我们使用ImageJ软件测量了2种甘油水液滴在5种网面上的静态接触角θeq, 如表3所示, μ = 5.01 mPa·s和μ = 80.16 mPa·s液滴在网面上的静态接触角最大差值分别为10.3°和9.0°, 该值相对各自的最小静态接触角均小于10%, 因此网面润湿性的影响可以忽略.
表 2 网面参数Table 2. The mesh parametersNm ω/μm d/μm 60 250 160 80 200 140 100 150 125 120 125 90 180 80 56 表 3 甘油水溶液在网面上的静态接触角Table 3. The equilibrium contact angles of aqueous glycerol solutions droplet on the meshNm θeq/(°)
(μ = 5.01 mPa·s)θeq/(°)
(μ = 80.16 mPa·s)60 125.1 128.8 80 123.0 124.2 100 119.5 124.8 120 114.8 126.7 180 125.3 119.8 如图1所示, 一台Photron Nova S12高速相机被用来以20 000帧/秒的时间分辨率和20.6 μm/pixel的空间分辨率拍摄液滴撞网过程. LED光源配合散光板可以为阴影成像提供均匀光照.
预湿网面是指网面上已经预先有液体存在, 普遍存在于实际应用中[34]. 本文预湿网面的制备如图1(c)所示, 使用微量移液器将液体轻放到网面上, 通过调整液体量实现预湿高度H的控制, H使用高速相机记录的液滴撞击前图片测量.
2. 结果与分析
2.1 液滴撞击干网
2.1.1 最大液指长度
液滴撞击网面后, 在网面下产生的液指如果未达到破裂产生二次液滴的临界长度, 就会在到达最大长度Lmax后回缩, 每一次撞击中产生的最长液指都出现在撞击点附近, 如图2 (a)所示.
图2(a) ~ 图2(c)为浓度分别为75% (μ = 27.73 mPa·s), 80% (μ = 45.37 mPa·s), 84% (μ = 71.06 mPa·s)的甘油水液滴以U = 1.97 m/s的速度撞击80目网面后的演化过程. 图中的无量纲时间t = T/[(D/2)/U] = 2TU/D, T为图像对应的实际时刻, 以液滴底部接触到网面或液膜的时刻为零点. 图2左上角标尺和本文其他图中的标尺长度均为2 mm. 液滴撞击网面后形成液指, 到达最大长度后逐渐缩短, 并融合成为整体. 在整个液滴撞网过程中, 网面下形成的液指长度先增大后减小; 图2(a)中75% (μ = 27.73 mPa·s)的甘油水液滴在t = 2.45时刻, 液指达到最大长度; 最大液指长度随黏性增加而减小, 说明液体黏性耗散了液滴撞击动能.
图3(a) ~ 图3(c)为浓度为80% (μ = 45.37 mPa·s)的甘油水液滴以速度U = 1.97 m/s分别撞击60目、80目、100目网面后的演化过程. 液指长度随时间的变化趋势与图2相同; 图3(a)中, t = 3.63时刻, 液指达到最大长度; 最大液指长度随网孔尺寸减小而减小.
图4为液滴撞击网面后形成的最大液指长度Lmax随液滴黏性μ和网目数Nm变化的曲线, 实验中液滴撞击速度U 均为 1.97 m/s. 图中数据是3次实验测得数据的平均值, 误差线表示标准偏差. 需要说明的是, 当液滴黏性较小和/或网孔尺寸较大时, 液滴撞击后形成完全穿透, 液指断裂形成二次液滴而无法测量液指最大长度, 所以没有数据, 如液滴黏性为5.01 mPa·s情况下, 撞击4种网面后均形成了完全穿透. 由图可知, 最大液指长度随液滴黏性的增大而减小, 如在Nm = 100情况下, 液滴黏性从11.54 mPa·s增加到80.16 mPa·s, 最大液指长度Lmax从0.52 mm降低至0.03 mm; 网线直径接近时, 最大液指长度随网孔宽度的增加而增加, 如在液滴黏性μ = 36.94 mPa·s情况下, 网孔宽度从125 μm(120目)增加到200 μm (80目), Lmax从0.10 mm增加至0.70 mm. 实验室空气流动、实验台微振动等随机扰动影响撞击点而产生的随机误差是图4中Nm = 60 和μ = 45.37 mPa·s及Nm = 80和μ = 27.73 mPa·s条件下Lmax标准偏差较大的主要原因, 该误差随网孔减小和液滴黏性增大而减小.
图5为浓度为78% (μ = 36.94 mPa·s)的甘油水液滴分别以U = 1.00, 1.54, 1.98 m/s的速度撞击80目网面后的演化过程. 由图可知, 最大液指长度随撞击速度的增加而增加, U = 1.00 m/s情况下Lmax = 0.12 mm; U = 1.54 m/s情况下Lmax = 0.30 mm; U = 1.98 m/s情况下Lmax = 0.59 mm. 3种撞击速度下, 液指均在t = 0.99时刻达到最大长度.
图6为液滴撞击网面后形成的最大液指长度Lmax随液滴黏性μ和速度U变化的曲线, 图中所有实验均使用80目网面. 与图4类似, 当液滴速度较大时, 撞网形成完全穿透而没有最大液指长度数据. 由图可知, 最大液指长度Lmax随液滴黏性的增大而减小, 如在U = 1.21 m/s情况下, 液滴黏性从5.01 mPa·s增加到80.16 mPa·s, 最大液指长度Lmax从0.85 mm降低至0.08 mm; Lmax随撞击速度增大而增大, 如在液滴黏性μ = 36.94 mPa·s情况下, 液滴速度从0.98 m/s增加到1.97 m/s, Lmax从0.13 mm增加至0.70 mm.
2.1.2 最大液指长度预测模型
Abouelsoud等[58]基于4种黏性的低表面张力液滴撞击1种网面(接近180目)的实验结果, 建立了一种预测液指最大长度的方法, 该方法中的经验系数有4个值, 是分别根据这4种黏性液滴的实验结果获得的, 因此仅能预测这4种黏性液滴撞击这1种网面形成的液指长度, 无法考虑其他黏性和网孔宽度变化对液指长度的影响, 应用受到限制. 使用该方法我们计算了9.35 mPa·s和19 mPa·s液滴撞击180目网形成液指的最大长度, 如图4中2个紫色实心圆所示, 预测与实验结果也有偏差. 为建立一个更为通用的预测模型, 下面我们对液滴撞网时刻和液指长度最大时刻的能量进行分析.
液指达到最大长度时, 所有能量转化为液体的表面能. 在撞击时刻, 初始撞击能量是动能 (ρDω2U2/2)和表面能(2ω2σ)的总和; 在撞击过程中, 重力对液指做正功(ρgω2L2/2), 黏性力对液指做负功(µωLU); 在液指达到最大长度时刻, 表面能为(σπωL + σπω2/4) , 其中g为重力加速度, L为液指的实时长度.
在撞击过程中, 部分液体穿透网面形成液指, 剩余液体在网面上铺展. 驱动液体铺展将消耗部分初始撞击能量, 因此初始撞击能量只有部分被用于驱动液体穿透网面[58], 而前述重力与黏性力做功则全部作用于液指. 定义取值范围为0 ~ 1的系数φ与初始撞击能量相乘, 表示驱动液体穿透网面的初始撞击能量比例, 即φ(ρDω2U2/2 + 2ω2σ). 由能量守恒可得
$$ \begin{split} & \varphi \left( {\rho D{\omega ^2}{U^2}/2 + 2{\omega ^2}\sigma } \right) + \rho g{\omega ^2}{L^2}/2= \\ &\qquad \sigma \text{π} \omega (L + \omega /4) + \mu \omega UL \end{split} $$ (2) 由实验结果可知, 液滴黏性和网孔宽度均影响φ取值, 取φ = φ1φ2 , φ1 ~ μ1, φ2 ~ ω1, 其中μ1 = μ/(1 mPa·s), ω1 = ω/(1 μm)分别为无量纲黏性系数和网孔宽度. 为根据图6实验数据计算各个黏性系数对应的φ1值, 需先确定ω1 = 200时的φ2值, 此处取φ2 = 1 (可取任意非零常数, 不影响最终φ的拟合公式), 由图6数据可得图7所示的φ1值, 根据这些值可拟合计算φ1的经验公式, 如图7右上所示. 根据φ1经验公式计算各个网孔宽度时的φ2值, 如图8中数据点所示, 根据这些值拟合计算φ2的经验公式, 如图8上部所示.
由式(2)可得最大液指长度
$$ \begin{split} & {L_{\max }} = \frac{{ {\mu U + \sigma \text{π} } }}{{\rho g\omega }} - \\ &\qquad \frac{{\sqrt {{{\left( {\mu U + \sigma \text{π} } \right)}^2} - \rho g{\omega ^2}\left( {\varphi \rho D{U^2} + 4\varphi \sigma - \dfrac{{\sigma \text{π} }}{2}} \right)} }}{{\rho g\omega }} \end{split} $$ (3) 由φ = φ1φ2结合φ1和φ2的经验公式得
$$ \varphi = 1.59\mu _1^{ - 0.68}(0.000\;1\omega _1^2 - 0.028{\omega _1} + 2.46) $$ (4) 由式(3)和式(4)可获得各种情况下的最大液指长度预测值, 如图4和图6中的实线所示. 由图可知, 该两式计算值与实验结果无论在数值还是变化趋势上都是较为吻合的, 为合理预测黏性液滴撞击网面形成的最大液指长度提供了一种方法.
2.1.3 黏性、网面尺寸对撞击产生二次液滴的影响
图9(a) ~ 图9(c)为浓度分别为50% (μ = 5.01 mPa·s), 64% (μ = 11.54 mPa·s), 75% (μ = 27.73 mPa·s)的甘油水液滴以U = 1.97 m/s的速度撞击80目网面后的演化过程. 浓度为50% (μ = 5.01 mPa·s)和64% (μ = 11.54 mPa·s)的甘油水液滴(图9(a)和图9(b))撞击网面后, 由于Plateau-Rayleigh不稳定性导致液指在网面下失稳破碎形成了二次液滴. 浓度为75%(μ = 27.73 mPa·s)的甘油水液滴(图9(c))撞击网面后形成液指, 到达最大长度后逐渐缩短, 并融合成为整体. 浓度为50% (μ = 5.01 mPa·s)的甘油水液滴(图9(a))在t = 5.93时刻, 网面下的液指完全破碎形成二次液滴; 从图9(a) ~ 图9(b), 二次液滴数量随黏性增大而减小, 直至图9(c)未出现二次液滴, 形象展示了液体黏性对液滴初始撞击能量的耗散.
图10(a) ~ 图10(c)为浓度为70% (μ = 18.07 mPa·s)的甘油水液滴以速度U = 1.97 m/s分别撞击60目、80目、100目网面后的演化过程. 液滴撞击后随时间的变化情况与图10相同; 图10(a)中t = 5.93时刻, 二次液滴完全形成; 二次液滴数量随网孔宽度减小而减少直至为0, 演示了毛细力在液指形成过程中的阻滞作用.
使用不同的黏性液滴撞击不同结构的干燥网面, 我们获得了临界撞击韦伯数We和临界毛细数Ca, 如图11所示, 临界韦伯数随临界毛细数增加而单调增加. 图中用不同的符号形状表示不同的网面, 不同的符号颜色表示不同的液滴黏性, 误差线表示不确定度, 实线为理论曲线.
2.1.4 产生二次液滴临界条件预测模型
液滴垂直撞击网面时, 最容易穿透产生液指和二次液滴的位置是撞击的中心点, 液滴的投影面积为 ~ πD2/4, 单个网孔结构大小为 ~ (ω + d)2, 因此液滴投影面积下的网格数目N = πD2/[4(ω + d)2], 实验中使用网面中网孔尺寸最大为60目, 最小为180目, 对应N = 34 ~ 265, 因此可以假设液滴在每次撞击时落点都能保持在其中一个网孔的中心而不是撞击在网线上.
前人已经建立了不考虑液滴黏性的临界韦伯数或者临界速度的预测模型. 在无黏情况下, 液滴能否穿透网面由液滴的动压 ~ ρU2和毛细压 ~ σГ/A之间的平衡决定[35,38]. 这里, A = w2和Γ = 4w分别表示网面单个开孔的开口面积和周长. 液体黏性作用可表示为 ~ μU/w [58]. 动压可看作是促进液滴穿透的因素, 而毛细压和黏性作用都在阻碍液滴穿透网面, 可得
$$ \rho {U^2} = m\frac{{\mu U}}{w} + n\frac{{\sigma \varGamma }}{A} = m\frac{{\mu U}}{w} + 4n\frac{\sigma }{w} $$ (5) $$ \frac{{\rho {U^2}\omega }}{\sigma } = m\frac{{\mu U}}{\sigma } + 4n $$ (6) $$ We = mCa + 4n $$ (7) 式中, m和n是两个常数, 需要通过实验确定. 根据图11所示实验结果拟合可得m = 8.22和n = 1.45. 图11中黑色实线即为将这两个常数代入式(7)获得的理论曲线, 由图可知, 理论曲线在数值和变化趋势上均与实验数据吻合较好, 可以用于合理预测黏性液滴撞击干燥网面时产生二次液滴的临界参数.
2.2 液滴撞击预湿网面
2.2.1 预湿对最大液指长度的影响
图12为浓度为70% (μ = 18.07 mPa·s)的甘油水液滴以U = 1.40 m/s的速度分别撞击无量纲预湿高度h = 0, 0.36, 0.98, 1.58的60目网面后的演化过程, 其中h = H/(D/2) = 2H/D. 液滴撞击网面后形成液指, 到达最大长度后逐渐缩短, 并融合成为整体. 从图12(a) ~ 图12(c), 网面下形成的最大液指长度随预湿高度的增加而减小, 直到图12(d)不再出现, 表明网面预湿可完全抑制不完全穿透现象的出现.
图13为液滴撞击预湿网面后形成的最大液指长度Lmax随液滴黏性μ和无量纲预湿高度h的变化, 液滴撞击速度U = 1.4 m/s, 被撞击网面为60目, 图中数据点均为3次实验的平均值, 误差线表示标准偏差, 实线为理论预测值. 与图3类似, 当预湿高度和液滴黏性较小时, 撞网形成完全穿透而没有最大液指长度数据. 由图可知, 最大液指长度随预湿高度的增大而减小, 如在液滴黏性μ = 18.07情况下, 预湿高度从0增加到1.18, 最大液指长度Lmax从0.83 mm降低至0.19 mm; 最大液指长度随液滴黏性的减小而增加, 如在h = 0.47情况下, 液滴黏性μ从80.16 mPa·s减小到5.01 mPa·s, Lmax从0.13 mm增加至2.56 mm.
2.2.2 最大液指长度预测模型
与低黏度液滴撞击预湿网面[34]类似, 黏性液滴撞击预湿网面的演化过程也可以分为如图14所示的两个步骤, 但不同之处是黏性的影响不再能被忽略. 图14(a)显示了液滴以速度U撞击预湿液膜, 此时, 液滴体积为Vd = πD3/6, 假设被撞击的预湿液膜为高度为H, 直径为D的圆柱体, 体积为Vunder = πD2H/4.
虽然液滴撞击液膜过程中局部流动异常复杂, 但是整体来看, 是液体流动的垂直分量决定了液滴穿透的类型, 因此在分析中可以仅考虑流动的垂直分量. 液滴与液滴下方液柱的总体积为Vtotal = Vd + Vunder, 融合后速度为U1, 如图14(b)所示. 撞击过程中液体黏性力为 ~ μDU [58], 作用时间由液膜高度H和液滴速度决定U, 即 ~ H/U, 由动量守恒可得
$$ {U_1} = \frac{{2D/3}}{{2D/3 + H}}U - \frac{{4\mu H{D^2}}}{{(2D/3 + H)\text{π} \rho {D^3}}} $$ (8) 式(8)化简得
$$ {U_1} = \frac{{4U}}{{4 + 3h}} - \frac{{12\mu h}}{{(4 + 3h)\text{π} \rho D}} $$ (9) 此后, 融合的液体可被认为以速度U1撞击干燥网面, 将U1代入式(3)可得黏性液滴撞击预湿网面形成的最大液指长度
$$\begin{split} &{L_{\max }} = \dfrac{ {\dfrac{{4\mu U}}{{4 + 3h}} - \dfrac{{12h{\mu ^2}}}{{(4 + 3h)\text{π} \rho D}} + \sigma \text{π} } }{{\rho g\omega }} - \\ &\qquad \dfrac{{\sqrt {\left[ {\dfrac{{4\mu U}}{{4 + 3h}} - \dfrac{{12h{\mu ^2}}}{{(4 + 3h)\text{π} \rho D}} + \sigma \text{π} } \right]^2 }}}{{\rho g\omega }}-\\ &\qquad \dfrac{\sqrt {\varphi g{\rho ^2}{\omega ^2}D{\left[ {\dfrac{{4U}}{{4 + 3h}} - \dfrac{{12\mu h}}{{(4 + 3h)\text{π} \rho D}}} \right]^2} }}{\rho g\omega}-\\ &\qquad \dfrac{\sqrt { g\rho {\omega ^2}\left( {4\varphi \sigma - \dfrac{{\sigma \text{π} }}{2}} \right)}}{\rho g\omega}\end{split} $$ (10) 将图13的实验条件代入式(10)可得黏性液滴撞击预湿网面后的最大液指长度理论曲线, 如图13中实线所示. 理论曲线的数值和变化趋势均与实验结果吻合良好, 说明式(10)可合理预测黏性液滴撞击预湿网面后的最大液指长度, 也进一步说明φ1和φ2两个拟合参数是合理的.
2.2.3 预湿对产生二次液滴的影响
图15为浓度为50% (μ = 5.01 mPa·s)的甘油水液滴以U = 1.54 m/s的速度分别撞击无量纲预湿高度h = 0, 0.35, 0.64, 1.63的60目网面后的演化过程. 液滴在撞击无量纲预湿高度h = 0, 0.35 (图15(a)和图15(b))的网面后, 液指在网面下破碎形成了二次液滴. 液滴在撞击无量纲预湿高度h = 0.64, 1.63 (图15(c)和图15(d))的网面后形成液指, 到达最大长度后回缩, 未形成二次液滴. 预湿高度为0 (图15(a))情况下, 在t = 4.13时刻, 网面下液指完全破碎形成二次液滴; 从图15(a)到图15(b), 二次液滴数量随预湿高度增大而减小直至图15(c) 未出现二次液滴, 表明网面预湿同样可抑制完全穿透现象.
图16 ~ 图18分别是浓度为50% (μ = 5.01 mPa·s)的甘油水液滴撞击120目、80目和60目网面后是否产生二次液滴的结果相图; 图18 ~ 图20分别是浓度为为50%, 78%和85% (μ = 5.01, 36.94, 80.16 mPa·s)的甘油水液滴撞击60目网面后是否产生二次液滴的结果相图, 纵坐标为撞击时液滴的雷诺数Re, 横坐标为无量纲预湿高度h, h = 0 代表干燥网面. 图中, 蓝色实心圆表示产生二次液滴, 红色实心圆表示不产生二次液滴, 黑色实线是理论预测的临界曲线. 由图可知, 完全穿透所需的雷诺数随预湿高度的增加而增大, 对于同样的撞击条件, 网面预湿可有效抑制完全穿透.
2.2.4 产生二次液滴临界条件预测模型
与最大液指长度预测模型类似, 首先引入液体黏性力的影响, 通过动量守恒获得液滴与预湿液膜融合后的撞击速度为U1 (式(9)); 此后的撞击过程, 可认为是液体以U1撞击干燥网面[34]. 将U1代入式(6)可得临界速度$U_P $
$$ {U_P} = \left(1 + \frac{{3h}}{4}\right)\left( {\frac{{m\mu }}{{2\rho \omega }} + \sqrt {\frac{{{m^2}{\mu ^2}}}{{4{\rho ^2}{\omega ^2}}} + \frac{{4n\sigma }}{{\rho \omega }}} } \right) $$ (11) 将式(11)带入雷诺数式可得临界雷诺数
$$ R{e_T} = \left(1 + \frac{{3h}}{4}\right)\left( {m + \sqrt {{m^2} + \frac{{16\rho \omega n\sigma }}{{{\mu ^2}}}} } \right) + \frac{{6\omega h}}{D} $$ (12) 将图16 ~ 图20实验条件带入式(12)可得产生二次液滴所需的临界雷诺数理论曲线, 如图中黑色实线所示, 理论曲线与实验结果合理吻合, 说明理论模型可以合理预测黏性液滴撞击预湿网面产生二次液滴的所需临界雷诺数.
3. 结 论
本文通过可视化实验研究了10种不同黏性的液滴撞击5种不同结构的干燥和预湿网面形成液指和二次液滴的演化规律, 通过理论分析了相应机理, 并建立了预测模型, 获得主要结论如下:
(1)液滴撞击干燥网面后形成的最大液指长度和液滴完全穿透干燥网面所需的临界速度都随网孔宽度减小、液滴黏性增加而减小;
(2)网面预湿可以完全抑制液滴撞击网面后的不完全穿透和完全穿透, 预湿液膜高度越高抑制效果越明显;
(3)建立了考虑液滴黏性、网孔宽度和网面预湿的液滴撞击网面后形成最大液指长度及出现完全穿透临界参数的理论预测模型, 模型预测结果与实验结果吻合良好.
-
表 1 甘油水液滴的物理属性(25 °C)
Table 1 The physical properties of aqueous glycerol solutions at 25 °C
ψ/wt% μ/(mPa·s) ρ/(kg·m−3) σ/(mN·m−1) 50 5.01 1113.41 68.05 64 11.54 1150.67 66.95 70 18.07 1167.41 66.85 75 27.73 1181.74 66.54 78 36.94 1190.50 66.16 80 45.37 1196.42 66.77 81.5 53.35 1200.90 65.43 83 63.16 1205.41 65.52 84 71.06 1208.43 66.20 85 80.16 1211.48 66.13 表 2 网面参数
Table 2 The mesh parameters
Nm ω/μm d/μm 60 250 160 80 200 140 100 150 125 120 125 90 180 80 56 表 3 甘油水溶液在网面上的静态接触角
Table 3 The equilibrium contact angles of aqueous glycerol solutions droplet on the mesh
Nm θeq/(°)
(μ = 5.01 mPa·s)θeq/(°)
(μ = 80.16 mPa·s)60 125.1 128.8 80 123.0 124.2 100 119.5 124.8 120 114.8 126.7 180 125.3 119.8 -
[1] Brunet P, Lapierre F, Zoueshtiagh F, et al. To grate a liquid into tiny droplets by its impact on a hydrophobic microgrid. Applied Physics Letters, 2009, 95: 254102 doi: 10.1063/1.3275709
[2] Hsu CF, Ashgriz N. Impaction of a droplet on an orifice plate. Physics of Fluids, 2004, 16: 400-411 doi: 10.1063/1.1637036
[3] Soto D, Girard HL, Helloco AL, et al. Droplet fragmentation using a mesh. Physical Review Fluids, 2018, 3: 083602 doi: 10.1103/PhysRevFluids.3.083602
[4] Kooij SA, Moqaddam AM, Goede TCD, et al. Sprays from droplets impacting a mesh. Journal of Fluid Mechanics, 2019, 871: 489-509 doi: 10.1017/jfm.2019.289
[5] Wang L, Wu X, Yu W, et al. Numerical study of droplet fragmentation during impact on mesh screens. Microfluidics and Nanofluidics, 2019, 23: 136
[6] Yarin AL. Novel nanofluidic and microfluidic devices and their applications. Current Opinion in Chemical Engineering, 2020, 29: 17-25 doi: 10.1016/j.coche.2020.02.004
[7] Brewer SA, Willis CR. Structure and oil repellency: Textiles with liquid repellency to hexane. Applied Surface Science, 2008, 254: 6450-6454 doi: 10.1016/j.apsusc.2008.04.053
[8] Kim H, Park Y, Kim H, et al. Critical heat flux enhancement by single layered metal wire mesh with micro and nano-sized pore structures. International Journal of Heat and Mass Transfer, 2017, 115: 439-449 doi: 10.1016/j.ijheatmasstransfer.2017.08.066
[9] Modak CD, Kumar A, Tripathy A, et al. Drop impact printing. Nature Communications, 2020, 11: 4327 doi: 10.1038/s41467-020-18103-6
[10] Lohse D. Fluid mech: Fundamental fluid dynamics challenges in inkjet printing. Annual Review of Fluid Mechanics, 2022, 54: 349-382
[11] He P, Cao J, Ding H, et al. Screen-printing of a highly conductive graphene ink for flexible printed electronics. ACS Applied Materials & Interfaces, 2019, 11: 32225-32234
[12] Al-Dughaither AS, Ibrahim AA, Al-Masry WA. Investigating droplet separation efficiency in wire-mesh mist eliminators in bubble column. Journal of Saudi Chemical Society, 2010, 14: 331-339 doi: 10.1016/j.jscs.2010.04.001
[13] Wang B, Guo ZG. Superhydrophobic copper mesh films with rapid oil/water separation properties by electrochemical deposition inspired from butterfly wing. Applied Physics Letters, 2013, 103: 063704 doi: 10.1063/1.4817922
[14] Dunderdale G, Urata C, Sato T, et al. Continuous, high-speed, and efficient oil/water separation using meshes with antagonistic wetting properties. ACS Applied Materials & Interfaces, 2015, 7(34): 18915-18919
[15] Pi P, Hou K, Zhou C, et al. A novel superhydrophilic-underwater superoleophobic Cu2 S coated copper mesh for efficient oil-water separation. Materials Letters, 2016, 182: 68-71 doi: 10.1016/j.matlet.2016.06.087
[16] Wen R, Xu S, Zhao D, et al. Sustaining enhanced condensation on hierarchical mesh-covered surfaces. National Science Review, 2018, 5: 878-887
[17] Tudu BK, Kumar A. Robust and durable superhydrophobic steel and copper meshes for separation of oil-water emulsions. Progress in Organic Coating, 2019, 133: 316-324.
[18] 张星, 刘金鑫, 张海峰等. 防护口罩用非织造滤料的制备技术与研究现状. 纺织学报, 2020, 41(3): 168-174 (Zhang Xing, Liu Jinxin, Zhang Haifeng, et al. Preparation technology and research status of nonwoven filtration materials for individual protective masks. Journal of Textile Research, 2020, 41(3): 168-174 (in Chinese) Zhang xing, Liu jinxin, Zhang haifeng, et al. Preparation technology and research status of nonwoven filtration materials for individual protective masks. Journal of Textile Research, 2020, 41(03): 168-174(in Chinese)
[19] Bagchi S, Basu S, Chaudhuri S, et al. Penetration and secondary atomization of droplets impacted on wet facemasks. Physical Review Fluids, 2021, 6: 110510 doi: 10.1103/PhysRevFluids.6.110510
[20] Sharma S, Pinto R, Saha A, et al. On secondary atomization and blockage of surrogate cough droplets in single- and multilayer face masks. Science Advances, 2021, 7: eabf0452 doi: 10.1126/sciadv.abf0452
[21] Solano T, Ni C, Mittal R, et al. Perimeter leakage of face masks and its effect on the mask's efficacy. Physics of Fluids, 2022, 34: 051902 doi: 10.1063/5.0086320
[22] Liu Y, Yan X, Wang Z. Droplet dynamics on slippery surfaces: Small droplet, big impact. Biosurface and Biotribology, 2019, 5: 35-45 doi: 10.1049/bsbt.2019.0004
[23] Xu L, Zhang WW, Nagel SR. Drop splashing on a dry smooth surface. Physical Review Letters, 2005, 94: 184505
[24] Yarin AL. Drop impact dynamics: Splashing, spreading, receding, bouncing. Annual Review of Fluid Mechanics, 2006, 38: 159-192
[25] Thoroddsen ST, Etoh TG, Takehara K. High-speed imaging of drops and bubbles. Annual Review of Fluid Mechanics, 2008, 40: 257-285
[26] Riboux G, Gordillo JM. Experiments of drops impacting a smooth solid surface: A model of the critical impact speed for drop splashing. Physical Review Letters, 2014, 113: 024507
[27] Liu Y, Moevius L, Xu X, et al. Pancake bouncing on superhydrophobic surfaces. Nature Physics, 2014, 10: 515-519 doi: 10.1038/nphys2980
[28] Laan N, de Bruin KG, Bartolo D, et al. Maximum diameter of impacting liquid droplets. Physical Review Applied, 2014, 2: 044018 doi: 10.1103/PhysRevApplied.2.044018
[29] Josserand C, Thoroddsen ST. Drop impact on a solid surface. Annual Review of Fluid Mechanics, 2016, 48: 365-391 doi: 10.1146/annurev-fluid-122414-034401
[30] Liang G, Mudawar I. Review of drop impact on heated walls. International Journal of Heat and Mass Transfer, 2017, 106: 103-126 doi: 10.1016/j.ijheatmasstransfer.2016.10.031
[31] Hao J, Lu J, Lee L, et al. Droplet splashing on an inclined surface. Physical Review Letters, 2019, 122: 054501
[32] 孙姣, 周维, 蔡润泽等. 垂直壁面附近上升单气泡的弹跳动力学研究. 力学学报, 2020, 52(1): 1-11 (Sun Jiao, Zhou Wei, Cai Runze, et al. The bounce dynamics of a rising single bubble near a vertical wall. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(1): 1-11 (in Chinese) doi: 10.6052/0459-1879-19-228 Sun Jiao, Zhou Wei, Cai Runze, et al. The bounce dynamics of a rising single bubble near a vertical wall. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(1): 1-11 (in Chinese) doi: 10.6052/0459-1879-19-228
[33] 万其文, 陈效鹏, 胡海豹等. 中性润湿平板上液膜的惯性收缩. 力学学报, 2022, 54(6): 1516-1522 (Wan Qiwen, Chen Xiaopeng, Hu Haibao, et al. Inertial retraction of liquid film on moderately wettable plate. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(6): 1516-1522 (in Chinese) Wan Qiwen, Chen Xiaopeng, Hu Haibao, et al. Inertial retraction of liquid film on moderately wettable plate. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(6): 1516-1522 (in Chinese)
[34] Xu L, Ji W, Lu J, et al. Droplet impact on a prewetted mesh. Physical Review Fluids, 2021, 6: L101602 doi: 10.1103/PhysRevFluids.6.L101602
[35] Ryu S, Sen P, Nam Y, et al. Water penetration through a superhydrophobic mesh during a drop impact. Physical Review Letters, 2017, 118: 014501 doi: 10.1103/PhysRevLett.118.014501
[36] Sen U, Roy T, Chatterjee S, et al. Post-impact behavior of a droplet impacting on a permeable metal mesh with a sharp wettability step. Langmuir, 2019, 35: 12711-12721 doi: 10.1021/acs.langmuir.9b02486
[37] Zhang G, Quetzeri-Santiago MA, Stone CA, et al. Droplet impact dynamics on textiles. Soft Matter, 2018, 14: 8182-8190 doi: 10.1039/C8SM01082J
[38] Lorenceau É, Quéré D. Drops impacting a sieve. Journal of Colloid and Interface Science, 2003, 263: 244-249 doi: 10.1016/S0021-9797(03)00126-7
[39] Sahu RP, Sinha-Ray S, Yarin A, et al. Drop impacts on electrospun nanofiber membranes. Soft Matter, 2012, 8: 3957-3970 doi: 10.1039/c2sm06744g
[40] Boscariol C, Chandra S, Sarker D, et al. Drop impact onto attached metallic meshes: liquid penetration and spreading. Experiments in Fluids, 2018, 59: 1-13 doi: 10.1007/s00348-017-2450-7
[41] Lembach AN, Tan HB, Roisman IV, et al. Drop impact, spreading, splashing, and penetration into electrospun nanofiber mats. Langmuir, 2010, 26: 9516-9523 doi: 10.1021/la100031d
[42] Sahu R, Sett S, Yarin A, et al. Impact of aqueous suspension drops onto non-wettable porous membranes: Hydrodynamic focusing and penetration of nanoparticles. Colloids and Surfaces A, 2015, 467: 31-45 doi: 10.1016/j.colsurfa.2014.11.023
[43] Bae C, Oh S, Han J, et al. Water penetration dynamics through a Janus mesh during drop impact. Soft Matter, 2020, 16: 6072-6081 doi: 10.1039/D0SM00567C
[44] Sun L, Lin S, Pang B, et al. Water sprays formed by impinging millimeter-sized droplets on superhydrophobic meshes. Physics of Fluids, 2021, 33(9): 092111 doi: 10.1063/5.0058512
[45] An T, Cho SJ, Choi W, et al. Preparation of stable superhydrophobic mesh with a biomimetic hierarchical Structure. Soft Matter, 2011, 7: 9867-9870 doi: 10.1039/c1sm06238g
[46] Kumar A, Tripathy A, Modak CD, et al. Designing assembly of meshes having diverse wettability for reducing liquid ejection at terminal velocity droplet impact. Journal of Microelectromechanical Systems, 2018, 27: 866-873 doi: 10.1109/JMEMS.2018.2850903
[47] Su M, Luo Y, Chu G, et al. Dispersion behaviors of droplet impacting on wire mesh and process intensification by surface micro/nano-structure. Chemical Engineering Science, 2020, 219: 115593 doi: 10.1016/j.ces.2020.115593
[48] Jamali M, Vahedi TH, Pourdeyhimi B, et al. Penetration of liquid droplets into hydrophobic fibrous materials under enhanced gravity. Journal of Applied Physics, 2019, 125: 145304
[49] Kumar A, Tripathy A, Nam Y, et al. Effect of geometrical parameters on rebound of impacting droplets on leaky superhydrophobic meshes. Soft Matter, 2018, 14: 1571-1580 doi: 10.1039/C7SM02145C
[50] de Goede TC, Moqaddam AM, Limpens K, et al. Droplet impact of Newtonian fluids and blood on simple fabrics: Effect of fabric pore size and underlying substrate. Physics of Fluids, 2021, 33: 033308 doi: 10.1063/5.0037123
[51] Tang Y, Su M, Chu G, et al. Impact phenomena of liquid droplet passing through stainless steel wire mesh units. Chemical Engineering Science, 2019, 198: 144-154 doi: 10.1016/j.ces.2018.12.035
[52] Xu J, Xie J, He X, et al. Water drop impacts on a single-layer of mesh screen membrane: Effect of water hammer pressure and advancing contact angles. Experimental Thermal and Fluid Science, 2017, 82: 83-93 doi: 10.1016/j.expthermflusci.2016.11.006
[53] Xu L, Zong S, Hao J, et al. Droplet penetration through an inclined Mesh. Physics of Fluids, 2022, 34: 122105 doi: 10.1063/5.0126982
[54] Blackwell BC, Nadhan AE, Ewoldt RH, et al. Impacts of yield-stress fluid drops on permeable mesh substrates. Journal of Non-Newtonian Fluid Mechanics, 2016, 238: 107-114 doi: 10.1016/j.jnnfm.2016.06.012
[55] Mehrizi A, Lin S, Sun L, et al. Spectacular behavior of a viscoelastic droplet impinging on a superhydrophobic mesh. Langmuir, 2022, 38: 6106-6115 doi: 10.1021/acs.langmuir.2c00385
[56] Wang G, Gao J, Luo KH. Droplet impacting a superhydrophobic mesh array: Effect of liquid properties. Physical Review Fluids, 2020, 5: 123605 doi: 10.1103/PhysRevFluids.5.123605
[57] Vontas K, Boscariol C, Andredaki M, et al. Droplet impact on suspended metallic meshes: effects of wettability, Reynolds and Weber numbers. Fluids, 2020, 5: 81 doi: 10.3390/fluids5020081
[58] Abouelsoud M, Kherbeche A, Thoraval, MJ. Drop impact on a mesh—Viscosity effect. Journal of Colloid and Interface Science, 2023, 648: 37-45 doi: 10.1016/j.jcis.2023.04.099
[59] Cheng N. Formula for the viscosity of a glycerol-water mixture. Chemical Engineering Research and Design, 2008, 47: 3285-3288 doi: 10.1021/ie071349z
[60] Zang D, Wang X, Geng X, et al. Impact dynamics of droplets with silica nanoparticles and polymer additives. Soft Matter, 2013, 9: 394-400 doi: 10.1039/C2SM26759D
[61] Lin K, Zang D, Geng X, et al. Revisiting the effect of hierarchical structure on the super hydrophobicity. The European Physical Journal E, 2016, 39: 15 doi: 10.1140/epje/i2016-16015-8
[62] Lin K, Zang D, Li X, et al. Superhydrophobic polytetrafluoroethylene surfaces by spray coating on porous and continuous substrates. The Royal Society of Chemistry, 2016, 6: 47096-47100
[63] Eral HB, ’tMannetje DJCM, Oh JM. Contact angle hysteresis: a review of fundamentals and applications. Colloid and Polymer Science, 2013, 291: 247-260 doi: 10.1007/s00396-012-2796-6
-
期刊类型引用(2)
1. 刘永洪,王意乐,代明璐,丁怡,靳亚康,陈龙泉. 黏弹性液滴撞击超疏水表面的反弹行为研究. 力学学报. 2024(06): 1552-1562 . 本站查看
2. 陈苗苗,王兆长,郑楠,张国涛,涂德浴,童宝宏. 微量润滑油滴冲击倾斜壁面黏附颗粒的实验研究. 力学学报. 2024(10): 2854-2864 . 本站查看
其他类型引用(1)