Evaluation of Biot-Savart integrals on tetrahedral meshes
For computational fluid dynamics researchers, this offers a faster alternative to analytical integration for vortex methods, though it is an incremental improvement.
The paper presents a method for evaluating Biot-Savart integrals on tetrahedral meshes using Gaussian quadrature and ray tracing, eliminating costly transcendental functions. The method achieves second-order convergence and near-linear parallel speedup.
An arithmetically simple method has been developed for the evaluation of Biot--Savart integrals on tetrahedralized distributions of vorticity. In place of the usual approach of analytical formulae for the velocity induced by a linear distribution of vorticity on a tetrahedron, the integration is performed using Gaussian quadrature and a ray tracing technique from computer graphics. This eliminates completely the need for the evaluation of square roots, logarithms and arc tangents, and almost completely eliminates the requirement for trigonometric functions, with no operation more costly than a division required during the main calculation loop. An assessment of the algorithm's performance is presented, demonstrating its accuracy, second order convergence and near-linear speedup on parallel systems.