Combined dimension reduction and tabulation strategy using ISAT–RCCE–GALI for the efficient implementation of combustion chemistry

Computations of turbulent combustion flows using detailed chemistry involving a large number of species and reactions are computationally prohibitive, even on a distributed computing system. Here, we present a new combined dimension reduction and tabulation methodology for the efficient implementation of combustion chemistry. In this study, the dimension reduction is performed using the rate controlled constrained-equilibrium (RCCE) method, and tabulation of the reduced space is performed using the in situ adaptive tabulation (ISAT) algorithm. The dimension reduction using RCCE is performed by specifying a set of represented (constrained) species, which in this study is selected using a new Greedy Algorithm with Local Improvement (GALI) (based on the greedy algorithm). This combined approach is found to be particularly fruitful in the probability density function (PDF) approach, wherein the chemical composition is represented by a large number of particles in the solution domain. In this work, the combined approach has been tested and compared to reduced and skeletal mechanisms using a partially-stirred reactor (PaSR) for premixed combustion of (i) methane/air (using the 31-species GRI-Mech 1.2 detailed mechanism and the 16-species ARM1 reduced mechanism) and (ii) ethylene/air (using the 111-species USC-Mech II detailed mechanism, a 38-species skeletal mechanism and a 24-species reduced mechanism). Results are presented to quantify the relative accuracy and efficiency of three different ways of representing the chemistry: (i) ISAT alone (with a detailed mechanism); (ii) ISAT (with a reduced or skeletal mechanism); and (iii) ISAT–RCCE with represented species selected using GALI. We show for methane/air: ISAT (with ARM1 reduced mechanism) incurs 6% error, while ISAT–RCCE incurs the same error using just 8 or more represented species, and less than 1% error using 11 or more represented species, with a twofold speedup relative to using ISAT alone with the GRI-Mech 1.2 detailed mechanism. And we show for ethylene/air: ISAT incurs 7% and 3% errors with the reduced and skeletal mechanisms, respectively, while ISAT–RCCE achieves the same levels of error 7% with just 18 and 3% with just 25 represented species, and also provides 15-fold speedup relative to using ISAT alone with the USC-Mech II detailed mechanism. With fewer species to track in the CFD code, this combined ISAT–RCCE–GALI reduction–tabulation algorithm provides an accurate and efficient way to represent combustion chemistry.

[1]  Stephen B. Pope,et al.  Implementation of combustion chemistry by in situ adaptive tabulation of rate-controlled constrained equilibrium manifolds , 2002 .

[2]  James C. Keck,et al.  Rate-controlled constrained-equilibrium method using constraint potentials , 1998 .

[3]  Stephen B. Pope,et al.  Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation , 1997 .

[4]  James C. Keck,et al.  Rate-controlled partial-equilibrium method for treating reacting gas mixtures , 1971 .

[5]  Stephen B. Pope,et al.  An improved algorithm for in situ adaptive tabulation , 2009, J. Comput. Phys..

[6]  S. Lam,et al.  The CSP method for simplifying kinetics , 1994 .

[7]  Samuel Paolucci,et al.  One-dimensional slow invariant manifolds for spatially homogenous reactive systems. , 2009, The Journal of chemical physics.

[8]  Michael Brunner,et al.  Mechanical and Aerospace Engineering , 2011 .

[9]  C. Westbrook,et al.  A comprehensive detailed chemical kinetic reaction mechanism for combustion of n-alkane hydrocarbons from n-octane to n-hexadecane , 2009 .

[10]  C. Law,et al.  A directed relation graph method for mechanism reduction , 2005 .

[11]  Harsha K. Chelliah,et al.  Explicit reduced reaction models for ignition, flame propagation, and extinction of C2H4/CH4/H2 and air systems , 2007 .

[12]  Harsha K. Chelliah,et al.  Skeletal reaction models based on principal component analysis: Application to ethylene–air ignition, propagation, and extinction phenomena , 2011 .

[13]  H. Pitsch,et al.  An efficient error-propagation-based reduction method for large chemical kinetic mechanisms , 2008 .

[14]  T. Turányi,et al.  Reduction of very large reaction mechanisms using methods based on simulation error minimization , 2009 .

[15]  Mitchell D. Smooke,et al.  Reduced Kinetic Mechanisms and Asymptotic Approximations for Methane-Air Flames: A Topical Volume , 1991 .

[16]  Robert W. Dibble,et al.  Pdf Modeling of Turbulent Nonpremixed Methane Jet Flames , 1989 .

[17]  Terese Løvås,et al.  Spectral uncertainty quantification, propagation and optimization of a detailed kinetic model for ethylene combustion , 2009 .

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

[19]  Zhuyin Ren,et al.  A greedy algorithm for species selection in dimension reduction of combustion chemistry , 2010 .

[20]  Stephen B. Pope.  Gibbs function continuation for the stable computation of chemical equilibrium , 2004 .

[21]  Thomas F. Edgar,et al.  Approximate nonlinear model predictive control with in situ adaptive tabulation , 2008, Comput. Chem. Eng..

[22]  S. Burgard.  ISR , 1999 .

[23]  Stephen B. Pope,et al.  An integrated PDF/neural network approach for simulating turbulent reacting systems , 1996 .

[24]  Hameed Metghalchi,et al.  Combustion modeling of mono-carbon fuels using the rate-controlled constrained-equilibrium method , 2009 .

[25]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[26]  Chung King Law,et al.  An augmented reduced mechanism for methane oxidation with comprehensive global parametric validation , 1998 .

[27]  Michael Frenklach,et al.  PRISM: piecewise reusable implementation of solution mapping. An economical strategy for chemical kinetics , 1998 .

[28]  Ulrich Maas,et al.  Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space , 1992 .

[29]  Stephen B. Pope,et al.  Probability density function calculations of local extinction and no production in piloted-jet turbulent methane/air flames , 2000 .

[30]  Sandip Mazumder,et al.  Modeling Full-Scale Monolithic Catalytic Converters: Challenges and Possible Solutions , 2007 .

[31]  Tamás Turányi,et al.  Parameterization of Reaction Mechanisms Using Orthonormal Polynomials , 1994, Comput. Chem..

[32]  J. Keck.  Rate-controlled constrained-equilibrium theory of chemical reactions in complex systems☆ , 1990 .

[33]  W. E. Stewart,et al.  Sensitivity analysis of initial value problems with mixed odes and algebraic equations , 1985 .