A flamelet analysis of the burning velocity of premixed turbulent expanding flames

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the corrugated flamelet regime are performed. The flamelet-generated manifold method is used to deal with detailed reaction kinetics. The numerical method is validated for both laminar and turbulent expanding flames. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. If the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is predicted accurately. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed to be constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate.     

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.     

Evaluation of models for flame stretch due to curvature in the thin reaction zones regime

Direct numerical simulation (DNS) was used to study modelling assumptions for the curvature-propagation component of flame stretch in the thin reaction zones regime of turbulent premixed combustion, a regime in which small eddies can penetrate the preheat zone but not the thinner fuel breakdown zone. Simulations of lean hydrogen–air and methane–air flames were conducted, and statistics of flame stretch due to curvature, henceforth referred to simply as stretch, were extracted from a species mass fraction iso-surface taken to represent the flame. The study focussed on investigating the modelling assumptions of Peters [J. Fluid Mech. 384 (1999) 107]. It was found that the mean stretch is dominated by stretch due to correlations of flame speed with curvature, and specifically the effects of tangential diffusion. The modelling suggestions of Peters were found to provide an improvement over the assumptions of a constant flame speed or a flame speed governed by the linear relationship with stretch at small and steady stretch. However for the conditions considered here, diffusive-thermal effects remain well into the thin reaction zones regime, and the suggestions of Peters generally over-predict the mean compressive stretch. An effective diffusivity for flame stretch was suggested and evaluated for the methane simulations. It was found that the effective diffusivity was comparable to the mass diffusivity for flames with a high ratio of flame time to eddy turnover time. The length scales contributing to stretch were investigated, and it was found that while most flame area has a radius of curvature greater than the laminar flame thickness, most stretch occurs in more tightly curved flame elements.     

How fast can we burn, 2.0

We review the state of the art in measurements and simulations of the behavior of premixed laminar and turbulent flames, subject to differential diffusion, stretch and curvature. The first part of the paper reviews the behavior of premixed laminar flames subject to flow stretch, and how it affects the accuracy of measurements of unstrained laminar flame speeds in stretched and spherically propagating flames. We then examine how flow field stretch and differential diffusion interact with flame propagation, promoting or suppressing the onset of thermodiffusive instabilities. Secondly, we survey the methodology for and results of measurements of turbulent flame speeds in the light of theory, and identify issues of consistency in the definition of mean flame speeds, and their corresponding mean areas. Data for methane at a single operating condition are compared for a range of turbulent conditions, showing that fundamental issues that have yet to be resolved for Bunsen and spherically propagating flames. Finally, we consider how the laminar flame scale response of flames to flow perturbations interacting with differential diffusion leads to very different outcomes to the overall sensitivity of the burning rate to turbulence, according to numerical simulations (DNS). The paper concludes with opportunities for future measurements and model development, including the perennial recommendation for robust archival databases of experimental and DNS results for future testing of models.     

湍流分层燃烧的直接数值模拟和小火焰模型研究

湍流分层燃烧广泛应用于工业燃烧装置,但是目前还比较缺乏适用于湍流分层燃烧的高精度数值模型。本文利用直接数值模拟数据库,对高Karlovitz数分层射流火焰的小火焰模型表现进行了先验性评估。考虑了两种小火焰模型,一种是基于自由传播层流预混火焰的小火焰模型M1,另一种是基于分层对冲小火焰的小火焰模型M2。研究发现M1和M2在c-Z空间的结果与直接数值模拟在定性上是一致的。在物理空间,M2对过程变量反应速率脉动值的预测结果要优于M1.     

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.     

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.     

Stratified turbulent Bunsen flames: flame surface analysis and flame surface density modelling

In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification.

All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1?φ?1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate.

In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1?φ?1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4?φ?1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1?φ?1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data, thereby rendering this approach promising for application in LES.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

Modeling curvature effects in turbulent autoigniting non-premixed flames using tabulated chemistry

Turbulent flames are intrinsically curved. In the presence of preferential diffusion, curvature effects either enhance or suppress molecular diffusion, depending on the diffusivity of the species and the direction of the flame curvature. When a tabulated chemistry type of modeling is employed, curvature-preferential diffusion interactions have to be taken into consideration in the construction of manifolds. In this study, we employ multistage stage flamelet generated manifolds (MuSt-FGM) method to model autoigniting non-premixed turbulent flames with preferential diffusion effects included. The conditions for the modeled flame are in MILD combustion regime. To model the above-mentioned curvature-preferential diffusion interactions, a new mixture fraction which has a non-unity Lewis number is defined and used as a new control variable in the manifold generation. 1D curved flames are simulated to create the necessary flamelets. The resulting MuSt-FGM tables are used in the simulation of 1D laminar flames, and then also applied to turbulent flames using 2D direct numerical simulations (DNS). It was observed that when the curvature effects are included in the manifold, the MuSt-FGM results agree well with the detailed chemistry results; whereas the results become unsatisfactory when the curvature effects are ignored.     

Recent developments in DNS of turbulent combustion

