Partition of Unity Extension of Functions on Complex Domains

Abstract

We introduce an efficient algorithm, called partition of unity extension or PUX, to construct an extension of desired regularity of a function given on a complex multiply connected domain in $2D$. Function extension plays a fundamental role in extending the applicability of boundary integral methods to inhomogeneous partial differential equations with embedded domain techniques. Overlapping partitions are placed along the boundaries, and a local extension of the function is computed on each patch using smooth radial basis functions; a trivially parallel process. A partition of unity method blends the local extrapolations into a global one, where weight functions impose compact support. The regularity of the extended function can be controlled by the construction of the partition of unity function. We evaluate the performance of the PUX method in the context of solving the Poisson equation on multiply connected domains using a boundary integral method and a spectral solver. With a suitable choice of parameters the error converges as a tenth order method down to $10^{-14}$.

References in corpus (5)

• Fast convolution with free-space Green's functions
• Immersed Boundary Smooth Extension (IBSE): A high-order method for solving incompressible flows in arbitrary smooth domains
• An adaptive fast multipole accelerated Poisson solver for complex geometries
• Fast Ewald summation for free-space Stokes potentials
• Estimation of quadrature errors in layer potential evaluation using quadrature by expansion

Cited by in corpus (9)

• A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations
• A Locally Corrected Multiblob Method with Hydrodynamically Matched Grids for the Stokes Mobility Problem
• Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension
• A correction function-based kernel-free boundary integral method for elliptic PDEs with implicitly defined interfaces
• A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
• The Smooth Forcing Extension Method: A High-Order Technique for Solving Elliptic Equations on Complex Domains
• Rapid evaluation of Newtonian potentials on planar domains
• Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation
• An integral equation method for the advection-diffusion equation on time-dependent domains in the plane