27 January 2020
SIAM Journal on Numerical Analysis
We consider first-kind weakly singular and hypersingular boundary integral operators for the Laplacian on screens in R3 and their Galerkin discretization by means of low-order piecewise polynomial boundary elements. For the resulting linear systems of equations we propose novel Calder´on-type preconditioners based on (i) new boundary integral operators, which provide the exact inverses of the weakly singular and hypersingular operators on flat disks, and (ii) stable duality pairings relying on dual meshes. On screens obtained as images of the unit disk under bi-Lipschitz transformations, this approach achieves condition numbers uniformly bounded in the meshwidth even on locally refined meshes. Comprehensive numerical tests also confirm its excellent pre-asymptotic performance.
Submitted to ORA: