Simple elastic finite-difference scheme

Jiri Zahradnik


This paper investigates which features of the elastic finite-difference schemes are essential for their accuracy, and which ones allow simplifications. It is shown that the schemes employing the geometrically averaged parameters are more accurate than those using local material parameters, mainly when a discontinuity passes between the grid lines. It is also shown that the accuracy of the mixed spatial derivatives at the internal grid points does not degrade when the number of the implicitly employed stress values and the geometrically averaged material parameters decreases from four to two (the co-called Full and Short Forms, respectively). The Short and Full Forms give the same numerical results, while 50 percent of the arithmetic operations are saved with the Short one. However, at the free-surface points such an simplification is not permitted, and the Full Form should be used. Based on these results a new simple elastic scheme (called PS2) is suggested.

Bull. Seism. Soc. Am., 85 (1995), 1879-1887.
SW3D - main page of consortium Seismic Waves in Complex 3-D Structures .