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 Citations (Scopus)

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

Fingerprint

Reaction Network
Runge Kutta methods
Rate Equations
Runge-Kutta Methods
Series
Alternatives
First-order
Sensitivity Coefficient
Matrix Exponential
Interval
Matrix Inversion
Diagonalization
Matrix Method
Runge-Kutta
Kinetic Equation
Rate Constant
Fourth Order
Numerical Solution
Differential equation
Rate constants

ASJC Scopus subject areas

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

Cite this

A matrix series method for the integration of rate equations in a reaction network. An alternative to Runge-Kutta methods. / Zamis, Thomas M.; Parkhurst, Lawrence J; Gallup, Gordon A.

In: Computers and Chemistry, Vol. 13, No. 3, 1989, p. 165-171.

Research output: Contribution to journalArticle

@article{5d535b2da0614478a485832134e52ff6,
title = "A matrix series method for the integration of rate equations in a reaction network. An alternative to Runge-Kutta methods",
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.",
author = "Zamis, {Thomas M.} and Parkhurst, {Lawrence J} and Gallup, {Gordon A.}",
year = "1989",
doi = "10.1016/0097-8485(89)85002-8",
language = "English (US)",
volume = "13",
pages = "165--171",
journal = "Computational Biology and Chemistry",
issn = "1476-9271",
publisher = "Elsevier Limited",
number = "3",

}

TY - JOUR

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

AU - Zamis, Thomas M.

AU - Parkhurst, Lawrence J

AU - Gallup, Gordon A.

PY - 1989

Y1 - 1989

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=0024915465&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0024915465&partnerID=8YFLogxK

U2 - 10.1016/0097-8485(89)85002-8

DO - 10.1016/0097-8485(89)85002-8

M3 - Article

AN - SCOPUS:0024915465

VL - 13

SP - 165

EP - 171

JO - Computational Biology and Chemistry

JF - Computational Biology and Chemistry

SN - 1476-9271

IS - 3

ER -