Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 918083, 9 pages
http://dx.doi.org/10.1155/2015/918083

Research Article

The Trapezoidal Rule for Computing Cauchy Principal Value Integral on Circle

Jin Li1,2

¹School of Science, Shandong Jianzhu University, Jinan 250101, China
²School of Mathematics, Shandong University, Jinan 250100, China

Received 27 July 2015; Accepted 20 September 2015

Academic Editor: Kishin Sadarangani

Copyright © 2015 Jin Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The composite trapezoidal rule for the computation of Cauchy principal value integral with the singular kernel is discussed. Our study is based on the investigation of the pointwise superconvergence phenomenon; that is, when the singular point coincides with some a priori known point, the convergence rate of the trapezoidal rule is higher than what is globally possible. We show that the superconvergence rate of the composite trapezoidal rule occurs at middle of each subinterval and obtain the corresponding superconvergence error estimate. Some numerical examples are provided to validate the theoretical analysis.

1. Introduction

Consider the Cauchy principal integral

where denotes a Cauchy principal value integral and is the singular point.

There are several different definitions which can be proved equally, such as the definition of subtraction of the singularity, regularity definition, and direct definition.

In this paper we adopt the following one:

Cauchy principal value integral has recently attracted a lot of attention; see, for example, [1]. The main reason for this interest is probably due to the fact that integral equations with Cauchy principal value integrals have been shown to be an adequate tool [1] for the modeling of many physical situations, such as acoustics, fluid mechanics, elasticity, fracture mechanics, and electromagnetic scattering problems. Numerous work has been devoted in developing efficient quadrature formulas, such as the Gaussian method [2–8], the Newton-Cotes methods [9–13], spline methods [14, 15], and some other methods [16–22].

It is the aim of this paper to investigate the superconvergence phenomenon of trapezoidal rule and, in particular, to derive error estimates. In this paper the density function is replaced by the approximation function while the singular kernel is computed in each subinterval, where is the trapezoidal rule. This method is different from the semidiscrete methods and the order of singularity kernel can be reduced somehow which was firstly presented by Linz in the paper to calculate the hypersingular integral on interval. He used the trapezoidal rule and Simpson rule to approximate the density function and the convergence rate was , , when the singular point is always located at the middle of certain subinterval.

The superconvergence of composite Newton-Cotes rules for Hadamard finite-part integrals was studied in [23, 24], where the superconvergence rate and the superconvergence point were presented, respectively. In [25, 26] the classical composite midpoint rectangle rule and classical composite trapezoidal rectangle rule for the computation of Cauchy principal value integrals are discussed. When the singular point coincides with some a priori known point, the convergence rate of the midpoint rectangle rule is higher than the global one, the same as the Riemann integral, which is called superconvergence phenomenon.

This paper focuses on the superconvergence of trapezoidal rule for Cauchy principle integrals in a circle. Based on the investigation of the pointwise superconvergence phenomenon, that is, when the singular point coincides with some a priori known point, the convergence rate of the trapezoidal rule is higher than what is globally possible. We show that the superconvergence rate of the composite trapezoidal rule occurs at middle of each subinterval and obtain the corresponding superconvergence error estimate. We prove both theoretically and numerically that the composite trapezoidal rule reaches the superconvergence rate when the local coordinate of the singular point is 0.

The rest of this paper is organized as follows. In Section 2, after introducing some basic formulas of the rectangle rule, we present the main results. In Section 3 we perform the proof. Finally, several numerical examples are provided to validate our analysis.

2. Main Result

Let be a uniform partition of the interval with mesh size . Define by the piecewise line interpolant for :

and a linear transformation

from the reference element to the subinterval . Replacing in (2) with gives the composite trapezoidal rule:

where denote the Cotes coefficients given by

We also define

Theorem 1. Assume . For the trapezoidal rule is defined as (5). Assuming that , there exists a positive constant , independent of and , such that

where

Proof. Let ; then we have . Thus

For the first part of (10), we have

For the second part of (10), we have

For the third part of (10), we have

Combining (11), (13), and (14) together, the proof is completed.

Set

Lemma 2. Assume with . Let be defined by (15); then there holds that

