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.     

Nonlinear analysis of the deformation and breakup of viscous microjets using the method of lines

We present a numerical model for predicting the instability and breakup of viscous microjets of Newtonian fluid. We adopt a one‐dimensional slender‐jet approximation and obtain the equations of motion in the form of a pair of coupled nonlinear partial differential equations (PDEs). We solve these equations using the method of lines, wherein the PDEs are transformed to a system of ordinary differential equations for the nodal values of the jet variables on a uniform staggered grid. We use the model to predict the instability and satellite formation in infinite microthreads of fluid and continuous microjets that emanate from an orifice. For the microthread analysis, we take into account arbitrary initial perturbations of the free‐surface and jet velocity, as well as Marangoni instability that is due to an arbitrary variation in the surface tension. For the continuous nozzle‐driven jet analysis, we take into account arbitrary time‐dependent perturbations of the free‐surface, velocity and/or surface tension as boundary conditions at the nozzle orifice. We validate the model using established computational data, as well as axisymmetric, volume of fluid (VOF) computational fluid dynamic (CFD) simulations. The key advantages of the model are its ease of implementation and speed of computation, which is several orders of magnitude faster than the VOF CFD simulations. The model enables rapid parametric analysis of jet breakup and satellite formation as a function of jet dimensions, modulation parameters, and fluid rheology. Copyright © 2010 John Wiley & Sons, Ltd.     

Temporal analysis of breakup for a power law liquid jet in a swirling gas

The breakup mechanism and instability of a power law liquid jet are investigated in this study. The power law model is used to account for the non-Newtonian behavior of the liquid fluid. A new theoretical model is established to explain the breakup of a power law liquid jet with axisymmetric and asymmetric disturbances, which moves in a swirling gas. The corresponding dispersion relation is derived by a linear approximation, and it is applicable for both shear-thinning and shear-thickening liquid jets. Analysis results are calculated based on the temporal mode. The analysis includes the effects of the generalized Reynolds number, the Weber number, the power law exponent, and the air swirl strength on the breakup of the jet. Results show that the shear-thickening liquid jet is more unstable than its Newtonian and shear-thinning counterparts when the effect of the air swirl is taken into account. The axisymmetric mode can be the dominant mode on the power law jet breakup when the air swirl strength is strong enough, while the non-axisymmetric mode is the domination on the instability of the power liquid jet with a high We and a low Re _n . It is also found that the air swirl is a stabilizing factor on the breakup of the power law liquid jet. Furthermore, the instability characteristics are different for different power law exponents. The amplitude of the power law liquid jet surface on the temporal mode is also discussed under different air swirl strengths.     

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.     

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.     

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.     

Instability of a compound liquid jet

In the linear Rayleigh theory [1] the degree of stability of a jet is determined by the viscosity and inertia characteristics of the fluids and the interphase surface tension. The stability of a jet in an infinite medium increases with increase in the viscosity of both the jet and the medium [2, 3]. The presence of two interfaces is responsible for various features of the development of instability in a liquid layer on the surface of a cylinder, and in particular a layer on the inner surface of a cylinder is more unstable than one on the outer surface [4]. In [5, 6] the breakup of a hollow jet in an external medium was investigated. In this paper we examine, in the linear approximation, the stability of a compound jet of nonmiscible liquids with respect to small axisynmetric perturbations of the interfaces. The instability characteristics are given for jets with inviscid and very viscous outer shells. The conditions governing the suppression of rapidly growing instabilities of the inner part (core) of the jet by a viscous shell are determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–8, July–August, 1985.     

The breakup and atomization of a viscous liquid jet

