1. Exterior Scattering Problem

We consider the standard exterior Helmholtz scattering problem. Let \(\mathcal O\subset\mathbb R^d\), \(d=2,3\), be a bounded scatterer with boundary \(\Gamma=\partial\mathcal O\), and let \(D = \mathbb R^d\setminus \overline{\mathcal O}\). For a sound-soft obstacle, we write the total field as \(U^{\rm tot}=U^{\rm inc}+U\), where \(U^{\rm inc}\) is a prescribed incident wave and \(U\) is the scattered field. In the simplest constant-index case,

\[ (\Delta+k^2)U=0 \qquad \text{in }D, \]

with the boundary condition

\[ U=-U^{\rm inc}\qquad \text{on }\Gamma. \]

For a plane wave incident from direction \(a\in\mathbb S^{d-1}\),

\[ U^{\rm inc}(x)=e^{ik a\cdot x}. \]

The scattering problem requires a boundary condition at infinity. The scattered field must be outgoing, which is expressed by the Sommerfeld radiation condition

\[ \lim_{r\to\infty} r^{\frac{d-1}{2}} \left(\partial_r U - ikU\right)=0, \qquad r=|x|. \]

An outgoing scattered field has the asymptotic expansion

\[ U(r,\omega) =r^{-\frac{d-1}{2}}e^{ikr}U_\infty(\omega) +\mathcal O\left(r^{-\frac{d+1}{2}}\right), \qquad \omega=\frac{x}{|x|}. \]

The function \(U_\infty\) is the far-field pattern. In conventional finite-element calculations, one solves for the scattered field \(U(r,\omega)\), and the domain is truncated at a finite radius. The artificially introduced boundary requires an approximate absorbing boundary condition, a Dirichlet-to-Neumann map, or a perfectly matched layer.

Our approach is to directly solve the unbounded problem for a suitably extended far-field variable by compactifying the exterior domain to a finite computational domain. The numerical solution gives both the near field around \(\mathcal O\) and \(U_\infty\) at the outer mesh boundary.

To put our approach in context, the main exterior-domain strategies differ in where they place the boundary and how they select the outgoing solution:

method	computational exterior	outgoing condition	access to \(U_\infty\)
absorbing boundary condition	truncated at a finite boundary	local approximation on that boundary	postprocessing from finite-radius data
perfectly matched layer	truncated after a complex absorbing layer	decay through coordinate stretching plus an outer termination	postprocessing from the interior physical region
Dirichlet-to-Neumann map	truncated at a finite boundary	exact or truncated nonlocal boundary operator	available from the boundary expansion
infinite elements	unbounded elements attached to a finite interface	outgoing asymptotics built into the exterior approximation space	available from exterior expansion coefficients
hyperboloidal compactification	infinity is a finite mesh boundary	incoming fields are excluded by regularity of the transformed unknown	boundary trace of the transformed solution