Premium
Fast and reliable numerical methods to simulate complex chemical kinetic mechanisms
Author(s) -
Olcese Luis E.,
Toselli Beatriz M.
Publication year - 1998
Publication title -
international journal of chemical kinetics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.341
H-Index - 68
eISSN - 1097-4601
pISSN - 0538-8066
DOI - 10.1002/(sici)1097-4601(1998)30:5<349::aid-kin5>3.0.co;2-s
Subject(s) - ode , solver , ordinary differential equation , computer science , problem solver , residual , mathematics , construct (python library) , stiff equation , numerical integration , explicit and implicit methods , l stability , algorithm , differential equation , computational science , differential algebraic equation , mathematical analysis , programming language
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