不同润湿性表面液滴铺展及蒸发过程的LBM数值模拟

基于格子玻尔兹曼方法（Lattice Boltzmann Method, LBM）对固着在加热基板上的液滴铺展及蒸发过程进行模拟,主要研究重力场、基板润湿性以及初始环境温度对液滴铺展及蒸发过程的影响。通过预测蒸发过程中液滴与基板的接触直径变化和液滴剩余质量变化,分析液滴形状及体积变化。研究结果发现,液滴形貌及蒸发过程受重力影响较大,重力作用下液滴铺展现象明显且蒸发加快。基板的接触角越小,液滴铺展现象越明显,其接触直径越大,蒸发越快。当环境温度与基板温度相差较大时,液滴内部出现涡流,强化换热使蒸发过程加快。     

非牛顿流体液滴铺展过程的格子Boltzmann模拟北大核心CSCD

格子Boltzmann可有效模拟气液两相流动,成功应用于牛顿流体动态湿润领域和非牛顿流体单相流动领域的研究。将牛顿流体的两相流动格子Boltzmann模型与非牛顿流体单相流动格子Boltzmann模型相结合,建立了幂指流体两相流动格子Boltzmann模型,模拟了牛顿流体和幂指流体液滴铺展过程,并将液滴半径随时间的变化规律模拟结果与液滴铺展理论模型进行了比较,对两类流体的铺展指数与流变指数变化关系进行了分析讨论。     

高温高压下机油液滴蒸发特性的数值模拟

建立了高温高压环境下双组分单液滴的一维非稳态蒸发模型.该模型可描述气/液两相质量及能量平衡、液相传热传质和相变过程.使用所建立的液滴蒸发模型,以正二十四烷(C_(24)H_(50))和正三十烷(C_(30)H_(62))作为机油的表征组分,分析了双组分机油液滴蒸发过程中液滴温度和组分摩尔分数分布的变化趋势,并对比了相同环境条件下机油与异辛烷液滴的不同蒸发特性.在此基础上,研究了环境压力、环境温度和液滴初始半径对机油液滴寿命、液滴蒸发百分数、液滴温度和组分摩尔分数等的影响.结果表明:在高温高压环境下,机油液滴能够留存较长时间,形成缸内高温早燃源的可能性较大.     

快速降压环境下盐水液滴蒸发过程的传热传质研究

采用VOF（Volume of Fluid）自由表面捕捉方法对盐水液滴蒸发过程中气液界面进行追踪,建立了降压环境下单个盐水液滴的蒸发模型,并通过盐水液滴蒸发的实验数据验证了此模型。通过对盐水液滴在相变过程中的形态变化以及传热传质特性的分析,研究了液滴内部温度、速度、蒸汽分布以及液滴形态等随时间的变化情况,分析了影响盐水液滴降压蒸发过程的主要因素。结果表明：在降压蒸发过程中液滴形态变化和环境中蒸汽的分布会随速度场的变化而变化;蒸发过程中初始盐组分质量浓度越大的液滴蒸发速率越缓慢,最终能达到的液滴最低中心温度越高,且液滴中心温度回升速度越慢、回升时间也越晚;液滴初始温度对蒸发速率影响较大,初始温度越高,表面蒸发速率越快,液滴中心温度回升速度越快。     

超音速分离管中水蒸气凝结相变的数值研究

基于超音速分离管中混合气体流动属于伴随凝结相变的可压缩、跨音速的特点,建立了考虑传质效应与非平衡凝结过程的数学模型,并采用数值方法对伴随水蒸气凝结的超音速分离管中的流动进行分析研究。以空气、水蒸气及液态水为流动介质,采用两相流动中的VOF模型结合凝结相变模型以及组分传输模型,研究不同进出口参数及不同水蒸气含量对凝结流场的影响。研究结果表明,所建立的分离管内部非平衡凝结相变模型可以较好的再现超音速流中的凝结成核及液滴生长过程;数值计算结果表明,入口压力、温度及水蒸气含量对分离管内流动凝结过程有直接且重要的影响。因此在进行超音速分离管设计时,考虑温度压力参数的同时,考虑水蒸气含量对分离管性能的影响也是非常重要的。     

非能动安全壳竖直平板降膜流动特性数值模拟

