**American Journal of Computational Mathematics**

Vol.06 No.02(2016), Article ID:67263,10 pages

10.4236/ajcm.2016.62011

The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle

Jin Li^{1,2}, Benxue Gong^{3}, Wei Liu^{4}

^{1}School of Science, Shandong Jianzhu University, Jinan, China

^{2}School of Mathematics, Shandong University, Jinan, China

^{3}School of Science, Shandong University of Technology, Zibo, China

^{4}School of Mathematics and Statistics Science, Ludong University, Yantai, China

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

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

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.

**Keywords:**

Cauchy Principal Value Integral, Extrapolation Method, Composite Rectangle Rule, Superconvergence, Error Expansion

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 f_{C}(x) while the singular kernel

is computed analysis in each subinterval, where f_{C}(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).

Cite this paper

Lianju Chen,Xue Su,Yuhong Li,Jin Li,1 1,Benxue Gong,Wei Liu, (2016) The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle. *American Journal of Computational Mathematics*,**06**,98-107. doi: 10.4236/ajcm.2016.62011

References

- 1. Ewing, R., Weber, W. and Clinard, F. (1995) Radiation Effects in Nuclear Waste Forms for High-Level Radioactive Waste. Progress in Nuclear Energy, 29, 63-127.

http://dx.doi.org/10.1016/0149-1970(94)00016-Y - 2. Ewing, C., Weber, W. and Lian, J. (2004) Nuclear Waste Disposal—Pyrochlore (A2B2O7): Nuclear Waste Form for the Immobilization of Plutonium and “Minor” Actinides. Journal of Applied Physics, 95, 5949-5971.

http://dx.doi.org/10.1063/1.1707213 - 3. Wang, S.X., Begg, B.D., Wang, L.M., Ewing, R.C., Weber, W.J. and Govidan Kutty, K.V. (1999) Radiation Stability of Gadolinium Zirconate: A Waste Form for Plutonium Disposition. Journal of Materials Research, 14, 4470-4473.

http://dx.doi.org/10.1557/JMR.1999.0606 - 4. Sickafus, K.E., Grimes, R.W., Valdez, J.A., Cleave, A., Tang, M., Ishimaru, M., Corish, S.M., Stanek, C.R. and Uberuaga, B.P. (2007) Radiation-Induced Amorphization Resistance and Radiation Tolerance in Structurally Related Oxides. Nature Materials, 6, 217-223.

http://dx.doi.org/10.1038/nmat1842 - 5. Begg, B.D., Hess, N.J., Weber, W.J., Devanathan, R., Icenhower, J.P., Thevuthasan, S. and McGrail, B.P. (2001) Heavy-Ion Irradiation Effects on Structures and Acid Dissolution of Pyrochlores. Journal of Nuclear Materials, 288, 208-216.

http://dx.doi.org/10.1016/S0022-3115(00)00708-X - 6. Vance, E.R., Lumpkin, G.R., Carter, M.L., Cassidy, D.J., Ball, C.J., Day, R.A. and Begg, B.D. (2002) Incorporation of Uranium in Zirconolite (CaZrTi2O7). Journal of the American Ceramic Society, 85, 1853-1859.

http://dx.doi.org/10.1111/j.1151-2916.2002.tb00364.x - 7. Lian, J., Chen, J., Wang, L., Ewing, R., Farmer, J., Boatner, L. and Helean, K. (2003) Radiation-Induced Amorphization of Rare-Earth Titanate Pyrochlores. Physical Review B, 68, Article ID: 134107.

http://dx.doi.org/10.1103/PhysRevB.68.134107 - 8. Lian, J., Zu, X., Kutty, K., Chen, J., Wang, L. and Ewing, R. (2002) Ion-Irradiation-Induced Amorphization of La2Zr2O7 Pyrochlore. Physical Review B, 66, Article ID: 054108.

http://dx.doi.org/10.1103/PhysRevB.66.054108 - 9. Lian, J., Helean, K., Kennedy, B., Wang, L., Navrotsky, A. and Ewing, R. (2006) Effect of Structure and Thermodynamic Stability on the Response of Lanthanide Stannate Pyrochlores to Ion Beam Irradiation. Journal of Physical Chemistry B, 110, 2343-2350.

http://dx.doi.org/10.1021/jp055266c - 10. Jiang, C., Stanek, C., Sickafus, K. and Uberuaga, B. (2009) First-Principles Prediction of Disordering Tendencies in Pyrochlore Oxides. Physical Review B, 79, Article ID: 104203.

http://dx.doi.org/10.1103/PhysRevB.79.104203 - 11. Marks, N.A., Thomas, B.S., Smith, K.L. and Lumpkin, G.R. (2008) Thermal Spike Recrystallisation: Molecular Dynamics Simulation of Radiation Damage in Polymorphs of Titania. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 266, 2665-2670.

