Galerkin–Legendre Spectral Method for the 3D Helmholtz Equation

Author: Auteri F.   Quartapelle L.  

Publisher: Academic Press

ISSN: 0021-9991

Source: Journal of Computational Physics, Vol.161, Iss.2, 2000-07, pp. : 454-483

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

A Galerkin–Legendre spectral method for the direct solution of Poisson and Helmholtz equations in a three-dimensional rectangular domain is presented. The method extends Jie Shen's algorithm for 2D problems by using the diagonalization of the three mass matrices in the three spatial directions and fully exploits the direct product nature of the spectral approximation. The Dirichlet boundary values are taken into account by means of a discrete lifting performed in three subsequent steps and built upon Gauss–Legendre quadrature points. A few numerical tests illustrate the accuracy and efficiency of the method.