液滴低速撞击润湿球面现象观测分析

采用高速摄像仪以10000 帧/s的拍摄速度对液滴低速撞击润湿球体表面过程进行了实验观测, 分析了液滴撞击后的反弹、局部反弹和铺展等现象, 考察了黏度对撞击过程的影响; 在此基础上, 定量讨论了液滴铺展特征参数随撞击速度、球体直径和黏度的变化规律. 观测发现: 黏度较大且撞击速度较低时, 撞击后可能出现反弹和局部反弹, 黏度较小时则不发生; 铺展面积随撞击速度的增大而增大; 黏度增大时, 铺展因子减小; 在球体直径为4–20 mm范围内, 随着球体直径的增加, 铺展因子呈上升趋势. 关键词： 液滴撞击 润湿球面 铺展 黏度     

单液滴正碰球面动态行为特性实验研究

在考虑空气阻力影响,确定液滴撞击球面速度的基础上,对较高韦伯数液滴撞击干燥球面动态行为过程进行了实验研究,分析了球面曲率与韦伯数对液滴撞击行为和铺展因子的影响,并与前人撞击平面结果进行了对比.实验表明,靠近撞击球面时,液滴降落速度出现明显波动;球面曲率对液滴撞击后行为影响明显,曲率较大时,液滴撞击后铺展液膜会超出球面直径并滑落,曲率较小时,液滴撞击后在球面上呈现明显的铺展、回缩、震荡、着附动态变化行为,此时最大铺展因子受曲率影响小,随曲率减小,逐渐趋向于撞击平面时的最大铺展因子;韦伯数对液膜铺展速率影响较小,但对液膜回缩时间影响明显,最大铺展因子随韦伯数增加逐渐增大,获得的关联式呈指数变化.     

液滴撞击移动表面动力学特性

试验观测了水滴撞击移动水平铝表面动力学行为,分析了液滴撞击移动表面铺展特征,探讨了液滴不同尾流形态的形成机理。试验研究发现,液滴撞击移动表面存在铺展及拉伸过程。表面剪切力对液滴动力学行为起决定作用,铺展阶段抑制液膜前端液指形成,拉伸阶段促进后端液指产生。表面剪切力对液滴纵向铺展因子影响较大,横向最大铺展因子随剪切韦伯数增大而减小。液滴拉伸后尾流可分为黏附沉积、液指断裂、单液指及多液指四种形态。     

液滴撞击固体表面铺展特性的实验研究

利用高速摄影仪记录了不同液滴撞击固体表面的形态变化过程, 并探讨了液滴撞击参数对撞击过程液滴形态的影响. 结果表明: 液滴黏度对液滴铺展过程起着决定性作用, 液滴表面张力对液滴铺展后的回缩起到主导作用, 两者的共同作用决定着液滴的震荡特性; 撞击速度的提高会增大液滴的最大铺展因子, 但达到最大铺展的时间因表面张力的不同呈现不同变化规律.     

盐水液滴撞击固体壁面接触特性实验研究

采用高速摄像仪对不同浓度盐水液滴撞击水平固体表面的现象进行了观测并记录,总结了接触角、铺展系数和无量纲高度的变化规律,并对其影响参数进行了分析。结果表明:在液滴铺展阶段,盐水液滴的动态接触角大于纯水液滴的,在回缩阶段,其动态接触角反而小于纯水液滴的;液滴铺展系数随盐浓度的增加而减小;液滴直径增大,其最大铺展系数增大,铺展系数变化的震荡频率减弱;对于较大直径液滴,在铺展阶段,盐浓度对液滴动态铺展特性影响较小;随液滴撞击速度增加,液滴的铺展系数增大,而回缩到最小铺展直径所用时间缩短,液滴震荡越易受到抑制;相对于速度而言,液滴直径大小更能影响液滴的震荡频率。     

小韦伯数下液滴撞击光滑壁面的数值模拟

采用界面追踪法数值模拟了液滴以小We撞击光滑壁面的动力学过程.通过数值模拟发现,当液滴初始动能较小时(We=12),液滴表面能主导了能量转换过程,因此最大铺展率βmax随Re增加而单调增加;当液滴初始动能继续增加到中等程度时(We=30),能量转换主导因素存在液滴黏性与液滴变形诱导强应变的竞争机制,最终导致βmax随Re增加而非单调变化,而液滴残余动能随Re增加而增大也进一步促进了非单调βmax;随着液滴初始动能的进一步增加(We=60),动能主导了液滴铺展的能量转换过程,最大铺展率βmax再次随Re增加而单调增加.当液滴处于"无气膜"铺展时,液滴底部表面张力的消失导致液滴可有更大的βmax.     

