Abstract:We present an energy- and potential enstrophy-conserving scheme for the non-traditional shallow water equations that include the complete Coriolis force and topography. These integral conservation properties follow from material conservation of potential vorticity in the continuous shallow water equations. The latter property cannot be preserved by a discretisation on a fixed Eulerian grid, but exact conservation of a discrete energy and a discrete potential enstrophy seems to be an effective substitute that prevents any distortion of the forward and inverse cascades in quasi-two dimensional turbulence through spurious sources and sinks of energy and potential enstrophy, and also increases the robustness of the scheme against nonlinear instabilities. We exploit the existing Arakawa–Lamb scheme for the traditional shallow water equations, reformulated by Salmon as a discretisation of the Hamiltonian and Poisson bracket for this system. The non-rotating, traditional, and our non-traditional shallow water equations all share the same continuous Hamiltonian structure and Poisson bracket, provided one distinguishes between the particle velocity and the canonical momentum per unit mass. We have determined a suitable discretisation of the non-traditional canonical momentum, which includes additional coupling between the layer thickness and velocity fields, and modified the discrete kinetic energy to suppress an internal symmetric computational instability that otherwise arises for multiple layers. The resulting scheme exhibits the expected second-order convergence under spatial grid refinement. We also show that the drifts in the discrete total energy and potential enstrophy due to temporal truncation error may be reduced to machine precision under suitable refinement of the timestep using the third-order Adams–Bashforth or fourth-order Runge–Kutta integration schemes.

A Numerical Scheme for Wave Turbulence: 3-Wave Kinetic Equations

Numerical Schemes for 3-Wave Kinetic Equations: A Complete Treatment of the Collision Operator

Numerical schemes for a fully nonlinear coagulation-fragmentation model coming from wave kinetic theory

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

NUMERICAL ANALYSIS OF AN ENERGY-CONSERVATION SCHEME FOR TWO-DIMENSIONAL HAMILTONIAN WAVE EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

A numerical method with properties of consistency in the energy domain for a class of dissipative nonlinear wave equations with applications to a Dirichlet boundary-value problem

An energy and potential enstrophy conserving numerical scheme for the multi-layer shallow water equations with complete Coriolis force

Strongly stable dual-pairing summation by parts finite difference schemes for the vector invariant nonlinear shallow water equations -- I: Numerical scheme and validation on the plane

A numerical method for wave-structure interactions in the Boussinesq regime

Efficient numerical approximations for a non-conservative Nonlinear Schrodinger equation appearing in wind-forced ocean waves

Two-dimensional Isotropic Inertia–gravity Wave Turbulence

An Energy-Preserving Algorithm for Nonlinear Hamiltonian Wave Equations with Neumann Boundary Conditions

Numerical solutions of the nonlinear wave equations with energy-preserving sixth-order finite difference schemes

An entropy stable discontinuous Galerkin method for the spherical thermal shallow water equations

Efficient numerical approximations for a nonconservative nonlinear Schrödinger equation appearing in wind‐forced ocean waves

A theory of turbulence. Part 1:Towards solutions of the Navier-stokes equations

Staggered finite volume scheme for solving cascade flow with a k-omega turbulence model

A finite difference scheme for smooth solutions of the general mKdV equation

A Bound Preserving Energy Stable Scheme for a Nonlocal Cahn-Hilliard Equation

Simulation of the Wave Turbulence of a Liquid Surface Using the Dynamic Conformal Transformation Method

A numerical method for two-phase flow with its application to the Kelvin Helmholtz instability problem