Authors: William W Symes Tetyana Vdovina
Publish Date: 2009/02/21
Volume: 13, Issue: 3, Pages: 363-371
Abstract
The numerical error associated with finitedifference simulation of wave propagation in discontinuous media consists of two components The first component is a higherorder error that leads to grid dispersion it can be controlled by higherorder methods The second component results from misalignment between numerical grids and material interfaces We provide an explicit estimate of the interface misalignment error for the second order in time and space staggered finitedifference scheme applied to the acoustic wave equation Our analysis confirmed by numerical experiments demonstrates that the interface error results in a firstorder time shift proportional to the distance between the interface and computational grids A 2D experiment shows that the interface error cannot be suppressed by higherorder methods and indicates that our 1D analysis gives a good prediction about the behavior of the numerical solution in higher dimensions
Keywords: