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.     

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

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

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

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.     

A multi-scale simulation method for high Reynolds number wall-bounded turbulent flows

We present an assessment and enhancement of the hybrid two-level large-eddy simulation method (A.G. Gungor and S. Menon, A new two-scale model for large eddy simulation of wall-bounded flows, Prog. Aerosp. Sci. 46 (2010), pp. 28–45), a multi-scale formulation for simulation of high Reynolds number wall-bounded turbulent flows. The assessment of the method is performed by examining role of static and dynamic blending functions used to perform hybridisation of two-level simulation (K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows, J. Comput. Phys. 220 (2006), pp. 290–311; K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows. Part 2: Non-homogeneous flows, J. Comput. Phys. 222 (2007), pp. 673–701) and large-eddy simulation methods. The sensitivity of first- and second-order turbulence statistics to the type of blending functions is investigated by simulating a fully developed turbulent flow in a channel at a friction Reynolds number Re_τ = 395 and comparing the results with those obtained using a direct numerical simulation. The first-order statistics do not show any significant differences for different blending functions, but the second-order statistics show some minor differences. The dynamic evaluation of the hybrid region and the blending function is necessary for non-equilibrium and complex flows where use of a static blending function can lead to inaccurate results. We propose two criteria for the dynamic evaluation; first evaluates extent of the hybrid region based on the subgrid turbulent kinetic energy and the second estimates the blending function based on a characteristic length scale. The computational efficiency of the method is enhanced by incorporating a hybrid programming paradigm where a standard domain decomposition by the message-passing-interface library is combined with the open multi-processing based parallelisation. A further enhancement of the method is achieved by incorporating a closure model for the unclosed hybrid terms in the governing equations, which appear due to hybridisation of two-level- and large-eddy-simulation methods. The model is based on an order of magnitude approximation and a preliminary assessment of the model shows improvement of turbulence statistics when used to simulate turbulent flow in a periodic channel. The assessment and improvements to the multi-scale method make it more suitable for simulation of practical wall-bounded turbulent flows at higher Reynolds number than a conventional large-eddy simulation. This is demonstrated by simulating two representative cases; turbulent flow at high Reynolds number in a periodic channel and flow over a bump placed on the lower surface of a channel, where a relatively coarser computational grid is found to be sufficient for reasonably accurate results.     

An assessment of large eddy simulations of premixed flames propagating past repeated obstacles

This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.     

A priori analysis of sub-grid variance of a reactive scalar using DNS data of high Ka flames

Direct numerical simulations (DNS) of low and high Karlovitz number (Ka) flames are analysed to investigate the behaviour of the reactive scalar sub-grid scale (SGS) variance in premixed combustion under a wide range of combustion conditions (regimes). An order of magnitude analysis is performed to assess the importance of various terms in the variance evolution equation and the analysis is validated using the DNS results. This analysis sheds light on the relative behaviour among turbulent transport and production, scalar dissipation and chemical processes involved in the evolution of the SGS variance at different Ka. The common expectation is that the variance equation shifts from a reaction-dissipation balance at low Ka to a production–dissipation balance at high Ka with diminishing reaction contribution. However, in large eddy simulation (LES), a high Ka alone does not make the reaction term negligible, as the relative importance of the reaction term has a concurrent increase with filter size. The filter size can be relatively large compared with the Kolmogorov length scale in practical LES of high Ka flames, and as a consequence a reaction–production–dissipation balance may prevail in the variance equation even in a high Ka configuration, and this possibility is quantified using the DNS analysis in this work. This has implications from modelling perspectives, and therefore two commonly used closures in LES for the SGS scalar dissipation rate are investigated a priori to estimate the importance of the above balance in LES modelling. The results are explained to highlight the interplay among turbulence, chemistry and dissipation processes as a function of Ka.     

Numerical simulation of turbulent combustion: Scientific challenges

Predictive simulation of engine combustion is key to understanding the underlying complicated physicochemical processes, improving engine performance, and reducing pollutant emissions. Critical issues as turbulence modeling, turbulence-chemistry interaction, and accommodation of detailed chemical kinetics in complex flows remain challenging and essential for high-fidelity combustion simulation. This paper reviews the current status of the state-of-the-art large eddy simulation (LES)/prob-ability density function (PDF)/detailed chemistry approach that can address the three challenging modelling issues. PDF as a subgrid model for LES is formulated and the hybrid mesh-particle method for LES/PDF simulations is described. Then the development need in micro-mixing models for the PDF simulations of turbulent premixed combustion is identified. Finally the different acceleration methods for detailed chemistry are reviewed and a combined strategy is proposed for further development.     

Large eddy simulation of a turbulent nonpremixed piloted methane jet flame (Sandia Flame D)

