The Numerical Solution of Some Kinetics Models with VODE and CHEMKIN II

Abstract We describe some numerical experiments we have performed in solving chemical kinetics models with CHEMKIN II and two high quality ordinary differential equations (ODE) packages—VODE and LSODE. For this mode of kinetics modeling, we give performance figures for these codes and for three chemical systems of varying size and behavior. In all cases, the Jacobian matrix is internally computed by the ODE solver, owing to the complexity of the models. Although these results are somewhat more conservative than others we have obtained, they do show that significant performance improvement can be obtained by using VODE instead of LSODE for these problems. We believe that the problems are typical of those that one of us (AMD) has been working with for several years. The main objective of this report is to make kineticists and other practitioners aware of VODE and its potential advantages in such applications. Performance results are given. For example, a small mechanism (17 species and 31 reactions) required 0.48 s CPU time with VODE vs 0.83 s with LSODE. A larger mechanism (135 species and 497 reactions) showed a more dramatic difference: 60 s CPU time for VODE vs 213 s with LSODE.