液滴撞击液膜的流热耦合界面追踪方法数值模拟

采用界面追踪方法研究蒸馏过程中液滴撞击高温薄液膜的流动和传热特性，将数值结果与解析解和实验进行比较验证模型的正确性，研究气液界面和热流分布的演变过程.同时，分析液滴We数和无量纲液膜厚度对传热的影响.液滴撞击后的热流密度分布显示：液膜可分为撞击区、过渡区和静态区.由于液滴的撞击作用，强制对流是撞击区内主要的传热机制.增大液滴的韦伯数或减小无量纲液膜厚度会加强热量传递.随着液滴韦伯数的增加，冲击引起的扰动增强，在动量和能量共同作用下，平均热流密度明显增大，撞击区冠状水花越明显.无量纲液膜厚度越小，平均热流密度越大，且有更长的时间保持高热流密度换热.     

双液滴撞击平面液膜的流动与传热特性

采用耦合的水平集-体积分数法(CLSVOF)对双液滴连续撞击恒定壁温壁面上的热液膜的流动和换热特性进行了数值模拟及分析,得到了双液滴撞击热液膜后形态演变的过程.分析了液滴垂直间距、撞击速度、液膜厚度以及液滴直径对双液滴撞击液膜后的流动与传热特性的影响,结果显示,壁面平均热流密度随液滴撞击速度的增大而增大,液滴垂直间距、液膜厚度和液滴直径对平均热流密度的影响较小,但会对热流密度在撞击区域和交界区的分布产生重要影响.     

双液滴同时垂直撞击壁面的数值研究

采用质量守恒的level set方法对双液滴同时垂直撞击干壁面后的流动过程进行了模拟研究,主要讨论了韦伯数(We)、壁面接触角(θ)以及双液滴水平间距(S)等物理参数对相界面流动过程的影响,分析了不同参数下射流高度和水平铺展半长随时间的变化规律.研究表明:We数较大时,中心射流液柱将产生二次液滴,随后液柱反弹至空中,且We数越大,中心射流产生的二次液滴次数越多,最大无量纲射流高度和最大无量纲铺展半长越大;随壁面接触角的增大,中心射流液柱出现反弹现象,水平铺展液流出现断裂的时间越早,最大无量纲射流高度和最大无量纲铺展半长越小;最大无量纲射流高度值与液滴水平间距的相关性不单调,铺展半长随水平间距的增大而增大.     

实验观测液滴撞击倾斜表面液膜的特殊现象

采用高速摄像仪以10000帧/s 的拍摄速度对液滴撞击倾斜表面液膜的过程进行了实验观测, 分析了液滴撞击倾斜表面液膜后的铺展、水花形成以及飞溅等现象, 考察了撞击角对液滴震荡变形过程的影响; 在此基础上, 定量讨论了液滴铺展速度随时间的变化规律, 揭示了液滴撞击速度和撞击角对前、后铺展因子及初始铺展速度的影响.观测发现, 在撞击角为28.0°–74.7°范围内, 随着撞击角的减小, 液滴在液膜表面的震荡变形程度增大; 前铺展因子随撞击速度的增大而增大, 随撞击角的减小而增大; 后铺展因子随撞击速度的增大几乎不发生变化, 但是随撞击角的增大而增大; 液滴初始铺展速度随撞击速度和撞击角的升高而增大. 关键词： 液滴撞击 倾斜液膜 铺展因子 铺展速度     

单液滴撞击液膜的颈部射流模拟及机理分析

本文采用DNS耦合ACLS方法对喷雾中常见的单液滴撞击水平液膜现象进行二维的数值模拟,通过撞击后形态演变及内部场信息对比,重点分析颈部射流的形成机理,随后将研究对象拓展到倾斜液膜。本文验证了ACLS方法的准确性,成功捕捉到液膜倾角对空气卷吸现象的促进作用。研究表明,局部压差是液滴撞击水平液膜后产生颈部射流的直接原因,而对于撞击倾斜液膜的案例,前后两端颈部射流成因不同:前端射流由局部压差产生而后端射流随着液膜倾斜角的增大,旋涡机理的影响比重逐渐增强,两者处于共存竞争关系。     

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

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

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

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

