Research Papers: Gas Turbines: Combustion, Fuels, and Emissions

A Detailed Validation Study of Multi-Environment Eulerian Probability Density Function Transport Method for Modeling Turbulent Nonpremixed Combustion

[+] Author and Article Information
Rakesh Yadav

ANSYS Fluent India Pvt. Ltd.,
Pune 411057, India
e-mail: rakesh.yadav@ansys.com

Abhijit Kushari

Department of Aerospace Engineering,
I.I.T Kanpur 208016, India
e-mail: akushari@iitk.ac.in

Vinayak Eswaran

Department of Mechanical Engineering,
I.I.T Hyderabad 502025, India
e-mail: eswar@iith.ac.in

Atul K. Verma

ANSYS Fluent India Pvt. Ltd.,
Pune 411057, India
e-mail: atul.verma@ansys.com

1Corresponding author.

Contributed by the Combustion and Fuels Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received January 11, 2014; final manuscript received February 1, 2014; published online March 17, 2014. Editor: David Wisler.

J. Eng. Gas Turbines Power 136(8), 081506 (Mar 17, 2014) (11 pages) Paper No: GTP-14-1021; doi: 10.1115/1.4026861 History: Received January 11, 2014; Revised February 01, 2014

The current work involves the validation of presumed shape multi-environment Eulerian probability density function (PDF) transport method (MEPDF) using direct quadrature method of moments (DQMOM)-interaction by exchange with mean (IEM) approach for modeling turbulence chemistry interactions in nonpremixed combustion problems. The joint composition PDF is represented as a collection of finite number of Delta functions. The PDF shape is resolved by solving the governing transport equations for probability of occurrence of each environment and probability-weighted mass fraction of species and enthalpy in Eulerian frame for each environment. A generic implementation of the MEPDF approach is carried out for an arbitrary number of environments. In the current work, the MEPDF approach is used for a series of problems to validate each component of MEPDF approach in an isolated manner as well as their combined effect. First of all, a nonreactive turbulent mixing problem with two different Reynolds numbers (Re = 7000 and 11,900) is used for validation of the mixing and correction terms appear in the MEPDF approach. The second problem studied is a diffusion flame with infinitely fast chemistry having an analytical solution. The reaction component is validated by considering a 1D premixed laminar flame. In order to validate the combined effect of mixing and turbulence chemistry interactions, two different turbulent nonpremixed problems using global one-step chemistry are used. The first reactive problem used is H2 combustion (DLR Flame H3), while the second reactive validation case is a pilot-stabilized CH4 flame. The current predictions for all validation problems are compared with experimental data or published results. The study is further extended by modeling a turbulent nonpremixed H2 combustion using finite-rate chemistry effects and radiative heat transfer. The current model predictions for different flame lengths as well as minor species are compared with experimental data. The current model gave excellent predictions of minor species like OH. The differences in the current predictions with experimental data are discussed.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Magnussen, B. F., and Hjertager, B. H., 1976, “On the Mathematical Modeling of Turbulent Combustion With Special Emphasis on Soot Formation and Combustion,” Symp. (Int.) Combust., 16(1), pp. 719–729. [CrossRef]
Peters, N., 1984, “Laminar Diffusion Flamelet Models in Non-Premixed Turbulent Combustion,” Prog. Energy Combust. Sci., 10(3), pp. 319–339. [CrossRef]
Pope, S. B., 2000, Turbulent Flows, Cambridge University, Cambridge, UK.
Xu, J., and Pope, S. B., 2000, “PDF Calculations of Turbulent Non-Premixed Flames With Local Extinction,” Combust. Flame, 123(3), pp. 281–307. [CrossRef]
Liu, K., Pope, S. B., and Caughey, D. A., 2005, “Calculation of Bluff Body Stabilized Flames Using a Joint Probability Density Function Model With Detailed Chemistry,” Combust. Flame, 141(1–2), pp. 89–117. [CrossRef]
Merci, B., Naud, B., and Roekaerts, D., 2005, “Flow and Mixing Fields for Transported Scalar PDF Simulations of a Piloted Jet Diffusion Flame (‘Delft Flame III’),” Flow Turb. Combust., 74(3), pp. 239–272. [CrossRef]
Pope, S. B., 1997, “Computationally Efficient Implementation of Combustion Chemistry Using In Situ Adaptive Tabulation,” Combust. Theory Modeling, 1(1), pp. 41–63. [CrossRef]
Fox, R. O., 2003, Computational Models for Turbulent Reacting Flows, Cambridge University, Cambridge, UK.
McGraw, R., 1997, “Description of Aerosol Dynamics by the Quadrature Method of Moments,” Aerosol Sci. Technol., 27(2), pp. 255–265. [CrossRef]
Marchisio, D. L., Pikturna, J. T., Fox, R. O., Vigil, R. D., and Barresi, A. A., 2003, “Quadrature Method of Moments for Population-Balance Equations,” AIChE J., 49(5), pp. 1266–1276. [CrossRef]
Tang, Q., Zhou, W., Bockelie, M., and Fox, R. O., 2007, “Multi-Environment Probability Density Function Method for Modeling Turbulent Combustion Using Realistic Chemical Kinetics,” Combust. Theory Modeling, 11(6), pp. 889–907. [CrossRef]
Jaishree, J., and Haworth, D. C., 2009, “Comparison of Lagrangian and Eulerian PDF Methods in Simulations of Non-Premixed Turbulent Jet Flames With Strong Turbulence Chemistry Interactions,” 6th U.S National Combustion Meeting, Ann Arbor, MI, May 17–20, pp. 2517–2525.
Akroyd, J., Smith, A. J., McGlashan, L. R., and Kraft, M., 2010, “Numerical Investigation of DQMoM-IEM as a Turbulent Reaction Closure,” Chem. Eng. Sci., 65(6), pp. 1915–1924. [CrossRef]
Akroyd, J., Smith, A. J.McGlashan, L. R., and Kraft, M., 2010, “Comparison of the Stochastic Fields Method and DQMoM-IEM as Turbulent Reaction Closures,” Chem. Eng. Sci., 65(20), pp. 5429–5441. [CrossRef]
ANSYS Fluent R14 User Guide, 2011, Ansys, Inc., Canonsburg, PA, www.ansys.com
Clean Combustion Research Group, 2011, “Aerospace, Mechanical and Mechatronic Engineering,” The University of Sydney, Sydney, Australia, http://www.aeromech.usyd.edu.au/thermofluids/main_frame.htm
Kuo, K. K., 1986, Principles of Combustion, Wiley-Interscience, New York.
Kee, R. J., Grcar, J., Smooke, M., and Miller, J., 1985, “Fortran Program for Modeling Steady Laminar One-Dimensional Premixed Flames,” Sandia National Laboratories, Livermore, CA, Technical Report No. SAND85.
“DLR Portal,” 2014, German Aerospace Center (DLR), Cologne, Germany, http://www.dlr.de
Li, G., and Modest, M. F., 2003, “Importance of Turbulence Radiation Interactions in Turbulent Diffusion Jet Flames,” ASME J. Heat Transfer, 125(5), pp. 831–838. [CrossRef]
Barlow, R. S., and Frank, J. H., 1998, “Effects of Turbulence on Species Mass Fractions in Methane/Air Jet Flames,” Symp. (Int.) Combust., 27(1), pp. 1087–1095. [CrossRef]
Westbrook, C. K., and Dryer, F. L., 1981, “Simplified Reaction Mechanisms for the Oxidation of Hydrocarbon Fuels in Flames,” Combust. Sci. Technol., 27(1–2), pp. 31–43. [CrossRef]
Yadav, R., Verma, A. K., and Bessette, D., 2009, “A Hybrid Solution Methodology for Efficient Solving of Radiative Heat Transfer in Non Gray Medium Using WSSG Model,” ASME Heat Transfer Summer Conference Collocated With the InterPACK09 and 3rd Energy Sustainability Conferences (HT2009), San Francisco, CA, July 19–23, ASME Paper No. HT2009-88482, pp. 903–910. [CrossRef]
Schlatter, M., Ferreira, J. C., Flury, M., and Gass, J., 1996, “Analysis of Turbulence-Chemistry Interaction With Respect to NO Formation in Turbulent Nonpremixed Hydrogen-Air Flames,” Symp. (Int.) Combust., 26(2), pp. 2215–2222. [CrossRef]
Barlow, R. S., 2003, “Sandia H2/He Flame: Scalar Data,” Sandia National Laboratories, Livermore, CA, http://www.sandia.gov/TNF/DataArch/Sandia/SandH2He.pdf
Li, J., Zhao, Z., Kazakov, A., and Dryer, F. L., 2003, “An Updated Comprehensive Kinetics Model of Hydrogen Combustion,” Fall Technical Meeting of the Eastern States Section of the Combustion Institute, Penn State University, University Park, PA, October 26–29, pp. 169–171.


