A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

Three-dimensional numerical simulations of cellular jet diffusion flames

Recent experimental investigations have demonstrated that the appearance of particular cellular states in circular non-premixed jet flames significantly depends on a number of parameters, including the initial mixture strength, reactant Lewis numbers, and proximity to the extinction limit (Damköhler number). For CO₂-diluted H₂/O₂ jet diffusion flames, these studies have shown that a variety of different cellular patterns or states can form. For given fuel and oxidizer compositions, several preferred states were found to co-exist, and the particular state realized was determined by the initial conditions. To elucidate the dynamics of cellular instabilities, circular non-premixed jet flames are modeled with a combination of three-dimensional numerical simulation and linear stability analysis (LSA). In both formulations, chemistry is described by a single-step, finite-rate reaction, and different reactant Lewis numbers and molecular weights are specified. The three-dimensional numerical simulations show that different cellular flames can be obtained close to extinction and that different states co-exist for the same parameter values. Similar to the experiments, the behavior of the cell structures is sensitive to (numerical) noise. During the transient blow-off process, the flame undergoes transitions to structures with different number of cells, while the flame edge close to the nozzle oscillates in the streamwise direction. For conditions similar to the experiments discussed, the LSA results reveal various cellular instabilities, typically with azimuthal wavenumber m = 1–6. Consistent with previous theoretical work, the propensity for the cellular instabilities is shown to increase with decreasing reactant Lewis number and Damköhler number.     

Response of non-premixed flames to bulk flow perturbations

This paper describes the dynamics of non-premixed flames responding to bulk velocity fluctuations, and compares the dynamics of the flame sheet position and spatially integrated heat release to that of a premixed flame. The space–time dynamics of the non-premixed flame sheet in the fast chemistry limit is described by the stoichiometric mixture fraction surface, extracted from the solution of the

-equation. This procedure has some analogies to premixed flames, where the premixed flame sheet location is extracted from the G = 0 surface of the solution of the G-equation. A key difference between the premixed and non-premixed flame dynamics, however, is the fact that the non-premixed flame sheet dynamics are a function of the disturbance field everywhere, and not just at the reaction sheet, as in the premixed flame problem. A second key difference is that the non-premixed flame does not propagate and so flame wrinkles are convected downstream at the axial flow velocity, while wrinkles in premixed flames convect downstream at a vector sum of the flame speed and axial velocity. With the exception of the flame wrinkle propagation speed, however, we show that that the solutions for the space–time dynamics of the premixed and non-premixed reaction sheets in high velocity axial flows are quite similar. In contrast, there are important differences in their spatially integrated unsteady heat release dynamics. Premixed flame heat release fluctuations are dominated by area fluctuations, while non-premixed flames are dominated by mass burning rate fluctuations. At low Strouhal numbers, the resultant sensitivity of both flames to flow disturbances is the same, but the non-premixed flame response rolls off slower with frequency. Hence, this analysis suggests that non-premixed flames are more sensitive to flow perturbations than premixed flames at O(1) Strouhal numbers.     

The differential diffusion effect of the intermediate species on the stability of premixed flames propagating in microchannels

The propagation of premixed flames in adiabatic and non-catalytic planar microchannels subject to an assisted or opposed Poiseuille flow is considered. The diffusive–thermal model and the well-known two-step chain-branching kinetics are used in order to investigate the role of the differential diffusion of the intermediate species on the spatial and temporal flame stability. This numerical study successfully compares steady-state and time-dependent computations to the linear stability analysis of the problem. Results show that for fuel Lewis numbers less than unity, Le_F < 1, and at sufficiently large values of the opposed Poiseuille flow rate, symmetry-breaking bifurcation arises. It is seen that small values of the radical Lewis number, Le_Z, stabilise the flame to symmetric shape solutions, but result in earlier flashback. For very lean flames, the effect of the radical on the flame stabilisation becomes less important due to the small radical concentration typically found in the reaction zone. Cellular flame structures were also identified in this regime. For Le_F > 1, flames propagating in adiabatic channels suffer from oscillatory instabilities. The Poiseuille flow stabilises the flame and the effect of Le_Z is opposite to that found for Le_F < 1. Small values of Le_Z further destabilise the flame to oscillating or pulsating instabilities.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

Flame structure analysis and flamelet progress variable modelling of strained coal flames

