skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A hybrid scalar model for sooting turbulent flames

Abstract

A Lagrangian Monte Carlo solution of the joint scalar pdf transport equation for mixture fraction and representative soot properties, coupled with an Eulerian solution for the turbulent flow field and here described as a 'hybrid model,' has been developed. The modeling of soot formation and destruction employs an existing description of the key processes based on two soot variables-the soot volume fraction (or mass concentration) and the particle number density. The gas-phase chemistry is introduced through flamelet-state relationships. The simulation strategy is based on tracing the evolution of reactive stochastic particles within the computational domain. The ensemble of these particles at a fixed location and time then describes the joint scalar pdf. Soot rate equations, represented as functions of mixture fraction, soot mass concentration, and number density, are solved exactly in terms of the scalar values of each individual stochastic particle and the associated gas-phase properties derived from laminar flamelet-state relationships. The solution for the turbulent flow field provides the mean velocity and mixing frequency required for the transport of the stochastic particles in both physical and compositional space, while the Monte Carlo simulation returns the computed mean density field to the CFD code. Density-weighted mean values are approximated bymore » ensemble averages over the scalar values of the stochastic particles in individual computational cells. The principal objective of the hybrid model is the improved treatment of nonlinear soot formation and oxidation, in particular, the capture of the intermittency in the oxidation process associated with the noncoexistence of soot and the principal oxidizing species. Significant computational economies accompany the adoption of the laminar flamelet approach for the source terms in the soot rate equations and the reduced number of scalars computed stochastically. (author)« less

Authors:
;  [1]
  1. School of Engineering, Cranfield University, Bedford (United Kingdom)
Publication Date:
OSTI Identifier:
20727302
Resource Type:
Journal Article
Resource Relation:
Journal Name: Combustion and Flame; Journal Volume: 145; Journal Issue: 1-2; Other Information: Elsevier Ltd. All rights reserved
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; SOOT; MONTE CARLO METHOD; FLAMES; MATHEMATICAL MODELS; TURBULENCE; OXIDATION; HYDROCARBONS

Citation Formats

