Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Large eddy simulation of soot evolution in turbulent nonpremixed bluff body flames

Large Eddy Simulation (LES) is utilized to investigate soot evolution in a series of turbulent nonpremixed bluff body flames featuring different bluff body diameters. The modeling framework relies on recent development in the soot subfilter Probability Density Function (PDF) model that can correctly account for the distribution of soot with respect to mixture fraction, correcting errors in previous soot subfilter PDF models that significantly overpredict soot oxidation. With the previous soot subfilter PDF model, no soot was observed outside of the recirculation zone in past studies on similar bluff body flames. Results obtained with the current LES modeling approach compare favorably with the experimental measurements of the flow field and the soot volume fraction. Notably, the current LES modeling approach correctly predicts large soot volume fractions in the recirculation zone, a decrease in the soot volume fraction through the high-strain neck region, and then an increase again in the downstream jet-like region. Consistent with the experimental measurements, the larger bluff body diameter, with its larger recirculation zone with longer residence times, generates more soot in the recirculation zone and also more soot in the high-strain neck region. Analysis of the soot volume fraction source terms lead to mechanistic understanding of soot evolution in the entirety of the bluff body flames. Most of the soot generated in the recirculation zone is oxidized but some escapes unoxidized and is passively transported through the neck region. Virtually no new soot forms in the downstream jet-like region, and the increase in the soot volume fraction in the jet-like region is due to acetylene-based surface growth of the soot transported through the neck region. This mechanism could not be predicted with the previous soot subfilter PDF model, with the recent soot subfilter PDF model being critical in the understanding of this basic mechanism.     

Subgrid combustion modeling of 3-D premixed flames in the thin-reaction-zone regime

Large eddy simulation (LES) is used to investigate three-dimensional (3D) lean premixed turbulent methane–air flames in the thin-reaction-zone regime. In this regime, the Kolmogorov scale is smaller than the preheat zone thickness, but larger than the reaction zone thickness. Past numerical studies of similar flames were primarily direct numerical simulation either in two-dimensions or using the artificially thickened flame approach in 3D. For an LES the effect of small (unresolved) scales on the scalar field must be, modeled accurately to capture the correct flame structure. A subgrid combustion model based on the linear-eddy-mixing (LEM) model is used within an LES framework (called LEM–LES hereafter) to capture the 3D flame-structure of the highly stretched premixed flames. A finite-rate, one-step methane–air chemistry with a non-unity Lewis number formulation is used in this study. The simulated flame structure resembles flames experimentally studied in the thin-reaction-zone regime. Even though the preheat zone is broadened by the penetration of small eddies, the chemical reaction zone remains thin and localized. This feature is captured properly in the current LEM–LES approach. The flame structure and other statistics such as the flame area evolution, curvature, and strain-rate statistics computed using the LEM–LES are also in good agreement with the past DNS studies.     

Large eddy simulation of multi-regime combustion with a two-progress variable approach for carbon monoxide

Simulations of two cases in a novel multi-regime burner configuration are undertaken using a presumed joint probability density function (PDF) approach with tabulated chemistry. The flame conditions are varied by changing the central jet equivalence ratio, which produces different multi-regime combustion modes in the non-premixed inner flame. An outer premixed flame and recirculation zone behind a bluff body are present to supply heat and combustion products to stabilise the inner flame. A two-progress variable approach is tested to improve predictions of carbon monoxide (CO) in the post-flame regions, where CO oxidation occurs. The large eddy simulation set-up and sub-grid combustion model are assessed through comparisons with time-averaged measurements for radial profiles at different streamwise locations. The jet break-up length, the shear layers and the mixture fraction distribution are well captured in both cases. The temperature distribution is well captured for the inner flame in each case but the temperature and mixture fraction are over predicted in the downstream regions of the outer premixed flame, which is due to increased dilatation that suppresses air entrainment. Improved predictions of the CO mass fraction are obtained for the outer premixed flames with the two-progress variable approach. Over predictions are seen in the upstream regions of the inner flame when the CO mass fraction is obtained from a look-up table, suggesting that the CO mass fraction should be transported to include the convection/diffusion balance in regions where there is no flame. Furthermore, transporting the CO mass fraction with a one-progress variable approach produces over predictions in the burnt regions, suggesting a two-progress variable model is needed to capture the consumption region of CO. The multi-regime combustion characteristics are observed to be stronger in flame MRB26b, where non-premixed and rich premixed combustion is present. For flame MRB18b, the non-premixed contribution is smaller and weak stratified combustion is observed.     

