Numerical implementation of mixing and molecular transport in LES/PDF studies of turbulent reacting flows

Probability Density Function (PDF) methods in combination with Large Eddy Simulations (LES) are a powerful tool for studying turbulent reacting flow problems and we are interested in the implementation of mixing and molecular transport in LES/PDF methods. The numerical methodology used for solution procedure is the hybrid particle/mesh method and a fractional step scheme is used to solve for transport, reaction and mixing sequentially. Mixing is modeled using the Interaction by Exchange with the Mean (IEM) model and the effects of molecular transport are incorporated as a mean drift term in the mixing step. This methodology avoids spurious production of scalar variance and also allows direct incorporation of differential diffusion effects. In this study, various numerical implementations of mixing and molecular transport are presented and evaluated, using the Method of Manufactured Solutions (MMS), for (1) accuracy, (2) detailed conservation, (3) realizability, and (4) stability. Moreover, the methodology is shown to be successful in capturing the effects of differential diffusion accurately with the additional property of ensuring realizability of species mass fractions. Finally and most importantly, we describe a new variance reduction technique by way of an implicit smoothing methodology. This smoothing scheme is shown to satisfy conservation, boundedness and regularity criteria. Moreover, for an appropriate choice of the smoothing length scale, significant improvements in accuracy can be achieved for an incremental increase in computational cost. Also, it is shown that with smoothing on a length scale greater than the grid size, the bias and statistical errors due to there being a finite number of particles in the Lagrangian Monte Carlo simulations scale as N"t"o"t^-^1 and N"t"o"t^-^1^/^2 respectively, where N"t"o"t is the total number of particles in the computational domain, whereas without smoothing these errors scale as N"p"c^-^1 and N"p"c^-^1^/^2, where N"p"c is the much smaller number of particles in a computational cell.

[1]  Dominique Pelletier,et al.  Verification of RANS solvers with manufactured solutions , 2007, Engineering with Computers.

[2]  Stephen B. Pope,et al.  The hybrid method for the PDF equations of turbulent reactive flows: consistency conditions and correction algorithms , 2001 .

[3]  César Dopazo,et al.  ''Relaxation'' of initial probability density functions in the turbulent convection of scalar fields , 1979 .

[4]  Stephen B. Pope,et al.  Molecular diffusion effects in LES of a piloted methane–air flame , 2011 .

[5]  James C. Sutherland,et al.  Quantification of differential diffusion in nonpremixed systems , 2004 .

[6]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[7]  Stephen B. Pope,et al.  A particle formulation for treating differential diffusion in filtered density function methods , 2006, J. Comput. Phys..

[8]  Stephen B. Pope,et al.  On the relationship between stochastic Lagrangian models of turbulence and second‐moment closures , 1994 .

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

[10]  Stephen B. Pope,et al.  Large eddy simulation/probability density function modeling of a turbulent CH4/H2/N2 jet flame , 2011 .

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

[12]  Heinz Pitsch,et al.  A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry , 2007 .

[13]  Christopher J. Roy,et al.  Verification of Euler/Navier–Stokes codes using the method of manufactured solutions , 2004 .

[14]  Haifeng Wang,et al.  Weak second-order splitting schemes for Lagrangian Monte Carlo particle methods for the composition PDF/FDF transport equations , 2010, J. Comput. Phys..

[15]  H. Pitsch.  LARGE-EDDY SIMULATION OF TURBULENT COMBUSTION , 2006 .

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

[17]  D. Haworth.  Progress in probability density function methods for turbulent reacting flows , 2010 .

[18]  Stephen B. Pope.  Self-conditioned fields for large-eddy simulations of turbulent flows , 2010, Journal of Fluid Mechanics.

[19]  Stephen B. Pope,et al.  Computations of turbulent combustion: Progress and challenges , 1991 .

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

[21]  J. Heimerl,et al.  A comparison of transport algorithms for premixed, laminar steady state flames , 1980 .

[22]  S. Pope,et al.  The effect of mixing models in PDF calculations of piloted jet flames , 2007 .

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

[24]  S. Pope.  An Improved Turbulent Mixing Model , 1982 .

[25]  Jinchao Xu,et al.  Assessment of Numerical Accuracy of PDF/Monte Carlo Methods for Turbulent Reacting Flows , 1999 .

[26]  Stephen B. Pope,et al.  The modeling of turbulent reactive flows based on multiple mapping conditioning , 2003 .

[27]  Alexander Y. Klimenko,et al.  A Generalised Multiple Mapping Conditioning Approach for Turbulent Combustion , 2009 .

[28]  Rs Cant,et al.  An Introduction to Turbulent Reacting Flows , 2007 .

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

[30]  Christopher J. Roy,et al.  Review of code and solution verification procedures for computational simulation , 2005 .

[31]  Stephen B. Pope,et al.  Diffusion Behind a Line Source in Grid Turbulence , 1985 .

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

[33]  Elaine S. Oran,et al.  Detailed modelling of combustion systems , 1981 .

[34]  P. Sagaut.  Large Eddy Simulation for Incompressible Flows , 2001 .

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

[36]  R. Curl.  Dispersed phase mixing: I. Theory and effects in simple reactors , 1963 .

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

[38]  Rodney O. Fox,et al.  On velocity‐conditioned scalar mixing in homogeneous turbulence , 1996 .

[39]  Stephen B. Pope,et al.  Paradigms in turbulent combustion research , 2005 .

[40]  Tianfeng Lu,et al.  Structure of a spatially developing turbulent lean methane–air Bunsen flame , 2007 .

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

[42]  Stephen B. Pope,et al.  The relationship between the probability approach and particle models for reaction in homogeneous turbulence , 1979 .

[43]  J Mathew,et al.  Filtered density function for large eddy simulation of turbulent reacting flows , 2008 .

[44]  S. Pope,et al.  Time-averaging strategies in the finite-volume/particle hybrid algorithm for the joint PDF equation of turbulent reactive flows , 2008 .