# Accurate computation of Galerkin double surface integrals in the 3-D boundary element method

 dc.contributor.author Adelman, Ross dc.contributor.author Gumerov, Nail A. dc.contributor.author Duraiswami, Ramani dc.date.accessioned 2015-06-23T07:44:43Z dc.date.available 2015-06-23T07:44:43Z dc.date.issued 2015-05-29 dc.identifier https://doi.org/10.13016/M23S5D dc.identifier.uri http://hdl.handle.net/1903/16394 dc.description.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 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. en_US 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. dc.language.iso en_US en_US dc.relation.ispartofseries UM Computer Science Department;CS-TR-5043 dc.relation.ispartofseries UMIACS;UMIACS-TR-2015-02 dc.title Accurate computation of Galerkin double surface integrals in the 3-D boundary element method en_US dc.type Other en_US
﻿