Hybrid Discontinuous Galerkin methods with relaxed H(div)-conformity for incompressible flows. Part II

The present work is the second part of a pair of papers, considering Hybrid Discontinuous Galerkin methods with relaxed H(div)-conformity. The first part mainly dealt with presenting a robust analysis with respect to the mesh size h and the introduction of a reconstruction operator to restore divergence-conformity and pressure robustness (pressure independent velocity error estimates) using a modified force discretization. The aim of this part is the presentation of a high order polynomial robust analysis for the relaxed H(div)-conforming Hybrid Discontinuous Galerkin discretization of the two dimensional Stokes problem. It is based on the recently proven polynomial robust LBB-condition for BDM elements, Lederer and Schöberl (IMA J. Numer. Anal. (2017)) and is derived by a direct approach instead of using a best approximation Céa like result. We further treat the impact of the reconstruction operator on the hp analysis and present a numerical investigation considering polynomial robustness. We conclude the paper presenting an efficient operator splitting time integration scheme for the Navier–Stokes equations which is based on the methods recently presented in Lehrenfeld and Schöberl (Comp. Methods Appl. Mech. Eng. 307 (2016) 339–361) and includes the ideas of the reconstruction operator.