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.     

A technique of flux reconstruction at the interfaces of nonconforming grids for aeroacoustic simulations

A flux reconstruction technique is presented to perform aeroacoustic computations using implicit high-order spatial schemes on multiblock structured grids with nonconforming interfaces. The use of such grids, with mesh spacing discontinuities across the block interfaces, eases local mesh refinements, simplifies the mesh generation process, and thus facilitates the computation of turbulent flows. In this work, the spatial discretization consists of sixth-order finite-volume implicit schemes with low-dispersion and low-dissipation properties. The flux reconstruction is based on the combination of noncentered schemes with local interpolations to define ghost cells and compute flux values at the grid interfaces. The flow variables in the ghost cells are calculated from the flow field in the grid cells using a meshless interpolation with radial basis functions. In this study, the flux reconstruction is applied to both plane and curved nonconforming interfaces. The performance of the method is first evaluated by performing two-dimensional simulations of the propagation of an acoustic pulse and of the convection of a vortex on Cartesian and wavy grids. No significant spurious noise is produced at the grid interfaces. The applicability of the flux reconstruction to a three-dimensional computation is then demonstrated by simulating a jet at a Mach number of 0.9 and a diameter-based Reynolds number of 4×10⁵ on a Cartesian grid. The nonconforming grid interface located downstream of the jet potential core does not appreciably affect the flow development and the jet sound field, while reducing the number of mesh points by a factor of approximately two.     

On The Penetration Of Liquid Contaminants Into Aquifers

The hydrodynamic aspects of the propagation of a liquid non-aqueous contaminant in aquifers and the development of a zone of primary contamination are considered. Special attention is paid to the interaction of the contaminant front with the interfaces between layers of different permeability as considerable lateral spreading of the jet of propagating contaminant occurs at such interfaces.     

Isothermal mass transfer in annular liquid jets with Sievert's solubility law

Summary A study of isothermal gas absorption by underpressurized, axisymmetric, thin, inviscid, incompressible, annular liquid jets which form enclosed volumes, where hazardous wastes may be burned, is presented. The study considers the nonlinear dynamical coupling between the fluid dynamics of, and the gases enclosed by, the annular liquid jet. It assumes equilibrium conditions at the interfaces, and employs Sievert's solubility law to determine the gas concentration at the gas-liquid interfaces. Both steady-state and transient conditions are considered. Under steady-state conditions, the fluid dynamics and mass transfer phenomena are uncoupled, and the rate of generation of combustion gases is equal to the mass absorption rate by the liquid. The transient behaviour of the annular jet is determined from initial conditions corresponding to steady-state operation, once there is no gas generation by the combustion of hazardous wastes. It is shown that, for most of the conditions considered in this paper, there is no leakage of gaseous combustion products through the jet's outer interface, and that the amount of gases dissolved in the liquid at the nozzle exit and the solubility ratio play a paramount role in determining the mass fluxes of hazardous combustion products at the annular jet's interfaces.The research reported in this paper was supported by Project PB91-0767 from the C.I.C.Y.T. of spain.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。     

Scalar interfaces in digital images of turbulent flows

A scalar interface is defined as the surface separating the scalar-marked regions of a turbulent flow from the rest. The problem of determining the two-dimensional intersections of scalar interfaces is examined, taking as a specific example digital images of an axisymmetric jet visualized by laser-induced fluorescence. The usefulness of gradient and Laplacian techniques for this purpose is assessed, and it is shown that setting a proper threshold on the pixel intensity works well if the signal/noise ratio is high. Two methods of determining the proper threshold are presented, and the results are discussed. As one application of the technique, the fractal dimension of the scalar interface is calculated.     

激波致伤的流体力学模型

本文从激波和人体组织相互作用的流体力学现象出发,提出和分析了激波致伤的四个流体动力学模型,即激波与不同密度的物质相互作用模型、气泡坍缩模型。激波作用下组织的表面形成射流模型和激波加速人体形成界面不稳定性模型。这几种模型是在激波和人体器官相互作用的复杂过程中引起原发性冲击伤的重要机制。     

铝-铝爆炸焊接界面的实验研究

对铝-铝同质金属爆炸焊接进行系列实验,研究不同装药比条件下焊接界面的变化特征。所得试件的金相照片表明,焊接界面均具有宏观波状,但其细观形貌又不同于文献报导的连续性波状界面。本次实验得到的三种界面随装药不同而分别呈正弦曲线型、点线型及纯点型,且后两者不再连续。最后利用爆轰波与金属壳表面相互作用理论进行分析,认为焊接界面细观形貌由射流厚度及堆积特征决定,其中界面连续性类似于水流连续性对流量的依赖,当射流厚度不足时难以形成连续界面。     

A hybrid finite-element–volume-of-fluid method for simulating free surface flows and interfaces

