A stable FDTD algorithm for non-diagonal, anisotropic dielectrics

A stable FDTD algorithm is developed for simulating Maxwell’s equations in anisotropic dielectric materials with principal axes not aligned with the grid. The algorithm is stable because the finite-difference operator that converts D to E is symmetric and positive semidefinite; for contrast, a previously developed asymmetric algorithm is shown to suffer from late-time instabilities. The presented algorithm has second-order error for continuous dielectric materials, and the error can be reduced to third-order by Richardson extrapolation. Applied to dielectrics with sharp interfaces, the algorithm has first-order error, even when averaging the dielectric within partially filled grid cells. However, averaging the dielectric permits Richardson extrapolation to obtain second-order error.

[1]  Miller Puckette,et al.  Accuracy of frequency estimates using the phase vocoder , 1998, IEEE Trans. Speech Audio Process..

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

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

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

[5]  T. I. Kosmanis,et al.  A systematic and topologically stable conformal finite-difference time-domain algorithm for modeling curved dielectric interfaces in three dimensions , 2003 .

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

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

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

[9]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[10]  Rolf Schuhmann,et al.  A uniformly stable conformal FDTD‐method in Cartesian grids , 2003 .

[11]  Salvador García,et al.  On the application of finite methods in time domain to anisotropic dielectric waveguides , 1996 .

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

[13]  J. Cary,et al.  VORPAL: a versatile plasma simulation code , 2004 .

[14]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[15]  Steven G. Johnson,et al.  Block-iterative frequency-domain methods for Maxwell's equations in a planewave basis. , 2001, Optics express.

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

[17]  A. Vegas,et al.  A treatment of magnetized ferrites using the FDTD method , 1993, IEEE Microwave and Guided Wave Letters.

[18]  Steven G. Johnson,et al.  Photonic Crystals: Molding the Flow of Light , 1995 .

[19]  Zhizhang Chen,et al.  The generalized TLM-based FDTD modeling of frequency-dependent and anisotropic media , 1998 .

[20]  K. Kunz,et al.  Finite-difference time-domain analysis of gyrotropic media. I. Magnetized plasma , 1992 .

[21]  Antti V. Räisänen,et al.  A generalized compact 2‐D FDTD model for the analysis of guided modes of anisotropic waveguides with arbitrary tensor permittivity , 1998 .

[22]  A. Taflove.  The Finite-Difference Time-Domain Method , 1995 .

[23]  V. Mandelshtam,et al.  Harmonic inversion of time signals and its applications , 1997 .