液滴撞击液膜的射流与水花形成机理分析

建立了单液滴撞击平面液膜的物理与数学模型,采用Coupled Level Set and Volume of Fluid方法对这种现象进行了数值模拟,探讨了黏度和表面张力对冠状水花形态的影响.通过分析撞击后液体内部的压力和速度分布,揭示了液滴颈部射流的产生机理,验证了Yarin和Weiss提出的运动间断理论.研究显示,表面张力对冠状水花形态的影响远大于黏度的影响.颈部射流的产生主要是由于撞击后颈部区域局部较大压差造成的,随着撞击过程的继续,压差作用减弱;液膜内流体的径向运动对射流发展成冠状水花具有推动作用.     

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

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

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

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

液体横向射流在气膜作用下的破碎过程

为了研究液体横向射流在气膜作用下的破碎过程,采用背景光成像技术及VOF TO DPM方法进行了实验研究和仿真研究,模拟介质为水和空气.研究结果表明,液体射流在气膜作用下主要存在两种破碎过程:柱状破碎和表面破碎.Rayleigh-Taylor(R-T)不稳定性产生的表面波是液体射流发生柱状破碎的主要原因,气流穿透表面波的波谷导致射流柱破碎,破碎后的液丝沿流向逐渐发展呈带状分布.Kelvin-Helmholtz(K-H)不稳定性产生的表面波是液体射流发生表面破碎的主要原因,液丝和液滴从射流表面剥离.局部动量比对液体横向射流的破碎过程具有重要影响,当局部动量比较低时,液体射流的破碎由K-H不稳定性主导;随着局部动量比的增大液体射流的破碎逐渐由R-T不稳定性主导.液体射流的破碎长度及穿透深度均随局部动量比的增大而增大.     

格子Boltzmann方法模拟双液滴撞击液膜的流动过程

采用单相自由面格子Boltzmann方法中流场计算,自由面条件和表面张力处理方法,研究双液滴撞击液膜的流动过程,实现对垂直间隔的两液滴相继撞击液膜流动过程的模拟.数值结果表明,液滴相继冲击液膜时,液滴间上下距离对液膜状态有很大的影响,液滴间距不同,液膜形状会出现很大差异.数值结果与实验结论定性一致.     

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.     

Air ingestion by a buckled viscous jet of silicone oil impacting the free surface of the same liquid

As frequently observed in common life, a jet of a viscous liquid impacting on a horizontal surface does not remain straight but instead buckles and folds periodically. We report experiments with planar (ribbonlike) jets of silicone oil impacting the free surface of the same liquid and describe the way in which jet folds incorporate air. It is shown that air ingestion proceeds through different modes, each of them acting as a source of monodisperse bubbles and featuring a threshold in jet height. These sources result from the breakup of remarkable cuspidal structures, produced by the recession of air domains within liquid folds.     

液滴撞击加热壁面传热实验研究

本文采用高速摄像仪对水滴和乙醇液滴撞击加热壁面后的蒸发过程进行了实验观测, 分析了液滴撞击加热壁面后的蒸发特性参数. 实验中, 两种液体初始温度均为20 ℃, 不锈钢壁面初始温度范围为68-126℃. 水滴初始直径为2.07 mm, 撞击壁面时Weber 数为2-44; 乙醇液滴初始直径为1.64 mm, Weber数为3-88. 结果表明, 液滴受到重力、表面张力及流动性的影响, 在蒸发过程的大部分时间内, 水滴高度持续降低而接触直径几乎不变; 蒸发后期, 液滴发生回缩, 水滴的接触直径、高度和接触角出现振荡现象. 乙醇液滴的接触角随时间的增加呈现先减小随后保持不变的趋势, 而接触直径和高度则持续减小, 直到液滴完全蒸发. 液滴蒸发总时长与液体物性和壁面温度有关, 随壁面温度的升高而减小, 与液滴撞击壁面时的Weber 数无关. 同时, 随着壁面温度的升高, 液滴显热部分占总换热量的比重增大, 显热部分能量不可忽略, 本文实验条件下得到水滴的平均热流密度为0.014-0.110 W·mm^-2.     