核电安全日益受到关注,非能动系统作为第三代核电系统具有很高的安全性。采用FLUENT流体体积分数(volume of fraction,VOF)模型和k-ε湍流模型对非能动安全壳冷却系统(passive containment cooling system,PCCS)三维平板降膜流动进行数值模拟。结果表明:1)在降膜过程中有波动现象,最终波动趋于平缓;2)水与空气逆流流动过程中发生轻微的液滴夹带;3)降膜流动受重力、表面张力与壁面黏滞力共同作用,液膜厚度沿横向分布均匀,沿高度方向平均液膜厚度越来越小,并且受进口水流速度与入口宽度影响,水流量一定时增加进口水流速度与入口宽度,平均液膜厚度增大,空气入口流速对水膜厚度影响相对较小。     

基于混合热格子玻尔兹曼模型的液滴蒸发数值模拟

采用格子玻尔兹曼方法(LBM)的单组分伪势模型与有限差分耦合的混合热格子玻尔兹曼模型(TLBM)对液滴蒸发过程进行了研究。首先,通过对液滴在方腔内蒸发过程进行模拟,验证了所采用计算方法及程序的有效性。随后,模拟了液滴撞击高温壁面后的蒸发过程,研究了壁面温度、液滴邦德数和液滴雷诺数对蒸发过程的影响。结果表明,壁面温度、液滴邦德数和液滴雷诺数的增加均会造成液滴撞击高温壁面后蒸发速率的增大。     

多壁碳纳米管水基纳米流体重力热管传热特性数值模拟

为揭示多壁碳纳米管水基纳米流体应用在重力热管中的传热特性,基于多相流模型(VOF)建立其重力热管数值模型,并将数值结果与实验数据进行对比验证。以热阻作为性能评价指标,改变加热功率和充液率,讨论二者对热管换热性能的具体影响。通过添加传热传质源项来编写用户自定义函数(UDF)完成内部流体蒸发冷凝过程中的相变模拟。模拟结果表明:该数值模型能够较好模拟多壁碳纳米管水基纳米流体应用重力热管内部复杂的流动与传热过程;在选定的加热功率及充液率参数范围内,该重力热管的整体热阻随蒸发段加热功率的增大而减小,随充液率的增大而增大。     

超临界工况下二甲醚蒸发的数值模拟

基于瞬态气液相控制方程,建立了超临界工况下二甲醚(DME)蒸发的数学模型,并进行了数值求解。详细考虑了高压下的非理想气体效应、液滴表面热动力相平衡、环境气体的可溶性、流体热物性的高压修正以及Soret和Dufour耦合效应等。计算结果表明,高压下气体溶解显著增强,相变焓明显偏离蒸发潜热,d~2律已不能准确描述液滴蒸发现象,液滴表面温度一直处于上升状态,甚至可以达到临界点;而且,随着压力的上升,耦合效应对蒸发的影响逐渐增强。     

掺混燃油在离心喷嘴内的流动与雾化特性研究

燃料性质的改变会导致雾化特性的变化,针对乙醇掺混航空煤油在离心式压力雾化喷嘴内的流动与雾化特性开展了研究。通过耦合流体体积法(VOF)和离散相模型(DPM),研究了不同乙醇掺混体积分数下掺混燃油在离心式喷嘴中的内部流动和外部雾化过程。研究结果表明：在压差不变时,喷嘴内空气芯直径随着掺混燃油内乙醇体积分数的增加而增大;而液膜厚度则与空气芯直径成反比,随着乙醇体积分数的增加而减小。喷嘴出口的速度随着乙醇体积分数的增加而增大;在油膜表面的波动及气动力的共同作用下,油膜失稳形成液滴,获得了不同比例下掺混燃油在喷雾外流场内的喷雾粒径分布特征,随着掺混乙醇体积分数的增加,液滴的平均直径逐渐减小。     

A numerical investigation of the evaporation process of a liquid droplet impinging onto a hot substrate

A numerical investigation of the evaporation process of n-heptane and water liquid droplets impinging onto a hot substrate is presented. Three different temperatures are investigated, covering flow regimes below and above Leidenfrost temperature. The Navier–Stokes equations expressing the flow distribution of the liquid and gas phases, coupled with the Volume of Fluid Method (VOF) for tracking the liquid–gas interface, are solved numerically using the finite volume methodology. Both two-dimensional axisymmetric and fully three-dimensional domains are utilized. An evaporation model coupled with the VOF methodology predicts the vapor blanket height between the evaporating droplet and the substrate, for cases with substrate temperature above the Leidenfrost point, and the formation of vapor bubbles in the region of nucleate boiling regime. The results are compared with available experimental data indicating the outcome of the impingement and the droplet shape during the impingement process, while additional information for the droplet evaporation rate and the temperature and vapor concentration fields is provided by the computational model.     

Thermal and solutal effects on convection inside a polymer solution droplet on a substrate

The convective flow inside polymer solution droplets drying on a lyophobic substrate is numerically studied. The evaporating droplet is presumed as a hemisphere shrinking with time at the constant contact angle. The thermal and solutal effects are simultaneously considered in the computation. The thermal Marangoni convection is induced due to the quick thermal diffusion, and this convection transports the solute resulting in the solutal Marangoni flow. The solutal dependence corresponds to our previous experimental work, but the flow pattern does not. Consideration of the pseudo evaporation rate distribution depending on the contact angle yields to the flow pattern correspondence.     

Onset of Marangoni instability of a two-component evaporating droplet

The temperature and solute concentration reductions across a thin boundary layer near the free surface of an evaporating droplet may induce cellular flow motion in the droplet because of Marangoni instability. The present study is aimed at investigating theoretically the onset of Marangoni instability due to the evaporation of a two-component evaporating droplet.

With the quasi-steady approximation which means that the surrounding gas motion is asymptotically steady, the size change of the droplet is negligible, and the temperature and concentration distributions of the droplet are temporarily frozen at each specified instant of interest, the onset condition for Marangoni instability is obtained through the linear stability analysis.

By assuming the surface tension is a monotonically decreasing function of both temperature and concentration of the higher-volatility substance, the thermocapillary and diffuso-capillary effects augment each other. Therefore, the theoretical analysis predicts a linear relation, with a negative slope, between the onset thermal Marangoni number, Ma_T, and the onset solute Marangoni number, Ma_S. Moreover, when liquid Lewis number Le_l>1, the critical wave number, l_c, may possess different values depending on the variation of the thermocapillary effect and diffuso-capillary effect. In addition, Le_l has a stronger effect on the critical solute Marangoni number Ma_S,C, than on the critical thermal Marangoni number Ma_T,C. That is, as Le_l decreases, Ma_T,C decreases mildly while Ma_S,C increases drastically.     

Numerical investigation on the evaporation of droplets depositing on heated surfaces at low Weber numbers

The evaporation of water droplets, impinging with low Weber number and gently depositing on heated surfaces of stainless steel is studied numerically using a combination of fluid flow and heat transfer models. The coupled problem of heat transfer between the surrounding air, the droplet and the wall together with the liquid vaporisation from the droplet’s free surface is predicted using a modified VOF methodology accounting for phase-change and variable liquid properties. The surface cooling during droplet’s evaporation is predicted by solving simultaneously with the fluid flow and heat transfer equations, the heat conduction equation within the solid wall. The droplet’s evaporation rate is predicted using a model from the kinetic theory of gases coupled with the Spalding mass transfer model, for different initial contact angles and substrate’s temperatures, which have been varied between 20–90° and 60–100 °C, respectively. Additionally, results from a simplified and computationally less demanding simulation methodology, accounting only for the heat transfer and vaporisation processes using a time-dependent but pre-described droplet shape while neglecting fluid flow are compared with those from the full solution. The numerical results are compared against experiments for the droplet volume regression, life time and droplet shape change, showing a good agreement.     

Internal flow in evaporating droplet on heated solid surface

This work numerically studies the evaporation process of a liquid droplet on a heated solid surface using a comprehensive model. The internal flow within the evaporating liquid droplet is elucidated, while considering the effects of buoyancy force, thermocapillary force, and viscous resistance. The evaporation process is modeled by simultaneously solving the Navier–Stokes equations and energy equation for the liquid domain and the heat conduction equation for the solid domain, while assuming the liquid–vapor interface is a free surface. Three dimensionless parameters are utilized to describe the contribution of individual driving forces to internal flow. Evolutions of the thermal and internal flows during evaporation are discussed. The volume evolution and experimental data are in good agreement.     

