Thomas F. Eibert

Simon B. Adrian, Francesco P. Andriulli, Thomas F. Eibert

We present a Calderón preconditioner for the electric field integral equation (EFIE), which does not require a barycentric refinement of the mesh and which yields a Hermitian, positive definite (HPD) system matrix allowing for the usage of the conjugate gradient (CG) solver. The resulting discrete equation system is immune to the low-frequency and the dense-discretization breakdown and, in contrast to existing Calderón preconditioners, no second discretization of the EFIE operator with Buffa-Christiansen (BC) functions is necessary. This preconditioner is obtained by leveraging on spectral equivalences between (scalar) integral operators, namely the single layer and the hypersingular operator known from electrostatics, on the one hand, and the Laplace-Beltrami operator on the other hand. Since our approach incorporates Helmholtz projectors, there is no search for global loops necessary and thus our method remains stable on multiply connected geometries. The numerical results demonstrate the effectiveness of this approach for both canonical and realistic (multi-scale) problems.

Simon B. Adrian, Francesco P. Andriulli, Thomas F. Eibert

This paper analyzes how hierarchical bases preconditioners constructed for the Electric Field Integral Equation (EFIE) can be effectively applied to the Combined Field Integral Equation (CFIE). For the case where no hierarchical solenoidal basis is available (e.g., on unstructured meshes), a new scheme is proposed: the CFIE is implicitly preconditioned on the solenoidal Helmholtz subspace by using a Helmholtz projector, while a hierarchical non-solenoidal basis is used for the non-solenoidal Helmholtz subspace. This results in a well-conditioned system. Numerical results corroborate the presented theory.