An efficient multigrid calculation of the far field map for Helmholtz and Schr\"odinger equations

TL;DR

This paper introduces a highly efficient multigrid-based method for calculating the far field pattern in scattering problems governed by Helmholtz and Schrödinger equations, utilizing a complex domain reformulation for improved scalability.

Contribution

It presents a novel complex contour reformulation combined with multigrid techniques for fast, scalable far field calculations in high-dimensional scattering problems.

Findings

Efficient computation of far field patterns in 2D and 3D.

Validation of method on Helmholtz and Schrödinger problems.

Demonstrated scalability with respect to wavenumber.

Abstract

In this paper we present a new highly efficient calculation method for the far field amplitude pattern that arises from scattering problems governed by the d-dimensional Helmholtz equation and, by extension, Schr\"odinger's equation. The new technique is based upon a reformulation of the classical real-valued Green's function integral for the far field amplitude to an equivalent integral over a complex domain. It is shown that the scattered wave, which is essential for the calculation of the far field integral, can be computed very efficiently along this complex contour (or manifold, in multiple dimensions). Using the iterative multigrid method as a solver for the discretized damped scattered wave system, the proposed approach results in a fast and scalable calculation method for the far field map. The complex contour method is successfully validated on Helmholtz and Schr\"odinger model problems in two and three spatial dimensions, and multigrid convergence results are provided to substantiate the wavenumber scalability and overall performance of the method.