Interaction of the transfer processes in semitransparent liquid droplets

The study presents the mathematical model of unsteady heat transfer in evaporating semitransparent droplets of non-isothermal initial state and the numerical research method, evaluating selective radiation absorption and its influence on the interaction of transfer processes. The relation of the transfer processes inside droplets and in their surroundings and the necessity of thorough research of these processes are substantiated. When modeling the combined energy transfer in water droplets, the evaluation of thermoconvective stability in evaporating semitransparent liquid droplets is presented; the influence of the droplet initial state on its heating and evaporation process is investigated. The influence of heat transfer peculiarities on the change of the evaporating droplet state is indicated. Main parameters, which decide the peculiarities of the interaction of unsteady transfer processes in droplets and their surroundings, are discussed. The results of the numerical research are compared to the known results of the experimental studies of water droplet temperature and evaporation rate.     

Heat convection within evaporating droplets in strong aerodynamic interactions

The temperature field within evaporating ethanol droplets is investigated, relying on the two-color laser induced fluorescence (LIF) measurement technique and on a Direct Numerical Simulation (DNS). The configuration studied corresponds to a monodisperse droplet stream in a diffusion flame sustained by the droplet vapor. An experimental probe volume, small compared to the droplet size, is used to characterize the temperature field within the droplets, whereas DNS takes into account key aspects of the droplet heating and evaporation such as the non-uniform and transient stress, and the mass and heat transfer coefficients at the droplet surface. These investigations reveal that the frictional stresses are strongly reduced due to the small spacing between the droplets. They also show that the Marangoni effect has a significant influence on the internal motion and hence on the internal temperature field.     

Free convection flow with thermal radiation and mass transfer past a moving vertical porous plate

The combined free convection boundary layer flow with thermal radiation and mass transfer past a permeable vertical plate is studied when the plate moves in its own plane. The plate is maintained at a uniform temperature with uniform species concentration and the fluid is considered to be gray, absorbing–emitting. The coupled unsteady non-linear momentum, energy and concentration equations governing the problem is obtained and made similar by introducing a time-dependent length scale. The similarity equations are solved numerically using superposition method. The resulting velocity, temperature and concentration distributions are shown graphically for different values of parameters entering into the problem. The numerical values of the local wall shear stress, local surface heat and mass flux are shown in tabular form.     

Motion of droplets induced by the Marangoni force on a wall with a temperature gradient

Motion of silicone oil and water droplets induced by the Marangoni force was numerically simulated by using two‐ and three‐dimensional second‐order finite difference methods with the CIP and the level set methods. The surface tension was introduced by the continuum surface force (CSF) method. The results clearly showed the flow induced by the Marangoni force and the dependence of droplet velocity on droplet size, contact angle, temperature gradient, and fluid properties. The Marangoni force balanced with the viscous force in the small contact angle case; on the other hand, in the large contact angle case, it balanced with the normal component of surface tension. As for the effect of fluid properties on droplet motion, the temperature coefficient of surface tension had a much larger effect than did viscosity, thermal diffusivity, or surface tension. © 2004 Wiley Periodicals, Inc. Heat Trans Asian Res, 33(2): 81–93, 2004; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.20004     

Numerical solution for film evaporation of a spherical liquid droplet on an isothermal and adiabatic surface

A numerical solution for the problem of film evaporation of a liquid droplet on a horizontal surface is presented. The droplets are small enough to be assumed spherical. Two principal cases are considered: (1) the horizontal surface is maintained at a constant temperature (case I), and (2) the surface is insulated while the ambience is hot (case II). The complete set of equations governing this problem were solved under the following assumptions : (1) evaporation is quasi-steady, (2) no internal liquid circulation, (3) constant properties, and (4) the droplet temperature is spatially uniform but temporally varying. The Lewis number is not assumed to be unity; gas phase viscous effects, Stefan type convection, and gas phase inertia are included in the analysis. The total droplet evaporation time was found to decrease with increasing plate (I) or ambient (II) temperature as expected, and the droplet progressively moves away from the plate as it evaporates. The numerical results agree with the analytical solution for film evaporation of a droplet above an adiabatic surface in a hot ambience in the limit of large effective Reynolds number (i.e. potential flow).