Based on the linear analysis of stability, a dispersion equation is deduced which delineates the evolution of a general 3-dimensional disturbance on the free surface of an incompressible viscous liquid jet. With respect to the spatial growing disturbance mode, the numerical results obtained from the solution of the dispersion equation reveal that a dimensionless parameterJ _e exists. AsJ _e>1, the axisymmetric disturbance mode is most unstable; and whenJ _e<1, the asymmetric disturbances come into being, their growth rate increases with the decrease, ofJ _e, till one of them becomes the most unstable disturbance. The breakup of a low-speed liquid jet results from the developing of axisymmetric disturbances, whose instability is produced by the surface tension; while the atomization of a high-speed liquid jet is brought about by the evolution of nonaxisymmetric disturbance, whose instability is caused by the aerodynamic force on the interface between the jet and the ambient gas. The project supported by the National Natural Science Foundation of China     

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.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Displacement of yield-stress fluids in a fracture

We consider a displacement of several yield-stress fluids in a Hele-Shaw cell. The topic is relevant to the development of a model for the flow of multiple phases inside a narrow fracture with application to hydraulically fracturing a hydrocarbon-bearing underground formation. Existing models for fracturing flows include only pure power-law models without yield stress, and the present work is aimed at filling this gap. The fluids are assumed to be immiscible and incompressible. We consider fluid advection in a plane channel in the presence of density gradients. Gravity is taken into account, so that there can be slumping and gravitational convection. We use the lubrication approximation so that governing equations are reduced to a 2D width-averaged system formed by the quasi-linear elliptic equation for pressure and transport equations for volume concentrations of fluids. The numerical solution is obtained using a finite-difference method. The pressure equation is solved using an iterative algorithm and the Multigrid method, while the transport equations are solved using a second-order TVD flux-limiting scheme with the superbee limiter. This numerical model is validated against three different sets of experiments: (i) gravitational slumping of fluids in a closed Hele-Shaw cell, (ii) viscous fingering of fluids with a high viscosity contrast due to the Saffman–Taylor (S–T) instability in a Hele-Shaw cell at microgravity conditions, (iii) displacement of Bingham fluids in a Hele-Shaw cell with the development of fingers due to the S–T instability. Good agreement is observed between simulations and laboratory data. The model is then used to investigate the joint effect of fingering and slumping. Numerical simulations show that the slumping rate of yield-stress fluid is significantly less pronounced than that of a Newtonian fluid with the same density and viscosity. If a low-viscosity Newtonian fluid is injected after a yield-stress one, the S–T instability at the interface leads to the development of fingers. As a result, fingers penetrating into a fluid with a finite yield stress locally decrease the pressure gradient and unyielded zones develop as a consequence.     

2D thermal/isothermal incompressible viscous flows

2D thermal and isothermal time‐dependent incompressible viscous flows are presented in rectangular domains governed by the Boussinesq approximation and Navier–Stokes equations in the stream function–vorticity formulation. The results are obtained with a simple numerical scheme based on a fixed point iterative process applied to the non‐linear elliptic systems that result after a second‐order time discretization. The iterative process leads to the solution of uncoupled, well‐conditioned, symmetric linear elliptic problems. Thermal and isothermal examples are associated with the unregularized, driven cavity problem and correspond to several aspect ratios of the cavity. Some results are presented as validation examples and others, to the best of our knowledge, are reported for the first time. The parameters involved in the numerical experiments are the Reynolds number Re, the Grashof number Gr and the aspect ratio. All the results shown correspond to steady state flows obtained from the unsteady problem. Copyright © 2005 John Wiley & Sons, Ltd.     

金属射流失稳断裂的理论分析

基于Hamilton原理,提出一个包含射流强度、剪切、应变率效应、动力学黏性、表面张力和速度梯度等多因素耦合的金属射流拉伸运动方程,具体分析了各种失稳因素,并由数值解定量给出其影响大小,以及最不稳定波长与初始应变率乘积\lambda_{m}\dot{\varepsilon}_{0} 值的变化范围;给出了射流断裂的时间判据和近似理论公式,计算得到的 t_{b}-\dot{\varepsilon}_{0}曲线与射流实验点、Chou\&Carleone拟合公式三者符合较好.      

Natural Convection Induced by a Centrifugal Force Field in a Horizontal Annular Porous Layer Saturated with a Binary Fluid

