非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

金属射流发散性理论分析及数值模拟

分析经典射流理论和相关文献,给出了在汇聚点坐标系中金属Cu对碰形成射流的汇聚射流区、无射流激波封锁区、无射流强度封锁区和形成发散射流区域。对金属超声速射流形成的发散性问题进行了理论分析,诠释了Walker基于实验提出射流发散理论,证明了金属射流超声速部分可形成发散射流也可形成汇聚射流,且超过1.23倍体声速的金属射流必定是发散的。最后,应用自编的欧拉计算程序MEPH对金属Cu以不同速度、倾角对碰射流形成过程进行数值模拟,得到了分叉射流、空洞射流和密实稀疏射流等的典型射流发散模式图像,印证了理论分析的结果。     

非对称射流形成的几何理论(一)

两股来流碰撞形成非对称射流问题的定常解具有不确定性，理论上尚未很好求解。从分析特殊射流形成入手，由定常出流所满足的最小动能原理出发，分析来流初始交汇位形和相互作用过程，建立了非对称射流形成的方程封闭条件，提出了非对称射流形成的几何理论。以平面二维不等厚度来流碰撞为例，给出了非对称射流出流方向、厚度随初始位形变化的理论预测曲线，与数值模拟结果吻合较好。     

Numerical investigation of the spray-mesh-turbulence interactions for high-pressure, evaporating sprays at engine conditions

This work presents a numerical methodology to simulate evaporating, high pressure Diesel sprays using the Eulerian-Lagrangian approach. Specific sub-models were developed to describe the liquid spray injection and breakup, and the influence of the liquid jet on the turbulence viscosity in the vicinity of the nozzle. To reduce the computational time and easily solve the problem of the grid dependency, the possibility to dynamically refine the grid where the fuel-air mixing process takes place was also included.The validity of the proposed approach was firstly verified simulating an evaporating spray in a constant-volume vessel at non-reacting conditions. The availability of a large quantity of experimental data allowed us to investigate in detail the effects of grid size, ambient diffusivity and used spray sub-models. In this way, different guidelines were derived for a successful simulation of the fuel-air mixture formation process. Finally, fuel injection and evaporation were simulated in an optical engine geometry and computed mixture fraction distributions were compared with experimental data.     

Multiscale simulations and experiments on water jet atomization

Numerical simulation of primary atomization at high Reynolds number is still a challenging problem. In this work a multiscale approach for the numerical simulation of liquid jet primary atomization is applied, using an Eulerian-Lagrangian coupling. In this approach, an Eulerian volume of fluid (VOF) method, where the Reynolds stresses are closed by a Reynolds stress model is applied to model the global spreading of the liquid jet. The formation of the micro-scale droplets, which are usually smaller than the grid spacing in the computational domain, is modelled by an energy-based sub-grid model. Where the disruptive forces (turbulence and surface pressure) of turbulent eddies near the surface of the jet overcome the capillary forces, droplets are released with the local properties of the corresponding eddies. The dynamics of the generated droplets are modelled using Lagrangian particle tracking (LPT). A numerical coupling between the Eulerian and Lagrangian frames is then established via source terms in conservation equations. As a follow-up study to our investigation in Saeedipour et al. (2016a), the present paper aims at modelling drop formation from liquid jets at high Reynolds numbers in the atomization regime and validating the simulation results against in-house experiments. For this purpose, phase-Doppler anemometry (PDA) was used to measure the droplet size and velocity distributions in sprays produced by water jet breakup at different Reynolds numbers in the atomization regime. The spray properties, such as droplet size spectra, local and global Sauter-mean drop sizes and velocity distributions obtained from the simulations are compared with experiment at various locations with very good agreement.     

A CORRELATION EQUATION FOR CALCULATING INCLINED JET PENETRATION LENGTH IN A GAS-SOLID FLUIDIZED BED

Numerical simulation of gas-solid flow in a two-dimensional fluidized bed with an inclined jet was performed. The numerical model is based on the two-fluid model of gas and solids phase in which the solids constitutive equations are based on the kinetic theory of granular flow. The improved ICE algorithm, which can be used for both low and high-velocity fluid flow, were used to solve the model equations. The mechanism of jet formation was analyzed using both numerical simulations and experiments. The emergence and movement of gas bubbles were captured numerically and experimentally. The influences of jet velocity, nozzle diameter, nozzle inclination and jet position on jet penetration length were obtained. A semi-empirical expression was derived and the parameters were correlated from experimental data. The correlation equation, which can be easily used to obtain the inclined jet penetration length, was compared with our experimental data and published correlation equations.     

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.     

Jet formation and penetration mechanism of W typed shaped charge

