FORCE FIELDS ON INVISCID,SLENDER, ANNULAR LIQUID JETS

Regular perturbation expansions are used to analyse the fluid dynamics of unsteady, inviscid, slender, thin, incompressible (constant density), axisymmetric, upward and downward, annular liquid jets subjected to non-homogeneous, conservative body forces when both the annular jets are very thin and the gases enclosed by and surrounding the jet are dynamically passive. Both inertia- and capillarity-dominated annular jets are considered. It is shown that, for inertia-dominated jets, closure of the leading-order equations is achieved at second order in the perturbation parameter, which is the slenderness ratio, whereas closure is achieved at first order for capillarity-dominated jets. The steady leading-order equations are solved numerically by means of both an adaptive finite difference method which maps the curvilinear geometry of the jet onto a unit square and a fourth-order-accurate Runge–Kutta technique. It is shown that the fluid dynamics of steady, annular liquid jets is very sensitive to the Froude and Weber numbers and nozzle exit angle in the presence of non-homogeneous, conservative body forces. For upward jets with inwardly or axially directed velocities at the nozzle exit the effect of the non-homogeneous, conservative body forces is to increase the leading-order axial velocity component, decrease the jet's mean radius and move the stagnation point downstream. For downward jets with radially outward velocity at the nozzle exit the axial velocity component decreases monotonically as the magnitude of the non-homogeneous, conservative body forces is increased.     

On the role of the liquid flow characteristics on low-Weber-number atomization processes

At low Weber numbers, the aerodynamic forces due to the interaction between gas and liquid do not influence liquid atomization processes. In these situations, atomization processes depend on issuing liquid flow characteristics only. According to the literature, the atomization efficiency is best when the issuing liquid flow shows a high turbulence level. Some injectors are based on this concept and promote the production of turbulence by imposing deflection of the flow inside the nozzle. However, many studies indicate that the level of turbulence does not solely control the atomization efficiency. By conducting a numerical and experimental study on the behavior of cavity nozzles, it is found that internal flow deflection to produce turbulence also produces a non-axial flow component at the nozzle exit whose effect on the atomization process is of paramount importance. Indeed, the results show that the surface energy produced during the atomization process is linearly dependent on the sum of the turbulent kinetic energy and the non-axial kinetic energy at the nozzle exit. This sum represents the energy available for the atomization process, and the influence of the injection pressure as well as of the nozzle geometry on this energy is investigated.     

A cross-correlation velocimetry technique for breakup of an annular liquid sheet

A novel experimental technique for studying the stability of the breakup of a two-phase flow is presented. High-speed backlit photography is used to capture realisations of the unstable mixing layer, and the edge velocity is derived in order to measure the effects of parameters such as gas/liquid momentum ratio, Reynolds and Weber Number. This has traditionally been an edge detection problem that introduces additional uncertainty. Here, a cross-correlation solution is presented, which overcomes the limitations of threshold techniques. Practical application is demonstrated for an atomising annular liquid sheet under several conditions. Sensitivity due to edge blurring and noise is quantified by artificial analysis. Sensitivity analysis shows accuracy and precision to permit sub-pixel precise velocity and stability measurements up to 0.6 sheet thicknesses from the nozzle exit at the conditions studied.     

Numerical investigation of a perturbed swirling annular two-phase jet

A swirling annular gas–liquid two-phase jet flow system has been investigated by solving the compressible, time-dependent, non-dimensional Navier–Stokes equations using highly accurate numerical methods. The mathematical formulation for the flow system is based on an Eulerian approach with mixed-fluid treatment while an adjusted volume of fluid method is utilised to account for the gas compressibility. Surface tension effects are captured by a continuum surface force model. Swirling motion is applied at the inlet while a small helical perturbation is also applied to initiate the instability. Three-dimensional spatial direct numerical simulation has been performed with parallelisation of the code based on domain decomposition. The results show that the flow is characterised by a geometrical recirculation zone adjacent to the nozzle exit and by a central recirculation zone further downstream. Swirl enhances the flow instability and vorticity and promotes liquid dispersion in the cross-streamwise directions. A dynamic precessing vortex core is developed demonstrating that the growth of such a vortex in annular configurations can be initiated even at low swirl numbers, in agreement with experimental findings. Analysis of the averaged results revealed the existence of a geometrical recirculation zone and a swirl induced central recirculation zone in the flow field.     

Oscillating annular liquid membranes

