- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Journal of Applied Mathematics
Volume 2013 (2013), Article ID 736834, 9 pages
The Modified Trapezoidal Rule for Computing Hypersingular Integral on Interval
1School of Science, Shandong Jianzhu University, Jinan 25010, China
2School 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.
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.
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 . 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  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  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  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.  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.
3. Proof of Main Results
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
Lemma 4. Assume that and ; there holds
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. ) 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.
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.
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).
- G. Monegato, “Numerical evaluation of hypersingular integrals,” Journal of Computational and Applied Mathematics, vol. 50, no. 1–3, pp. 9–31, 1994.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- D. H. Yu, Natural Boundary Integral Method and Its Applications, Kluwer Academic Publishers, 2002.
- P. Linz, “On the approximate computation of certain strongly singular integrals,” Computing, vol. 35, no. 3-4, pp. 345–353, 1985.
- D. H. Yu, “The approximate computation of hypersingular integrals on interval,” Numerical Mathematics, vol. 1, no. 1, pp. 114–127, 1992.
- 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.
- 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.
- L. C. Andrews, Special Functions of Mathematics for Engineers, McGraw-Hill, New York, NY, USA, 2nd edition, 1992.