Numerical Simulation of Isotropic 2D Elastic Waves in Homogeneous Media

Numerical simulation of 2D elastic waves in homogeneous isotropic media represents a crucial geophysical technique that enables better understanding of Earth's internal structure and material composition. This implementation employs Cerjan boundary conditions to confine wave propagation domains, ensuring computational accuracy. The code incorporates velocity models and Ricker wavelet functions to effectively simulate elastic wave propagation through Earth's interior. The simulation utilizes a finite-difference time-domain (FDTD) method with second-order spatial derivatives and second-order time integration. The Cerjan boundary condition implementation applies exponential damping coefficients along the grid edges to minimize artificial reflections. Before executing the 2D elastic wave numerical simulation, users need to extract the compressed files and install required software dependencies. The package includes pre-configured velocity models (homogeneous media properties) and Ricker wavelet source functions defined by: f(t) = (1 - 2π²f₀²t²)exp(-π²f₀²t²) where f₀ represents the dominant frequency. The simulation can be run directly after setup completion. The algorithm's precision and stability have been extensively validated, with applications spanning seismic exploration and geological research. For deeper understanding of 2D elastic wave numerical simulation principles and applications, refer to the following literature: - E. Haber, and I. Stadler, "A multigrid method for elliptic equations with highly oscillatory coefficients," Journal of Computational Physics, vol. 215, no. 2, pp. 631-652, 2006. - S. Zhang, and Y. Chen, "A stable and efficient finite-difference method for the seismic wave equation in second-order form," Geophysical Journal International, vol. 192, no. 2, pp. 901-918, 2013. We hope this information proves valuable for your research and implementation needs!