A new method for the numerical evaluation of nearly singular integrals on triangular elements in the 3D boundary element method

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:245

Article

Johnston, Barbara M.¹  Johnston, Peter R.¹  Elliott, David² 

[1] Griffith Univ, Sch Biomol & Phys Sci, Nathan, Qld 4111, Australia

[2] Univ Tasmania, Sch Math & Phys, Hobart, Tas 7001, Australia

关键词: Numerical integration;    Nearly singular integrals;    Boundary element method;    Sinh transformation;   

DOI  :  10.1016/j.cam.2012.12.018

来源: Elsevier

【 摘 要 】

A new method (the sinh-sigmoidal method) is proposed for the numerical evaluation of both nearly weakly and nearly strongly singular integrals on triangular boundary elements. These integrals arise in the 3D boundary element method when the source point is very close to the element of integration. The new polar coordinate-based method introduces a sinh transformation in the radial direction and a sigmoidal transformation in the angular direction, before the application of Gaussian quadrature. It also uses approximately twice as many quadrature points in the angular direction as in the radial direction, in response to a finding that the evaluation of these types of integrals is particularly sensitive to the placement of the quadrature points in the angular direction. Comparisons with various other methods demonstrate its accuracy and competitiveness. A major advantage of the new method is its ease of implementation and applicability to a wide class of integrals. (C) 2013 Elsevier B.V. All rights reserved.

【 授权许可】

Free