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

Large Eddy Simulation of an Enclosed Turbulent Reacting Methane Jet With the Tabulated Premixed Conditional Moment Closure Method

[+] Author and Article Information
Carlos Velez

Center for Advanced Turbomachinery and
Energy Research,
University of Central Florida,
Orlando, FL 32816
e-mail: velezcar@yahoo.com

Scott Martin

Eagle Flight Research Center,
Embry-Riddle Aeronautical University,
Daytona Beach, FL 32114
e-mail: martis38@erau.edu

Aleksander Jemcov

Institute for Flow Physics and Control,
University of Notre Dame,
South Bend, IN 46556
e-mail: ajemcov@nd.edu

Subith Vasu

Center for Advanced Turbomachinery and
Energy Research,
University of Central Florida,
Orlando, FL 32816
e-mail: subith@ucf.edu

Contributed by the Combustion and Fuels Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received July 21, 2015; final manuscript received January 28, 2016; published online April 12, 2016. Assoc. Editor: Song-Charng Kong.

J. Eng. Gas Turbines Power 138(10), 101501 (Apr 12, 2016) (9 pages) Paper No: GTP-15-1357; doi: 10.1115/1.4032846 History: Received July 21, 2015; Revised January 28, 2016

The tabulated premixed conditional moment closure (T-PCMC) method has been shown to provide the capability to model turbulent, premixed methane flames with detailed chemistry and reasonable runtimes in Reynolds-averaged Navier–Stokes (RANS) environment by Martin et al. (2013, “Modeling an Enclosed, Turbulent Reacting Methane Jet With the Premixed Conditional Moment Closure Method,” ASME Paper No. GT2013-95092). Here, the premixed conditional moment closure (PCMC) method is extended to large eddy simulation (LES). The new model is validated with the turbulent, enclosed reacting methane backward facing step data from El Banhawy et al. (1983, “Premixed, Turbulent Combustion of a Sudden-Expansion Flow,” Combust. Flame, 50, pp. 153–165). The experimental data have a rectangular test section at atmospheric pressure and temperature with an inlet velocity of 10.5 m/s and an equivalence ratio of 0.9 for two different step heights. Contours of major species, velocity, and temperature are provided. The T-PCMC model falls into the class of table lookup turbulent combustion models in which the combustion model is solved offline over a range of conditions and stored in a table that is accessed by the computational fluid dynamic (CFD) code using three controlling variables: the reaction progress variable (RPV), variance, and local scalar dissipation rate. The local scalar dissipation rate is used to account for the affects of the small-scale mixing on the reaction rates. A presumed shape beta function probability density function (PDF) is used to account for the effects of subgrid scale (SGS) turbulence on the reactions. SGS models are incorporated for the scalar dissipation and variance. The open source CFD code OpenFOAM is used with the compressible Smagorinsky LES model. Velocity, temperature, and major species are compared to the experimental data. Once validated, this low “runtime” CFD turbulent combustion model will have great utility for designing the next generation of lean premixed (LPM) gas turbine combustors.

FIGURES IN THIS ARTICLE
<>
Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.

References