http://dx.doi.org/10.1016/j.nimb.2008.03.098 - 12. Weber, W., Ewing, R., Catlow, C., Rubia, T.D., Hobbs, L., Kinoshita, C., Matzke, H., Motta, A., Nastasi, M. and Salje, E. (1998) Radiation Effects in Crystalline Ceramics for the Immobilization of High-Level Nuclear Waste and Plutonium. Journal of Materials Research, 13, 1434-1484.

http://dx.doi.org/10.1557/JMR.1998.0205 - 13. Chen, Z.J., Xiao, H.Y., Zu, X.T. and Gao, F. (2008) First-Principles Calculation of Defect Formation Energies and Electronic Properties in Stannate Pyrochlores. Journal of Applied Physics, 104, Article ID: 093702.
- 14. Xiao, H.Y., Zu, X.T., Gao, F. and Weber, W.J. (2008) First-Principles Study of Energetic and Electronic Properties of A2Ti2O7 (A = Sm, Gd, Er) Pyrochlore. Journal of Applied Physics, 104, Article ID: 073503.

http://dx.doi.org/10.1063/1.2986156 - 15. Zhang, Z., Xiao, H., Zu, X.T., Gao, F. and Weber, W.J. (2009) First-Principles Calculation of Structural and Energetic Properties for A2Ti2O7 (A = Lu, Er, Y, Gd, Sm, Nd, La). Journal of Materials Research, 24, 1335-1341.

http://dx.doi.org/10.1557/jmr.2009.0152 - 16. Panero, W.R., Stixrude, L. and Ewing, R. (2004) First-Principles Calculation of Defect-Formation Energies in the Y2(Ti,Sn,Zr)2O7 Pyrochlore. Physical Review B, 70, Article ID: 054110.

http://dx.doi.org/10.1103/PhysRevB.70.054110 - 17. Li, N., Xiao, H.Y., Zu, X.T., Wang, L.M., Ewing, R.C., Lian, J. and Gao, F. (2007) First-Principles Study of Electronic Properties of La2Hf2O7 and Gd2Hf2O7. Journal of Applied Physics, 102, Article ID: 063704.

http://dx.doi.org/10.1063/1.2779262 - 18. Minervini, L., Grimes, R.W. and Sickafus, K.E. (2000) Disorder in Pyrochlore Oxides. Journal of the American Ceramic Society, 83, 1873-1878.

http://dx.doi.org/10.1111/j.1151-2916.2000.tb01484.x - 19. Lian, J., Weber, W.J., Jiang, W., Wang, L.M., Boatner, L.A. and Ewing, R.C. (2006) Radiation-Induced Effects in Pyrochlores and Nanoscale Materials Engineering. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 250, 128-136.

http://dx.doi.org/10.1016/j.nimb.2006.04.157 - 20. Sickafus, K., Minervini, L., Grimes, R.W., Valdez, J.A., Ishimaru, M., Li, F., McClellan, K.J. and Hartmann, T. (2000) Radiation Tolerance of Complex Oxides. Science, 289, 748-751.

http://dx.doi.org/10.1126/science.289.5480.748 - 21. Lumpkin, G.R., Smith, K.L. and Blackford, M.G. (2001) Heavy Ion Irradiation Studies of Columbite, Brannerite, and Pyrochlore Structure Types. Journal of Nuclear Materials, 289, 177-187.

http://dx.doi.org/10.1016/S0022-3115(00)00695-4 - 22. Lumpkin, G.R., Pruneda, M., Rios, S., Smith, K.L., Trachenko, K., Whittle, K.R. and Zaluzec, N.J. (2007) Nature of the Chemical Bond and Prediction of Radiation Tolerance in Pyrochlore and Defect Fluorite Compounds. Journal of Solid State Chemistry, 180, 1512-1518.

http://dx.doi.org/10.1016/j.jssc.2007.01.028 - 23. Xiao, H.Y., Gao, F. and Weber, W.J. (2010) Threshold Displacement Energies and Defect Formation Energies in Y2Ti2O7. Journal of Physics: Condensed Matter, 22, Article ID: 415801.

http://dx.doi.org/10.1088/0953-8984/22/41/415801 - 24. Li, Y.H., Uberuaga, B., Jiang, C., Choudhury, S., Valdez, J., Patel, M., Won, J., Wang, Y.Q., Tang, M., Safarik, D., Byler, D., McClellan, K., Usov, I.O., Hartmann, T., Baldinozzi, G. and Sickafus, K.E. (2012) Role of Antisite Disorder on Preamorphization Swelling in Titanate Pyrochlores. Physical Review Letters, 108, Article ID: 195504.