高炉渣熔融液滴撞击壁面行为特性实验研究

针对高温熔渣粒化技术的开发,本文通过可视化实验研究了高炉渣熔融液滴撞击不锈钢壁面的动态行为特性。结果表明:随液滴雷诺数增大,撞击过程液滴形态演变模式依次由铺展-回缩转变为铺展-回缩破碎和铺展-破碎-凝固;增大壁面粗糙度可减弱与壁面换热,抑制液滴铺展;减小壁面粗糙度促进液滴发生破碎;液膜回缩过程普遍出现回卷现象,壁面倾角越大,液膜回卷和液滴滚动现象越显著,且液滴铺展面积越大,在壁面停留时间越长;减小液滴雷诺数、减小壁面粗糙度并采用垂直粒化仓壁面有利于壁面防黏结。     

Simulation of impact of a hollow droplet on a flat surface

Despite many theoretical and experimental works dealing with the impact of dense continuous liquid droplets on a flat surface, the dynamics of the impact of hollow liquid droplets is not well addressed. In an effort to understand dynamics of the hollow droplet impingement, a numerical study for the impact of a hollow droplet on a flat surface is presented. The impingement model considers the transient flow dynamics during impact and spreading of the droplet using the volume of fluid surface tracking method (VOF) coupled with the momentum transport model within a one-domain continuum formulation. The model is used to simulate the hydrodynamic behaviour of the impact of glycerin hollow droplet. It is found that the impact and spreading of the hollow droplet on a flat surface is distinctly different from the conventional dense droplet and has some new hydrodynamic features. A phenomenon of formation of a central counter jet of the liquid is predicted. With the help of simulations the cause of this phenomenon is discussed. Comparison of the predicted length of the central counter jet and the velocity of the counter jet front shows good agreements with the experimental data. The influence of the droplet initial impact velocity and the hollow droplet shell thickness on the impact behaviour is highlighted.     

A multicomponent multiphase lattice Boltzmann model with large liquid–gas density ratios for simulations of wetting phenomena

A multicomponent multiphase(MCMP) pseudopotential lattice Boltzmann(LB) model with large liquid–gas density ratios is proposed for simulating the wetting phenomena. In the proposed model, two layers of neighboring nodes are adopted to calculate the fluid–fluid cohesion force with higher isotropy order. In addition, the different-time-step method is employed to calculate the processes of particle propagation and collision for the two fluid components with a large pseudoparticle mass contrast. It is found that the spurious current is remarkably reduced by employing the higher isotropy order calculation of the fluid–fluid cohesion force. The maximum spurious current appearing at the phase interfaces is evidently influenced by the magnitudes of fluid–fluid and fluid–solid interaction strengths, but weakly affected by the time step ratio.The density ratio analyses show that the liquid–gas density ratio is dependent on both the fluid–fluid interaction strength and the time step ratio. For the liquid–gas flow simulations without solid phase, the maximum liquid–gas density ratio achieved by the present model is higher than 1000:1. However, the obtainable maximum liquid–gas density ratio in the solid–liquid–gas system is lower. Wetting phenomena of droplets contacting smooth/rough solid surfaces and the dynamic process of liquid movement in a capillary tube are simulated to validate the proposed model in different solid–liquid–gas coexisting systems. It is shown that the simulated intrinsic contact angles of droplets on smooth surfaces are in good agreement with those predicted by the constructed LB formula that is related to Young's equation. The apparent contact angles of droplets on rough surfaces compare reasonably well with the predictions of Cassie's law. For the simulation of liquid movement in a capillary tube, the linear relation between the liquid–gas interface position and simulation time is observed, which is identical to the analytical prediction. The simulation results regarding the wetting phenomena of droplets on smooth/rough surfaces and the dynamic process of liquid movement in the capillary tube demonstrate the quantitative capability of the proposed model.