Search

The finite-difference time-domain (FDTD) method has been widely used for the electromagnetic analysis of complex dispersive media. The shift-operator (SO)-FDTD or auxiliary differential equation (ADE)-FDTD formulation has been mainly employed for dusty plasma. Each FDTD formulation has its pros and cons. SO-FDTD is accurate but not computationally efficient. ADE-FDTD needs fewer computational resources, but its accuracy is poor. Here, we propose an accurate and efficient FDTD formulation for dusty plasma, based on the bilinear transform (BT). We perform a comprehensive study on the numerical permittivity and the computational efficiency for three FDTD formulations. Numerical examples are employed to illustrate that the proposed BT-FDTD outperforms the previously reported FDTD formulations for dusty plasma. In addition, based on the proposed BT-FDTD simulations, the effect of dust particles on EM wave propagation is investigated in the GHz band and the THz band. [ABSTRACT FROM AUTHOR]

The complex-conjugate pole-residue (CCPR) model has been popularly adopted because CCPR-finite-difference time domain (FDTD) can reduce the memory requirement with the help of complex conjugate property of auxiliary variables. To fully utilize CCPR-FDTD, it is of great necessity to investigate its numerical stability since the FDTD method is conditionally stable. Nonetheless, the numerical stability conditions of CCPR-FDTD have not been studied because its derivation is not straightforward. In this communication, the numerical stability conditions of CCPR-FDTD are systematically derived by combining the von Neumann method with Routh–Hurwitz criterion. It is found that the numerical stability conditions of CCPR-FDTD are the same as those of the modified Lorentz-FDTD with bilinear transform. Moreover, the numerical accuracy of CCPR-FDTD is studied, and numerical examples are employed to validate this work. [ABSTRACT FROM AUTHOR]

Finite-difference time domain (FDTD) has been widely used to analyze electromagnetic wave interaction with dispersive media. It is of great necessity to incorporate a dispersion model into FDTD formulation for electromagnetic wave analysis of dispersive media. Recently, it was reported that the modified Lorentz model can cover Debye, Drude, Lorentz, critical point, and quadratic complex rational function models. In this work, it is illustrated that the modified Lorentz model can also cover the complex-conjugate pole-residue model which is one of the most popular dispersion models. Modified Lorentz-based dispersive FDTD has not been thoroughly studied, especially for numerical aspects. In this work, we investigate auxiliary differential equation (ADE)-FDTD formulations for the modified Lorentz model based on electric flux density (D), current (J), or polarization (P). We perform a comprehensive study on memory requirement, the number of arithmetic operations, numerical stability, and numerical permittivity for the above three ADE-FDTD formulations. In addition, the bilinear transformation (BT) is incorporated into modified Lorentz-based FDTD formulations and it will be shown that the utilization of the BT can lead to better performance in terms of numerical stability and numerical accuracy. Numerical examples are used to demonstrate our work. [ABSTRACT FROM AUTHOR]

In this work, we present a dispersive finite-difference time-domain (FDTD) algorithm using a four-pole complex rational function (CRF). For the sake of a better curve fitting of the four-pole CRF dispersion model, we use a particle swarm optimization technique. We also discuss an efficient memory storage strategy using a state-space approach. The numerical aspects of four-pole CRF-FDTD, the numerical accuracy and the numerical stability, are investigated in detail. Numerical examples are used to validate four-pole CRF-FDTD and numerical stability issues are discussed in detail. We also discuss the computational accuracy and the computational efficiency of an arbitrary $N$-pole CRF-FDTD. [ABSTRACT FROM PUBLISHER]