A nodal discontinuous Galerkin finite element method for the poroelastic wave equation
We use the nodal discontinuous Galerkin method with a Lax-Friedrich flux to model the wave propagation in transversely isotropic and poroelastic media. The effect of dissipation due to global fluid flow causes a stiff relaxation term, which is incorporated in the numerical scheme through an operator splitting approach. The well-posedness of the poroelastic system is proved by adopting an approach based on characteristic variables. An error analysis for a plane wave propagating in poroelastic media shows a convergence rate of O(hn+1). Computational experiments are shown for various combinations of homogeneous and heterogeneous poroelastic media.
NDG KShukla 2018.pdf
openaccess
4.13 MB
Adobe PDF
95b5b047b90583005e5ebc581c650058