Swirl number and nozzle confinement effects in a flat-vane axial swirler

This paper presents the results of a parametric experimental study of free swirling flow at the exit of a flat-vane axial swirler. A total of 16 data sets were acquired by combining four swirler vane angles (22°, 29°, 50.5°, and 58.3°) and four exit nozzles of different diameters (30, 40, 52, and 76 mm). Sophisticated pressure probes consisting of precise microphones and a two-component LDV system were used to investigate the effect of these geometrical parameters on swirling flow regimes characterized by the swirl number. Particular attention was paid to the precessing vortex core (PVC) phenomenon observed at the exit of the swirler nozzle. It has been shown that by varying the vane angle and the diameter of the exit nozzle, it is possible to independently control the swirl number value and the occurrence of a PVC. A distinct correlation has been found between the PVC-induced pressure pulsations detected by acoustic probes and the tangential velocity fluctuations measured by LDV. The use of microphones provides a quick way to measure the frequency response of swirl flow in a wide range of geometries and flow configurations. The PVC effect does not occur at low subcritical values of the integral swirl number (S < 0.5) and in the case of strong swirl flow (S_g = 0.9 and 1.2) in the absence of constriction by the nozzle (D_e/D₀ = 1). The disappearance of the PVC effect for strong swirl flow without constriction is due to the extreme displacement of the flow to the nozzle walls. The absence of a PVC in the flow was inferred not only from measurements of the frequency response of the flow over a wide range of Re numbers, but also from the absence of specific markers in velocity RMS distributions. Measurement results are used to derive an empirical correlation of the integral swirl number and the Strouhal number with a modified geometric swirl number. This allows a generalization of the frequency characteristics of swirling flows with a PVC for flat-vane axial swirlers, which are widely used in engineering.     

Fluid Dynamical Analysis of Atmospheric Reacting and Isothermal Swirling Flows

A series of unconfined swirling premixed natural gas/air flames was investigated. Reynolds-numbers spanned from 10,000 to 42,300. Respective isothermal flows were studied additionally to gain insight into changes of fluid dynamical features caused by combustion. Statistical moments, Reynolds-stresses, temporal time scales, spatial length scales, and power spectral densities were deduced from one- and two-point laser Doppler velocimetry (LDV) data. Properties of the turbulent flows and dependencies on Reynolds-number, swirl number, and chemical reactions are discussed. Most distinct differences between combusting and isothermal flows were precessing vortex cores (PVC) occurring only for the latter cases. The study is aimed to serve as a database of a generic flame geometry featuring important characteristics of industrial applications for validation of numerical simulations. Therefore, nozzle exit profiles as important inlet conditions to numerical simulations are thoroughly documented.     

Studies on stability and dynamics of a swirling jet

The temporal instability and nonlinear evolution of the swirling jet near a nozzle exit are studied by both normal-mode method and three dimensional direct numerical simulation (3D DNS). It is found that the swirl enhances the maximum linear growth rates for negative helical modes, while decreases the growth rate for axisymmetric mode. Numerical simulations show that the evolution in early stage is compared well with the linear stability theory. In nonlinear stage, the swirl promotes the breakup of 3D large scale organized structures in the flow into small eddies. The project supported by the National Natural Science Foundation of China (19772052)     

Computation of a drastic flow pattern change in an annular swirling jet caused by a small decrease in inlet swirl

This paper investigates the flow pattern change in an annular jet caused by a sudden change in the level of inlet swirl. The jet geometry consists of an annular channel followed by a specially designed stepped‐conical nozzle, which allows the existence of four different flow patterns as a function of the inlet swirl number. This paper reports on the transition between two of them, called the ‘open jet flow high swirl’ and the ‘Coanda jet flow.’ It is shown that a small sudden decrease of 4% in inlet swirl results in a drastic and irreversible change in flow pattern. The objective of this paper is to reveal the underlying physical mechanisms in this transition by means of numerical simulations. The flow is simulated using the unsteady Reynolds‐averaged Navier–Stokes (URANS) approach for incompressible flow with a Reynolds stress turbulence model. The analysis of the numerical results is based on a study of different forces on a control volume, which consists of the jet boundaries. The analysis of these forces shows that the flow pattern change consists of three different regimes: an immediate response regime, a quasi‐static regime and a Coanda regime. The simulation reveals that the pressure–tangential velocity coupling during the quasi‐static regime and the Coanda effect at the nozzle outlet during the Coanda regime are the driving forces behind the flow pattern change. These physical mechanisms are validated with time‐resolved stereo‐PIV measurements, which confirm the numerical simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental investigation of helical structures in swirling flows