Large eddy simulation (LES) is conducted of the Sandia Flame D [Proc. Combust. Inst. 27 (1998) 1087, Sandia National Laboratories (2004)], which is a turbulent piloted nonpremixed methane jet flame. The subgrid scale (SGS) closure is based on the scalar filtered mass density function (SFMDF) methodology [J. Fluid Mech. 401 (1999) 85]. The SFMDF is basically the mass weighted probability density function (PDF) of the SGS scalar quantities [Turbulent Flows (2000)]. For this flame (which exhibits little local extinction), a simple flamelet model is used to relate the instantaneous composition to the mixture fraction. The modelled SFMDF transport equation is solved by a hybrid finite-difference/Monte Carlo scheme. This is the first LES of a realistic turbulent flame using the transported PDF method as the SGS closure. The results via this method capture important features of the flame as observed experimentally.     

LES/PDF for premixed combustion in the DNS limit

We investigated the behaviour of the composition probability density function (PDF) model equations used in a large-eddy simulation (LES) of turbulent combustion in the direct numerical simulation (DNS) limit; that is, in the limit of the LES resolution length scale Δ (and the numerical mesh spacing h) being small compared to the smallest flow length scale, so that the resolution is sufficient to perform a DNS. The correct behaviour of a PDF model in the DNS limit is that the resolved composition fields satisfy the DNS equations, and there are no residual fluctuations (i.e. the PDF is everywhere a delta function). In the DNS limit, the treatment of molecular diffusion in the PDF equations is crucial, and both the ‘random-walk’ and ‘mean-drift’ models for molecular diffusion are investigated. Two test cases are considered, both of premixed laminar flames (of thickness δ_L). We examine the solutions of the model PDF equations for these test cases as functions of Δ/δ_L and h/δ_L. Each of the two PDF models has advantages and disadvantages. The mean-drift model behaves correctly in the DNS limit, but it is more difficult to implement and computationally more expensive. The random-walk model does not have the correct behaviour in the DNS limit in that it produces non-zero residual fluctuations. However, if the specified mixing rate Ω normalised by the reaction timescale τ_c is sufficiently large (Ωτ_c ? 1), then the residual fluctuations are less than 10% and the observed flame speed and thickness are close to their laminar values. Away from the DNS limit (i.e. h/δ_L ? 1), the observed flame thickness scales with the mesh spacing h, and the flame speed scales with Ωh. For this case it is possible to construct a non-general specification of the mixing rate Ω such that the flame speed matches the laminar flame speed.     

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.     

DNS assessment of relation between mean reaction and scalar dissipation rates in the flamelet regime of premixed turbulent combustion

The linear relation between the mean rate of product creation and the mean scalar dissipation rate, derived in the seminal paper by K.N.C. Bray [‘The interaction between turbulence and combustion’, Proceedings of the Combustion Institute, Vol. 17 (1979), pp. 223–233], is the cornerstone for models of premixed turbulent combustion that deal with the dissipation rate in order to close the reaction rate. In the present work, this linear relation is straightforwardly validated by analysing data computed earlier in the 3D Direct Numerical Simulation (DNS) of three statistically stationary, 1D, planar turbulent flames associated with the flamelet regime of premixed combustion. Although the linear relation does not hold at the leading and trailing edges of the mean flame brush, such a result is expected within the framework of Bray's theory. However, the present DNS yields substantially larger (smaller) values of an input parameter c_m (or K₂ = 1/(2c_m ? 1)), involved by the studied linear relation, when compared to the commonly used value of c_m = 0.7 (or K₂ = 2.5). To gain further insight into the issue and into the eventual dependence of c_m on mixture composition, the DNS data are combined with the results of numerical simulations of stationary, 1D, planar laminar methane–air flames with complex chemistry, with the results being reported in terms of differently defined combustion progress variables c, i.e. the normalised temperature, density, or mole fraction of CH₄, O₂, CO₂ or H₂O. Such a study indicates the dependence of c_m both on the definition of c and on the equivalence ratio. Nevertheless, K₂ and c_m can be estimated by processing the results of simulations of counterpart laminar premixed flames. Similar conclusions were also drawn by skipping the DNS data, but invoking a presumed beta probability density function in order to evaluate c_m for the differently defined c's and various equivalence ratios.     

Direct numerical simulations of rich premixed turbulent n-dodecane/air flames at diesel engine conditions

