Full numerical quadrature of weakly singular double surface integrals in Galerkin boundary element methods

D'ELIA, J.; LAURA BATTAGLIA; ALBERTO CARDONA; STORTI, MARIO ALBERTO

Communications in Numerical Methods in Engineering

John Wiley & Sons

Lugar: Chichester; Año: 2008

1069-8299

When a Galerkin discretization of a boundary integral equation with a weakly singular kernel isperformed over triangles, a double surface integral must be evaluated for each pair of them. If these pairs are not contiguous or not coincident the kernel is regular and a Gauss{Legendre quadrature can be employed. When the pairs have a common edge or a common vertex, then edge and vertex weak singularities appear. If the pairs have both facets coincident the whole integration domain is weakly singular. D. J. Taylor (IEEE Trans. on Antennas and Propagation, 51(7):1630{1637 (2003)) proposed a systematic evaluation based on a reordering and partitioning of the integration domain, together with a use of the Duy transformations in order to remove the singularities, in such a way that a Gauss{Legendre quadrature was performed on three coordinates with an analytic integration in the fourth coordinate. Since this scheme is a bit restrictive because it was designed for electromagnetic kernels, a full numerical quadrature is proposed in order to handle kernels with a weak singularity with a general framework. Numerical tests based on modications of the one proposed by W. Wang and N. Atalla (Comm. in Num. Meth. Eng., 13(11):885{890 (1997)) are included.