Exponential versus IMEX high-order time integrators for thermal convection in rotating spherical shells

Keywords

SC4I/Digitization, Innovation, and Competitiveness of the Production System
Code:
36/2013
Title:
Exponential versus IMEX high-order time integrators for thermal convection in rotating spherical shells
Date:
Sunday 25th August 2013
Author(s):
Ferran Garcia, Luca Bonaventura, Marta Net, Juan Sanchez
Download link:
Abstract:
We assess the accuracy and efficiency of several exponential time integration methods coupled to a spectral discretization of the three-dimensional Boussinesq thermal convection equations in rotating spherical shells. Exponential methods are compared to implicit-explicit (IMEX) multi-step methods. The results of a wide range of numerical simulations highlight the superior accuracy of exponential methods for a given time step, especially when employed with large time steps and at low Ekman number. However, presently available implementations of exponential methods appear to be in general computationally more expensive than those of IMEX methods and further research is needed to reduce their computational cost per time step. A physically justified extrapolation argument suggests that some exponential methods could be the most efficient option for integrating flows near Earth’s outer core conditions.
This report, or a modified version of it, has been also submitted to, or published on
Journal of Computational Physics, Vol. 264, pp. 41-54, 2014