Because of its high computational cost, we needed to develop an efficient numerical scheme for the frequency-domain 3D elastic wave equation. In addition, the numerical scheme should be applicable to media with a liquid-solid interface. To address these two issues, we have developed a new average-derivative optimal 27-point scheme with arbitrary directional grid intervals and a corresponding numerical dispersion analysis for the frequency-domain 3D elastic wave equation. The novelty of this scheme is that its optimal coefficients depend on the ratio of the directional grid intervals and Poisson’s ratio. In this way, this scheme can be applied to media with a liquid-solid interface and a computational grid with arbitrary directional grid intervals. For media with a variable Poisson’s ratio, we have developed an effective and stable interpolation method for optimization coefficients. Compared with the classic 19-point scheme, this new scheme reduces the required number of grid points per wavelength for equal and unequal directional grid intervals. The reduction of the number of grid points increases as the Poisson’s ratio becomes larger. In particular, the numerical S-wave phase velocity of this new scheme becomes zero, whereas the classic 19-point scheme produces a spurious numerical S-wave phase velocity, as Poisson’s ratio reaches 0.5. We have performed numerical examples to develop the theoretical analysis.
Mode-conversion of the scattering kernel for the elastic wave equation
