High-Order Multirate Explicit Time-Stepping Schemes for the Baroclinic-Barotropic Split Dynamics in Primitive Equations

05/27/2021
by   Rihui Lan, et al.
0

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. In this paper, we propose second and third-order multirate explicit time-stepping schemes for such split systems based on the strong stability-preserving Runge-Kutta (SSPRK) framework. Our method allows for a large time step to be used for advancing the three-dimensional (slow) baroclinic mode and a small time step for the two-dimensional (fast) barotropic mode, so that each of the two mode solves only need satisfy their respective CFL condition to maintain numerical stability. It is well known that the SSPRK method achieves high-order temporal accuracy by utilizing a convex combination of forward-Euler steps. At each time step of our method, the baroclinic velocity is first computed by using the SSPRK scheme to advance the baroclinic-barotropic system with the large time step, then the barotropic velocity is specially corrected by using the same SSPRK scheme with the small time step to advance the barotropic subsystem with a barotropic forcing interpolated based on values from the preceding baroclinic solves. Finally, the fluid thickness and the sea surface height perturbation is updated by coupling the predicted baroclinic and barotropic velocities. Temporal truncation error analyses are also carried out for the proposed schemes. Two benchmark tests drawn from the -MPAS-Ocean" platform are used to numerically demonstrate the accuracy and parallel performance of the proposed schemes.

READ FULL TEXT

page 18

page 19

page 21

page 22

research
04/01/2021

Stability and error analysis of IMEX SAV schemes for the magneto-hydrodynamic equations

We construct and analyze first- and second-order implicit-explicit (IMEX...
research
09/07/2017

Adaptive PCA for Time-Varying Data

In this paper, we present an online adaptive PCA algorithm that is able ...
research
02/20/2022

Adomian Decomposition Based Numerical Scheme for Flow Simulations

Time efficiency is one of the more critical concerns in computational fl...
research
08/21/2020

Exact targeting of Gibbs distributions using velocity-jump processes

This work introduces and studies a new family of velocity jump Markov pr...
research
10/10/2019

Positivity-preserving finite difference WENO scheme for Ten-Moment equations with source term

We develop a positivity-preserving finite difference WENO scheme for the...
research
10/05/2019

Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models

The tracer equations are part of the primitive equations used in ocean m...

Please sign up or login with your details

Forgot password? Click here to reset