Abstract:Dispersive media refer to a class of natural substances, encompassing living organisms, composite materials, plasma and water, among others, with diverse applications in areas such as biomedicine, microwave sensing, electromagnetic protection, and stealth technology. In the pursuit of investigating the electromagnetic properties of these media, time-domain numerical methods, including finite difference in time domain (FDTD), finite element method (FEM), and time domain boundary integral equation method (TD-BEM), have been widely utilized. Time-domain numerical methods have been preferred over their frequency-domain counterparts due to their capability to handle nonlinear and wideband problems, as well as various material properties. The FDTD method, in particular, is a highly adaptable, robust, and easy-to-use numerical method that directly solves the Maxwell equations while also simulating the reflection, transmission, and scattering of electromagnetic waves in complex dispersion media. Nonetheless, the traditional FDTD method suffers from low computational efficiency arising from the CFL (Courant-Friedrichs-Lewy) stability condition. To address the issue of low computational efficiency, a new method, the complying divergence implicit finite-difference time-domain (CDI-FDTD) method with a one-step leapfrog scheme, is introduced for lossy Debye dispersive media. The Maxwell equations in the frequency domain form the starting point, and the Fourier transform is utilized to transform the electromagnetic field components from the frequency domain to the time domain. To approximate the integral terms arising from the frequency-to-time domain transformation, a recursive integration (RI) method is employed. Subsequently, the time-domain Maxwell equations and auxiliary variables are discretized with a one-step leapfrog implicit scheme. The iterative formula of the RI-CDI-FDTD algorithm for lossy Debye dispersive media is then derived. The RI-CDI-FDTD method does not alter the formulas of the traditional CDI-FDTD method while only requiring to add auxiliary variables for updating field components in the dispersive medium region. The numerical implementation is straightforward, and the electromagnetic modeling is flexible. Moreover, the unconditional stability of the RI-CDI-FDTD algorithm is proven using the von Neumann method. Finally, some numerical examples are presented to demonstrate the effectiveness and efficiency of the proposed method. In conclusion, our work contributes a crucial numerical simulation tool for accurately modeling complex dispersive media while providing a systemic stability analysis method for time-domain numerical methods.

A Symplectic FDTD Algorithm for the Simulations of Lossy Dispersive Materials

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

The Stability and Numerical Dispersion Analyses of High-Order Symplectic FDTD Scheme for Solving Time-Dependent Schrödinger Equation

A Conformal Scheme for Modeling Curved Dispersive Medium Objects Compatible with ADE-FDTD (2, 4) Method

A Modified PML Implementation in FDTD for Handling Dispersive Media

A Stability Extensible Symplectic FDTD(4,4) Method and Its Application in Subgrid Technique for Arbitrary Grid Ratios

Improved HIE-FDTD Method for Simulating Thin Layer Structure in Linear Dispersive Media

Application of the high-order symplectic FDTD scheme to the curved three-dimensional perfectly conducting objects

Application of the High‐order Symplectic FDTD Scheme to the Curved Three‐dimensional Perfectly Conducting Objects

A Novel Eigenvalue Method for Calculating the Band Structure of Lossy and Dispersive Photonic Crystals

Thin Absorber with Dispersive Medium Simulated by Conformal Hybrid Explicit Implicit FDTD Method

The HIE-FDTD Method for Simulating Dispersion Media Represented by Drude, Debye, and Lorentz Models

A Modular Implementation of Dispersive Materials for Time-Domain Simulations with Application to Gold Nanospheres at Optical Frequencies.

Ade-Wcs-Fdtd Method for General Dispersive Materials and Pml Implementation

Modeling and Simulation of Active Plasmonics with the FDTD Method by Using Solid State and Lorentz–Drude Dispersive Model

Time-domain Simulation for Plasmonic Devices

RI-CDI-FDTD method and program implementation for electromagnetic characteristics simulation of lossy Debye dispersive medium

Single Field HIE FDTD Solver with Nonuniform Grid Stratagem for More Fine Structure

ADE-Laguerre-FDTD Method for Wave Propagation in General Dispersive Materials