Previously, we have applied the method of difference potentials to solving several variable coefficient interior Helmholtz problems with fourth and sixth order accuracy. We have employed compact finite difference schemes as a core discretization methodology. Those schemes enable high order accuracy on narrow stencils and hence require only as many boundary conditions as needed for the underlying differential equation itself. Numerical experiments corroborate the high order accuracy of our method for variable coefficients, regular grids, and non-conforming boundaries.
In the current paper, we extend the previously developed methodology to exterior problems. We present a complete theoretical analysis of the algorithm, as well as the results of a series of numerical simulations. Specifically, we study the scattering of time-harmonic waves about smooth shapes, subject to various boundary conditions. We also solve the transmission/scattering problems, in which not only do the waves scatter off a given shape but also propagate through the interface and travel across the heterogeneous medium inside. In all the cases, our methodology guarantees high order accuracy for regular grids and non-conforming boundaries and interfaces.