An interpolation‐based fast multipole method for higher‐order boundary elements on parametric surfaces

Abstract

SUMMARYIn this article, a black-box higher order fast multipole method for solving boundary integral equations on parametric surfaces in three spatial dimensions is proposed. Such piecewise smooth surfaces are the topic of recent studies in isogeometric analysis. Due to the exact surface representation, the rate of convergence of higher order methods is not limited by approximation errors of the surface. An element-wise clustering strategy yields a balanced cluster tree and an efficient numerical integration scheme for the underlying Galerkin method. By performing the interpolation for the fast multipole method directly on the reference domain, the cost complexity in the polynomial degree is reduced by one order. This gain is independent of the application of either H-or H 2 -matrices. In fact, several simplifications in the construction of H 2 -matrices are pointed out, which are a by-product of the surface representation. Extensive numerical examples are provided in order to quantify and qualify the proposed method.

Isogeometric Boundary Elements in Electromagnetism: Rigorous Analysis, Fast Methods, and Examples

“…In this paper, we follow the approach of [20,43], which allows for a fast method with element-wise quadrature and fully avoids redundant evaluations of geometry and fundamental solution. Our method exploits the isogeometric structure and yields a simplified implementation based on interpolation on the unit square.…”

A Numerical Comparison of an Isogeometric and a Parametric Higher Order Raviart–Thomas Approach to the Electric Field Integral Equation

“…To showcase the differences, each of the Figures 5, 6 and 8 shows the accuracy of the solution w.r.t. the number of iterations required for solving the linear system (8).…”

Bembel: The fast isogeometric boundary element C++ library for Laplace, Helmholtz, and electric wave equation

“…The development of the software started in the context of wavelet Galerkin methods on parametric surfaces, see [14], where the integration routines for the Green's function of the Laplacian have been developed and implemented. It was then extended to hierarchical matrices (H-matrices) in [15] and to H 2 -matrices and higher-order B-splines in [16]. With support of B-splines and NURBS for the geometry mappings, the Laplace and Helmholtz code became isogeometric in [7].…”