The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle

Abstract

The classical composite rectangle (constant) rule for the computation of Cauchy principle value integral with the singular kernel  is discussed. We show that the superconvergence rate of the composite midpoint rule occurs at certain local coodinate of each subinterval and obtain the corresponding superconvergence error estimate. Then collation methods are presented to solve certain kind of Hilbert singular integral equation. At last, some numerical examples are provided to validate the theoretical analysis.

Share and Cite:

Li, J. , Gong, B. and Liu, W. (2016) The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle. American Journal of Computational Mathematics, 6, 98-107. doi: 10.4236/ajcm.2016.62011.

Received 8 March 2016; accepted 10 June 2016; published 13 June 2016

1. Introduction

Consider the Cauchy principle integral

(1)

where denotes a Cauchy principle value integral and s 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, direct definition and so on. In this paper we adopt the following one

(2)

Cauchy principal value integrals have recently attracted a lot of attention [1] - [5] . The main reason for this interest is probably due to the fact that integral equations with Cauchy principal value integrals have shown to be an adequate tool for the modeling of many physical situations, such as acoustics, fluid mechanics, elasticity, fracture mechanics and electromagnetic scattering problems and so on. It is the aim of this paper to investigate the superconvergence phenomenon of rectangle rule for it and, in particular, to derive error estimates.

The superconvergence of composite Newton-Cotes rules for Hadamard finite-part integrals was studied in [6] - [8] , where the superconvergence rate and the superconvergence point were presented, respectively. Lyness [9] derived the Euler-Maclaurin formula for Cauchy principal value integrals. Elliott and Venturino [2] employed sigmoidal transformations to obtain better approximation to Cauchy principal value integrals. In the reference Avram Sidi [10] [11] and [12] presented high-accuracy numerical quadrature methods for integrals of singular periodic functions. The classical Euler-Maclaurin summation formula [13] expressed the difference between a definite integral over and its approximation using the trapezoidal rule with step length as an asymptotic expansion in powers of h together with a remainder term.

The extrapolation method for the computation of Hadamard finite-part integrals on the interval and in a circle is studied in [14] and [15] which focus on the asymptotic expansion of error function. Based on the asymptotic expansion of the error functional, algorithm with theoretical analysis of the generalized extrapolation is given.

In this paper, the density function f(x) is replaced by the approximation function fC(x) while the singular kernel

is computed analysis in each subinterval, where fC(x) is the midpoint rectangle rule. This methods

may be considered as the semi-discrete methods and the order of singularity kernel can be reduced somehow. This idea was firstly presented by Linz [16] in the paper to calculated the hypersingular integral on interval. He used the trapezoidal rule and Simpson rule to approximate the density function f(x) and the convergence rate was when the singular point was always located at the middle of certain subinterval. This paper focuses on the superconvergence of mid-rectangle rule for Cauchy principle integrals. We prove both theoreti- cally and numerically that the composite mid-rectangle rule reaches the superconvergence rate when the local coordinate of the singular point s is. Then a collation methods is presented to solve certain kind of Hilbert singular integral equation.

The rest of this paper is organized as follows. In Sect. 2, after introducing some basic formulas of the rectangle rule, we present the main resluts. In Sect. 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 constant interpolant for

(3)

and a linear transformation

(4)

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

(5)

where denotes the error functional and is the Cote coefficients given by

(6)

We also define

(7)

Theorem 1: Assume. For the rectangle rule defined as (5). Assume that, there exist a positive constant C, independent of h and s, such that

(8)

where

(9)

Proof: Let, then we have As

(10)

For the first part of (10), we have

(11)

For the second part of (10),

(12)

(13)

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

Setting

(14)

Lemma 1: Assume with. Let be defined by (14), then there holds that

(15)

Proof: For, by the definition of cauchy principal value integral, we have

(16)

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

(17)

Now, by using the well-known identity

(18)

and

(19)

The proof is completed.

By the identity in [17]

(20)

then we get

(21)

and

(22)

For, by the definition of cauchy principal value integral, we have

(23)

Let be the function of the second kind associated with the Legendre polynomial, defined by (cf. [17] )

(24)

We also define

(25)

Then, by the definition of W,

it follows that

(26)

Theorem 2: Assume. For the rectangle rule defined as (5). Assume that, there exist a positive constant C, independent of h and s, such that

(27)

where

(28)

is defined as (9).

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

3. Proof of the Theorem

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

Preliminaries

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

Lemma 2: Under the same assumptions of Theorem 2, it holds that

where.

Proof: Performing Taylor expansion of at the point x, we have

(29)

and

(30)

Combining (29) and (30) together we get the results.

Proof of Theorem 2: we have

(31)

For, we have

(32)

Putting (31) and (32) together yields

(33)

Here

with the linear transformation from to the identity interval. As for the last part of

which can be considered as the error estimate of left rectangle rule for the definite integral. Obviously,by the Theorem, it can be expanded by the Euler-Maclaurin expansions and we have

(34)

It is easy to see that there are not relation with the singular point which can be written as

(35)

The proof is complete.

We actually obtain the error expansion of the rectangle rule and moreover, get the explicit expression of the first order term. So it is easy for us to get the superconvergence point with, which means that is the superconvergence point in subinterval not near the end of the interval.

Based on the theorem 1, we present the modify rectangle rule

(36)

4. Numerical Example

In this section, computational results are reported.

Example 1: We consider the Hilbert singular integral with. with is the superconvergence point.

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

For the modify classical rectangle rule, from Table 3 and Table 4, for the non-superconvergence point and the supersonvergence point, we all get the supercocergence phenomenon.

In this section, we consider the integral equation

(37)

