Hybrid multigrid methods for high-order discontinuous Galerkin discretizations

Abstract

The present work develops hybrid multigrid methods for high-order discontinuous Galerkin discretizations of elliptic problems. Fast matrix-free operator evaluation on tensor product elements is used to devise a computationally efficient PDE solver. The multigrid hierarchy exploits all possibilities of geometric, polynomial, and algebraic coarsening, targeting engineering applications on complex geometries. Additionally, a transfer from discontinuous to continuous function spaces is performed within the multigrid hierarchy. This does not only further reduce the problem size of the coarse-grid problem, but also leads to a discretization most suitable for state-of-the-art algebraic multigrid methods applied as coarse-grid solver. The relevant design choices regarding the selection of optimal multigrid coarsening strategies among the various possibilities are discussed with the metric of computational costs as the driving force for algorithmic selections. We find that a transfer to a continuous function space at highest polynomial degree (or on the finest mesh), followed by polynomial and geometric coarsening, shows the best overall performance. The success of this particular multigrid strategy is due to a significant reduction in iteration counts as compared to a transfer from discontinuous to continuous function spaces at lowest polynomial degree (or on the coarsest mesh). The coarsening strategy with transfer to a continuous function space on the finest level leads to a multigrid algorithm that is robust with respect to the penalty parameter of the SIPG method. Detailed numerical investigations are conducted for a series of examples ranging from academic test cases to more complex, practically relevant geometries. Performance comparisons to state-of-the-art methods from the literature demonstrate the versatility and computational efficiency of the proposed multigrid algorithms.

References in corpus (15)

• Fast matrix-free evaluation of discontinuous Galerkin finite element operators
• A performance comparison of continuous and discontinuous Galerkin methods with fast multigrid solvers
• A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow
• FFT, FMM, or Multigrid? A comparative Study of State-Of-the-Art Poisson Solvers for Uniform and Nonuniform Grids in the Unit Cube
• On the Strong Scaling of the Spectral Element Solver Nek5000 on Petascale Systems
• Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods
• Efficiency of high-performance discontinuous Galerkin spectral element methods for under-resolved turbulent incompressible flows
• A p-Multigrid Strategy with Anisotropic p-Adaptation Based on Truncation Errors for High-Order Discontinuous Galerkin Methods
• On thin plate spline interpolation
• Matrix-free multigrid block-preconditioners for higher order Discontinuous Galerkin discretisations
• GPU accelerated spectral finite elements on all-hex meshes
• High-performance Implementation of Matrix-free High-order Discontinuous Galerkin Methods
• Robust multigrid for high-order discontinuous Galerkin methods: A fast Poisson solver suitable for high-aspect ratio Cartesian grids
• Scaling to the stars -- a linearly scaling elliptic solver for $p$-multigrid
• FFT, FMM, and Multigrid on the Road to Exascale: performance challenges and opportunities

Cited by in corpus (13)

• A Flexible, Parallel, Adaptive Geometric Multigrid method for FEM
• A matrix-free high-order solver for the numerical solution of cardiac electrophysiology
• High-order arbitrary Lagrangian-Eulerian discontinuous Galerkin methods for the incompressible Navier-Stokes equations
• Algorithms for Parallel Generic $hp$-adaptive Finite Element Software
• Matrix-Free Higher-Order Finite Element Methods for Hyperelasticity
• Smoothers with localized residual computations for geometric multigrid methods
• Geometric Uncertainty of Patient-Specific Blood Vessels and its Impact on Aortic Hemodynamics
• A discontinuous Galerkin discretization of elliptic problems with improved convergence properties using summation by parts operators
• SynthAorta: A 3D Mesh Dataset of Parametrized Physiological Healthy Aortas
• An hp Multigrid Approach for Tensor-Product Space-Time Finite Element Discretizations of the Stokes Equations
• Deep Medial Voxels: Learned Medial Axis Approximations for Anatomical Shape Modeling
• On efficient multigrid solvers for mixed-degree local discontinuous Galerkin multiphase Stokes problems
• Higher-Order Discontinuous Galerkin Splitting Schemes for Fluids with Variable Viscosity