An adaptive fast multipole accelerated Poisson solver for complex geometries
arXiv:1610.00823 · doi:10.1016/j.jcp.2017.04.063
Abstract
We present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the sum of a volume integral to account for the source distribution and a layer potential to enforce the desired boundary condition. The volume integral is computed by applying the FMM on a square box that encloses the domain of interest. For the sake of efficiency and convergence acceleration, we first extend the source distribution (the right-hand side in the Poisson equation) to the enclosing box as a function using a fast, boundary integral-based method. We demonstrate on multiply connected domains with irregular boundaries that this continuous extension leads to high accuracy without excessive adaptive refinement near the boundary and, as a result, to an extremely efficient "black box" fast solver.
31 pages
References in corpus (2)
Cited by in corpus (14)
- Adaptive quadrature by expansion for layer potential evaluation in two dimensions
- Partition of Unity Extension of Functions on Complex Domains
- A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations
- Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension
- Integral equation formulation of the biharmonic Dirichlet problem
- A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
- Rapid evaluation of Newtonian potentials on planar domains
- Fast integral equation methods for linear and semilinear heat equations in moving domains
- Higher-order accurate diffuse-domain methods for partial differential equations with Dirichlet boundary conditions in complex, evolving geometries
- High-order Finite Element--Integral Equation Coupling on Embedded Meshes
- An adaptive fast Gauss transform in two dimensions
- Second-kind integral equations for the Laplace-Beltrami problem on surfaces in three dimensions
- Distributed and Adaptive Fast Multipole Method In Three Dimensions
- Boundary Integral Analysis for the Non-homogeneous 3D Stokes Equation