Lieuwen, T. , Torres, H. , Johnson, C. , and Zinn, B. , 2001, “ A Mechanism of Combustion Instability in Lean Premixed Gas Turbine Combustors,” ASME J. Eng. Gas Turbines Power, 123(1), pp. 182–189. [CrossRef]
Amzin, S. , 2012, “ Conditional Moment Closure Method for Turbulent Premixed Flames,” Ph.D. thesis, Cambridge University, Cambridge, UK.
Poinsot, T. , and Veynante, D. , 2005, Theoretical and Numerical Combustion, Edwards, Philadelphia, PA.
Driscoll, J. , 2008, “ Turbulent Premixed Combustion: Flamelet Structure and Its Effect on Turbulent Burning Velocities,” Prog. Energy Combust. Sci., 34(1), pp. 91–134. [CrossRef]
Martin, S. , Kramlich, J. , Kosaly, G. , and Riley, J. , 2001, “ The Conditional Moment Closure Method For Premixed Turbulent Combustion,” Eastern States Section Fall Technical Meeting of the Combustion Institute, Hilton Head, SC, Dec. 3–5, pp. 245–248.
Klimenko, A. , 1990, “ Multicomponent Diffusion of Various Admixtures in Turbulent Flow,” Fluid Dyn., 25(3), pp. 327–334. [CrossRef]
Bilger, R. W. , 1991, “ Conditional Moment Methods for Turbulent Reacting Flow Using Crocco Variable Conditions,” The University of Sydney, Sydney, Australia, Charles Kolling Report No. TN F-99.
Klimenko, A. Y. , 2001, “ On the Relation Between the Conditional Moment Closure and Unsteady Flamelets,” Combust. Theory Modell., 5(3), pp. 275–294. [CrossRef]
Bolla, M. , Wright, Y. , Boulouchos, K. , Borghesi, G. , and Mastorakos, E. , 2013, “ Soot Formation Modeling of n-Heptane Sprays Under Diesel Engine Conditions Using the Conditional Moment Closure Approach,” Combust. Sci. Technol., 185(5), pp. 766–793. [CrossRef]
Garmory, A. , and Mastorakos, E. , 2011, “ Capturing Localised Extinction in Sandia Flame F With LES–CMC,” Proc. Combust. Inst., 33(1), pp. 1673–1680. [CrossRef]
Bottone, F. , Kronenburg, A. , Gosman, D. , and Marquis, A. , 2012, “ The Numerical Simulation of Diesel Spray Combustion With LES-CMC,” Flow, Turbul. Combust., 89(4), pp. 651–673. [CrossRef]
Navarro-Martinez, S. , and Kronenburg, A. , 2009, “ LES–CMC Simulations of a Lifted Methane Flame,” Proc. Combust. Inst., 32(1), pp. 1509–1516. [CrossRef]
Ma, M. , and Devaud, C. B. , 2014, “ A Conditional Moment Closure (CMC) Formulation Including Differential Diffusion Applied to a Non-Premixed Hydrogen–Air Flame,” Combust. Flame, 162(1), pp. 144–158. [CrossRef]
Thornber, B. , Bilger, R. W. , Masri, A. R. , and Hawkes, E. R. , 2011, “ An Algorithm for LES of Premixed Compressible Flows Using the Conditional Moment Closure Model,” J. Comput. Phys., 230(20), pp. 7687–7705. [CrossRef]
Bilger, R. W. , 1993, “ Conditional Moment Closure Modelling and Advanced Laser Measurements,” Turbulence and Molecular Processes in Combustion, Vol. 138, T. Takeno, ed., Elsevier, Amsterdam, pp. 267–285.
Martin, S. , Kramlich, J. C. , Kosaly, G. , and Riley, J. , 2003, “ The Premixed Conditional Moment Closure Method,” ASME J. Eng. Gas Turbines Power, 125(4), pp. 895–900. [CrossRef]
Martin, S. M. , 2003, “ The Conditional Moment Closure Method for Modeling Lean Premixed Turbulent Combustion,” Ph.D. thesis, The University of Washington, Seattle, WA.
Martin, S. , Jemcov, A. , and de Ruijter, B. , 2013, “ Modeling an Enclosed, Turbulent Reacting Methane Jet With the Premixed Conditional Moment Closure Method,” ASME Paper No. GT2013-95092.
Smith, N. S. , 1994, “ Development of the Conditional Moment Closure Method for Modelling Turbulent Combustion,” Ph.D. thesis, The University of Sydney, Sydney, Australia.
Devaud, C. B. , and Bray, N. C. , 2003, “ Assessment of the Applicability of Conditional Moment Closure to a Lifted Turbulent Flame: First Order Model,” Combust. Flame, 132(1), pp. 102–114. [CrossRef]
Klimenko, A. Y. , and Bilger, R. W. , 1999, “ Conditional Moment Closure for Turbulent Combustion,” Prog. Energy Combust. Sci., 25(6), pp. 595–687. [CrossRef]
Jimenex, C. , Ducros, F. , and Bedat, B. , 2001, “ Subgrid Scale Variance and Dissipation of a Scalar Field in Large Eddy Simulations,” Phys. Fluids, 13(6), p. 1748. [CrossRef]
Donnini, A. , Martin, S. , Bastiaans, R. , van Oijen, J. , and de Goey, L. , 2013, “ Numerical Simulations of a Premixed Turbulent Confined Jet Flame Using the Flamelet Generated Manifold Approach With Heat Loss Inclusion,” ASME Paper No. GT2013-94363.
El Banhawy, Y. , Sivasegaram, S. , and Whitelaw, J. , 1983, “ Premixed, Turbulent Combustion of a Sudden-Expansion Flow,” Combust. Flame, 50, pp. 153–165. [CrossRef]
Smith, G. P. , Golder, D. M. , Frenklach, M. , Moriarty, N. W. , Eiteneer, B. , Goldenberg, B. , Bowmann, C. T. , Hanson, R. K. , Song, S. , Gardiner, W. C., Jr. , Lissianski, V. V. , and Qin, Z. , 2002, “ GRI-Mech 3.0,” University of California, Berkeley, Berkeley, CA, http://www.me.berkeley.edu/gri_mech/
Pitsch, H. , 2002, “ Improved Pollutant Predictions in Large-Eddy Simulations of Turbulent Non-Premixed Combustion by Considering Scalar Dissipation Rate Fluctuations,” Proc. Combust. Inst., 29(2), pp. 1971–1978. [CrossRef]
Tanahashi, M. , Nada, Y. , Ito, Y. , and Miauchi, T. , 2002, “ Local Flame Structure in the Well-Stirred Reactor Regime,” Proc. Combust. Inst., 29(2), pp. 2041–2049. [CrossRef]
Fureby, C. , Grinstein, F. , Menon, S. , and Weller, H. , 2001, “ A Comparitive Study of LES Combustion Models for Premixed Turbulent Combustion,” Second Symposium on Turbulence and Shear Flow Phenomena (TSFP2), Stockholm, Sweden, June 26–29, Vol. 3, pp. 309–314.
Weller, H. , Tabor, G. , Gosman, D. , and Fureby, C. , 1998, “ Application of a Flame Wrinkling LES Combustion Model to a Turbulent Mixing Layer,” Combust. Symp., 27(1), pp. 899–907. [CrossRef]
Park, N. , Kobayashi, T. , and Taniguchi, N. , 2001, “ Application of Dynamic Subgrid G-Equation Model to LES of Turbulent Premixed Flame Over Backward Facing Step,” Turbulence and Shear Flow Phenomena (TSFP2), Stockholm, Sweden, June 26–29, Vol. 3, pp. 291–296.
Jasak, H. , Jemcov, A. , and Tukovic, Z. , 2007, “ OpenFOAM: A C++ Library for Complex Physics Simulation,” International Workshop on Coupled Methods in Numerical Dynamics, (CMND 2007), Dubrovnik, Croatia, Sept. 19–21.
Weller, H. , Tabor, G. , Jasak, H. , and Fureby, C. , 1998, “ A Tensorial Approach to Computational Continuum Mechanics Using Object Oriented Techniques,” Comput. Phys., 12(6), pp. 620–631. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Flow chart of the T-PCMC CFD model

