Fast Converging Single Trace Quasi-local PMCHWT Equation for the Modelling of Composite Systems

TL;DR

This paper introduces a new single trace quasi-local PMCHWT integral equation that models composite systems with junctions efficiently, converges rapidly, and avoids interior resonances, improving computational performance in electromagnetic scattering problems.

Contribution

A novel single trace quasi-local PMCHWT equation is developed, enabling faster convergence and better handling of junctions in composite systems compared to existing methods.

Findings

Converges with a slow increase in iterations as mesh size decreases

Demonstrates correctness and convergence through numerical experiments

Free from interior resonances, enhancing solution stability

Abstract

The PMCHWT integral equation enables the modelling of scattering of time-harmonic fields by penetrable, piecewise homogeneous, systems. They have been generalised to include the modelling of composite systems that may contain junctions, i.e. lines along which three or more materials meet. Linear systems resulting upon discretisation of the PMCHWT are, because of their large dimension, typically solved by Krylov iterative methods. The number of iterations required for this solution critically depends on the eigenvalue distribution of the system matrix. For systems that do not contain junction lines, Calder\'on preconditioning, which was first applied to the electric field integral equation, has been generalised to the PMCHWT equation. When junctions are present, this approach cannot be applied. Alternative approaches, such as the global multi-trace method, conceptually remove the junction lines and as a result are amenable to Calder\'on preconditioning. This approach entails a doubling of the degrees of freedom, and the solution that is produced only approximately fulfils the continuity conditions at interfaces separating domains. In this contribution, a single trace quasi-local PMCHWT equation is introduced that requires a number of iterations for its solution that only slowly increases as the mesh size tends to zero. The method is constructed as a generalisation of the classic PMCHWT, and its discretisation is thoroughly discussed. A comprehensive suite of numerical experiments demonstrates the correctness, convergence behaviour, and efficiency of the method. The integral equation is demonstrated to be free from interior resonances.