Existing classical shaped charges are well known for their longer jets capable of achieving large hole depth to hole diameter ratios in metallic targets. However, in some situations, there arises demand to obtain 1：1 ratio for hole depth to hole diameter which is beyond normal shaped charges capability. A new variant of shape charge, named W typed shape charge （WSC）, is proposed in this paper, which can meet the demand of 1：1 ratio, and is based on the geometry that can produce annular jets upon proper initiation scheme. In this paper, we present formation and penetration results of WSC based on three different schemes. We also show that not all WSC designs can form annular jets, only annularly initiated WSC, which also fulfils the ＂Internal-External Liners Equal-Impulse＂ criterion, has the capability to form annular jet. The experimental and numerical results show that when the ratio between annular initiation ring diameter and the charge diameter is 0.75, an annular jet is formed, which was also supported by high speed photographs performed in vacuum. 2D numerical simulations are performed with indigenously developed simulation software, where Eulerian approach with multi-material interface tracking algorithm is utilized, to find various mechanisms involved during jet formation process. The calculation results are found in good agreement with the experimental results, indicating that the interface treatment algorithm proposed in this paper can not only deal with large deformation problem, but also depict clearly the variation of materials interface. It is especially suitable for simulation of the process from liner collapse to formation of shaped charge jet.     

Drag of a Plate Planing on the Shallow Water with Formation of Waves

The plane problem of the plate planing at a constant velocity on the surface of a heavy, ideal, incompressible, finite-depth fluid is considered. The approximate, depth-independent expression for the force acting on the plate is derived from the linear distribution of the fluid velocity along the plate and the height of the flow stagnation point, without regard for jet formation near the leading edge. In this approximate formulation the plate drag depends on its velocity and the trailing edge immersion and does not depend on the planing angle. Experiments and numerical calculations in the exact formulation are performed in the near-critical flow regimes. It is shown that the wave patterns in the experiments and numerical calculations coincide, the formula for the drag being in agreement with the numerical experiments. An approximate criterion of the formation of waves going away from the plate in the forward direction is proposed.     

Bubble collapse in compressible fluids using a spectral element marker particle method. Part 1. Newtonian fluids

This paper is concerned with the development of a high‐order numerical scheme for the modelling of two‐phase Newtonian flows. The companion paper, herein referred to as Part 2, extends the scheme to two‐phase viscoelastic flows. The particular problem of the collapse of a two‐dimensional bubble in the vicinity of a rigid boundary is considered. The governing equations are discretized using the spectral element method, and the two phases are modelled using a marker particle method. The marker particle scheme is validated using the Zalesak slotted disk rotation test problem. A comprehensive set of results is presented for the problem of bubble collapse near a rigid wall, and qualitative agreement is obtained with other numerical studies and experimental observations. Viscous effects are shown to inhibit bubble collapse and prevent jet formation and are therefore likely to have a mitigating effect on cavitation damage.Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical and Experimental Investigation of the Influence of the Wake Behind an Injector Frame on Jet Dilution in a Crossflow

In this paper, a mixing of gases through square Jets issuing normally Into a CrossFlow (JICF) is investigated by means of both numerical simulation and experiment. The jets are emitted by two injectors mounted at the top and bottom of an Injector Frame (IF) which is installed at the center of an Eiffel type wind-tunnel. This jet configuration makes it possible to approximate an industrial gas mixer placed at the center of a pipe. Large Eddy Simulation based on the Smagorinsky model is used, enabling characterization of the mean and fluctuating velocities as well as the oscillating flow frequencies. Different diagnostic techniques, such as Laser Doppler Anemometry and Particle Image Velocimetry are employed for validating the numerical models, and a good agreement between prediction and experiment is obtained. In the numerical simulation, introduction of a passive scalar through the jet makes it possible to show three dilution phenomena. They are generated respectively by the wake of the IF, the jet/wake assemblage and the jets alone in function of the momentum flux ratio between jet and crossflow. Influence of the various parameters on the mixing process between the jets and the crossflow is identified. The numerical results show that if the IF wake is suppressed with the presence of a trailing edge behind the IF, classical formation of Counter-rotating Vortex Pair is found.     

Vortical structures behind a transverse jet in a supersonic flow at high jet to crossflow pressure ratios

A three-dimensional supersonic turbulent flow with symmetric normal injection of circular jets from the channel walls is numerically simulated. The initial Favre-averaged Navier–Stokes equations closed by the k–ω turbulence model are solved by an algorithm based on an ENO scheme. The mechanism of the formation of vortical structures due to the interaction of the jet with the free stream is studied for jet to crossflow total pressure ratios ranging from 3 to 50. It is known from experiments reported in the literature that, for n ? 10, mixing of the jet with the high-velocity flow leads to the formation of a pair of vortices and of an additional separation zone near the wall behind the jet. It is demonstrated that the present numerical results are consistent with such findings and that the pressure distribution on the wall ahead of the jet in the plane of symmetry is also in reasonable agreement with available experimental data.     

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.     

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

