Christian Merdon

Philip L. Lederer, Christoph Lehrenfeld, Christian Merdon et al.

This paper studies pressure-robustness for the axisymmetric Stokes problem. The transformation to cylindrical coordinates requires that the radially weighted velocity is divergence-free in the classical sense. Consequently, traditional divergence-free finite element methods from the Cartesian setting -- even if inf-sup stable -- are in general not divergence-free in the axisymmetric formulation. We therefore explore the approach that restores pressure-robustness via reconstruction operators for a low-order Bernardi--Raugel discretization. We show that an application of standard interpolation operators from the Cartesian setting to radially weighted test functions works in principle, but it lacks properties needed to derive optimal consistency error estimates. To address this, we introduce a reconstruction operator into a finite element space spanned by Raviart--Thomas functions that are modified such that they vanish on the rotation axis. This vanishing-on-axis property is the key to obtain optimal consistency error estimates. Numerical examples demonstrate the overall feasibility of the approach and include cases where the vanishing-on-axis property yields significantly better results.

Philip L. Lederer, Alexander Linke, Christian Merdon et al.

Classical inf-sup stable mixed finite elements for the incompressible (Navier-)Stokes equations are not pressure-robust, i.e., their velocity errors depend on the continuous pressure. However, a modification only in the right hand side of a Stokes discretization is able to reestablish pressure-robustness, as shown recently for several inf-sup stable Stokes elements with discontinuous discrete pressures. In this contribution, this idea is extended to low and high order Taylor-Hood and mini elements, which have continuous discrete pressures. For the modification of the right hand side a velocity reconstruction operator is constructed that maps discretely divergence-free test functions to exactly divergence-free ones. The reconstruction is based on local $H(\mathrm{div})$-conforming flux equilibration on vertex patches, and fulfills certain orthogonality properties to provide consistency and optimal a-priori error estimates. Numerical examples for the incompressible Stokes and Navier-Stokes equations confirm that the new pressure-robust Taylor-Hood and mini elements converge with optimal order and outperform significantly the classical versions of those elements when the continuous pressure is comparably large.