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.     

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.     

The role of tangential diffusion in evaluating the performance of flamelet models

For non-premixed combustion, the steady laminar flamelet model (SLFM) and flamelet/progress variable approach (FPVA) are two popular methods for tabulating flamelet manifolds. Even if the two methods are used to tabulate and parameterize the same flamelet database, their results sometimes differ in the subsequent simulation. In this work, a novel perspective is provided to assess the performance of the SLFM and FPVA. Both approaches are compared with respect to their capabilities to capture tangential diffusion (TD) of the thermochemical state variables along iso-surfaces of mixture fraction. The relevance of TD effects is identified using generalized flamelet equations and regimes by comparing flamelet solutions with and without TD terms to a FTC (full transport and chemistry) solution of a well-known non-premixed coflow flame. It is found that TD effects can play an important role in entire mixture fraction space, even in the classical flamelet regime. This suggests that the ability to characterize TD effects is an important performance indicator for tabulation strategies. Thereafter, an a priori analysis is conducted comparing the results from the FPVA and SLFM (using the same flamelet database) with the FTC results. The results show that the FPVA is able to more accurately describe the thermochemical state and the flame structure than the SLFM. For a more detailed assessment of the two tabulation strategies, the TD terms reconstructed from the FPVA and SLFM are compared to those from the FTC results. It is found that the FPVA can capture a significant portion of TD effects, while the SLFM can hardly characterize TD effects. This particular capability allows the FPVA to describe chemistry-transport interaction and flame structure more accurately than the SLFM.     

A three mixture fraction flamelet model for multi-stream laminar pulverized coal combustion

A three mixture fraction flamelet model is proposed for multi-stream laminar pulverized coal combustion. The technique of coordinate transformation is utilized to map the flamelet solutions from a unit pyramid space into a unit cubic space to improve the stability of the simulation. The validity of the three mixture fraction flamelet model was assessed on different configurations, including a laminar counterflow pulverized coal/methane flame and a laminar piloted pulverized coal jet flame. The flamelet predictions were compared to the reference results of the detailed chemistry solutions. For the counterflow flame, it was found that the flame temperature and major species mass fractions are correctly predicted by the three mixture fraction flamelet model. However, discrepancies are observed for combustion-mode-sensitive species such as CO and H₂ in the premixed combustion region. The thermo-chemical quantities in the char surface reaction zone cannot be correctly predicted if the mixing between the char off-gas stream and other streams is neglected. For the piloted jet flame, it was shown that the stable thermo-chemical variables can be correctly predicted at the upper and middle stream locations. However, at the downstream location, discrepancies can be observed in certain regions. Overall, the validity of the three mixture fraction flamelet model for multi-stream pulverized coal combustion is confirmed and its performance in turbulent pulverized coal combustion will be tested in future work.     

Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.     

Comparison of well-mixed and multiple representative interactive flamelet approaches for diesel spray combustion modelling

The application of detailed chemistry to the computational fluid dynamics simulation of combustion process in diesel engines has many potentials, including the possibility to predict auto-ignition, diffusion flame structure, stabilisation and soot formation in a wide range of operating conditions, also taking into account the effects of different fuel types. Among the approaches that were proposed over the years, the ones that are mostly used in practical calculations can be divided into two main categories: the first assumes each cell to be a well-stirred reactor, while the second employs the flamelet assumption to describe both auto-ignition and turbulent diffusion flame propagation. Despite the fact that both types of model have been widely validated over the years, a detailed comparison between them appears to be very useful in order to understand better the relevant parameters governing auto-ignition, flame stabilisation and the formation of pollutant emissions. This work is focused on a comparison of two different combustion models that were recently implemented by the authors in an open-source code. The first assumes each cell to be a homogeneous reactor and neglects interaction between turbulence and chemistry, while in the second, multiple laminar flamelets are used to represent the structure of a turbulent diffusion flame. Suitable techniques for online reaction rate tabulation and chemical mechanism reduction are also incorporated, to make the use of bigger mechanisms possible (up to 150 species). The two models are compared and validated by simulating constant-volume diesel combustion in a wide range of operating conditions, including variations of ambient temperature and oxygen concentration. Comparison between the computed and experimental data on flame structure, auto-ignition and flame lift-off enables an understanding of the main relevant differences between the models in the way both auto-ignition and flame stabilisation processes are predicted.     

Combustion and emission modelling of a direct-injection spark-ignition engine by combining flamelet models for premixed and diffusion flames

The combustion and emission production processes of a DISI (direct-injection spark-ignition) engine were modelled by combining flamelet models for premixed and diffusion flames. A new surrogate fuel was proposed to approximate the complicated composition of real gasoline. In contrast to simpler conventional models, the fuel was modelled as a ternary mixture of three hydrocarbons: iso-octane, n-heptane and toluene. Turbulent flame propagation in a partially premixed field was modelled by a premixed flamelet model. The mass fractions of the detailed composition of species in burnt gas were predicted by a diffusion flamelet model. For the pollutant formation modelling, a two-step oxidation of CO and H₂ was used to simulate the secondary diffusion flame. The extended Zeldovich mechanism was used to model NOx formation, while a phenomenological model was used to model soot formation. This model was initially applied to a simple geometry to investigate the fundamentals of the model's behaviour, after which three-dimensional computational fluid dynamic (CFD) simulations were performed in a realistic engine geometry.     

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.     

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

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

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.     

