Computation of Galerkin Double Surface Integrals in the 3-D Boundary Element Method

Abstract

Many boundary element integral equation kernels are based on the Green's functions of the Laplace and Helmholtz equations in three dimensions. These include, for example, the Laplace, Helmholtz, elasticity, Stokes, and Maxwell's equations. Integral equation formulations lead to more compact, but dense linear systems. These dense systems are often solved iteratively via Krylov subspace methods, which may be accelerated via the fast multipole method. There are advantages to Galerkin formulations for such integral equations, as they treat problems associated with kernel singularity, and lead to symmetric and better conditioned matrices. However, the Galerkin method requires each entry in the system matrix to be created via the computation of a double surface integral over one or more pairs of triangles. There are a number of semi-analytical methods to treat these integrals, which all have some issues, and are discussed in this paper. We present novel methods to compute all the integrals that arise in Galerkin formulations involving kernels based on the Laplace and Helmholtz Green's functions to any specified accuracy. Integrals involving completely geometrically separated triangles are non-singular and are computed using a technique based on spherical harmonics and multipole expansions and translations, which results in the integration of polynomial functions over the triangles. Integrals involving cases where the triangles have common vertices, edges, or are coincident are treated via scaling and symmetry arguments, combined with automatic recursive geometric decomposition of the integrals. Example results are presented, and the developed software is available as open source.

Efficient Fast Multipole Accelerated Boundary Elements via Recursive Computation of Multipole Expansions of Integrals

“…It is remarkable that the present method can be also used in computations of the Galerkin integrals (e.g., see [2]). Indeed, the no-touch cases can be computed using the multipole expansions.…”

“…When used with the FMM accelerated iterative methods (FMMBEM), the matrix vector products are evaluated approximately, and these integrals are only needed as far as their elements contribute to the final product. While there exist methods to evaluate such integrals exactly or approximately with controlled accuracy (e.g., [2]), this approximation may not be consistent with that needed in the FMM. We need exact integrals for some elements (those in the near field) and for the others (in the far field) the integrals are only needed as far as their contribution to the consolidated far-field multipole expansions are concerned.…”

Improving Efficiency of the Secondary Sources Method for Modeling of the Three-Dimensional Electromagnetic Field of Eddy Currents

“…The disadvantages of the collocation method are known and are noted by many authors who are investigating the method of boundary integral equations. In some works, it is proposed to solve the problem of collocation points using the Galerkin method to obtain a SLAE from the corresponding boundary integral equation [17][18][19]. Mathematically, this method leads to SLAE, the coefficients of which are expressed in terms of double integrals.…”

GPU accelerated fast multipole boundary element method for simulation of 3D bubble dynamics in potential flow

“…There exist extensive literature for accurate numerical and analytical evaluation of the integrals of the Green's function and its derivatives over triangles (e.g., [49,50]). However, efficient use of the FMM for large N requires numerically inexpensive quadratures and approximations, which brings forward strategies, such as described [51].…”