Impact of swirl and bluff-body on the transfer function of premixed flames

The frequency response of three lean methane/air flames submitted to flowrate perturbations is analyzed for flames featuring the same equivalence ratio and thermal power, but a different stabilization mechanism. The first flame is stabilized by a central bluff body without swirl, the second one by the same bluff body with the addition of swirl and the last one only by swirl without central insert. In the two last cases, the swirl level is roughly the same. These three flames feature different shapes and heat release distributions, but their Flame Transfer Function (FTF) feature about the same phase lag at low frequencies. The gain of the FTF also shows the same behavior for the flame stabilized by the central insert without swirl and the one fully aerodynamically stabilized by swirl. Shedding of vortical structures from the injector nozzle that grow and rollup the flame tip controls the FTF of these flames. The flame stabilized by the swirler-plus-bluff-body system features a peculiar response with a large drop of the FTF gain around a frequency at which large swirl number oscillations are observed. Velocity measurements in cold flow conditions reveal a strong reduction of the size of the vortical structures shed from the injector lip at this forcing condition. The flame stabilized aerodynamically only by swirl and the one stabilized by the bluff body without swirl do not exhibit any FTF gain drop at low frequencies. In the former case, large swirl number oscillations are still identified, but large vortical structures shed from the nozzle also persist at the same forcing frequency in the cold flow response. These different flame responses are found to be intimately related to the dynamics of the internal recirculation region, which response strongly differs depending upon the injector used to stabilize the flame.     

Large eddy simulations of turbulent flames using the filtered density function model

Large eddy simulations (LES) of the Sandia/Sydney swirl burners (SM1 and SMA1) and the Sandia/Darmstadt piloted jet diffusion flame (Flame D) are performed. These flames are part of the database of turbulent reacting flows widely considered as benchmark test cases for validating turbulent-combustion models. In the simulations presented in this paper, the subgrid scale (SGS) closure model adopted for turbulence-chemistry interactions is based on the transport filtered density function (FDF) model. In the FDF model, the transport equation for the joint probability density function (PDF) of scalars is solved. The main advantage of this model is that the filtered reaction rates can be exactly computed. However, the density field, computed directly from the FDF solver and needed in the hydrodynamic equations, is noisy and causes numerical instability. Two numerical approaches that yield a smooth density field are examined. The two methods are based on transport equations for specific sensible enthalpy (h_s) and RT, where R is the gas constant and T is the temperature. Consistency of the two methods is assessed in a bluff-body configuration using Reynolds averaged Navier-Stokes (RANS) methodology in conjunction with the transported PDF method. It is observed that the h_s method is superior to the RT method. Both methods are used in LES of the SM1 burner. In the near-field region, the h_s method produces better predictions of temperature. However, in the far-field region, both methods show deviation from data. Simulations of the SMA1 burner and Flame D are also presented using the h_s method. Some deficiencies are seen in the predictions of the SMA1 burner that may be related to the simple chemical kinetics model and mixing model used in the simulations. Simulations of Flame D show good agreement with data. These results indicate that, while further improvements to the methodology are needed, the LES/FDF method has the potential to accurately predict complex turbulent flames.     

Effect of Lewis number on premixed laminar lean-limit flames stabilized on a bluff body

