Prediction of local extinction and re-ignition effects in non-premixed turbulent combustion using a flamelet/progress variable approach

The flamelet/progress variable approach (FPVA) has been proposed by Pierce and Moin as a model for turbulent non-premixed combustion in large-eddy simulation. The filtered chemical source term in this model appears in unclosed form, and is modeled by a presumed probability density function (PDF) for the joint PDF of the mixture fraction Z and a flamelet parameter λ. While the marginal PDF of Z can be reasonably approximated by a beta distribution, a model for the conditional PDF of the flamelet parameter needs to be developed. Further, the ability of FPVA to predict extinction and re-ignition has also not been assessed. In this paper, we address these aspects of the model using the DNS database of Sripakagorn et al. It is first shown that the steady flamelet assumption in the context of FPVA leads to good predictions even for high levels of local extinction. Three different models for the conditional PDF of the flamelet parameter are tested in an a priori sense. Results obtained using a delta function to model the conditional PDF of λ lead to an overprediction of the mean temperature, even with only moderate extinction levels. It is shown that if the conditional PDF of λ is modeled by a beta distribution conditioned on Z, then FPVA can predict extinction and re-ignition effects, and good agreement between the model and DNS data for the mean temperature is observed.     

Tabulated chemistry approaches for laminar flames: Evaluation of flame-prolongation of ILDM and flamelet methods

The present study considers the performance of tabulation methods for numerical simulation of complex chemical kinetics in laminar combusting flows and compares their predictions to results obtained by direct calculation. Two tabulation methods are considered: the Flame Prolongation of Intrinsic low-dimensional manifold (FPI) method and Steady Laminar Flamelet Model (SLFM). The FPI method is of current interest as it is a potentially unifying approach capable of dealing with both premixed and non-premixed flames for gaseous fuels. SLFM tabulation methods are popular for non-premixed flames and form a good basis for comparing the performance of the FPI approach. The performance of each method is also evaluated by comparing the results to the direct simulation of the laminar flames using two chemical kinetic schemes: simplified chemistry involving five species and one reaction and detailed chemistry involving 53 species and 325 reaction steps. As part of the evaluation process, the computational cost of each method is also assessed. The laminar flames considered in this study include: freely propagating laminar premixed flames, a two-dimensional axisymmetric methane–air opposed-jet diffusion flame, and a two-dimensional axisymmetric methane–air co-flow diffusion flame. Both tabulation methods are implemented in a parallel adaptive mesh refinement (AMR) framework for solving the complete set of governing partial differential equations. These equations are solved using a fully-coupled finite-volume formulation on body-fitted multi-block quadrilateral mesh. Significant improvements in terms of reduced computational requirements, as measured by both storage and processing time, are demonstrated for the tabulated methods.     

Thermal and chemical effects of differential diffusion in turbulent non-premixed H2 flames

As a renewable fuel, hydrogen (H₂) may play an increasingly important role in the development and control of piston and gas turbine engines to achieve zero carbon emissions. Predictive modeling of H₂-fueled combustion processes requires a clear understanding of differential diffusion (DD) due to the high diffusivity of H₂. On the assumption that turbulent mixing is a far more dominant process than molecular mixing, DD effects are typically neglected in turbulent combustion simulations to reduce modeling complications. While this assumption is reasonable for hydrocarbon fuels, it is less valid for H₂ combustion, where DD is significant. In this work, two three-dimensional direct numerical simulations of temporally evolving turbulent H₂ jet flames with and without considering DD are performed and compared with laminar flamelet solutions to assess DD effects under turbulent conditions. The emphasis is placed on assessing the suitability of classical mixture fraction Z and Bilger mixture fraction Z_Bilger as conditioning variables for non-premixed turbulent combustion modeling through analyzing DD effects on flame structure, chemical reactions, and tangential diffusion (TD). Furthermore, the persistence of DD effects under turbulent conditions and the suitability of a conventional DD parameter are investigated by comparing the turbulent flames to laminar flamelet solutions. It is found that conditioning the thermochemical state on Z_Bilger helps to capture DD effects and mitigate the relative contribution of TD, which gives Z_Bilger advantages over Z when employing flamelet modeling. Due to close coupling between DD and local chemical reactions, DD can affect the turbulent/laminar flames in the form of thermal effects due to the change in flame temperature, chemical effects due to the change in chemical reactions, and transport effects due to multiple species with varying diffusivities that could result in the difference between Z and Z_Bilger. While the transport effects are suppressed, significant chemical and thermal effects of DD still persist under turbulent conditions, which indicates that the DD parameter is probably unsuitable for comprehensively characterizing and assessing DD effects on the structure of turbulent non-premixed flames.     

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.     

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

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

