Search citation statements
Paper Sections
Citation Types
Year Published
Publication Types
Relationship
Authors
Journals
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.
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.
The employment of uniaxial/biaxial bianisotropic materials in the field of optical device manufacturing is widespread due to the distinctive electromagnetic response characteristics exhibited by such materials. To effectively analyze the electromagnetic properties of uniaxial/biaxial bianisotropic materials, rapid-transfer matrix method (R-TMM) is proposed to investigate the propagation process of plane waves in the media. Starting from the Maxwell’s equations in the time domain, the homogeneous differential equation about the electric field is constructed by processing the matrix containing dielectric and magnetic conductivity, electric and magnetic loss, tellegen and chirality carrier parameters, and the complex matrix operation is applied to that equation to get the booker quartic equation, and then the formulae method are utilized to obtain the eigenvalues in the uniaxial/biaxial bianisotropic media. Subsequently, the tangential continuity of layered media at the interface is employed to establish a transfer matrix for single-layered media. In the case of multi-layered media, the transfer matrix of plane waves propagating in multi-layered uniaxial/biaxial bianisotropic media can be obtained by means of a continuous iteration process based on the transfer matrix of single-layered media. The formula for calculating the propagation coefficient of uniaxial/biaxial bianisotropic materials can be derived based on the different upward and downward waves in the reflection/transmission region. Finally, the reliability and efficiency of R-TMM is verified from two numerical experiments with the plane waves incident at different angles on a uniaxial/biaxial bianisotropic media. The first experiment is designed as a single-layered biaxial bianisotropic model with more general electromagnetic parameters, and the second experiment is designed as a double-layered uniaxial and biaxial bianisotropic model consisting of common optical materials, which are composed of two non-magnetic materials, lithium niobate (LiNbLO<sub>3</sub>) and cadmium sulfide (CdS). The experimental results demonstrate that, in comparison with the conventional conventional-transfer matrix method (C-TMM), the R-TMM reduces the computational memory and CPU time required for calculating the reflection and transmission coefficients of the uniaxial/ biaxial bianisotropic model by over 98%, while maintaining the accuracy of the reflection and transmission coefficients calculation. Therefore, R-TMM provides an efficient and dependable approach for the design of complex optical devices and the analysis for uniaxial/biaxial bianisotropic propagation characteristics.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2025 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.