edoc

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

Dölz, Jürgen and Harbrecht, Helmut and Peters, Michael. (2016) An interpolation-based fast multipole method for higher order boundary elements on parametric surfaces. International Journal for Numerical Methods in Engineering, 108 (13). pp. 1705-1728.

[img] PDF - Accepted Version
3247Kb

Official URL: http://edoc.unibas.ch/44933/

Downloads: Statistics Overview

Abstract

In 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 (mathcal{H})- or (mathcal{H}^2)-matrices. In fact, several simplifications in the construction of (mathcal{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. In 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 reducedby one order. This gain is independent of the application of either H - or H 2- matrices. In fact, several simplificationsin 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.
Faculties and Departments:05 Faculty of Science > Departement Mathematik und Informatik > Mathematik > Computational Mathematics (Harbrecht)
UniBasel Contributors:Harbrecht, Helmut and Dölz, Jürgen and Peters, Michael
Item Type:Article, refereed
Article Subtype:Research Article
Publisher:Wiley
ISSN:0029-5981
e-ISSN:1097-0207
Note:Publication type according to Uni Basel Research Database: Journal article -- The final publication is available at Wiley via DOI link.
Language:English
Identification Number:
edoc DOI:
Last Modified:29 Nov 2016 15:43
Deposited On:29 Nov 2016 15:24

Repository Staff Only: item control page