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.     

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.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.     

Self-similar behavior and chemistry tabulation of burnt-gas diluted premixed flamelets including heat-loss

In many combustion systems, the reactive gases feeding the reaction zones are diluted by burnt products, to favor flame stabilization, homogenize the temperature distribution and reduce pollutant emission. The objective of this paper is to discuss a premixed flamelet detailed chemistry tabulation strategy for vitiated and non-adiabatic combustion. Dilution by burnt products is parameterized here with two controlling quantities: the amplitude of the heat-loss in the burnt gases, for instance at walls, and the level of reactant vitiation. The chemical response of premixed flames to variations of these parameters is studied and it is shown that most chemical properties of burnt-gas diluted flames feature self-similar behavior, which can be used to dramatically downsize chemical tables based on canonical flamelets. The self-similar behavior of the flamelets is studied for both molecular diffusion and chemical source budgets in a progress variable composition space. It is found that two different scaling relations are needed to ensure self-similar behavior of both major and radical species.     

湍流扩散火焰的非稳态火焰面模拟

采用稳态的和非稳态的火焰面模型同时对一个湍流甲烷射流扩散火焰进行了数值模拟,比较了两者对湍流平均火焰结构、活性自由基和污染物(氮氧化物)排放的模拟效果。速度场采用κ-ε模型计算,守恒标量混合物分数的分布通过其概率密度函数(PDF)输运方程的求解得到。稳态的火焰面结构由查询火焰面数据库得到,而非稳态的火焰面结构由火焰面方程和流场方程耦合求解来计算。采用详细的GRI—Mech 3.0机理描述甲烷的氧化和氮氧化物的形成。数值模拟结果和实验数据作了广泛的对比,验证了火焰面模型对湍流扩散燃烧的定量模拟能力。     

Unstrained and strained flamelets for LES of premixed combustion

The unstrained and strained flamelet closures for filtered reaction rate in large eddy simulation (LES) of premixed flames are studied. The required sub-grid scale (SGS) PDF in these closures is presumed using the Beta function. The relative performances of these closures are assessed by comparing numerical results from large eddy simulations of piloted Bunsen flames of stoichiometric methane–air mixture with experimental measurements. The strained flamelets closure is observed to underestimate the burn rate and thus the reactive scalars mass fractions are under-predicted with an over-prediction of fuel mass fraction compared with the unstrained flamelet closure. The physical reasons for this relative behaviour are discussed. The results of unstrained flamelet closure compare well with experimental data. The SGS variance of the progress variable required for the presumed PDF is obtained by solving its transport equation. An order of magnitude analysis of this equation suggests that the commonly used algebraic model obtained by balancing source and sink in this transport equation does not hold. This algebraic model is shown to underestimate the SGS variance substantially and the implications of this variance model for the filtered reaction rate closures are highlighted.     

A self-contained progress variable space solution method for thermochemical variables and flame speed in freely-propagating premixed flamelets

Flamelet models for premixed combustion, which are based on equations formulated and solved in progress variable space, have been proposed in the past, but have not been adopted for chemistry reduction methods. This is due to one limitation of these models: they need a closure for both the magnitude and the shape of the gradient (or scalar dissipation rate) of the progress variable, which is essential for an accurate prediction of the flame displacement speed. So far, solution methods for the aforementioned models require gradient information as an input, which is either modelled and non-generic, or extracted from a previous physical space flame solution for the analogous problem. The objective of this work is to provide a self-contained solution method for freely-propagating premixed flamelets in progress variable space, by solving an additional flamelet equation for the gradient of the progress variable. With this, the novel method provides both magnitude and shape of the gradient. Studying hydrogen-air and methane-air configurations, it is demonstrated that an accurate prediction of the laminar flame speed without the necessity for further input parameters can be obtained.     

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.     

Regularization of reaction progress variable for application to flamelet-based combustion models

Many combustion models that are based on the flamelet paradigm employ a reaction progress variable. While such a progress variable is well defined for one-step reaction kinetics, this is typically not the case for complex chemical mechanisms. Consequently, several expressions for a progress variable have been utilized. In this paper a formal method for the generation of a reaction progress variable is proposed that is optimal with respect to a set of constraints. The potential of the method is demonstrated in applications to partially premixed and diffusion flames, and the extension to premixed combustion is discussed. It is shown that the proposed method can lead to significant improvements in the definition of an optimal progress variable over conventional formulations, essentially eliminating the expert knowledge previously required in identifying such quantities.     

A DNS assessment of linear relations between filtered reaction rate,flame surface density,and scalar dissipation rate in a weakly turbulent premixed flame

Linear relations between (i) filtered reaction rate and filtered flame surface density (FSD) and (ii) filtered reaction rate and filtered scalar dissipation rate (SDR), which are widely used in Large Eddy Simulation (LES) research into premixed turbulent combustion, are examined by processing DNS data obtained from a statistically 1D planar flame under weakly turbulent conditions that are most favourable for the two approaches (flamelet combustion regime, single-step chemistry, equidiffusive mixture, adiabatic burner, and low Mach number). The analysis well supports the former approach provided that the filtered reaction rate is combined with filtered molecular transport term. In such a case, both the RANS and LES FSD approaches are based on local relations valid within weakly perturbed flamelets. Accordingly, simply recasting RANS expressions to a filtered form works well. On the contrary, while the FSD and SDR approaches appear to be basically similar at first glance, the analysis does not support the latter one, but shows that a ratio of the filtered reaction rate to the filtered SDR is strongly scattered within the studied flame brush, with its conditionally mean value varying significantly with Favre-filtered combustion progress variable. As argued in the paper, these limitations of the LES SDR approach stem from the fact that it is based on a relation valid after integration over weakly perturbed flamelets, but this relation does not hold locally within such flamelets. Consequently, when a sufficiently small filter is applied to instantaneous fields, the filter may contain only a part of the local flamelet, whereas the linear relation holds solely for the entire flamelet and may not hold within the filtered flamelet volume. Thus, the present study implies that straightforwardly recasting well-established RANS equations to a filtered form is a flawed approach if the equations are based on integral features of local burning.     

