Comparison of exponential integrators and traditional time integration schemes for the shallow water equations
- University of Grenoble Alpes, Grenoble (France)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
We report the time integration scheme is probably one of the most fundamental choices in the development of an ocean model. In this paper, we investigate several time integration schemes when applied to the shallow water equations. This set of equations is accurate enough for the modeling of a shallow ocean and is also relevant to study as it is the one solved for the barotropic (i.e. vertically averaged) component of a three dimensional ocean model. We analyze different time stepping algorithms for the linearized shallow water equations. High order explicit schemes are accurate but the time step is constrained by the Courant-Friedrichs-Lewy stability condition. Implicit schemes can be unconditionally stable but, in practice lack accuracy when used with large time steps. In this paper we propose a detailed comparison of such classical schemes with exponential integrators. The accuracy and the computational costs are analyzed in different configurations.
- Research Organization:
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); French National Research Agency; European Union’s Horizon 2020
- Grant/Contract Number:
- NA0003525; ANR-18-CE46-0008; 821926
- OSTI ID:
- 1872031
- Report Number(s):
- SAND2022-6879J; 706614
- Journal Information:
- Applied Numerical Mathematics, Vol. 180; ISSN 0168-9274
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Development of a diffusive wave shallow water model with a novel stability condition and other new features
A high order characteristic discontinuous Galerkin scheme for advection on unstructured meshes