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.

Mie scattering calculation by FDTD employing a modified Debye model for Gold material

An Area-Averaged Conformal Technique for the Electromagnetic Simulations of Debye Medium

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

Development of a Novel FDTD (2, 4)-Compatible Conformal Scheme for Electromagnetic Computations of Complex Curved PEC Objects

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

The Validation of the Simulated HIRF Effects in Metallic Cases Using Parallel FDTD Solver on a High Performance Computer

Hybrid FDTD Method Based on Conformal Technique for Modeling Wedges Composed of Conductive and Dielectric Materials

A Modified PML Implementation in FDTD for Handling Dispersive Media

Finite-Element-Based Generalized Impedance Boundary Condition for Modeling Plasmonic Nanostructures

A modified high-order FDTD conformal technique for modeling perfectly conducting objects

Optimized Conformal FDTD (2, 4) Method for Calculating Reflected and Diffracted Electromagnetic Fields of Perfectly Conducting Wedges

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

FDTD modeling of 3D metal-LTCC structures for RF(MM)ICs

FDTD far-field scattering amplitudes: Comparison of surface and volume integration methods

Modified Debye Model Parameters Of Metals Applicable For Broadband Calculations

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

Electromagnetic scattering analysis using a combined magnetic field integral equation for small objects with flat surfaces.

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

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

COMPUTATIONAL INVESTIGATION OF NANOSCALE SEMICONDUCTOR DEVICES AND OPTOELECTRONIC DEVICES FROM THE ELECTROMAGNETICS AND QUANTUM PERSPECTIVES BY THE FINITE DIFFERENCE TIME DOMAIN METHOD (INVITED REVIEW)