A weight-adjusted discontinuous Galerkin method for the poroelastic wave equation: penalty fluxes and micro-heterogeneities

Abstract

We introduce a high-order weight-adjusted discontinuous Galerkin (WADG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in anisotropic porous media. We use a coupled first-order symmetric stress-velocity formulation. Careful attention is directed at (a) the derivation of an energy-stable penalty-based numerical flux, which offers high-order accuracy in presence of material discontinuities, and (b) proper treatment of micro-heterogeneities (sub-element variations) in the numerical scheme. The use of a penalty-based numerical flux avoids the diagonalization of Jacobian matrices into polarized wave constituents necessary when solving element-wise Riemann problems. Micro-heterogeneities are accurately and stably incorporated in the numerical scheme using easily-invertible weight-adjusted mass matrices. The convergence of the proposed numerical scheme is proven and verified by using convergence studies against analytical plane wave solutions. The proposed method is also compared against an existing implementation using the spectral element method to solve the poroelastic wave equation.

References in corpus (5)

• Nodal Discontinuous Galerkin Methods on Graphics Processors
• GPU-accelerated discontinuous Galerkin methods on hybrid meshes
• A Discontinuous Galerkin method with a modified penalty flux for the propagation and scattering of acousto-elastic waves
• Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media
• Bernstein-Bezier weight-adjusted discontinuous Galerkin methods for wave propagation in heterogeneous media

Cited by in corpus (3)

• A Memory-efficient Implementation of Perfectly Matched Layer with Smoothly-varying Coefficients in Discontinuous Galerkin Time-Domain Method
• A weight-adjusted discontinuous Galerkin method for wave propagation in coupled elastic-acoustic media
• A Discontinuous Galerkin method for three-dimensional elastic and poroelastic wave propagation: forward and adjoint problems