SIAM JOURNAL ON SCIENTIFIC COMPUTING, v.32, no.2, pp.684 - 702
Abstract
The domain decomposition method based on overlapping grids is developed to solve the two-dimensional Maxwell equations in the time domain. The finite difference schemes for rectangular and polar coordinate systems are presented. Since interpolation plays a crucial role in our method, the Newton and the Fourier interpolation methods are surveyed in detail. The computational studies of the electromagnetic wave propagation in free space and the back-scattering by a perfect electric conducting object of a circular shape are performed to test the accuracy, the convergence, and the efficiency of our method. Moreover, we give a methodology to model dispersive media in time domain simulations by introducing Drude conductivity in the constitutive equations. The problem of light scattering by metallic nanoparticles is solved, and its results show that our algorithm is efficient and reliable in capturing the small scale phenomena.