with the compatibility condition

(38)

As in [5] , under the condition of (38), there exists a unique solution for the integral Equation (37). In order to get a unique solution, we adopt the following condition

(39)

By choosing the middle points, we get the composite rectangle rule to approximate the Hilbert singular integral in (37), then the following linear system is obtained

(40)

Table 1. An error estimate of the rectangle rule.

Table 2. An error estimate of the rectangle rule.

Table 3. An error estimate of the modify rectangle rule.

Table 4. An error estimate of the modify rectangle rule.

and written as the matrix expression as

(41)

where

(42)

here denotes the numerical solution of f at. By directly calculation, we get that is not only a symmetric Toeplitz matrix but also a circulant matrix. As for any,

(43)

from (43), we know that is singular matrix, then we cannot use system (40) or (41) to solve the integral Equation (37).

In order to get a well-conditioned definite system, we introduce a regularizing factor in (40), which leads to linear system

(44)

where defined by

(45)

Then the matrix form of system (44) can be presented as

(46)

where

(47)

and.

Example 2: Now we consider an example of solving Hilbert integral equation by collocation scheme. Let, the exact solution is.

We examine the maximal nodal error, defined by

(48)

where denotes the approximation of at. Numerical results presented in Table 5 show that both the maximal nodal errors are as follow.

Acknowledgements

The work of Jin Li was supported by National Natural Science Foundation of China (Grant No. 11471195 and

Table 5. Errors for the solution of the Hilbert integral equation of first kind.

Grant No. 11101247), China Postdoctoral Science Foundation (Grant No. 2013M540541 and 2015T80703). The work of Wei Liu was supported by National Natural Science Foundation of China (Grant No. 11401289).

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Diethdm, K. (1995) Asymptotically Sharp Error Bounds for a Quadrature Rule for Cauchy Principal Value Integrals Based on Piecewise Linear Interpolation, Approx. Theory of Probability and Its Applications, 11, 78-89. [2] Elliott, D. and Venturino, E. (1997) Sigmoidal Transformations and the Euler-Maclaurin Expansion for Evaluating Certain Hadamard Finite-Part Integrals. Numerische Mathematik, 77, 453-465. http://dx.doi.org/10.1007/s002110050295 [3] Hasegawa, T. (2004) Uniform Approximations to Finite Hilbert Transform and Its Derivative. Journal of Computational and Applied Mathematics, 163, 127-138. http://dx.doi.org/10.1016/j.cam.2003.08.059 [4] Ioakimidis, N.I. (1985) On the Uniform Convergence of Gaussian Quadrature Rules for Cauchy Principal Value Integrals and Their Derivatives. Mathematics of Computation, 44, 191-198. http://dx.doi.org/10.1090/S0025-5718-1985-0771040-8 [5] Yu, D.H. (2002) Natural Boundary Integrals Method and Its Applications. Kluwer Academic Publishers, Dordrecht, 50. [6] Yu, D.H. (1992) The Approximate Computation of Hypersingular Integrals on Interval. Numerical Mathematics: A Journal of Chinese Universities (English Series), 1, 114-127. [7] Zhang, X.P., Wu, J.M. and Yu, D.H. (2010) The Superconvergence of Composite Trapezoidal Rule for Hadamard Finite-Part Integral on a Circle and Its Application. International Journal of Computer Mathematics, 87, 855-876. http://dx.doi.org/10.1080/00207160802226517 [8] Wu, J.M. and Sun, W.W. (2008) The Superconvergence of Newton-Cotes Rules for the Hadamard Finite-Part Integral on an Interval. Numerische Mathematik, 109, 143-165. http://dx.doi.org/10.1007/s00211-007-0125-7 [9] Lyness, J.N. and Ninhan, B.W. (1967) Numerical Quadrature and Asymptotic Expansions. Mathematics of Computation, 21, 162-178. http://dx.doi.org/10.1090/S0025-5718-1967-0225488-X [10] Sidi, A. and Israeli, M. (1988) Quadrature Methods for Periodic Singular and Weakly Singular Fredholm Integral Equations. Journal of Scientific Computing, 3, 201-231. http://dx.doi.org/10.1007/BF01061258 [11] Sidi, A. (2003) Practical Extrapolation Methods: Theory and Applications. Cambridge University Press, Cambridge. http://dx.doi.org/10.1017/CBO9780511546815 [12] Zeng, G., Lei, L. and Huang, J. (2015) A New Construction of Quadrature Formulas for Cauchy Singular Integral. Journal of Computational Analysis and Applications, 17, 426-436. [13] Sidi, A. (2014) Analysis of Errors in Some Recent Numerical Quadrature Formulas for Periodic Singular and Hypersingular Integrals via Regularization. Applied Numerical Mathematics, 81, 30-39. http://dx.doi.org/10.1016/j.apnum.2014.02.011 [14] Li, J., Wu, J.M. and Yu, D.H. (2009) Generalized Extrapolation for Computation of Hypersingular Integrals in Boundary Element Methods. CMES: Computer Modeling in Engineering & Sciences, 42, 151-175. [15] Li, J., Zhang, X.P. and Yu, D.H. (2013) Extrapolation Methods to Compute Hypersingular Integral in Boundary Element Methods. Science China Mathematics, 56, 1647-1660. http://dx.doi.org/10.1007/s11425-013-4593-1 [16] Linz, P. (1985) On the Approximate Computation of Certain Strongly Singular Integrals. Computing, 35, 345-353. http://dx.doi.org/10.1007/BF02240199 [17] Andrews, L.C. (1992) Special Functions of Mathematics for Engineers. 2nd Edition, McGraw-Hill, Inc., Cambridge, 258.