Subgrid modeling of intrinsic instabilities in premixed flame propagation

This work is devoted to the investigation and subgrid-scale modeling of intrinsic flame instabilities occurring in the propagation of a deflagration wave. Such instabilities, of hydrodynamic and thermodiffusive origin, are expected to be of particular relevance in recent technological trends such as in the use of hydrogen as a clean energy carrier or as a secondary fuel in hydrogen enriched combustion. A dedicated set of direct numerical simulations is presented and used, in conjunction with coherent literature results, in order to develop scaling arguments for the propagation speed of self-wrinkled flames which are also supported by the outcomes of a weakly non-linear model, namely the Sivashinsky equation. The observed scaling is based on the definition of the number of unstable wavelengths in a reference hydrodynamic lengthscale, in other words the ratio between the neutral or cutoff lengthscale of intrinsic instabilities and the lateral domain of a planar flame. The scalings are then employed to develop an algebraic model for the wrinkling factor in the context of a flame surface density closure approach. An a-priori analysis shows that the model correctly captures the flame wrinkling caused by intrinsic instability at sub grid level. A strategy to include the developed self-wrinkling model in the context of a turbulent combustion model is finally discussed on the basis of the turbulence induced cut-off concept.     

Large scale effects in weakly turbulent premixed flames

In this study we numerically investigate large scale premixed flames in weakly turbulent flow fields. A large scale flame is classified as such based on a reference hydrodynamic lengthscale being larger than a neutral (cutoff) lengthscale for which the hydrodynamic or Darrieus–Landau (DL) instability is balanced by stabilizing diffusive effects. As a result, DL instability can develop for large scale flames and is inhibited otherwise. Direct numerical simulations of both large scale and small scale three-dimensional, weakly turbulent flames are performed at constant Karlovitz and turbulent Reynolds number, using two paradigmatic configurations, namely a statistically planar flame and a slot Bunsen flame. As expected from linear stability analysis, DL instability induces its characteristic cusp-like corrugation only on large scale flames. We therefore observe significant morphological and topological differences as well as DL-enhanced turbulent flame speeds in large scale flames. Furthermore, we investigate issues related to reaction rate modeling in the context of flame surface density closure. Thicker flame brushes are observed for large scale flames resulting in smaller flame surface densities and overall larger wrinkling factors.     

Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames

Lean hydrogen/air flames are prone to hydrodynamic and thermodiffusive instabilities. In this work, the contribution of each instability mechanism is quantified separately by performing detailed simulations of laminar planar lean hydrogen/air flames with different diffusivity models and equations of state to selectively suppress the hydrodynamic or thermodiffusive instability mechanism.From the analysis of the initial phase of the simulations, the thermodiffusive instability is shown to dominate the flame dynamics. If differential diffusion and, hence, the thermodiffusive instability is suppressed, the flame features a strong reduction of the instability growth rates, whereas if present, a wide range of unstable wave numbers is observed due to the strong destabilizing nature of differential diffusion. When instabilities are fully developed, lean hydrogen/air flames feature the formation of small-scale cellular structures and large-scale flame fingers. While the size of the former is known to be close to the most unstable wave length of a linear stability analysis, this work shows that flame fingers also originate from the thermodiffusive instability and most noteworthy, are not linked to an interaction of the two instability mechanisms. They are stable with respect to external perturbations and feature an enhanced flame propagation as the formation of a central cusp at their tip enables the co-existence of two strongly curved leading edges with high reactivity. The thermodiffusive instability is shown to significantly affect the flames’ consumption speed, while the consumption speed enhancement caused by the hydrodynamic instability is significantly smaller. Further, the surface area increase due to wrinkling is strongly diminished if one of the two instability mechanisms is missing. This is linked to a synergistic interaction between the two mechanisms, as the propagation of flame fingers is enhanced by the presence of the hydrodynamic instability due to a widening of the streamlines ahead of the flame fingers.     