Grahic Jump Location
Fig. 2

Density as a function of RPV and variance for a scalar dissipation rate of 200

Grahic Jump Location
Fig. 3

RPV source term as a function of RPV and scalar dissipation rate at a variance equal to 1/3Cvar,Max

Grahic Jump Location
Fig. 4

Borghi diagram: computed from CFD in flame zone (c = 0.05–0.95)

Grahic Jump Location
Fig. 5

Experimental (bottom) and T-PCMC (top) temperature contours

Grahic Jump Location
Fig. 6

Compares the contours of the measured and predicted mean axial velocity

Grahic Jump Location
Fig. 7

T-PCMC (top) and experimental (bottom) RMS of axial velocity contours

Grahic Jump Location
Fig. 8

T-PCMC (top) and experimental (bottom) mean CO2 % volume concentrations

Grahic Jump Location
Fig. 9

Experimental (top) and T-PCMC (bottom) mean CO% volume concentrations

Grahic Jump Location
Fig. 10

Experimental (top) and T-PCMC (bottom) mean UHC% volume concentrations

Grahic Jump Location
Fig. 11

Experimental versus T-PCMC profiles for mean CO2 % volume concentrations. Lines—CFD and symbols—experiment.

Grahic Jump Location
Fig. 12

Experimental versus T-PCMC profiles for mean UHC% volume concentrations. Lines—CFD and symbols—experiment.

Tables

Errata

Discussions

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