A frequency-dependent finite-difference time-domain formulation for transient propagation in plasma

Previous FDTD (finite-difference time-domain) formulations were not capable of analyzing plasmas for two reasons. First, FDTD requires that at each time step the permittivity and conductivity be specified as constants that do not depend on frequency, while even for the simplest plasmas these parameters vary with frequency. Second, the permittivity of a plasma can be negative, which can cause terms in FDTD expressions to become singular. A novel FDTD formulation for frequency-dependent materials (FD)/sup 2/TD has been developed. It is shown that (FD)/sup 2/TD can be applied to compute transient propagation in plasma when the plasma can be characterized by a complex frequency-dependent permittivity. While the computational example presented is for a pulse normally incident on an isotropic plasma slab, the (FD)/sup 2/TD formulation is fully three-dimensional. It can accommodate arbitrary transient excitation, with the limitation that the excitation pulse must have no zero frequency energy component. Time-varying electron densities and/or collision frequencies could also be included. The formulation presented is for an isotropic plasma, but extension to anisotropic plasma should be fairly straightforward. >

[1]  A. Bayliss,et al.  Radiation boundary conditions for wave-like equations , 1980 .

[2]  J. Lee Pulse distortion through a lossy plasma medium , 1979 .

[3]  R. M. Hyde An experimental study of the absorption of energy from a large amplitude electromagnetic pulse by a collisionless plasma , 1979 .

[4]  K. Tschu,et al.  Propagation of electromagnetic pulse with a Gaussian envelope through an inhomogeneous anisotropic plasma layer , 1987 .

[5]  K. Kunz,et al.  Finite-Difference Analysis of EMP Coupling to Lossy Dielectric Structures , 1980, IEEE Transactions on Electromagnetic Compatibility.

[6]  G. Mur Absorbing Boundary Conditions for the Finite-Difference Approximation of the Time-Domain Electromagnetic-Field Equations , 1981, IEEE Transactions on Electromagnetic Compatibility.

[7]  Norinobu Yoshida,et al.  Transient analysis of a magnetized plasma in three-dimensional space , 1988 .

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

[9]  Karl Kunz,et al.  A Three-Dimensional Finite-Difference Solution of the External Response of an Aircraft to a Complex Transient EM Environment: Part I-The Method and Its Implementation , 1978, IEEE Transactions on Electromagnetic Compatibility.

[10]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[11]  Allen Taflove,et al.  A numerical technique for analyzing electromagnetic wave scattering from moving surfaces in one and two dimensions , 1989 .

[12]  Andreas C. Cangellaris,et al.  Conformal time domain finite difference method , 1984 .

[13]  A. Papoulis,et al.  The Fourier Integral and Its Applications , 1963 .

[14]  A. Taflove,et al.  Numerical Solution of Steady-State Electromagnetic Scattering Problems Using the Time-Dependent Maxwell's Equations , 1975 .

[15]  J. Wait Electromagnetic wave theory , 1985 .

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

[17]  Richard Holland,et al.  Finite-Difference Analysis of EMP Coupling to Thin Struts and Wires , 1981, IEEE Transactions on Electromagnetic Compatibility.

[18]  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.