A matrix series method for the integration of rate equations in a reaction network. An alternative to Runge-Kutta methods

Thomas M. Zamis, Lawrence J. Parkhurst, Gordon A. Gallup

Research output: Contribution to journalArticle

13 Scopus citations

Abstract

A matrix method is presented for obtaining numerical solutions to the differential equations pertaining to reaction networks. Neither matrix inversions nor diagonalizations are required. The kinetic equations are east in matrix and vector form and solved as an exponential matrix, which is approximated by a truncated series. The series matrix. S(t) is in terms of the rate constant matrix K that characterizes the reaction network. Accuracy for a given step size is determined by the order of the approximating series. Optimum coefficients for the series are given through the fifth order. Concentrations of all species and sensitivity coefficients are calculated for each time interval. The procedure for first order networks is particularly simple when all times are equally spaced, since S remains fixed throughout the problem and only the concentration vector changes. A second order process is approximated over a small time interval by a first order expression, but both S and the concentration vector change at each step. The procedure is approx. 2-3 times more rapid than conventional fourth order Runge-Kutta integration to the same level of accuracy.

Original languageEnglish (US)
Pages (from-to)165-171
Number of pages7
JournalComputers and Chemistry
Volume13
Issue number3
DOIs
StatePublished - 1989

ASJC Scopus subject areas

  • Biotechnology
  • Applied Microbiology and Biotechnology
  • Chemical Engineering(all)

Fingerprint Dive into the research topics of 'A matrix series method for the integration of rate equations in a reaction network. An alternative to Runge-Kutta methods'. Together they form a unique fingerprint.

  • Cite this