Fredrik Fryklund

Fredrik Fryklund, Erik Lehto, Anna-Karin Tornberg

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}$.

Fredrik Fryklund

Integral equation methods provide an effective framework for solving partial differential equations, but their applicability typically relies on the availability of explicit free-space Green's functions. For coupled systems arising in multiphysics applications, such Green's functions are generally not available, limiting the scope of classical potential theory-based approaches. In this work, we introduce a formulation of potential theory that avoids explicit use of Green's functions entirely, relying instead on the Fourier symbol of the governing operator. The central idea is a parabolic regularization of the symbol, which yields a decomposition of the solution into a smooth, nonlocal component and a spatially localized residual. For the localized component, we derive explicit asymptotic expansions for volume, single layer, and double layer potentials in powers of a length scale parameter $\varepsilon$. The coefficients are expressed in terms of local geometric quantities and derivatives of the source data. The derivation is carried out entirely in the Fourier domain and applies to the Poisson equation in two and three dimensions, as well as to a class of coupled strongly elliptic systems.