In this paper, we present a study on the effect of Lewis number, Le, on the stabilization and blow-off of laminar lean limit premixed flames stabilized on a cylindrical bluff body. Numerical simulations and experiments are conducted for propane, methane and two blends of hydrogen with methane as fuel gases, containing 20% and 40% of hydrogen by volume, respectively. It is found that the Le?>?1 flame blows-off via convection from the base of the flame (without formation of a neck) when the conditions for flame anchoring are not fulfilled. Le?≤?1 flames exhibit a necking phenomenon just before lean blow-off. This necking of the flame front is a result of the local reduction in mass burning rates causing flame merging and quenching of the thin flame tube formed. The structure of these flames at the necking location is found to be similar to tubular flames. It is found that extinction stretch rates for tubular flames closely match values at the neck location of bluff-body flames of corresponding mixtures, suggesting that excessive flame stretch is directly responsible for blow-off of the studied Le?≤?1 flames. After quenching of the neck, the upstream part forms a steady and stable residual flame in the wake of the bluff body while the downstream part is convected away.     

Partially premixed flamelet in LES of acetone spray flames

An effective partially premixed flamelet model for large eddy simulation (LES) of turbulent spray combustion is formulated. Different flame regimes are identified with a flame index defined by budget terms in a 2-D multi-phase flamelet formulation, and the application in LES of partially pre-vaporized spray flames shows a favorable agreement with experiments. Simulations demonstrate that, compared to the conventional single-regime flamelets, the present partially premixed flamelet formulation shows its ability in capturing the subgrid regime transitions, yielding a well prediction of peak gas temperature and the downstream flame spreading. A propagating premixed flame front is found coupled with a trailing diffusion burning through the spray evaporation, and the spray effect on regime discrimination is manifested with transport budget analysis. A two-phase regime indicator is then proposed, by which the evaporation-dictated regime is properly described. Its intended use will rely on both gas and spray flamelet structures.     

Effects of CO2 dilution on partially premixed CH4-air flames in swirl and bluff body stabilized combustor

The effect of CO₂ dilution on the flame characteristics and pollutant emission of a partially premixed CH₄-air flame in a confined bluff body and swirl influenced flowfield is investigated using optical and laser diagnostic methods. The non-premixed burner produced a converging-diverging flowfield at the burner exit and a lifted flame is produced at all test cases, with an upstream movement of the flame with decreasing global equivalence ratios (?_g). Based on variations in ?_g, two flame stabilization modes – bluff body influenced and swirl stabilized – with a transition mode in-between is observed for the cases with (flame FB) and without dilution (flame FM). The characteristics of the heat release zone are influenced by dilution, with the FB flames being longer and also less intense when compared to FM flames. Pollutant measurement at 30 mm downstream from the combustor exit highlighted the ultra-low NO_x capability of the IIST-GS2 burner. CO₂ dilution leads to a reduction in NO_x emission due to both thermal and chemical effects. For ?_g ≥ 0.7 extreme low levels of CO and unburned hydrocarbons (UHC) are observed for both cases. For ?_g ≤ 0.6 the dramatic increase of both CO and UHC maybe due to the lower flame temperatures and shorter flame zone residence times, respectively.     

Flame stabilization regimes for premixed flames anchored behind cylindrical flame holders

The objective of this study is to construct a regime diagram for laminar flames stabilized behind flame holders with respect to the presence of a recirculation zone (RZ), trend of heat loss to the burner, and flow strain and flame curvature effects. This is achieved by varying the radius of the cylindrical flame holder and the mixture velocity between the flashback limit and the blow-off limit at a fixed equivalence ratio. It is found that for all flame holders, a RZ vortex is not present near the flashback limit. At flashback, flow strain is almost zero and the flame curvature is found to be the main contributor to flame stretch. With increasing mixture velocity, the heat loss to the flame holder decreases for smaller radii and a RZ does not appear till blow-off occurs. For flame holders with radii greater than twice the flame thickness, the heat loss to the flame holder first decreases with increasing mixture velocity without a RZ. A further increase in the mixture velocity does not result in blow-off but instead, a RZ appears behind the flame holder reversing the heat loss trend. In this scenario, flow strain is found to increase significantly and becomes the major contributor to flame stretch, although curvature effects are still present. With the RZ present, the blow-off limits are significantly extended and the stabilization mechanism is altered. The RZ vortex shields the flame base from intense pre-heating resulting from the increase in heat loss to the flame-holder while it provides support to the flame leading edge by recirculation of hot products. The results obtained from this study are used to construct a regime diagram, which offers a broader view of the whole flame stabilization process and its mechanisms.     

