Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

TL;DR

This paper presents the scaled boundary cubature scheme for efficient and accurate numerical integration over complex planar regions, including those with curved boundaries and singularities, demonstrating broad applicability and superior performance.

Contribution

The paper introduces the SBC scheme that generalizes numerical integration over polygons and curved regions, including methods for handling singularities and non-homogeneous functions.

Findings

The SBC scheme effectively integrates over convex and nonconvex regions.

It achieves high accuracy and efficiency compared to existing methods.

The scheme handles singularities and non-homogeneous functions well.

Abstract

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by $m$ curves to integration over $m$ regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. The SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.