A numerical technique is reported for the evaluation of improper integrals associated with the self-interaction terms arising in the application of linear (Rao-Wilton-Glisson) current basis functions, defined on planar triangular patches, to three-dimensional electromagnetic surface scattering problems. The two-dimensional numerical integration arising in more conventional approaches, which follow Graglia [1993], is replaced by a one-dimensional integration by means of a suitable change of the local coordinate system, and analytical expressions for the functions to be numerically integrated are derived. Numerical results obtained using Graglia's method, our alternative method, and a reliable reference solution are compared for accuracy and computational complexity. The alternative technique appears to be conceptually simpler than the conventional method, is easier to implement, and causes no degradation in accuracy; in fact, it seems to more efficiently achieve a slightly specified level of accuracy.
[1]
D. Wilton,et al.
Electromagnetic scattering by surfaces of arbitrary shape
,
1980
.
[2]
D. J. Shippy,et al.
Alternative integration formulae for triangular finite elements
,
1981
.
[3]
Roberto D. Graglia,et al.
On the numerical integration of the linear shape functions times the 3-D Green's function or its gradient on a plane triangle
,
1993
.
[4]
I. S. Gradshteyn,et al.
Table of Integrals, Series, and Products
,
1976
.
[5]
Roland Klees,et al.
Numerical calculation of weakly singular surface integrals
,
1996
.
[6]
D. Wilton,et al.
Potential integrals for uniform and linear source distributions on polygonal and polyhedral domains
,
1984
.