PDF 52999.pdf | Request a copy Full text restricted Available under University of Tasmania Standard License.

Abstract

In the two-dimensional boundary element method, one often needs to evaluate numerically integrals of the form where j2 is a quadratic, g is a polynomial and f is a rational, logarithmic or algebraic function with a singularity at zero. The constants a and b are such that -1a1 and 0<b1 so that the singularities of f will be close to the interval of integration. In this case the direct application of Gauss–Legendre quadrature can give large truncation errors. By making the transformation x=a+bsinh(μu-η), where the constants μ and η are chosen so that the interval of integration is again [-1,1], it is found that the truncation errors arising, when the same Gauss–Legendre quadrature is applied to the transformed integral, are much reduced. The asymptotic error analysis for Gauss–Legendre quadrature, as given by Donaldson and Elliott [A unified approach to quadrature rules with asymptotic estimates of their remainders, SIAM J. Numer. Anal. 9 (1972) 573–602], is then used to explain this phenomenon and justify the transformation.

Item Type:	Article
Authors/Creators:	Elliott, D and Johnston, PR
Keywords:	Numerical integration; Boundary element method; Nearly singular integrals; Nonlinear coordinate transformation; Sinh function
Journal or Publication Title:	Journal of Computational and Applied Mathematics
Publisher:	Elsevier Science BV
ISSN:	0377-0427
DOI / ID Number:	https://doi.org/10.1016/j.cam.2006.03.012
Item Statistics:	View statistics for this item

Actions (login required)

Item Control Page