A DG Method for the Stokes Equations on Tensor Product Meshes with $[P_k]^d-P_{k-1}$ Element

Abstract

We consider the mixed discontinuous Galerkin (DG) finite element approximation of the Stokes equation and provide the analysis for the $[P_k]^d-P_{k-1}$ element on the tensor product meshes. Comparing to the previous stability proof with $[Q_k]^d-Q_{k-1}$ discontinuous finite elements in the existing references, our first contribution is to extend the formal proof to the $[P_k]^d-P_{k-1}$ discontinuous elements on the tensor product meshes. Numerical inf-sup tests have been performed to compare $Q_k$ and $P_k$ types of elements and validate the well-posedness in both settings. Secondly, our contribution is to design the enhanced pressure-robust discretization by only modifying the body source assembling on $[P_k]^d-P_{k-1}$ schemes to improve the numerical simulation further. The produced numerical velocity solution via our enhancement shows viscosity and pressure independence and thus outperforms the solution produced by standard discontinuous Galerkin schemes. Robustness analysis and numerical tests have been provided to validate the scheme’s robustness.

Funding

Simons Foundation(964995)