Communications on Applied Mathematics and Computation ›› 2022, Vol. 4 ›› Issue (3): 904-944.doi: 10.1007/s42967-021-00155-0

• ORIGINAL PAPER • Previous Articles     Next Articles

A Space-Time Interior Penalty Discontinuous Galerkin Method for the Wave Equation

Poorvi Shukla, J. J. W. van der Vegt   

  1. Department of Applied Mathematics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
  • Received:2020-12-20 Revised:2021-06-11 Online:2022-09-20 Published:2022-07-04
  • Contact: J. J. W. van der Vegt,E-mail:j.j.w.vandervegt@utwente.nl;Poorvi Shukla,E-mail:p.shukla@utwente.nl E-mail:j.j.w.vandervegt@utwente.nl;p.shukla@utwente.nl
  • Supported by:
    This research was funded by the Shell-NWO Computational Sciences for Energy Research Program.

Abstract: A new higher-order accurate space-time discontinuous Galerkin (DG) method using the interior penalty flux and discontinuous basis functions, both in space and in time, is presented and fully analyzed for the second-order scalar wave equation. Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method. The theoretical analysis shows that the DG discretization is stable and converges in a DG-norm on general unstructured and locally refined meshes, including local refinement in time. The space-time interior penalty DG discretization does not have a CFL-type restriction for stability. Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time step $ \Delta t $ satisfy $ h\cong C\Delta t $, with C a positive constant. The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems. These calculations also show that for pth-order tensor product basis functions the convergence rate in the $ L^\infty $ and $ L^2 $-norms is order $ p+1 $ for polynomial orders $ p=1 $ and $ p=3 $ and order p for polynomial order $ p=2 $.

Key words: Wave equation, Space-time methods, Discontinuous Galerkin methods, Interior penalty method, A priori error analysis

CLC Number: