Eulerian transported probability density function sub-filter model for large-eddy simulations of turbulent combustion

Reactive flow simulations using large-eddy simulations (LES) require modelling of sub-filter fluctuations. Although conserved scalars like mixture fraction can be represented using a beta-function, the reactive scalar probability density function (PDF) does not follow an universal shape. A one-point one-time joint composition PDF transport equation can be used to describe the evolution of the scalar PDF. The high-dimensional nature of this PDF transport equation requires the use of a statistical ensemble of notional particles and is directly coupled to the LES flow solver. However, the large grid sizes used in LES simulations will make such Lagrangian simulations computationally intractable. Here we propose the use of a Eulerian version of the transported-PDF scheme for simulating turbulent reactive flows. The direct quadrature method of moments (DQMOM) uses scalar-type equations with appropriate source terms to evolve the sub-filter PDF in terms of a finite number of delta-functions. Each delta-peak is characterized by a location and weight that are obtained from individual transport equations. To illustrate the feasibility of the scheme, we compare the model against a particle-based Lagrangian scheme and a presumed PDF model for the evolution of the mixture fraction PDF. All these models are applied to an experimental bluff-body flame and the simulated scalar and flow fields are compared with experimental data. The DQMOM model results show good agreement with the experimental data as well as the other sub-filter models used.

[1]  Bassam B. Dally,et al.  Instantaneous and Mean Compositional Structure of Bluff-Body Stabilized Nonpremixed Flames , 1998 .

[2]  S. Pope,et al.  A Consistent Hybrid Finite-Volume / Particle Method for the PDF Equations of Turbulent Reactive Flows , 2001 .

[3]  Rodney O. Fox,et al.  A multienvironment conditional probability density function model for turbulent reacting flows , 2004 .

[4]  Bassam B. Dally,et al.  Flow and mixing fields of turbulent bluff-body jets and flames , 1998 .

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

[6]  J. Riley,et al.  A subgrid model for equilibrium chemistry in turbulent flows , 1994 .

[7]  Laurent Falk,et al.  A generalized mixing model for initial contacting of reactive fluids , 1994 .

[8]  Rodney O. Fox,et al.  Comparison of micromixing models for CFD simulation of nanoparticle formation , 2004 .

[9]  C. Pierce,et al.  Progress-variable approach for large-eddy simulation of turbulent combustion , 2001 .

[10]  Juan Pedro Mellado,et al.  Reconstruction subgrid models for nonpremixed combustion , 2003 .

[11]  Daniele Marchisio,et al.  Solution of population balance equations using the direct quadrature method of moments , 2005 .

[12]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[13]  Stephen B. Pope,et al.  PDF modeling of a bluff-body stabilized turbulent flame , 2003 .

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

[15]  Stephen B. Pope,et al.  PDF calculations of piloted turbulent nonpremixed flames of methane , 1990 .

[16]  Heinz Pitsch,et al.  Large-eddy simulation of a bluff-body-stabilized non-premixed flame using a recursive filter-refinement procedure , 2005 .

[17]  R. Fox.  Computational Models for Turbulent Reacting Flows , 2003 .

[18]  Paul E. Dimotakis,et al.  Mixing and combustion with low heat release in a turbulent shear layer , 1984, Journal of Fluid Mechanics.

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

[20]  P. Moin,et al.  Progress-variable approach for large-eddy simulation of non-premixed turbulent combustion , 2004, Journal of Fluid Mechanics.

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

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

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

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

[25]  S. Pope,et al.  A Hybrid Algorithm for the Joint PDF Equation of Turbulent Reactive Flows , 2001 .

[26]  Kang Y. Huh,et al.  Use of the conditional moment closure model to predict NO formation in a turbulent CH4/H2 flame over a bluff-body , 2002 .

[27]  Robert W. Bilger,et al.  Conditional moment closure for turbulent reacting flow , 1993 .

[28]  Parviz Moin,et al.  An evaluation of the assumed beta probability density function subgrid-scale model for large eddy simulation of nonpremixed, turbulent combustion with heat release , 2000 .