Turbulent propagation of premixed flames in the presence of Darrieus–Landau instability

We investigate the role played by hydrodynamic instability in the wrinkled flamelet regime of turbulent combustion, where the intensity of turbulence is small compared to the laminar flame speed and the scale large compared to the flame thickness. To this end the Michelson–Sivashinsky (MS) equation for flame front propagation in one and two spatial dimensions is studied in the presence of uncorrelated and correlated noise representing a turbulent flow field. The combined effect of turbulence intensity, integral scale, and an instability parameter related to the Markstein length are examined and turbulent propagation speed monitored for both stable planar flames and corrugated flames for which the planar conformation is unstable. For planar flames a particularly simple scaling law emerges, involving quadratic dependence on intensity and a linear dependence on the degree of instability. For corrugated flames we find the dependence on intensity to be substantially weaker than quadratic, revealing that corrugated flames are more resilient to turbulence than planar flames. The existence of a threshold turbulence intensity is also observed, below which the corrugated flame in the presence of turbulence behaves like a laminar flame. We also analyze the conformation of the flame surface in the presence of turbulence, revealing primary, large-scale wrinkles of a size comparable to the main corrugation. When the integral scale is much smaller than the characteristic corrugation length we observe, in addition to primary wrinkles, secondary small-scale wrinkles contaminating the surface. The flame then acquires a multi-scale, self-similar conformation, with a fractal dimension, for one-dimensional flames, plateauing at 1.23 for large intensities. The existence of an intermediate integral scale is also found at which the turbulent speed is maximized. When two-dimensional flames are subject to turbulence, the primary wrinkling patterns give rise to polyhedral-cellular structures which bear a very close resemblance to those observed in experiments on hydrodynamically unstable expanding spherical flames.     

A numerical study of auto-ignition in turbulent lifted flames issuing into a vitiated co-flow

This paper presents a numerical study of auto-ignition in simple jets of a hydrogen–nitrogen mixture issuing into a vitiated co-flowing stream. The stabilization region of these flames is complex and, depending on the flow conditions, may undergo a transition from auto-ignition to premixed flame propagation. The objective of this paper is to develop numerical indicators for identifying such behavior, first in well-known simple test cases and then in the lifted turbulent flames. The calculations employ a composition probability density function (PDF) approach coupled to the commercial CFD code, FLUENT. The in-situ-adaptive tabulation (ISAT) method is used to implement detailed chemical kinetics. A simple k–ε turbulence model is used for turbulence along with a low Reynolds number model close to the solid walls of the fuel pipe.

The first indicator is based on an analysis of the species transport with respect to the budget of convection, diffusion and chemical reaction terms. This is a powerful tool for investigating aspects of turbulent combustion that would otherwise be prohibitive or impossible to examine experimentally. Reaction balanced by convection with minimal axial diffusion is taken as an indicator of auto-ignition while a diffusive–reactive balance, preceded by a convective–diffusive balanced pre-heat zone, is representative of a premixed flame. The second indicator is the relative location of the onset of creation of certain radical species such as HO₂ ahead of the flame zone. The buildup of HO₂ prior to the creation of H, O and OH is taken as another indicator of autoignition.

The paper first confirms the relevance of these indicators with respect to two simple test cases representing clear auto-ignition and premixed flame propagation. Three turbulent lifted flames are then investigated and the presence of auto-ignition is identified. These numerical tools are essential in providing valuable insights into the stabilization behaviour of these flames, and the demarcation between processes of auto-ignition and premixed flame propagation.     

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.     

Temporal evolution of flame stretch due to turbulence and the hydrodynamic instability

