Seismic modeling with RBF-FD – a simplified treatment of interfaces Public Deposited

Downloadable Content

Download PDF
  • In a previous study of seismic modeling with radial basis function-generated finite differences (RBF-FD), we outlined a numerical method for solving 2-D wave equations in domains with material interfaces between different regions. The method was applicable on a mesh-free set of data nodes. It included all information about interfaces within the weights of the stencils (allowing the use of traditional time integrators), and was shown to solve problems of the 2-D elastic wave equation to 3rd-order accuracy. In the present paper, we discuss a refinement of that method that makes it simpler to implement. It can also improve accuracy for the case of smoothly-variable model parameter values near interfaces. We give several test cases that demonstrate the method solving 2-D elastic wave equation problems to 4th-order accuracy, even in the presence of smoothly-curved interfaces with jump discontinuities in the model parameters.
Date Issued
  • 2017-01-01
Additional Information
  • This is a post-print version of an article published in the Journal of Computational Physics. - Journal of Computational Physics
Academic Affiliation
Journal Title
Journal Volume
  • 335.0
Last Modified
  • 2020-01-06
Resource Type
Rights Statement
  • 1090-2716