Large Eddy Simulations of gaseous flames in gas turbine combustion chambers

Recent developments in numerical schemes, turbulent combustion models and the regular increase of computing power allow Large Eddy Simulation (LES) to be applied to real industrial burners. In this paper, two types of LES in complex geometry combustors and of specific interest for aeronautical gas turbine burners are reviewed: (1) laboratory-scale combustors, without compressor or turbine, in which advanced measurements are possible and (2) combustion chambers of existing engines operated in realistic operating conditions. Laboratory-scale burners are designed to assess modeling and funda- mental flow aspects in controlled configurations. They are necessary to gauge LES strategies and identify potential limitations. In specific circumstances, they even offer near model-free or DNS-like LES computations. LES in real engines illustrate the potential of the approach in the context of industrial burners but are more difficult to validate due to the limited set of available measurements. Usual approaches for turbulence and combustion sub-grid models including chemistry modeling are first recalled. Limiting cases and range of validity of the models are specifically recalled before a discussion on the numerical breakthrough which have allowed LES to be applied to these complex cases. Specific issues linked to real gas turbine chambers are discussed: multi-perforation, complex acoustic impedances at inlet and outlet, annular chambers.. Examples are provided for mean flow predictions (velocity, temperature and species) as well as unsteady mechanisms (quenching, ignition, combustion instabil- ities). Finally, potential perspectives are proposed to further improve the use of LES for real gas turbine combustor designs

Using LES to Study Reacting Flows and Instabilities in Annular Combustion Chambers

Great prominence is put on the design of aeronautical gas turbines due to increasingly stringent regulations and the need to tackle rising fuel prices. This drive towards innovation has resulted sometimes in new concepts being prone to combustion instabilities. In the particular field of annular combustion chambers, these instabilities often take the form of azimuthal modes. To predict these modes, one must compute the full combustion chamber, which remained out of reach until very recently and the development of massively parallel computers. Since one of the most limiting factors in performing Large Eddy Simulation (LES) of real combustors is estimating the adequate grid, the effects of mesh resolution are investigated by computing full annular LES of a realistic helicopter combustion chamber on three grids, respectively made of 38, 93 and 336 million elements. Results are compared in terms of mean and fluctuating fields. LES captures self-established azimuthal modes. The presence and structure of the modes is discussed. This study therefore highlights the potential of LES for studying combustion instabilities in annular gas turbine combustors

LES evaluation of the effects of equivalence ratio fluctuations on the dynamic flame response in a real gas turbine combustion chamber

Large Eddy Simulations (LES) of a lean swirl-stabilized gas turbine burner are used to analyze mechanisms triggering combustion instabilities. To separately study the effect of velocity and equivalence ratio fluctuations, two LES of the same geometry are performed: one where the burner operates in a “technically” premixed mode (methane is injected by holes in the vanes located in the diagonal passage upstream of the chamber) and the second one where the flow is fully premixed in the diagonal passage. The inlet is acoustically modulated and the mechanisms affecting the dynamic flame response are identified. LES reveals that both cases provide similar averaged (non-)pulsated flame shapes. However, even though the mean flames are only slightly modified, the delays change when mixing is not perfect. LES fields and a simple model for the methane jets trajectories show that mixing in the diagonal passage is not sufficient to damp heterogeneities induced by unsteady fuel flow rate and varying fuel jet trajectories. These mixing fluctuations are phased with velocity oscillations and modify the flame response to forcing. Local fields of delays and interaction indices are obtained, showing that the flame is not compact and is affected by fluctuations of mixing

Acoustic and Large Eddy Simulation studies of azimuthal modes in annular combustion chambers

The objectives of this paper are the description of azimuthal instability modes found in annular combus- tion chambers using two numerical tools: (1) Large Eddy Simulation (LES) methods and (2) acoustic solv- ers. These strong combustion instabilities are difficult to study experimentally and the present study is based on a LES of a full aeronautical combustion chamber. The LES exhibits a self-excited oscillation at the frequency of the first azimuthal eigenmode. The mesh independence of the LES is verified before ana- lysing the nature of this mode using various indicators over more than 100 cycles: the mode is mostly a pure standing mode but it transitions from time to time to a turning mode because of turbulent fluctu- ations, confirming experimental observations and theoretical results. The correlation between pressure and heat release fluctuations (Rayleigh criterion) is not verified locally but it is satisfied when pressure and heat release are averaged over sectors. LES is also used to check modes predicted by an acoustic Helmholtz solver where the flow is frozen and flames are modelled using a Flame Transfer Function (FTF) as done in most present tools. The results in terms of mode structure compare well confirming that the mode appearing in the LES is the first azimuthal mode of the chamber. Moreover, the acoustic solver provides stability maps suggesting that a reduction of the time delay of the FTF would be enough to sta- bilise the mode. This is confirmed with LES by increasing the flame speed and verifying that this modi- fication leads to a damped mode in a few cycles

