Abstract
Simulation of large networks of chemical reactions via the numerical integration of large systems of ordinary differential equations is of growing importance in real-world problems. We propose an attractive novel numerical integration method, that is largely independent from ill-conditioning and is suitable for any nonlinear problem; moreover, the method, being exact for linear problems, is especially precise for quasi-linear problems, the most frequent kind in the real world. The method is based on a new approach to the computation of a matrix exponential, includes an automatic correction of rounding errors, is not too expensive computationally, and lends itself to a short and robust software implementation that can be easily inserted in large simulation packages. A preliminary numerical verification has been performed, with encouraging results, on two sample problems. The full source listing (in standard C language) of an academic version of the algorithm is freely available on request (e-mail address: Valerio.Parisi@roma2.infn.it), together with a very simple but very stiff chemical problem.
Similar content being viewed by others
References
J.D. Lambert and S.T. Sigurdsson, SIAM J. Numer. Anal. 9 (1972) 715.
G.R. Gavalas, Nonlinear Differential Equations of Chemically Reacting Systems (Springer-Verlag, Berlin, 1968).
E.H. Kerner, J. Math. Phys. 22 (1981) 1366.
M. Hochbruck, C. Lubich and H. Selhofer, SIAM J. Sci. Comput. 19 (1998) 1552.
C. Moler and C. van Loan, SIAM Review 20 (1978) 801.
G.H. Golub and C.F. van Loan, Matrix Computations (Johns Hopkins University Press, Baltimore, 1996).
Y. Saad, SIAM J. Numer. Anal. 29 (1992) 209.
R.B. Sidje, ACM TOMS 24 (1998) 130.
D.E. Knuth, The Art of Computer Programming, Vol. 2, Seminumerical Algorithms, 2nd ed. (Addison-Wesley, Reading, MA, 1981).
W. Kahan, CACM 8 (1965) 40.
W. Kahan, in: Proceedings of the IFIP Congress, Information Processing 71, Ljubljana, August 1971, eds. C.V. Freiman, J.E. Griffith and J.L. Rosenfeld (North-Holland, Amsterdam, 1972) pp. 1214–1239.
E. Schäfer, J. Math. Biol. 2 (1975) 41.
E. Hairer and G. Wanner, Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2nd ed. (Springer-Verlag, New York, 1996).
M. Tomita, K. Hashimoto, K. Takahashi, T.S. Shimizu, Y. Matsuzaki, F. Miyoshi, K. Saito, S. Tanida, K. Yugi, J.C. Venter and C.A. Hutchison, III, Bioinformatics 15 (1999) 72.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Aluffi-Pentini, F., De Fonzo, V. & Parisi, V. A Novel Algorithm for the Numerical Integration of Systems of Ordinary Differential Equations Arising in Chemical Problems. Journal of Mathematical Chemistry 33, 1–15 (2003). https://doi.org/10.1023/A:1023295712640
Issue Date:
DOI: https://doi.org/10.1023/A:1023295712640