Summary The response of annular liquid membranes to sinusoidal mass flow rate fluctuations at the nozzle exit is analyzed as a function of the amplitude and frequency of the axial velocity fluctuations at the nozzle exit and thermodynamic compression of the gas enclosed by the membrane. It is shown that both the pressure of the gases enclosed by the annular membrane and the axial distance at which the annular membrane merges on the symmetry axis are periodic functions of time which have the same period as that of the mass flow rate fluctuations at the nozzle exit. They are also nearly sinusoidal functions of time for small amplitudes of the mass flow rate fluctuations at the nozzle exit, and exhibit delay and lag times with respect to the sinusoidal axial velocity fluctuations at the nozzle exit. Both the delay and the lag times are functions of the amplitude and frequency of the mass flow rate fluctuations at the nozzle exit and the polytropic exponent. The amplitudes of both the pressure of the gases enclosed by the annular liquid membrane and the convergence length increase and decrease, resp., as the amplitude and frequency of the mass flow rate fluctuations at the nozzle exit, resp., are increased. They also increase as the polytropic exponent is increased.

---

Oszilierende rotationsschalenförmige Flüssigkeitsmembranen

Übersicht Untersucht wird das Verhalten rotationsschalenförmiger Flüssigkeitsmembranen in Abhängigkeit von der Amplitude und Frequenz der axialen Geschwindigkeitsschwankungen an der Düsenmündung und der thermodynamischen Verdichtung des eingeschlossenen Gases, wenn sich der Massestrom an der Düsenmündung sinusförmig ändert. Es wird gezeigt, daß der Druck des eingeschlossenen Gases und der axiale Mündungsabstand des Scheitels der geschlossenen Membran periodische Zeitfunktionen mit der Frequenz der Masseflußschwankung am Düsenaustritt sind. Für kleine Amplituden des Massestroms ist ihr Zeitverhalten ebenfalls fast-sinusförmig, wobei sie bezüglich der sinusförmigen axialen Geschwindigkeitsschwankungen an der Düsenmündung eine Ansprechzeit und Phasenverschiebung aufweisen. Ansprechzeit und Phasenverschiebung sind Funktionen von Amplitude und Frequenz der Massestromschwankung sowie des polytropen Exponenten. Die Amplitude von Gasdruck und Abstand des Membranscheitels von der Düse wächst bzw. fällt mit wachsender Amplitude und Frequenz des Massestroms. Beide nehmen außerdem mit dem polytropen Exponenten zu.

     

Fluid dynamics of slender,thin, annular liquid jets

Perturbation methods are used to obtain the one-dimensional, asymptotic equations that govern the fluid dynamics of slender, thin, inviscid, incompressible, axisymmetric, irrotational, annular liquid jets from the Euler equations. It is shown that, depending on the magnitude of the Weber number, two flow regimes are possible: an inertia-dominated one corresponding to large Weber numbers, and a capillary regime for Weber numbers of the order of unity. The steady equations governing these two regimes have analytical solutions for the liquid's axial velocity component and require a numerical integration to determine the jet's mean radius for inertia-dominated jets. The one-dimensional equations derived in this paper are shown to be particular cases of a hydraulic model for annular liquid jets, and this model is used to determine the effects of gravity modulation on the unsteady fluid dynamics of annular liquid jets in the absence of mass injection into the volume enclosed by the jet and mass absorption. It is shown that both the convergence length and the pressure coefficient are periodic functions of time which have the same period as that of the gravity modulation, but undergo large variations as the amplitude, frequency and width of gravitational pulses is varied.     

Experimental analysis of liquid–gas interface at low Weber number: interface length and fractal dimension

The present paper reports an experimental investigation on atomizing liquid flows produced by simplified cavity nozzles. The Weber number being kept low, the sprays produced by these injectors depend on the liquid flow characteristics only, and more precisely, on the non-axial kinetic energy and of the turbulent kinetic energy at the nozzle exit. The investigation reported here concentrates on the characterization of liquid flows during atomization by measuring the spatial variation of the local interface length and of the local interface fractal dimension. Both parameters were found representative of the physics of atomization process: they depend on the characteristics of the flow issuing from the nozzle and they are related to the subsequent drop size distribution. The local interface length is representative of the amount of liquid–gas interface surface area, and is a function of both the non-axial and the turbulent kinetic energies at the nozzle exit. The fractal dimension is representative of the tortuosity of the liquid–gas interface and, as expected, is mainly related to the turbulent kinetic energy at the nozzle exit. As far as the drop size distribution is concerned, it is found that the local interface length at the instant of break-up determines a representative drop diameter of some kind, whereas the fractal dimension at the same instant controls the dispersion of the distribution.     

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.     

High-speed visualisation of primary break-up of an annular liquid sheet

