New method for the efficient summation of double infinite series arising from the spectral domain analysis of frequency selective surfaces

When the method of moments (MoM) in the spectral domain is applied to the analysis of frequency selective surfaces, the entries of the MoM matrix are slowly convergent double infinite series. In this paper, a two-step acceleration technique is developed which makes it possible the fast and accurate computation of these double series in the particular case where subsectional rooftops are used as basis functions. The technique is based on a combination of the use of Kummer's transformation, the use of Poisson's transformation, and the determination of judicious Chebyshev polynomial interpolations of some of the spectral discrete functions involved in the infinite series. The results obtained show that when all the double series of the MoM matrix are to be computed with an accuracy of three significant figures, the new acceleration technique turns out to be about one thousand times faster than brute-force computation, and a few times faster than the acceleration technique based on fast Fourier transform.