The temporal evolution of the strain rate on a turbulent premixed flame was measured experimentally using cinema-stereoscopic particle image velocimetry. Turbulence strains a flame due to velocity gradients associated both directly with the turbulence and those caused by the hydrodynamic instability, which are initiated by the turbulence. The development of flame wrinkles caused by both of these mechanisms was observed. Wrinkles generated by the turbulence formed around vortical structures, which passed through the flame and were attenuated. After the turbulent structures had passed, the hydrodynamic instability flow pattern developed and caused additional strain. The hydrodynamic instability also caused the growth of small flame front perturbations into large wrinkles. In the moderately turbulent flame investigated, it was found that the evolution of the strain rate caused by turbulence–flame interactions followed a common pattern involving three temporal regimes. In the first, the turbulence exerted extensive (positive) strain on the flame, creating a wrinkle that had negative curvature (concave towards the reactants). This was followed by a transition period, leading into the third regime in which the flow pattern and strain rate were dominated by the hydrodynamic instability mechanism. It was also found that the magnitudes of the strain rate in the first and third regimes were similar. Hence, the hydrodynamic instability mechanism caused significant strain on a flame and should be included in turbulent combustion models.     

Numerical simulation of flames as gas-dynamic discontinuities

The dynamics of thin premixed flames is computationally studied within the context of a hydrodynamic theory. A level-set method is used to track down the flame, which is treated as a free-boundary interface. The flow field is described by the incompressible Navier–Stokes equations, with different densities for the burnt and unburnt gases, supplemented by singular source terms that properly account for thermal expansion effects. The numerical scheme has been tested on several benchmark problems and was shown to be stable and accurate. In particular, the propagation of a planar flame front and the dynamics of hydrodynamically unstable flames were successfully simulated. This includes recovering the planar front in narrow domains, the Darrieus–Landau linear growth rate for long waves of small amplitude, and the nonlinear development of cusp-like structures predicted by the Michelson–Sivashinsky equation for a small density change. The stationary flame of a Bunsen burner with uniform and parabolic outlet flows were also simulated, showing in particular a careful mapping of the flow field. Finally, the evolution of a hydrodynamically unstable flame was studied for finite amplitude disturbances and realistic values of thermal expansion. These results, which constitute one of the main objectives of this study, elucidate the effect of thermal expansion on flame dynamics.     

Spherical turbulent flame propagation of pulverized coal particle clouds in an O2/N2 atmosphere

The present study aims to clarify the effects of turbulence intensity and coal concentration on the spherical turbulent flame propagation of a pulverized coal particle cloud. A unique experimental apparatus was developed in which coal particles can be dispersed homogeneously in a turbulent flow field generated by two fans. Experiments on spherical turbulent flame propagation of pulverized coal particle clouds in a constant volume spherical chamber in various turbulence intensities and coal concentrations were conducted. A common bituminous coal was used in the present study. The flame propagation velocity was obtained from an analysis of flame propagation images taken using a high-speed camera. It was found that the flame propagation velocity increased with increasing flame radius. The flame propagation velocity increases as the turbulence intensity increases. Similar trends were observed in spherical flames using gaseous fuel. The coal concentration has a weak effect on the flame propagation velocity, which is unique to pulverized coal combustions in a turbulent field. These are the first reports of experimental results for the spherical turbulent flame propagation behavior of pulverized coal particle clouds. The results obtained in the present study are obviously different from those of previous pulverized coal combustion studies and any other results of gaseous fuel combustion research.     

Problems of predicting turbulent burning rates

Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, P_b, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When P_b is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate P_b and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.     

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.     

On polyhedral structures of lean methane/hydrogen Bunsen flames: Combined experimental and numerical analysis

