Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations

We present optimal preconditioners for a recently introduced hybridized discontinuous Galerkin finite element discretization of the Stokes equations. Typical of hybridized discontinuous Galerkin methods, the method has degrees-of-freedom that can be eliminated locally (cell-wise), thereby significantly reducing the size of the global problem. Although the linear system becomes more complex to analyze after static condensation of these element degrees-of-freedom, the pressure Schur complement of the original and reduced problem are the same. Using this fact, we prove spectral equivalence of this Schur complement to two simple matrices, which is then used to formulate optimal preconditioners for the statically condensed problem. Numerical simulations in two and three spatial dimensions demonstrate the good performance of the proposed preconditioners

Rellich-type Discrete Compactness for Some Discontinuous Galerkin FEM

We deduce discrete compactness of Rellich type for some discontinuous Galerkin finite element methods (DGFEM) including hybrid ones, under fairly general settings on the triangulations and the finite element spaces. We make use of regularity of the solutions to an auxiliary second-order elliptic boundary value problem as well as the error estimates of the associated finite element solutions. The present results can be used for analyzing DGFEM applied to some boundary value and eigenvalue problems, and also to derive the discrete Poincar´e-Friedrichs inequalities.discontinuous Galerkin FEM, polygonal element, discrete compactness, Rellich’s selection theorem

An embedded-hybridized discontinuous Galerkin method for the coupled Stokes-Darcy system

We introduce an embedded-hybridized discontinuous Galerkin (EDG-HDG) method for the coupled Stokes-Darcy system. This EDG-HDG method is a pointwise mass-conserving discretization resulting in a divergence-conforming velocity field on the whole domain. In the proposed scheme, coupling between the Stokes and Darcy domains is achieved naturally through the EDG-HDG facet variables. \emph{A priori} error analysis shows optimal convergence rates, and that the velocity error does not depend on the pressure. The error analysis is verified through numerical examples on unstructured grids for different orders of polynomial approximation

High-order DG solvers for under-resolved turbulent incompressible flows: A comparison of L2L^2 and HH(div) methods

The accurate numerical simulation of turbulent incompressible flows is a challenging topic in computational fluid dynamics. For discretisation methods to be robust in the under-resolved regime, mass conservation as well as energy stability are key ingredients to obtain robust and accurate discretisations. Recently, two approaches have been proposed in the context of high-order discontinuous Galerkin (DG) discretisations that address these aspects differently. On the one hand, standard $L^2$-based DG discretisations enforce mass conservation and energy stability weakly by the use of additional stabilisation terms. On the other hand, pointwise divergence-free $H(\operatorname{div})$-conforming approaches ensure exact mass conservation and energy stability by the use of tailored finite element function spaces. The present work raises the question whether and to which extent these two approaches are equivalent when applied to under-resolved turbulent flows. This comparative study highlights similarities and differences of these two approaches. The numerical results emphasise that both discretisation strategies are promising for under-resolved simulations of turbulent flows due to their inherent dissipation mechanisms.Comment: 24 pages, 13 figure