Further development of efficient and accurate time integration schemes for meteorological models
arXiv:1805.02144 · doi:10.1016/j.jcp.2018.10.018
Abstract
In this paper, we investigate the use of higher-order exponential Rosenbrock time integration methods on the shallow water equations on the sphere. This stiff, nonlinear model provides a testing ground for accurate and stable time integration methods in weather modeling, serving as the focus for exploration of novel methods for many years. We therefore identify a candidate set of three recent exponential Rosenbrock methods of orders four and five (exprb42, pexprb43 and exprb53) for use on this model. Based on their multi-stage structure, we propose a set of modifications to the phipm_IOM2 algorithm for efficiently calculating the matrix phi-functions. We then investigate the performance of these methods on a suite of four challenging test problems, comparing them against the epi3 method investigated previously in [1, 2] on these problems. In all cases, the proposed methods enable accurate solutions at much longer time-steps than epi3, proving considerably more efficient as either the desired solution error decreases, or as the test problem nonlinearity increases.
44 pages, 11 figures, 7 tables (accepted version), Journal of Computational Physics
References in corpus (1)
Cited by in corpus (5)
- High-order numerical solutions to the shallow-water equations on the rotated cubed-sphere grid
- Exponential time differencing for mimetic multilayer ocean models
- LeXInt: Package for Exponential Integrators employing Leja interpolation
- Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models
- A scalable exponential-DG approach for nonlinear conservation laws: with application to Burger and Euler equations