CHEMEQ2: A Solver for the Stiff Ordinary Differential Equations of Chemical Kinetics

Abstract : This report describes and documents the subroutine CHEMEQ2, used to integrate stiff ordinary differential equations arising from reaction kinetics. This is a second generation improvement of CHEMEQ using a new quasi-steady-state predictor-corrector method that is A-stable for linear equations and second-order accurate. A single integration method can now be used for all species, regardless of the timescales of the individual equations. Start-up costs and memory requirements are low, so CHEMEQ2 is ideal for multi-dimensional reacting-flow simulations which require the solution of a process-split, initial-value problem in every computational cell for every global timestep. The algorithm, its implementation, the FORTRAN code, the internal variables and the argument lists are presented, along with several test problem results.

[1]  Michael Frenklach,et al.  Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method—combustion of methane , 1992 .

[2]  Colin J. Aro.  CHEMSODE: a stiff ODE solver for the equations of chemical kinetics , 1996 .

[3]  W. H. Vandevender,et al.  The SLATEC mathematical subroutine library , 1982, SGNM.

[4]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[5]  J. G. Verwer,et al.  Explicit method for stiff ODEs from atmospheric chemistry , 1995 .

[6]  Elaine S. Oran,et al.  Numerical simulations of loops heated to solar flare temperatures. I - Gasdynamics. II - X-ray and UV spectroscopy , 1983 .

[7]  Elaine S. Oran,et al.  Load balancing and performance issues for the data parallel simulation of stiff chemical nonequilibrium flows , 1995 .

[8]  Gopal Patnaik,et al.  COMPUTATIONAL STUDIES OF PULSE DETONATION ENGINES: A STATUS REPORT , 1999 .

[9]  Elaine S. Oran,et al.  Application of time-dependent numerical methods to the description of reactive shocks , 1979 .

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

[11]  J. Boris,et al.  A numerical technique for solving stiff ordinary differential equations associated with the chemical kinetics of reactive-flow problems , 1977 .

[12]  P. Deuflhard,et al.  A semi-implicit mid-point rule for stiff systems of ordinary differential equations , 1983 .

[13]  Elaine S. Oran,et al.  Weak and strong ignition. I. Numerical simulations of shock tube experiments , 1982 .

[14]  Elaine S. Oran,et al.  A theoretical study of the ignition of premixed gases , 1982 .

[15]  B. V. Leer,et al.  A quasi-steady state solver for the stiff ordinary differential equations of reaction kinetics , 2000 .

[16]  Elaine S. Oran,et al.  Numerical simulations of loops heated to solar flare temperatures. I , 1983 .

[17]  David T. Pratt,et al.  Fast Algorithm for Calculating Chemical Kinetics in Turbulent Reacting Flow , 1988 .

[18]  T. Young,et al.  CHEMEQ - A Subroutine for Solving Stiff Ordinary Differential Equations , 1980 .

[19]  Peter Deuflhard,et al.  Uniqueness Theorems for Stiff ODE Initial Value Problems. , 1987 .

[20]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[21]  Elaine S. Oran,et al.  A Numerical Study of a Two-Dimensional H2-O2-Ar Detonation Using a Detailed Chemical Reaction Model , 1998 .

[22]  Krishnan Radhakrishnan,et al.  New Integration Techniques for Chemical Kinetic Rate Equations. I. Efficiency Comparison , 1986 .

[23]  Colin J. Aro.  CHEMSODE: A stiff ODE solver for the equations of chemical kinetics , 1995 .

[24]  Jan Verwer,et al.  An evaluation of explicit pseudo-steady-state approximation schemes for stiff ODE systems from chemical kinetics , 1993 .

[25]  Adrian Sandu,et al.  Improved Quasi-Steady-State-Approximation Methods for Atmospheric Chemistry Integration , 1997, SIAM J. Sci. Comput..

[26]  C HindmarshAlan.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980 .

[27]  R. Lorenzini,et al.  Test of numerical methods for the integration of kinetic equations in tropospheric chemistry , 1999 .

[28]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[29]  D. Mott,et al.  New Quasi-Steady-State And Partial-Equilibrium Methods For Integrating Chemically Reacting Systems , 1999 .

[30]  Elaine S. Oran,et al.  Ignition of flamelets behind incident shock waves and the transition to detonation. Memorandum report , 1983 .

[31]  Elaine S. Oran,et al.  LCPFCT-A Flux-Corrected Transport Algorithm for Solving Generalized Continuity Equations , 1993 .

[32]  P. Deuflhard.  Order and stepsize control in extrapolation methods , 1983 .