Flame propagation in a small-scale parallel flow

We consider the propagation of laminar premixed flames in the presence of a parallel flow whose scale is smaller than the laminar flame thickness. The study addresses fundamental aspects with relevance to flame propagation in narrow channels, to the emerging micro-combustion technology, and to the understanding of the effect of small scales in a (turbulent) flow on the flame structure. In part, the study extends the results of a previous analytical study carried out in the thick flame asymptotic limit which has in particular addressed the validity of Damköhler's second hypothesis in the context of laminar steady parallel flows. Several new contributions are made here.

Analytical contributions include the derivation of an explicit formula for the effective speed of a premixed flame U _T in the presence of an oscillatory parallel flow whose scale ? (measured with the laminar flame thickness δ _L ) is small and amplitude A (measured with the laminar flame speed U _L ) is (1). The formula shows a quadratic dependence on both the amplitude and the scale of the flow. The validity of the formula is established analytically in two distinguished limits corresponding to (1) frequencies of oscillations (measured with the natural frequency of the flame U _L /δ _L ), and to higher frequencies of (A/?) (the natural frequency of the flow). The analytical study yields partial support of Damköhler's second hypothesis in that it shows that the flame behaves as a planar flame (to leading order) with an increased propagation speed which depends on both the scale and amplitude of the velocity fluctuation. However our formula for U _T contradicts the formula given by Damköhler in his original paper where U _T has a square root dependence on the scale and amplitude.

Numerical contributions include a significant set of two-dimensional calculations which determine the range of validity of the asymptotic findings. In particular, these account for volumetric heat loss and differential diffusion effects. Good agreement between the numerics and asymptotics is found in all cases, both for steady and oscillatory flows, at least in the expected range of validity of the asymptotics. The effect of the frequency of oscillation is also discussed. Additional related aspects such as the difference in the response of thin and thick flames to the combined effect of heat loss and fluid flow are also addressed. It is found for example that the sensitivity of thick flames to volumetric heat loss is negligibly affected by the parallel flow intensity, in marked contrast to the sensitivity of thin flames. Interestingly, and somewhat surprisingly, thin flames are found to be more resistant to heat loss when a flow is present, even for unit Lewis number; this ceases to be the case, however, when the Lewis number is large enough.     

Turbulent scalar flux transport in head-on quenching of turbulent premixed flames: a direct numerical simulations approach to assess models for Reynolds averaged Navier Stokes simulations

The statistical behaviour and the modelling of turbulent scalar flux transport have been analysed using a direct numerical simulation (DNS) database of head-on quenching of statistically planar turbulent premixed flames by an isothermal wall. A range of different values of Damköhler, Karlovitz numbers and Lewis numbers has been considered for this analysis. The magnitudes of the turbulent transport and mean velocity gradient terms in the turbulent scalar flux transport equation remain small in comparison to the pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation when the flame is away from the wall but the magnitudes of all these terms diminish and assume comparable values during flame quenching before vanishing altogether. It has been found that the existing models for the turbulent transport, pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation do not adequately address the respective behaviours extracted from DNS data in the near-wall region during flame quenching. Existing models for transport equation-based closures of turbulent scalar flux have been modified in such a manner that these models provide satisfactory prediction both near to and away from the wall.     

Direct numerical simulations of exhaust gas recirculation effect on multistage autoignition in the negative temperature combustion regime for stratified HCCI flow conditions by using H2O2 addition

Direct numerical simulations (DNSs) of a stratified flow in a homogeneous compression charge ignition (HCCI) engine are performed to investigate the exhaust gas recirculation (EGR) and temperature/mixture stratification effects on the autoignition of synthetic dimethyl ether (DME) in the negative temperature combustion region. Detailed chemistry for a DME/air mixture is employed and solved by a hybrid multi-time scale (HMTS) algorithm to reduce the computational cost. The effect of to mimic the EGR effect on autoignition are studied. The results show that adding enhances autoignition by rapid OH radical pool formation (34–46% reduction in ignition delay time) and changes the ignition heat release rates at different ignition stages. Sensitivity analysis is performed and the important reactions pathways affecting the autoignition are specified. The DNS results show that the scales introduced by thermal and mixture stratifications have a strong effect after the low temperature chemistry (LTC) ignition especially at the locations of high scalar dissipation rates. Compared to homogenous ignition, stratified ignitions show similar first autoignition delay times, but 18% reduction in the second and third ignition delay times. The results also show that molecular transport plays an important role in stratified low temperature ignition, and that the scalar mixing time scale is strongly affected by local ignition in the stratified flow. Two ignition-kernel propagation modes are observed: a wave-like, low-speed, deflagrative mode and a spontaneous, high-speed, ignition mode. Three criteria are introduced to distinguish these modes by different characteristic time scales and Damkhöler numbers using a progress variable conditioned by an ignition kernel indicator. The low scalar dissipation rate flame front is characterized by high displacement speeds and high mixing Damkhöler number. The proposed criteria are applied successfully at the different ignition stages and approximate characteristic values are identified to delineate between the different ignition propagation modes.     