剪切变稀液滴撞击不同浸润性壁面的数值模拟研究

基于有限元法,采用水平集方法捕捉相界面的移动,构建了液滴撞击固体壁面的数值模型.通过修正的幂律模型描述流体的非牛顿剪切变稀特性,探讨了剪切变稀特性对液滴撞击固体壁面后铺展行为的影响,分析了撞击不同浸润性壁面时剪切变稀特性对液滴撞击壁面行为的影响差异.研究结果表明:随着幂律指数m的减小,液滴撞击过程中的黏性耗散减小,液滴的形貌变化及无量纲参数变化更为显著.接触角为55°的情况下:当m降低至0.85时,液滴铺展过程中开始出现显著区别于牛顿流体液滴的振荡现象;当m降低至0.80时,液滴在回缩过程中会出现中心液膜断裂的情况.接触角为100°时,剪切变稀液滴均会出现振荡行为,振荡幅度随着m的减小而增大.接触角为160°时,牛顿流体液滴与剪切变稀液滴均会在回缩过程中弹起,但剪切变稀液滴的弹起速度更快.此外,基于数值计算结果,本文提出了接触角为55°情况下剪切变稀液滴撞击壁面后的最大无量纲铺展直径预测模型.     

液滴撞击微尺度矩形沟槽表面的数值研究

采用CLSVOF方法建立了液滴撞击微尺度矩形沟槽表面的三维数值模型,对撞击过程中的动态特性及传热特性进行了数值研究。分析了液滴在微尺度矩形沟槽表面的润湿状态,给出了液滴润湿状态转变的临界速度。研究了表面接触角及撞击速度对液滴铺展特性的影响。数值研究结果表明:微尺度沟梢结构会使液滴在撞击后产生横纵差异。液滴的最大铺展因子随着撞击表面接触角的增大而减小,随着撞击速度的增大而增大。液滴撞击微尺度矩形沟槽表面后,撞击表面的热流量先增大后减小。撞击表面的最大热流量受到表面接触角与沟槽深度的耦合作用,撞击表面的最大热流量随着表面接触角的增大而减小。     

液滴撞击超亲水表面冷液膜的研究

液滴撞击过程因具有较强的传热传质性能被广泛应用于工业领域中。本文利用高速摄影机和红外热像仪,研究了液滴撞击超亲水表面冷液膜的水力学特性和温度分布,探讨了撞击We数和液膜温度对撞击过程中水力学特性以及液膜温度分布的影响.总结了液滴在超亲水表面的液膜上的水力学特征的变化规律。实验结果表明,液滴铺展速度和最大铺展直径随撞击We数的增大而增大。同时,液膜温度会影响液滴撞击薄液膜后的水力学特性.在低We数下液滴撞击低温薄液膜后液膜的温度呈高低相间的环状分布,随着撞击We的增大,该环状温度分布消失。这对要求精确喷雾控温的工业过程起到了十分重要的意义.     

含活性剂液滴在倾斜粗糙壁面上的铺展稳定性

针对含非溶性活性剂液滴在倾斜粗糙壁面上的铺展过程, 应用润滑理论推导出基态和扰动态下液膜厚度和活性剂浓度的演化方程组, 基于非模态理论研究了液滴铺展的稳定性特征, 探讨了相关参数的影响及其内在机理. 研究表明: 液膜厚度和活性剂浓度扰动量均呈现双驼峰型变化, 且峰值位于液滴底部凹陷处; 随扰动波数k增加, 最易失稳区域由液滴底部右侧凹陷处移至左侧凹陷, 壁面结构的影响逐渐减弱, 液滴铺展历程趋于稳定; 增加壁面倾角θ 将导致液滴铺展不稳定性加剧, 增大壁面高度D和壁面波数k₀均导致液滴铺展稳定性先增强后减弱; 随毛细数C减小, 液滴铺展稳定性下降, 重力的影响逐渐突显, 扰动量最大值呈现先增大后减小的变化趋势.     

单液滴撞击倾斜液膜飞溅过程的耦合Level Set-VOF模拟

