Development of a Multi-Phase Flamelet Generated Manifold for Spray Combustion Simulations

Author(s):  
Xu Zhang ◽  
Ran Yi ◽  
C. P. Chen

Abstract In this study, a model flame of quasi-1D counterflow spray flame has been developed. The two-dimensional multiphase convection-diffusion-reaction (CDR) equations have been simplified to one dimension using similarity reduction under the Eulerian framework. This model flame is able to directly account for non-adiabatic heat loss as well as multiple combustion regimes present in realistic spray combustion processes. A spray flamelet library was generated based on the model flame. To retrieve data from the spray flamelet library, the enthalpy was used as an additional controlling variable to represent the interphase heat transfer, while the mixing and chemical reaction processes were mapped to the mixture fraction and the progress variable. The spray-flamelet/progress-variable (SFPV) approach was validated against the results from the direct integration of finite-rate chemistry as a benchmark. The SFPV approach gave a better performance in terms of temperature predictions, while the conventional gas-phase flamelet/progress-variable (FPV) approach over-predicted by nearly 20%. In terms of species mass fractions, there was no significant difference between the two, both showing good agreements with the direct integration of chemistry (DIC) model.

2021 ◽  
Vol 143 (6) ◽  
Author(s):  
Xu Zhang ◽  
Ran Yi ◽  
C. P. Chen

Abstract In this study, a model flame of quasi-one-dimensional (1D) counterflow spray flame has been developed. The two-dimensional (2D) multiphase convection-diffusion-reaction equations have been simplified to one dimension using similarity reduction under the Eulerian framework. This model flame is able to directly account for nonadiabatic heat loss, preferential evaporation, as well as multiple combustion regimes present in realistic spray combustion processes. A spray flamelet library was generated based on the model flame. To retrieve data from the spray flamelet library, the enthalpy was used as an additional controlling variable to represent the interphase heat transfer, while the mixing and chemical reaction processes were mapped to the mixture fraction and the progress variable. The spray flamelet generated manifolds (SFGM) approach was validated against the results from the direct integration of finite rate chemistry as a benchmark. The SFGM approach was found to give a better performance in terms of predictions of temperature and species mass fractions.


Author(s):  
Graham Goldin ◽  
Yongzhe Zhang

The Flamelet Generated Manifold (FGM) model requires a reaction progress variable which is usually defined as a weighted sum of species mass fractions. This progress variable should increase monotonically as flamelet states progress from unburnt to chemical equilibrium. A favorable attribute of the progress variable is that the flamelet species should change gradually with the progress variable, which reduces sensitivity of these species to any predicted errors in the progress variable. Previous publications have presented optimization algorithms for specific flamelet operating conditions, including fuel and oxidizer compositions and temperatures, and pressures. This work applies the HEEDS optimization software to find optimal species weights for a range of fuels and operating conditions. The fuels included are methane, methane-hydrogen, n-dodecane and n-heptane, at fuel-oxidizer temperatures of 293K and 1000K, and pressures of 1 and 30 atmospheres. For manifolds modeled by constant pressure ignition reactors, the optimal progress variable weights using four species weights are {αCO2 = 1, αCO = 0.91, αH2O = 0.52, αH2 = 1}, and for eight species weights are {αCO2 = 1, αCO = 0.91, αH2O = 0.51, αH2 = 1, αC2H2 = 0.16, αOH = −0.66, αH = −0.38, αO = 0.4}.


Author(s):  
Chitralkumar V. Naik ◽  
Hossam Elasrag ◽  
Rakesh Yadav ◽  
Ahad Validi ◽  
Ellen Meeks