Numerical analyses of deflagration initiation by a hot jet

Numerical simulations of axisymmetric reactive jets with one-step Arrhenius kinetics are used to investigate the problem of deflagration initiation in a premixed fuel–air mixture by the sudden discharge of a hot jet of its adiabatic reaction products. For the moderately large values of the jet Reynolds number considered in the computations, chemical reaction is seen to occur initially in the thin mixing layer that separates the hot products from the cold reactants. This mixing layer is wrapped around by the starting vortex, thereby enhancing mixing at the jet head, which is followed by an annular mixing layer that trails behind, connecting the leading vortex with the orifice rim. A successful deflagration is seen to develop for values of the orifice radius larger than a critical value a _c in the order of the flame thickness of the planar deflagration δ_L. Introduction of appropriate scales provides the dimensionless formulation of the problem, with flame initiation characterised in terms of a critical Damköhler number Δ_c=(a _c/δ_L)², whose parametric dependence is investigated. The numerical computations reveal that, while the jet Reynolds number exerts a limited influence on the criticality conditions, the effect of the reactant diffusivity on ignition is much more pronounced, with the value of Δ_c increasing significantly with increasing Lewis numbers . The reactant diffusivity affects also the way ignition takes place, so that for reactants with the flame develops as a result of ignition in the annular mixing layer surrounding the developing jet stem, whereas for highly diffusive reactants with Lewis numbers sufficiently smaller than unity combustion is initiated in the mixed core formed around the starting vortex. The analysis provides increased understanding of deflagration initiation processes, including the effects of differential diffusion, and points to the need for further investigations incorporating detailed chemistry models for specific fuel–air mixtures.     

A posteriori testing of algebraic flame surface density models for LES

In the application of Large Eddy Simulation (LES) to premixed combustion, the unknown filtered chemical source term can be modelled by the generalised flame surface density (FSD) using algebraic models for the wrinkling factor Ξ. The present study compares the behaviour of the various models by first examining the effect of sub-grid turbulent velocity fluctuation on Ξ through a one-dimensional analysis and by the LES of the ORACLES burner (Nguyen, Bruel, and Reichstadt, Flow, Turbulence and Combustion Vol. 82 [2009], pp. 155–183) and the Volvo Rig (Sjunnesson, Nelsson, and Max, Laser Anemometry, Vol. 3 [1991], pp. 83–90; Sjunnesson, Henrikson, and Löfström, AIAA Journal, Vol. 28 [1992], pp. AIAA–92–3650). Several sensitivity studies on parameters such as the turbulent viscosity and the grid resolution are also carried out. A statistically 1-D analysis of turbulent flame propagation reveals that counter gradient transport of the progress variable needs to be accounted for to obtain a realistic flame thickness from the simulations using algebraic FSD based closure. The two burner setups are found to operate mainly within the wrinkling/corrugated flamelet regime based on the premixed combustion diagram for LES (Pitsch and Duchamp de Lageneste, Proceedings of the Combustion Institute, Vol. 29 [2002], pp. 2001–2008) and this suggests that the models are operating within their ideal range. The performance of the algebraic models are then assessed by comparing velocity statistics, followed by a detailed error analysis for the ORACLES burner. Four of the tested models were found to perform reasonably well against experiments, and one of these four further excels in being the most grid-independent. For the Volvo Rig, more focus is placed upon the comparison of temperature data and identifying changes in flame structure amongst the different models. It is found that the few models which largely over-predict velocities in the ORACLES case and volume averaged in a previous a priori DNS analysis (Chakraborty and Klein, Physics of Fluids, Vol. 20 [2008], p. 085108), deliver satisfactory agreement with experimental observations in the Volvo Rig, whereas a few of the other models are only able to capture the experimental data of the Volvo Rig either quantitatively or qualitatively.     

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.     

Quasi steady state and partial equilibrium approximations: their relation and their validity

The quasi steady state and partial equilibrium approximations are analysed in the context of a system of nonlinear differential equations exhibiting multiscale behaviour. Considering systems in the most general and dimensional form , it is shown that both approximations are limiting cases of leading-order asymptotics. Algorithmic conditions are established which guarantee that the accuracy and stability delivered by the two approximations are equivalent to those obtained with leading-order asymptotics. It is shown that the quasi steady state approximation is a limiting case of the partial equilibrium approximation. Algorithms are reported for the identification of the variables in quasi steady state and/or of the processes in partial equilibrium.     

