Positive and conservative schemes for mass action kinetics

Abstract Some implicit schemes for the discretization of mass action kinetics are presented and discussed. Although all the schemes are classical, the solution procedure takes great advantage of the peculiar structure of the ODE. All the schemes maintain linear first integrals such as atomic mass conservation. The main features of the first-order scheme is positivity preservation for arbitrarily large time steps, while for the higher order schemes, there is numerical evidence that this property is maintained with fairly large time steps. For the first-order scheme, the existence of a nonnegative solution of the algebraic system arising in the advancing step is proved for arbitrarily large time steps, while uniqueness is proved for limited time steps. Finally, an efficient procedure is presented for solving the nonlinear systems involved in the time step. These systems are solved without the use of Jacobian matrices or their approximations, but by the repeated inversion of M -matrices, a procedure which is both easy and fast, considerably simplifying and accelerating the computer implementation of the schemes.

[1]  A. Bellen,et al.  Contractivity of waveform relaxation Runge-Kutta iterations and related limit methods for dissipative systems in the maximum norm , 1994 .

[2]  Philip W. Sharp,et al.  Two-stage and Three-stage Diagonally Implicit Runge-Kutta Nyström Methods of Orders Three and Four , 1990 .

[3]  G. J. Cooper,et al.  Some schemes for the implementation of implicit Runge-Kutta methods , 1993 .

[4]  R. Alexander,et al.  Stability of Runge-Kutta methods for stiff ordinary differential equations , 1994 .

[5]  L. Collatz.  The numerical treatment of differential equations , 1961 .

[6]  E. Zeidler.  Nonlinear functional analysis and its applications , 1988 .

[7]  Kevin Burrage,et al.  On smoothing and order reduction effects for implicit Runge-Kutta formulae , 1993 .

[8]  M. N. Spijker,et al.  On the algebraic equations in implicit Runge-Kutta methods , 1987 .

[9]  R. Vermiglio.  Multistep high-order interpolants of Runge-Kutta methods , 1993 .

[10]  Kevin Burrage,et al.  Efficiently Implementable Algebraically Stable Runge–Kutta Methods , 1982 .

[11]  R. Scherer,et al.  Complete Algebraic Characterization of A-Stable Runge–Kutta Methods , 1994 .

[12]  K. Burrage.  The dichotomy of stiffness: pragmatism versus theory , 1989, Conference on Numerical Ordinary Differential Equations.

[13]  Alston S. Householder,et al.  The Theory of Matrices in Numerical Analysis , 1964 .

[14]  Winfried Auzinger,et al.  Modern convergence theory for stiff initial-value problems , 1993 .

[15]  Péter Érdi,et al.  Mathematical models of chemical reactions , 1989 .

[16]  Y. Chang.  Solving STIFF systems by Taylor series , 1989, Conference on Numerical Ordinary Differential Equations.

[17]  Peter Albrecht,et al.  A new theoretical approach to Runge-Kutta methods , 1987 .

[18]  R. Alexander.  Diagonally implicit runge-kutta methods for stiff odes , 1977 .