Continuity properties of the inf-sup constant for the divergence

The inf-sup constant for the divergence, or LBB constant, is explicitly known for only few domains. For other domains, upper and lower estimates are known. If more precise values are required, one can try to compute a numerical approximation. This involves, in general, approximation of the domain and then the computation of a discrete LBB constant that can be obtained from the numerical solution of an eigenvalue problem for the Stokes system. This eigenvalue problem does not fall into a class for which standard results about numerical approximations can be applied. Indeed, many reasonable finite element methods do not yield a convergent approximation. In this article, we show that under fairly weak conditions on the approximation of the domain, the LBB constant is an upper semi-continuous shape functional, and we give more restrictive sufficient conditions for its continuity with respect to the domain. For numerical approximations based on variational formulations of the Stokes eigenvalue problem, we also show upper semi-continuity under weak approximation properties, and we give stronger conditions that are sufficient for convergence of the discrete LBB constant towards the continuous LBB constant. Numerical examples show that our conditions are, while not quite optimal, not very far from necessary

A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow

We present an efficient discontinuous Galerkin scheme for simulation of the incompressible Navier-Stokes equations including laminar and turbulent flow. We consider a semi-explicit high-order velocity-correction method for time integration as well as nodal equal-order discretizations for velocity and pressure. The non-linear convective term is treated explicitly while a linear system is solved for the pressure Poisson equation and the viscous term. The key feature of our solver is a consistent penalty term reducing the local divergence error in order to overcome recently reported instabilities in spatially under-resolved high-Reynolds-number flows as well as small time steps. This penalty method is similar to the grad-div stabilization widely used in continuous finite elements. We further review and compare our method to several other techniques recently proposed in literature to stabilize the method for such flow configurations. The solver is specifically designed for large-scale computations through matrix-free linear solvers including efficient preconditioning strategies and tensor-product elements, which have allowed us to scale this code up to 34.4 billion degrees of freedom and 147,456 CPU cores. We validate our code and demonstrate optimal convergence rates with laminar flows present in a vortex problem and flow past a cylinder and show applicability of our solver to direct numerical simulation as well as implicit large-eddy simulation of turbulent channel flow at $Re_{\tau}=180$ as well as $590$.Comment: 28 pages, in preparation for submission to Journal of Computational Physic