A probability density function Eulerian Monte Carlo field method for large eddy simulations: Application to a turbulent piloted methane/air diffusion flame (Sandia D)

Abstract The Eulerian stochastic field method is applied to the solution of the modeled evolution equation for the subgrid joint probability density function (JPDF) of the reacting scalars in a large eddy simulation (LES) of a piloted methane/air diffusion flame (Sandia Flame D). A simple model for subgrid scale (SGS) stresses and fluxes and a global four-step mechanism for combustion are combined in the formulation. Test cases with varying mesh sizes and numbers of stochastic fields were completed. The differences in the results obtained with the two grids were very small and this indicates that the mesh resolution was sufficient. However, incorporation of the JPDF, via the stochastic field solution method, improved the quality of predictions significantly, particularly those quantities related to reaction, such as temperature. Eight stochastic fields were shown to be enough to characterize the influence of SGS fluctuations on filtered species formation rate to reasonable accuracy and at moderate computational cost. With the exceptions of H2 and CO, good agreement between measured and computed mean and RMS profiles of velocity, composition, and temperature was achieved. The discrepancies in H2 and CO concentrations are attributable to limitations in the global chemistry mechanism used in the LES. Overall the results serve to highlight the potential of the Eulerian stochastic field method in LES.

[1]  Stephen B. Pope,et al.  Filtered mass density function for large-eddy simulation of turbulent reacting flows , 1999, Journal of Fluid Mechanics.

[2]  C. Dopazo.  Probability density function approach for a turbulent axisymmetric heated jet. Centerline evolution , 1975 .

[3]  F. C. Lockwood,et al.  The prediction of the fluctuations in the properties of free, round-jet, turbulent, diffusion flames , 1975 .

[4]  ShuChi-Wang,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes, II , 1989 .

[5]  R. Fox,et al.  Hybrid finite-volume/transported PDF simulations of a partially premixed methane–air flame , 2004 .

[6]  W. Jones,et al.  Pdf modeling of finite-rate chemistry effects in turbulent nonpremixed jet flames , 1998 .

[7]  G. Hauke,et al.  Computing reactive flows with a field Monte Carlo formulation and multi-scale methods , 2004 .

[8]  Stephen B. Pope,et al.  Large eddy simulation of a turbulent nonpremixed piloted methane jet flame (Sandia Flame D) , 2004 .

[9]  M. Rogers,et al.  A priori testing of subgrid models for chemically reacting non-premixed turbulent shear flows , 1997, Journal of Fluid Mechanics.

[10]  Robert W. Bilger,et al.  Second-order conditional moment closure for the autoignition of turbulent flows , 1998 .

[11]  R. Barlow,et al.  Effects of turbulence on species mass fractions in methane/air jet flames , 1998 .

[12]  M. Germano,et al.  Turbulence: the filtering approach , 1992, Journal of Fluid Mechanics.

[13]  J. Janicka,et al.  Flow field measurements of stable and locally extinguishing hydrocarbon-fuelled jet flames , 2003 .

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

[15]  S. Pope.  Lagrangian PDF Methods for Turbulent Flows , 1994 .

[16]  W. P. Jones,et al.  Large Eddy simulation of a turbulent non-premixed flame , 2001 .

[17]  W. P. Jones,et al.  Global reaction schemes for hydrocarbon combustion , 1988 .

[18]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

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

[20]  L. Valiño,et al.  A Field Monte Carlo Formulation for Calculating the Probability Density Function of a Single Scalar in a Turbulent Flow , 1998 .

[21]  R. Lindstedt,et al.  Joint scalar probability density function modeling of pollutant formation in piloted turbulent jet diffusion flames with comprehensive chemistry , 2000 .

[22]  R. Bilger,et al.  Conditional moment closure (CMC) predictions of a turbulent methane-air jet flame , 2001 .

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

[24]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[25]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[26]  T. Grundy,et al.  Progress in Astronautics and Aeronautics , 2001 .

[27]  Ulrich Schumann,et al.  Coherent structure of the convective boundary layer derived from large-eddy simulations , 1989, Journal of Fluid Mechanics.

[28]  P. Moin,et al.  A dynamic model for subgrid-scale variance and dissipation rate of a conserved scalar , 1998 .

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

[30]  H. Pitsch,et al.  Large-eddy simulation of a turbulent piloted methane/air diffusion flame (Sandia flame D) , 2000 .

[31]  Ugo Piomelli,et al.  Large-eddy simulation of rotating channel flows using a localized dynamic model , 1995 .

[32]  W. Jones,et al.  Large eddy simulation of a model gas turbine combustor , 2004 .

[33]  Yiliang Chen,et al.  PDF modelling of turbulent non-premixed combustion with detailed chemistry , 2004 .

[34]  C. Dopazo,et al.  Functional formulation of nonisothermal turbulent reactive flows , 1974 .

[35]  A. Klimenko,et al.  Conditional moment closure for turbulent combustion , 1999 .

[36]  N. Peters.  Laminar diffusion flamelet models in non-premixed turbulent combustion , 1984 .

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

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

[39]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[40]  N. Fueyo,et al.  Statistical Description of the Turbulent Mixing of Scalar Fields , 1997 .

[41]  B. E. Launder,et al.  Closure Strategies for Turbulent and Transitional Flows , 2002 .

[42]  Stephen B. Pope,et al.  A Monte Carlo Method for the PDF Equations of Turbulent Reactive Flow , 1981 .

[43]  Andrew W. Cook,et al.  A laminar flamelet approach to subgrid-scale chemistry in turbulent flows , 1997 .