The simulation of turbulent flames fully resolving the smallest flow scales and the thinnest reaction zones goes along with specific requirements, which are discussed from dimensionless numbers useful to introduce the generic context in which direct numerical simulation (DNS) of turbulent flames is performed. Starting from this basis, the evolution of the DNS landscape over the past five years is reviewed. It is found that the flow geometries, the focus of the studies and the overall motivations for performing DNS have broadened, making DNS a standard tool in numerical turbulent combustion. Along these lines, the emerging DNS of laboratory burners for turbulent flame modeling development is discussed and illustrated from DNS imbedded in Large Eddy Simulation (LES) and flow resolved simulation of bluff-body flames. The literature shows that DNS generated databases constitute a fantastic playground for developing and testing a large spectrum of promising machine learning methods for the control and the optimisation of combustion systems, including novel numerical approaches based on the training of neural networks and which can be evaluated in DNS free from sub-model artefacts. The so-called quasi-DNS is also progressively entering the optimisation loop of combustion systems, with the application of techniques to downsize real combustion devices in order to perform fully resolved simulations of their complex geometries. An example of such study leading to the improvement of an incinerator efficiency is reported. Finally, numbers are given relative to the carbon footprint of the generation of DNS databases, motivating the crucial need for community building around database sharing.     

Dynamics of triple-flames in ignition of turbulent dual fuel mixture: A direct numerical simulation study

Pilot-ignited dual fuel combustion involves a complex transition between the pilot fuel autoignition and the premixed-like phase of combustion, which is challenging for experimental measurement and numerical modelling, and not sufficiently explored. To further understand the fundamentals of the dual fuel ignition processes, the transient ignition and subsequent flame development in a turbulent dimethyl ether (DME)/methane-air mixing layer under diesel engine-relevant conditions are studied by direct numerical simulations (DNS). Results indicate that combustion is initiated by a two-stage autoignition that involves both low-temperature and high-temperature chemistry. The first stage autoignition is initiated at the stoichiometric mixture, and then the ignition front propagates against the mixture fraction gradient into rich mixtures and eventually forms a diffusively-supported cool flame. The second stage ignition kernels are spatially distributed around the most reactive mixture fraction with a low scalar dissipation rate. Multiple triple flames are established and propagate along the stoichiometric mixture, which is proven to play an essential role in the flame developing process. The edge flames gradually get close to each other with their branches eventually connected. It is the leading lean premixed branch that initiates the steady propagating methane-air flame. The time required for the initiation of steady flame is substantially shorter than the autoignition delay time of the methane-air mixture under the same thermochemical condition. Temporal evolution of the displacement speed at the flame front is also investigated to clarify the propagation characteristics of the combustion waves. Cool flame and propagation of triple flames are also identified in this study, which are novel features of the pilot-ignited dual fuel combustion.     

Towards a classification of models for the numerical simulation of premixed combustion based on a generalized regime diagram

The inner structure, and the physical behaviour of turbulent premixed flames are usually described, and classified by means of the regime diagram introduced by Borghi and Peters. Thereby properties related to both the flame and the (turbulent) flow are considered. In this work a diagram valid for all physical regimes, comprising suitable requirements for laminar simulations, direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaging based numerical simulation (RANS) is proposed. In particular the diagram describes essential situations within the validity limits of the “Borghi, Peters diagram” which physical phenomena are resolved by the simulation, and which have to be modelled. This information is used for systematic classification of various models by suggesting specific models that are appropriate depending on the regime and numerical resolution, and may provide guidance for numerical simulation methods and model development in turbulent premixed combustion. This might help users as a guideline in choosing appropriate models for a given device, and numerical effort available. The regime diagram suggested by Pitsch and Duchamp de Lageneste, which includes DNS and LES by explicitely accounting for the numerical related variable filterwidth, emerges here as one of the special two-dimensional cases possible. In contrast to the generalized regime diagram, their diagram does not include laminar simulations, and RANS based considerations, while transition between wrinkled and corrugated flamelets is not clearly established.     

An algorithm for LES of premixed compressible flows using the Conditional Moment Closure model

A novel numerical method has been developed to couple a recent high order accurate fully compressible upwind method with the Conditional Moment Closure combustion model. The governing equations, turbulence modelling and numerical methods are presented in full. The new numerical method is validated against direct numerical simulation (DNS) data for a lean premixed methane slot burner. Although the modelling approaches are based on non-premixed flames and hence not expected to be valid for a wide range of premixed flames, the predicted flame is just 10% longer than that in the DNS and excellent agreement of mean mass fractions, conditional mass fractions and temperature is demonstrated. This new numerical method provides a very useful framework for future application of CMC to premixed as well as non-premixed combustion.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

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.     

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.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

Structural response of different Lewis number premixed flames interacting with a toroidal vortex

Simultaneous measurements of temperature, CH* and OH* chemiluminescent species are carried out to explore the impact of stretch rate and curvature on the structure of premixed flames. The configuration of an initially flat premixed flame interacting with a toroidal vortex is selected for the present study and reasons for this choice are discussed. Lewis number effects are assessed by comparing methane and propane flames. It is emphasized that the flame structure experiences very strong variations. In particular, the flame is shrunk (broadened) in the initial (final) period of the interaction with the vortex where strain rate (curvature) contribution of the stretch rate is predominant. By further analysing independently the thickness of the preheat and reaction zones, it is shown that for propane flames, not only the former but also the latter is significantly altered in zones where the flame curvature is negative. Changes in the reaction zone properties are further emphasized using CH* and OH* radicals. It is demonstrated that higher thermal diffusivity plays a significant role around curved regions, in which the enhanced diffusion of heat leads to a strong increase of CH* compared to OH* intensity. As an overall conclusion, this study suggests that it would be interesting to reassess the internal flame structure at lower and moderate Karlovitz numbers since changes might appear for a moderate vortex intensity with typical size much larger than the flame thickness.