Iterative Solution of the Helmholtz Equation By a Second-Order Method

Loading...
Thumbnail Image

Files

CS-TR-3727.ps (1.06 MB)
No. of downloads: 341
CS-TR-3727.pdf (1010.98 KB)
No. of downloads: 2153

Publication or External Link

Date

1998-10-15

Advisor

Citation

DRUM DOI

Abstract

The numerical solution of the Helmholtz equation subject to nonlocal radiation boundary conditions is studied. The specific problem is discretized with a second-order accurate finite-difference method, resulting in a linear system of equations. To solve the system of equations, a preconditioned Krylov subspace method is employed. The preconditioner is based on fast transforms, and yields a direct fast Helmholtz solver for rectangulay domains. Numerical experiments for curved ducts demonstrate that the rate of convergence is high. Compared with band Gaussian elimination the preconditioned iterative method shows a significant gain in both storage requirement and arithmetic complexity.

(Also cross-referenced as UMIACS-TR-96-95)

Notes

Rights