Numerical quadrature over smooth surfaces with boundaries

Abstract This paper describes a high order accurate method to calculate integrals over curved surfaces with boundaries. Given data locations that are arbitrarily distributed over the surface, together with some functional description of the surface and its boundary, the algorithm produces matching quadrature weights. This extends on the authors' earlier methods for integrating over the surface of a sphere and over arbitrarily shaped smooth closed surfaces by also considering domain boundaries. The core approach consists again of combining RBF-FD (radial basis function-generated finite difference) approximations for curved surface triangles, which together make up the full surface. The provided examples include both curved and flat domains. In the highly special case of equi-spaced nodes over a regular interval in 1-D, the method provides a new opportunity for improving on the classical Gregory enhancements of the trapezoidal rule.

[1]  Johan M. De Villiers,et al.  Gregory type quadrature based on quadratic nodal spline interpolation , 2000, Numerische Mathematik.

[2]  Bengt Fornberg,et al.  On spherical harmonics based numerical quadrature over the surface of a sphere , 2014, Advances in Computational Mathematics.

[3]  Bengt Fornberg,et al.  On spherical harmonics based numerical quadrature over the surface of a sphere , 2014, Advances in Computational Mathematics.

[4]  Lloyd N. Trefethen,et al.  Euler–Maclaurin and Gregory interpolants , 2016, Numerische Mathematik.

[5]  Lloyd N. Trefethen,et al.  The Exponentially Convergent Trapezoidal Rule , 2014, SIAM Rev..

[6]  Louis J. Wicker,et al.  Enhancing finite differences with radial basis functions: Experiments on the Navier-Stokes equations , 2016, J. Comput. Phys..

[7]  Joseph D. Ward,et al.  Kernel based quadrature on spheres and other homogeneous spaces , 2012, Numerische Mathematik.

[8]  Bengt Fornberg,et al.  A primer on radial basis functions with applications to the geosciences , 2015, CBMS-NSF regional conference series in applied mathematics.

[9]  J. Halton On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals , 1960 .

[10]  Bengt Fornberg,et al.  Solving PDEs with radial basis functions * , 2015, Acta Numerica.

[11]  Bengt Fornberg,et al.  Fast generation of 2-D node distributions for mesh-free PDE discretizations , 2015, Comput. Math. Appl..

[12]  J A Reeger,et al.  Numerical quadrature over smooth, closed surfaces , 2016, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[13]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs , 2017, J. Comput. Phys..