Isoparametric singularity extraction technique for 3D potential problems in BEM

Abstract

To solve boundary integral equations for potential problems using collocation Boundary Element Method (BEM) on smooth curved 3D geometries, an analytical singularity extraction technique is employed. By adopting the isoparametric approach, curved geometries that are represented by mapped rectangles or triangles from the parametric domain are considered. The singularity extraction on the governing singular integrals can be performed either as an operation of subtraction or division, each having some advantage.A particular series expansion of a singular kernel about a source point is investigated. The series in the parametric coordinates consists of functions of a type Rpxqyr, where R is a square root of a quadratic bivariate homogeneous polynomial with coefficients corresponding to the first fundamental form of a smooth surface, and p,q,r are integers, satisfying p≤−1 and q,r≥0. By extracting more terms from the series expansion of the singular kernel, the smoothness of the regularized kernel at the source point can be increased. Analytical formulae for integrals of such terms are obtained from antiderivatives of Rpxqyr, using recurrence formulae, and by evaluating them at the edges of rectangular or triangular parametric domains.Numerical tests demonstrate that the singularity extraction technique can be a useful prerequisite for a numerical quadrature scheme to obtain high order of convergence for the governing singular integrals in 3D collocation BEM with respect to the number of quadrature nodes or with respect to size of the discretization mesh cells.