An enhanced finite difference time domain method for two dimensional Maxwell's equations

An efficient finite-difference time-domain (FDTD) algorithm is built to solve the transverse electric 2D Maxwell's equations with inhomogeneous dielectric media where the electric fields are discontinuous across the dielectric interface. The new algorithm is derived based upon the integral version of the Maxwell's equations as well as the relationship between the electric fields across the interface. It is an improvement over the contour-path effective-permittivity algorithm by including some extra terms in the formulas. The scheme is validated in solving the scattering of a dielectric cylinder with exact solution from Mie theory and is then compared with the above contour-path method, the usual staircase and the volume-average method. The numerical results demonstrate that the new algorithm has achieved significant improvement in accuracy over the other methods. Furthermore, the algorithm has a simple structure and can be merged into any existing FDTD software package very easily.

[1]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[2]  P. Barber.  Absorption and scattering of light by small particles , 1984 .

[3]  John R. Cary,et al.  A stable FDTD algorithm for non-diagonal, anisotropic dielectrics , 2007, J. Comput. Phys..

[4]  Ahmad Mohammadi,et al.  Contour-path effective permittivities for the two-dimensional finite-difference time-domain method. , 2005, Optics express.

[5]  Hung-chun Chang,et al.  The Second-order Condition of FDTD Method at Sloped Dielectric Interfaces by Averaging the Permittivities along Surface Normal , 2012, IEEE Transactions on Antennas and Propagation.

[6]  R. Mittra,et al.  A conformal finite-difference time-domain technique for modeling cylindrical dielectric resonators , 1999 .

[7]  Raj Mittra,et al.  A study of the nonorthogonal FDTD method versus the conventional FDTD technique for computing resonant frequencies of cylindrical cavities , 1992 .

[8]  K. Yee,et al.  A subgridding method for the time-domain finite-difference method to solve Maxwell's equations , 1991 .

[9]  Steven G. Johnson,et al.  Accurate finite-difference time-domain simulation of anisotropic media by subpixel smoothing. , 2009, Optics letters.

[10]  Ahmad Mohammadi,et al.  Dispersive contour-path algorithm for the two-dimensional finite-difference time-domain method. , 2008, Optics express.

[11]  A. Cangellaris,et al.  Effective permittivities for second-order accurate FDTD equations at dielectric interfaces , 2001, IEEE Microwave and Wireless Components Letters.

[12]  Moysey Brio,et al.  Overlapping Yee FDTD Method on Nonorthogonal Grids , 2009, J. Sci. Comput..

[13]  John R. Cary,et al.  A more accurate, stable, FDTD algorithm for electromagnetics in anisotropic dielectrics , 2012, J. Comput. Phys..

[14]  M. Brio,et al.  Subpixel smoothing finite-difference time-domain method for material interface between dielectric and dispersive media. , 2012, Optics letters.

[15]  Yuzo Yoshikuni,et al.  The second-order condition for the dielectric interface orthogonal to the Yee-lattice axis in the FDTD scheme , 2000 .

[16]  Shan Zhao,et al.  A second order dispersive FDTD algorithm for transverse electric Maxwell's equations with complex interfaces , 2016, Comput. Math. Appl..

[17]  Steven G. Johnson,et al.  Improving accuracy by subpixel smoothing in the finite-difference time domain. , 2006, Optics letters.

[18]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[19]  Y. Yoshikuni,et al.  Effective permittivities with exact second-order accuracy at inclined dielectric interface for the two-dimensional finite-difference time-domain method. , 2010, Applied optics.

[20]  Tatsuo Itoh,et al.  FDTD analysis of dielectric resonators with curved surfaces , 1997 .

[21]  John R. Cary,et al.  A second-order 3D electromagnetics algorithm for curved interfaces between anisotropic dielectrics on a Yee mesh , 2011, J. Comput. Phys..