In this experimental study, a thin annular moving water sheet is placed between two annular co-flowing air streams. The shear at the interface gives rise to Kelvin–Helmholtz type instabilities and promotes development of a sinuous surface wave at the gas–liquid interface. The amplitude of the surface wave is amplified as it travels downstream of the nozzle exit until it ruptures forming spanwise and streamwise ligaments. The liquid sheet is illuminated with high-powered halogen lamps. High-speed imaging is used in this study to qualitatively visualise the structure of the spray—of particular interest is the evolution of the spray into a ligament structure during the primary break-up and the role the outer air stream plays in this process. Sequences of images with high temporal resolution (∼2,000 fps) are recorded for image processing and analysis of the surface waves and ligament formation. A preliminary analysis of the waveform of the outer gas–liquid interface of the annular liquid sheet over a range of conditions shows the sheet Strouhal number to increase with increasing gas to liquid momentum ratio. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Free Surface of a High Speed Capillary Jet

A free surface shape of a viscous liquid jet is investigated at large Reynolds and Weber numbers. The jet is ejected into a vacuum from a cylindrical nozzle with a flat exterior surface. The liquid is completely wetting the nozzle material (zero contact angle). Free jet surface is non-cylindrical near the nozzle. There is a smooth connection between the flat external surface of the nozzle and the cylindrical surface of the jet away from the nozzle. The size of the connection region is estimated by means of the boundary layer technique.     

Swirling and Impinging Effects in an Annular Nonpremixed Jet Flame

The effects of swirl and downstream wall confinement on an annular nonpremixed flame were investigated using direct numerical simulation (DNS). Fully three-dimensional parallel DNS was performed employing high-order numerical methods and high-fidelity boundary conditions to solve governing equations for variable-density flow and finite-rate Arrhenius chemistry. Three swirl numbers have been examined: 0 (without swirl), 0.4 and 0.8, while the effects of downstream wall confinement have been examined for swirl numbers 0 and 0.4. Results have been presented in terms of instantaneous and time-averaged flow quantities, which have also been analysed using energy spectra and proper orthogonal decomposition (POD). Effects of swirl on the fluid dynamic behaviour of the annular nonpremixed flame were found to be significant. The fluid dynamic behaviour of the flame is greatly affected by the interaction between the geometrical recirculation zone (GRZ) near the jet nozzle exit due to the annular configuration, the central recirculation zone (CRZ) associated with swirl, the unsteady vortical structures in the jet column due to the shear instability, and the downstream wall confinement. Depending on the degree of swirl, the GRZ near the burner mouth and the CRZ may co-exist or one zone may be overwhelmed by another. At a moderate swirl number, the co-existence leads to a flame with strong reaction attached to the burner mouth; while at a high swirl number, the CRZ dominates over the GRZ. The precessing vortex core was observed to exist in the swirling flow fields. The Nusselt number distribution of the annular impinging flames differs from that of round impinging jets. The POD analysis revealed that wall effects on the flow field are mainly associated with the higher mode numbers.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Experimental and analytical investigation of liquid sheet breakup characteristics

The main objective of this research is to study analytically and experimentally the liquid sheet breakup of a flat fan jet nozzle resulting from pressure-swirling. In this study the effects of nozzle shape and spray pressure on the liquid sheet characteristics were investigated for four nozzles with different exit widths (1.0, 1.5, 2.0 and 2.5 mm). The length of liquid sheet breakup, liquid sheet velocity and the size of formed droplets were measured by a digital high speed camera. The breakup characteristics of plane liquid sheets in atmosphere are analytically investigated by means of linear and nonlinear hydrodynamic instability analyses. The liquid sheet breakup process was studied for initial sinuous and also varicose modes of disturbance. The results presented the effect of the nozzle width and the spray pressure on the breakup length and also on the size of the formed droplets. Comparing the experimental results with the theoretical ones for all the four types of nozzles, gives a good agreement with difference ranges from 4% to 12%. Also, the comparison between the obtained results and the results due to others shows a good agreement with difference ranged from 5% to 16%. Empirical correlations have been deduced describing the relation between the liquid sheet breakup characteristics and affecting parameters; liquid sheet Reynolds number, Weber number and the nozzle width.     

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.     

Shock wave configurations and flow structures in non-axisymmetric underexpanded sonic jets

An experimental and numerical study of underexpanded free sonic jet flows issuing from rectangular, elliptical and slot nozzles has been undertaken. Aspect ratios (AR) of 1, 2, and 4 are described at pressure ratios (exit plane pressure to ambient pressure), of 2 and 3. There is good qualitative agreement between the experimental observations and the numerical predictions. In the case of rectangular jets, a complex system of shock waves forming the incident shock system is identified. This shock wave system originates at the corners of the nozzle exits, and proceeds downstream. Mach reflections are found to occur on the incident shock wave surface as well as the presence of a Mach disk terminating the first jet cell. This Mach disk has the shape of a square, a hexagon, or an octagon depending on the nozzle shape. For slot and elliptical jets, the formation of the incident shock wave was not observed along the minor axis plane of the nozzle for AR > 2. The incident shock wave was observed to originate downstream of the nozzle exit in the major axis plane. This wave system undergoes a transition to Mach reflection as it propagates downstream of the nozzle exit. In all cases tested, the shape of the jet boundary is significantly distorted. In rectangular jets, the narrowing of the jet boundary along the diagonal axis of the nozzle exit is observed, and in the case of the elliptical and slot jets axis switching is noted.     

