Approximation of Integrals for Boundary Element Methods

A new method for approximating two-dimensional integrals ∫ B f(x)μ (dx) over surfaces B ⊂ R³ is introduced where μ is the standard measure of surface area. Such integrals typically occur in boundary element methods. The algorithm is based on triangulations T : = ∪ Tᵢ approximating B. Under the assumption that the surface B is given implicitly by an equation $H(x) = 0$, a retraction P:U ⊃ B → B is used to obtain a curved subdivision B = ∪ Bᵢ via Bᵢ : = PTᵢ. Except in very special cases, this retraction is not analytically accessible, but is generated by a subroutine. Hence standard multiple integral techniques are not available. Thus, the approach given here differs from the usual panel method. It is shown how to calculate the integrals as precisely as wished. Two numerical examples are given. The first integrand f(x) ≡ 1 is regular, and it is shown that a very accurate extrapolation method can be used. The second integrand f(x) ~ ||x - x₀ ||- 1 is singular, and an adaptive refinement procedure is displayed.