Numerical dispersion analysis for three-dimensional Laplace-Fourier-domain scalar wave equation

摘要

Based on the phase velocity and attenuation propagation velocity, a method for performing numerical dispersion analysis of three-dimensional Laplace-Fourier-domain scalar wave equation is presented. This method is applied to a 27-point average-derivative optimal scheme and a 27-point finite-element scheme. Within the relative error of 1%, the 27-point average-derivative optimal scheme requires seven grid points per wavelength and pseudo-wavelength while the 27-point finite-element scheme requires 23 grid points per wavelength and pseudo-wavelength for equal and unequal directional sampling intervals. Numerical examples show that the 27-point Laplace-Fourier-domain average-derivative optimal scheme is more accurate than the 27-point Laplace-Fourier-domain finite-element scheme for the same computational cost. By using larger directional sampling intervals while maintaining accuracy, the 27-point Laplace-Fourier-domain average-derivative optimal scheme can greatly reduce the computational cost of three-dimensional Laplace-Fourier-domain modelling.