A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
arXiv:2203.05933 · doi:10.1016/j.jcp.2022.111688
Abstract
This article presents a new high-order accurate algorithm for finding a particular solution to a linear, constant-coefficient partial differential equation (PDE) by means of a convolution of the volumetric source function with the Green's function in complex geometries. Utilizing volumetric domain decomposition, the integral is computed over a union of regular boxes (lending the scheme compatibility with adaptive box codes) and triangular regions (which may be potentially curved near boundaries). Singular and near-singular quadrature is handled by converting integrals on volumetric regions to line integrals bounding a reference volume cell using cell mappings and elements of the Poincaré lemma, followed by leveraging existing one-dimensional near-singular and singular quadratures appropriate to the singular nature of the kernel. The scheme achieves compatibility with fast multipole methods (FMMs) and thereby optimal asymptotic complexity by coupling global rules for target-independent quadrature of smooth functions to local target-dependent singular quadrature corrections, and it relies on orthogonal polynomial systems on each cell for well-conditioned, high-order and efficient (with respect to number of required volume function evaluations) approximation of arbitrary volumetric sources. Our domain discretization scheme is naturally compatible with standard meshing software such as Gmsh, which are employed to discretize a narrow region surrounding the domain boundaries. We present 8th-order accurate results, demonstrate the success of the method with examples showing up to 12-digit accuracy on complex geometries, and, for static geometries, our numerical examples show well over of evaluation time of the particular solution is spent in the FMM step.
arXiv comment33 pages, 9 figures
References in corpus (10)
- Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
- A state redistribution algorithm for finite volume schemes on cut cell meshes
- An adaptive fast multipole accelerated Poisson solver for complex geometries
- Adaptive quadrature by expansion for layer potential evaluation in two dimensions
- Partition of Unity Extension of Functions on Complex Domains
- Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme
- A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations
- Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D
- Optimization of Fast Algorithms for Global Quadrature by Expansion Using Target-Specific Expansions
- General-purpose kernel regularization of boundary integral equations via density interpolation