Grid Resolution Effects on VSFMDF/LES - Springer Link

3 downloads 74 Views 731KB Size Report
Jul 9, 2010 - Submitted for the Special Issue Dedicated to S. B. Pope. Mehdi B. Nik·Server L. Yilmaz·. Mohammad Reza H. ..... VSFMDF methodology [30, 31] is significantly more powerful as it also accounts for the effects of SGS ...
Flow Turbulence Combust (2010) 85:677–688 DOI 10.1007/s10494-010-9272-5

Grid Resolution Effects on VSFMDF/LES Submitted for the Special Issue Dedicated to S. B. Pope Mehdi B. Nik · Server L. Yilmaz · Mohammad Reza H. Sheikhi · Peyman Givi

Received: 22 January 2010 / Accepted: 11 June 2010 / Published online: 9 July 2010 © Springer Science+Business Media B.V. 2010

Abstract The grid dependence of LES/VSFMDF is studied on a series of grids with progressively increased resolution reaching over 10 million grids for simulation of a turbulent piloted nonpremixed methane jet flame (Sandia D). In VSFMDF, the effects of the subgrid scale chemical reaction and convection appear in closed forms. The modeled transport equation for the VSFMDF is solved by a hybrid finite-difference/Monte Carlo scheme. A flamelet model is employed to relate the instantaneous composition to the mixture fraction. The simulated results are assessed via comparison with laboratory data. In addition, the dependence of predicted statistics on the grid size of the simulation is studied. The first order moments converge for the finest grid, but the higher order statistics including the PDFs are more sensitive to the grid resolution.

Keywords Large eddy simulation · Filter density function · Reacting flow

M. B. Nik (B) · P. Givi Department of Mechanical Engineering and Material Science, University of Pittsburgh, Pittsburgh, PA 15261, USA e-mail: [email protected] S. L. Yilmaz Center for Simulation and Modeling, University of Pittsburgh, Pittsburgh, PA 15213, USA M. R. H. Sheikhi Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA 02115, USA

678

Flow Turbulence Combust (2010) 85:677–688

1 Introduction The filtered density function (FDF) [1–3] is now regarded as one of the most effective means of conducting large eddy simulation (LES) in turbulent combustion. In its initial form, the marginal scalar FDF (SFDF) [4]; and its mass weighted scalar filtered mass density function (SFMDF) [5] provided the first demonstration of a “transported” FDF in reacting flows. The primary advantage of SFDF (SFMDF) is that it accounts for the subgrid scale (SGS) chemical reaction in a closed form. This closure is one of the reasons for SFMDF’s popularity and its widespread recent applications [6–27]. Inclusion of the “velocity” in the FDF accounts for the effects of “convection” in a closed form as well. This is demonstrated in the velocity-FDF (VFDF) [28], the joint velocity-scalar FDF (VSFDF) [29] and its density weighted VSFMDF [30]. Currently, this is the most comprehensive form of the FDF for hydrocarbon flame simulations [31]. The objective of the present work is to assess the grid dependence of the VSFMDF for LES of turbulent flames. For that, we consider the piloted non-premixed methane jet flame as studied in the experiments at the Sandia National Laboratories [32, 33] and at the TU-Darmstadt [34]. Flame D is considered in this work; this flame has been the subject of extensive previous LES via SFMDF by several investigators [12–14, 19–24]. With utilization of progressively larger number of grid points, the effects of resolution on various statistics of the flow variables as predicted by VSFMDF are investigated.

