The New Jersey Institute of Technology's
Electronic Theses & Dissertations Project

Title:	Asymptotic and numerical analysis of time-dependent wave propagation in dispersive dielectric media that exhibit fractional relaxation
Author:	Causley, Matthew Frank
View Online:	njit-etd2011-048 (xvi, 104 pages ~ 2.7 MB pdf)
Department:	Department of Mathematical Sciences
Degree:	Doctor of Philosophy
Program:	Mathematical Sciences
Document Type:	Dissertation
Advisory Committee:	Petropoulos, Peter G. (Committee chair) Kriegsmann, Gregory A. (Committee member) Luke, Jonathan H.C. (Committee member) Jiang, Shidong (Committee member) Whitman, Gerald Martin (Committee member)
Date:	2011-05
Keywords:	Time-dependent Maxwell's equations FDTD method Cole-Cole model Dispersive dielectric Havrilink-Negami model Fractional derivative
Availability:	Unrestricted
Abstract:	This dissertation addresses electromagnetic pulse propagation through anomalously dispersive dielectric media. The Havriliak-Negami (H-N) and Cole-Cole (C-C) models capture the non-exponential nature of such dielectric relaxation phenomena, which is manifest in a variety of dispersive dielectric media. In the C-C model, the dielectric polarization is coupled to the time-dependent Maxwell's equations by a fractional differential equation involving the electric field. In the H-N case, a more general pseudo-fractional differential operator describes the polarization. The development and analysis of a robust method for implementing such models is presented, with an emphasis on algorithmic efficiency. Separate numerical schemes are presented for C-C and H-N media. A straightforward numerical implementation of these models using finite-difference time-domain (FD-TD) techniques is expected to be second order accurate in both space and time. However due to the singular nature of the kernels appearing in the fractional convolution operators, the standard C-C implementation, produces first order accuracy in time. As we show, this method is equivalent to most approaches presented in the current literature, which implies that they are also first order. The desired accuracy is instead achieved by applying multistep methods to the fractional differential equation; however multistep methods are unnecessary in the H-N implementation to preserve the accuracy. Furthermore, the C-C model is a specific case of the H-N model and can therefore be constructed using the latter of these approaches. The FD-TD methods are validated by evaluating the electric field for the signaling problem, using numerical quadrature to evaluate the integral form of the solution. This is accomplished using the Green's function of the dispersive medium; in addition, the behavior of pulse propagation is studied asymptotically using the Green's function, which further validates the observed results of the numerical experiments.