Results of the numerical solution of the problem of impingement of an overexpanded supersonic jet onto an obstacle are reported. The mass-flow-rate mechanism of self-oscillations is revealed. This mechanism consists of periodic changes in the regimes of gas inflow and outflow from the separation region to the jet around this region. It is shown that the shock-wave structure of the impinging supersonic jet exerts a significant effect on the amplitude of self-oscillations.     

Filling of channels with a non-Newtonian fluid in the gravity field

The flow of a pseudoplastic fluid with a free surface realized during the filling of a channel at a constant flow rate in the entry section is modeled. The fluid flow direction coincides with that of gravity. The mathematical problem is formulated on the basis of the complete equations of motion. A numerical technique for solving the problem in the two-dimensional formulation is developed on the basis of the SIMPLE algorithm and the method of invariants. A parametric investigation of the free surface evolution is carried out for different values of the governing criteria and rheological parameters. Three filling regimes are found to exist, namely, a complete filling regime, a transitional regime characterized by air cavity formation, and a jet regime. Criterial dependences separating these regimes are presented.     

Numerical simulations of bouncing jets

Bouncing jets are fascinating phenomenon occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non‐Newtonian fluids when the jets fall in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier–Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is carried out with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing. Copyright © 2015 John Wiley & Sons, Ltd.     

On the problem of funnel formation during liquid outflow from vessels

The process of funnel formation during unsteady outflow of initially still water under gravity through one or two circular holes in the bottom of a prismatic vessel rectangular in plan is analyzed experimentally and numerically. The process is studied for various locations of the sinkholes, for vessel bottoms and walls with areas that differ in roughness, and for several values of the initial water depth h ₀. A sufficient, but not necessary, condition for the absence of vortex funnels is formulated. On the basis of a method proposed for determining the fluid rotation direction in vortex funnels, regions are found such that above sinkholes located in them vortex funnels with a given rotation direction are formed, as well as lines such that for sinkhole centers positioned on them stable vortex funnels are not formed. These results are confirmed by numerical calculations and experiments. At different times t the structure of the circumferential velocity field in the sinkhole cross-section is investigated. In the sinkhole crosssection, at several levels above it, and in the outflowing jet the dependence of the integral moment of the circumferential velocities on time and the initial depth is found. The time dependences of the flow-rate characteristics (flow-rate coefficient and volume flow-rate) are obtained for various initial depths and compared with results known in hydraulics. On the range of determining parameters investigated it is established that the fluid depth in the vessel at which a vortex funnel begins to be formed is independent of the initial water depth. The numerical calculations were carried out using STAR-CD software (license of the Institute of Mechanics of Moscow State University) tested in experiments on water outflow from prismatic vessels circular and rectangular in plan.     

Surface instability and primary atomization characteristics of straight liquid jet sprays

Using the detailed numerical simulation data of primary atomization, the liquid surface instability development that leads to atomization is characterized. The numerical results are compared with a theoretical analysis of liquid–gas layer for a parameter range close to high-speed Diesel jet fuel injection. For intermittent and short-duration Diesel injection, the aerodynamic surface interaction and transient head formation play an important role. The present numerical setting excludes nozzle disturbances to primarily investigate this interfacial instability mechanism and the role of jet head. The first disturbed area is the jet head region, and the generated disturbances are fed into the upstream region through the gas phase. This leads to the viscous boundary layer instability development on the liquid jet core. By temporal tracking of surface pattern development including the phase velocity and stability regime and by the visualization of vortex structures near the boundary layer region, it is suggested that the instability mode is the Tollmien–Schlichting (TS) mode similar to the turbulent transition of solid-wall boundary layer. It is also demonstrated that the jet head and the liquid core play an interacting role, thus the jet head cannot be neglected in Diesel injection. In this study, this type of boundary layer instability has been demonstrated as a possible mechanism of primary atomization, especially for high-speed straight liquid jets. The effect of nozzle turbulence is a challenging but important issue, and it should be examined in the future.     

Solution to the problem of a swirling jet from an annular orifice

The problem of the propagation of a laminar immersed fan jet with swirling was considered in [1–3]. In [1], the jet source scheme was used to find a self-similar solution for a weakly swirling jet. An attempt to solve by an integral method the analogous problem for a jet emanating from a slit of finite size was made in [2]. In [3], the equations of motion for a jet with arbitrary swirling were reduced under a number of assumptions to the equations that describe the flow of a flat immersed jet. This paper gives the numerical solution to the problem of the propagation of a radial jet emanating with arbitrary swirling from a slit of finite size and an analytic solution for the main section of the jet.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 49–54, March–April, 1991.