2 Simulations Sandia Flame D consists of a main jet with a mixture of 25% methane and 75% air by volume. The nozzle is placed in a coflow of air and the flame is stabilized by a substantial pilot. The Reynolds number for the main jet is Re = 22400 based on the nozzle diameter D = 7.2mm and the bulk jet velocity 49.6m/s. Details of the VSFMDF formulation for simulation of this flame is given in Ref. [31]. In this formulation, the primary transport variables are the density ρ(x, t), the velocity vector ui (x, t) (i = 1, 2, 3), the pressure p(x, t), the enthalpy h(x, t) and the species’ mass fractions Yα (x, t) (α = 1, 2, . . . , Ns ), where x ≡ xi , i = 1, 2, 3, t and Ns denote the space coordinate, the time and the total number of species involved in the chemical reaction. Large  +∞ eddy simulation involves the spatial filtering operation [35, 36] Q(x, t) = −∞ Q(x , t)G(x , x)dx , where G(x , x) denotes a filter function, and Q(x, t) is the filtered value of the transport variable Q(x, t). In variable-density flows it is convenient to use the Favre-filtered quantity Q(x, t) L = ρ Q / ρ . We consider a filter function that is spatially and tempoinvariant and localized, thus: G(x , x) ≡ G(x − x) with the properties G(x) ≥ 0, rally +∞ −∞ G(x)dx = 1. For closure of the VSFMDF transport equation, we consider the general diffusion process, [37] given by the system of stochastic differential equations (SDEs). In this context developed in Refs. [4, 28, 29, 38, 39] we utilize the simplified Langevin model (SLM) and the linear mean square estimation (LMSE) model [40].

Flow Turbulence Combust (2010) 85:677–688

679

 dX+ i

=

U i+ dt 

dU i+

=

+

2γ dWi , ρ

(1a)

 

    ∂ uj L ∂ uj L 2 1 ∂ ∂ ui  L 1 ∂ 2 ∂ γ − γ dt γ + ρ ∂ xj ρ ∂ xj ∂ xj ∂ xi 3 ρ ∂ xi ∂ xj

  1 ∂  p + Gij U +j − u j L dt + C0  dWi + − ρ ∂ xi



  dφα+ = −Cφ ω φα+ − φα  L dt + Sα (φ + ) dt,

2γ ∂ ui  L dWj, ρ ∂ xj

(1b) (1c)

where 

 1 3 Gij = −ω + C0 δij 2 4 k3/2  = C L

ω=

 k

(2)

1 k = τ L (ui , ui ) . 2

HereXi+ , U i+ , φα+ are probabilistic representations of position, velocity vector, and scalar variables, respectively. W terms denote the Wiener-Lévy processes [41, 42], γ denotes the molecular diffusion coefficients, and the chemical reaction source terms   Sα (φ(x, t)) are functions of compositional scalars (φ ≡ φ1 , φ2 , . . . , φ Ns +1 ). In Eq. 2, ω is the SGS mixing frequency,  is the dissipation rate, k is the SGS kinetic energy, and L is the LES filter size. The SGS correlations are dented by τ L (a, b) = ab L − a L b L , and the model parameters are the same as those suggested by Sheikhi et al. [30]: C0 = 2.1, Cφ = 1.0 and C = 1.0. No attempt is made to optimize the values of these parameters. Numerical solution of the modeled VSFMDF transport equation is obtained by a hybrid finite-difference (FD)/Monte Carlo (MC) procedure. The computational domain is discretized on equally spaced finite-difference grid points and the FMDF is represented by an ensemble of statistically identical MC particles which carry information pertaining to the velocity and the scalar values. Statistical information is obtained by considering an ensemble of N E computational particles residing within an ensemble domain of volume V E centered around each of the FD grid points. To reduce the computational cost, a procedure involving the use of non-uniform weights is also considered. This procedure allows a smaller number of particles in regions where a low degree of variability is expected. Conversely, in regions of high variability, a large number of particles is allowed. The sum of weights within the ensemble domain is related to filtered fluid density [43]. The FD solver is fourth-order accurate in space and second-order accurate in time [44]. To avoids the commutation error in the filtered and VSFMDF transport equations, all of the FD operations are conducted on fixed and equally spaced grid points. The transfer of information from the FD points to the MC particles is accomplished via a linear interpolation. The inverse transfer is accomplished via

680

