Spherical polar coordinate transformation for integration of singular functions on tetrahedra

arXiv2110.07270MaRDI QIDQ6380315FDOQ6380315

Publication date: 14 October 2021

Abstract: A method is presented for the evaluation of integrals on tetrahedra where the integrand has an integrable singularity at one vertex. The approach uses a transformation to spherical polar coordinates which explicitly eliminates the singularity and facilitates the evaluation of integration limits. The method is also implemented in an adaptive form which gives convergence to a required tolerance. Results from the method are compared to the output from an exact analytical method for one tetrahedron and show high accuracy. In particular, when the adaptive algorithm is used, highly accurate results are found for poorly conditioned tetrahedra which normally present difficulties for numerical quadrature techniques. The approach is also demonstrated for evaluation of the Biot-Savart integral on an unstructured mesh in combination with a fixed node quadrature rule and demonstrates good convergence and accuracy.

Has companion code repository: https://github.com/mjcarley/tetquad

This page was built for publication: Spherical polar coordinate transformation for integration of singular functions on tetrahedra

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6380315)