Closure of the scalar dissipation rate in the spray flamelet equations through a transport equation for the gradient of the mixture fraction

2019 ◽  
Vol 208 ◽  
pp. 330-350 ◽  
Author(s):  
Hernan Olguin ◽  
Arne Scholtissek ◽  
Sebastian Gonzalez ◽  
Felipe Gonzalez ◽  
Matthias Ihme ◽  
...  
Author(s):  
M. Balthasar ◽  
F. Mauss ◽  
M. Pfitzner ◽  
A. Mack

The modelling of soot formation and oxidation under industrially relevant conditions has made significant progress in recent years. Simplified models introducing a small number of transport equations into a CFD code have been used with some success in research configurations simulating a reciprocating diesel engine. Soot formation and oxidation in the turbulent flow is calculated on the basis of a laminar flamelet library model. The gas phase reactions are modelled with a detailed mechanism for the combustion of heptane containing 89 species and 855 reactions developed by Frenklach and Warnatz and revised by Mauss. The soot model is divided into gas phase reactions, the growth of polycyclic aromatic hydrocarbons (PAH) and the processes of particle inception, heterogeneous surface growth, oxidation and condensation. The first two are modelled within the laminar flamelet chemistry, while the soot model deals with the soot particle processes. The time scales of soot formation are assumed to be much larger than the turbulent time scales. Therefore rates of soot formation are tabulated in the flamelet libraries rather than the soot volume fraction itself. The different rates of soot formation, e.g. particle inception, surface growth, fragmentation and oxidation, computed on the basis of a detailed soot model, are calculated in the mixture fraction / scalar dissipation rate space and further simplified by fitting them to simple analytical functions. A transport equation for the mean soot mass fraction is solved in the CFD-code. The mean rate in this transport equation is closed with the help of presumed probability density functions for the mixture fraction and the scalar dissipation rate. Heat loss due to radiation can be taken into account by including a heat loss parameter in the flamelet calculations describing the change of enthalpy due to radiation, but was not used for the results reported here. The soot model was integrated into an existing commercial CFD code as a post-processing module to existing combustion CFD flow fields and is very robust with high convergence rates. The model is validated with laboratory flame data and using a realistic 3-D BMW Rolls-Royce combustor configuration, where test data at high pressure are available. Good agreement between experiment and simulation is achieved for laboratory flames, whereas soot is overpredicted for the aeroengine combustor configuration by 1–2 orders of magnitude.


2000 ◽  
Vol 124 (1) ◽  
pp. 66-74 ◽  
Author(s):  
M. Balthasar ◽  
F. Mauss ◽  
M. Pfitzner ◽  
A. Mack

The modeling of soot formation and oxidation under industrially relevant conditions has made significant progress in recent years. Simplified models introducing a small number of transport equations into a CFD code have been used with some success in research configurations simulating a reciprocating diesel engine. Soot formation and oxidation in the turbulent flow is calculated on the basis of a laminar flamelet library model. The gas phase reactions are modeled with a detailed mechanism for the combustion of heptane containing 89 species and 855 reactions developed by Frenklach and Warnatz and revised by Mauss. The soot model is divided into gas phase reactions, the growth of polycyclic aromatic hydrocarbons (PAH) and the processes of particle inception, heterogeneous surface growth, oxidation, and condensation. The first two are modeled within the laminar flamelet chemistry, while the soot model deals with the soot particle processes. The time scales of soot formation are assumed to be much larger than the turbulent time scales. Therefore rates of soot formation are tabulated in the flamelet libraries rather than the soot volume fraction itself. The different rates of soot formation, e.g., particle inception, surface growth, fragmentation, and oxidation, computed on the basis of a detailed soot model, are calculated in the mixture fraction/scalar dissipation rate space and further simplified by fitting them to simple analytical functions. A transport equation for the mean soot mass fraction is solved in the CFD code. The mean rate in this transport equation is closed with the help of presumed probability density functions for the mixture fraction and the scalar dissipation rate. Heat loss due to radiation can be taken into account by including a heat loss parameter in the flamelet calculations describing the change of enthalpy due to radiation, but was not used for the results reported here. The soot model was integrated into an existing commercial CFD code as a post-processing module to existing combustion CFD flow fields and is very robust with high convergence rates. The model is validated with laboratory flame data and using a realistic three-dimensional BMW Rolls-Royce combustor configuration, where test data at high pressure are available. Good agreement between experiment and simulation is achieved for laboratory flames, whereas soot is overpredicted for the aeroengine combustor configuration by 1–2 orders of magnitude.


Author(s):  
Hernan Olguin ◽  
Philip Hindenberg ◽  
Eva Gutheil

The paper presents a combined theoretical and numerical study of laminar counterflow mono-disperse spray flames. The numerical model includes a similarity transformation of the two-dimensional governing gas phase equations into a one-dimensional formulation. The reduced computational time enables the use of detailed chemical reaction mechanisms to study the spray flame structure. In particular, the effect of spray evaporation on combustion is investigated by means of numerical simulations. For this purpose, the transport equation of the scalar dissipation rate of the mixture fraction is derived, where the spray evaporation source term is included. Numerical simulations of laminar liquid and gaseous ethanol and combustion products mono disperse spray flames under fuel-rich conditions are presented and discussed. The parametric dependence of the flame structures on strain rate is studied with emphasis on the spray evaporation. Droplet reversal and oscillation are found to dominate the flame structure, and they determine the location of the main reaction zone as well as the profile of the scalar dissipation rate. The study aims to develop a novel spray flamelet model for use in the numerical simulations of turbulent spray combustion with particular emphasis on flameless conditions.


2016 ◽  
Vol 167 ◽  
pp. 60-71 ◽  
Author(s):  
Christoph M. Arndt ◽  
Michael J. Papageorge ◽  
Frederik Fuest ◽  
Jeffrey A. Sutton ◽  
Wolfgang Meier ◽  
...  

Author(s):  
Norio Arai ◽  
Takahisa Yamamoto ◽  
Tomohiko Furuhata

In this study, in order to obtain fundamental data for designing the practical gas turbine combustors in the chemical gas turbine (Ch/GT) combined cycle system and/or other gas turbine systems which utilize low heating value fuel such as coal and biomass gasification syngas, we have simulated low heating value fuel–air turbulent diffusion combustion. The simulated results for the profiles of temperature and species concentrations have been compared with the measured ones. As a reaction model, the flamelet model has been applied to predict the turbulent diffusion combustion characteristics. A flamelet library which indicates relationships for scalar dissipation rate, mixture fraction, and mass fraction of chemical species were given by the calculated results of counterflow diffusion flame using 16 chemical species and 25 elementary reactions. Clipped Gaussian and a log-normal distribution were employed as the PDFs of the mixture fraction and the scalar dissipation rate, respectively. In this model, the k-ε two equations turbulence model is used to describe the turbulent flow field and the radiative heat transfer is calculated by the six-flux model. The governing equations were solved iteratively by the SIMPLE algorithm. By the comparisons of the simulated and measured results, it was clarified that the simulated profiles conformed fairly well to the measured ones in the turbulent combustor, and that this combustion simulation model could predict the low heating value fuel turbulent diffusion combustion characteristics.


Sign in / Sign up

Export Citation Format

Share Document