On Long-Time Instabilities in Staggered Finite Difference Simulations of the Seismic Acoustic Wave Equations on Discontinuous Grids
KAUST DepartmentApplied Mathematics and Computational Science Program
Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
Extreme Computing Research Center
KAUST Grant NumberCCF/CAF/URF/1/2596
Online Publication Date2017-10-27
Print Publication Date2018-02-01
Permanent link to this recordhttp://hdl.handle.net/10754/626096
MetadataShow full item record
AbstractWe consider the long-time instability issue associated with finite difference simulation of seismic acoustic wave equations on discontinuous grids. This issue is exhibited by a prototype algebraic problem abstracted from practical application settings. Analysis of this algebraic problem leads to better understanding of the cause of the instability and provides guidance for its treatment. Specifically, we use the concept of discrete energy to derive the proper solution transfer operators and design an effective way to damp the unstable solution modes. Our investigation shows that the interpolation operators need to be matched with their companion restriction operators in order to properly couple the coarse and fine grids. Moreover, to provide effective damping, specially designed diffusive terms are introduced to the equations at designated locations and discretized with specially designed schemes. These techniques are applied to simulations in practical settings and are shown to lead to superior results in terms of both stability and accuracy.
CitationGao L, Ketcheson D, Keyes D (2017) On Long-Time Instabilities in Staggered Finite Difference Simulations of the Seismic Acoustic Wave Equations on Discontinuous Grids. Geophysical Journal International. Available: http://dx.doi.org/10.1093/gji/ggx470.
SponsorsThe authors acknowledge the support of the King Abdullah University of Science and Technology through the project CCF/CAF/URF/1/2596.
PublisherOxford University Press (OUP)