Flow Turbulence Combust (2010) 85:677–688

Table 1 Grid Parameters and the corresponding filter size

Grid

Resolution (x,y,z)

Cells ×106

G1 G2 G3 G4

101 × 81 × 81 150 × 123 × 123 200 × 161 × 161 270 × 215 × 215

0.662 2.260 5.180 12.48

ensemble averaging. The FD transport equations include unclosed second order moments which are obtained from the MC. Simulations are conducted on a three-dimensional Cartesian mesh with uniform spacings in each of the three directions. With available resources, simulations are afforded within a domain spanning 20D × 16D × 16D in streamwise (x), and the two lateral (y, z) directions, respectively. The number of grid points for different test cases in the x, y and z directions are presented in Table 1. The filter size is set equal to L = 2( x y z)(1/3) where x, y and z denote the grid spacings in the corresponding directions. The MC particles are supplied in the inlet region and are free to move within the domain due to combined actions of convection and diffusion. The volume of the ensemble domain is V E = ( x y z) within which N E ≈ 40. The MC simulation results are monitored to ensure the particles fully encompass and extend well beyond regions of non-zero vorticity and reaction. For efficient parallel simulations, the domain is partitioned into equally sized partitions and each partition is assigned to a different processor at the onset of the simulations. This provides an effective parallelization and is relatively easy to implement. Each CPU is assigned relatively equal number of grid points and (approximately) equal number of particles. Simulations are conducted in conjunction with MPI and the PETSC [45–47] library. The flow variables at the inflow are set the same as those in the experiments, including the inlet profiles of the velocity and the mixture fraction. The inlet condition for the velocity is presented in Fig. 1. The flow is excited by superimposing oscillating axisymmetric perturbations at the inflow. The procedure is similar to

(a)

(b)

Fig. 1 Radial distribution of the mean and RMS values of the filtered axial velocity. Ucl denotes the mean axial velocity at the centerline at the inlet, the symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean axial velocity at x/D = 0.138, b RMS of the axial velocity at x/D = 0.138

Flow Turbulence Combust (2010) 85:677–688

(a)

681

(b)

Fig. 2 Radial distribution of the mean and RMS values of the filtered axial velocity. Ucl denotes the mean axial velocity at the centerline at the inlet, the symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean axial velocity at x/D = 15, b RMS of the axial velocity at x/D = 15

that in Ref. [48], but the amplitude of forcing is set in such a way to match the experimentally measured turbulent intensity of the streamwise velocity at the inlet. Standard characteristic boundary conditions [49] are implemented in all of the FD simulations. The methane-air reaction mechanism is taken into account via the “f lamelet” model. This model considers a laminar, one-dimensional counterflow (opposed jet) flame configuration [50]. The detailed kinetics mechanism of the Gas Research Institute (GRI2.11) [51] is employed to describe combustion. The flamelet table at strain rate of a = 100 1/s is used to relate the thermo-chemical variables to the mixture fraction. The predictive capability of VSFMDF is demonstrated by comparing the flow statistics with the experimental data [32–34]. These statistics are obtained by long-time averaging of the filtered field during 6 flow through times. The notations

(a)

(b)

Fig. 3 Radial distribution of the mean and RMS values of the filtered mixture fraction. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean mixture fraction at x/D = 15, b RMS of the mixture fraction at x/D = 15

682

Flow Turbulence Combust (2010) 85:677–688

(a)

(b)

Fig. 4 Radial distribution of the mean and RMS values of the filtered temperature values. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denote the RMS values. a Mean temperature (K) at x/D = 15, b RMS of the temperature (K) at x/D = 15

Q and RMS(Q) denote the time-averaged mean and root mean square values of the variable Q, respectively.

3 Results To assess the grid dependency of the VSFMDF simulations, we consider four grid configurations (G1 through G4) with a progressively increasing resolution from about 0.6 to 12.2 million grid points (Table 1). To compare with the experimental data, we consider the Reynolds-averaged moments of the thermo-chemical variables. In the presentation below, the overline denotes the Reynolds-averaged operator.

(a)

(b)

Fig. 5 Radial distribution of the mean and RMS values of filtered CH4 mass fractions. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean CH4 mass fraction at x/D = 15, b RMS of CH4 mass fraction at x/D = 15

Flow Turbulence Combust (2010) 85:677–688

(a)

683

(b)

Fig. 6 Radial distribution of the mean and RMS values of the filtered O2 mass fractions. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean O2 mass fraction at x/D = 15, b RMS of O2 mass fraction at x/D = 15

The root mean square (RMS) includes the contributions from both the resolved and the residual fields. The capability of the method in predicting the hydrodynamics field is demon strated by examining some of the (reported) radial (r = z2 + y2 ) distributions of the flow statistics in Fig. 2. The results at x/D = 15 are shown in this and all of the subsequent figures as the profiles at other streamwise location portray similar behaviors. The VSFMDF predicts the peak value of mean axial velocity profile and the spread of the jet reasonably well. The RMS values, however, is underpredicted. The radial distribution of the mixture fraction is also shown to compare well with data (Fig. 3). The mean fields of velocity and mixture fraction both converge to the results corresponding to the G4 grid. The convergence for the RMS field is slower, but the results for G4 compare best with experimental data.

(a)

(b)

Fig. 7 Radial distribution of the mean and RMS values of the filtered CO mass fractions. The symbols denote experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean CO mass fraction at x/D = 15, b RMS of CO mass fraction at x/D = 15

684

Flow Turbulence Combust (2010) 85:677–688

(a)

(b)

Fig. 8 Radial distribution of the mean and RMS values of the filtered CO2 mass fractions. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean CO2 mass fraction at x/D = 15, b RMS of CO2 mass fraction at x/D = 15

The radial distribution of the mean temperature and its corresponding RMS values are presented in Fig. 4. Similar to hydrodynamic quantities, the mean profiles show a relatively good agreement with measured data while the RMS values are underpredicted at the outerlayer. The statistics of the mass fractions (denoted by Y) of several of the species at different streamwise locations are compared with data in Figs. 5, 6, 7, 8 and 9. The mean profiles of the species show a close agreements with measurements. The RMS values show close agreements with measured data at the inner layer, but not as good at the outer layer. These disagreements can be attributed, in part, to the shortcoming of the flamelet model in relating the thermo-chemical variables to the mixture fraction. The mean fields of temperature and species mass fractions both converge to the results corresponding to the G4 grid.

(a)

(b)

Fig. 9 Radial distribution of the mean and RMS values of the filtered H2 O mass fractions. The symbols denote the experimental data. The line denotes the mean values and the thick solid line denotes the RMS values. a Mean H2 O mass fraction at x/D = 15, b RMS values of H2 O mass fraction at x/D = 15

Flow Turbulence Combust (2010) 85:677–688

685

In Fig. 10, the PDFs of the resolved mixture fraction as predicted by the VSFMDF are compared with those measured experimentally at several locations throughout the domain. In general, both the peak and the spreads of the PDFs are predicted

(a)

(b)

(c)

(d)

(e)

(f)

Fig. 10 PDF of the resolved filtered mixture fraction at x/D = 15 and different radial locations. The symbols and the thick lines denote the experimental data and LES predictions, respectively. a r = 2mm, b r = 4mm, c r = 6mm, d r = 8mm, e r = 10mm, f r = 12mm

686

Flow Turbulence Combust (2010) 85:677–688

well. The PDF’s converge to the results corresponding to the G4 grid; in which the best agreement with experimental data is observed.

4 Summary and Conclusions Since its original development, the SFMDF [4, 5] has experiences widespread applications for LES of a variety of reacting flows [6–8, 10, 11, 13–26, 29]. The extended VSFMDF methodology [30, 31] is significantly more powerful as it also accounts for the effects of SGS convection in an exact manner. In this work, we asses the grid dependence of VSFMDF in prediction of turbulent flyames. For that, we consider the piloted, nonpremixed, turbulent, methane jet flame (Sandia D). For this flame, the thermo-chemical variables are related to the mixture fraction via a flamelet table. A modeled transport equation for the mass weighted joint FDF of the velocity and the mixture fraction is considered [30]. This equation is solved by a hybrid finite-difference/Monte Carlo method. The predictive capability of the overall scheme is assessed by comparison of the ensemble (long time Reynolds-averaged) values of the thermo-chemical variables. The grid dependency of the predicted results is assessed via consideration of progressively larger number of grid points. In the case with the finest resolution, there are over 10 million grid points. The results of this grid resolution assessment indicate that the first order moments converges quickly; but the rate of convergence for the second order moments and the PDFs of the resolved field is slower. All of the mean quantities are, generally, in good comparison with experimental data. For the RMS values, the predicted values agree closely with the experimental data in the inner layer, but not as good in the outer layer. This discrepancy is attributed, in part, to the use of the flamelet model in relating thermo-chemical variables to the mixture fraction. Acknowledgements This work is sponsored by NASA under Grant NNX08AB36A. The computations support is provided by the National Science Foundation through TeraGrid resources at the Pittsburgh Supercomputing Center.

References 1. Pope, S.B.: Computations of turbulent combustion: progress and challenges. Proc. Combust. Inst. 23, 591–612 (1990) 2. Givi, P.: Model free simulations of turbulent reactive flows. Prog. Energy Combust. Sci. 15, 1–107 (1989) 3. Givi, P.: Filtered density function for subgrid scale modeling of turbulent combustion. AIAA J. 44(1), 16–23 (2006) 4. Colucci, P.J., Jaberi, F.A., Givi, P., Pope, S.B.: Filtered density function for large eddy simulation of turbulent reacting flows. Phys. Fluids 10(2), 499–515 (1998) 5. Jaberi, F.A., Colucci, P.J., James, S., Givi, P., Pope, S.B.: Filtered mass density function for large eddy simulation of turbulent reacting flows. J. Fluid Mech. 401, 85–121 (1999) 6. Garrick, S.C., Jaberi, F.A., Givi, P.: Large eddy simulation of scalar transport in a turbulent jet flow. In: Knight, D., Sakell, L. (eds.) Recent Advances in DNS and LES, Fluid Mechanics and its Applications, vol. 54, pp. 155–166. Kluwer Academic Publishers, The Netherlands (1999) 7. James, S., Jaberi, F.A.: Large scale simulations of two-dimensional nonpremixed methane jet flames. Combust. Flame 123, 465–487 (2000)

Flow Turbulence Combust (2010) 85:677–688

687

8. Zhou, X.Y., Pereira, J.C.F.: Large eddy simulation (2D) of a reacting plan mixing layer using filtered density function. Flow Turbul. Combust. 64, 279–300 (2000) 9. Zhou, X., Mahalingam, S.: A flame surface density based model for large eddy simulation of turbulent nonpremixed combustion. Phys. Fluids 14(11), 77–80 (2002) 10. Heinz, S.: On Fokker-Planck equations for turbulent reacting flows. Part 2. Filter density function for large eddy simulation. Flow Turbul. Combust. 70(1–4), 153–181 (2003) 11. Cha, C.M., Troullet, P.: A subgrid-scale mixing model for large-eddy simulations of turbulent reacting flows using the filtered density function. Phys. Fluids 15(6), 1496–1504 (2003) 12. Sheikhi, M.R.H., Drozda, T.G., Givi, P., Jaberi, F.A., Pope, S.B.: Large eddy simulation of a turbulent nonpremixed piloted methane jet flame (Sandia Flame D). Proc. Combust. Inst. 30, 549–556 (2005) 13. Raman, V., Pitsch, H., Fox, R.O.: Hybrid large-eddy simulation/lagrangian filtered density function approach for simulating turbulent combustion. Combust. Flame 143(1–2), 56–78 (2005) 14. Raman, V., Pitsch, H.: Large-Eddy simulation of a bluff-body-stabilized nonpremixed flame using a recursive filter-refinement procedure. Combust. Flame 142, 329–347 (2005) 15. van Vliet, E., Derksen, J.J., van den Akker, H.E.A.: Turbulent mixing in a tubular reactor: assessment of an FDF/LES approach. AIChE J. 51(3), 725–739 (2005) 16. Carrara, M.D., DesJardin, P.E.: A filtered mass density function approach to modeling separated two-phase flows using LES I: mathematical formulation. Int. J. Multiph. Flow 32, 365–384 (2006) 17. Mustata, R., Valiño, L., Jiménez, C., Jones, W.P., Bondi, S.: A probability density function Eulerian Monte Carlo field method for large eddy simulations: application to a turbulent piloted methane/air diffusion flame (Sandia D). Combust. Flame 145(1–2), 88–104 (2006) 18. Drozda, T.G., Sheikhi, M.R.H., Madnia, C.K., Givi, P.: Developments in formulation and application of the filtered density function. Flow Turbul. Combust. 78, 35–67 (2007) 19. Jones, W.P., Navarro-Martinez, S., Röhl, O.: Large eddy simulation of hydrogen auto-ignition with a probability density function method. Proc. Combust. Inst. 31, 1765–1771 (2007) 20. Jones, W.P., Navarro-Martinez, S.: Large eddy simulation of autoignition with a subgrid probability density function method. Combust. Flame 150, 170–187 (2007) 21. James, S., Zhu, J., Anand, M.S.: Large eddy simulations of turbulent flames using the filtered density function model. Proc. Combust. Inst. 31, 1737–1745 (2007) 22. Chen, J.Y.: A Eulerian PDF scheme for LES of nonpremixed turbulent combustion with secondorder accurate mixture fraction. Combust. Theor. Model. 11(5), 675–695 (2007) 23. McDermott, R., Pope, S.B.: A particle formulation for treating differential diffusion in filtered density function methods. J. Comput. Phys. 226(1), 947–993 (2007) 24. Raman, V., Pitsch, H.: A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry. Proc. Combust. Inst. 31, 1711–1719 (2007) 25. Afshari, A., Jaberi, F.A., Shih, T.I.P.: Large-eddy simulations of turbulent flows in an axisymmetric dump combustor. AIAA J. 46(7), 1576–1592 (2008) 26. Drozda, T.G., Wang, G., Sankaran, V., Mayo, J.R., Oefelein, J.C., Barlow, R.S.: Scalar filtered mass density functions in nonpremixed turbulent jet flames. Combust. Flame 155(1–2), 54–69 (2008) 27. Givi, P., Sheikhi, M.R.H., Drozda, T.G., Madnia, C.K.: Large scale simulation of turbulent combustion. Combust. Plasma Chem. 6(1), 1–9 (2008) 28. Gicquel, L.Y.M., Givi, P., Jaberi, F.A., Pope, S.B.: Velocity filtered density function for large eddy simulation of turbulent flows. Phys. Fluids 14(3), 1196–1213 (2002) 29. Sheikhi, M.R.H., Drozda, T.G., Givi, P., Pope, S.B.: Velocity-scalar filtered density function for large eddy simulation of turbulent flows. Phys. Fluids 15(8), 2321–2337 (2003) 30. Sheikhi, M.R.H., Givi, P., Pope, S.B.: Velocity-scalar filtered mass density function for large eddy simulation of turbulent reacting flows. Phys. Fluids 19(9), 095196 1–21 (2007) 31. Nik, M.B., Yilmaz, S.L., Givi, P., Sheikhi, M.R.H., Pope, S.B.: Simulation of sandia flame D using velocity-scalar filtered density function. AIAA J. 48(7), 1513–1522 (2010) 32. Nooren, P.A., Versiuis, M., Van der Meer, T.H., Barlow, R.S., Frank, J.H.: Raman-RayleighLIF measurements of temperature and species concentrations in the delft piloted turbulent jet diffusion flame. Appl. Phys. B71, 95–111 (2000) 33. Sandia National Laboratories, Combustion Research Facility website. http://www.ca.sandia. gov/crf/ (2003). Accessed 2 July 2010 34. Schneider, C., Dreizler, A., Janicka, J., Hassel, E.P.: Flow field measurements of stable and locally extinguishing hydrocarbon-fuelled jet flames. Combust. Flame 135(1), 185–190 (2003) 35. Pope, S.B.: Turbulent Flows. Cambridge University Press, Cambridge, UK (2000)

