| Abstract: | Models that simulate atmospheric photochemistry require the use of a stiff ordinary differential equations (ODEs) solver. Since the simulation of the chemical transformations taking place in the system takes up to 80 percent of the CPU time, the numerical solver must be computationally fast. Also, the residual error from the solver must be small. Because most accurate solvers are relatively slow, modelers continue to search for timely, yet accurate integration methods. Over the past years an extensive number of articles have been dedicated to this subject. One of the highly debated questions is whether one should construct specialized algorithms or instead use general methods for stiff ODEs. In the present article we use the second alternative. We apply three linearly (semi-)implicit methods from the classical stiff ODE literature which we modified to implement the sparse routines to solve the system of equations describing a complex kinetic mechanism. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet 30: 349–358, 1998 |