Aksit, I.M., and Moss, J.B.. A hybrid scalar model for sooting turbulent flames. United States: N. p., 2006. Web. doi:10.1016/j.combustflame.2005.10.010.
Aksit, I.M., & Moss, J.B.. A hybrid scalar model for sooting turbulent flames. United States. doi:10.1016/j.combustflame.2005.10.010.
Aksit, I.M., and Moss, J.B.. Sat . "A hybrid scalar model for sooting turbulent flames". United States. doi:10.1016/j.combustflame.2005.10.010.
@article{osti_20727302,
title = {A hybrid scalar model for sooting turbulent flames},
author = {Aksit, I.M. and Moss, J.B.},
abstractNote = {A Lagrangian Monte Carlo solution of the joint scalar pdf transport equation for mixture fraction and representative soot properties, coupled with an Eulerian solution for the turbulent flow field and here described as a 'hybrid model,' has been developed. The modeling of soot formation and destruction employs an existing description of the key processes based on two soot variables-the soot volume fraction (or mass concentration) and the particle number density. The gas-phase chemistry is introduced through flamelet-state relationships. The simulation strategy is based on tracing the evolution of reactive stochastic particles within the computational domain. The ensemble of these particles at a fixed location and time then describes the joint scalar pdf. Soot rate equations, represented as functions of mixture fraction, soot mass concentration, and number density, are solved exactly in terms of the scalar values of each individual stochastic particle and the associated gas-phase properties derived from laminar flamelet-state relationships. The solution for the turbulent flow field provides the mean velocity and mixing frequency required for the transport of the stochastic particles in both physical and compositional space, while the Monte Carlo simulation returns the computed mean density field to the CFD code. Density-weighted mean values are approximated by ensemble averages over the scalar values of the stochastic particles in individual computational cells. The principal objective of the hybrid model is the improved treatment of nonlinear soot formation and oxidation, in particular, the capture of the intermittency in the oxidation process associated with the noncoexistence of soot and the principal oxidizing species. Significant computational economies accompany the adoption of the laminar flamelet approach for the source terms in the soot rate equations and the reduced number of scalars computed stochastically. (author)},
doi = {10.1016/j.combustflame.2005.10.010},
journal = {Combustion and Flame},
number = 1-2,
volume = 145,
place = {United States},
year = {Sat Apr 15 00:00:00 EDT 2006},
month = {Sat Apr 15 00:00:00 EDT 2006}
}
  • The scalar structure of the overfire (fuel-lean) region of sooting turbulent diffusion flames was investigated, considering ethylene and acetylene burning in air. Measurements and predictions are reported of the mean concentrations of major gas species and mean soot volume fractions. Predictions were based on the conserved-scalar formalism in conjunction with the laminar flamelet approximation. The comparison between predictions and measurements was encouraging, suggesting that state relationships for major gas species, found in laminar diffusion flames, were preserved in the overfire region of the turbulent flames. Measurements also indicated nearly constant soot generation efficiencies from point to point in the overfiremore » region for sufficiently long characteristic residence times to yield nearly universal soot volume fraction state relationships at the same conditions. However, effects attributed to finite-rate chemistry were observed at shorter characteristic residence times, causing spatial variations of soot generation efficiencies in the overfire region, with associated loss of universal soot volume fraction state relationships.« less
  • A comprehensive model for luminous turbulent flames is presented. The model features detailed chemistry, radiation and soot models and state-of-the-art closures for turbulence-chemistry interactions and turbulence-radiation interactions. A transported probability density function (PDF) method is used to capture the effects of turbulent fluctuations in composition and temperature. The PDF method is extended to include soot formation. Spectral gas and soot radiation is modeled using a (particle-based) photon Monte Carlo method coupled with the PDF method, thereby capturing both emission and absorption turbulence-radiation interactions. An important element of this work is that the gas-phase chemistry and soot models that have beenmore » thoroughly validated across a wide range of laminar flames are used in turbulent flame simulations without modification. Six turbulent jet flames are simulated with Reynolds numbers varying from 6700 to 15,000, two fuel types (pure ethylene, 90% methane-10% ethylene blend) and different oxygen concentrations in the oxidizer stream (from 21% O{sub 2} to 55% O{sub 2}). All simulations are carried out with a single set of physical and numerical parameters (model constants). Uniformly good agreement between measured and computed mean temperatures, mean soot volume fractions and (where available) radiative fluxes is found across all flames. This demonstrates that with the combination of a systematic approach and state-of-the-art physical models and numerical algorithms, it is possible to simulate a broad range of luminous turbulent flames with a single model. (author)« less
  • The behaviour of the turbulent scalar flux in premixed flames has been studied using Direct Numerical Simulation (DNS) with emphasis on the effects of Lewis number in the context of Reynolds-averaged closure modelling. A database was obtained from DNS of three-dimensional freely propagating statistically planar turbulent premixed flames with simplified chemistry and a range of global Lewis numbers from 0.34 to 1.2. Under the same initial conditions of turbulence, flames with low Lewis numbers are found to exhibit counter-gradient transport, whereas flames with higher Lewis numbers tend to exhibit gradient transport. The Reynolds-averaged transport equation for the turbulent scalar fluxmore » is analysed in detail and the performance of existing models for the unclosed terms is assessed with respect to corresponding quantities extracted from DNS data. Based on this assessment, existing models which are able to address the effects of non-unity Lewis number on turbulent scalar flux transport are identified, and new or modified models are suggested wherever necessary. In this way, a complete set of closure models for the scalar flux transport equation is prescribed for use in Reynolds-Averaged Navier-Stokes simulations. (author)« less
  • Velocity and scalar statistics of premixed flames stabilized in a novel weak-swirl burner have been studied under a broad range of mixture and turbulence conditions. The nonreacting flows are also investigated to determine the influence of weak-swirl on isotropic turbulence. Deduced from the data are the flame speeds, combustion induced flow acceleration, flame generated turbulence, and flame crossing frequencies. They are compared with previous results obtained in other flame configurations. Three series of experiments are conducted to investigate the evolution of velocity statistics with decreasing equivalence ratios. Under very lean conditions ({phi} = 0.6), the velocity profiles have no characteristicmore » features and are not very different from the corresponding nonreacting profiles. An indicator parameter, K{sub T}, derived from fundamental mixture and flow properties is developed to estimate the relative influences of the turbulence stretch and the mean aerodynamic stretch on combustion induced flow acceleration. The analysis shows that intense turbulence and bulk stretch (i.e., large values of K{sub T}) can obscure and counteract the dynamic effects of combustion.« less
  • The objective of the present work is twofold: (1) extend the coupled Monte Carlo probability density function (PDF)/computational fluid dynamics (CFD) computations to the modeling of turbulent spray flames, and (2) extend the PDF/SPRAY/CFD module to parallel computing in order to facilitate large-scale combustor computations. In this approach, the mean gas phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The PDF transport equation is solved by a Montemore » Carlo method, and the mean gas phase velocity and turbulence fields together with the liquid phase equations are solved by existing state-of-the-art numerical representations. The application of the method to both open as well as confined axisymmetric swirl-stabilized spray flames shows good agreement with the measured data. Preliminary estimates indicate that it is well within reach of today`s modern parallel computer to do a realistic gas turbine combustor simulation within a reasonable turnaround time. The article provides complete details of the overall algorithm, parallelization, and other numerical issues related to coupling between the three solvers.« less