Flame characterisation of gas-assisted pulverised coal combustion using FPV-LES

A multiphase flamelet/progress variable (FPV) model for the large eddy simulation (LES) of gas-assisted pulverised coal combustion (PCC) is developed. The target of the simulation is the Darmstadt turbulent gas-assisted swirling solid fuel combustion chamber. The coal particles are treated as Lagrangian point particles, the position, momentum and energy of which are tracked. The gas phase is described by the low-Mach Navier-Stokes equations alongside the Eulerian transport equations of the governing variables for the FPV model. The set of chemical states of the PCC flame is pre-tabulated in a six-dimensional flamelet table and determined by the mixing of the primary fuel stream, volatiles and char off-gases with the oxidising air, the progress of chemical reactions, the interphase heat transfer, as well as sub-grid scale variations. A presumed $\beta$-PDF approach for the total mixture fraction is applied to capture sub-grid scale effects. The discrete ordinate method (DOM) with the weighted sum of grey gases model (WSGGM) is employed to model radiation. The FPV-LES results are validated against the experimental evidence and a good agreement of the predicted mean and RMS velocities, as well as the mean gas temperature between experiments and simulations is obtained. The contributions of the pilot, volatile and char off-gas fuel streams to the coal flame are analysed. It is found that most regions of the furnace are dominated by either pilot or volatile combustion, while char conversion only occurs in the far downstream and outer furnace regions. The pilot gas dominates the near-wall region inside the quarl, whereas the volatile gas mainly released from small particles dominates a first volatile combustion zone in the interior of the internal recirculation zone. Larger particles heat up more slowly and release their volatile content further downstream, leading to a secondary volatile combustion zone.     

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.     

Prediction of NO in turbulent diffusion flames using Eulerian particle flamelet model

The combustion characteristics for the turbulent diffusion flames using the unsteady flamelet concept have been numerically investigated. The Favre-averaged Navier–Stokes equations are solved by a finite volume method of SIMPLE type that incorporates the laminar flamelet concept with a modified k ? ε turbulence model. The NO formation is estimated by solving the Eulerian particle transport equations in a postprocessing mode. Two test problems are considered: CH₄/H₂/N₂ jet flame and CH₄/H₂ stabilised bluff body flame. The temperature and species profiles are well captured by the flamelet model. Two different chemical mechanisms (GRI 2.11 and 3.0) give nearly identical results for temperature and species except NO. The GRI 3.0 gives significantly higher NO levels compared to the GRI 2.11. This is mainly attributed to the difference in NO formation by the prompt mechanism. The NO formation is sensitive to the number of flamelet particles. The NO levels for two test flames do not change when the flamelet particle number exceeds six.     

A hybrid EDC/Flamelet approach for modelling biomass combustion of grate-firing furnace

This paper proposes a new combustion model for the simulation of biomass combustion. It is developed based on the framework of the well-known Eddy Dissipation Concept (EDC) approach, which has the ability to incorporate chemical kinetics in turbulent reacting flows and thus makes it suitable for modelling gas-phase combustion. However, its high computational cost when using detailed chemistry has made it impractical for modelling large/industrial setups. To address this handicap, the proposed approach decouples the real-time calculation of chemical and mixing processes by importing a pre-calculated steady laminar flamelet library into EDC. The development of this new model is performed based on a modified version of EDC (called Extended EDC), which is capable of modelling the gas-phase of biomass combustion over a wide range of turbulent flow conditions. The proposed model is validated by simulating the well-documented experiment of the piloted jet flames of Barlow and Frank. The performance of the model is then evaluated by simulating a small-scale grate firing biomass furnace. The results show that, overall, the proposed model can be used to model biomass combustion at substantially low computational cost.     

A comparison of different approaches to integrate flamelet tables with presumed-shape PDF in flamelet models for turbulent flames

Flamelet models for turbulent combustion modelling make use of presumed-shape probability density functions (PDFs) for integrating laminar flamelet solutions to obtain an integrated flamelet table that can readily be used for turbulent flame calculations. The existence of non-unique approaches for such an integration has rarely been investigated before. For the first time, this work studies systematically the non-uniqueness of the flamelet table integration approaches. A flamelet model called the flamelet/progress variable model is used in the study, although the issue exists generally in many other flamelet models. Two classes of table integration approaches are investigated, one preserving the laminar flamelet structures during integration and the other not. Three different table integration approaches are examined and compared in detail to provide a thorough understanding of the different approaches. A partially stirred reactor is used as a test case for examining the different approaches. A method based on the transported PDF method is also employed to provide a reference for the assessment of the different flamelet table integration approaches. It is found in general that the flamelet preserving integration approach yields a more reasonable joint PDF of the mixture fraction and the progress variable, and the prediction results are closer to the referenced transported PDF results.     

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.     

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.