Some quadrature formulae for the numerical evaluation of singular integrals of arbitrary order are established and both the estimate of remainder and the convergence of each quadrature formula derived here are also gi...Some quadrature formulae for the numerical evaluation of singular integrals of arbitrary order are established and both the estimate of remainder and the convergence of each quadrature formula derived here are also given.展开更多
In this paper, we introduce a quadrature rule for the numerical evaluation of certain hyper singular integrals. The rule is obtained by using Hermite interpolation polynomial. Error bound is also made.
Finite part integrals introduced by Hadamard in connection with hyperbolic partial differential equations,have been useful in a number of engineering applications.In this paper we investigate some numerical methods fo...Finite part integrals introduced by Hadamard in connection with hyperbolic partial differential equations,have been useful in a number of engineering applications.In this paper we investigate some numerical methods for computing finite-part integrals.展开更多
A new family of numerical integration formula is presented, which uses the function evaluation at the midpoint of the interval and odd derivatives at the endpoints. Because the weights for the odd derivatives sum to z...A new family of numerical integration formula is presented, which uses the function evaluation at the midpoint of the interval and odd derivatives at the endpoints. Because the weights for the odd derivatives sum to zero, the derivative calculations cancel out for the interior points in the composite form, so that these derivatives must only be calculated at the endpoints of the overall interval of integration. When using N subintervals, the basic rule which uses the midpoint function evaluation and the first derivative at the endpoints achieves fourth order accuracy for the cost of N/2 function evaluations and 2 derivative evaluations, whereas the three point open Newton-Cotes method uses 3N/4 function evaluations to achieve the same order of accuracy. These derivative-based midpoint quadrature methods are shown to be more computationally efficient than both the open and closed Newton-Cotes quadrature rules of the same order. This family of derivative-based midpoint quadrature rules are derived using the concept of precision, along with the error term. A theorem concerning the order of accuracy of quadrature rule using the concept of precision is provided to justify its use to determine the leading order error term.展开更多
In this paper, the author discusses some singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, and obtain some useful p...In this paper, the author discusses some singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, and obtain some useful properties for them. These results improve both the classical theory of singular integral equation and the classical theory of singular quadrature.展开更多
This study presents the determination of the stress intensity factors （SIFs） at the edges of the cracks in an elastic strip weakened by N-collinear cracks. The problem of an orthotropic elastic strip is reduced to a...This study presents the determination of the stress intensity factors （SIFs） at the edges of the cracks in an elastic strip weakened by N-collinear cracks. The problem of an orthotropic elastic strip is reduced to a system of Cauchy type singular integral equations. The system of singular integral equations is approached by a Quadrature technique. Under two different loading conditions, the results are obtained for the different cases of crack numbers. The resistance of the strip is examined by considering the orthotropic properties of the strip material. Finally, the crack interactions are clarified during the analysis.展开更多
The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element ...The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.展开更多
We address the evaluation of highly oscillatory integrals,with power-law and logarithmic singularities.Such problems arise in numerical methods in engineering.Notably,the evaluation of oscillatory integrals dominates ...We address the evaluation of highly oscillatory integrals,with power-law and logarithmic singularities.Such problems arise in numerical methods in engineering.Notably,the evaluation of oscillatory integrals dominates the run-time for wave-enriched boundary integral formulations for wave scattering,and many of these exhibit singularities.We show that the asymptotic behaviour of the integral depends on the integrand and its derivatives at the singular point of the integrand,the stationary points and the endpoints of the integral.A truncated asymptotic expansion achieves an error that decays faster for increasing frequency.Based on the asymptotic analysis,a Filon-type method is constructed to approximate the integral.Unlike an asymptotic expansion,the Filon method achieves high accuracy for both small and large frequency.Complex-valued quadrature involves interpolation at the zeros of polynomials orthogonal to a complex weight function.Numerical results indicate that the complex-valued Gaussian quadrature achieves the highest accuracy when the three methods are compared.However,while it achieves higher accuracy for the same number of function evaluations,it requires signi cant additional cost of computation of orthogonal polynomials and their zeros.展开更多
基金Supported by NNSF and RFDP of Higher Education of China.
文摘Some quadrature formulae for the numerical evaluation of singular integrals of arbitrary order are established and both the estimate of remainder and the convergence of each quadrature formula derived here are also given.
文摘In this paper, we introduce a quadrature rule for the numerical evaluation of certain hyper singular integrals. The rule is obtained by using Hermite interpolation polynomial. Error bound is also made.
文摘Finite part integrals introduced by Hadamard in connection with hyperbolic partial differential equations,have been useful in a number of engineering applications.In this paper we investigate some numerical methods for computing finite-part integrals.
文摘A new family of numerical integration formula is presented, which uses the function evaluation at the midpoint of the interval and odd derivatives at the endpoints. Because the weights for the odd derivatives sum to zero, the derivative calculations cancel out for the interior points in the composite form, so that these derivatives must only be calculated at the endpoints of the overall interval of integration. When using N subintervals, the basic rule which uses the midpoint function evaluation and the first derivative at the endpoints achieves fourth order accuracy for the cost of N/2 function evaluations and 2 derivative evaluations, whereas the three point open Newton-Cotes method uses 3N/4 function evaluations to achieve the same order of accuracy. These derivative-based midpoint quadrature methods are shown to be more computationally efficient than both the open and closed Newton-Cotes quadrature rules of the same order. This family of derivative-based midpoint quadrature rules are derived using the concept of precision, along with the error term. A theorem concerning the order of accuracy of quadrature rule using the concept of precision is provided to justify its use to determine the leading order error term.
文摘In this paper, the author discusses some singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, and obtain some useful properties for them. These results improve both the classical theory of singular integral equation and the classical theory of singular quadrature.
文摘This study presents the determination of the stress intensity factors （SIFs） at the edges of the cracks in an elastic strip weakened by N-collinear cracks. The problem of an orthotropic elastic strip is reduced to a system of Cauchy type singular integral equations. The system of singular integral equations is approached by a Quadrature technique. Under two different loading conditions, the results are obtained for the different cases of crack numbers. The resistance of the strip is examined by considering the orthotropic properties of the strip material. Finally, the crack interactions are clarified during the analysis.
基金supported by the National Natural Science Foundation of China(11304344,11404364)the Project of Hubei Provincial Department of Education(D20141803)+1 种基金the Natural Science Foundation of Hubei Province(2014CFB378)the Doctoral Scientific Research Foundation of Hubei University of Automotive Technology(BK201604)
文摘The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.
基金The work is supported by Royal Society International Exchanges(grant IE141214)the Projects of International Cooperation and Exchanges NSFC-RS(Grant No.11511130052)+1 种基金the Key Science and Technology Program of Shaanxi Province of China(Grant No.2016GY-080)the Fundamental Research Funds for the Central Universities.
文摘We address the evaluation of highly oscillatory integrals,with power-law and logarithmic singularities.Such problems arise in numerical methods in engineering.Notably,the evaluation of oscillatory integrals dominates the run-time for wave-enriched boundary integral formulations for wave scattering,and many of these exhibit singularities.We show that the asymptotic behaviour of the integral depends on the integrand and its derivatives at the singular point of the integrand,the stationary points and the endpoints of the integral.A truncated asymptotic expansion achieves an error that decays faster for increasing frequency.Based on the asymptotic analysis,a Filon-type method is constructed to approximate the integral.Unlike an asymptotic expansion,the Filon method achieves high accuracy for both small and large frequency.Complex-valued quadrature involves interpolation at the zeros of polynomials orthogonal to a complex weight function.Numerical results indicate that the complex-valued Gaussian quadrature achieves the highest accuracy when the three methods are compared.However,while it achieves higher accuracy for the same number of function evaluations,it requires signi cant additional cost of computation of orthogonal polynomials and their zeros.