Proof. For , by the definition of Cauchy principal value integral, we have

For , taking integration by parts on the correspondent Riemann integral, we have

Now, by using the well-known identity,

The proof is completed.

Lemma 3. Under the same assumptions of Lemma 2, there holds that

Proof. By (15), we have

The proof is completed.

Theorem 4. Assume . For the trapezoidal rule is defined as (5). Assuming that , there exists a positive constant , independent of and , such that

where

is defined as (9).

It is known that the global convergence rate of the composite trapezoidal rule is lower than Riemann integral.

3. Proof of the Theorem

In this section, we study the superconvergence of the composite trapezoidal rule for Cauchy principle integrals.

3.1. Preliminaries

In the following analysis, will denote a generic constant that is independent of and and it may have different values in different places.

Lemma 5. Under the same assumptions of Theorem 4, it holds that

where

Proof. Performing Taylor expansion of at the point , we have

Similarly, we have

We set

Lemma 6. Let , denote to be the error functional for the composite trapezoidal rule, and assume for any ; then there holds

Proof. As , we get . We have

For the first part of (30), by the definition of

then we have

then we get

For the second part of (30), we have

And the proof is completed.

Proof of Theorem 4. By Lemma 5, we have

By the definition of , we have

Putting (35), (36) together yields

where

Now we estimate term by term. For the first part of , we have

For the second part of , there is no singularity and we have

The third and fourth parts of can be similarly obtained. For , by Lemmas 5 and 6, then we get

and the proof is completed.

Remark 7. By the identity in [27],

then we get

We can prove that the error expansion of this paper has the same error expansion of [25] for the Cauchy principal integral defined on the interval. Based on the error expansion, we get the same superconvergence point with local coordinate point equal to zero.
Based on Theorem 4, we present the modified trapezoidal rule

4. Numerical Example

In this section, computational results are reported to confirm our theoretical analysis.

Example. We consider the singular integral with . and with are the superconvergence point.

From Tables 1 and 2, we know that the superconvergence point is 0 with the coordinate location of singular point equal to zero, while, for the local coordinate of singular point not equal to zero, it is not convergence in general which coincides with our analysis.

Table 1: An error estimate of the trapezoidal rule with .

Table 2: An error estimate of the trapezoidal rule with .

For the modified classical trapezoidal rule, from Tables 3 and 4, for the nonsuperconvergence point and the superconvergence point, we all get the superconvergence phenomenon.

Table 3: An error estimate of the modified trapezoidal rule with .

Table 4: An error estimate of the modified trapezoidal rule with .

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work of Jin Li was supported by the National Natural Science Foundation of China (Grant no. 11471195, Grant no. 11101247, and Grant no. 11201209) and China Postdoctoral Science Foundation (Grant nos. 2013M540541 and 2015T80703).

References

D. H. Yu, Natural Boundary Integral Method and Its Applications, Kluwer Academic Publishers, 2002.

G. Criscuolo and G. Mastroianni, “On the uniform convergence of Gaussian quadrature rules for Cauchy principal value integrals,” Numerische Mathematik, vol. 54, no. 4, pp. 445–461, 1989.

G. Criscuolo and G. Mastroianni, “Gaussian formulas for computing Cauchy principal value integrals, and their convergence,” Calcolo, vol. 22, no. 3, pp. 391–411, 1985.

K. Diethelm, “Gaussian quadrature formulae of the third kind for Cauchy principal value integrals: basic properties and error estimates,” Journal of Computational and Applied Mathematics, vol. 65, no. 1–3, pp. 97–114, 1995.

K. Diethelm, “Asymptotically sharp error bounds for a quadrature rule for Cauchy principal value integrals based on piecewise linear interpolation,” Approximation Theory and Its Applications, vol. 11, no. 4, pp. 78–89, 1995.

N. I. Ioakimidis, “On the uniform convergence of Gaussian quadrature rules for Cauchy principal value integrals and their derivatives,” Mathematics of Computation, vol. 44, no. 169, pp. 191–198, 1985.

T. Hasegawa, “Uniform approximations to finite Hilbert transform and its derivative,” Journal of Computational and Applied Mathematics, vol. 163, no. 1, pp. 127–138, 2004.