Air flow and pressure inside a pressure-swirl spray and their effects on spray development

Air flow and pressure inside a pressure-swirl spray for direct injection (DI) gasoline engines and their effects on spray development have been analyzed at different injector operating conditions. A simulation tool was utilized and the static air pressure at the centerline of the spray was measured to investigate the static pressure and flow structure inside the swirl spray. To investigate the effect of static air pressure on swirl spray development, a liquid film model was applied and the Mie-scattered images were captured. The simulation and experiment showed that recirculation vortex and air pressure drop inside the swirl spray were observable and the air pressure drop was greater at high injection pressure. At high fuel temperature, the air pressure at the nozzle exit showed higher value compared to the atmospheric pressure and then continuously decreased up to few millimeters distance from the nozzle exit. The pressure drop at high fuel temperatures was more than that of atmospheric temperature. This reduced air pressure was recovered to the atmospheric pressure at further downstream. The results from the liquid film model and macroscopic spray images showed that the air pressure started to affect the liquid film trajectory about 3 mm from the nozzle exit and this effect was sustained until the air pressure recovered to the atmospheric pressure. However, the entrained air motion and droplet size have more significant influence on the spray development after the most of the liquid sheet is broken-up and the spray loses its initial momentum.     

Investigation of high frequency noise generation in the near-nozzle region of a jet using large eddy simulation

This paper reports on the simulation of the near-nozzle region of an isothermal Mach 0.6 jet at a Reynolds number of 100,000 exhausting from a round nozzle geometry. The flow inside the nozzle and the free jet outside the nozzle are computed simultaneously by a high-order accurate, multi-block, large eddy simulation (LES) code with overset grid capability. The total number of grid points at which the governing equations are solved is about 50 million. The main emphasis of the simulation is to capture the high frequency noise generation that takes place in the shear layers of the jet within the first few diameters downstream of the nozzle exit. Although we have attempted to generate fully turbulent boundary layers inside the nozzle by means of a special turbulent inflow generation procedure, an analysis of the simulation results supports the fact that the state of the nozzle exit boundary layer should be characterized as transitional rather than fully turbulent. This is believed to be most likely due to imperfections in the inflow generation method. Details of the computational methodology are presented together with an analysis of the simulation results. A comparison of the far field noise spectrum in the sideline direction with experimental data at similar flow conditions is also carried out. Additional noise generation due to vortex pairing in the region immediately downstream of the nozzle exit is also observed. In a second simulation, the effect of the nozzle exit boundary layer thickness on the vortex pairing Strouhal frequency (based on nozzle diameter) and its harmonics is demonstrated. The limitations and deficiencies of the present study are identified and discussed. We hope that the lessons learned in this study will help guide future research activities towards resolving the pending issues identified in this work.

---

Presented as AIAA Paper 2006-2499 at 12th AIAA/CEAS Aeroacoustics Conference, 8–10 May 2006, Cambridge, MA, USA.     

Experimental Investigation of Nozzle Exit Reflector Effect on Supersonic Jet

The present study describes an experimental work to investigate the effect of a nozzle exit reflector on a supersonic jet that is discharged from a convergent–divergent nozzle with a design Mach number of 2.0. An annular reflector is installed at the nozzle exit and its diameter is varied. A high-quality spark schlieren optical system is used to visualize detailed jet structures with and without the reflector. Impact pressure measurement using a pitot probe is also carried out to quantify the reflector’s effect on the supersonic jet which is in the range from an over-expanded to a moderately under-expanded state. The results obtained show that for over-expanded jets, the reflector substantially increases the jet spreading rate and reduces the supersonic length of the jet, compared with moderately under-expanded jets. The reflector’s effect appears more significant in imperfectly expanded jets that have strong shock cell structures, but is negligible in correctly expanded jet.     

A study of the flow characteristics in micro-abrasive jets

An experimental study of particle velocities in micro-abrasive jets by using the particle image velocimetry (PIV) technique is presented. It has been found that the particle jet flow has a nearly linear expansion downstream. The particle velocities increase with air pressure, and the increasing rate increases with nozzle diameter within the range considered. The instantaneous velocity profile of the particle flow field in terms of the particle velocity distribution along the axial and radial directions of the jets is discussed. For the axial profile in the jet centerline downstream, there exists an extended acceleration stage, a transition stage, and a deceleration stage. For the radial velocity profiles, a relatively flat shape is observed at a jet cross-section near the nozzle exit. Mathematical models for the particle velocities in the air jet are then developed. It is shown that the results from the models agree well with experimental data in both the variation trend and magnitude.