2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

Yder Masson,Chao Lyu,Peter Moczo,Yann Capdeville,Barbara Romanowicz,Jean Virieux
DOI: https://doi.org/10.1093/gji/ggae025
IF: 3.352
2024-01-18
Geophysical Journal International
Abstract:SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.
geochemistry & geophysics
What problem does this paper attempt to address?