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.

IFDDA, an Easy-to-use Code for Simulating the Field Scattered by 3D Inhomogeneous Objects in a Stratified Medium: Tutorial

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

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

A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

A Modified PML Implementation in FDTD for Handling Dispersive Media

Fast and accurate solution of scattering from 3-D object with large aspect ratio

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

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

Efficient Implementation of the Invariant Imbedding T-matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Simulation of the Inhomogeneous Medium with a Self-Adapting Grid.

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

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

A Novel Method to Simulate the Electromagnetic Susceptibility Problem

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

An enhanced finite difference time domain method for two dimensional Maxwell's equations

DPO-SDF: Differentiable Physical Optics Method with Signed Distance Field for 3-D Reconstruction of Perfect Electric Conductors

A 3D Finite Volume Method Solver for Scattering from Lossy Objects in Layered Media

Accelerated 3D Maxwell Integral Equation Solver using the Interpolated Factored Green Function Method

Fast finite difference solver for optical microscopy in deep biological tissue

The 3DV-FMM Method to Calculate Scattering from Inhomogeneous Dielectric Objects