A FREQUENCY-DEPENDENT FINITE-DIFFERENCE TIME-DOMAIN FORMULATION FOR GENERAL DISPERSIVE MEDIA

A weakness of the finite-difference time-domain (FDTD) method is that dispersion of the dielectric properties of the scattering/absorbing bo dy is often ignored and frequency-independent properties are generally taken. While this is not a disadvantage for CW or narrow-band irradia tion, the results thus obtained may be highly erroneous for short puls es where ultrawide bandwidths are involved. In some recent publication s, procedures based on a convolution integral describing D(t) in terms of E(t) are given for media for which the complex permittivity epsilo n(omega) may be described by a single-order Debye relaxation equation or a modified version thereof. Procedures are, however, needed for ge neral dispersive media for which epsilon(omega) and mu*(omega) may be expressible in terms of rational functions, or for human tissues wher e multiterm Debye relaxation equations must generally be used. We desc ribe a new differential equation approach, which can be used for gener al dispersive media. In this method D(t) is expressed in terms of E(t) by means of a differential equation involving D, E. and their time de rivatives. The method is illustrated by means of one- and three-dimens ional examples of media for which epsilon(omega) is given by a multit erm Debye equation, and for an approximate two-thirds muscle-equivalen t model of the human body.