Continuous interior penalty stabilization for divergence-free finite element methods

In this paper, we propose, analyze and test numerically a pressure-robust stabilized finite element for a linearized problem in incompressible fluid mechanics, namely, the steady Oseen equation with low viscosity. Stabilization terms are defined by jumps of different combinations of derivatives for the convective term over the element faces of the triangulation of the domain. With the help of these stabilizing terms, and the fact the finite element space is assumed to provide a point-wise divergence-free velocity, an $\mathcal O\big(h^{k+\frac 12}\big)$ error estimate in the $L^2$-norm is proved for the method (in the convection-dominated regime), and optimal order estimates in the remaining norms of the error. Numerical results supporting the theoretical findings are provided.