The Darcy Model with the Boussinesq approximation is used to study natural convection in a horizontal annular porous layer filled with a binary fluid, under the influence of a centrifugal force field. Neumann boundary conditions for temperature and concentration are applied on the inner and outer boundary of the enclosure. The governing parameters for the problem are the Rayleigh number, Ra, the Lewis number, Le, the buoyancy ratio, j{\varphi } , the radius ratio of the cavity, R, the normalized porosity, e{\varepsilon } , and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in a thin annular layer (R → 1), analytical solutions for the stream function, temperature and concentration fields are obtained using a concentric flow approximation and an integral form of the energy equation. The critical Rayleigh number for the onset of supercritical convection is predicted explicitly by the present model. Also, results are obtained from the analytical model for finite amplitude convection for which the flow and heat and mass transfer are presented in terms of the governing parameters of the problem. Numerical solutions of the full governing equations are obtained for a wide range of the governing parameters. A good agreement is observed between the analytical model and the numerical simulations.     

水动力作用下低速射流碎化的数值模拟

为研究射流在水动力作用下的碎化特性,采用有限体积法对轴对称坐标下Navier-Stokes方程进行了求解,考虑重力和表面张力的影响,并通过Volume-of-Fluid法与Level-Set法成功捕捉到界面的不稳定发展、变形及射流碎化过程,分析了流场内部速度场和压力场分布,结果表明,射流碎化长度随Re/We“5数呈指数型增加,最后探讨了射流速度、直径及周围流体密度、粘性等参量对射流的碎化过程的影响规律.     

An Euler system source term that develops prototype Z‐pinch implosions intended for the evaluation of shock‐hydro methods

In this paper, a phenomenological model for a magnetic drive source term for the momentum and total energy equations of the Euler system is described. This body force term is designed to produce a Z‐pinch like implosion that can be used in the development and evaluation of shock‐hydrodynamics algorithms that are intended to be used in Z‐pinch simulations. The model uses a J × B Lorentz force, motivated by a 0‐D analysis of a thin shell (or liner implosion), as a source term in the equations and allows for arbitrary current drives to be simulated. An extension that would include the multi‐physics aspects of a proposed combined radiation hydrodynamics (rad‐hydro) capability is also discussed. The specific class of prototype problems that are developed is intended to illustrate aspects of liner implosions into a near vacuum and with idealized pre‐fill plasma effects. In this work, a high‐resolution flux‐corrected‐transport method implemented on structured overlapping meshes is used to demonstrate the application of such a model to these idealized shock‐hydrodynamic studies. The presented results include an asymptotic solution based on a limiting‐case thin‐shell analytical approximation in both (x, y) and (r, z). Additionally, a set of more realistic implosion problems that include density profiles approximating plasma pre‐fill and a set of perturbed liner geometries that excite a hydro‐magnetic like Rayleigh–Taylor instability in the implosion dynamics are demonstrated. Finally, as a demonstration of including and evaluating multiphysics effects in the Euler system, a simple radiation model is included and self‐convergence results for two types of (r, z) implosions are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite‐element/level‐set/operator‐splitting (FELSOS) approach for computing two‐fluid unsteady flows with free moving interfaces

The present work is devoted to the study on unsteady flows of two immiscible viscous fluids separated by free moving interface. Our goal is to elaborate a unified strategy for numerical modelling of two‐fluid interfacial flows, having in mind possible interface topology changes (like merger or break‐up) and realistically wide ranges for physical parameters of the problem. The proposed computational approach essentially relies on three basic components: the finite element method for spatial approximation, the operator‐splitting for temporal discretization and the level‐set method for interface representation. We show that the finite element implementation of the level‐set approach brings some additional benefits as compared to the standard, finite difference level‐set realizations. In particular, the use of finite elements permits to localize the interface precisely, without introducing any artificial parameters like the interface thickness; it also allows to maintain the second‐order accuracy of the interface normal, curvature and mass conservation. The operator‐splitting makes it possible to separate all major difficulties of the problem and enables us to implement the equal‐order interpolation for the velocity and pressure. Diverse numerical examples including simulations of bubble dynamics, bifurcating jet flow and Rayleigh–Taylor instability are presented to validate the computational method. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.