High-Order Meshfree Surface Integration, Including Singular Integrands

Abstract.

We develop and test high-order methods for integration on surface point clouds. The task of integrating a function on a surface arises in a range of applications in engineering and the sciences, particularly those involving various integral methods for partial differential equations. Mesh-based methods require a curved mesh for high-order convergence, which can be difficult to reliably obtain on many surfaces, and most meshfree methods require the ability to integrate a set of functions (such as radial basis functions) exactly on the domain of interest; these integrals are generally not known in closed form on most surfaces. We describe two methods for integrating on arbitrary, piecewise-smooth surfaces with or without boundary. Our approaches do not require a particular arrangement of points or an initial triangulation of the surface, making them completely meshfree. We also show how the methods can be extended to handle singular integrals while maintaining high accuracy without changing the point density near singularities.

Reproducibility of computational results. This paper has been awarded the “SIAM Reproducibility Badge: Code and Data Available” as a recognition that the authors have followed reproducibility principles valued by SISC and the scientific computing community. Code and data that allow readers to reproduce the results in this paper are available at https://github.com/venndaniel/meshfree_integration and in the supplementary materials ( M181158_Supplementary_Materials.zip [19.4KB]), linked from the main article webpage. [Formula: see text]