The effects of inlet conditions and expansion ratio on the onset of flow reversal in swirling flow in a sudden expansion

Experiments are reported in which the minimum swirl intensity required to produce a central recirculation zone in a swirling sudden expansion flow is determined as a function of expansion ratio and inlet conditions. Using a swirl generator which allows for independent variation of velocity profile shape and swirl number, it is shown that an inlet tangential velocity distribution with a large solid body vortex core or an axial velocity profile with a maximum on the axis will lead to a higher critical swirl.     

Reynolds stress modelling of jet and swirl interaction inside a gas turbine combustor

Influences of the inlet swirl levels on the interaction between the dilution air jets and the swirling cross‐flow to the interior flow field inside a gas turbine combustor were investigated numerically by Reynolds stress transport model (RSTM). Due to the intense swirl and jet interaction, a high level of swirl momentum is transported to the centreline and hence, an intense vortex core is formed. The strength of the centreline vortex core was found to depend on the inlet swirl levels. For the higher swirling inlet, the decay of the swirling motion causes strong streamline variation of pressure; and consequently leads to an elevated level of deceleration of its axial velocity. Predictions contrasted with measurements indicate that the stress model reproduces the flow correctly and is able to reflect the influences of inlet swirl levels on the interior flow structure. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Large Eddy Simulation of Low Swirl Flames Under External Flow Excitations

Low swirl flame characteristics under external flow excitations are numerically investigated using large eddy simulations with a dynamically thickened flame combustion model. A finite volume scheme on a Cartesian grid with a dynamic one equation eddy viscosity subgrid scale model is used for large eddy simulations. The excitations are imposed on inlet velocity profiles by a sinusoidal forcing function over a wide range of amplitudes and frequencies. Present investigation shows that although, the swirling motion of the low swirl flame is not intense enough to induce a recirculation zone in ensemble averaged results, external flow excitations increase the local swirl number upstream of the flame front. Such increase in the local swirl number is at maximum value when the low swirl flame is excited at the dominant frequency of the flow field, which in turn induces a vortex breakdown and hence a central recirculation zone. The strength and size of the time averaged recirculation zone depend on both the amplitude and frequency of the excitations. Moreover, phase-locked results indicate that external flow excitations induce local swirl fluctuations ahead of the flame front which alter the strength of the recirculation zone at different phase angles of the excitations.     

Dispersion of Entropy Perturbations Transporting through an Industrial Gas Turbine Combustor

In the context of combustion noise and combustion instabilities, the transport of entropy perturbations through highly simplified turbulent flows has received much recent attention. This work performs the first systematic study into the transport of entropy perturbations through a realistic gas turbine combustor flow-field, exhibiting large-scale hydrodynamic flow features in the form of swirl, separation, recirculation zones and vortex cores, these being ubiquitous in real combustor flows. The reacting flow-field is simulated using low Mach number large eddy simulations, with simulations validated by comparison to available experimental data. A generic artificial entropy source, impulsive in time and spatially localized at the flame-front location, is injected. The conservation equation describing entropy transport is simulated, superimposed on the underlying flow-field simulation. It is found that the transport of entropy perturbations is dominated by advection, with both thermal diffusion and viscous production being negligible. It is furthermore found that both the mean flow-field and the large-scale unsteady flow features contribute significantly to advective dispersion — neither can be neglected. The time-variation of entropy perturbation amplitude at combustor exit is well-modelled by a Gaussian profile, whose dispersion exceeds that corresponding to a fully-developed pipe mean flow profile roughly by a factor of three. Finally, despite the attenuation in entropy perturbation amplitude caused by advective dispersion, sufficient entropy perturbation strength is likely to remain at combustor exit for entropy noise to make a meaningful contribution at low frequencies.     

Large Eddy Simulations of Unconfined Non-reacting and Reacting Turbulent Low Swirl Jets

