Preconditioned iterative methods for inhomogeneous acoustic scattering applications
Doctor of Philosophy
This thesis develops and analyzes efficient iterative methods for solving discretizations of the Lippmann--Schwinger integral equation for inhomogeneous acoustic scattering. Analysis and numerical illustrations of the spectral properties of the scattering problem demonstrate that a significant portion of the spectrum is approximated well on coarse grids. To exploit this, I develop a novel restarted GMRES method with adaptive deflation preconditioning based on spectral approximations on multiple grids. Much of the literature in this field is based on exact deflation, which is not feasible for most practical computations. This thesis provides an analytical framework for general approximate deflation methods and suggests a way to rigorously study a host of inexactly-applied preconditioners. Approximate deflation algorithms are implemented for scattering through thin inhomogeneities in photonic band gap problems. I also develop a short term recurrence for solving the one dimensional version of the problem that exploits the observation that the integral operator is a low rank perturbation of a self-adjoint operator. This method is based on strategies for solving Schur complement problems, and provides an alternative to a recent short term recurrence algorithm for matrices with such structure that we show to be numerically unstable for this application. The restarted GMRES method with adaptive deflation preconditioning over multiple grids, as well as the short term recurrence method for operators with low rank skew-adjoint parts, are very effective for reducing both the computational time and computer memory required to solve acoustic scattering problems. Furthermore, the methods are sufficiently general to be applicable to a wide class of problems.
Applied mathematics; Physics; Acoustics