Assessment of the presumed mapping function approach for the stationary laminar flamelet modelling of reacting double scalar mixing layers

This paper assesses the Presumed Mapping Function (PMF) approach in the context of the Stationary Laminar Flamelet Modelling (SLFM) of a reacting Double Scalar Mixing Layer (DSML). Starting from a prescribed Gaussian reference field, the PMF approach provides a presumed Probability Density Function (PDF) for the mixture fraction that is subsequently employed to close the Conditional Scalar Dissipation Rate (CSDR) upon doubly-integrating the homogeneous PDF transport equation. The PMF approach is unique in its ability to yield PDF and CSDR distributions that capture the effect of multiple fuel injections of different composition. This distinct feature overcomes the shortcomings of the classical SLFM closures (the β-distribution for the PDF and the counterflow solution for the CSDR). The current study analyses the impact of the binary (two-stream) and trinary (three-stream) PMF approaches on the structure of laminar flamelets in a DSML formed by the mixing of a fuel stream and an oxidiser stream separated by a pilot. The conditions of a partially-premixed methane/air piloted jet flame are considered. A parametric assessment is performed by varying the local mixing statistics and the findings are compared to those of the classical SLFM approach. Further, the influence of the PMF approach on flamelet extinction and transport by means of differential diffusion is thoroughly investigated. It is shown that the trinary PMF approach captures the influence of the pilot stream as it is capable of yielding bimodal CSDR and trimodal PDF distributions. It is further demonstrated that, when the influence of the pilot is significant, flamelets generated using the trinary CSDR closure extinguish at higher strain levels compared to flamelets generated using the binary and counterflow closures. Lastly, it is shown that the trinary PMF approach can be critical for accurate SLFM computations of DSMLs when differential diffusion effects are important.     

Flame structure analysis and flamelet progress variable modelling of strained coal flames

Strained two-phase pulverised coal flames in a counterflow configuration are investigated numerically. Three operating conditions with different coal-to-primary-air ratios and inlet velocities were evaluated in order to establish different flame regimes. At first, the two-phase flow of the fully resolved reference cases is calculated solving the transport equation for the species and directly evaluating the reaction rates. Different flame structures are identified using the heat release rate and the chemical explosive mode as markers, showing that complex structures with a combination of lean premixed and non-premixed flames can be observed in strained counterflow coal flames. In addition to the fully resolved simulation, the suitability of the Flamelet-Progress Variable (FPV) model is investigated. Both premixed and non-premixed tables are employed. At first, the suitability of the look-up tables is evaluated by means of an a priori analysis, using the fully resolved simulations as reference solutions, showing that the non-premixed flamelet table correctly predicts the structure of the strained coal flames, while the premixed table shows sensible deviations in terms of temperature and species, especially at rich conditions. Finally, the a posteriori analysis shows that the fully coupled FPV model with a non-premixed flamelet look-up table can accurately predict strained coal flames.     

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. These observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.     

Non-premixed turbulent combustion modeling based on the filtered turbulent flamelet equation

In turbulent combustion simulations, the flow structure at the unresolved scale level needs to be reasonably modeled. Following the idea of turbulent flamelet equation for the non-premixed flame case, which was derived based on the filtered governing equations(L. Wang, Combust. Flame 175, 259(2017)), the scalar dissipation term for tabulation can be directly computed from the resolved flowing quantities, instead of solving species transport equations. Therefore, the challenging source term closure for the scalar dissipation or any assumed probability density functions can be avoided;meanwhile the chemical sources are closed by scaling relations. The general principles are discussed in the context of large eddy simulation with case validation. The new model predictions of the bluff-body flame show sufficiently improved results, compared with these from the classic progress-variable approach.     

Evaluation of a flamelet/progress variable approach for pulverized coal combustion in a turbulent mixing layer