Rich premixed turbulent n-dodecane/air flames at diesel engine conditions are analyzed using direct numerical simulations. The conditions correspond to a parametric variation of the Engine Combustion Network Spray A (pressure 60 atm; oxidizer oxygen level and temperature 21% and 900 K, respectively; fuel temperature 363 K). Three simulations with equivalence ratios of 3, 5, and 7 are performed with a Karlovitz number (Ka, based on flame time) of order 100 to match the estimated Ka of the rich premixed combustion region in Spray A. At these conditions, the reference laminar flames exhibit a complex structure which involves both low-temperature chemistry (LTC) and high-temperature chemistry over a wide range of length scales. In the presence of turbulence, the flame structure is strongly affected in physical space and the reaction zone exhibits a very complex structure in which broken, distributed, and thin regions co-exist, especially for the leanest case. However, the contribution of the LTC pathway is only weakly affected by turbulence. In progress variable space, the mean flame structure, including the chemical source terms, is found to match remarkably well that of the corresponding unity Lewis number laminar flame, particularly for the $?=$ 3 and 5 cases. This behavior is attributed to the strong turbulent mixing occurring throughout the flames/reaction zones, which suppresses differential diffusion effects. Nevertheless, large conditional fluctuations around the mean chemical source terms are identified. These are found to correlate very well with radical species mass fractions such as OH. In addition, a similar functional dependence is obtained from counterflow laminar flames. As such, it appears from these results that laminar flame models have a potential to be used to represent the thermochemical state of rich premixed turbulent flames under diesel engine conditions.     

A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry

A hybrid large-Eddy simulation/filtered-density function (LES–FDF) methodology is formulated for simulating variable density turbulent reactive flows. An indirect feedback mechanism coupled with a consistency measure based on redundant density fields contained in the different solvers is used to construct a robust algorithm. Using this novel scheme, a partially premixed methane/air flame is simulated. To describe transport in composition space, a 16-species reduced chemistry mechanism is used along with the interaction-by-exchange with the mean (IEM) model. For the micro-mixing model, typically a constant ratio of scalar to mechanical time-scale is assumed. This parameter can have substantial variations and can strongly influence the combustion process. Here, a dynamic time-scale model is used to prescribe the mixing time-scale, which eliminates the time-scale ratio as a model constant. Two different flame configurations, namely, Sandia flames D and E are studied. Comparison of simulated radial profiles with experimental data show good agreement for both flames. The LES–FDF simulations accurately predict the increased extinction near the inlet and re-ignition further downstream. The conditional mean profiles show good agreement with experimental data for both flames.     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

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.     

Conditional Source-term Estimation using dynamic ensemble selection and parallel iterative solution

A modified version of the Least-Square QR-factorisation (LSQR) algorithm has been implemented in conjunction with Conditional Source-term Estimation (CSE) for lean, turbulent premixed methane–air combustion via Large Eddy Simulation (LES). The iterative solver can reduce computational times by an order of magnitude during the inversion phase of CSE in comparison with the conventional LU-decomposition method. The advantages of iterative and parallel iterative solvers become more prominent as the size of the system increases. The ensemble selection procedure for computing averages within localised regions of the simulation domain has also been updated to a dynamic routine. This allows for more flexible and efficient allocation of computational resources along with reduced input from the user, especially for complex geometries. Preliminary LES calculations have shown that the implementation of an iterative solver and a dynamic ensemble selection algorithm will reduce computational times significantly with negligible error contribution for one-condition CSE, which is applicable to purely premixed or non-premixed turbulent combustion problems. In addition, these algorithms provide the foundation for exceptional computational cost savings for the inversion in two-condition CSE, or Doubly Conditional Source-term Estimation (DCSE), which has shown promise for predicting partially-premixed combustion. Parallel computation of the inverse solution is particularly beneficial to DCSE as the computational cost of the inversion process is considerably larger than in one-condition CSE.     

Comparative analysis of methods for heat losses in turbulent premixed flames using physically-derived reduced-order manifolds

Heat losses have the potential to substantially modify turbulent combustion processes, especially the formation of pollutants such as nitrogen oxides. The chemistry governing these species is strongly temperature sensitive, making heat losses critical for an accurate prediction. To account for the effects of heat loss in large eddy simulation (LES) using a precomputed reduced-order manifold approach, thermochemical states must be precomputed not only for adiabatic conditions but also over a range of reduced enthalpy states. However, there are a number of methods for producing reduced enthalpy states, which invoke different implicit assumptions. In this work, a set of a priori and a posteriori LES studies have been performed for turbulent premixed flames considering heat losses within a precomputed reduced-order manifold approach to determine the sensitivity to the method by which reduced enthalpy states are generated. Two general approaches are explored for generating these reduced enthalpy states and are compared in detail to assess any effects on turbulent flame structure and emissions. In the first approach, the enthalpy is reduced at the boundary of the one-dimensional (1D) premixed flame solution, resulting in a single enthalpy deficit for a single premixed flame solution. In the second approach, a variable heat loss source term is introduced into the 1D flame solutions by mimicking a real heat loss to reduce the post-flame enthalpy. The two approaches are compared in methane–air piloted turbulent premixed planar jet flames with different diluents that maintain a constant adiabatic flame temperature but experience different radiation heat losses. Both a priori and a posteriori results, as well as a chemical pathway analysis, indicate that the manner by which the heat loss is accounted for in the manifold is of secondary importance compared to other model uncertainties such as the chemical mechanism, except in situations where heat loss is unphysically fast compared to the flame time scale. A new theoretical framework to explain this insensitivity is also proposed, and its validity is briefly assessed.     

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.     

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.