A second-order 3D electromagnetics algorithm for curved interfaces between anisotropic dielectrics on a Yee mesh

A new frequency-domain electromagnetics algorithm is developed for simulating curved interfaces between anisotropic dielectrics embedded in a Yee mesh with second-order error in resonant frequencies. The algorithm is systematically derived using the finite integration formulation of Maxwell's equations on the Yee mesh. Second-order convergence of the error in resonant frequencies is achieved by guaranteeing first-order error on dielectric boundaries and second-order error in bulk (possibly anisotropic) regions. Convergence studies, conducted for an analytically solvable problem and for a photonic crystal of ellipsoids with anisotropic dielectric constant, both show second-order convergence of frequency error; the convergence is sufficiently smooth that Richardson extrapolation yields roughly third-order convergence. The convergence of electric fields near the dielectric interface for the analytic problem is also presented.

[1]  M. Clemens,et al.  Discrete Electromagnetism With the Finite Integration Technique - Abstract , 2001 .

[2]  J. Jackson.  Classical Electrodynamics, 3rd Edition , 1998 .

[3]  J. Hesthaven,et al.  Staircase-free finite-difference time-domain formulation for general materials in complex geometries , 2001 .

[4]  Bo Zhang,et al.  An explicit fourth-order staggered finite-difference time-domain method for Maxwell's equations , 2002 .

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

[6]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[7]  N. Myung,et al.  Locally tensor conformal FDTD method for modeling arbitrary dielectric surfaces , 1999 .

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

[9]  Steven G. Johnson,et al.  Perturbation theory for anisotropic dielectric interfaces, and application to subpixel smoothing of discretized numerical methods. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Steven G. Johnson,et al.  Accurate FDTD simulation of anisotropic media by subpixel smoothing , 2009 .

[11]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[12]  R. Mittra,et al.  A conformal finite difference time domain technique for modeling curved dielectric surfaces , 2001, IEEE Microwave and Wireless Components Letters.

[13]  Xianliang Wu,et al.  Application of the symplectic finite-difference time-domain scheme to electromagnetic simulation , 2007, J. Comput. Phys..

[14]  J. Hesthaven,et al.  Convergent Cartesian Grid Methods for Maxwell's Equations in Complex Geometries , 2001 .

[15]  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.

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

[17]  Allen Taflove,et al.  Finite-difference time-domain modeling of curved surfaces (EM scattering) , 1992 .

[18]  Peter Deuflhard,et al.  A 3-D tensor FDTD-formulation for treatment of sloped interfaces in electrically inhomogeneous media , 2003 .

[19]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[20]  R. Holland,et al.  Pitfalls of staircase meshing , 1993 .

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

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

[23]  T. Weiland.  A discretization model for the solution of Maxwell's equations for six-component fields , 1977 .

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

[25]  Wei Cai,et al.  An upwinding embedded boundary method for Maxwell's equations in media with material interfaces: 2D case , 2003 .

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

[27]  R. Mittra,et al.  A locally conformal finite-difference time-domain (FDTD) algorithm for modeling three-dimensional perfectly conducting objects , 1997 .

[28]  Amir Yefet,et al.  A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations , 2001 .

[29]  John R. Cary,et al.  Application of Dey-Mittra conformal boundary algorithm to 3D electromagnetic modeling , 2009, J. Comput. Phys..

[30]  Shan Zhao,et al.  High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces , 2004 .

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

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