Grahic Jump Location
Fig. 1

Schematic of the computational domain having distinct fuel and air inlets

Grahic Jump Location
Fig. 2

Radial profiles of mean and rms fluctuations of mixture fraction at axial locations of (a) x/D = 20, (b) x/D = 40, (c) x/D = 60 for bluff body case with Re = 7000 (left) and Re = 11,900 (right)

Grahic Jump Location
Fig. 4

Comparison of velocity, temperature, density, and species mole fraction for 1D premixed methane flame (lines—current results, symbols—reference results [18])

Grahic Jump Location
Fig. 3

Profiles of mass fraction of fuel and oxidizer and the nondimensional temperature with mixture fraction for infinitely fast chemistry case (lines—present simulation, symbols—analytical results)

Grahic Jump Location
Fig. 9

Axial profiles at centerline for H2 jet flame

Grahic Jump Location
Fig. 8

Schematic of turbulent H2 jet flame geometry

Grahic Jump Location
Fig. 7

Axial profiles of temperature at centerline for piloted methane-air flame

Grahic Jump Location
Fig. 6

Axial profiles of mean mixture fraction, H2 mass fraction, H2O mass fraction, and temperature at centerline for DLR Flame H3 [19]

Grahic Jump Location
Fig. 5

Minor species profiles with different level of ISAT tolerances for 1D premixed reaction problem

Grahic Jump Location
Fig. 10

Radial profiles at axial locations of x/Lvis = 1/4 and 1/2 for H2 flame

Grahic Jump Location
Fig. 11

Initialization of the different flamelets for multiple unsteady radial profiles of mean and rms of mixture fraction for x/D (a) 20, (b) 40, and (c) 60 with different environments for the nonreactive mixing problem

Grahic Jump Location
Fig. 12

Axial profiles of species mass fraction and temperature at centerline for DLR Flame H3 with different number of environments



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In