In this investigation the flow in a generic swirl tube with a tangential double-inlet swirl generator and variable exit orifices was experimentally investigated. Using magnetic resonance velocimetry (MRV) three-dimensional, three-component velocity fields were measured for two different Reynolds numbers: 10,000 and 15,000, and for three different exit orifices. The swirl generator had a fixed geometry producing an initial swirl number of 1.6 for all cases. One major observation is the occurrence of a three-layered flow structure. An annular main flow was surrounded by a recirculation zone, as reported in previous literature. However, this recirculation zone – also of an annular shape – exhibited a third layer inside: a thin, high speed jet in the center of the tube with the same flow direction as the main flow. Therefore, the conventional classification of swirling flows into ring and recirculation zone, has to be extended by a core zone. This three-layered flow structure develops independently of the exit configuration. Helical structures were observed in the near-wall region for all cases investigated. Applying an eccentric exit orifice results in the occurrence of strong stationary helical structures not only in the near-wall region but also in the center of the tube. The results, deviating significantly from previous results in the literature, underline the need for more detailed research on the topic of cyclone type flows.     

Formation of Coanda jet from sharp-edged swirl nozzle with base plate

Experiments were conducted using tufts and PIV to determine the conditions for which a swirled gas jet issuing from a sharp-edge nozzle, in flush with a base plate, would form a Coanda jet. The flow field was also simulated. The inception of the Coanda jet was observed to be associated with the formation of a recirculation bubble at the nozzle exit. A threshold value of swirl number, which increased monotonically with Reynolds number, was required for the formation of the Coanda jet. The Coanda jet was associated with hysteresis. The flow features and transition from a diverging jet to a Coanda jet are discussed.     

Analysis of Unsteady Motion with Respect to Noise Sources in a Gas Turbine Combustor: Isothermal Flow Case

In order to evaluate the direct and indirect contributions to the total combustion noise emission, a combustion chamber consisting of a swirl burner and an exit nozzle of Laval-shape, representing a gas turbine combustor, is investigated by means of experiments and large eddy simulation. Focused on the isothermal flow case first and encouraged by a good overall agreement between the LES and the experimental data for the flow field, a first characterisation of the flow with respect to noise sources is performed. To analyse acoustic properties of the flow, time and length scales are evaluated inside the combustor. Furthermore, the evidence for the existence of a precessing vortex core (PVC), typical for configurations with swirl, is revealed. Finally, the effect of the PVC on the flow inside the Laval nozzle is discussed.     

Some features of the flow around rapidly rotating bodies made of cellular-porous materials

Two types of gas flows arising near a rapidly rotating cellular-porous disk are studied numerically and experimentally. Steady-state limits for the flow around a disk rotating in free space and the type and scenario of the loss of stability are determined. Transitional flows are characterized by formation of a vortex sheet at the boundary of the exhausting jet. Numerical simulations of the flow around a cellular-porous disk rotating near a flat screen show that it is possible to form a closed swirl flow responsible for redistribution of swirl in the gap between the disk and the flat screen. The computed results offer an explanation for the experimentally observed excess of tangential velocity of the flow in the gap over the velocity of disk rotation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 86–96, January–February, 2007.     

Numerical predictions of the bifurcation of confined swirling flows

The bifurcation of confined swirling flows was numerically investigated by employing both the k-? and algebraic stress turbulence models. Depending upon the branch solution examined, dual flow patterns were predicted at certain swirl levels. In the lower-branch solution which is obtained by gradually increasing the swirl level from a low-swirl flow, the flow changes with increasing swirl number from the low-swirl flow pattern to a high-swirl flow pattern. In the upper-branch solution which is acquired by gradually decreasing the swirl level from a high-swirl flow, on the other hand, the flow can maintain itself in the high-swirl flow pattern at the swirl levels where it exhibits the low-swirl flow pattern in the lower branch. The bifurcation of confined swirling flows was predicted with either the k-? model or the algebraic stress model being employed. Both the k-? and algebraic stress models result in comparable and sufficiently good predictions for confined swirling flows if high-order numerical schemes are used. The reported poor performance of the k-? model was clarified to be mainly attributable to the occurrence of the bifurcation and the use of low-order numerical schemes.     

Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.     

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.     

Ballistic imaging in the near-field of an effervescent spray

