GPU-accelerated discontinuous Galerkin methods on hybrid meshes

JOURNAL OF COMPUTATIONAL PHYSICS

卷:318

Article

Chan, Jesse¹  Wang, Zheng¹  Modave, Axel¹  Remacle, Jean-Francois²  Warburton, T.¹ 

[1] Rice Univ, Dept Computat & Appl Math, 6100 Main St, Houston, TX 77005 USA

[2] Catholic Univ Louvain, Sch Engn, Ave Georges Lemaitre 4, B-1348 Louvain, Belgium

关键词: Discontinuous Galerkin;    GPU;    High order;    Hybrid mesh;    Timestep restriction;    Wave equation;   

DOI  :  10.1016/j.jcp.2016.04.003

来源: Elsevier

【 摘 要 】

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. (C) 2016 Published by Elsevier Inc.

【 授权许可】

Free