An elementary model for the validation of flamelet approximations in non-premixed turbulent combustion

The fundamental soundness of three flamelet models for non-premixed turbulent combustion is examined on the basis of their performance in an idealized model problem that merges ideas from the laminar asymptotic theory for non-premixed flames and rigorous homogenization theory for the diffusion of a passive scalar. The overall flame configuration is stabilized by a mean gradient in the passive scalar: large Damköhler number asymptotics results are available for the laminar case to quantify the finite-rate effects that cause the flame to depart from its equilibrium state; the same results can also be used to incorporate higher-order corrections in the approximation of the reactive variables in terms of the passive scalar. The use of such flamelet approximations has been extended well beyond the laminar regime as they lie at the core of practical strategies to simulate non-premixed flames in the turbulent regime: the flamelet representation avoids the problem of turbulence closure for the reactive variables by replacing it by the presumably much simpler closure problem for a passive scalar. It is precisely the validity of this substitution outside the laminar regime that is addressed here in the idealized context of a class of small-scale periodic flows for which extensive rigorous results are available for the passive scalar statistics. Results for this simplified problem are reported here for significant wide ranges of Peclet and Damköhler numbers. Asymptotic convergence is observed in terms of the Damköhler number, with a convergence rate that is found to match the laminar predictions and appears relatively insensitive to the Peclet number. The passive scalar dissipation plays a key role in achieving higher-order corrections for the finite-rate case: replacing its pointwise value by an averaged value is convenient practically and can be rigorously motivated for the class of flows studied here, but while it does achieve an overall improvement over the lower-order equilibrium model, the simplification compromises the higher asymptotic convergence observed with the original finite-rate flamelet model with exact local dissipation.(Some figures in this article are in colour only in the electronic version; see www.iop.org)     

Turbulent scalar fluxes in H2-air premixed flames at low and high Karlovitz numbers

The statistical behaviour and closure of sub-grid scalar fluxes in the context of turbulent premixed combustion have been assessed based on an a priori analysis of a detailed chemistry Direct Numerical Simulation (DNS) database consisting of three hydrogen-air flames spanning the corrugated flamelets (CF), thin reaction zones (TRZ) and broken reaction zones (BRZ) regimes of premixed turbulent combustion. The sub-grid scalar fluxes have been extracted by explicit filtering of DNS data. It has been found that the conventional gradient hypothesis model is not appropriate for the closure of sub-grid scalar flux for any scalar in the context of a multispecies system. However, the predictions of the conventional gradient hypothesis exhibit a greater level of qualitative agreement with DNS data for the flame representing the BRZ regime. The aforementioned behaviour has been analysed in terms of the properties of the invariants of the anisotropy tensor in the Lumley triangle. The flames in the CF and TRZ regimes are characterised by a pronounced two-dimensional anisotropy due to strong heat release whereas a three-dimensional and more isotropic behaviour is observed for the flame located in the BRZ regime. Two sub-grid scalar flux models which are capable of predicting counter-gradient transport have been considered for a priori DNS assessment of multispecies systems and have been analysed in terms of both qualitative and quantitative agreements. By combining the latter two sub-grid scalar flux closures, a new modelling strategy is suggested where one model is responsible for properly predicting the conditional mean accurately and the other model is responsible for the correlations between model and unclosed term. Detailed physical explanations for the observed behaviour and an assessment of existing modelling assumptions have been provided. Finally, the classical Bray–Moss–Libby theory for the scalar flux closure has been extended to address multispecies transport in the context of large eddy simulations.     

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.     

Modelling of turbulent scalar flux in turbulent premixed flames based on DNS databases

A transport equation for scalar flux in turbulent premixed flames was modelled on the basis of DNS databases. Fully developed turbulent premixed flames were obtained for three different density ratios of flames with a single-step irreversible reaction, while the turbulent intensity was comparable to the laminar burning velocity. These DNS databases showed that the countergradient diffusion was dominant in the flame region. Analyses of the Favre-averaged transport equation for turbulent scalar flux proved that the pressure related terms and the velocity–reaction rate correlation term played important roles on the countergradient diffusion, while the mean velocity gradient term, the mean progress variable gradient term and dissipation terms suppressed it. Based on these analyses, modelling of the combustion-related terms was discussed. The mean pressure gradient term and the fluctuating pressure term were modelled by scaling, and these models were in good agreement with DNS databases. The dissipation terms and the velocity–reaction rate correlation term were also modelled, and these models mimicked DNS well.