688

Flow Turbulence Combust (2010) 85:677–688

36. Sagaut, P.: Large Eddy Simulation for Incompressible Flows, 3rd edn. Springer-Verlag, New York, NY (2005) 37. Karlin, S., Taylor, H.M.: A Second Course in Stochastic Processes. Academic Press, New York, NY (1981) 38. Haworth, D.C., Pope, S.B.: A generalized Langevin model for turbulent flows. Phys. Fluids 29(2), 387–405 (1986) 39. Dreeben, T.D., Pope, S.B.: Probability density function and Reynolds-stress modeling of nearwall turbulent flows. Phys. Fluids 9(1), 154–163 (1997) 40. O’Brien, E.E.: The probability density function (PDF) approach to reacting turbulent flows. In: Libby, P.A., Williams, F.A. (eds.) Turbulent Reacting Flows, Topics in Applied Physics, vol. 44, pp. 185–218. Springer-Verlag, Heidelberg (1980) 41. Wax, N.: Selected Papers on Noise and Stochastic Processes. Dover (1954) 42. Gardiner, C.W.: Handbook of Stochastic Methods. Springer-Verlag, New York, NY (1990) 43. Madnia, C.K., Jaberi, F.A., Givi, P.: Large eddy simulation of heat and mass transport in turbulent flows. In: Minkowycz, W.J., Sparrow, E.M., Murthy, J.Y. (eds.) Handbook of Numerical Heat Transfer, 2nd edn, pp. 167–189. John Wiley & Sons, Inc., New York, NY (2006) 44. Kennedy, C.A., Carpenter, M.H.: Several new numerical methods for compressible shear-layer simulations. Appl. Numer. Math. 14, 397–433 (1994) 45. Balay, S., Buschelman, K., Gropp, W.D., Kaushik, D., Knepley, M.G., McInnes, L.C., et al.: PETSc Web page. http://www.mcs.anl.gov/petsc (2009). Accessed 2 July 2010 46. Balay, S., Buschelman, K., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M.G., et al.: PETSc users manual. Technical Report ANL-95/11-Revision 2.1.5, Argonne National Laboratory (2004) 47. Balay, S., Gropp, W.D., McInnes, L.C., Smith, B.F.: Efficient management of parallelism in object oriented numerical software libraries. In: Arge, E., Bruaset, A.M., Langtangen, H.P. (eds.) Modern Software Tools in Scientific Computing, pp. 163–202. Birkhäuser Press (1997) 48. Danaila, I., Boersma, B.J.: Direct numerical simulation of bifurcating jets. Phys. Fluids 12(5), 1255–1257 (2000) 49. Poinsot, T.J., Lele, S.K.: Boundary conditions for direct simulations of compressible viscous flows. J. Comp. Phys. 101, 104–129 (1992) 50. Peters, N.: Turbulent Combustion. Cambridge University Press, Cambridge, UK (2000) 51. Smith, G.P., Golden, D.M., Frenklach, M., Moriarty, N.W., Eiteneer, B., Goldenberg, M., et al.: http://www.me.berkeley.edu/gri_mech. Accessed 2 July 2010