Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Details and results from an extension to heterogeneous media are also included. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the proposed method. A novel local high order ABC for elastic waves based on farfield expansions is constructed and preliminary results applying it to elastic scattering problems are presented.
College and Department
Physical and Mathematical Sciences
BYU ScholarsArchive Citation
Grundvig, Dane Scott, "High Order Numerical Methods for Problems in Wave Scattering" (2020). Theses and Dissertations. 8617.
acoustic scattering, elastic scattering, Helmholtz equation, high order numerical methods, variable wave number, heterogeneous media, deferred corrections