Large-Eddy Simulation of combustion instabilities in a variable-length combustor.

This article presents a simulation of a model rocket combustor with continuously variable acoustic properties thanks to a variable-length injector tube. Fully compressible Large-Eddy Simulations are conducted using the AVBP code. An original flame stabilization mechanism is uncovered where the recirculation of hot gases in the corner recirculation zone creates a triple flame structure. An unstable operating point is then chosen to investigate the mech- anism of the instability. The simulations are compared to experimental results in terms of frequency and mode structure. Two-dimensional axi-symmetric computations are com- pared to full 3D simulations in order to assess the validity of the axi-symmetry assumption for the prediction of mean and unsteady features of this flow. Despite the inaccuracies in- herent to the 2D description of a turbulent flow, for this configuration and the particular operating point investigated, the axi-symmetric simulation qualitatively reproduces some features of the instability

Triple flame structure and diffusion flame stabilization

The stabilization of diffusion ñames is studied using asymptotic techniques and numerical tools. The configuration studied corresponda to parallel streams of cold oxidizer and fuel initially separated by a splitter píate. It is shown that stabilization of a diffusion flame may only occur in this situation by two processes. First, the flame may be stabilized behind the flame holder in the wake of the splitter píate. For this case, numerical simulations confirm scalings previously predicted by asymptotic analysis. Second, the flame may be lifted. In this case a triple flame is found at longer distanees downstream of the flame holder. The structure and propagation speed of this flame are studied by using an actively controlled numerical technique in which the triple flame is tracked in its own reference frame. It is then possible to investigate the triple flame structure and velocity. It is shown, as suggested from asymptotic analysis, that heat reléase may induce displacement speeds of the triple flame larger than the laminar flame speed corresponding to the stoichiometric conditions prevailing in the mixture approaching the triple flame. In addition to studying the characteristics of triple flames in a uniform flow, their re-sistance to turbulence is investigated by subjecting triple flames to different vortical configurations

LES-CMC simulations of different auto-ignition regimes of hydrogen in a hot turbulent air Co-flow

Large-Eddy Simulation (LES) results in combination with first-order Conditional Moment Closure (CMC) are presented for a hydrogen jet, diluted with nitrogen, issued into a turbulent co-flowing hot air stream. The fuel mixes with the co-flow air, ignites and forms a lifted-like flame. Global trends in the experimental observations are in general well reproduced: the auto-ignition length decreases with increase in co-flow temperature and increases with increase in co-flow velocity. In the experiments, the co-flow temperature was varied, so that different auto-ignition regimes, including low Damkohler number situations, were obtained (no ignition, random spots, flashback and lifted flame). All regimes are recovered in the simulations. Auto-ignition is found to be the stabilizing mechanism. The impact of different detailed chemistry mechanisms on the auto-ignition predictions is discussed. With increasing air temperature, the differences between the mechanisms considered diminish. The evolution of temperature, H2O, H, HO2 and OH from inert to burning conditions is discussed in mixture fraction space

Prediction of combustion instability limit cycle oscillations by combining flame describing function simulations with a thermoacoustic network model

Accurate prediction of limit cycle oscillations resulting from combustion instability has been a long-standing challenge. The present work uses a coupled approach to predict the limit cycle characteristics of a combustor, developed at Cambridge University, for which experimental data are available (Balachandran, Ph.D. thesis, 2005). The combustor flame is bluff-body stabilised, turbulent and partially-premixed. The coupled approach combines Large Eddy Simulation (LES) in order to characterise the weakly non-linear response of the flame to acoustic perturbations (the Flame Describing Function (FDF)), with a low order thermoacoustic network model for capturing the acoustic wave behaviour. The LES utilises the open source Computational Fluid Dynamics (CFD) toolbox, OpenFOAM, with a low Mach number approximation for the flow-field and combustion modelled using the PaSR (Partially Stirred Reactor) model with a global one-step chemical reaction mechanism for ethylene/air. LES has not previously been applied to this partially-premixed flame, to our knowledge. Code validation against experimental data for unreacting and partially-premixed reacting flows without and with inlet velocity perturbations confirmed that both the qualitative flame dynamics and the quantitative response of the heat release rate were captured with very reasonable accuracy. The LES was then used to obtain the full FDF at conditions corresponding to combustion instability, using harmonic velocity forcing across six frequencies and four forcing amplitudes. The low order thermoacoustic network modelling tool used was the open source OSCILOS (http://www.oscilos.com). Validation of its use for limit cycle prediction was performed for a well-documented experimental configuration, for which both experimental FDF data and limit cycle data were available. The FDF data from the LES for the present test case was then imported into the OSCILOS geometry network and limit cycle oscillations of frequency 342 Hz and normalised velocity amplitude of 0.26 were predicted. These were in good agreement with the experimental values of 348 Hz and 0.21 respectively. This work thus confirms that a coupled numerical prediction of limit cycle behaviour is possible using an entirely open source numerical framework