We have investigated liquid breakup mechanisms in the near nozzle region of a high-pressure effervescent atomizer using ballistic imaging. This technique has revealed various breakup regimes depending upon total flow rate and the gas-to-liquid ratio (GLR). At low total speeds, the jet does not exhibit the wide spread angle and rapid breakup for which effervescent sprays are known, even at high GLR. Above a distinct threshold value for total flow rate, the jet passes through several recognizable flow regimes depending on GLR and it does achieve the expected wide spread angle and rapid breakup. Intermediate GLR’s produce interesting flow patterns that seem to be generated by surging at the nozzle exit, and this surging can probably be attributed to the flow pattern just at the nozzle exit. Indeed, specific interior flows seem to generate the most rapid breakup and should be investigated further.     

On coherent structures and mixing characteristics in the near field of a rotating-pipe jet

Mixing characteristics and coherent structures populating the near-nozzle area of a rotating-pipe jet at the Reynolds number of 5300 were studied by Large-eddy simulation (LES). The swirl rate, defined as the ratio of the tangential velocity of the inner pipe wall to the bulk axial velocity, varied from 0 to 1, corresponding to a weak-to-moderate swirl intensity, insufficient to induce reverse flow near the nozzle. The visualization shows that for the non-swirling jet the near-wall streaky structures generated in the pipe interact with the shear layer, evolving into hairpin-like structures that become tilted at low rotation rates. For higher swirl, they cannot be recognized as they are destroyed at the nozzle exit. No large-scale coherent structures akin to Kelvin–Helmholtz vortical rings in the ‘top-hat’ jets are identifiable close to the nozzle. Using the single and joint probability density functions of velocity and passive scalar (temperature) fields we quantify the events responsible for the intensive entrainment at various swirl numbers. The isosurface of the temperature field indicates the meandering and precessing motion of the rotating jet core at the axial distance of 6D downstream, where D is the diameter of the pipe. The Fourier analysis with respect to the azimuthal angle and time reveals an interplay between the co- and counter-rotating modes. These findings explain the previously detected but not fully clarified phenomenon of the weakly counter-rotating jet core at low swirl rates.     

The spray characteristics of a pressure-swirl injector with various exit plane tilts

The gasoline spray characteristics of a pressure-swirl injector were investigated with various exit plane tilts. The analysis focused on the correlation between tilt angle and flow angle. Mie-scattering technique and phase Doppler anemometry were employed to analyze the macroscopic spray development and droplet size distribution of the spray. An analytical method for mass flux estimation was applied to understand the velocity distribution at the nozzle exit. The results showed that the spray shape and velocity distribution of the spray were more asymmetrical at high tilt angles. In particular, an opened hollow cone spray was formed when the tilt angle is greater than the complementary flow angle. The pressure drop inside the spray, one of the crucial factors for the swirl spray collapse at various surrounding conditions, was attenuated in this opened hollow cone spray since the pressure inside the spray was assimilated to the surrounding air pressure. The spray collapse at high fuel temperature and back pressure conditions did not appear when the tilt angle is larger than the complementary flow angle due to the reduced pressure drop inside the spray. However, tilt angle should be optimized to fulfill the requirements of spray robustness and avoid the locally rich area. The droplet size of 70° tilted nozzle spray shows a value similar to that of the original swirl spray in the plane that includes nozzle axis and the major axis of exit surface ellipse (Major axis plane) while it shows an increased value in the plane that includes nozzle axis and the minor axis of exit surface ellipse (Minor axis plane).     

Shock and wave dynamics in cavitating compressible liquid flows in injection nozzles

Due to the exceptional high inlet pressures up to 2,000 bar flow dynamics and efficiency of modern injection systems are controlled by high frequency wave dynamics of the compressible liquid flow. Corresponding to alternating shock and expansion waves the liquid fluid evaporates and recondenses instantaneously. Here we present CFD simulations of the time accurate evolution of cavitating flows in 2-D plane and in six-hole injection nozzles with focus on the wave dynamics just after initialisation of the flow and within the time scale Δt ≤ 10^?4 s of pilot and multi-point injection. Due to shock reflections at the bottom of the sack hole the instantaneous maximum pressure increases more than three times higher as compared with the prescribed pressure at the nozzle inlet. For instance, in case of an inlet pressure of 600 bar the maximum pressure in the sack and therefore ahead of the nozzle bore holes reaches about 2,100 bar. It is quite reasonable that this amplification of the pressure affects the evolution of the convective flow and therefore the mass flow through the nozzle bore holes.     

Effects of confinement,geometry, inlet velocity profile,and Reynolds number on the asymmetry of opposed-jet flows

