In this paper, we consider a high-order discontinuous Galerkin (DG) method for modelling wave propagation in coupled poroelastic–elastic media. The upwind numerical flux is derived as an exact solution for the Riemann problem including the poroelastic–elastic interface. Attenuation mechanisms in both Biot's low- and high-frequency regimes are considered. The current implementation supports non-uniform basis orders which can be used to control the numerical accuracy element by element. In the numerical examples, we study the convergence properties of the proposed DG scheme and provide experiments where the numerical accuracy of the scheme under consideration is compared to analytic and other numerical solutions.
Bibliographical note© 2017 Elsevier B.V. or its licensors or contributors. This is an author-produced version of the published paper. Uploaded in accordance with the publisher’s self-archiving policy.
- Discontinuous Galerkin method
- Non-uniform basis order
- Poroelastic waves