HILBERT - a MATLAB implementation of adaptive 2D-BEM

被引：0

作者：

Aurada, Markus ^{[1
]}

Ebner, Michael ^{[1
]}

Feischl, Michael ^{[1
]}

Ferraz-Leite, Samuel ^{[3
]}

Fuehrer, Thomas ^{[1
]}

Goldenits, Petra ^{[1
]}

Karkulik, Michael ^{[2
]}

Mayr, Markus ^{[1
]}

Praetorius, Dirk ^{[1
]}

机构：

[1] Vienna Univ Technol, A-1040 Vienna, Austria

[2] Pontificia Univ Catolica Chile, Santiago, Chile

[3] Max Planck Inst Math Sci, D-04103 Leipzig, Germany

来源：

NUMERICAL ALGORITHMS | 2014年 / 67卷 / 01期

基金：

奥地利科学基金会;

关键词：

Boundary element methods; Adaptive mesh-refinement; A posteriori error estimation; MATLAB implementation; HYPERSINGULAR INTEGRAL-EQUATION; OPTIMAL CONVERGENCE RATE; AVERAGING TECHNIQUES; 2-LEVEL METHODS; BEM;

D O I：

10.1007/s11075-013-9771-2

中图分类号：

O29 [应用数学];

学科分类号：

070104 ;

摘要：

We report on the Matlab program package HILBERT. It provides an easily-accessible implementation of lowest order adaptive Galerkin boundary element methods for the numerical solution of the Poisson equation in 2D. The library was designed to serve several purposes: The stable implementation of the integral operators may be used in research code. The framework of Matlab ensures usability in lectures on boundary element methods or scientific computing. Finally, we emphasize the use of adaptivity as general concept and for boundary element methods in particular. In this work, we summarize recent analytical results on adaptivity in the context of BEM and illustrate the use of HILBERT. Various benchmarks are performed to empirically analyze the performance of the proposed adaptive algorithms and to compare adaptive and uniform mesh-refinements. In particular, we do not only focus on mathematical convergence behavior but also on the usage of critical system resources such as memory consumption and computational time. In any case, the superiority of the proposed adaptive approach is empirically supported.

引用

页码：1 / 32

页数：32