A mapping method for numerical evaluation of two-dimensional integrals with 1/r singularity

Singular integrals occur commonly in applications of the boundary element method (BEM). A simple mapping method is presented here for the numerical evaluation of two-dimensional integrals in which the integrands, at worst, are O(1/r) (r being the distance from a source to a field point). This mapping transforms such integrals over general curved triangles into regular 2-D integrals. Over flat and curved quadratic triangles, regular line integrals are obtained, and these can be easily evaluated by standard Gaussian quadrature. Numerical tests on some typical singular integrals, encountered in BEM applications, demonstrate the accuracy and efficacy of the method.