A posteriori testing of the flame surface density transport equation for LES

Flame Surface Density (FSD) models for Large Eddy Simulation (LES) are implemented and tested for a canonical configuration and a practical bluff body stabilised burner, comparing common algebraic closures with a transport equation closure in the context of turbulent premixed combustion. The transported method is expected to yield advantages over algebraic closures, as the equilibrium of subgrid production and destruction of FSD is no longer enforced and resolved processes of strain, propagation and curvature are explicitly accounted for. These advantages might have the potential to improve the ability to capture large-scale unsteady flame propagation in situations with combustion instabilities or situations where the flame encounters progressive wrinkling with time. The initial study of a propagating turbulent flame in wind-tunnel turbulence shows that the Algebraic Flame Surface Density (FSDA) method can predict an excessively wrinkled flame under fine grid conditions, potentially increasing the consumption rate of reactants to artificially higher levels. In contrast, the Flame Surface Density Transport (FSDT) closure predicts a smooth flame front and avoids the formation of artificial flame cusps when the grid is refined. Five FSDA models and the FSDT approach are then applied to the LES of the Volvo Rig. The predicted mean velocities are found to be relatively insensitive to the use of the FSDT and FSDA approaches, whereas temperature predictions exhibit appreciable differences for different formulations. The FSDT approach yields very similar temperature predictions to two of the tested FSDA models, quantitatively capturing the mean temperature. Grid refinement is found to improve the FSDT predictions of the mean flame spread. Overall, the paper demonstrates that the apparently complicated FSD transport equation approach can be implemented and applied to realistic, strongly wrinkled flames with good success, and opens up the field for further work to improve the models and the overall FSDT approach.     

Pressure gradient tailoring effects on vorticity dynamics in the near-wake of bluff-body premixed flames

We investigate the role of mean streamwise pressure gradients in the development of a bluff-body-stabilized premixed flame in the near-wake of the bluff body. To this end, a triangular prism flame holder is situated in three different channel geometries: a nominal case with straight walls, a nozzle with a stronger mean pressure gradient, and a diffuser with a comparatively weaker mean pressure gradient. All geometries are implemented using embedded boundaries, and adaptive mesh refinement is used to locally resolve all relevant thermal (i.e., flame) and fluid-mechanical (i.e., vorticity) scales. A premixed propane flame, modeled using a 66-step skeletal mechanism, interacts with vorticity in the boundary layer of the triangular bluff body in the presence of each mean pressure gradient. Analysis of flame-related enstrophy budget terms reveals key differences in the behavior of baroclinic torque between cases, the specifics of which are tied to larger variations in the mean flow structure, recirculation zone structure, and confinement effects. Our results show that the baroclinic torque changes significantly among the configurations, with the nozzle exhibiting the largest baroclinic torque production. However, these differences are shown to be only a secondary consequence of the background pressure gradient, with the primary consequence being the change in the recirculation zone length resulting from the different channel configurations. These results are relevant for flame stabilization with bluff bodies, where clear understanding of the sensitivities to global mean pressure gradient is important to engineering design.     

LES of the Sydney swirl flame series: A study of vortex breakdown in isothermal and reacting flows

