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.
Key words.Maxwell's equations, finite difference schemes, computational electromagnetics, domain decomposition methods, interpolation, dispersive media, metallic nanoparticles AMS subject classifications. 65M06, 78M20 DOI. 10.1137/070705374 1. Introduction. One of the main difficulties in the development of reliable tools for computational electromagnetics (CEM) based on conventional algorithms like Yee's method [20,21] (commonly known as the finite difference time domain method (FDTD)) is the treatment of curved interfaces. Brute force discretization using only a rectangular grid requires excessively small grid sizes to reasonably resolve the curvature of the interface, and this frequently forces the use of large amounts of memory and unreasonably long computation times. Furthermore, numerical artifacts (known as staircasing errors) can be introduced and corrupt the solutions. Other algorithms for CEM are discussed in [3,4,7,12,15,20].Recently, domain decomposition methods (DDMs) [18] based on overlapping (or composite) grids have been successfully applied to CEM problems with complex geometries [6,8]. Henshaw [8] solved Maxwell's equations as a second-order vector wave equation rather than as a first-order system and developed high-order schemes by discretizing the Laplace operator. The jump (or transmission) condition at the interface of two different constant media are derived using boundary-(or interface-) fitted grids.In this paper, we suggest a more general method. We solve the two-dimensional Maxwell equations as a first-order system with constitutive equations. By doing this, we separate numerical stability issues from material modeling issues (i.e., stable schemes for Maxwell's equations can be used with various materials by simply