Calderón preconditioning of PMCHWT boundary integral equations for scattering by multiple absorbing dielectric particles

Abstract

We consider the simulation of electromagnetic scattering by single and multiple isotropic homogeneous dielectric particles using boundary integral equations. Galerkin discretizations of the classical Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) boundary integral equation formulation provide accurate solutions for complex particle geometries, but are well-known to lead to ill-conditioned linear systems. In this paper we carry out an experimental investigation into the performance of Calderón preconditioning techniques for single and multiple absorbing obstacles, which involve a squaring of the PMCHWT operator to produce a well-conditioned second-kind formulation. For single-particle scattering configurations we find that Calderón preconditioning is actually often outperformed by simple "mass-matrix" preconditioning, i.e. working with the strong form of the discretized PMCHWT operator. In the case of scattering by multiple particles we find that a significant saving in computational cost can be obtained by performing block-diagonal Calderón preconditioning in which only the self-interaction blocks are preconditioned. Using the boundary element software library Bempp (www.bempp.com) the numerical performance of the different methods is compared for a range of wavenumbers, particle geometries and complex refractive indices relevant to the scattering of light by atmospheric ice crystals.

An Introduction to Operator Preconditioning for the Fast Iterative Integral Equation Solution of Time-Harmonic Scattering Problems

“…In electromagnetism, contributions to the Calderón preconditioning introduced by Christiansen and Nédélec [52,53] for the EFIE are given e.g. in [2,3,7,8,21,67,83,97,109]. Domain decomposition preconditioners were recently proposed in [103,104] for the EFIE.…”

Helmholtz scattering by random domains: first-order sparse boundary element approximation

“…1 The induced linear systems are preconditioned by strong-form multiplicative Calderón preconditioning (cf. [24]). Tests are executed on a 32 core, 4 GB RAM per core, 64-bit Linux server using Python 2.7.6.…”

“…To finish, we introduce M the mass and A the impedance matrices. In Table 13, we plot the eigenvalues distributions of the resulting linear systems (in strong form, such as done in [24]). We remark that the spectra present some clustering at one and have similar patterns.…”

RCS Estimation of Singly Curved Dielectric Shell Structure with PMCHWT Method and Experimental Verification

“…Among them, the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) formulation has been popularly used, because of the accuracy and stability of the results at arbitrarily shaped objects [8]. In addition, the Carderon preconditioner has been used with it for the acceleration of the iterative solver [14,15].…”