Many problems that involve the propagation of time-harmonic waves are naturally posed in unbounded domains. For instance, a common problem in the are a of acoustic scattering is the determination of the sound field that is generated when an incoming time-harmonic wave (which is assumed to arrive ``from infinity'') impinges onto a solid body (the scatterer). The boundary
conditions to be applied on the surface of the scatterer (most often of Dirichlet, Neumann or Robin type) tend to be easy to enforce in most numerical solution schemes. Conversely, the imposition of a suitable decay condition (typically a variant of the Sommerfeld radiation condition), which is required to ensure the well-posedness of the solution, is considerably more involved. As a result, many numerical schemes generate spurious reflections from the outer boundary of the finite computational domain.
Perfectly matched layers (PMLs) are in this context a versatile alternative to the usage of classical approaches such as employing absorbing boundary conditions or Dirichlet-to-Neumann mappings, but unfortunately most PML formulations contain adjustable parameters which have to be optimised to give the best possible performance for a particular problem. In this talk I will present a parameter-free PML formulation for the case of the two-dimensional Helmholtz equation. The performance of the proposed method is demonstrated via extensive numerical experiments, involving domains with smooth and polygonal boundaries, different solution types (smooth and singular, planar and non-planar waves), and a wide range of wavenumbers (R. Cimpeanu, A. Martinsson and M.Heil, J. Comp. Phys., 296, 329-347 (2015)). Possible extensions and generalisations will also be touched upon.