Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations
arXiv:1312.1703 · doi:10.1109/TAP.2014.2367492
Abstract
We present a generic technique, automated by computer-algebra systems and available as open-source software \cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solvers has often required the aggregation of multiple disparate integral-evaluation schemes to treat all of the distinct types of integrals needed for a given BEM formulation; in contrast, our technique allows many different types of integrals to be handled by the \emph{same} algorithm and the same code implementation. Our method is a significant generalization of the Taylor--Duffy approach \cite{Taylor2003,Duffy1982}, which was originally presented for just a single type of integrand; in addition to generalizing this technique to a broad class of integrands, we also achieve a significant improvement in its efficiency by showing how the \emph{dimension} of the final numerical integral may often be reduced by one. In particular, if is the number of common vertices between the two triangles, in many cases we can reduce the dimension of the integral from to , obtaining a closed-form analytical result for (the common-triangle case).
Revised version containing a rearrangement of the text and updates to figures
Cited by in corpus (5)
- PCMSolver: an Open-Source Library for Solvation Modeling
- General-purpose kernel regularization of boundary integral equations via density interpolation
- On the Generalization of DIRECTFN for Singular Integrals Over Quadrilateral Patches
- The scattering of a scalar beam from isotropic and anisotropic two-dimensional randomly rough Dirichlet or Neumann surfaces: The full angular intensity distributions
- Isoparametric singularity extraction technique for 3D potential problems in BEM