Polynomial robust stability analysis for $H(\text{div})$-conforming finite elements for the Stokes equations

In this work we consider a discontinuous Galerkin method for the discretization of the Stokes problem. We use $H(\text{div})$-conforming finite elements as they provide major benefits such as exact mass conservation and pressure-independent error estimates. The main aspect of this work lies in the analysis of high order approximations. We show that the considered method is uniformly stable with respect to the polynomial order $k$ and provides optimal error estimates $\mathrm{\parallel }\mathit{u}-{\mathit{u}}_h{\mathrm{\parallel }}_{1_h}+\mathrm{\parallel }{\mathrm{\Pi }}^{Q_h}p-p_h\mathrm{\parallel }\le c{\left(h\mathrm{/}k\right)}^s\mathrm{\parallel }\mathit{u}{\mathrm{\parallel }}_{s+1}$. To derive those estimates, we prove a $k$-robust LBB condition. This proof is based on a polynomial $H^2$-stable extension operator. This extension operator itself is of interest for the numerical analysis of $C^0$-continuous discontinuous Galerkin methods for $4^{th}$ order problems.