A numerical technique is developed for the simulation of free surface flows and interfaces. This technique combines the strength on the finite element method (FEM) in calculating the field variables for a deforming boundary and the versatility of the volume-of-fluid (VOF) technique in advection of the fluid interfaces. The advantage of the VOF technique is that it allows the simulation of interfaces with large deformations, including surface merging and breaking. However, its disadantage is that is solving the flow equations, it cannot resolve interfaces smaller than the cell size, since information on the subgrid scale is lost. Therefore the accuracy of the interface reconstruction and the treatment of the boundary conditions (i.e. viscous stresses and surface tension forces) become grid-size-dependent. On the other hand, the FEM with deforming interface mesh allows accurate implementation of the boundary conditions, but it cannot handle large surface deformations occurring in breaking and merging of liquid regions. Combining the two methods into a hybrid FEM-VOF method eliminates the major shortcomings of both. The outcome is a technique which can handle large surface deformations with accurate treatment of the boundary conditions. For illustration, two computational examples are presented, namely the instability and break-up of a capillary jet and the coalescence collision of two liquid drops.     

Mach number effects on shock-bubble interaction

An investigation of Mach number effects on the interaction of a shock wave with a cylindrical bubble, is presented. We have conducted simulations with the Euler equations for various incident shock Mach numbers () in the range of , using high-resolution Godunov-type methods and an implicit solver. Our results are found in a very good agreement with previous investigations and further reveal additional gasdynamic features with increasing the Mach number. At higher Mach numbers larger deformations of the bubble occur and a secondary-reflected shock wave arises upstream of the bubble. Negative vorticity forms at all Mach numbers, but the “c-shaped” vortical structure appeared at gives its place to a circular-shaped structure at higher Mach numbers. The computations reveal that the (instantaneous) displacements of the upstream, downstream and jet interfaces are not significantly affected by the incident Mach number for values (approximately) greater than . With increasing the incident Mach number, the speed of the jet (arising from the centre of the bubble during the interaction) also increases. Received 21 December 2000 / Accepted 23 April 2001     

激波诱导环形SF6气柱演化的机理

基于可压缩多组分Navier-Stokes控制方程,结合5阶加权本质无振荡格式以及网格自适应加密技术和level-set方法,数值模拟了平面激波(Ma=1.23)与环形SF₆气柱(内外半径分别为8和17.5 mm)界面的相互作用过程。相比于之前的实验结果,数值模拟结果揭示了入射激波在界面内4次透射过程中的复杂波系结构,观察到透射激波在内部界面传播时形成自由前导折射结构并向自由前导冯诺依曼折射结构转换的波系演变过程;另外,界面内的复杂激波结构诱导内部下游界面上的涡量发生了3次反向;在界面演化后期,内部界面形成的“射流”结构与下游界面相互作用,诱导界面形成一对主涡、一对次级涡以及一个反向“射流”结构。定量分析了环形界面长度、宽度、位移、环量以及混合率的变化情况,结果表明,内部气柱的存在减弱了前期小涡结构合并形成大涡结构过程中对界面高度与长度的影响,同时提高了重质气体与环境气体的混合率。     

用改进的ghost fluid方法模拟强激波与界面的相互作用

为了解决原来的ghost fluid方法在计算强激波和界面相互作用时界面附近出现的速度和压力振荡问题,对原来的ghost fluid方法进行了改进,通过在界面处构造Riemann问题并求出界面的压力和速度,ghost fluid流体的压力和速度分别用界面的压力和速度代替,ghost流体的密度通过熵常数外推得到。改进的ghost fluid保持了原来的ghost fluid的简单性,对一维强激波与气-气、气-液界面的相互作用问题以及射流问题进行了数值计算,得到了分辨率较高的计算结果。     

Mechanism of self-oscillations in a supersonic jet impact onto an obstacle 1. Obstacle with a spike

Results of numerical simulations and experimental investigations of self-oscillations arising in the case of impingement of an overexpanded or underexpanded jet onto an obstacle with a spike are reported. The mechanisms of the emergence and maintaining of self-oscillations for overexpanded and underexpanded jets are elucidated. It is demonstrated that self-oscillations are caused by disturbances in a supersonic jet, which induce mass transfer between the supersonic flow and the region between the shock wave and the obstacle. The feedback is ensured by acoustic waves generated by the radial jet on the obstacle. These waves propagate in the gas surrounding the jet, impinge onto the nozzle exit, and initiate disturbances of the supersonic jet parameters. In the overexpanded jet, these disturbances penetrate into the jet core, where they are amplified in oblique shock waves.     

Computations of breakup modes in laminar compound liquid jets in a coflowing fluid

