Numerical computation of 2D Sommerfeld integrals -- a novel asymptotic extraction technique (Q1184616)

It is important to develop efficient techniques for the calculation of the elements in the impedance matrix. In a previous paper [ibid. 98, No. 2, 199-216 (1992; Zbl 0744.65018)] the same authors have shown that the angular integral can be decomposed into a finite number of incomplete Lipschitz-Hankel integrals (ILHI). The decomposition of the angular integral in terms of ILHIs can also be used as the starting point for the derivation of a novel asymptotic extraction technique (AET) for the outer infinite integral. The general idea is to use the series expansions for the ILHIs to find an asymptotic expansion for the integrand which holds for large values of the spectral variable \(\lambda\), and which when integrated from some lower limit \(L\) to infinity can be efficiently calculated. If this can be accomplished, then numerical integration will only be required over the range from 0 to \(L\), thereby significantly improving the efficiency of the algorithm. A novel AET is derived using the above procedure.