高Bond数下黏性液滴的Rayleigh-Taylor不稳定性

给出了高Bond数下黏性液滴表面Rayleigh-Taylor线性不稳定性的分析解，这种不稳定性对于超音速气流作用下液滴破碎的早期阶段起着至关重要的作用．基于稳定性分析的结果，导出了用于估算稳定液滴的最大直径及液滴无量纲初始破碎时间的计算式，这些计算式与相关文献给出的实验和分析结果比较显示了良好的一致．     

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

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

轴对称液体环的形成、变形和破碎的研究

为了实现液体环的轴对称径向扩展运动,提出了一种新的实验方法及设备。这种新的实验设备是由垂直的无膜激波管和液体环发生器所组成,可以用来观察轴对称液体环的形成、变形和破碎过程。轴对称液体环抛撒过程的实验研究已在这种实验设备上完成,并得到了相应的一系列流场照片。对实验结果的分析表明,液体表面的不稳定性是造成液体环破碎的主要原因。     

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

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

在轴向气流作用下液体轴对称抛撒的近场研究

在轴向气流作用下液体轴对称抛撒的研究,是以飞行物体在运动状态下向大气抛撒液体燃料所导致的破碎和雾化为背景的。为了研究在轴向气流作用下液体轴对称抛撒所产生的雾场特性,本文提出了一种新的组合型实验设备。该设备由两台激波管、一套电子同步控制系统组成,可以观察在轴向气流作用下液体轴对称抛撒、破碎和雾化的过程。通过在此设备上的一系列实验,获得了在不同压力和不同气流速度下液体轴对称抛撒的近场纹影照片。通过对照片的研究发现,液体轴对称抛撒具有两个明显的阶段:液核生长阶段和液核稳定阶段。另外近场云雾区轮廓有明显的转折点,此转折点即为液核发生首次破碎的位置。进一步分析表明,轴向气流能促使液体轴对称抛撒首次破碎发生的时间缩短。     

液体环二次破碎所形成云雾颗粒尺寸测量和测试系统的标定

介绍了激光散射法测量颗粒尺寸系统的工作原理和标定结果,并对液体环轴对称抛撒进行了光学测量。实验结果表明,液体环二次破碎产生云雾区的液滴Ｓａｕｔｅｒ平均直径在固定点随时间的增加呈减小的趋势,而云雾区的宽度和云雾区前缘的液滴颗粒的Ｓａｕｔｅｒ平均直径则随测量的距离增加均有所增加。     

液体环轴对称抛撒首次破碎的实验研究和稳定性分析

首次破碎是液体抛撒、破碎过程中一个非常重要的阶段．本文提出了一种新的实验设备,并通过在这套实验设备上的一系列实验,得到了在不同激波马赫数、不同的液体种类和不同的抛撒容量下液体首次破碎过程的照片．实验结果表明,在确定的抛撒条件下,液体环的失稳,即不稳定的发生、发展并最终导致液体环破碎的过程,是由于液体环运动的加速度改变方向而引起的     

激波作用下液体喷射雾化的实验研究

液体在气流作用下的喷射抛撒过程是流体力学研究中令人感兴趣的重要领域,笔者采用改进后的激波管实现了激波作用下液体的喷射抛撒,并通过阴影照相和激光散射法分别对液体的抛撒状态和抛撒液滴的直径进行了测量,研究表明,在喷射过程中,流场中固定位置所测得的液滴Ｓａｒｔｅｒ平均直径随时间的发展而逐渐减小,在开始时时刻液滴直三小较快而最终渐趋平缓;在对不同抛撒距离雾化场的测量中发现,没位置测得的颗粒最大直 隧测得位     

运动状态下液体轴对称抛撒首次破碎的实验研究

在实验室条件下利用组合式激波管设备,对运动状态下液体轴对称抛撒进行了实验研究。通过纹影装置获得其所形成雾化场的外形轮廓照片,测量获得了液核发生首次破碎的位置与对称轴之间的距离。通过对抛撒过程中R-T不稳定性与K-H不稳定性的分析认为,轴向气流作用下液体轴对称抛撒的首次破碎点与对称轴的距离主要由轴向气流的速度、轴向气流的密度、液体轴对称抛撒的出口速度、抛撒液体的表面张力系数、环形喷口的宽度等参数所决定。在此基础上,利用相似性理论和无量纲分析,获得了运动状态下液体轴对称抛撒首次破碎点与对称轴之间的距离与相关参数的无量纲关系式。     

激波诱导气流与液幕、液柱相互作用的实验研究

利用激波管对激波诱导气流与液幕、液柱的相互作用进行了实验研究。通过比较发现,这种相互作用下的液体块变形破碎过程与以往对于液滴进行的研究结果很不相同。当激波与液幕相互作用时,阴影照片和直接照相都表明,液幕的变形破碎行为有很强的三维性,较之液滴的变形破坏机理更为复杂,并且在局部区域,初始时刻液幕破碎抛撒的速度相较激波诱导气流速度为快,本文应用一维变截面激波管理论对这一现象进行了理论分析。     

Application of maximum entropy method for droplet size distribution prediction using instability analysis of liquid sheet