We present a numerical investigation of breakup modes of an axisymmetric, laminar compound jet of immiscible fluids, which flows in a coflowing immiscible outer fluid. We use a front-tracking/finite difference method to track the unsteady evolution and breakup of the compound jet, which is governed by the Navier–Stokes equations for incompressible Newtonian fluids. Numerical results show that depending on parameters such as the Reynolds number Re (in the range of 5–30) and Weber Number We (in the range of 0.1–0.7), based on the inner jet radius and inner fluid properties, the compound jet can break up into drops in various modes: inner dripping–outer dripping (dripping), inner jetting–outer jetting (jetting), and mixed dripping–jetting. Decreasing Re or increasing We promotes the jetting mode. The transition from dripping to jetting is also strongly affected by the velocity ratios, U₂₁ (intermediate to inner velocities) and U₃₁ (outer to inner velocities). Increasing U₂₁ makes the inner jet thinner and stretches the outer jet and thus promotes jetting. In contrast, increasing U₃₁ thins the outer jet, and thus, when the inner jet is dripping, the outer jet can break up into drops in the mixed dripping–jetting mode. Continuously increasing U₃₁ results in thinning both inner and outer jets and thus produces small drops in the jetting mode. In addition, starting from dripping, a decrease in the interfacial tension ratio of the outer to inner interfaces results in the mixed dripping–jetting and jetting modes. These modes produce various types of drops: simple drops, and compound drops with a single inner drop (single-core compound drops) or a few inner drops (multi-core compound drops).     

Numerical treatment of the steady flow of a liquid compound jet

A numerical treatment of the outflow of a two-layer laminar jet into a non-viscous continuous phase is performed. The dispersed phases (i.e. the central core and the concentric layer) are immiscible, incompressible and Newtonian fluids. The method of solution allows for the simultaneous determination of the shape of both interfaces, as well as of the corresponding velocity profiles. The equations of motion of both phases are obtained in a boundary layer approximation. The pressure jump in the radial direction, owing to forces of interfacial tension, is taken into account. Also studied is how the initial velocity profiles at the nozzle exit and some dimensionless parameters affect the interaction between the primary and secondary flow. Numerical results agree qualitatively with some experimental evidence. The approach can also be employed to predict the flow within a viscous continuous phase.     

Symmetric shapes of contact of a jet with the surface of a heavy liquid

The plane linearized problem of oblique impingement of a weightless jet of an ideal incompressible fluid on the surface of a heavy fluid is considered. Flows are sought with symmetric forms of the contact region. Mathematically we arrive at the problem of the eigenvalues and eigenfunctions of an integral equation; solving this problem we obtain various contact forms. The fundamental result for the infinitesimally thin jet of finite intensity is derived by passing to the limit, yielding a result analagous with the forms of free vibrations of a string. Some results are presented for the problem under consideration in the nonlinear formulation.The two-dimensional problem on (vertical) impingement of a jet on a liquid was solved by Olmstead and Raynor [1]. Some results for oblique impingement of a sufficiently thin, slightly curved jet are presented by Frolov [2], Information on other studies, primarily experimental, is presented in [3].This problem is related to the model of a jet curtain of an air-cushion vehicle; in this regard we note the study of Stepanov [4] in which, in particular, a result is obtained for an infinitesimally thin jet curtain.     

The preheated Airy wall jet

The Airy jet is a wall-bounded flow belonging to the similarity class of the well known free jet but, in contrast to the latter, its far field behavior is an algebraically decaying rotational flow. The velocity and temperature distributions of a preheated Airy jet flowing over an insulated wall are investigated using both analytical and numerical methods, and are compared with those of the classical (preheated) exponentially decaying wall jet. For the same value of the dimensionless skin friction parameter, the maximum of the similar velocity profile of the Airy jet exceeds that of the classical wall jet by approximately 20%. The dimensionless temperature along the insulated wall scales for large values of the Prandtl number with Pr^2/3 for both jets, while for small values of the Prandtl number the temperature scales with Pr^1/3 for the Airy jet and goes to 1 for the classical wall jet.This work is dedicated to Michael B. Glauert who passed away on June 14, 2004     

Flow in a viscous jet escaping through a supersonic nozzle into a semi-infinite ambient space

The problem of an axisymmetric gas flow in a supersonic nozzle and in the jet escaping from the nozzle to a quiescent gas is solved within the framework of Navier-Stokes equations. The calculated pressure distribution is compared with that measured in the jet by a Pitot tube. The influence of the jet pressure ratio, Reynolds number, and half-angle of the supersonic part of the nozzle on nozzle flow and jet flow parameters is studied. It is shown that the distributions of gas-dynamic parameters at the nozzle exit are nonuniform, which affects the jet flow. The flow pattern for an overexpanded jet shows that jet formation begins inside the nozzle because of boundary-layer displacement from the nozzle walls. This result cannot be obtained with the inviscid formulation of the problem.     

Color schlieren deflectometry study of jet mixing: effect of buoyancy and perforation

Mixing of jets is crucial for optimal performance of many industrial applications and there is a need to optimize both nozzle geometry and flow conditions. The present study reports the influence of buoyancy and perforation on mixing between a jet and its environment. Optical techniques are ideal for the study of jet mixing due to their non-intrusive and inertia free properties. The present study gives an account of mixing between helium jet and the ambient fluid using a combination of color schlieren deflectometry and radial tomographic mathematics. Four different perforation sizes have been used and the experiments are performed for Reynolds numbers 21–676 and Richardson numbers 3.27–0.0015. Color schlieren images show distinct influence of perforation and flow conditions (Richardson number). Oxygen concentration and jet width quantify effectiveness of jet mixing. Buoyancy plays an important role in mixing at high Richardson number. Perforation improves jet mixing i.e. there is about 120% increase in jet width and the size of perforation plays an important role.