Partition of unity extension of functions on complex domains

DOI10.1016/J.JCP.2018.08.012zbMATH Open1416.65475arXiv1712.08461OpenAlexW2779258471WikidataQ129422657 ScholiaQ129422657MaRDI QIDQ2002228FDOQ2002228

Authors: Fredrik Fryklund, Erik Lehto, Anna-Karin Tornberg 

Publication date: 11 July 2019

Published in: Journal of Computational Physics (Search for Journal in Brave)

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

Full work available at URL: https://arxiv.org/abs/1712.08461

Recommendations

• Partition of unity finite element method implementation for Poisson equation
• THE PARTITION OF UNITY METHOD
• Fourier embedded domain methods: Extending a function defined on an irregular region to a rectangle so that the extension is spatially periodic and \(C^{\infty}\)
• The piecewise polynomial partition of unity functions for the generalized finite element methods
• Fourier embedded domain methods: Periodic and \(C^{\infty}\) extension of a function defined on an irregular region to a rectangle via convolution with Gaussian kernels

zbMATH Keywords

radial basis functionboundary integral methodlinear elliptic partial differential equationpartition of unityfunction extensionembedded domain

Mathematics Subject Classification ID

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Fundamental solutions, Green's function methods, etc. for boundary value problems involving PDEs (65N80) Boundary value problems for second-order elliptic equations (35J25) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35)

Cites Work

• Stable Computations with Gaussian Radial Basis Functions
• Meshfree approximation methods with Matlab. With CD-ROM.
• Spectral Methods in MATLAB
• Title not available (Why is that?)
• Layer potentials and regularity for the Dirichlet problem for Laplace's equation in Lipschitz domains
• Error estimates for interpolation by compactly supported radial basis functions of minimal degree
• The Numerical Solution of Integral Equations of the Second Kind
• Fast integral equation methods for the modified Helmholtz equation
• Fast integral equation methods for Rothe's method applied to the isotropic heat equation
• Impossibility of fast stable approximation of analytic functions from equispaced samples
• Accelerating the Nonuniform Fast Fourier Transform
• Sampling inequalities for infinitely smooth functions, with applications to interpolation and machine learning
• Laplace's equation and the Dirichlet-Neumann map in multiply connected domains
• On the evaluation of layer potentials close to their sources
• Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions
• Title not available (Why is that?)
• Title not available (Why is that?)
• Spectral domain embedding for elliptic PDEs in complex domains
• High-order unconditionally stable FC-AD solvers for general smooth domains. I: Basic elements
• A spectral FC solver for the compressible Navier-Stokes equations in general domains. I: Explicit time-stepping
• A high-order 3D boundary integral equation solver for elliptic PDEs in smooth domains
• High-order unconditionally stable FC-AD solvers for general smooth domains. II: Elliptic, parabolic and hyperbolic PDEs; theoretical considerations
• An Integral Equation Approach to the Incompressible Navier--Stokes Equations in Two Dimensions
• Fast convolution with free-space Green's functions
• Fast Ewald summation for free-space Stokes potentials
• Immersed boundary smooth extension (IBSE): a high-order method for solving incompressible flows in arbitrary smooth domains
• A sharp-interface active penalty method for the incompressible Navier-Stokes equations
• An adaptive fast multipole accelerated Poisson solver for complex geometries
• Error estimation for quadrature by expansion in layer potential evaluation
• A least squares radial basis function partition of unity method for solving PDEs

Cited In (19)

• A fast integral equation method for the two-dimensional Navier-Stokes equations
• A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
• Quadrature by fundamental solutions: kernel-independent layer potential evaluation for large collections of simple objects
• Efficient least squares approximation and collocation methods using radial basis functions
• Extension and division on complex manifolds
• An integral equation method for the advection-diffusion equation on time-dependent domains in the plane
• Direct RBF-PU method combined with the tangent plane approach for parabolic equation on surface
• The smooth forcing extension method: a high-order technique for solving elliptic equations on complex domains
• Rapid Evaluation of Newtonian Potentials on Planar Domains
• Numerical analytic continuation
• A correction function-based kernel-free boundary integral method for elliptic PDEs with implicitly defined interfaces
• Immersed hybrid difference methods for elliptic boundary value problems by artificial interface conditions
• Stable and accurate least squares radial basis function approximations on bounded domains
• An FMM Accelerated Poisson Solver for Complicated Geometries in the Plane Using Function Extension
• Two-dimensional Fourier continuation and applications
• AAA interpolation of equispaced data
• Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension
• A locally corrected multiblob method with hydrodynamically matched grids for the Stokes mobility problem
• An integral equation-based numerical method for the forced heat equation on complex domains

Uses Software

• Matlab
• rbf_qr

This page was built for publication: Partition of unity extension of functions on complex domains

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2002228)