In premixed flame propagation of lean hydrogen or hydrogen-enriched blends, both hydrodynamic and thermo-diffusive instabilities are governing the flame front shape and affect its propagation velocity. As a result, different types of cellular patterns can occur along the flame front in a laminar scenario. In this context, an interesting phenomenon is the formation of polyhedral flames which can be observed in a Bunsen burner. It is the objective of this work to systematically characterize the polyhedral structures of premixed methane/hydrogen Bunsen flames in a combined experimental and numerical study. A series of lean flames with hydrogen content varying between 20 and 85% at two equivalence ratios is investigated. The experiments encompass chemiluminescence imaging together with Planar Laser-induced Fluorescence (PLIF) measurements of the OH radical. Characteristic cell sizes are quantified from the experiments and related to the characteristic length scales obtained from a linear stability analysis. In the experiments, it is observed that the cell sizes at the base of the polyhedral Bunsen flames decrease almost linearly with hydrogen addition and only a weak dependence on the equivalence ratio is noted. These trends are well reflected in the numerical results and the length scale comparison further shows that the wavelength with the maximum growth rate predicted by the linear stability analysis is comparable to the cell size obtained from the experiment. The correlation between the experimental findings and the linear stability analysis is discussed from multiple perspectives considering the governing time and length scales, furthermore drawing relations to previous studies on cellular flames.     

CO2 addition and pressure effects on laminar and turbulent lean premixed CH4 air flames

An experimental study on CH₄–CO₂–air flames at various pressures is conducted by using both laminar and turbulent Bunsen flame configurations. The aim of this research is to contribute to the characterization of fuel lean methane/carbon dioxide/air premixed laminar and turbulent flames at different pressures, by studying laminar and turbulent flame propagation velocities, the flame surface density and the instantaneous flame front wrinkling parameters. PREMIX computations and experimental results indicate a decrease of the laminar flame propagation velocities with increasing CO₂ dilution rate. Instantaneous flame images are obtained by Mie scattering tomography. The image analysis shows that although the height of the turbulent flame increases with the CO₂ addition rate, the flame structure is quite similar. This implies that the flame wrinkling parameters and flame surface density are indifferent to the CO₂ addition. However, the pressure increase has a drastic effect on both parameters. This is also confirmed by a fractal analysis of instantaneous images. It is also observed that the combustion intensity S_T/S_L increases both with pressure and the CO₂ rate. Finally, the mean fuel consumption rate decreases with the CO₂ addition rate but increases with the pressure.     

Propagation of Darrieus–Landau unstable laminar and turbulent expanding flames

The propagation of laminar and turbulent expanding flames subjected to Darrieus–Landau (DL), hydrodynamic instability was experimentally studied by employing stoichiometric H₂/O₂/N₂ flames under quiescent and turbulent conditions performed in a newly developed medium-scale, fan-stirred combustion chamber. In quiescent environment, DL unstable laminar flame exhibits three-stage propagation, i.e. smooth expansion, transition acceleration, and self-similar acceleration. The self-similar acceleration is characterized by a power-law growth of acceleration exponent, α, with normalized Peclet number, which is different from the usually suggested self-similar propagation with a constant α. The imposed turbulence advances the onset of both transition acceleration and self-similar acceleration stages and promotes the strength of flame acceleration as additional wrinkles are invoked by turbulence eddies. A DL–turbulent interaction regime is confirmed to be the classical corrugated flamelets regime. Furthermore, the DL instability significantly facilitates the propagation of expanding flames in medium and even intense turbulence. The development of DL cells is not suppressed by turbulence eddies, and it needs to be considered in turbulent combustion.     

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.     

Conditional moment closure and transient flamelet modelling for detailed structure and NO x formation characteristics of turbulent nonpremixed jet and recirculating flames

This study has been mainly motivated to assess computationally and theoretically the conditional moment closure (CMC) model and the transient flamelet model for the simulation of turbulent nonpremixed flames. These two turbulent combustion models are implemented into the unstructured grid finite volume method that efficiently handles physically and geometrically complex turbulent reacting flows. Moreover, the parallel algorithm has been implemented to improve computational efficiency as well as to reduce the memory load of the CMC procedure. Example cases include two turbulent CO/H₂/N₂ jet flames having different flow timescales and the turbulent nonpremixed H₂/CO flame stabilized on an axisymmetric bluff-body burner. The Lagrangian flamelet model and the simplified CMC formulation are applied to the strongly parabolic jet flame calculation. On the other hand, the Eulerian particle flamelet model and full conservative CMC formulation are employed for the bluff-body flame with flow recirculation. Based on the numerical results, a detailed discussion is given for the comparative performances of the two combustion models in terms of the flame structure and NO _x formation characteristics.     

