A robust solver for elliptic PDEs in 3D complex geometries
arXiv:2002.04143 · doi:10.1016/j.jcp.2021.110511
Abstract
We develop a boundary integral equation solver for elliptic partial differential equations on complex \threed geometries. Our method is efficient, high-order accurate and robustly handles complex geometries. A key component is our singular and near-singular layer potential evaluation scheme, \qbkix: a simple extrapolation of the solution along a line to the boundary. We present a series of geometry-processing algorithms required for \qbkix to run efficiently with accuracy guarantees on arbitrary geometries and an adaptive upsampling scheme based on a iteration-free heuristic for quadrature error. We validate the accuracy and performance with a series of numerical tests and compare our approach to a competing local evaluation method.
References in corpus (24)
- Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
- Stability and convergence of the Method of Fundamental Solutions for Helmholtz problems on analytic domains
- Theory and implementation of -matrix based iterative and direct solvers for Helmholtz and elastodynamic oscillatory kernels
- Fast algorithms for Quadrature by Expansion I: Globally valid expansions
- High-volume fraction simulations of two-dimensional vesicle suspensions
- Efficient numerical solution of acoustic scattering from doubly-periodic arrays of axisymmetric objects
- A Fast Algorithm for Quadrature by Expansion in Three Dimensions
- A high-order integral solver for scalar problems of diffraction by screens and apertures in three dimensional space
- A Simple Method for Computing Singular or Nearly Singular Integrals on Closed Surfaces
- Newton's method in practice: finding all roots of polynomials of degree one million efficiently
- A local target specific quadrature by expansion method for evaluation of layer potentials in 3D
- Adaptive quadrature by expansion for layer potential evaluation in two dimensions
- Estimation of quadrature errors in layer potential evaluation using quadrature by expansion
- Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme
- Taylor States in Stellarators: A Fast High-order Boundary Integral Solver
- A Fast Algorithm with Error Bounds for Quadrature by Expansion
- Regularized Single and Double Layer Integrals in 3D Stokes Flow
- Asymptotic analysis for close evaluation of layer potentials
- Adaptive BEM with optimal convergence rates for the Helmholtz equation
- Scalable Simulation of Realistic Volume Fraction Red Blood Cell Flows through Vascular Networks
- Optimization of Fast Algorithms for Global Quadrature by Expansion Using Target-Specific Expansions
- Adaptive boundary element methods for the computation of the electrostatic capacity on complex polyhedra
- On the Approximation of Local Expansions of Laplace Potentials by the Fast Multipole Method
- Asymptotic approximations for the close evaluation of double-layer potentials