A sinh transformation for evaluating nearly singular boundary element integrals.

*(English)*Zbl 1119.65318Summary: An implementation of the boundary element method requires the accurate evaluation of many integrals. When the source point is far from the boundary element under consideration, a straightforward application of Gaussian quadrature suffices to evaluate such integrals. When the source point is on the element, the integrand becomes singular and accurate evaluation can be obtained using the same Gaussian points transformed under a polynomial transformation which has zero Jacobian at the singular point. A class of integrals which lies between these two extremes is that of nearly singular integrals. Here, the source point is close to, but not on, the element and the integrand remains finite at all points. However, instead of remaining flat, the integrand develops a sharp peak as the source point moves closer to the element, thus rendering accurate evaluation of the integral difficult.

This paper presents a transformation, based on the sinh function, which automatically takes into account the position of the projection of the source point onto the element, which we call the nearly singular point, and the distance from the source point to the element. The transformation again clusters the points towards the nearly singular point, but does not have a zero Jacobian. Implementation of the transformation is straightforward and could easily be included in existing boundary element method software. It is shown that, for the two-dimensional boundary element method, several orders of magnitude improvement in relative error can be obtained using this transformation compared to a conventional implementation of Gaussian quadrature. Asymptotic estimates for the truncation errors are also quoted.

This paper presents a transformation, based on the sinh function, which automatically takes into account the position of the projection of the source point onto the element, which we call the nearly singular point, and the distance from the source point to the element. The transformation again clusters the points towards the nearly singular point, but does not have a zero Jacobian. Implementation of the transformation is straightforward and could easily be included in existing boundary element method software. It is shown that, for the two-dimensional boundary element method, several orders of magnitude improvement in relative error can be obtained using this transformation compared to a conventional implementation of Gaussian quadrature. Asymptotic estimates for the truncation errors are also quoted.

##### MSC:

65D32 | Numerical quadrature and cubature formulas |

41A55 | Approximate quadratures |

65N38 | Boundary element methods for boundary value problems involving PDEs |

##### Keywords:

non-linear co-ordinate transformation; boundary element method; nearly singular integrals; numerical integration; sinh function; error estimates; Gaussian quadrature
PDF
BibTeX
XML
Cite

\textit{P. R. Johnston} and \textit{D. Elliott}, Int. J. Numer. Methods Eng. 62, No. 4, 564--578 (2005; Zbl 1119.65318)

Full Text:
DOI

##### References:

[1] | Boundary Elements An Introductory Course (2nd edn). Computational Mechanics Publications: Southampton, 1992. · Zbl 0780.73002 |

[2] | Telles, International Journal for Numerical Methods in Engineering 24 pp 959– (1987) |

[3] | Sanz-Serna, Communications in Numerical Methods in Engineering 6 pp 121– (1990) |

[4] | Telles, Engineering Analysis with Boundary Elements 13 pp 135– (1994) |

[5] | Johnston, Journal of Computational and Applied Mathematics 131 pp 223– (2001) |

[6] | Yun, SIAM Journal for Scientific Computing 24 pp 1203– (2003) |

[7] | Yun, International Journal for Numerical Methods in Engineering 57 pp 1883– (2003) |

[8] | Yun, SIAM Journal for Scientific Computing 25 pp 284– (2003) |

[9] | Yun, Communications in Numerical Methods in Engineering 20 pp 401– (2003) |

[10] | Johnston, International Journal for Numerical Methods in Engineering 45 pp 1333– (1999) |

[11] | Johnston, International Journal for Numerical Methods in Engineering 47 pp 1709– (2000) |

[12] | Johnston, International Journal for Numerical Methods in Engineering 48 pp 949– (2000) |

[13] | Johnston, Journal for Computational and Applied Mathematics 146 pp 231– (2002) |

[14] | Sladek, International Journal for Numerical Methods in Engineering 36 pp 1609– (1993) |

[15] | Sladek, Applied Mathematical Modelling 25 pp 901– (2001) |

[16] | Ma, Engineering Analysis with Boundary Elements 26 pp 329– (2002) |

[17] | Zhang, Computational Mechanics 23 pp 389– (1999) |

[18] | Dirgantara, International Journal of Fracture 105 pp 27– (2000) |

[19] | Aliabadi, International Journal for Numerical Methods in Engineering 48 pp 995– (2000) |

[20] | Liu, International Journal for Numerical Methods in Engineering 41 pp 541– (1998) |

[21] | Calculus (5th edn). Thomson: CA, USA, 2003. |

[22] | Calculus with Analytic Geometry (7th edn). Houghton Mifflin: Boston, 2002. |

[23] | Donaldson, SIAM Journal of Numerical Analysis 9 pp 573– (1972) |

[24] | Error analysis for a sinh transformation used in evaluating nearly singular boundary element integrals, submitted. · Zbl 1116.65032 |

[25] | Mathematical Handbook of Formulas and Tables. Schaum’s Outline Series (International edn) McGraw-Hill: New York, 1990. |

[26] | Table of Integrals, Series and Products (5th edn). Academic Press: Boston, 1994. |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.