An Eulerian–Lagrangian hybrid model for the coarse-grid simulation of turbulent liquid jet breakup

In this paper we present a numerical model for the coarse-grid simulation of turbulent liquid jet breakup using an Eulerian–Lagrangian coupling. To picture the unresolved droplet formation near the liquid jet interface in the case of coarse grids we considered a theoretical model to describe the unresolved flow instabilities leading to turbulent breakup. These entrained droplets are then represented by an Eulerian–Lagrangian hybrid concept. On the one hand, we used a volume of fluid method (VOF) to characterize the global spreading and the initiation of droplet formation; one the other hand, Lagrangian droplets are released at the liquid–gas interface according to the theoretical model balancing consolidating and disruptive energies. Here, a numerical coupling was required between Eulerian liquid core and Lagrangian droplets using mass and momentum source terms. The presented methodology was tested for different liquid jets in Rayleigh, wind-induced and atomization regimes and validated against literature data. This comparison reveals fairly good qualitative agreement in the cases of jet spreading, jet instability and jet breakup as well as relatively accurate size distribution and Sauter mean diameter (SMD) of the droplets. Furthermore, the model was able to capture the regime transitions from Rayleigh instability to atomization appropriately. Finally, the presented sub-grid model predicts the effect of the gas-phase pressure on the droplet sizes very well.     

液体环轴对称抛洒首次破碎的理论分析

对于水和无水乙醇环形轴对称抛洒实验的首次破碎结果的分析表明，在给定实验条件下，液体密度大约是气体的1000倍，液体首次破碎的过程分2个阶段:首先，界面不稳定性非线性发展，液体尖钉在惯性力和空气阻力的作用下，拉伸变细变长；然后在射流不稳定性的作用下，断裂成液体珠串。这个理论分析给出的首次破碎的平均液体直径与实验结果基本一致。理论分析表明，无水乙醇破碎液滴的平均直径比水液滴小，是无水乙醇表面张力值仅为水的1/3的一个必然结果。     

A unified spray model for engine spray simulation using dynamic mesh refinement

This study is based on dynamic mesh refinement and uses spray breakup models to simulate engine spray dynamics. It is known that the Lagrangian discrete particle technique for spray modeling is sensitive to gird resolution. An adequate spatial resolution in the spray region is necessary to account for the momentum and energy coupling between the gas and liquid phases. This study uses a dynamic mesh refinement algorithm that is adaptive to spray particles to increase the accuracy of spray modeling. On the other hand, the accurate prediction of the spray structure and drop vaporization requires accurate physical models to simulate fuel injection and spray breakup. The present primary jet breakup model predicts the initial breakup of the liquid jet due to the surface instability to generate droplets. A secondary breakup model is then responsible for further breakup of these droplets. The secondary breakup model considers the growth of the unstable waves that are formed on the droplet surface due to the aerodynamic force. The simulation results are compared with experimental data in gasoline spray structure and liquid penetration length. Validations are also performed by comparing the liquid length of a vaporizing diesel spray and its variations with different parameters including the orifice diameter, injection pressure, and ambient gas temperature and density. The model is also applied to simulate a direct-injection gasoline engine with a realistic geometry. The present spray model with dynamic mesh refinement algorithm is shown to predict the spray structure and liquid penetration accurately with reasonable computational cost.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.     

圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布相关性研究

