Numerical Solution of the Helmholtz Equation in 2D and 3D Using a High-Order Nyström Discretization

Author: Canino L.F.   Ottusch J.J.   Stalzer M.A.   Visher J.L.   Wandzura S.M.  

Publisher: Academic Press

ISSN: 0021-9991

Source: Journal of Computational Physics, Vol.146, Iss.2, 1998-11, pp. : 627-663

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

We show how to solve time-harmonic scattering problems by means of a high-order Nyström discretization of the boundary integral equations of wave scattering in 2D and 3D. The novel aspect of our new method is its use of local corrections to the discretized kernel in the vicinity of the kernel singularity. Enhanced by local corrections, the new algorithm has the simplicity and speed advantages of the traditional Nyström method, but also enjoys the advantages of high-order convergence for controlling solution error. We explain the practical details of implementing a scattering code based on a high-order Nyström discretization and demonstrate by numerical example that a scattering code based on this algorithm can achieve high-order convergence to the correct answer. We also demonstrate its performance advantages over a high-order Galerkin code.