Abstract

Implicit integrators are very useful in efficiently solving stiff systems of ODEs arising from atmospheric chemistry kinetics, provided that they are modified to take full advantage of the structure of the problem. A systematic way of treating sparsity for reducing the linear algebra cost is presented. 3 Efficient implementation of fully implicit methods 4 1 Introduction It is well known that the equations arising from chemical kinetics comprise a system of stiff ordinary differential equations (ODE). For solving these equations numerically, implicit integrators with infinite stability regions are likely to work with relatively large step-sizes when the accuracy requirements are not too stringent. However, at each integration step, a nonlinear system of equations has to be solved. This involves the repeated evaluation of Jacobians and the solution of linear algebraic systems of dimension n, the number of species considered in the model. General stiff ODE solvers do not take advanta...

Citations