Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers

A. Sandu, J.G. Verwer, J.G. Blom, E.J. Spee, G.R. Carmichael, F.A. Potra
1997 Atmospheric Environment  
In the numerical simulation of atmospheric transport-chemistry processes, a major task is the integration of the stiff systems of ordinary differential equations describing the chemical transformations. It is therefore of interest to systematically search for stiff solvers which can be identified as close to optimal for atmospheric applications. In this paper we continue our investigation from Sandu et al. (1996, CWI Report NM-R9603 and Report in Comput. Math., No. 85) and compare eight solvers
more » ... on a set of seven box-models used in present day models. The focus is on Rosenbrock solvers. These tum out to be very well suited for our application when they are provided with highly efficient sparse matrix techniques to economize on the linear algebra. Two of the Rosenbrock solvers tested are from the literature, viz. RODAS and ROs4, and two are new and specially developed for air quality applications, viz. RODAS3 and ROS3. (f) 1997 Elsevier Science Ltd. Key word index: Atmospheric chemistry, air pollution modeling, numerical stiff ODEs, Rosenbrock methods, sparsity. VODE (Brown et al., 1989) and LSODES (Hindmarsh, 1983), the Runge-Kutta solver SDIRK4 (Hairer and Wanner, 1991) and the Rosenbrock solver RODAS (Hairer and Wanner, 1991) .
doi:10.1016/s1352-2310(97)83212-8 fatcat:gov5e3wj2bfldcf5qth4az3sv4