A spectral boundary integral equation method for the 2D Helmholtz equation
- Old Dominion Univ., Norfolk, VA (United States)
In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the non-smoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the non-smoothness of the integral kernels in the spectral implementation. The present method is robust for a general smooth boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. Numerical examples of wave scattering are given in which the exponential accuracy of the present numerical method is demonstrated. 15 refs., 3 figs., 4 tabs.
- OSTI ID:
- 110897
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 2 Vol. 120; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Galerkin Boundary Integral Analysis for the 3D Helmholtz Equation
A new FFT-based algorithm to compute Born radii in the generalized Born theory of biomolecule solvation