Journal of Applied Mathematics

Volume 2013, Article ID 736834, 9 pages

http://dx.doi.org/10.1155/2013/736834

## The Modified Trapezoidal Rule for Computing Hypersingular Integral on Interval

^{1}School of Science, Shandong Jianzhu University, Jinan 25010, China^{2}School of Mathematics, Shandong University, Jinan 250100, China

Received 3 August 2013; Revised 12 October 2013; Accepted 12 October 2013

Academic Editor: Kazutake Komori

Copyright © 2013 Jin Li and Xiuzhen 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 modified trapezoidal rule for the computation of hypersingular integrals in boundary element methods is discussed. When the special function of the error functional equals zero, the convergence rate is one order higher than the general case. A new quadrature rule is presented and the asymptotic expansion of error function is obtained. Based on the error expansion, not only do we obtain a high order of accuracy, but also a posteriori error estimate is conveniently derived. Some numerical results are also reported to confirm the theoretical results and show the efficiency of the algorithms.

#### 1. Introduction

Consider the following integral: where denotes a Hadamard finite-part integral ( is called hypersingular integral and is called supersingular integral) and is the singular point. The formulation of these classes of boundary value problems in terms of hypersingular integral equations has drawn lots of interest. Many scientific and engineering problems, such as acoustics, electromagnetic scattering, and fracture mechanics, can be reduced to boundary integral equations with hypersingular kernels. There exist several definitions, equivalent mathematically, for such kind of integrals in some literatures.

We mention the following one:

Accurate calculation of boundary element methods (BEM) arising in boundary integral equations has been a subject of intensive research in recent years. The hypersingular integrals have certain properties different from regular and weak singular integrals. One of the major problems arising from boundary element method, for solving such integral equations, is how to evaluate the hypersingular integrals on the interval or on the circle efficiently.

Hypersingular integral must be considered in Hadamard finite-part sense. Numerous works [1–18] have been devoted towards developing efficient quadrature formulas. In 1983, the series expansion of hypersingular integral kernel on circle was firstly suggested by Yu [19]. He solved the harmonic and biharmonic natural boundary integral equations successfully. The Newton-Cotes methods to compute the hypersingular integral on interval were firstly studied by Linz [20] with generalized trapezoidal and Simpson rules which fail altogether when the singular point is close to a mesh point. In order to make the mesh be selected in such a way that falls near the center of a subinterval, two shorter subintervals at the end of the interval were allowed. Then Yu [21] gave new quadrature formulae to compute the case of singular point coinciding with the mesh point which presented that the error estimate is . In 1999, Wu and Yu [22] presented simple, easy to be implemented methods not affected by the location of singular point with calculation of double. In recent years, the case of singular point coincided with the mesh point, and Wu et al. [23] presented a modified trapezoidal rule and proved the convergence rate.

In this paper, for the case of singular point coinciding with the mesh point a new quadrature rule is introduced. Based on the expansion of the error functional, the error estimate is presented and a posteriori error estimate is given. Then not only do we obtain a high order of accuracy, but also a posteriori error estimate is conveniently derived.

The rest of this paper is organized as follows. In Section 2, after introducing some basic formulas of the general (composite) trapezoidal rule and notations, we present our main result. In Section 3, the corresponding theoretical analysis is given. Finally, several numerical examples are given to validate our analysis.

#### 2. Main Result

Let be a uniform partition of the interval with mesh size and set then we get the new partition:

We define , the modified trapezoidal interpolation for , as and a linear transformation from the reference element to the subinterval . For the two subintervals and near the end of the interval , values in and , respectively.

Replacing in (1) with gives the new composite trapezoidal rule: where is the Cotes coefficients: , denotes the Kronecker delta, and denotes the error functional.

Theorem 1. *Assume , . For the trapezoidal rule defined in (7), there exists a positive constant , independent of and , such that
**
where
*

*Proof. *Let ; then we have , as
For the first part of (11), since , by Taylor expansion, we have

By the definition of finite-part integral,
we have
Now, we estimate the right hand side of (14) term by term. Since , we have
For the second part of (11),

Combining (14), (15), and (16) leads to (9) and the proof is completed.

Firstly, we set then we have By straight calculation, we get Let and the operator be defined as Obviously, the operator is linear operator. Then we set

Now we present our main results below. The proof will be given in the next section.

Theorem 2. *Assume . For the trapezoidal rule defined in (7), there exists a positive constant , independent of and , such that
**
where , , and
**
( defined as (10)) and
*

#### 3. Proof of Main Results

##### 3.1. Preliminaries

In the following section, denotes certain constant independent of and , and its value varies with places.

Lemma 3. *Assume that and are defined by (5); there holds
**
where
**,,, and
*

*Proof. *Taking Taylor expansion for , at , there holds
by applying Taylor expansion to in (28) at , we have
Therefore, (25) can be obtained directly from (28) and (29).

Lemma 4. *Assume that and ; there holds
*

*Proof. *According to (2) and the linear transformation (6) for , we have
where we have used . The second identity can be similarly obtained.

Lemma 5. *Suppose and are defined by (18) and (24), respectively; then one has
**
where denotes that the first interval is certain part of the reference element.*

*Proof. *By straightforward calculation, we have
Noting that , we have and
On the other hand, since , we have and
Combining (34) and (35), we get (32).

Set

Lemma 6. *Under the same assumptions of Theorem 2, for in (36), there holds that
**
where is defined in (10).*

*Proof. *Since , by Taylor expansion, we have

By the definition of finite-part integral (13), we have
Now, we estimate the right hand side of (39) term by term. Since , we have

Combining (40), (41), and (42) leads to (37) and the proof is completed.

*Proof of Theorem 2. *Consider

By the definition of , we have

Putting (43) and (44) together, we have
where
For the first part of ,
where and . We have also used the identity = , , and .

For the second part of ,
By Lemmas 5 and 6, we have
Then the proof is completed.

##### 3.2. The Calculation of

Let be the function of the second kind associated with the Legendre polynomial , defined by (cf. [24]) We also define Then, by the definition of , It follows that which means

What remains is to determine the constant . By using the identities (cf. [24, Chapter 1, Section ]), where denote the Bernoulli numbers, we have Setting gives Then we have where we have used the formulae (cf. [24, Chapter 1, Section ]) Then we have

By Theorem 2, we get the following error expansion:

Let where is the starting meshes, is the refining numbers, and . Then we have the following.

Corollary 7. *Under the same assumption of Theorem 2 and (62), there holds
*

#### 4. Numerical Examples

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

*Example 1. *Consider the hypersingular integral

We examine the dynamic point , in Table 1 show that when the local coordinate of singular point , the quadrature reach the convergence rate of as for the nonsupersingular point, there are no convergence rate which agree with our theorematically analysis. For the case of , Table 2 shows that there is no superconvergence phenomenon because of the influence of which coincides with our theoretical analysis.

*Example 2. *Consider the hypersingular integral

The numerical results show that the convergence rate reaches when the singular point coincides with the mesh point in Table 3. In Table 4, a posteriori error estimate is presented and the convergence rate is also which agrees with our theoretical analysis.

*Example 3. *Now we consider an example of less regularity. Let , , and
Obviously, . The exact value of the integral is

The numerical results are presented in Tables 5 and 6. When the density function is smooth enough , the error bound is , and if the density function has less regularity , there is no hyperconvergence phenomenon, which means the regularity of density function cannot be reduced.

#### Acknowledgments

The work of Jin Li was supported by the National Natural Science Foundation of China (nos. 11101247, 11201209, and 11101246), China Postdoctoral Science Foundation Fund Project (no. 2013M540541), the Shandong Provincial Natural Science Foundation of China (no. ZR2011AQ020), and a Project of Shandong Province Higher Educational Science and Technology Program (no. J11LE8).

#### References