Strained two-phase pulverised coal flames in a counterflow configuration are investigated numerically. Three operating conditions with different coal-to-primary-air ratios and inlet velocities were evaluated in order to establish different flame regimes. At first, the two-phase flow of the fully resolved reference cases is calculated solving the transport equation for the species and directly evaluating the reaction rates. Different flame structures are identified using the heat release rate and the chemical explosive mode as markers, showing that complex structures with a combination of lean premixed and non-premixed flames can be observed in strained counterflow coal flames. In addition to the fully resolved simulation, the suitability of the Flamelet-Progress Variable (FPV) model is investigated. Both premixed and non-premixed tables are employed. At first, the suitability of the look-up tables is evaluated by means of an a priori analysis, using the fully resolved simulations as reference solutions, showing that the non-premixed flamelet table correctly predicts the structure of the strained coal flames, while the premixed table shows sensible deviations in terms of temperature and species, especially at rich conditions. Finally, the a posteriori analysis shows that the fully coupled FPV model with a non-premixed flamelet look-up table can accurately predict strained coal flames.     

Diffusion flame instabilities in a 0.75 mm non-premixed microburner

We examine the cellular instabilities of laminar non-premixed diffusion flames that arise in a polycrystalline alumina microburner with a channel wall gap of dimension 0.75 mm. Changes in the flame structure are observed as a function of the fuel type (H₂, CH₄, and C₃H₈) and diluent. The oxidizer is O₂/inert. In contrast to previous observations on laminar diffusion flame instabilities, the current instabilities occur in the direction of flow above the splitter plate, and only occur for the heavier fuel types. They are not observed in a H₂–O₂ mixture, which will only support a continuous laminar flame inside our burner, regardless of the initial mixture strength and whether or not the flame is in near-quenching conditions. The only exception is when helium is added to the H₂–O₂ mixture, raising the effective Lewis numbers of both components.     

Reaction zone structures and mixing characteristics of partially premixed swirling CH4/air flames in a gas turbine model combustor

The mixing, reaction progress, and flame front structures of partially premixed flames have been investigated in a gas turbine model combustor using different laser techniques comprising laser Doppler velocimetry for the characterization of the flow field, Raman scattering for simultaneous multi-species and temperature measurements, and planar laser-induced fluorescence of CH for the visualization of the reaction zones. Swirling CH₄/air flames with Re numbers between 7500 and 60,000 have been studied to identify the influence of the turbulent flow field on the thermochemical state of the flames and the structures of the CH layers. Turbulence intensities and length scales, as well as the classification of these flames in regime diagrams of turbulent combustion, are addressed. The results indicate that the flames exhibit more characteristics of a diffusion flame (with connected flame zones) than of a uniformly premixed flame.     

On transient behavior of non-premixed counter-flow diffusion flames within the REDIM based model reduction concept

In the present work non-stationary behavior of the counter-flow diffusion flame is examined in the context of the recently developed approach of model reduction called REaction–DIffusion Manifolds (REDIM) method. It is a natural extension of the ILDM approach which takes into account both the chemical reaction and the diffusion processes. It has been developed to treat both premixed and non-premixed regimes of combustion. In this work we investigate the ability of the concept to describe transient processes of extinction and re-ignition. A very simple flame configuration and transport model are considered in this current study for the sake of transparency because the main focus is on the transient and non-stationary behavior of flames. H₂/O₂/N₂ combustion system is considered in a non-premixed counter-flow diffusion 1D flame configuration. This study shows how the REDIM concept performs in the transient regimes; it interprets the effect of local extinction and reigniting phenomena using detailed and reduced models. It shows how the unstable/transient behavior of a detailed system can be successfully accounted with the help of the REDIM based reduced model.     

An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

Effect of Soret diffusion on lean hydrogen/air flames at normal and elevated pressure and temperature

The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H₂ while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H₂ is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H₂ on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H₂ caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.     

On the stability of one-dimensional diffusion flames established between plane,parallel, porous walls

A one-dimensional, non-premixed flame stability analysis is undertaken.Oscillatory and cellular flame instabilities are identified by a careful studyof the numerically calculated eigenvalues of the linearized system of equations. The numerical investigation details the critical locations for changes in flame behaviour, as well as the critical values of variousparameters that affect flame stability. A critical Lewis number, greaterthan unity, is identified as the value where unstable oscillations maybegin to appear (Le?>?Le _c) and for which cellular flames can exist(Le?<?Le _c). Some prior discussions are clarified regarding theaforementioned critical values, as well as the role of convection inproducing flame instabilities. The methodology of the stability analysis isdiscussed in detail.     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

Effects of fuel Lewis number on flame spread over solids

Using a detailed two-dimensional numerical model, a systematic investigation has been made to study the effect of fuel Lewis number (Le_F = α/D_F) and mass transfer on flame spread over thin solids. The fuel Lewis number affects the flame spread rates for both concurrent and opposed flames over thin fuels. The dependence of the flame spread rate on Le_F for these two spreading modes is, however, not the same. In opposed flame spreads (zero-gravity, self-propagation, and normal gravity downward propagation), the flame spread rate vs. Le_F curve is non-monotonic with a maximum value occurring at an intermediate value of Le_F = 0.5. In steady, concurrent spread in zero-gravity with low-speed flow and a constant flame length, the flame spread rate decreases with Le_F in a monotonic manner. By using the computational model as a tool, the effects of fuel mass diffusion perpendicular to and parallel with the solid surface are isolated to obtain more physical insight on the two-dimensional aspect of fuel mass transfer on flame spread. In addition, the model has also been used to decouple the solid evaporation process so that the fuel diffusion effect in the gas-phase can be isolated. Both of these theoretical exercises contribute to the understanding of mass transfer effects on the flame spreading phenomena over solids.     

