Adaptive boundary element methods for the computation of the electrostatic capacity on complex polyhedra
arXiv:1901.08393 · doi:10.1016/j.jcp.2019.07.036
Abstract
The accurate computation of the electrostatic capacity of three dimensional objects is a fascinating benchmark problem with a long and rich history. In particular, the capacity of the unit cube has widely been studied, and recent advances allow to compute its capacity to more than ten digits of accuracy. However, the accurate computation of the capacity for general three dimensional polyhedra is still an open problem. In this paper, we propose a new algorithm based on a combination of ZZ-type a posteriori error estimation and effective operator preconditioned boundary integral formulations to easily compute the capacity of complex three dimensional polyhedra to 5 digits and more. While this paper focuses on the capacity as a benchmark problem, it also discusses implementational issues of adaptive boundary element solvers, and we provide codes based on the boundary element package Bempp to make the underlying techniques accessible to a wide range of practical problems.
References in corpus (7)
- On 2D Newest Vertex Bisection: Optimality of Mesh-Closure and H1-Stability of L2-Projection
- An abstract analysis of optimal goal-oriented adaptivity
- Adaptive Boundary Element Methods
- Local inverse estimates for non-local boundary integral operators
- Adaptive BEM with optimal convergence rates for the Helmholtz equation
- Adaptive BEM with inexact PCG solver yields almost optimal computational costs
- ZZ-type aposteriori error estimators for adaptive boundary element methods on a curve