A mixed discontinuous/continuous finite element pair for shallow-water ocean modelling

Abstract

We introduce a mixed discontinuous/continuous finite element pair for ocean modelling, with continuous quadratic pressure/layer depth and discontinuous velocity. We investigate the finite element pair applied to the linear shallow-water equations on an f-plane. The element pair has the property that all geostrophically balanced states which strongly satisfy the boundary conditions have discrete divergence equal to exactly zero and hence are exactly steady states of the discretised equations. This means that the finite element pair has excellent geostrophic balance properties. We illustrate these properties using numerical tests and provide convergence calculations which show that the discretisation has quadratic errors, indicating that the element pair is stable.

Cited by in corpus (12)

• Thetis coastal ocean model: discontinuous Galerkin discretization for the three-dimensional hydrostatic equations
• Mixed finite elements for numerical weather prediction
• A framework for automated PDE-constrained optimisation
• Numerical wave propagation for the triangular $P{1}_{DG}$-$P2$ finite element pair
• Solving the Poisson equation on small aspect ratio domains using unstructured meshes
• An implicit wetting and drying approach for non-hydrostatic baroclinic flows in high aspect ratio domains
• Well balanced finite volume schemes for shallow water equations on manifolds
• Flexible, Scalable Mesh and Data Management using PETSc DMPlex
• Reconstructing wave profiles from inundation data
• Geostrophic balance preserving interpolation in mesh adaptive shallow-water ocean modelling
• The representation of boundary currents in a finite element shallow water model
• A family of mixed finite element pairs with optimal geostrophic balance