Abstract Combustion models can have a significant impact on flame simulations. While solving finite rate chemistry typically yields more accurate predictions, they depend significantly on the detailed kinetics mechanism used. To demonstrate the effect, Large Eddy Simulation (LES) of Sandia Flame D [1] has been performed using various combustion models. Four different detailed kinetics mechanisms have been considered. They include DRM mechanism with 22 species, GRI-mech 2.11 with 49 species, GRI-mech 3.0 with 53 species [2], and Model Fuel Library (MFL) mechanism with 29 species [3]. In addition to the mechanisms, two modeling approaches considered are direct integration of finite rate kinetics (FR) and Flamelet Generated Manifold (FGM). The performance is compared between combinations of the mechanisms and combustion-modeling approaches for prediction of the flame structure and pollutants, including NO and CO. The mesh contains about half a million hexahedral cells and LES statistics were collected over ten flow throughs. Advanced solvers including dynamic cell clustering using the Chemkin-CFD solver in Fluent have been used for faster simulation time. Based on comparison of simulation results to the measurements at various axial and radial positions, we find that the results using the FGM approach were comparable to those using direct integration of FR chemistry, except for NO. In general, the simulation results are in good agreement with the experiment in terms of aerodynamics, mixture fraction and temperature profiles. However, kinetics mechanisms were found to have the most pronounced effect on emissions predictions. NO was especially more sensitive to the kinetics mechanism. Both versions of the GRI-mech fell short in predicting emissions. Overall, the MFL mechanism was found to yield the closest match with the data for flame structure, CO, and NO.


Author(s):  
George Mallouppas ◽  
Graham Goldin ◽  
Yongzhe Zhang ◽  
Piyush Thakre ◽  
Jim Rogerson

Abstract Three Flamelet Generated Manifold reaction source term closure options and two different reactor types are examined with Large Eddy Simulation of an industrial gas turbine combustor operating at 3 bar. This work presents the results for the SGT-100 Dry Low Emission (DLE) gas turbine provided by Siemens Industrial Turbomachinery Ltd. The related experimental study was performed at the German Aerospace Centre, DLR, Stuttgart, Germany. The FGM model approximates the thermo-chemistry in a turbulent flame as that in a simple 0D constant pressure ignition reactors and 1D strained opposed-flow premixed reactors, parametrized by mixture fraction, progress variable, enthalpy and pressure. The first objective of this work is to compare the flame shape and position predicted by these two FGM reactor types. The Kinetic Rate (KR) model, studied in this work, uses the chemical rate from the FGM with assumed shapes, which are a Beta function for mixture fraction and delta functions for reaction progress variable and enthalpy. Another model investigated is the Turbulent Flame-Speed Closure (TFC) model with Zimont turbulent flame speed, which propagates premixed flame fronts at specified turbulent flame speeds. The Thickened Flame Model (TFM), which artificially thickens the flame to sufficiently resolve the internal flame structure on the computational grid, is also explored. Therefore, a second objective of this paper is to compare KR, TFC and TFM with the available experimental data.


Open Physics ◽  
2019 ◽  
Vol 17 (1) ◽  
pp. 905-915
Author(s):  
Mijo Tvrdojevic ◽  
Milan Vujanovic ◽  
Peter Priesching ◽  
Ferry A. Tap ◽  
Anton Starikov ◽  
...  

Abstract Soot prediction for diesel engines is a very important aspect of internal combustion engine emissions research, especially nowadays with very strict emission norms. Computational Fluid Dynamics (CFD) is often used in this research and optimisation of CFD models in terms of a trade-off between accuracy and computational efficiency is essential. This is especially true in the industrial environment where good predictivity is necessary for engine optimisation, but computational power is limited. To investigate soot emissions for Diesel engines, in this work CFD is coupled with chemistry tabulation framework and semi-empirical soot model. The Flamelet Generated Manifold (FGM) combustion model precomputes chemistry using detailed calculations of the 0D homogeneous reactor and then stores the species mass fractions in the table, based on six look-up variables: pressure, temperature, mixture fraction, mixture fraction variance, progress variable and progress variable variance. Data is then retrieved during online CFD simulation, enabling fast execution times while keeping the accuracy of the direct chemistry calculation. In this work, the theory behind the model is discussed as well as implementation in commercial CFD code. Also, soot modelling in the framework of tabulated chemistry is investigated: mathematical model and implementation of the kinetic soot model on the tabulation side is described, and 0D simulation results are used for verification. Then, the model is validated using real-life engine geometry under different operating conditions, where better agreement with experimental measurements is achieved, compared to the standard implementation of the kinetic soot model on the CFD side.


2012 ◽  
Vol 2012 ◽  
pp. 1-11 ◽  
Author(s):  
Kazui Fukumoto ◽  
Yoshifumi Ogami