The low swirl flow is a novel method for stabilizing lean premixed combustion to achieve low emissions of nitrogen oxides. Understanding the characteristics of low swirl flows is of both practical and fundamental interest. In this paper, in order to gain better insight into low swirl stabilized combustion, large eddy simulation and dynamically thickened flame combustion modeling are used to characterize various features of non-reacting and reacting low swirl flows including vortex breakdown, shear layers’ instability, and coherent structures. Furthermore, four test cases with different equivalence ratios are studied to evaluate the effects of equivalence ratio on the flame and flow characteristics. A finite volume scheme on a Cartesian grid with a dynamic one equation eddy viscosity subgrid model is used for large eddy simulations. The obtained results show that the combustion heat release and increase in equivalence ratio toward the stoichiometric value decrease the local swirl number of the flow field, while increasing the flow spreading at the burner outlet. Results show that the flame becomes W shaped as the equivalence ratio increases. Moreover, the combination of the swirling motion and combustion heat release temporally imposes a vortex breakdown in the post-flame region, which leads to occurrence of a transient recirculation zone. The temporal recirculation zone disappears downstream of the burner outlet due to merging of the inner shear layer from all sides at the centerline. Also, various analyses of shear layers’ wavy and vortical structures show that combustion heat release has the effect of decreasing the instability amplitude and vortex shedding frequency.     

Influence of flow field scaling on flashback of swirl flames

In this paper, the effect of geometrical scaling on the onset of flashback into a cylindrical premixing zone of a swirl flame is investigated. We discriminate two types of flashback. In the first type of flashback the flame propagates upstream inside an already present axial recirculation zone. This flashback is caused by turbulent burning along the vortex axis (TBVA¹) and is controlled by flame extinction inside the recirculation zone. The second type of flashback is caused by combustion induced vortex breakdown (CIVB²). This type of flashback is characterised by the aerodynamic influence of the combustion heat release that leads to propagation of the axial recirculation zone and the flame in upstream direction.To study the effects of geometrical scaling on the flow fields and the two types of flashback, the operation of two geometrically scaled burners are compared at equal Reynolds number. By this method it is possible to observe the flashback phenomena in similar swirl flow fields but with different turbulent scales affecting the combustion process. To check flow field similarity and to indentify the flashback type, the non-reacting and reacting flow fields have been examined by planar particle imaging velocimetry and simultaneous recording of the flame luminescence.It is shown that geometrical scaling of the burner shifts the equivalence ratio at which flashback occurs and that this shift is different for the two types of flashback. Consistency and inconsistency with known scaling and stability criterions is discussed. Analysing the fluid dynamics and turbulent combustion gives a first explanation of why CIVB and TBVA are affected differently by geometrical scaling at constant Reynolds number which is in good agreement with the experimental observations.     

模型燃烧室内瞬态紊流的大涡模拟

本文采用三种不同亚网格尺度模型对带有V型稳定器的模型燃烧室二维瞬态紊流流动进行了大涡模拟。并在交错网格系下用SIMPLE算法和混合差分格式求解离散方程。数值研究拟不同型式入口速度分布和不同亚网格尺度模型下模型燃烧室二维瞬态紊流流场。计算结果表明不同入口速度分布和不同亚网格尺度模型对瞬态流场和出口速度分布有一定的影响。本文通过数值模拟,揭示了V型稳定器后旋涡的产生和脱落过程。通过计算结果及实验数据的比较可知,本文采用的亚网格尺度模型可以用来模拟模型燃烧室紊流流场及稳定器后面回流区的流动情况。     

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.     

Flow in a simple swirl chamber with and without controlled inlet forcing

Results are presented from a swirl chamber with and without controlled inlet forcing. The controlled inlet forcing is induced using arrays of vortex generators placed along one wall of the swirl chamber inlet duct. Flow visualization results are given, along with surveys of circumferential mean velocity, static pressure, and total pressure, at Reynolds numbers (based on inlet duct characteristics) as high as 8000. The controlled inlet forcing provides means to alter and control: (i) the spacing and number of Görtler vortices across the span of the swirl chamber, (ii) the amount of vortex development at a particular Reynolds number and circumferential location, (iii) the circumferential location and Reynolds number of initial Görtler vortex development, and (iv) the circumferential location and Reynolds number of Görtler vortex breakup into more chaotic flow.     

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.     

