×

zbMATH — the first resource for mathematics

A quasi-steady-state solver for the stiff ordinary differential equations of reaction kinetics. (English) Zbl 0967.65083
The presented quasi-steady-state method integrates stiff ordinary differential equations arising from reaction kinetics. This predictor-corrector method is A-stable for linear equations and second-order accurate. It works well for problems typical for hydrocarbon combustion. The algorithm is described, and error analysis and linear stability analysis are included. The algorithm is applied to test problems, and the results are compared to those of the stiff integrator CHEMEQ. The method proposed in this paper is more stable and more accurate than CHEMEQ.

MSC:
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
80A32 Chemically reacting flows
65L20 Stability and convergence of numerical methods for ordinary differential equations
65L05 Numerical methods for initial value problems involving ordinary differential equations
80A25 Combustion
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Hoffman, J.D., Numerical methods for engineers and scientists, (1992) · Zbl 0823.65006
[2] Lambert, J.D., Numerical methods for ordinary differential systems; the initial value problem, (1991) · Zbl 0745.65049
[3] Oran, E.S.; Boris, J.P., Numerical simulation of reactive flow, (1987) · Zbl 0762.76098
[4] Young, T.R.; Boris, J.P., A numerical technique for solving ordinary differential equations associated with the chemical kinetics of reactive-flow problems, J. phys. chem., 81, 2424, (1977)
[5] Young, T.R., CHEMEQ—A subroutine for solving stiff ordinary differential equations, (1980)
[6] Oran, E.S.; Young, T.R.; Boris, J.P., Application of time-dependent numerical methods to the description of reactive shocks, Proceedings of 17th international symp. on combustion, 43, (1979)
[7] Kailasanath, K.; Oran, E.S.; Boris, J.P., A theoretical study of the ignition of premixed gasses, Combust. flame, 47, 173, (1982)
[8] Oran, E.S.; Young, T.R.; Boris, J.P.; Cohen, A., Weak and strong ignition. I. numerical simulations of shock tube experiments, Combust. flame, 48, 135, (1982)
[9] Oran, E.S.; Boris, J.P., Weak and strong ignition. II. sensitivity of the hydrogen – oxygen system, Combust. flame, 48, 149, (1982)
[10] Kailasanath, K.; Oran, E.S., Ignition of flamelets behind incident shock waves and the transition to detonation, Combust. sci. technol., 34, 345, (1983)
[11] Weber, J.W.; Oran, E.S.; Anderson, J.D.; Patnaik, G., Load balancing and performance issues for the data parallel simulation of stiff chemical nonequilibrium flow, Aiaa j., 35, 147, (1998)
[12] Oran, E.S.; Weber, J.W.; Stefaniw, E.I.; Lefebvre, M.H.; Anderson, J.D., A numerical study of a two-dimensional H_{2}-O2-ar detonation using a detailed chemical reaction model, Combust. flame, 113, 147, (1998)
[13] Doschek, G.A.; Boris, J.P.; Cheng, C.-C.; Mariska, J.T.; Oran, E.S., Numerical simulation of cooling coronal flare plasma, Astrophys. J., 258, 373, (1982)
[14] Cheng, C.-C.; Oran, E.S.; Doschek, G.A.; Boris, J.P.; Mariska, J.T., Numerical simulations of loops heated to solar flare temperatures. I. gasdynamics, Astrophys. J., 265, 1090, (1983)
[15] Cheng, C.-C.; Oran, E.S.; Doschek, G.A.; Boris, J.P.; Mariska, J.T., Numerical simulations of loops heated to solar flare temperatures. II. X-ray and UV spectroscopy, Astrophys. J., 265, 1103, (1983)
[16] Avro, C.J., CHEMSODE: A stiff ODE solver for the equations of chemical kinetics, Comput. phys. commun., 97, 304, (1996) · Zbl 0926.65072
[17] Lorenzini, R.; Passoni, L., Test of numerical methods for the integration of kinetic equations in tropospheric chemistry, Comput. phys. commun., 117, 241, (1999) · Zbl 0980.65068
[18] G. Patnaik, and, K. Kailasanath, Numerical predictions of the cell-split limit in lean premixed hydrogen – air flames in microgravity, presented at the Joint Meeting of the, United States Sections of the Combustion Institute, Washington, DC, March 1999.
[19] K. Kailasanath, G. Patnaik, and, C. Li, Computational studies of pulse detonation engines: A status report, AIAA Paper 99-2634, presented at the 35th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Los Angeles, 1999.
[20] Mott, D.R., New quasi-steady-state and partial-equilibrium methods for integrating chemically reacting systems., (1999)
[21] D. R. Mott, and, E. S. Oran, CHEMEQ2: A Solver for the Stiff Differential Equations of Reaction Kinetics, NRL Memorandum Report, Naval Research Laboratory, Washington, DC, in preparation. · Zbl 0967.65083
[22] Verwer, J.G.; van Loon, M., An evaluation of explicit psuedo-steady-state approximation schemes for stiff ODE systems from chemical kinetics, J. comput. phys., 113, 347, (1994) · Zbl 0810.65068
[23] Verwer, J.G.; Simpson, D., Explicit methods for stiff ODEs from atmospheric chemistry, Appl. numer. math., 18, 413, (1995) · Zbl 0839.65078
[24] Jay, L.O.; Sandu, A.; Porta, F.A.; Carmichael, G.R., Improved quasi-steady-state-approximation methods for atmospheric chemistry integration, SIAM J. sci. comput., 18, 182, (1997) · Zbl 0869.65045
[25] Radhakrishnan, K.; Pratt, D.T., Fast algorithm for calculating chemical kinetics in turbulent reacting flow, Combust. sci. technol., 58, 155, (1988)
[26] Hindmarsh, A.C., LSODE and LSODEI, two new initial value ordinary differential equation solvers, ACM SIGNUM newslett., 15, 10, (1980)
[27] A. C. Hindmarsh, ODEPACK: A systemized collection of ODE solvers, in Scientific Computing, edited by R. S. Steplemanet al.North-Holland, Amsterdam, 1983, pp. 55-64. [Also UCRL-880-07, Lawrence Livermore Laboratory, Livermore, CA, 1993]
[28] Gear, C.W., Numerical initial value problems in ordinary differential equations, (1971) · Zbl 0217.21701
[29] Radhakrishnan, K., New integration techniques for chemical kinetic rate equations. I. efficiency comparison, Combust. sci. technol., 46, 59, (1986)
[30] Williams, F.A., Combustion theory, the fundamental theory of chemically reacting systems, (1985)
[31] Lam, S.H.; Goussis, D.A., The CSP method for simplifying kinetics, Int. J. chem. kinet., 26, 461, (1994)
[32] Frenklach, M.; Wang, H.; Rabinowitz, M.J., Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method—combustion of methane, Prog. energy combust. sci., 18, 47, (1992)
[33] U. Uphoff, D. Hänel, and, P. Roth, A grid refinement study for detonation simulation with detailed chemistry, in, Proc. of 6th Int. Conf. Numerical Combustion, New Orleans, March 4-6, 1996.
[34] P. Deuflhard, U. Nowak, and, U. Poehle, Scientific Software Group, Konrad-Zuse-Zentrum fuer Informationstechnik, Berlin, 1989.
[35] Deuflhard, P., A semi-implicit midpoint rule for stiff systems of ordinary differential equations, Numer. math., 41, 373, (1983) · Zbl 0522.65050
[36] Deuflhard, P., Order and stepsize control in extrapolation methods, Numer. math., 41, 399, (1983) · Zbl 0543.65049
[37] P. Deuflhard, Uniqueness Theorems for Stiff ODE Initial Value Problems, Preprint SC-87-3, Konrad-Zuse-Zentrum fuer Informationstechnik, Berlin, 1987.
[38] Patnaik, G.; Laskey, K.J.; Kailasanath, K.; Oran, E.S.; Brun, T.A., FLIC—A detailed, two-dimensional flame model, (1989)
[39] Boris, J.P.; Landsberg, Alexandra M.; Oran, Elaine S.; Gardner, John H., LCPFCT—A flux-corrected transport algorithm for solving generalized continuity equations, (1993)
[40] Vandevender, W.H.; Haskell, K.H., The SLATEC mathematical subroutine library, SIGNUM newslett., 17, 16, (1982)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.