2019/03/15 by Longfei Gao, David E. Keyes, Gao, Longfei +1
Earth and Planetary Sciences · Engineering · #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #FOS: Mathematics #Numerical Analysis (math.NA) #Seismic Imaging and Inversion Techniques
paper · pdf · doi:10.48550/arxiv.1903.06795
openalex publication_date 2019/03/15 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We consider the finite difference discretization of isotropic elastic wave equations on nonuniform grids. The intended applications are seismic studies, where heterogeneity of the earth media can lead to severe oversampling for simulations on uniform grids. To address this issue, we demonstrate how to properly couple two non-overlapping neighboring subdomains that are discretized uniformly, but with different grid spacings. Specifically, a numerical procedure is presented to impose the interface conditions weakly through carefully designed penalty terms, such that the overall semi-discretization conserves a discrete energy resembling the continuous energy possessed by the elastic wave system.