An Implementation-Friendly SDG Scheme based on Cartesian Grids for Stokes Equations with Pressure Robustness and Superconvergence
Abstract
This paper develops a staggered discontinuous Galerkin (SDG) scheme based on Cartesian grids for Stokes equations that is simple to implement, intrinsically pressure-robust, and superconvergent for all variables. Instead of the composite meshes used in standard SDG, we construct staggered quadrilateral meshes from Cartesian grids. The scheme takes velocity, pressure, and velocity gradient as unknowns, and employs piecewise-constant spaces with carefully designed staggered continuity. The additional gradient unknowns are locally eliminated via static condensation and can be further removed by mass lumping without loss of accuracy. An explicit pointwise formulation of the scheme is derived, which facilitates implementation and enables a detailed pointwise analysis. We rigorously prove pressure robustness and second-order superconvergence, which holds on general non-uniform Cartesian grids. The scheme is further extended to Navier-Stokes equations by introducing a novel discrete convection term with second-order consistency. Combined with the scalar auxiliary variable (SAV) approach and the Crank-Nicolson (CN) scheme, this yields an unconditionally energy-stable and second-order accurate scheme. Numerical experiments validate the theory and demonstrate accuracy and robustness.
Related
- Wavenumber-explicit stability of three-dimensional Cartesian UPMLs with piecewise-constant profiles for Helmholtz and Maxwell equations
- A ghost-cell immersed boundary method for the shallow water equations on Cartesian grids: mass conservation and wetting and drying over variable topography
- A compact high-order and positivity-preserving active flux method for compressible Navier--Stokes equations on Cartesian meshes
- Polynomial preserving recoveries of edge element method on Cartesian grids for the time-harmonic Maxwell equations with large wave number
- Admissible characterization in delay equations and robustness of exponential dichotomies against small-delay perturbations
BibTeX
@misc{yang2026implementationfriendly,
title = {{An Implementation-Friendly SDG Scheme based on Cartesian Grids for Stokes Equations with Pressure Robustness and Superconvergence}},
author = {Bohan Yang and Eric Chung},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.24232},
eprint = {2609.24232},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.24232},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - An Implementation-Friendly SDG Scheme based on Cartesian Grids for Stokes Equations with Pressure Robustness and Superconvergence
AU - Bohan Yang
AU - Eric Chung
PY - 2026
DA - 2026/09/21
PB - arXiv
DO - 10.48550/arxiv.2609.24232
UR - https://arxiv.org/abs/2609.24232
AB - This paper develops a staggered discontinuous Galerkin (SDG) scheme based on Cartesian grids for Stokes equations that is simple to implement, intrinsically pressure-robust, and superconvergent for all variables. Instead of the composite meshes used in standard SDG, we construct staggered quadrilateral meshes from Cartesian grids. The scheme takes velocity, pressure, and velocity gradient as unknowns, and employs piecewise-constant spaces with carefully designed staggered continuity. The additional gradient unknowns are locally eliminated via static condensation and can be further removed by mass lumping without loss of accuracy. An explicit pointwise formulation of the scheme is derived, which facilitates implementation and enables a detailed pointwise analysis. We rigorously prove pressure robustness and second-order superconvergence, which holds on general non-uniform Cartesian grids. The scheme is further extended to Navier-Stokes equations by introducing a novel discrete convection term with second-order consistency. Combined with the scalar auxiliary variable (SAV) approach and the Crank-Nicolson (CN) scheme, this yields an unconditionally energy-stable and second-order accurate scheme. Numerical experiments validate the theory and demonstrate accuracy and robustness.
ER -