We present a new combustion simulation technique based on a lookup table approach. In the proposed technique, a flow solver extracts the reaction rates from the look-up table using the mixture fraction, progress variable, and reaction time. Look-up table building and combustion simulation are carried out simultaneously. The reaction rates of the chemical species are recorded in the look-up table according to the mixture fraction, progress variable, and time scale of the reaction. Once the reaction rates are recorded, a direct integration to solve the chemical equations becomes unnecessary; thus, the time for computing the reaction rates is shortened. The proposed technique is applied to an eddy dissipation concept (EDC) model and it is validated through a simulation of a CO-H2-air nonpremixed flame. The results obtained by using the proposed technique are compared with experimental and computational data obtained by using the EDC model with direct integration. Good agreement between our method and the EDC model and the experimental data was found. Moreover, the computation time for the proposed technique is approximately 99.2% lower than that of the EDC model with direct integration.


Author(s):  
O. Colin ◽  
S. Chevillard ◽  
J. Bohbot ◽  
P. K. Senecal ◽  
E. Pomraning ◽  
...  

The current work presents a recent development of the Extended Coherent Flamelet Model (ECFM) for 3D combustion modeling in spark-ignited gasoline engines. The reference-based ECFM model, originally published in 2003, computes the conditional unburned and burned gas species mass fractions from both real species and species tracers. This current work is motivated by two limitations of the reference-based model. First, the difference between convection of species tracers and convection of real species leads to small discrepancies between the two, due to high velocity gradients during gas exchange. This can lead to inaccurate estimation of the progress variable and consequently to negative conditional mass fractions in the burned gases after ignition. Second, the reference-based ECFM model assumes implicitly that the unburned and burned states correspond to the same mixture fraction. This assumption is valid for low stratification cases, but it can lead to substantial conditioning errors for highly stratified systems like gasoline direct injection (GDI) engines. To address these shortcomings, a new species-based ECFM (SB-ECFM) implementation is presented. In this species-based model, the unburned and burned gas states are entirely defined by the transported species in each zone. It is shown that SB-ECFM more reliably defines conditional quantities and the progress variable. This enhancement allows the use of a second-order central scheme in space when using full decoupling of auto-ignition and premixed flame progress variables as proposed in Robert et al., Proc. Comb. Inst, 2015, while the reference model is limited to the first-order upwind scheme in this case. Finally, simulations of a GDI engine are presented at different loads and rpm conditions. It is shown that, with the higher order scheme, SB-ECFM demonstrates very good agreement with measured pressure.


Author(s):  
Rakesh Yadav ◽  
Shaoping Li ◽  
Ellen Meeks

In this work, a scale separation method has been proposed and implemented in the framework of Flamelet Generated Manifold (FGM) model. In this approach, first a list of slow evolving species like NO, N2O etc., are identified. Then, a separate transport equation for each of these species (called FGM scalars) is solved in addition to the mixture fraction and progress variable equations. The forward and reverse reaction rates of these slow forming species are computed in two-dimensional FGM flamelets and pre-tabulated as a function of progress variable, mixture fraction and their respective variances. At run time, the pre-tabulated probability density function (PDF) averaged production rates of these FGM scalars are used, while their tabulated reverse rates are modified with a linear scaling based on the ratio of tabulated values of the FGM scalar and the prevailing values of the FGM scalars from three dimensional CFD solution. This mechanism allows the reverse rates to provide continuous feedback and respond to the slow evolution of scalar. Other than the list of selected scalars, all other species and temperature are still computed as a function of the main progress variable and mixture fraction. Since, a small set of scalars can be used to track key species, this methodology remains computationally efficient. The current approach has been implemented into commercial CFD solver, ANSYS Fluent, and has been validated for two lab scale turbulent flames, the first one is Sandia Flame D, while the second one is a lifted turbulent methane flame in vitiated co-flow. In the current work, two additional FGM scalar transport equations are solved for CO and NO and comparisons have been made against the tabulated values as well as the experimental data. It has been seen that the scale separation methodology of these scalars leads ∼10–15% improvements in the CO mass fraction, while it reduces the peak NO formation up to 4 times leading to better agreement with experimental data compared to tabulated values. The quality of predictions from the current method is also evaluated against finite rate chemistry-based model as well as reduced order NO model. It is found that the current model has consistent results, and is an improvement over current reduced order modeling approach.


Sign in / Sign up

Export Citation Format

Share Document