http://dx.doi.org/10.1103/PhysRevLett.108.195504 - 25. Kresse, G. and Hafner, J. (1993) Ab Initio Molecular Dynamics for Liquid Metals. Physical Review B, 47, 558-561.

http://dx.doi.org/10.1103/PhysRevB.47.558 - 26. Kresse, G., Hafner, J. and Needs, R. (1992) Optimized Norm-Conserving Pseudopotentials. Journal of Physics: Condensed Matter, 4, 7451.

http://dx.doi.org/10.1088/0953-8984/4/36/018 - 27. Kresse, G. and Furthmüller, J. (1996) Efficiency of ab Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Computational Materials Science, 6, 15-50.

http://dx.doi.org/10.1016/0927-0256(96)00008-0 - 28. Kresse, G. and Furthmüller, J. (1996) Efficient Iterative Schemes for ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Physical Review B, 54, Article ID: 11169.

http://dx.doi.org/10.1103/PhysRevB.54.11169 - 29. Blöchl, P.E. (1994) Projector Augmented-Wave Method. Physical Review B, 50, Article ID: 17953.

http://dx.doi.org/10.1103/PhysRevB.50.17953 - 30. Kresse, G. and Joubert, D. (1999) From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Physical Review B, 59, 1758.

http://dx.doi.org/10.1103/PhysRevB.59.1758 - 31. Perdew, J.P., Chevary, J., Vosko, S., Jackson, K.A., Pederson, M.R., Singh, D. and Fiolhais, C. (1992) Atoms, Molecules, Solids, and Surfaces: Applications of the Generalized Gradient Approximation for Exchange and Correlation. Physical Review B, 46, 6671.

http://dx.doi.org/10.1103/PhysRevB.46.6671 - 32. Lian, J., Wang, L., Wang, S., Chen, J., Boatner, L. and Ewing, R. (2001) Nanoscale Manipulation of Pyrochlore: New Nanocomposite Ionic Conductors. Physical Review Letters, 87, Article ID: 145901.

http://dx.doi.org/10.1103/PhysRevLett.87.145901 - 33. Subramanian, M., Aravamudan, G. and Subba Rao, G. (1983) Oxide Pyrochlores—A Review. Progress in Solid State Chemistry, 15, 55-143.

http://dx.doi.org/10.1016/0079-6786(83)90001-8 - 34. Kennedy, B.J. (1997) Structural Trends in Pyrochlore-Type Oxides. Physica B: Condensed Matter, 241-243, 303-310.

http://dx.doi.org/10.1016/S0921-4526(97)00570-X - 35. Brisse, F. and Knop, O. (1968) Pyrochlores. III. X-Ray, Neutron, Infrared, and Dielectric Studies of A2Sn2O7 Stannates. Canadian Journal of Chemistry, 46, 859-873.

http://dx.doi.org/10.1139/v68-148 - 36. Chen, Z.J., Xiao, H., Zu, X., Wang, L., Gao, F., Lian, J. and Ewing, R. (2008) Structural and Bonding Properties of Stannate Pyrochlores: A Density Functional Theory Investigation. Computational Materials Science, 42, 653-658.

http://dx.doi.org/10.1016/j.commatsci.2007.09.019 - 37. Xiao, H.Y., Wang, L.M., Zu, X.T., Lian, J. and Ewing, R.C. (2007) Theoretical Investigation of Structural, Energetic and Electronic Properties of Titanate Pyrochlores. Journal of Physics: Condensed Matter, 19, Article ID: 346203.

http://dx.doi.org/10.1088/0953-8984/19/34/346203 - 38. 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.
- 39. 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 - 40. 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 - 41. 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 - 42. Yu, D.H. (2002) Natural Boundary Integrals Method and Its Applications. Kluwer Academic Publishers, Dordrecht, 50.
- 43. Yu, D.H. (1992) The Approximate Computation of Hypersingular Integrals on Interval. Numerical Mathematics: A Journal of Chinese Universities (English Series), 1, 114-127.
- 44. 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 - 45. 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 - 46. 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 - 47. 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 - 48. Sidi, A. (2003) Practical Extrapolation Methods: Theory and Applications. Cambridge University Press, Cambridge.

http://dx.doi.org/10.1017/CBO9780511546815 - 49. 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.
- 50. 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 - 51. 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.
- 52. 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 - 53. Linz, P. (1985) On the Approximate Computation of Certain Strongly Singular Integrals. Computing, 35, 345-353.

http://dx.doi.org/10.1007/BF02240199 - 54. Andrews, L.C. (1992) Special Functions of Mathematics for Engineers. 2nd Edition, McGraw-Hill, Inc., Cambridge, 258.