A discontinuous Galerkin method for three-dimensional poroelastic wave propagation: forward and adjoint problems

From MaRDI portal
Publication:2129488

DOI10.1007/S40315-021-00395-4zbMATH Open1484.76046arXiv2001.09478OpenAlexW3176675718WikidataQ116902704 ScholiaQ116902704MaRDI QIDQ2129488FDOQ2129488


Authors: Nick Dudley Ward, Timo Lähivaara, Simon P. Eveson Edit this on Wikidata


Publication date: 22 April 2022

Published in: Computational Methods and Function Theory (Search for Journal in Brave)

Abstract: We develop a numerical solver for three-dimensional wave propagation in coupled poroelastic-elastic media, based on a high-order discontinuous Galerkin (DG) method, with the Biot poroelastic wave equation formulated as a first order conservative velocity/strain hyperbolic system. To derive an upwind numerical flux, we find an exact solution to the Riemann problem, including the poroelastic-elastic interface; we also consider attenuation mechanisms both in Biot's low- and high-frequency regimes. Using either a low-storage explicit or implicit-explicit (IMEX) Runge-Kutta scheme, according to the stiffness of the problem, we study the convergence properties of the proposed DG scheme and verify its numerical accuracy. In the Biot low frequency case, the wave can be highly dissipative for small permeabilities; here, numerical errors associated with the dissipation terms appear to dominate those arising from discretisation of the main hyperbolic system. We then implement the adjoint method for this formulation of Biot's equation. In contrast with the usual second order formulation of the Biot equation, we are not dealing with a self-adjoint system but, with an appropriate inner product, the adjoint may be identified with a non-conservative velocity/stress formulation of the Biot equation. We derive dual fluxes for the adjoint and present a simple but illuminating example of the application of the adjoint method.


Full work available at URL: https://arxiv.org/abs/2001.09478




Recommendations




Cites Work


Cited In (2)





This page was built for publication: A discontinuous Galerkin method for three-dimensional poroelastic wave propagation: forward and adjoint problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2129488)