Large-Eddy Simulation of Lean Premixed Turbulent Flames of Three Different Combustion Configurations using a Novel Reaction Closure

This large eddy simulation (LES) study is applied to three different premixed turbulent flames under lean conditions at atmospheric pressure. The hierarchy of complexity of these flames in ascending order are a simple Bunsen-like burner, a sudden-expansion dump combustor, and a typical swirl-stabilized gas turbine burner–combustor. The purpose of this paper is to examine numerically whether the chosen combination of the Smagorinsky turbulence model for sgs fluxes and a novel turbulent premixed reaction closure is applicable over all the three combustion configurations with varied degree of flow and turbulence. A quality assessment method for the LES calculations is applied. The cold flow data obtained with the Smagorinsky closure on the dump combustor are in close proximity with the experiments. It moderately predicts the vortex breakdown and bubble shape, which control the flame position on the double-cone burner. Here, the jet break-up at the root of the burner is premature and differs with the experiments by as much as half the burner exit diameter, attributing the discrepancy to poor grid resolution. With the first two combustion configurations, the applied subgrid reaction model is in good correspondence with the experiments. For the third case, a complex swirl-stabilized burner–combustor configuration, although the flow field inside the burner is only modestly numerically explored, the level of flame stabilization at the junction of the burner–combustor has been rather well captured. Furthermore, the critical flame drift from the combustor into the burner was possible to capture in the LES context (which was not possible with the RANS plus k–ɛ model), however, requiring tuning of a prefactor in the reaction closure.     

Characterization of particle-laden, confined swirling flows by phase-doppler anemometry and numerical calculation

The particle dispersion characteristics in a confined swirling flow with a swirl number of approx. 0.5 were studied in detail by performing measurements using phase-Doppler anemometry (PDA) and numerical predictions. A mixture of gas and particles was injected without swirl into the test section, while the swirling airstream was provided through a co-flowing annular inlet. Two cases with different primary jet exit velocities were considered. For these flow conditions, a closed central recirculation bubble was established just downstream of the inlet.

The PDA measurements allowed the correlation between particle size and velocity to be obtained and also the spatial change in the particle size distribution throughout the flow field. For these results, the behaviour of different size classes in the entire particle size spectrum, ranging from about 15 to 80 μm, could be studied, and the response of the particles to the mean flow and the gas turbulence could be characterized. Due to the response characteristics of particles with different diameters to the mean flow and the flow turbulence, a considerable separation of the particles was observed which resulted in a streamwise increase in the particle mean number diameter in the core region of the central recirculation bubble. For the lower particle inlet velocity (i.e. low primary jet exit velocity), this effect is more pronounced, since here the particles have more time to respond to the flow reversal and the swirl velocity component. This also gave a higher mass of recirculating particle material.

The numerical predictions of the gas flow were performed by solving the time-averaged Navier-Stokes equations in connection with the well known kε turbulence model. Although this turbulence model is based on the assumption of isotropic turbulence, the agreement of the calculated mean velocity profiles compared to the measured gas velocities is very good. The gas-phase turbulent kinetic energy, however, is considerably underpredicted in the initial mixing region. The particle dispersion characteristics were calculated by using the Lagrangian approach, where the influence of the particulate phase on the gas flow could be neglected, since only very low mass loadings were considered. The calculated results for the particle mean velocity and the mass flux are also in good agreement with the experiments. Furthermore, the change in the particle mean diameter throughout the flow field was predicted approximately, which shows that the applied simple stochastic dispersion model also gives good results for such very complex flows. The variation of the gas and particle velocity in the primary inlet had a considerable impact on the particle dispersion behaviour in the swirling flow and the particle residence time in the central recirculation bubble, which could be determined from the numerical calculations. For the lower particle inlet velocity, the maximum particle size-dependence residence time within the recirculation region was considerably shifted towards larger particles.     

LES of Chemical and Acoustic Forcing of a Premixed Dump Combustor

This paper describes the first steps in the development of a large eddy simulation (LES) code able to compute combustion instabilities in gas turbines. This code was used to compute the forcing of an experimentally investigated premixed dump combustor. It is shown that the main effect of acoustic waves entering the combustion chamber is to create large vortices and unsteady heat release when these vortices burn. Another effect of waves entering the combustor is to modulate the fuel and air flow rates produced by the feeding lines. In this case the equivalence ratio of the mixture entering the combustor may also vary. This was investigated in a “chemical effect” simulation where the inlet equivalence ratio fluctuates but the total flow rate remains constant. For perturbations from stoichiometric burning, this mechanism was shown to induce less destabilizing effects than the purely aerodynamical mechanism due to vortex formation and combustion. It is shown that the LES methodology developed is able to reproduce the experimentally observed phase shift between acoustic excitation and total reaction rate in the chamber. This revised version was published online in July 2006 with corrections to the Cover Date.     

双燃式超燃发动机冷态内流场的数值研究

研究了双燃式一体化通道(包含进气道、双燃式燃烧室和尾喷管)的冷态内流场特性.首次在激波风洞中对内流场进行纹影照相,用TVD格式求解三维全N-S方程对喷管和一体化通道进行分区数值模拟,并考察了几何参数对内流场的影响.结果表明对典型工况(h     

Mass and momentum transport in the near field of swirling turbulent jets. Effect of swirl rate

Local transport of the flow momentum and scalar admixture in the near-field of turbulent swirling jets (Re = 5,000) has been investigated by using a combination of the particle image velocimetry and planar laser-induced fluorescence methods. Advection and turbulent and molecular diffusions are evaluated based on the measured distributions of the mean velocity and concentration and the Reynolds stresses and fluxes. As has been quantified from the data, the flow swirl intensifies the entrainment of the surrounding fluid and promotes mass and momentum exchange in the outer mixing layer. A superimposed swirl results in the appearance of a wake/recirculation region at the jet axis and, consequently, the formation of an inner shear layer. In contrast to the scalar admixture, the momentum exchange in the inner shear layer is found to be strongly intensified by the swirl. For the jet with the highest considered swirl rate, a substantial portion of the surrounding fluid is found to enter the unsteady central recirculation zone, where it mixes with the jet that is issued from the nozzle. The contribution of the coherent velocity fluctuations, which are induced by large-scale vortex structures, to the turbulent transport has been evaluated based on triple decomposition, which was based on proper orthogonal decomposition analysis of the velocity data sets. For the considered domain of the jet with the highest swirl rate and vortex breakdown, the contributions of detected helical vortex structures, inducing pressing vortex core, to the radial fluxes of the flow momentum and the scalar admixture are found to locally exceed 65% and 80%, respectively.     

Inviscid evolution of incompressible swirling flows in pipes: The dependence of the flow structure upon the inlet velocity field

This paper analyses the influence of the inlet swirl on the structure of incompressible inviscid flows in pipes. To that end, the inviscid evolution along a pipe of varying radius with a central body situated inside the pipe is studied for three different inlet swirling flows by solving the Bragg–Hawthorne equation both asymptotically and numerically. The downstream structure of the flow changes abruptly above certain threshold values of the swirl parameter (L). In particular, there exist a value L_r above which a near-wall region of flow reversal is formed downstream, and a critical value L_f above which the axial vortex flow breaks down. It is shown that the dependence upon the pipe geometry of these critical values of the swirl parameter varies strongly with the inlet azimuthal velocity profile considered. An excellent agreement between asymptotic and numerical results is found.     

Experiments on the formation of a recirculation zone in swirling coaxial jets

 The near flow field of coaxial air jets, with swirl in the outer one, was studied using flow visualization and hot-wire anemometry. The flow is sensitive to both the swirl number and the mass flow ratio between the outer and inner jets. A necessary condition for the formation of an internal recirculation zone (IRZ) is that the swirl number must exceed a minimum value which depends on the mass flow ratio. Spectral analysis of the velocity fluctuations indicates that the formation of an IRZ in the present flow does not appear to be related to the growth of convective flow instabilities. Analysis of the flow visualization and X-wire data indicates that the vorticity dynamics model for vortex breakdown proposed by Brown and Lopez [J Fluid Mech (1990) 222: 553] provides a plausible mechanism for the formation of an IRZ in this flow. Received: 14 June 1999/Accepted: 7 December 1999     

Dynamics of annular gas–liquid two-phase swirling jets

The dynamics of annular gas–liquid two-phase swirling jets have been examined by means of direct numerical simulation and proper orthogonal decomposition. An Eulerian approach with mixed-fluid treatment, combined with an adapted volume of fluid and a continuum surface force model, was used to describe the two-phase flow system. The unsteady, compressible, three-dimensional Navier–Stokes equations have been solved by using highly accurate numerical methods. Two computational cases have been performed to examine the effects of liquid-to-gas density ratio on the flow development. It was found that the higher density ratio case is more vortical with larger spatial distribution of the liquid, in agreement with linear theories. Proper orthogonal decomposition analysis revealed that more modes are of importance at the higher density ratio, indicating a more unstable flow field. In the lower density ratio case, both a central and a geometrical recirculation zone are captured while only one central recirculation zone is evident at the higher density ratio. The results also indicate the formation of a precessing vortex core at the high density ratio, indicating that the precessing vortex core development is dependent on the liquid-to-gas density ratio of the two-phase flow, apart from the swirl number alone.     

Interaction of heat release and vortex breakdown during flame flashback driven by combustion induced vortex breakdown

The interaction of heat release by chemical reaction and the flow dominates flame transition in swirling flows caused by combustion induced vortex breakdown (CIVB). The simultaneous application of 1 kHz high-speed particle imaging velocimetry (PIV) for the analysis of the flow field and OH planar laser-induced fluorescence for the detection of the flame front is particularly useful for the improvement of the understanding of the observed fast CIVB driven flame propagation. For the first time, the combination of both techniques was successfully applied to confined swirling flows. In the study, the flow field characteristics of an aerodynamically stabilized burner system with CIVB are analyzed in great depth. The influence of geometric parameters of the swirl generator was investigated and conclusions concerning the proper burner design of vortex breakdown premix burners are drawn from the experimental results. In particular, the effect of the vortex core with respect to the stability of the swirl stabilized burner is analyzed. The contribution of combustion to vortex breakdown is shown comparing isothermal and reacting flows. The presented data reveals that at the onset of CIVB driven flame transition, the azimuthal vorticity leads to the formation of a closed recirculation bubble at the tip of the internal recirculation zone. Once this bubble propagates upstream, the flame is able to follow and propagate relative to the bulk flow velocity with a velocity far beyond the turbulent flame speed. The interaction of reaction and flow was observed for different volumetric heat releases. The experiments confirm the CIVB theory of the authors, which was initially developed on the basis of a CFD study alone. Both the volume expansion and the baroclinic torque have an effect on whether fast flame propagation occurs. Whereas the volume expansion caused by the heat release stabilizes the flow field and the reaction, the baroclinic torque stimulates flame transition. For upstream propagation the flame tip has to have a position downstream of the stagnation point of the bubble. Else, the required transition inducing force is insufficient and the flame remains stable. In case the flame reaches positions too close or even upstream of the stagnation point, the fast propagation is interrupted or even prohibited. The key finding that the relative position of flame and stagnation bubble governs CIVB is discussed on the basis of high-speed LIF/PIV data as well as chemiluminescence. Since essentially the same behavior has been observed before in tests of a totally different swirler design and flow field, the conclusion can be made that the root cause for CIVB independent of the special geometry has been found.