Weight‐adjusted discontinuous Galerkin methods: Matrix‐valued weights and elastic wave propagation in heterogeneous media. (28th November 2017)

• Record Type: Journal Article
• Authors: Chan, Jesse
• Is Part Of: International journal for numerical methods in engineering. Volume 113:Number 12(2018)
• Page Start: 1779
• Page End: 1809
• Publication Date: 2017-11-28
• Subjects: Numerical analysis -- Periodicals; Engineering mathematics -- Periodicals; 620.001518
• DOI: 10.1002/nme.5720
• Languages: English
• ISSNs: 0029-5981

Abstract

Summary: Weight‐adjusted inner products are easily invertible approximations to weighted L 2 inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time‐domain method for wave propagation which is low storage, energy stable, and high‐order accurate for arbitrary heterogeneous media and curvilinear meshes. In this work, we extend weight‐adjusted DG methods to the case of matrix‐valued weights, with the linear elastic wave equation as an application. We present a DG formulation of the symmetric form of the linear elastic wave equation, with upwind‐like dissipation incorporated through simple penalty fluxes. A semidiscrete convergence analysis is given, and numerical results confirm the stability and high‐order accuracy of weight‐adjusted DG for several problems in elastic wave propagation.