PDF modeling and simulation of premixed turbulent combustion

Abstract The use of probability density function (PDF) methods for turbulent combustion simulations is very attractive because arbitrary finite-rate chemistry can be exactly taken into account. PDF methods are well developed for non-premixed turbulent combustion. However, many real flames involve a variety of mixing regimes (non-premixed, partially-premixed and premixed turbulent combustion), and the development of PDF methods for partially-premixed and premixed turbulent combustion turned out to be a very challenging task. The paper shows a promising way to overcome this problem by extending existing PDF methods such that a variety of mixing regimes can be covered. The latter is done by a generalization of the standard scalar mixing frequency model. The generalized scalar mixing frequency model accounts for several relevant processes in addition to velocity-scalar correlations that are represented by the standard model for the mixing of scalars. The suitability of the new mixing frequency model is shown by applications to several premixed turbulent Bunsen flames which cover various regimes ranging from flamelet to distributed combustion. Comparisons to existing concepts focused on the inclusion of reaction effects in mixing frequency models for non-reacting scalars reveal the advantages of the new mixing frequency model. It is worth noting that the same methodology can be used in corresponding filter density function (FDF) methods.

[1]  Stephen B. Pope,et al.  TURBULENT PREMIXED FLAMES , 1987 .

[2]  Stephen B. Pope,et al.  Numerical integration of stochastic differential equations: weak second-order mid-point scheme for application in the composition PDF method , 2003 .

[3]  S. Pope.  PDF methods for turbulent reactive flows , 1985 .

[4]  T. Shih,et al.  A new k-ϵ eddy viscosity model for high reynolds number turbulent flows , 1995 .

[5]  B. Merci,et al.  Flow and Mixing Fields for Transported Scalar PDF Simulations of a Piloted Jet Diffusion Flame (‘Delft Flame III’) , 2005 .

[6]  C. Dopazo,et al.  An approach to the autoignition of a turbulent mixture , 1974 .

[7]  B. Merci,et al.  INTERACTION BETWEEN CHEMISTRY AND MICRO-MIXING MODELING IN TRANSPORTED PDF SIMULATIONS OF TURBULENT NON-PREMIXED FLAMES , 2007 .

[8]  Stephen B. Pope,et al.  The influence of chemical mechanisms on PDF calculations of nonpremixed piloted jet flames , 2005 .

[9]  Roland Borghi,et al.  A unified pdf-flamelet model for turbulent premixed combustion , 2003 .

[10]  Peyman Givi,et al.  Filtered Density Function for Subgrid Scale Modeling of Turbulent Combustion , 2006 .

[11]  Stephen B. Pope,et al.  A mixing model for turbulent reactive flows based on Euclidean minimum spanning trees , 1998 .

[12]  N. Peters.  Laminar flamelet concepts in turbulent combustion , 1988 .

[13]  Stephen B. Pope,et al.  Local Time-Stepping Algorithm for Solving Probability Density Function Turbulence Model Equations , 2002 .

[14]  R. P. Lindstedt,et al.  Transported PDF modeling of high-Reynolds-number premixed turbulent flames , 2006 .

[15]  P. Jenny,et al.  A mixing model for turbulent flows based on parameterized scalar profiles , 2006 .

[16]  B. Merci,et al.  Study of the performance of three micromixing models in transported scalar PDF simulations of a piloted jet diffusion flame ( Delft Flame III ) , 2006 .

[17]  Chong M. Cha,et al.  A subgrid-scale mixing model for large-eddy simulations of turbulent reacting flows using the filtered density function , 2003 .

[18]  M. Herrmann.  Numerical simulation of turbulent Bunsen flames with a level set flamelet model , 2006 .

[19]  Feng Gao,et al.  A large‐eddy simulation scheme for turbulent reacting flows , 1993 .

[20]  W. P. Jones,et al.  Closure of the Reynolds stress and scalar flux equations , 1988 .

[21]  Jay P. Gore,et al.  An evaluation of flame surface density models for turbulent premixed jet flames , 1999 .

[22]  A. Masri,et al.  Pdf calculations of turbulent lifted flames of H2/N2 fuel issuing into a vitiated co-flow , 2004 .

[23]  T. Shih,et al.  A New K-epsilon Eddy Viscosity Model for High Reynolds Number Turbulent Flows: Model Development and Validation , 1994 .

