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 condition (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. Split-field and unsplit (anisotropic medium) formulations of the cylindrical PML-PLRC-FDTD schemes are implemented and compared in the time domain. The comparison includes the late-time stability properties of the update schemes. Numerical simulations of subsurface 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 parallel computers.

[1]  J. Calhoun.  A finite difference time domain (FDTD) simulation of electromagnetic wave propagation and scattering in a partially conducting layered Earth , 1997, IGARSS'97. 1997 IEEE International Geoscience and Remote Sensing Symposium Proceedings. Remote Sensing - A Scientific Vision for Sustainable Development.

[2]  R. B. Standler,et al.  A frequency-dependent finite-difference time-domain formulation for dispersive materials , 1990 .

[3]  Peter Monk,et al.  The Perfectly Matched Layer in Curvilinear Coordinates , 1998, SIAM J. Sci. Comput..

[4]  A Taflove,et al.  Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses. , 1991, Optics letters.

[5]  Weng Cho Chew,et al.  Finite-difference time-domain simulation of ground penetrating radar on dispersive, inhomogeneous, and conductive soils , 1998, IEEE Trans. Geosci. Remote. Sens..

[6]  Allen Taflove,et al.  Theory and application of radiation boundary operators , 1988 .

[7]  Dennis M. Sullivan,et al.  Frequency-dependent FDTD methods using Z transforms , 1992 .

[8]  Carey M. Rappaport,et al.  A general method for FDTD modeling of wave propagation in arbitrary frequency-dispersive media , 1997 .

[9]  Weng Cho Chew,et al.  Modeling of the dielectric logging tool at high frequencies: theory , 1988 .

[10]  Om P. Gandhi,et al.  A frequency-dependent finite-difference time-domain formulation for general dispersive media , 1993 .

[11]  Thomas Weiland,et al.  A CONSISTENT SUBGRIDDING SCHEME FOR THE FINITE DIFFERENCE TIME DOMAIN METHOD , 1996 .

[12]  Fernando L. Teixeira,et al.  General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media , 1998 .

[13]  Weng Cho Chew,et al.  Modeling Of The Subsurface Interface Radar , 1990, 10th Annual International Symposium on Geoscience and Remote Sensing.

[14]  D. Sullivan,et al.  Three dimensional optical fiber simulation , 1999, IEEE Antennas and Propagation Society International Symposium. 1999 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.99CH37010).

[15]  E. Michielssen,et al.  Complex coordinate system as a generalized absorbing boundary condition , 1997, IEEE Antennas and Propagation Society International Symposium 1997. Digest.

[16]  C. J. Railton,et al.  Passive equivalent circuit of FDTD: an application to subgridding , 1997 .

[17]  Stephen D. Gedney,et al.  An Anisotropic PML Absorbing Media for the FDTD Simulation of Fields in Lossy and Dispersive Media , 1996 .

[18]  3D electromagnetic modeling using staggered finite differences , 1997, IGARSS'97. 1997 IEEE International Geoscience and Remote Sensing Symposium Proceedings. Remote Sensing - A Scientific Vision for Sustainable Development.

[19]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, IMS 1996.

[20]  P. Monk,et al.  Optimizing the Perfectly Matched Layer , 1998 .

[21]  W. Chew,et al.  Differential Forms, Metrics, and the Reflectionless Absorption of Electromagnetic Waves , 1999 .

[22]  Weng Cho Chew,et al.  Analytical derivation of a conformal perfectly matched absorber for electromagnetic waves , 1998 .

[23]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, 1996 IEEE MTT-S International Microwave Symposium Digest.

[24]  Glenn S. Smith,et al.  A fully three-dimensional simulation of a ground-penetrating radar: FDTD theory compared with experiment , 1996, IEEE Trans. Geosci. Remote. Sens..

[25]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[26]  Carey M. Rappaport.  Interpreting and improving the PML absorbing boundary condition using anisotropic lossy mapping of space , 1996 .

[27]  Qing Huo Liu,et al.  Inversion of induction tool measurements using the distorted Born iterative method and CG-FFHT , 1994, IEEE Trans. Geosci. Remote. Sens..

[28]  M. Stuchly,et al.  Hybrid PEE-FDTD method for efficient field modelling in cylindrical co-ordinates , 1996 .

[29]  Jian-Ming Jin,et al.  Perfectly Matched Layers in the Discretized Space: An Analysis and Optimization , 1996 .

[30]  Liang C. Shen,et al.  Numerical simulation of subsurface radar for detecting buried pipes , 1991, IEEE Trans. Geosci. Remote. Sens..

[31]  James R. Wait,et al.  Transient Electromagnetic Fields , 1976 .

[32]  W. Chew,et al.  Lattice electromagnetic theory from a topological viewpoint , 1999 .

[33]  Q.H. Liu,et al.  A PML-FDTD algorithm for general dispersive media in GPR and plasma applications , 1998, IEEE Antennas and Propagation Society International Symposium. 1998 Digest. Antennas: Gateways to the Global Network. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.98CH36.

[34]  J. E. Hipp.  Soil electromagnetic parameters as functions of frequency, soil density, and soil moisture , 1974 .

[35]  Qing Huo Liu,et al.  Computation of transient electromagnetic waves in inhomogeneous media , 1991 .

[36]  W. Chew,et al.  Response of a Point Source in a Multicylindrcally Layered Medium , 1987, IEEE Transactions on Geoscience and Remote Sensing.

[37]  W. Chew,et al.  Systematic derivation of anisotropic PML absorbing media in cylindrical and spherical coordinates , 1997 .

[38]  M. Fusco,et al.  FDTD algorithm in curvilinear coordinates (EM scattering) , 1990 .

[39]  Hans Wenzel,et al.  The effective frequency method in the analysis of vertical-cavity surface-emitting lasers , 1997 .

[40]  Qing Huo Liu,et al.  A nonuniform cylindrical FDTD algorithm with improved PML and quasi-PML absorbing boundary conditions , 1999, IEEE Trans. Geosci. Remote. Sens..

[41]  Weng Cho Chew,et al.  Modeling of the dielectric logging tool at high frequencies: applications and results , 1988 .

[42]  W. Chew.  Waves and Fields in Inhomogeneous Media , 1990 .

[43]  Weng Cho Chew,et al.  Application of perfectly matched layers to the transient modeling of subsurface EM problems , 1997 .

[44]  E. Miller,et al.  Pole extraction from real-frequency information , 1980, Proceedings of the IEEE.

[45]  Dennis M. Sullivan,et al.  A frequency-dependent FDTD method for biological applications , 1992 .

[46]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[47]  R. J. Luebbers,et al.  Piecewise linear recursive convolution for dispersive media using FDTD , 1996 .

[48]  Weng Cho Chew,et al.  Response ofaPoint Source inaMulticylindrcall y Layered Medium , 1987 .

[49]  Weng Cho Chew,et al.  Electromagnetic borehole fields in a layered, dipping bed environment with invasion , 1991 .

[50]  Weng Cho Chew,et al.  PML-FDTD in cylindrical and spherical grids , 1997 .

[51]  Dennis M. Sullivan,et al.  Comparison of measured and simulated data in an annular phased array using an inhomogeneous phantom , 1992 .

[52]  D. Katz,et al.  Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FD-TD meshes , 1994, IEEE Microwave and Guided Wave Letters.

[53]  Jin-Fa Lee,et al.  A perfectly matched anisotropic absorber for use as an absorbing boundary condition , 1995 .

[54]  David Gottlieb,et al.  A Mathematical Analysis of the PML Method , 1997 .

[55]  Ian J Craddock,et al.  Derivation and application of a passive equivalent circuit for the finite difference time domain algorithm , 1996 .

[56]  Ichiro Fukai,et al.  A finite-difference time-domain formulation for transient propagation in dispersive media associated with Cole-Cole's circular arc law , 1990 .

[57]  Qing Huo Liu,et al.  Electromagnetic field generated by an off-axis source in a cylindrically layered medium with an arbitrary number of horizontal discontinuities , 1993 .

[58]  Yiwei He,et al.  Polarization effects on two-dimensional active imaging of conducting objects buried in a dielectric half-space , 1993, Proceedings of IEEE Antennas and Propagation Society International Symposium.

[59]  Weng Cho Chew,et al.  Transient response of some borehole mandrel tools , 1989 .

[60]  Radiation patterns of higher azimuthal order spatial modes from a concentric-circle-grating waveguide cavity using the volume-current method , 1998 .