This work examines the flow and mixing in selected non-premixed cases of the Sydney Swirl Flame series by Large Eddy Simulation. A mixture fraction approach with a steady flamelet model, based on a detailed chemical mechanism, is applied to determine the chemical state in the flame. The isothermal case N29S054 is simulated to provide insights into the flow field, the resolution requirements for the simulation, and to allow for various measures of validation and verification. For the reactive case, the high-speed hydrogen/methane flame SMH1 is chosen for its similarity to the non-swirling Sydney flames and its good stability. In experiments carried out previously at Sydney University, vortex breakdown has been observed and in the isothermal case, this is clearly predicted by the LES. However, no vortex breakdown is observed in the simulations of flame SMH1, which necessitates further studies on this and similar flames to investigate this phenomenon. Studies of the low-velocity flames SM1 and SM2 show that reactive vortex breakdown can be predicted successfully. This difficulty in the prediction of vortex breakdown is another indication that the Sydney Swirl Flame series, especially at high velocities of the central jet, is an interesting and challenging test-case for the development of combustion LES.     

Ultra-lean hydrogen-enriched oscillating flames behind a heat conducting bluff-body: Anomalous and normal blow–off

This paper presents a numerical study of ultra-lean hydrogen-methane flames stabilized behind a rectangular, highly conducting metallic bluff body acting as a flame holder. Using high fidelity numerical simulations, we show that lean inverted steady flames exist below normal flammability limits. They have distinct stabilization mechanism from pure methane flames. These flames are blown-off for sufficiently small velocities, a phenomenon called anomalous blow-off. At even leaner conditions oscillating ultra–lean hydrogen-methane flames can be established. These oscillating flames exist within a rather small range of equivalence ratios and inflow velocities, and move to mean locations closer to the flame holder as the reactant flow is increased. We show that the oscillations are associated with the shedding of flame balls from the downstream end of a “residual flame” that remains attached. Unlike their steady counterparts, the oscillating flames exhibit blow-off at both low velocities (anomalous blow-off) and at sufficiently high inflow velocities (normal blow-off). We show that normal blow-off is linked to heat losses to the flame holder.     

A physical relationship between consumption and displacement speed for premixed flames with finite thickness

The characterization of premixed flames by a flame speed has been a subject that has occupied much interest in the literature in many systematic studies on combustion phenomena. Consumption and displacement speeds are two such flame speeds that are understood to describe the flame dynamics under the effect of flame curvature, flow non-uniformities, Lewis number and turbulence effects along with heat transfer with flame holders and cold walls. As such, much work has been done in the past where either one of these two speeds has been employed along with a linear sensitivity coefficient (Markstein length) for describing different sensitivities to stretch effects. However, despite recent attempts using the asymptotic theory, the relationship between these two quantities has only been clarified in a limited manner for flames of finite thickness. In this study, we use flame stretch theory that takes into account changes of stretch, curvature, heat transfer and Lewis number effects throughout the pre-heat zone and its integral effect on the flame reaction zone. A sound mathematical and physical basis is provided for understanding the two speeds that is valid for weak as well as strong stretch effects. Understanding from theory is further demonstrated by analysing several example 1D stretched flames along with a 2D bluff body flame near extinction.     

LES-CMC simulations of a turbulent bluff-body flame

The large Eddy simulations (LES)-conditional moment closure (CMC) method with detailed chemistry is applied to a bluff-body stabilized flame. Computations of the velocity and mixture fraction fields show good agreement with the experiments. Temperature and major species are well-predicted throughout the flame with the exception of the flow regions in the outer shear layer close to the nozzle where the pure mixing between hot recirculating products and fresh oxidizer cannot be captured. LES-CMC generally improves on results obtained with RANS-CMC and on LES that uses one representative flamelet to model the dependence of reactive species on mixture fraction. Simulated CO mass fractions are generally in good agreement with the experimental data although a 10% overprediction can be found at downstream positions. NO predictions show a distinct improvement over the flamelet approach, however, simulations overpredict NO mass fractions at all downstream locations due to an overprediction of temperature close to the nozzle. The potential of LES-CMC to predict unsteady finite rate effects is demonstrated by the prediction of endothermic—or “flame cooling”—regions close to the neck of the recirculation zone that favours ethylene production via the methane fuel decomposition channel.     

Application of the Lagrangian CMC method in elliptic bluff body and swirl flames

