Fast Simulation of Scattering Problem for Magnetodielectric Materials With General Anisotropy in Layered Media

In this paper, the mixed-order stabilized biconjugate gradient fast Fourier transform (mixed-order BCGS-FFT) method is presented to solve the scattering problem of magnetodielectric materials with general anisotropy in layered media. While the volumetric roof-top functions are used as the testing functions for the coupled field volume integral equation and the basis functions for flux densities, the second-order curl conforming basis functions are applied to expand the vector potentials with the aim of both preserving the continuity of their tangential components and avoiding the zero terms that might otherwise be caused by the divergence operator. The layered medium Green's function (LMGF) is efficiently evaluated through the recursive matrix method along with an interpolation technique. Several numerical experiments are presented to demonstrate the high accuracy and efficiency of the method. Different from the previously published works that aim to solve the similar problem, the new contribution of this work is to extend the mixed-order BCGS-FFT method to accommodate the layered background medium. Therefore, the 3-D FFT acceleration for integral kernels associated with the LMGF as well as the interpolation technique has been implemented and combined with the mixed-order BCGS-FFT method.