DFT MODAL ANALYSIS OF SPECTRAL ELEMENT METHODS FOR ACOUSTIC WAVE PROPAGATION

Abstract

Spectral element methods are now widely used for wave propagation simulations. They are appreciated for their high order of accuracy, but are still used on a heuristic basis. In this work we present the numerical dispersion of spectral elements, which allows us to assess their error limits and to devise efficient numerical simulations, particularly for 2D and 3D problems. We propose a novel approach based on a discrete Fourier transform of both the probing plane waves and the discrete wave operators. The underlying dispersion relation is estimated by the Rayleigh quotients of the plane waves with respect to the discrete operator. Together with the Kronecker product properties, this approach delivers numerical dispersion estimates for 1D operators as well as for 2D and 3D operators, and is well suited for spectral element methods, which use nonequidistant collocation points such as Gauss–Lobatto–Chebyshev and Gauss–Lobatto–Legendre points. We illustrate this methodology with dispersion and anisotropy graphs for spectral elements with polynomial degrees from 4 to 12. These graphs confirm the rule of thumb that spectral element methods reach a safe level of accuracy at about four grid points per wavelength.