Finite-difference computation of transient electromagnetic waves for cylindrical geometries in complex media

We present two novel, fully; three-dimensional (3-D) finite-difference time -domain (FDTD) schemes in cylindrical coordinates for transient simulation of electromagnetic wave propagation in complex (inhomogeneous, dispersive, and conductive) and unbounded media. The proposed FDTD schemes incorporate an extension of the perfectly matched layer (PML) absorbing boundary condit ion (ABC) to three-dimensional (3-D) cylindrical coordinates. Dispersion on the media is modeled by using the piecewise-linear recursive Convolution ( PLRC) algorithm, accounting for multiterm Lorentz and/or Debye models. Spli t-field and unsplit (anisotropic medium) formulations of the cylindrical PM L-PLRC-FDTD schemes are implemented and compared in the time domain. The co mparison includes the late-time stability properties of the update schemes. Numerical simulations of susbsurface electromagnetic problems are included . Because the proposed schemes retain the nearest-neighbor property of the ordinary FDTD, they are well suited for implementation on massively paralle l computers.