A steady flamelet/progress variable (FPV) approach for pulverized coal flames is employed to simulate coal particle burning in a turbulent shear and mixing layer. The configuration consists of a carrier-gas stream of air laden with coal particles that mixes with an oxidizer stream of hot products from lean combustion. Carrier-phase DNS (CP-DNS) are performed, where the turbulent flow field is fully resolved, whereas the coal is represented by Lagrangian point particles. CP-DNS with direct chemistry integration is performed first and provides state-of-the-art validation data for FPV modeling. In a second step the control variables for FPV are extracted from the CP-DNS and used to test if the tabulated manifold can correctly describe the reacting flow (a priorianalysis). Finally a fully coupled a posteriori FPV simulation is performed, where only the FPV control variables are transported, and the chemical state is retrieved from the table and fed back to the flow solver. The a priori results show that the FPV approach is suitable for modeling the complex reacting multiphase flow considered here. The a posteriori data is similarly in good agreement with the reference CP-DNS, although stronger deviations than a priori can be observed. These discrepancies mainly appear in the upper flame (of the present DNS), where premixing and highly unsteady extinction and re-ignition effects play a role, which are difficult to capture by steady non-premixed FPV modeling. However, the present FPV model accurately captures the lower, more stable flame that burns in non-premixed mode.     

An a priori study of different tabulation methods for turbulent pulverised coal combustion

In many practical pulverised coal combustion systems, different oxidiser streams exist, e.g. the primary- and secondary-air streams in the power plant boilers, which makes the modelling of these systems challenging. In this work, three tabulation methods for modelling pulverised coal combustion are evaluated through an a priori study. Pulverised coal flames stabilised in a three-dimensional turbulent counterflow, consisting of different oxidiser streams, are simulated with detailed chemistry first. Then, the thermo-chemical quantities calculated with different tabulation methods are compared to those from detailed chemistry solutions. The comparison shows that the conventional two-stream flamelet model with a fixed oxidiser temperature cannot predict the flame temperature correctly. The conventional two-stream flamelet model is then modified to set the oxidiser temperature equal to the fuel temperature, both of which are varied in the flamelets. By this means, the variations of oxidiser temperature can be considered. It is found that this modified tabulation method performs very well on prediction of the flame temperature. The third tabulation method is an extended three-stream flamelet model that was initially proposed for gaseous combustion. The results show that the reference gaseous temperature profile can be overall reproduced by the extended three-stream flamelet model. Interestingly, it is found that the predictions of major species mass fractions are not sensitive to the oxidiser temperature boundary conditions for the flamelet equations in the a priori analyses.     

钝体驻定湍流扩散火焰的数值研究——燃烧模型比较

分别采用标量联合的概率密度函数方法、稳态火焰面模型、Euler非稳态火焰面模型和基于有限体积/Monte Carlo混合算法的完备PDF模型对钝体驻定的Sydney湍流扩散火焰HM1进行数值模拟,以比较不同燃烧模型的性能,并比较标量联合的概率密度函数方法和Euler非稳态火焰面模型对氮氧化物排放预测的差异.计算结果和实验数据的比较表明,采用概率密度函数方法计算化学反应可以得到更好的结果但计算量较大,而用火焰面模型求解计算量较小,在接近完全燃烧的情形下,其计算结果比较合理.     

Modelling compression ignition engines by incorporation of the flamelet generated manifolds combustion closure

Tabulated chemistry models allow to include detailed chemistry effects at low cost in numerical simulations of reactive flows. Characteristics of the reactive fluid flows are described by a reduced set of parameters that are representative of the flame structure at small scales so-called flamelets. For a specific turbulent combustion configuration, flamelet combustion closure, with proper formulation of the flame structure can be applied. In this study, flamelet generated manifolds (FGM) combustion closure with progress variable approach were incorporated with OpenFOAM® source code to model combustion within compression ignition engines. For IC engine applications, multi-dimensional flamelet look-up tables for counter flow diffusive flame configuration were generated. Source terms of non-premixed combustion configuration in flamelet domain were tabulated based on pressure, temperature of unburned mixture, mixture fraction, and progress variable. A new frozen flamelet method was introduced to link one dimensional reaction diffusion space to multi-dimensional Computational Fluid Dynamics (CFD) physical space to fulfill correct modelling of thermal state of the engine at expansion stroke when charge composition was changed after combustion and reaction rates were subsided. Predictability of the developed numerical framework were evaluated for Sandia Spray A (constant volume vessel), Spray B (light duty optical Diesel engine), and a heavy duty Diesel engine experiments under Reynolds averaged Navier Stokes turbulence formulation. Results showed that application of multi-dimensional FGM combustion closure can comprehensively predict key parameters such as: ignition delay, in-cylinder pressure, apparent heat release rate, flame lift-off , and flame structure in Diesel engines.     