压力旋流喷嘴被广泛应用于航空发动机、船用发动机、车用汽油缸内直喷发动机、燃气轮机等动力机械的燃油喷射系统中.以压力旋流喷嘴射流为研究对象,开展了圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布相关性问题研究.对于液体射流,以往的研究往往对破碎液滴粒径数量密度分布或速度数量密度分布进行单独研究,对于这两种数量密度分布之间关系的研究较少;从相关性的角度对圆环旋转黏性液体射流破碎液滴粒径与速度数量密度分布之间的关系进行研究.采用最大熵原理方法建立了圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数.对圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数进行了讨论,对圆环旋转黏性液体射流破碎液滴粒径数量密度分布与速度数量密度分布的相关性问题进行了研究.研究结果表明,为了给出正确的圆环旋转黏性液体射流破碎液滴粒径与速度联合概率密度函数,射流守恒约束条件中必须同时包括质量守恒定律、动量守恒定律以及能量守恒定律;破碎液滴粒径的数量密度分布与速度数量密度分布密切相关;射流旋转强度对破碎液滴粒径数量密度与速度数量密度分布结构影响不大,对破碎液滴粒径数量密度和速度数量密度的分布区域影响较大.      

爆炸驱动液体分层现象和液滴尺寸模拟

爆炸驱动液体介质外界面的分散和破碎是气溶胶云团形成的重要过程。采用基于维数分裂的欧拉程序和Youngs混合界面处理方法,对中心药爆炸驱动甘油和水介质流场的液体分层现象进行了数值模拟。结合试验结果推断提出了液滴形成过程的三种并存机制:外层射流破碎、内层R-T失稳和中间液层"空化"破碎,分别建立了不同液层破碎液滴的尺寸模拟方法。对比给出抛撒甘油和水装置初级液滴的尺寸分布及最外层理论射流量。     

Bag breakup of nonturbulent liquid jets in crossflow

An experimental investigation of the bag breakup of round nonturbulent liquid jets in gaseous crossflow at room temperature and pressure is described. Pulsed photography, pulsed shadowgraphy, and high-speed imaging were used to observe the column and surface waves along the liquid jet and the formation and breakup of bags. Measurements included: wavelengths of column and surface waves, jet velocities, the number of bags along the liquid jet, the number of nodes per bag, droplets sizes and velocities, and trajectories of droplets. Present results show that the column waves of a nonturbulent liquid jet in crossflow within bag breakup regime can be explained based on Rayleigh–Taylor instability. The number of nodes per bag affected the breakup mechanism of the bags. Three distinctive sizes of droplets were produced due the breakup of the bag membrane, the ring strings and the ring nodes. The size of the droplets resulting from the breakup of the bag membrane was constant independent of the crossflow Weber number. Finally different trajectories were observed for the three groups of droplets.     

液滴在激波冲击下的破裂过程

对液滴在入射激波作用下的变形破碎过程进行了实验研究和数值模拟,得知数值模拟结果与实验结果基本吻合,以及在什么情况下两者出现分歧。结果显示,液滴在激波的作用下要经历从压缩变形、RM不稳定性变形、细小液雾剥离到全部雾化破碎等过程。结果还表明,不同液滴直径、入射激波马赫数和液滴介质等参数下的液滴变形破碎的发展趋势是一致的,而其发展速度明显则不同。其中Weber数的增加加速了液滴的破碎,而Ohnesorge数和黏性的增加则抑制了液滴的破碎。     

The effect of surfactants on the instability of a rotating liquid jet

It has long been known that the presence of surfactants on the free surface of a liquid jet can create surface tension gradients along the interface. The resulting formation of tangential stresses along the surface lead to Marangoni type flows and greatly affect the resulting dynamics of rupture. In this way surfactants can be used to manipulate the breakup of a liquid jet and control the size of droplets produced. In this paper we investigate the effects of insoluble surfactants on the breakup of rotating liquid jets with applications to industrial prilling. Using a long wavelength approximation we reduce the governing equations into a set of one-dimensional equations. We use an asymptotic theory to find steady solutions and then carry out a linear instability analysis on these solutions. We show that steady state centreline solutions are independent of viscosity to leading order and that the most unstable wavenumber and growth rate of disturbances decrease as the effectiveness of surfactants is increased. We also numerically solve these equations using a finite difference scheme to investigate the effects of changing the initial surfactant concentration and other fluid parameters. Our results show that differences in breakup lengths between rotating surfactant-laden jets and surfactant-free jets increase with the rate of rotation. Moreover, we find that satellite droplet sizes decrease as the rate of rotation is decreased with the effect of surfactants amplifying the reduction in sizes. Furthermore, the presence of surfactants at fixed rotation rates is shown to produce larger main droplets at low disturbance wavenumbers whilst satellite droplets are smaller for moderate disturbance wavenumbers κ≈0.7.     

