GPU-accelerated discontinuous Galerkin methods on hybrid meshes
august, 2016
Publication type:
Paper in peer-reviewed journals
Journal:
Journal of Computational Physics, vol. 318, pp. 142 - 168
Publisher:
Elsevier
HAL:
arXiv:
Keywords :
Discontinuous Galerkin; GPU; High order; Hybrid mesh; Timestep restriction; Wave equation;
Abstract:
We present a time-explicit discontinuous Galerkin (DG) solver for the time-domain acoustic wave equation on hybrid meshes containing vertex-mapped hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable formulations are presented for both Gauss–Legendre and Gauss–Legendre–Lobatto (Spectral Element) nodal bases for the hexahedron. Stable timestep restrictions for hybrid meshes are derived by bounding the spectral radius of the DG operator using order-dependent constants in trace and Markov inequalities. Computational efficiency is achieved under a combination of element-specific kernels (including new quadrature-free operators for the pyramid), multi-rate timestepping, and acceleration using Graphics Processing Units.
BibTeX:
@article{Cha-Wan-Mod-Rem-War-2016, author={Jesse Chan and Zheng Wang and Axel Modave and Jean-Francois Remacle and Tim Warburton }, title={GPU-accelerated discontinuous Galerkin methods on hybrid meshes }, doi={10.1016/j.jcp.2016.04.003 }, journal={Journal of Computational Physics }, year={2016 }, month={8}, volume={318 }, pages={142--168}, }