- G. Monegato, “Numerical evaluation of hypersingular integrals,”
*Journal of Computational and Applied Mathematics*, vol. 50, no. 1–3, pp. 9–31, 1994. View at Publisher · View at Google Scholar · View at MathSciNet - U. J. Choi, S. W. Kim, and B. I. Yun, “Improvement of the asymptotic behaviour of the Euler-Maclaurin formula for Cauchy principal value and Hadamard finite-part integrals,”
*International Journal for Numerical Methods in Engineering*, vol. 61, no. 4, pp. 496–513, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. T. Chen and H. K. Hong, “Review of dual boundary element methods with emphasis on hypersingular integrals and divergent series,”
*Applied Mechanics Reviews*, vol. 52, no. 1, pp. 17–32, 1999. View at Google Scholar · View at Scopus - A. Frangi and M. Guiggiani, “A direct approach for boundary integral equations with high-order singularities,”
*International Journal for Numerical Methods in Engineering*, vol. 49, no. 7, pp. 871–898, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - C. Y. Hui and D. Shia, “Evaluations of hypersingular integrals using Gaussian quadrature,”
*International Journal for Numerical Methods in Engineering*, vol. 44, no. 2, pp. 205–214, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Kim and U. J. Choi, “Two trigonometric quadrature formulae for evaluating hypersingular integrals,”
*International Journal for Numerical Methods in Engineering*, vol. 56, no. 3, pp. 469–486, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li, J. M. Wu, and D. H. Yu, “Generalized extrapolation for computation of hypersingular integrals in boundary element methods,”
*Computer Modeling in Engineering & Sciences*, vol. 42, no. 2, pp. 151–175, 2009. View at Google Scholar · View at MathSciNet - J. Li and D. H. Yu, “The superconvergence of certain two-dimensional Cauchy principal value integrals,”
*Computer Modeling in Engineering & Sciences*, vol. 71, no. 4, pp. 331–346, 2011. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Google Scholar · View at MathSciNet - J. Li and D. H. Yu, “The error estimate of Newton-Cotes methods to compute hypersingular integrals,”
*Mathematica Numerica Sinica*, vol. 33, no. 1, pp. 77–86, 2011. View at Google Scholar · View at MathSciNet - J. Li, X. Zhang, and D. Yu, “Superconvergence and ultraconvergence of Newton-Cotes rules for supersingular integrals,”
*Journal of Computational and Applied Mathematics*, vol. 233, no. 11, pp. 2841–2854, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li, X. Zhang, and D. Yu, “Extrapolation methods to compute hypersingular integral in boundary element methods,”
*Science China Mathematics*, vol. 56, no. 8, pp. 1647–1660, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. N. J. Moore, L. J. Gray, and T. Kaplan, “Evaluation of supersingular integrals: second-order boundary derivatives,”
*International Journal for Numerical Methods in Engineering*, vol. 69, no. 9, pp. 1930–1947, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. K. Du, “Evaluations of certain hypersingular integrals on interval,”
*International Journal for Numerical Methods in Engineering*, vol. 51, no. 10, pp. 1195–1210, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Wu and W. Sun, “The superconvergence of the composite trapezoidal rule for Hadamard finite part integrals,”
*Numerische Mathematik*, vol. 102, no. 2, pp. 343–363, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Zhang, J. Wu, and D. Yu, “Superconvergence of the composite Simpson's rule for a certain finite-part integral and its applications,”
*Journal of Computational and Applied Mathematics*, vol. 223, no. 2, pp. 598–613, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. H. Yu,
*Natural Boundary Integral Method and Its Applications*, Kluwer Academic Publishers, 2002. View at MathSciNet - P. Linz, “On the approximate computation of certain strongly singular integrals,”
*Computing*, vol. 35, no. 3-4, pp. 345–353, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. H. Yu, “The approximate computation of hypersingular integrals on interval,”
*Numerical Mathematics*, vol. 1, no. 1, pp. 114–127, 1992. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. M. Wu and D. H. Yu, “An approximate computation of hypersingular integrals on an interval,”
*Journal on Numerical Methods and Computer Applications*, vol. 19, no. 2, pp. 118–126, 1998. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Wu, Y. Wang, W. Li, and W. Sun, “Toeplitz-type approximations to the Hadamard integral operator and their applications to electromagnetic cavity problems,”
*Applied Numerical Mathematics*, vol. 58, no. 2, pp. 101–121, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. C. Andrews,
*Special Functions of Mathematics for Engineers*, McGraw-Hill, New York, NY, USA, 2nd edition, 1992. View at MathSciNet