可压缩旋转气体中超空化射流的热稳定性

液体射流热稳定性研究是对射流稳定性问题的更深层次的探讨,可以进一步加深对液体射流分裂与雾化机理的认识,具有重要的学术意义和工程应用价值. 基于射流稳定性理论,在同时考虑射流周围气体旋转、射流和周围气体可压缩性以及射流液体中含空化气泡的条件下,建立了描述可压缩旋转气体中超空化射流热稳定性的数学模型,并对数学模型及其求解方法进行了验证分析;在此基础上,分析了液体射流表面与周围气体间温差及射流内部温度梯度同时作用下对射流稳定性的影响;并进一步探讨了超空化射流的热稳定性. 结果表明,射流表面扰动波的最大扰动增长率、最不稳定频率以及最大扰动波数皆随气液温差的增大呈近似线性增大趋势;射流内部温度梯度的存在使得气液温差对射流的失稳作用更加显著;射流内部温度梯度会抑制超空化对射流稳定性的影响,但气液温差会在一定程度上促进超空化对射流的失稳作用.      

Flow visualization of shock/water column interactions

The breakup of a liquid droplet induced by a high speed gas stream is a typical multiphase flow problem. The shock/droplet interaction is the beginning stage of the droplet breakup. Therefore, investigation of the shock/droplet interactions would be a milestone for interpreting the mechanism of the droplet breakup. In this study, a compressible multiphase solver with a five-equation model is successfully developed to study shock/water column interactions. For code validation, interface-only, gas–gas shock tube, and gas–liquid shock tube problems are first computed. Subsequently, a planar shock wave interacting with a water column is simulated. The transmitted wave and the alternative appearances of local high- and low-pressure regions inside the water column are observed clearly. Finally, a planar shock wave interacting with two water columns is investigated. In this work, both horizontal and vertical arrangements of two water columns are studied. It is found that different arrangements can result in the diversity of the interacting process. The complex flow structures generated by shock/water column interactions are presented by flow-visualization techniques.      

竖直振动无黏液滴的法拉第不稳定性分析

由于外部周期性的振动而在液滴表面产生的Faraday不稳定效应广泛存在于超声雾化、喷涂加工等应用中,对Faraday不稳定性进行分析对研究振动液滴的表面动力学有着重要意义.本文将Faraday不稳定性问题从径向振动拓展到竖直振动,研究了竖直振动无黏液滴表面波的不稳定性.竖直方向的振动使得液滴动量方程为含有空间相关项和时...     

Droplet impact dynamics for two liquids impinging on anisotropic superhydrophobic surfaces

Droplet impingement experiments were performed on grooved hydrophobic surfaces with cavity fractions of 0, 80, and 93?% using droplets of water and a 50?%/50?% water/glycerol mixture. The influence of liquid viscosity, cavity fraction, and spreading direction, relative to the surface grooves, is explored qualitatively and quantitatively. The maximum droplet spread diameter, velocity of the rebounding jet, and the time delay between droplet impact and jet emission were characterized for Weber numbers, We, based on droplet impact speed and diameter, up to 500. The unequal shear stresses and contact angles influence the maximum spread diameters in the two primary spread directions. At We?>?100, the ratio of the spread diameter along the direction of the grooves to the spread diameter perpendicular to the grooves increases above unity with increasing We. The maximum droplet spread diameter is compared to recent predictive models, and the data reveal differing behavior for the two fluids considered. The results also reveal the existence of very high relative jet velocities in the range 5????We????15 for water droplets, while such jets were not observed for the more viscous mixture. Further, in the range 115????We????265, the water/glycerol jet formation dynamics are radically different from the water behavior. Most evident is the existence of two-pronged jets, which arise from the anisotropy of the surface and the unequal shear stresses and contact angles that prevail on the surfaces. It is these influences that give rise to differences in the maximum spread diameters in the two primary spread directions. Similar two-pronged jet emission was observed for water over the very narrow range of We from 91 to 96. The issuing jet velocities were also observed to increase with increasing cavity fraction for both fluids and over the entire range of We explored. Lastly, the elapsed time between droplet impact and jet emission decreased with increasing cavity fraction.     