The Lagrangian CMC method was implemented in the open source programme OpenFOAM and applied to turbulent nonpremixed bluff body and swirl flames. Lagrangian CMC is more efficient than Eulerian CMC with the number of Lagrangian flame groups much less than the number of computational cells for Eulerian CMC equations in general. It is based on the conditional flame structure depending on the residence time of the fuel of fixed Lagrangian identity from the nozzle. According to sensitivity study the injected fuel was divided into ten flame groups according to the injection sequence with the resulting conditional profiles between those by ISR and Eulerian CMC. Minor deviation from Eulerian CMC was attributed to the flame structure that is difficult to be characterised by the residence time only in elliptic recirculating flows of the bluff body and swirl flames. The Eulerian and Lagrangian CMC showed the same trend of deviation from measurements for conditional temperature, H₂O, OH, CO and H₂ mass fractions. The significant deviation of H₂ was due to uncertainty in the reaction chemistry, as observed in the previous works based on other reaction mechanisms for methane and methanol.     

Large eddy simulation of piloting effects on turbulent swirling flames

Pilot flames, created by additional injectors of pure fuel, are often used in turbulent burners to enhance flame stabilization and reduce combustion instabilities. The exact mechanisms through which these additional rich zones modify the flame anchoring location and the combustion dynamics are often difficult to identify, especially when they include unsteady hydrodynamic motion. This study presents Large Eddy Simulations (LES) of the reacting flow within a large-scale gas turbine burner for two different cases of piloting, where either 2 or 6% of the total methane used in the burner is injected through additional pilot flame lines. For each case, LES shows how the pilot fuel injection affects both flame stabilization and flame stability. The 6% case leads to a stable flame and limited hydrodynamic perturbations in the initial flame zone. The 2% case is less stable, with a small-lift-off of the flame and a Precessing Vortex Core (PVC) in the cold stabilization zone. This PVC traps some of the lean cold gases issuing from the pilot passage stream, changes the flame stabilization point and induces instability.     

Effect of swirl on combustion dynamics in a lean-premixed swirl-stabilized combustor

The effect of inlet swirl on the flow development and combustion dynamics in a lean-premixed swirl-stabilized combustor has been numerically investigated using a large-eddy-simulation (LES) technique along with a level-set flamelet library approach. Results indicate that when the inlet swirl number exceeds a critical value, a vortex-breakdown-induced central toroidal recirculation zone is established in the downstream region. As the swirl number increases further, the recirculation zone moves upstream and merges with the wake recirculation zone behind the centerbody. Excessive swirl may cause the central recirculating flow to penetrate into the inlet annulus and lead to the occurrence of flame flashback. A higher swirl number tends to increase the turbulence intensity, and consequently the flame speed. As a result, the flame surface area is reduced. The net heat release, however, remains almost unchanged because of the enhanced flame speed. Transverse acoustic oscillations often prevail under the effects of strong swirling flows, whereas longitudinal modes dominate the wave motions in cases with weak swirl. The ensuing effect on the flow/flame interactions in the chamber is substantial.     

Effect of fuel mixture fraction and velocity perturbations on the flame transfer function of swirl stabilized flames

A novel methodology is developed to decompose the classic Flame Transfer Function (FTF) used in the thermo-acoustic stability analysis of lean premix combustors into contributions of different types. The approach is applied, in the context of Large Eddy Simulation (LES), to partially-premixed and fully-premixed flames, which are stabilized via a central recirculation zone as a result of the vortex breakdown phenomenon. The first type of decomposition is into contributions driven by fuel mixture fraction and dynamic velocity fluctuations. Each of these two contributions is further split into the components of turbulent flame speed and flame surface area. The flame surface area component, driven by the pure dynamic velocity fluctuation, which is shown to be a dominant contribution to the overall FTF, is also additionally decomposed over the coherent flow structures using proper orthogonal decomposition. Using a simplified model for the dynamic response of premixed flames, it is shown that the distribution of the FTF, as obtained from LES, is closely related to the characteristics of the velocity field frequency response to the inlet perturbation. Initially, the proposed method is tested and validated with a well characterized laboratory burner geometry. Subsequently, the method is applied to an industrial gas turbine burner.