A non-overlapping domain decomposition method with perfectly matched layer transmission conditions for the Helmholtz equation

It is well-known that the convergence rate of non-overlapping domain decomposition methods (DDMs) applied to the parallel finite-element solution of large-scale time-harmonic wave problems strongly depends on the transmission condition enforced at the interfaces between the subdomains. Transmission operators based on perfectly matched layers (PMLs) have proved to be well-suited for configurations with layered domain partitions. They are shown to be a good compromise between basic impedance conditions, which can lead to slow convergence, and computational expensive conditions based on the exact Dirichlet-to-Neumann (DtN) map related to the complementary of the subdomain. Unfortunately, the extension of the PML-based DDM for more general partitions with cross-points (where more than two subdomains meet) is rather tricky and requires some care.