A Multilevel Approach for Trace System in HDG Discretizations
arXiv:1903.11045 · doi:10.1016/j.jcp.2020.109240
Abstract
We propose a multilevel approach for trace systems resulting from hybridized discontinuous Galerkin (HDG) methods. The key is to blend ideas from nested dissection, domain decomposition, and high-order characteristic of HDG discretizations. Specifically, we first create a coarse solver by eliminating and/or limiting the front growth in nested dissection. This is accomplished by projecting the trace data into a sequence of same or high-order polynomials on a set of increasingly coarser edges/faces. We then combine the coarse solver with a block-Jacobi fine scale solver to form a two-level solver/preconditioner. Numerical experiments indicate that the performance of the resulting two-level solver/preconditioner depends only on the smoothness of the solution and is independent of the nature of the PDE under consideration. While the proposed algorithms are developed within the HDG framework, they are applicable to other hybrid(ized) high-order finite element methods. Moreover, we show that our multilevel algorithms can be interpreted as a multigrid method with specific intergrid transfer and smoothing operators. With several numerical examples from Poisson, pure transport, and convection-diffusion equations we demonstrate the robustness and scalability of the algorithms.
References in corpus (6)
- A hybridizable discontinuous Galerkin method for the Navier--Stokes equations with pointwise divergence-free velocity field
- A performance comparison of continuous and discontinuous Galerkin methods with fast multigrid solvers
- Analysis of a hybridized/interface stabilized finite element method for the Stokes equations
- A Robust Multilevel Method for Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation
- Analysis of Schwarz methods for a hybridizable discontinuous Galerkin discretization
- An Improved Iterative HDG Approach for Partial Differential Equations
Cited by in corpus (6)
- HDGlab: An open-source implementation of the hybridisable discontinuous Galerkin method in MATLAB
- A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
- Preconditioning for a pressure-robust HDG discretization of the Stokes equations
- An unfitted high-order HDG method for two-fluid Stokes flow with exact NURBS geometries
- Uniform auxiliary space preconditioning for HDG methods for elliptic operators with a parameter dependent low order term
- Multigrid preconditioners for the hybridized Discontinuous Galerkin discretisation of the shallow water equations