G. Monegato, “Convergence of product formulas for the numerical evaluation of certain two-dimensional Cauchy principal value integrals,” Numerische Mathematik, vol. 43, no. 2, pp. 161–173, 1984.

S. Amari, “Evaluation of Cauchy principal-value integrals using modified Simpson rules,” Applied Mathematics Letters, vol. 7, no. 3, pp. 19–23, 1994.

K. Diethelm, “Modified compound quadrature rules for strongly singular integrals,” Computing, vol. 52, no. 4, pp. 337–354, 1994.

P. Köhler, “On the error of quadrature formulae for Cauchy principal value integrals based on piecewise interpolation,” Approximation Theory and Its Applications, vol. 13, no. 3, pp. 58–69, 1997.

J. Li and D. H. Yu, “The superconvergence of certain two-dimensional Cauchy principal value integrals,” CMES, vol. 71, no. 4, pp. 331–346, 2011.

J. Li and D. H. Yu, “The superconvergence of certain two-dimensional Hilbert singular integrals,” Computer Modeling in Engineering & Sciences, vol. 82, no. 3-4, pp. 233–252, 2011.

C. Dagnino and E. Santi, “On the evaluation of one-dimensional Cauchy principal value integrals by rules based on cubic spline interpolation,” Computing, vol. 43, no. 3, pp. 267–276, 1990.

A. P. Orsi, “Spline approximation for Cauchy principal value integrals,” Journal of Computational and Applied Mathematics, vol. 30, no. 2, pp. 191–201, 1990.

G. Behforooz, “Approximation of Cauchy principal value integrals by piecewise Hermite quartic polynomials by spline,” Applied Mathematics Letters, vol. 5, no. 1, pp. 75–78, 1992.

P. Kim and U. J. Choi, “A quadrature rule of interpolatory type for Cauchy integrals,” Journal of Computational and Applied Mathematics, vol. 126, no. 1-2, pp. 207–220, 2000.

R. Kress, I. H. Sloan, and F. Stenger, “A sinc quadrature method for the double-layer integral equation in planar domains with corners,” Journal of Integral Equations & Applications, vol. 10, no. 3, pp. 291–317, 1998.

P. Junghanns and B. Silbermann, “Numerical analysis for one-dimensional Cauchy singular integral equations,” Journal of Computational and Applied Mathematics, vol. 125, no. 1-2, pp. 395–421, 2000.

A. Natarajan and N. Mohankumar, “A comparison of some quadrature methods for approximating Cauchy principal value integrals,” Journal of Computational Physics, vol. 116, no. 2, pp. 365–368, 1995.

F. Stenger, “Sinc approximation of Cauchy-type integrals over arcs,” in Proceedings of the Meeting in Honor of David Elliott's 65th Birthday, 1998.

C. F. Xu and Y. F. Yao, “Numerical calculation of Cauchy principal value integrals using wavelet functions,” Numerical Mathematics. A Journal of Chinese Universities, vol. 20, no. 2, pp. 97–105, 1998.

D. J. Liu, J. M. Wu, and D. H. Yu, “The superconvergence of the Newton-Cotes rule for Cauchy principal value integrals,” Journal of Computational and Applied Mathematics, vol. 235, no. 3, pp. 696–707, 2010.

J. M. Wu and W. W. Sun, “The superconvergence of Newton-Cotes rules for the Hadamard finite-part integral on an interval,” Numerische Mathematik, vol. 109, no. 1, pp. 143–165, 2008.

J. Li and D. H. Yu, “Error expansion of classical trapezoidal rule for computing Cauchy principal value integral,” Computer Modeling in Engineering & Sciences, vol. 93, no. 1, pp. 47–67, 2013.

J. Li, J. E. Yang, and D. H. Yu, “Error expansion of classical mid-point rectangle rule for computing Cauchy principal value integrals on an interval,” International Journal of Computer Mathematics, vol. 91, no. 10, pp. 2294–2306, 2014.

L. C. Andrews, Special Functions of Mathematics for Engineers, McGraw-Hill, 2nd edition, 1992.