High order corrected trapezoidal rules for a class of singular integrals

Federico Izzo,Olof Runborg,Richard Tsai
DOI: https://doi.org/10.1007/s10444-023-10060-0
2023-07-26
Abstract:We present a family of high order trapezoidal rule-based quadratures for a class of singular integrals, where the integrand has a point singularity. The singular part of the integrand is expanded in a Taylor series involving terms of increasing smoothness. The quadratures are based on the trapezoidal rule, with the quadrature weights for Cartesian nodes close to the singularity judiciously corrected based on the expansion. High order accuracy can be achieved by utilizing a sufficient number of correction nodes around the singularity to approximate the terms in the series expansion. The derived quadratures are applied to the Implicit Boundary Integral formulation of surface integrals involving the Laplace layer kernels.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the precision in the numerical calculation of singular integrals of a class of integrands with point singularities. Specifically, when using the traditional Trapezoidal Rule to deal with these singular integrals, due to the influence of singularities, the calculation precision will decrease significantly. To solve this problem, the author proposes a high - order correction method based on the Trapezoidal Rule. By adding appropriate correction nodes near the singularity, the calculation precision is improved. ### Specific Problem Description 1. **Limitations of the Traditional Trapezoidal Rule**: - The Trapezoidal Rule is a simple and robust algorithm for approximate integration and has second - order precision in general cases. - But for integrands with singularities (such as singular integrals in boundary integral equations), the precision of the Trapezoidal Rule will drop sharply, making it unsuitable for this type of integrals. 2. **Challenges in Singularity Handling**: - The integrand \( f \) has a singularity at a point \( x_0 \), that is, \( f(x)=s(x)v(x) \), where \( s(x) \) is the singular part and \( v(x) \) is the smooth part. - Direct application of the Trapezoidal Rule will lead to low - precision results, so a systematic method is needed to correct this error. 3. **Objective**: - Develop a systematic method to use the high - order corrected Trapezoidal Rule to deal with a class of integrals with singularities, so as to achieve high - precision numerical calculations. ### Solution The method proposed by the author includes the following key steps: 1. **Taylor Expansion**: - Perform a Taylor expansion on the singular part \( s(x) \) near the singularity point to obtain a series of gradually smoothed terms. - For example, in the two - dimensional case, \( s(x)=\frac{1}{|x|}\ell(|x|,\frac{x}{|x|}) \), where \( \ell \) is a smooth function. 2. **Correction Weights**: - For each expansion term, introduce correction weights to compensate for the error of the Trapezoidal Rule near the singularity. - These weights are calculated by the difference between the integral and the discrete sum of smooth test functions. 3. **Gradually Improving Precision**: - By gradually increasing the number of correction nodes, the precision of the integral can be improved. For example, start from the first - order correction and gradually expand to higher - order correction methods. - Eventually, by combining multiple correction terms, a high - order precision Trapezoidal Rule correction formula can be constructed. ### Application Scenarios This method is especially suitable for Implicit Boundary Integral Methods (IBIM), which are very useful in dealing with boundary integral equations of elliptic problems. Through this method, surface integrals, integral equations and variational problems can be accurately calculated without the need for explicit surface parameterization. In conclusion, this paper aims to improve the numerical calculation precision of integrals with singularities by improving the Trapezoidal Rule, thereby providing a more accurate tool for scientific computing involving singular integrals.