This paper describes the implementation of the instability analysis of wave growth on liquid jet surface, and maximum entropy principle (MEP) for prediction of droplet diameter distribution in primary breakup region. The early stage of the primary breakup, which contains the growth of wave on liquid–gas interface, is deterministic; whereas the droplet formation stage at the end of primary breakup is random and stochastic. The stage of droplet formation after the liquid bulk breakup can be modeled by statistical means based on the maximum entropy principle. The MEP provides a formulation that predicts the atomization process while satisfying constraint equations based on conservations of mass, momentum and energy. The deterministic aspect considers the instability of wave motion on jet surface before the liquid bulk breakup using the linear instability analysis, which provides information of the maximum growth rate and corresponding wavelength of instabilities in breakup zone. The two sub-models are coupled together using momentum source term and mean diameter of droplets. This model is also capable of considering drag force on droplets through gas–liquid interaction. The predicted results compared favorably with the experimentally measured droplet size distributions for hollow-cone sprays.     

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.     

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.     

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.     

基于移动粒子半隐式法的表面张力模拟

采用移动粒子半隐式法(MPS)模拟了受表面张力影响的自由面流动。表面张力的计算采取了一种较适合于MPS方法的表面自由能模型。方形液滴振荡和射流断裂的模拟结果分别与理论分析和试验结果一致,同时进行了三维射流注水模拟,从而验证了MPS方法结合该表面张力模型可以有效、方便地进行自由面流动中表面张力现象的模拟。     

An experimental study of a water droplet impinging on a liquid surface

An experimental study is presented for water droplet impingement on a liquid surface. The impaction process was recorded using a high-speed digital camera at 1,000 frames/s. The initial droplet diameter was fixed at 3.1 mm ± 0.1 mm, and all experiments were performed in atmospheric air. The impact velocity was varied from 0.36 m/s to 2.2 m/s thus varying the impact Weber number from 5.5 to 206. The impacted liquid surface consisted of two fluids, namely water and methoxy-nonafluorobutane, C₄F₉OCH₃ (HFE7100). The depth of the water and HFE7100 pool was varied from 2 mm to 25 mm. The collision dynamics of water in the HFE7100 pool was observed to be drastically different from that observed for the water droplet impingement on a water pool. The critical impact Weber number for jet breakup was found to be independent of liquid depth. Water–HFE7100 impact resulted in no jet breakup over the range of velocities studied. Therefore, no critical impact Weber number can be defined for water–HFE7100 impact. Received: 27 June 2001/Accepted: 29 November 2001     

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.      

Surface breakup of a non-turbulent liquid jet injected into a high pressure gaseous crossflow

We employ detailed numerical simulations to understand the physical mechanism underlying the surface breakup of a non-turbulent liquid jet injected transversely into a high pressure gaseous crossflow under isothermal conditions. The numerical observations reveal the existence of shear instability on the jet periphery as the primary destabilization mechanism. The temporal growth of such azimuthal instabilities leads to the formation of interface corrugations, which are eventually sheared off of the jet surface as sheet-like structures. The sheets next undergo disintegration into ligaments and drops during the surface breakup process. The proposed instability mechanism is inherently an inviscid mechanism, contrary to the previously suggested mechanism of surface breakup (known as “boundary layer stripping”), which is relied on a viscous interpretation. The numerically obtained length and time scales of the shear instabilities on the jet laterals are compared with the results of Behzad et al. (2015) on temporal linear stability analyses of a jet in crossflow at near the nozzle. The stability characteristics of the most amplified modes (i.e., the wavenumber and the corresponding growth rate) obtained from the numerical simulations and the stability analyses are in good agreement.     

Nonlinear instability of an annular liquid sheet exposed to gas flow

Nonlinear instability and breakup of an annular liquid sheet has been modeled in this paper. The liquid sheet is considered to move axially and is exposed to co-flowing inner and outer gas streams. Also, the effect of outer gas swirl on sheet breakup has been studied. In the developed model a perturbation expansion method has been used with the initial magnitude of the disturbance as the perturbation parameter. This is a comprehensive model in that other geometries of planar sheet and a coaxial jet can be obtained as limiting cases of very large inner radius and inner radius equal to zero, respectively. In this temporal analysis, the effect of liquid Weber number, initial disturbance amplitude, inner gas-to-liquid velocity ratio, outer gas-to-liquid velocity ratio and outer gas swirl strength on the breakup time is investigated. The model is validated by comparison with earlier analytical studies for the limiting case of a planar sheet as well as with experimental data of sheet breakup length available in literature. It is shown that the linear theory cannot predict breakup of an annular sheet and the developed nonlinear model is necessary to accurately determine the breakup length. In the limiting case of a coaxial jet, results show that gas swirl destabilizes the jet, makes helical modes dominant compared to the axisymmetric mode and decreases jet breakup length. These results contradict earlier linear analyses and agree with experimental observations. For an annular sheet, it is found that gas flow hastens the sheet breakup process and shorter breakup lengths are obtained by increasing the inner and the outer gas velocity. Axially moving inner gas stream is more effective in disintegrating the annular sheet compared to axially moving outer gas stream. When both gas streams are moving axially, the liquid sheet breakup is quicker compared to that with any one gas stream. In the absence of outer gas swirl, the axisymmetric mode is the dominant instability mode. However, when outer gas flow has a swirl component higher helical modes become dominant. With increasing outer gas swirl strength, the maximum disturbance growth rate increases and the most unstable circumferential wave number increases resulting in a highly asymmetric sheet breakup with shorter breakup lengths and thinner ligaments.