The unstable behavior of cellular premixed flames induced by intrinsic instability

The unstable behavior of cellular premixed flames induced by intrinsic instability is studied by two-dimensional unsteady calculations of reactive flows. In the present numerical simulation, the compressible Navier–Stokes equation including a one-step irreversible chemical reaction is employed. We consider two basic types of phenomena to account for the intrinsic instability of premixed flames, i.e., hydrodynamic and diffusive-thermal effects. The hydrodynamic effect is caused by the thermal expansion through the flame front; the diffusive-thermal effect is caused by the preferential diffusion of mass versus heat. A disturbance with several wavelength components is superimposed on a planar flame, and the formation of a cellular flame induced by hydrodynamic and diffusive-thermal effects is numerically simulated. After the cellular-flame formation, the combination and division of cells are observed. The behavior of cellular-flame fronts becomes more unstable when the Lewis number is lower than unity, since the diffusive-thermal effect has a great influence on the unstable behavior. The cell size changes with time, and its average is greater than the critical wavelength and becomes smaller by decreasing the Lewis number. The flame velocity of cellular flames depends strongly on the length of computational domain in the direction tangential to the flame front. As the length of computational domain increases, the flame velocity becomes larger. This is because the long-wavelength components of disturbances play an important role in the shape of cellular flames, i.e., in the flame-surface area.     

Modelling of lifted turbulent diffusion flames in a channel mixing layer by the flame hole dynamics

The partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics as an effort to develop a prediction model for the turbulent flame lift off. The essence of the flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of the quenching holes initially created by the local quenching events. The numerical simulation for the flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for a constant-density fuel–air channel mixing layer to obtain the background turbulent flow and mixing fields, from which a time series of two-dimensional scalar-dissipation-rate array is extracted. Subsequently, a Lagrangian simulation of the flame hole random walk mapping, projected to the scalar dissipation rate array, yields a temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. In particular, the probability of encountering the reacting state, while conditioned with the instantaneous scalar dissipation rate, is examined to reveal that the conditional probability has a sharp transition across the crossover scalar dissipation rate, at which the flame edge changes its direction of propagation. This statistical characteristic implies that the flame edge propagation instead of the local quenching event is the main mechanism controlling the partial quenching events in turbulent flames. In addition, the conditional probability can be approximated by a heavyside function across the crossover scalar dissipation rate.     

Effects of Lewis number,density ratio and gravity on burning velocity and conditional statistics in stagnating turbulent premixed flames

DNS is performed to analyse the effects of Lewis number (Le), density ratio and gravity in stagnating turbulent premixed flames. The results show good agreement with those of Lee and Huh (Combustion and Flame, Vol. 159, 2012, pp. 1576–1591) with respect to the turbulent burning velocity, S_T, in terms of turbulent diffusivity, flamelet thickness, mean curvature and displacement speed at the leading edge. In all four stagnating flames studied, a mean tangential strain rate resulting in a mean flamelet thickness smaller than the unstretched laminar flame thickness leads to an increase in S_T. A flame cusp of positive curvature involves a superadiabatic burned gas temperature due to diffusive–thermal instability for an Le less than unity. Wrinkling tends to be suppressed at a larger density ratio, not enhanced by hydrodynamic instability, in the stagnating flow configuration. Turbulence is produced, resulting in highly anisotropic turbulence with heavier unburned gas accelerating through a flame brush by Rayleigh–Taylor instability. Results are also provided on brush thickness, flame surface density and conditional velocities in burned and unburned gas and on flame surfaces to represent the internal brush structures for all four test flames.