A stable fourth-order FDTD method for modeling electrically long dielectric waveguides

Accurate time-domain modeling of electrically long structures such as dielectric waveguides in integrated optics often demands prohibitive computational cost. In this paper, an efficient fourth-order finite-difference time-domain (FDTD) method is investigated in a two-dimensional (2-D) transverse electric (TE) case. The method is a combination of fourth-order staggered backward-differentiation time integrator and fourth-order staggered spatial discretization. A rigorous stability and dispersion analysis is conducted to show its useful numerical characteristics. Fourth-order convergence of the numerical scheme is demonstrated by monitoring numerical errors in the L/sub 2/ norm in 2-D cavities. Numerical efficiency of the fourth-order method is validated through its applications for full-wave time-domain simulations of long 2-D optical waveguide structures.

[1]  Kyu-Pyung Hwang.  A fourth-order accurate FDTD scheme with long-time stability , 2005 .

[2]  S. Obayya.  Novel finite element analysis of optical waveguide discontinuity problems , 2004, Journal of Lightwave Technology.

[3]  C. E. Rubio-Mercedes,et al.  Pade/spl acute/ boundary conditions for the finite-element modeling of arbitrary planar junctions , 2004, Journal of Lightwave Technology.

[4]  C. Balanis,et al.  A hybrid fourth-order FDTD utilizing a second-order FDTD subgrid , 2001, IEEE Microwave and Wireless Components Letters.

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

[6]  Angel Vegas,et al.  Analyzing the stability of the FDTD technique by combining the von Neumann method with the Routh-Hurwitz criterion , 2001 .

[7]  David W. Zingg,et al.  Numerical Solution of the Time-Domain Maxwell Equations Using High-Accuracy Finite-Difference Methods , 2000, SIAM J. Sci. Comput..

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

[9]  Nikolaos V. Kantartzis,et al.  A generalized methodology based on higher‐order conventional and non‐standard FDTD concepts for the systematic development of enhanced dispersionless wide‐angle absorbing perfectly matched layers , 2000 .

[10]  Tobin A. Driscoll,et al.  Staggered Time Integrators for Wave Equations , 2000, SIAM J. Numer. Anal..

[11]  J. Shang,et al.  High-Order Compact-Difference Schemes for Time-Dependent Maxwell Equations , 1998 .

[12]  Jeffrey L. Young,et al.  Toward the construction of a fourth-order difference scheme for transient EM wave simulation: staggered grid approach , 1997 .

[13]  Melinda Piket-May,et al.  A modified FDTD (2, 4) scheme for modeling electrically large structures with high-phase accuracy , 1997 .

[14]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

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

[16]  Peter G. Petropoulos,et al.  Phase error control for FD-TD methods of second and fourth order accuracy , 1994 .

[17]  Robert J. Lee,et al.  On the Accuracy of Numerical Wave Simulations Based on Finite Methods , 1992 .

[18]  Sai T. Chu,et al.  A finite-difference time-domain method for the design and analysis of guided-wave optical structures , 1989 .

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

[20]  Benjamin C. Kuo,et al.  AUTOMATIC CONTROL SYSTEMS , 1962, Universum:Technical sciences.