Shock–Free Breakup of Droplets. Temporal Characteristics

In the present paper, we consider the shock–free breakup of droplets in their encounter with a layer (sheet) of a moving gas in the absence of pressure perturbations when the droplets are affected by a short U–shaped pulse of aerodynamic forces. Under a high pressure of the ambient gas medium p₀ = 20—80 bar, the droplets (ethanol or liquid oxygen) have a chance to break up after stay in a thing (2—5 mm thick) gas layer (jet) moving with a velocity of 1—10 m/sec. A distinctive feature of the process is that the characteristic time of droplet deformation and the period of natural oscillations coincide with the residence time for the droplets in the region of their interaction with the gas stream. Empirical formulas are proposed for determination of the total breakup time and the duration of the droplet disintegration stage in shock–free breakup.     

受径向振荡激励的黏弹性液滴稳定性分析

当液滴受到外部周期性的径向激励时, 在其表面会形成驻波模式的不稳定, 这就是在球面上的Faraday不稳定问题. 不稳定的表面波的振荡频率根据流体物性参数和所施加激励条件的不同呈现为谐波或是亚谐波模式的振荡. 本文基于线性小扰动理论, 研究了受径向振荡体积力的黏弹性液滴表面波的不稳定性. 振荡的体积力导致动量方程为含有时间周期系数的Mathieu方程, 系统因此变成参数不稳定问题, 采用Floquet理论进行求解. 本模型中将黏弹性的特征处理为与流变模型参数相关的等效黏度, 从而简化了问题的求解. 基于对中性稳定曲线及增长率的分析, 研究了黏弹性参数对液滴稳定性的影响. 结果表明零剪切黏度和应变驰豫时间的增加具有抑制液滴表面波增长的作用, 提高了使液滴表面发生谐波不稳定的激励幅值. 随着振荡幅值的增加, 增长率不稳定的区域减少, 且随着振荡频率的增加, 液滴表面波增长率减小. 通过对增长率的分析可以得出, 应力松弛时间的增加使得增长率增加, 从而促进了液滴表面波的增长.      

Breakup of Newtonian and non-Newtonian fluids in air jets

The breakup of droplets of non-Newtonian fluids has been investigated by high speed photography and impaction following preliminary results of Newtonian fluids, which confirmed the suitability of the measurement techniques. Single droplets with diameters from 2.4 to 3.3 mm, were arranged to fall under gravity into a jet of air with velocities up to 36o m/s. The droplets of Newtonian fluids, water and Diesel oil, were atomised in the expected manner within three main regimes characterised by the Weber number of the droplet and air jet conditions, while similar droplets of non-Newtonian fluids were found not to atomise but to develop under shear and stretching into ligaments of fluid separated from a local region of their surface; these ligaments were elongated until breakup occurred, though not into small droplets as with the Newtonian fluids. Some of the non-Newtonian fluids (TEP with 7.5% and lo% K125, with and without water) were found not to break up at the maximum speed of the tests and they will be re-examined at higher jet velocities. Increase in the concentration of K125 in TEP resulted in higher critical speed for a given droplet diameter.The authors would like to thank Prof. J. H. Whitelaw of Imperial College for many discussions and useful suggestions during the course of this work, and Dr. G. Cambray of CBDE for his valuable administrative support     

Large eddy simulation and experiment of shear breakup in liquid-liquid jet: Formation of ligaments and droplets

