Second kind boundary integral equation for multi-subdomain diffusion problems

We consider isotropic scalar diffusion boundary value problems whose diffusion coefficients are piecewise constant with respect to a partition of space into Lipschitz subdomains. We allow so-called material junctions where three or more subdomains may abut. We derive a boundary integral equation of the second kind posed on the skeleton of the subdomain partition that involves, as unknown, only one trace function at each point of each interface. We prove the well-posedness of the corresponding boundary integral equations. We also report numerical tests for Galerkin boundary element discretisations, in which the new approach proves to be highly competitive compared to the well-established first kind direct single-trace boundary integral formulation. In particular, GMRES seems to enjoy fast convergence independent of the mesh resolution for the discrete second kind BIE.