Numerical Analysis of a Stable Discontinuous Galerkin Scheme for the Hydrostatic Stokes Problem

We propose a Discontinuous Galerkin (DG) scheme for the Hydrostatic Stokes equations. These equations, related to large-scale PDE models in Oceanography, are characterized by the loss of ellipticity of the vertical momentum equation. This fact provides some interesting challenges, such as the design of stable numerical schemes. The new scheme proposed here is based on the SIP penalty technique, with a particular treatment of the vertical velocity. It is well-known that stability of the mixed formulation of Primitive Equations requires, besides the LBB inf-sup condition, an additional hydrostatic inf-sup restriction relating pressure and vertical velocity. This hydrostatic inf-sup condition invalidates stability of usual Stokes stable continuous finite elements like Taylor-Hood 𝒫2/𝒫1 or bubble 𝒫1b/𝒫1. Here we prove stability for our 𝒫k/𝒫k DG scheme. Some novel numerical tests are provided which are in agreement with the previous analysis.