Investigation of lengthscales, scalar dissipation, and flame orientation in a piloted diffusion flame by LES

This work investigates the structure of a diffusion flame in terms of lengthscales, scalar dissipation, and flame orientation by using large eddy simulation. This has been performed for a turbulent, non-premixed, piloted methane/air jet flame (Flame D) at a Reynolds-number of 22,400. A steady flamelet model, which was represented by artificial neural networks, yields species mass fractions, density, and viscosity as a function of the mixture fraction. This will be shown to suffice to simulate such flames. To allow to examine scalar dissipation, a grid of 1.97 × 10⁶ nodes was applied that resolves more than 75% of the turbulent kinetic energy. The accuracy of the results is assessed by varying the grid-resolution and by comparison to experimental data by Barlow, Frank, Karpetis, Schneider (Sandia, Darmstadt), and others. The numerical procedure solves the filtered, incompressible transport equations for mass, momentum, and mixture fraction. For subgrid closure, an eddy viscosity/diffusivity approach is applied, relying on the dynamic Germano model. Artificial turbulent inflow velocities were generated to feature proper one- and two-point statistics. The results obtained for both the one- and two-point statistics were found in good agreement to the experimental data. The PDF of the flame orientation shows the tilting of the flame fronts towards the centerline. Finally, the steady flamelet approach was found to be sufficient for this type of flame unless slowly reacting species are of interest.     

Investigation of the ignition processes of a multi-injection flame in a Diesel engine environment using the flamelet model

In this work, the first flamelet analysis is conducted of a highly resolved DNS of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction. A novel method is proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. The state-of-the-art DNS database generated by Rieth et al. (US National Combustion Meeting, 2019) for a multi-injection flame in a Diesel engine environment is investigated. Three-dimensional flamelets are extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process are investigated by varying the species composition boundary conditions of the transient flamelet equations. Budget analyses of the generalized flamelet equations show that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.     

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.     

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)     

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.     

Unsteady flamelet/progress variable approach for non-premixed turbulent lifted flames

The unsteady flamelet/progress variable approach has been developed for the prediction of a lifted flame to capture the extinction and re-ignition physics. In this work inclusion of the time variant behavior in the flamelet generation embedded in the large eddy simulation technique, allows better understanding of partially premixed flame dynamics. In the process sufficient simulations to generate unsteady laminar flamelets are performed, which are a function of time. These flamelets are used for the generation of the look-up table and the flamelet library is produced. This library is used for the calculation of temperature and other species in the computational domain as the solution progresses. The library constitutes filtered quantities of all the scalars as a function of mean mixture fraction, mixture fraction variance and mean progress variable. Mixture fraction and progress variable distributions are assumed to be β-PDF and δ-PDF respectively. The technique used here is known as the unsteady flamelet progress variable (UFPV) approach. One of the well known lifted flames is considered for the present modeling which shows flame lift-off. The results are compared with the experimental data for the mixture fraction and temperature. Lift off height is predicted from the numerical calculations and compared with the experimentally given value. Comparisons show a reasonably good agreement and the UFPV combustion model appears to be a promising technique for the prediction of lifted and partially premixed flames.     

Large eddy simulation of turbulent premixed combustion using tabulated detailed chemistry and presumed probability density function

A method of chemistry tabulation combined with presumed probability density function (PDF) is applied to simulate piloted premixed jet burner flames with high Karlovitz number using large eddy simulation. Thermo-chemistry states are tabulated by the combination of auto-ignition and extended auto-ignition model. To evaluate the predictive capability of the proposed tabulation method to represent the thermo-chemistry states under the condition of different fresh gases temperature, a-priori study is conducted by performing idealised transient one-dimensional premixed flame simulations. Presumed PDF is used to involve the interaction of turbulence and flame with beta PDF to model the reaction progress variable distribution. Two presumed PDF models, Dirichlet distribution and independent beta distribution, respectively, are applied for representing the interaction between two mixture fractions that are associated with three inlet streams. Comparisons of statistical results show that two presumed PDF models for the two mixture fractions are both capable of predicting temperature and major species profiles, however, they are shown to have a significant effect on the predictions for intermediate species. An analysis of the thermo-chemical state-space representation of the sub-grid scale (SGS) combustion model is performed by comparing correlations between the carbon monoxide mass fraction and temperature. The SGS combustion model based on the proposed chemistry tabulation can reasonably capture the peak value and change trend of intermediate species. Aspects regarding model extensions to adequately predict the peak location of intermediate species are discussed.