[24]  S. Pope,et al.  Comparative study of micromixing models in transported scalar PDF simulations of turbulent nonpremixed bluff body flames , 2006 .

[25]  F. Williams,et al.  Turbulent Reacting Flows , 1981 .

[26]  Stefan Heinz,et al.  Unified turbulence models for LES and RANS, FDF and PDF simulations , 2007 .

[27]  C. Cha,et al.  A model for the mixing time scale of a turbulent reacting scalar , 2001 .

[28]  Frediano V. Bracco,et al.  A Discussion of Turbulent Flame Structure in Premixed Charges , 1985 .

[29]  S. Heinz,et al.  Reynolds number effects on mixing and reaction in a turbulent pipe flow , 2001 .

[30]  Johannes Janicka,et al.  Large Eddy Simulation of Turbulent Combustion Systems , 2005 .

[31]  M. Sheikhia,et al.  Large eddy simulation of a turbulent nonpremixed piloted methane jet flame (Sandia Flame D) , 2004 .

[32]  Robert W. Bilger,et al.  Scalar gradient and related quantities in turbulent premixed flames , 1997 .

[33]  James J. Riley,et al.  Testing of mixing models for Monte Carlo probability density function simulations , 2005 .

[34]  I. Gökalp,et al.  Scalar dissipation rate modelling in variable density turbulent axisymmetric jets and diffusion flames , 1998 .

[35]  Stephen B. Pope,et al.  Direct numerical simulation of a passive scalar with imposed mean gradient in isotropic turbulence , 1996 .

[36]  Norbert Peters,et al.  The detailed flame structure of highly stretched turbulent premixed methane-air flames , 1996 .

[37]  Nedunchezhian Swaminathan,et al.  Effect of dilatation on scalar dissipation in turbulent premixed flames , 2005 .

[38]  T. Poinsot,et al.  Theoretical and numerical combustion , 2001 .

[39]  Stephen B. Pope,et al.  Velocity-scalar filtered density function for large eddy simulation of turbulent flows , 2003 .

[40]  Otto Zeman,et al.  Modeling Buoyancy Driven Mixed Layers , 1976 .

[41]  H. Pitsch,et al.  Large-eddy simulation of premixed turbulent combustion using a level-set approach , 2002 .

[42]  Stephen B. Pope,et al.  Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation , 1997 .

[43]  Roland Borghi,et al.  Towards an extended scalar dissipation equation for turbulent premixed combustion , 2003 .

[44]  Stephen B. Pope,et al.  PDF calculations of turbulent nonpremixed flames with local extinction , 2000 .

[45]  Heinz Pitsch,et al.  Hybrid large-eddy simulation/Lagrangian filtered-density-function approach for simulating turbulent combustion , 2005 .

[46]  E. O'brien,et al.  The probability density function (pdf) approach to reacting turbulent flows , 1980 .

[47]  Rodney O. Fox,et al.  Computational Models for Turbulent Reacting Flows: References , 2003 .

[48]  J. Janicka,et al.  Closure of the Transport Equation for the Probability Density Funcfion of Turbulent Scalar Fields , 1979 .

[49]  S. Pope,et al.  Filtered density function for large eddy simulation of turbulent reacting flows , 1998 .

[50]  Thierry Mantel,et al.  A new model of premixed wrinkled flame propagation based on a scalar dissipation equation , 1994 .

[51]  Stefan Heinz,et al.  On Fokker–Planck Equations for Turbulent Reacting Flows. Part 2. Filter Density Function for Large Eddy Simulation , 2003 .

[52]  S. Pope,et al.  Velocity filtered density function for large eddy simulation of turbulent flows , 2000 .

[53]  R. Lindstedt,et al.  Finite Rate Chemistry Effects in Turbulent Reacting Flows , 2004 .

[54]  Fachgebiet Strömungsmechanik,et al.  On Fokker – Planck Equations for Turbulent Reacting Flows . Part 2 . Filter Density Function for Large Eddy Simulation , 2003 .

[55]  Stefan Heinz,et al.  On Fokker–Planck Equations for Turbulent Reacting Flows. Part 1. Probability Density Function for Reynolds-Averaged Navier–Stokes Equations , 2003 .

[56]  R. Borghi.  Turbulent combustion modelling , 1988 .

[57]  Stephen B. Pope,et al.  Calculations of premixed turbulent flames by PDF methods , 1987 .