Plane waves in FDTD simulations and a nearly perfect total-field/scattered-field boundary

The total-field/scattered-field (TFSF) boundary has been successfully used for a number of years to introduce energy into finite-difference time-domain (FDTD) grids. If the propagation of the incident field is grid-aligned, a perfect TFSF implementation can be realized by using an auxiliary one-dimensional FDTD simulation which models propagation of the incident field. Here "perfect" implies the incident field propagation exactly matches the way in which the field propagates in the FDTD grid. However, for propagation which is not grid-aligned, no similarly perfect implementation has previously been presented. This work provides a framework for a perfect TFSF boundary for pulsed plane waves which do not propagate in a grid-aligned fashion. To achieve this, homogeneous plane-wave propagation is rigorously quantified. Using this knowledge and a specification of the desired incident field, the dispersion relation is used to ascertain the incident field at any point in the grid. It is required to account for, unlike in the continuous world, the electric field, the magnetic field, and the wavenumber vector not forming a mutually orthogonal set. Group velocity is also considered because of its relevance to the implementation.

[1]  Allen Taflove,et al.  Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures , 1988 .

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

[3]  C. Guiffaut,et al.  A perfect wideband plane wave injector for FDTD method , 2000, IEEE Antennas and Propagation Society International Symposium. Transmitting Waves of Progress to the Next Millennium. 2000 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (C.

[4]  L. Trefethen.  Group velocity in finite difference schemes , 1981 .

[5]  Allen Taflove,et al.  Efficient modeling of infinite scatterers using a generalized total-field/scattered-field FDTD boundary partially embedded within PML , 2002 .

[6]  John B. Schneider,et al.  FDTD dispersion revisited: faster-than-light propagation , 1999 .

[7]  Wojciech Gwarek,et al.  On the nature of solutions produced by finite difference schemes in time domain , 1999 .

[8]  K. Demarest,et al.  FDTD modeling of scatterers in stratified media , 1995 .

[9]  R. J. Kruhlak,et al.  Plane waves and planar boundaries in FDTD simulations , 2000 .

[10]  Steven G. Johnson,et al.  FFTW: an adaptive software architecture for the FFT , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[11]  D. Merewether,et al.  On Implementing a Numeric Huygen's Source Scheme in a Finite Difference Program to Illuminate Scattering Bodies , 1980, IEEE Transactions on Nuclear Science.

[12]  Levent Gurel,et al.  Interpolation techniques to improve the accuracy of the plane wave excitations in the finite difference time domain method , 1997 .

[13]  Orhan Arikan,et al.  An efficient and accurate technique for the incident-wave excitations in the FDTD method , 1998 .

[14]  J. E. Roy,et al.  The dispersion characteristics of the FD-TD method , 1989, Digest on Antennas and Propagation Society International Symposium.

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

[16]  김덕영.  [신간안내] Computational Electrodynamics (the finite difference time - domain method) , 2001 .

[17]  R. J. Kruhlak,et al.  Dispersion of homogeneous and inhomogeneous waves in the Yee finite-difference time-domain grid , 2001 .