采用耦合水平集--体积分数法(CLSVOF)对液滴撞击倾斜表面液膜后液膜的形态演化及飞溅过程进行数值模拟, 并对液滴撞击液膜过程中形成的空气卷吸现象进行研究并探讨了撞击角对此的影响, 分析了液滴撞击后液体内部的压力和速度分布, 对液滴撞击倾斜表面液膜的飞溅过程进行讨论, 并与实验结果进行了对比, 验证了CLSVOF方法研究液滴撞击倾斜液膜的可行性. 结果表明, 液滴撞击倾斜液膜时前后两部分飞溅现象产生的机理不同, 前半部分飞溅是由于压差引起的颈部射流, 而后半部分则是由液膜径向流动产生的飞溅现象. 随着撞击角的增大, 空气卷吸气泡数量减少.     

液滴撞击圆柱内表面的数值研究

针对液滴撞击圆柱内表面的过程,利用基于相场的格子Boltzmann方法模拟液滴以不同初速度、从不同初始高度、撞击不同大小的圆柱内表面时液滴的形态变化,分析了液滴自身物性(如密度和黏性等)和圆柱内表面润湿性等因素对撞击现象的具体影响.研究发现:撞击韦伯数、密度比及动力黏性比、圆柱半径等对液滴撞击后沿圆柱内表面的铺展均有一定影响,较高的韦伯数下液滴可能会发生分裂;液滴初始高度对大密度比和动力黏性比的撞击影响较小;液滴反弹现象可能出现在接触角较大时;重力作用会抑制撞击后液滴的振荡.     

液滴振荡模型及与数值模拟的对比

由于碰撞壁面后液滴内部流动的复杂性, 以及气-液-固三相间的相互作用, 对液滴碰撞壁面形态变化的数学理论研究有较大的难度, 因此所见者多为实验和数值模拟. 本文通过对液滴受力状态的分析, 得到了惯性力、黏性力和表面张力带经验系数的表达式, 并进一步建立了液滴碰撞壁面振荡模型, 得到了液滴铺展半径的振荡表达式, 以及表面张力、黏性系数等参数对液滴铺展的影响. 最后通过与液滴衰减振荡数值模拟结果的对比, 确定了液滴振荡模型中的修正系数, 验证了模型的可行性. 关键词： 液滴碰撞 振荡 铺展半径/高度 数值模拟     

倾斜粗糙壁面上含不溶性活性剂溶液的动力学特性

针对倾斜随机粗糙壁面上含不溶性活性剂溶液的流动过程, 采用润滑理论建立了液膜厚度和浓度的时空演化模型, 通过PDECOL程序数值求解得到了液膜流/液滴铺展的动力学特性及壁面结构参数的影响. 研究表明: 在重力分量和Marangoni效应共同作用下, 液膜流/液滴铺展速度加快, 液膜边缘和液滴中心出现毛细隆起, 液膜/液滴底部出现凹陷, 同时受粗糙壁面影响, 液膜表面变形更显著. 增加壁面倾角θ具有使重力分量和Marangoni效应增强, 导致隆起和凹陷程度均有所增加的作用. 增大壁面高度D可使液膜流/液滴铺展速度加快, 表面变形放大. 而壁面波数k₀则使液膜流/液滴铺展过程减缓, 抑制隆起和凹陷产生. 与液膜流相比, D和k₀对液滴铺展速度的影响相对较小. 关键词： 随机粗糙壁面 液膜 Marangoni效应 倾斜流动     

Numerical simulation of two droplets impacting upon a dynamic liquid film

The impact of droplets on the liquid film is widely involved in industrial and agricultural fields. In recent years, plenty of works are limited to dry walls or stationary liquid films, and the research of multi-droplet impact dynamic films is not sufficient. Based on this, this paper employs a coupled level set and volume of fluid (CLSVOF) method to numerically simulate two-droplet impingement on a dynamic liquid film. In our work, the dynamic film thickness, horizontal central distance between the droplets, droplets' initial impact speed, and simultaneously the flow velocity of the moving film are analyzed. The evolution phenomenon and mechanism caused by the collision are analyzed in detail. We find that within a certain period of time, the droplet spacing does not affect the peripheral crown height; when the droplet spacing decreases or the initial impact velocity increases, the height of the peripheral crown increases at the beginning, and then, because the crown splashed under Rayleigh-Plateau instability, this results in the reduction of the crown height. At the same time, it is found that when the initial impact velocity increases, the angle between the upstream peripheral jet and the dynamic film becomes larger. The more obvious the horizontal movement characteristics, the more restrained the crown height; the spread length increases with the increase of the dynamic film speed, droplet spacing and the initial impact velocity. When the liquid film is thicker, more fluid enters the crown, due to the crown being unstable, the surface tension is not enough to overcome the weight of the rim at the end of the crown, resulting in droplets falling off.