Effects of turbulence-chemistry interactions on auto-ignition and flame structure for n-dodecane spray combustion

The Engine Combustion Network (ECN) spray A under diesel engine conditions is investigated with a non-adiabatic 5D Flamelet Generated Manifolds (FGM) model with the consideration of detailed chemical kinetic mechanisms. The enthalpy deficit due to droplet vapourisation is considered by employing an additional controlling parameter in the FGM library. In this FGM model, β-PDF is used for the PDF integration over the control variable space. Validation results in non-reacting conditions indicate relatively good agreement between the predicted and experimental data in terms of liquid and vapour penetrations and mixture fraction spatial distribution. In reacting conditions, the effects of variance of mixture fraction and progress variable were examined. The ignition delay time and the quasi-steady flame structure are both affected by the variances. The variance of mixture fraction delays the ignition process and the variance of progress variable accelerates it. For mixture fraction, the ignition process is quicker at any stage in the case of neglecting variance. While things are more complex for progress variable, the ignition process is advanced in the case of neglecting variance at early times, but surpassed by the case of β-PDF later and until auto-ignition. When variance of mixture fraction is considered, the OH mass fraction shows a wide spatial distribution. While if not, a very thin flame is observed with a higher peak in OH, and a very large lift-off length. The variance of progress variable has little impact on the global flame structure, but makes the flame lift-off length much shorter. This study confirms the general observation, that the variance of mixture fraction is of higher importance in high temperature non-premixed combustion, however, we found that the variance of progress variable is far from negligible.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

Influences of stoichiometry on steadily propagating triple flames in counterflows

Most studies of triple flames in counterflowing streams of fuel and oxidizer have been focused on the symmetric problem in which the stoichiometric mixture fraction is 1/2. There then exist lean and rich premixed flames of roughly equal strengths, with a diffusion flame trailing behind from the stoichiometric point at which they meet. In the majority of realistic situations, however, the stoichiometric mixture fraction departs appreciably from unity, typically being quite small. With the objective of clarifying the influences of stoichiometry, attention is focused on one of the simplest possible models, addressed here mainly by numerical integration. When the stoichiometric mixture fraction departs appreciably from 1/2, one of the premixed wings is found to be dominant to such an extent that the diffusion flame and the other premixed flame are very weak by comparison. These curved, partially premixed flames are expected to be relevant in realistic configurations. In addition, a simple kinematic balance is shown to predict the shape of the front and the propagation velocity reasonably well in the limit of low stretch and low curvature.     

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

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

Capturing multi-regime combustion in turbulent flames with a virtual chemistry approach

Multiple flame regimes are encountered in industrial combustion chambers, where premixed, stratified and non-premixed flame regions may coexist. To obtain a predictive tool for pollutant formation predictions, chemical flame modeling must take into account the influence of such complex flame structure. The objective of this article is to apply and compare two reduced chemistry models on both laminar and turbulent multi-regime flame configurations in order to analyze their capabilities in predicting flame structure and CO formation. The challenged approaches are (i) a premixed flamelet-based tabulated chemistry method, whose thermochemical variables are parameterized by a mixture fraction and a progress variable, and (ii) a virtual chemical scheme which has been optimized to retrieve the properties of canonical premixed and non-premixed 1-D laminar flames. The methods are first applied to compute a series of laminar partially-premixed methane-air counterflow flames. Results are compared to detailed chemistry simulations. Both approaches reproduced the thermal flame structure but only the virtual chemistry captures the CO formation in all ranges of equivalence ratio from stoichiometry premixed flame to pure non-premixed flame. Finally, the two chemical models combined with the Thickened Flame model for LES are challenged on a piloted turbulent jet flame with inhomogeneous inlet, the Sydney inhomogeneous burner. Mean and RMS of temperature and CO mass fraction radial profiles are compared to available experimental data. Scatter data in mixture fraction space and Wasserstein metric of numerical and experimental data are also studied. The analyses confirm again that the virtual chemistry approach is able to account for the impact of multi-regime turbulent combustion on the CO formation.     

A quantitative method for a priori evaluation of combustion reaction models

In recent years, direct numerical simulations have been used increasingly to evaluate the validity and performance of combustion reaction models. This study presents a new, quantitative method to determine the ideal model performance attainable by a given parameterization of the state variables. Data from direct numerical simulation (DNS) of unsteady CO/H₂–air jet flames is analysed to determine how well various parameterizations represent the data, and how well specific models based on those parameterizations perform. Results show that the equilibrium model performs poorly relative to an ideal model parameterized by the mixture fraction. The steady laminar flamelet model performs quite well relative to an ideal model parameterized by mixture fraction and dissipation rate in some cases. However, at low dissipation rates or at dissipation rates exceeding the steady extinction limit, the steady flamelet model performs poorly. Interestingly, even in many cases where the steady flamelet model fails (particularly at low dissipation rate), the DNS data suggests that the state may be parameterized well by the mixture fraction and dissipation rate. A progress variable based on the CO₂ mass fraction is proposed, together with a new model based on the CO₂ progress variable. This model performs nearly ideally, and demonstrates the ability to capture extinction with remarkable accuracy for the CO/H₂ flames considered.