skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Comparison of exponential integrators and traditional time integration schemes for the shallow water equations

Journal Article · · Applied Numerical Mathematics

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

References (31)

KIOPS: A fast adaptive Krylov subspace solver for exponential integrators journal November 2018
Approximation of the matrix exponential for matrices with a skinny field of values journal May 2020
GPU Accelerated Algorithms for Computing Matrix Function Vector Products with Applications to Exponential Integrators and Fractional Diffusion journal January 2016
On the time-splitting scheme used in the Princeton Ocean Model journal May 2009
Nineteen Dubious Ways to Compute the Exponential of a Matrix journal October 1978
Stability constraints for oceanic numerical models: implications for the formulation of time and space discretizations journal August 2015
Exponential time differencing for the tracer equations appearing in primitive equation ocean models journal June 2020
On Krylov Subspace Approximations to the Matrix Exponential Operator journal October 1997
An efficient exponential time integration method for the numerical solution of the shallow water equations on the sphere journal October 2016
On the use of exponential time integration methods in atmospheric models journal May 2013
Runge-Kutta methods with minimum error bounds journal January 1962
A Lawson-type exponential integrator for the Korteweg–de Vries equation journal October 2019
Useful time-stepping methods for the Coriolis term in a shallow water model journal January 2009
Exponential time differencing for mimetic multilayer ocean models journal December 2019
Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow. Part I journal August 1966
Bedload transport in shallow water models: Why splitting (may) fail, how hyperbolicity (can) help journal August 2011
Exponential integrators with parallel-in-time rational approximations for the shallow-water equations on the rotating sphere journal July 2019
Exponential versus IMEX high-order time integrators for thermal convection in rotating spherical shells journal May 2014
Exponential integrators journal May 2010
Exponential Rosenbrock-Type Methods journal January 2009
Domain decomposition based exponential time differencing method for fluid dynamics problems with smooth solutions journal November 2019
Preconditioned implicit-exponential integrators (IMEXP) for stiff PDEs journal April 2017
Compact Finite Difference Schemes for Ocean Models journal November 2000
Algorithm 919: A Krylov Subspace Algorithm for Evaluating the ϕ-Functions Appearing in Exponential Integrators journal April 2012
Analysis of Some Krylov Subspace Approximations to the Matrix Exponential Operator journal February 1992
Exponential integrators preserving local conservation laws of PDEs with time-dependent damping/driving forces journal May 2019
Generalized Runge-Kutta Processes for Stable Systems with Large Lipschitz Constants journal September 1967
The scaling and modified squaring method for matrix functions related to the exponential journal March 2009
The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model journal January 2005
Stability analysis of split-explicit free surface ocean models: Implication of the depth-independent barotropic mode approximation journal December 2019
A standard test set for numerical approximations to the shallow water equations in spherical geometry journal September 1992