A high-order Discontinuous Galerkin method for the seismic wave propagation
Résumé
We are interested in the simulation of P-SV seismic wave propagation by a high-order Discontinuous Galerkin method based on centered fluxes at the interfaces combined with a leap-frog time-integration. This non-diffusive method, previously developed for the Maxwell equations, is particularly well adapted to complex topographies and fault discontinuities in the medium. We prove that the scheme is stable under a CFL type condition and that a discrete energy is preserved on an infinite domain. Convergence properties and efficiency of the method are studied through numerical simulations in two and three dimensions of space.