Understanding the shear breakup in jet flows and the formation of droplets from ligaments is important to determine the final droplet size distribution (DSD). The initial droplet size, which affects the final DSD, is considered to be generated by the shear breakup. Large eddy simulation (LES) was performed to investigate the shear breakup in liquid-liquid jet flows. The explicit Volume of Fluid (VOF) model with the geometric reconstruction scheme was used to capture the oil-water interface. The estimated oil distribution including wave peaks, ligaments, droplets and water streaks were compared to the experiments with a good agreement. The estimated DSD matched with the measurements favorably well. In the simulation, the formation of droplets with a smooth and curved surface from ligaments or sheet-like structures was obtained. Different mechanisms were observed along with the shear layer including the formation of droplets from ligament through the capillary forces, breakage of a droplet into smaller ones and attachment of a droplet to a ligament. The destructive shear forces and resisting surface tension forces were quantified on stretching and retracting ligaments. The influence of internal viscous force was found to be negligible due to low oil viscosity. The critical capillary number was found to be larger than 5.0 for ligaments breaking with the shear breakup. The capillary number was below unity for retracting ligaments. The coalescence of two equal-sized droplets was obtained in the shear breakup region. The shear stress magnitude at the contact region increased more than two folds. The total surface area decreased nearly 20% after the coalescence.     

Modeling of droplet collision-induced breakup process

The present article proposes a new droplet collision model considering droplet collision-induced breakup process with the formation of satellite droplets. The new model consists of several equations to investigate the post-collision characteristics of colliding droplets and satellite droplets. These equations are derived from the conservations of droplet mass, momentum, and energy between before and after collision, and make it possible to predict the number of satellite droplets, and the droplet size and velocity in the analytical way. To validate the new collision model, numerical calculations are performed and their results are compared with experimental data published earlier for binary collision of water droplets. It is found from the results that the new model shows good agreement with experimental data for the number of satellite droplets. It can be also shown that the predicted mean diameter by the new model decrease with increasing the Weber number because of the collision-induced breakup, whereas the O’Rourke model fails to predict the size reduction via the binary droplet collision.     

Study of the shock-induced acceleration of hexane droplets

An experimental study of the interaction of a shock wave with a hexane droplet is presented. The main goal of the experiments was to record images of the process and measure basic parameters describing movement, dispersion and evaporation of the droplets engulfed by a shock wave propagating in air. A shock tube with a visualization section was used for this research. Photography of the process allowed one to measure the positions, velocities and sizes of mist clouds created by the interaction processes. Analysis of the pictures shows that there is no qualitative difference between cases for different size droplets, but shock Mach number had a significant effect on the process. Quantitative analysis shows that under certain conditions, a catastrophic breakup mechanism of dispersion occurred. The droplets are shattered into a mist cloud before they achieve mechanical equilibrium with the surrounding gas. The approximate time for the complete dispersion and acceleration of the fuel droplet varies from 300 to 500 μs, and depends both on the droplet diameter and shock velocity. The dispersion time is controlled principally by the droplet diameter, and to a lesser extent, the shock Mach number. This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems, Montreal, Canada, July 31–August 5, 2005.     

An experiment on the breakup of impinging droplets on a hot surface

Breakup characteristics of liquid droplets impinging on a hot surface are investigated experimentally with the wall temperatures in the Leidenfrost temperature range of 220–330°C for n-decane fuel. Factors influencing droplet breakup are wall temperature, impinging velocity, droplet diameter and impinging angle. The 50% breakup probability shows that the impinging velocity decreases linearly with the droplet diameter increase and there exists an optimum impinging angle near 80° having the minimum value in the impinging velocity for given wall temperature and droplet size. Near the wall temperature of 250°C corresponding to the Leidenfrost temperature, a peculiar nonlinear behavior in the breakup probability is observed.This work was supported by the Turbo and Power Machinery Research Center, Seoul National University.