The effect of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame

The influence of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame was investigated by numerical simulation using three different transport property calculation methods. One simulation included preferential diffusion and the other two neglected preferential diffusion. The results show that the neglect of preferential diffusion or the use of unity Lewis number for all species results in a significant underprediction of soot volume fraction. The peak soot volume fraction is reduced from 8.0 to 2.0 ppm for the studied flame when preferential diffusion is neglected in the simulation. Detailed examination of numerical results reveals that the underprediction of soot volume fraction in the simulation neglecting preferential diffusion is due to the slower diffusion of some species from main reaction zone to PAH and soot formation layer. The slower diffusion of these species causes lower PAH formation rate and thus results in lower soot inception rate and smaller particle surface area. The smaller surface area further leads to smaller surface growth rate. In addition, the neglect of preferential diffusion also leads to higher OH concentration in the flame, which causes the higher specific soot oxidation rate. The lower inception rate, smaller surface growth rate and higher specific oxidation rate results in the lower soot volume fraction when preferential diffusion is neglected. The finding of the paper implies the importance of preferential diffusion for the modeling of not only laminar but maybe also some turbulent flames.     

Simulations of laminar non-premixed flames of kerosene with hot combustion products as oxidiser

The effects of hot combustion product dilution in a pressurised kerosene-burning system at gas turbine conditions were investigated with laminar counterflow flame simulations. Hot combustion products from a lean (φ = 0.6) premixed flame were used as an oxidiser with kerosene surrogate as fuel in a non-premixed counterflow flame at 5, 7, 9 and 11 bar. Kerosene-hot product flames, referred to as ‘MILD’, exhibit a flame structure similar to that of kerosene–air flames, referred to as ‘conventional’, at low strain rates. The Heat Release Rate (HRR) of both conventional and MILD flames reflects the pyrolysis of the primary and intermediate fuels on the rich side of the reaction zone. Positive HRR and OH regions in mixture fraction space are of similar width to conventional kerosene flames, suggesting that MILD flames are thin fronts. MILD flames do not exhibit typical extinction behaviour, but gradually transition to a mixing solution at very high rates of strain (above A = 160, 000 s^?1 for all pressures). This is in agreement with literature that suggests heavily preheated and diluted flames have a monotonic S-shaped curve. Despite these differences in comparison with kerosene–air flames, MILD flames follow typical trends as a function of both strain and pressure. Further still, the peak locations of the overlap of OH and CH₂O mass fractions in comparison with the peak HRR indicate that the pixel-by-pixel product of OH- and CH₂O-PLIF signals is a valid experimental marker for non-premixed kerosene MILD and conventional flames.     

LIF measurements and chemical kinetic analysis of methylidyne formation in high-pressure counter-flow partially premixed and non-premixed flames

We report quantitative, spatially resolved, linear laser-induced fluorescence (LIF) measurements of methylidyne concentration ([CH]) in laminar, methane–air, counter-flow partially premixed and non-premixed flames using excitation near 431.5 nm in the A–X (0,0) band. For partially premixed flames, fuel-side equivalence ratios (_B) of 1.45, 1.6 and 2.0 are studied at pressures of 1, 3, 6, 9 and 12 atm. For non-premixed flames, the fuel-side mixture consists of 25% CH₄ and 75% N₂; measurements are obtained at pressures of 1, 2, 3, 4, 5, 6, 9 and 12 atm. The quantitative CH measurements are compared with predictions from an opposed-flow flame code utilizing two GRI chemical kinetic mechanisms (versions 2.11 and 3.0). LIF measurements of [CH] are corrected for variations in the quenching rate coefficient by using major species concentrations and temperatures generated by the code along with suitable quenching cross sections for CH available from the literature. A pathway analysis provides relative contributions from important elementary reactions to the total amount of CH produced at various pressures. Key reactions controlling peak CH concentrations are also identified by using a sensitivity analysis. For the partially premixed flames, measured CH profiles are reproduced reasonably well by GRI 3.0, although some quantitative disagreement exists at all pressures. Two CH radical peaks are observed for _B=1.45 and _B=1.6 at pressures above 3 atm. Peak CH concentrations for the non-premixed flames are significantly underpredicted by GRI 3.0. The latter agrees with previously reported NO concentrations, which are also underpredicted in these same high-pressure counter-flow diffusion flames. PACS 07.35.+k; 42.62.Fi; 82.33.Vx