The opposed-jet counterflow configuration is widely used to measure fundamental flame properties that are essential targets for validating chemical kinetic models. The main and key assumption of the counterflow configuration in laminar flame experiments is that the flow field is steady and quasi-one-dimensional. In this study, experiments and numerical simulations were carried out to investigate the behavior and controlling parameters of counterflowing isothermal air jets for various nozzle designs, Reynolds numbers, and surrounding geometries. The flow field in the jets’ impingement region was analyzed in search of instabilities, asymmetries, and two-dimensional effects that can introduce errors when the data are compared with results of quasi-one-dimensional simulations. The modeling involved transient axisymmetric numerical simulations along with bifurcation analysis, which revealed that when the flow field is confined between walls, local bifurcation occurs, which in turn results in asymmetry, deviation from the one-dimensional assumption, and sensitivity of the flow field structure to boundary conditions and surrounding geometry. Particle image velocimetry was utilized and results revealed that for jets of equal momenta at low Reynolds numbers of the order of 300, the flow field is asymmetric with respect to the middle plane between the nozzles even in the absence of confining walls. The asymmetry was traced to the asymmetric nozzle exit velocity profiles caused by unavoidable imperfections in the nozzle assembly. The asymmetry was not detectable at high Reynolds numbers of the order of 1000 due to the reduced sensitivity of the flow field to boundary conditions. The cases investigated computationally covered a wide range of Reynolds numbers to identify designs that are minimally affected by errors in the experimental procedures or manufacturing imperfections, and the simulations results were used to identify conditions that best conform to the assumptions of quasi-one-dimensional modeling.     

Flows through plane sudden‐expansions

A calculation method has been developed and used to represent flows downstream of plane symmetric expansions with dimensions and velocities encompassing laminar and turbulent flows. Except for very low Reynolds numbers, the flows are time‐dependent and asymmetric and the calculated results are appraised first in relation to published measurements of laminar flows and then to new measurements obtained at a Reynolds number of 26 500. The time‐dependent laminar simulations indicate that the critical Reynolds numbers are predicted with excellent accuracy for different expansion ratios and the details of the asymmetric velocity profiles are in good agreement with experimental measurements. The laminar flow calculations also show that increasing the thickness of the separating boundary layer leads to longer regions of separation and no dominant frequency for Reynolds numbers up to those at which the third separation region was observed. The turbulent flow simulations made use of the k–ε turbulence model and provided a satisfactory representation of measurements, except in regions close to the wall and within the recirculation regions. Also, the longer reattachment length was underestimated. Limitations are discussed in relation to these and higher‐order assumptions. Copyright © 2000 John Wiley & Sons, Ltd.     

Sensitivity of an asymmetric 3D diffuser to vortex-generator induced inlet condition perturbations

Modifications of the turbulent separated flow in an asymmetric three-dimensional diffuser due to inlet condition perturbations were investigated using conventional static pressure measurements and velocity data acquired using magnetic resonance velocimetry (MRV). Previous experiments and simulations revealed a strong sensitivity of the diffuser performance to weak secondary flows in the inlet. The present, more detailed experiments were conducted to obtain a better understanding of this sensitivity. Pressure data were acquired in an airflow apparatus at an inlet Reynolds number of 10,000. The diffuser pressure recovery was strongly affected by a pair of longitudinal vortices injected along one wall of the inlet channel using either dielectric barrier discharge plasma actuators or conventional half-delta wing vortex generators. MRV measurements were obtained in a water flow apparatus at matched Reynolds number for two different cases with passive vortex generators. The first case had a pair of counter-rotating longitudinal vortices embedded in the boundary layer near the center of the expanding wall of the diffuser such that the flow on the outsides of the vortices was directed toward the wall. The MRV data showed that the three-dimensional separation bubble initially grew much slower causing a rapid early reduction in the core flow velocity and a consequent reduction of total pressure losses due to turbulent mixing. This produced a 13% increase in the overall pressure recovery. For the second case, the vortices rotated in the opposite sense, and the image vortices pushed them into the corners. This led to a very rapid initial growth of the separation bubble and formation of strong swirl at the diffuser exit. These changes resulted in a 17% reduction in the overall pressure recovery for this case. The results emphasize the extreme sensitivity of 3D separated flows to weak perturbations.     

Experimental analysis of flashback in lean premixed swirling flames: upstream flame propagation

Simultaneous 10-kHz OH-PLIF and 20-kHz two-component PIV were made in conjunction with wide-field 20-kHz flame luminescence imaging of an unconfined, swirling, lean premixed, bluff-body stabilized flame during flashback. Flashback was induced by increasing the stoichiometry or swirl number or reducing the Reynolds number. A detailed stability regime was prepared and compared to predictions. Analysis of the time-correlated flame history inside the exit nozzle during flashback and non-flashback flame events led to a new hypothesis for the flashback mechanism.