高阶色散介质的改进移位算子时域有限差分方法

提出一种适用于高阶Debye,Drude,Lorentz及其混合模型的改进移位算子的时域有限差分(SO-FDTD)方法。从介质极化率函数出发,将其写成一阶或二阶有理分式求和的形式,并在随时间步推进计算的过程中,通过引入中间变量和设置临时变量,克服了常规SO-FDTD将高阶模型直接转化为有理分式所导致的计算复杂性和内存占用量大的问题。同时,改进SO-FDTD方法的时域推进计算步骤具有通用性,克服了常规递归卷积(RC-FDTD)方法对各种高阶模型具有不同计算公式,因而不能形成通用计算程序的问题。最后,通过空气-高阶色散介质界面的反射系数计算验证了算法的有效性和通用性。     

高阶色散介质的改进移位算子时域有限差分方法

提出一种适用于高阶Debye,Drude,Lorentz及其混合模型的改进移位算子的时域有限差分(SO-FDTD)方法。从介质极化率函数出发,将其写成一阶或二阶有理分式求和的形式,并在随时间步推进计算的过程中,通过引入中间变量和设置临时变量,克服了常规SO-FDTD将高阶模型直接转化为有理分式所导致的计算复杂性和内存占用量大的问题。同时,改进SO-FDTD方法的时域推进计算步骤具有通用性,克服了常规递归卷积(RC-FDTD)方法对各种高阶模型具有不同计算公式,因而不能形成通用计算程序的问题。最后,通过空气-高阶色散介质界面的反射系数计算验证了算法的有效性和通用性。     

基于半解析递归卷积的通用色散介质FDTD方法

基于数字信号处理中的半解析递归卷积（SARC）算法,提出了一种用于色散介质电磁特性分析的半解析递归卷积时域有限差分方法（SARC FDTD）. 该方法既保持了FDTD方法处理复杂目标的灵活性,又吸取了线性系统SARC算法的绝对稳定性和高精度、低内存、高效率等优点,只需给出色散介质模型的极点和对应系数,即可运用SARC FDTD递推公式和通用程序进行计算,具有较好的通用性. 通过Debye,Drude和Lorentz三种常用色散介质模型对算法进行了验证. 关键词： 时域有限差分 色散介质 半解析递归卷积     

Calculation of the Effect on the Reflection of the Plane Electromagnetic Wave for Non-Magnetized Plasma with Different Electron Density Distributions

In this paper, shift operator finite-difference time-domain (SO-FDTD) method is applied for the calculation of the dispersive medium. The high efficiency and accuracy of this method is verified by calculating the reflection of the plane electromagnetic wave impinging on a non-magnetized plasma slab with different electron density distributions. The results show that the average electron density only determines overall trends of the reflection, and different distributions affect the oscillating process of the reflection. If the average electron density maintains the same, the distribution of electron density with homogenous or inhomogeneous alternation will sharply take effects on the reflection. And magnitude of alternation of electron density affects the incident frequency directly when the reflection tends to uniform.     

等离子体的分段线性电流密度递推卷积FDTD算法

给出了一种新的色散介质的时域有限差分(FDTD)算法，称为分段线性电流密度递推卷积FDTD算法.利用电流密度和电场强度的卷积关系，给出了该算法的计算公式.通过计算碰撞冷均匀等离子体平板对电磁波的反射系数和透射系数，验证了该算法的高效性和高精度. 关键词： 色散介质 FDTD算法 电磁波     

SO-FDTD Analysis on Magnetized Plasma

In this paper, a novel finite-difference time-domain (FDTD) method with recursive relationships among operators is developed for magnetized dispersive medium, named as the shift operator FDTD(SO-FDTD). The dielectric property of magnetized dispersive medium is written as rational polynomial function, the relationship between D and E is deduced in time-domain. And its high accuracy and efficiency are verified by calculating the reflection and transmission coefficients of electromagnetic waves through a collision plasma slab.     

Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions

In this paper, we will present a unified formulation of discontinuous Galerkin method (DGM) for Maxwell's equations in linear dispersive and lossy materials of Debye type and in the artificial perfectly matched layer (PML) regions. An auxiliary differential equation (ADE) method is used to handle the frequency-dependent constitutive relations with the help of auxiliary polarization currents in the computational and PML regions. The numerical flux for the dispersive lossy Maxwell's equations with the auxiliary polarization current variables is derived. Various numerical results are provided to validate the proposed formulation.     

Interior penalty DG methods for Maxwell’s equations in dispersive media

In this paper, we develop a fully-discrete interior penalty discontinuous Galerkin method for solving the time-dependent Maxwell’s equations in dispersive media. The model is described by a vector integral–differential equation. Our scheme is proved to be unconditionally stable and achieve optimal error estimates in both L² norm and energy norm. The scheme is implemented and numerical results supporting our analysis are presented.     

FDTD analysis on the effect of plasma parameters on the reflection coefficient of the electromagnetic wave

In this paper, the dielectric constant of dispersive medium is written as rational polynomial function, and the relationship between D and E is derived in time-domain. It is named shift operator FDTD (SO-FDTD) method. Compared to the analytical solution, the high accuracy and efficiency of this method is verified by calculating the reflection coefficient of the electromagnetic wave through a cold plasma slab. The effect on reflection coefficient is calculated by using the SO-FDTD method. The result shows that some factors effect on reflection coefficient. They are as follows: plasma thickness, electron density, electron distribution and incident frequency. And on most conditions, parabola distribution helps reduce reflection coefficient more effectively than homogeneous distribution.     

The comparison of the CPU time between SO-FDTD method and PLCDRC method for dispersive media

The precision and computational efficiency of the shift operator finite-difference time-domain (SO-FDTD) method are compared with those of the piecewise linear current density recursive convolution (PLCDRC) method by simulating the electromagnetic wave at different frequency propagation in homogeneous unmagnetized plasma in this paper. The results show that the two methods almost have the same precision, but the CPU time of the SO-FDTD method is less than that of the PLCDRC method, this is because the SO-FDTD method eliminates many convolution operations of the PLCDRC method.     

Numerical analysis of propagation characteristics of electromagnetic wave in lossy left-handed material media

The propagation characteristics of electromagnetic wave in lossy left-handed materials (LHM) are studied using finite-difference time-domain (FDTD) method base on auxiliary differential equation (ADE) technology. The LHM medium is realized with lossy Drude models for both the negative electric permittivity and the negative magnetic permeability. The discretized ADE-FDTD equations are derived in detail. The incident wave used in the simulation is a multiple cycle m-n-m pulses source. The term of Poynting's vector E_xH_y was calculated. These numerical results demonstrate conclusively that the phase velocity direction of electromagnetic wave propagation and the direction of the Poynting vectors are anti-parallel in LHM. The amplitude of electric field is reduced with the enhancive distance of LHM slab. It is also demonstrated that the energy of electromagnetic wave in the LHM slab is obviously attenuated, and the attenuation of energy becomes stronger with the angular plasma frequency ω_p increasing. These results indicate that LHM stealth is effective in theory, and reasonable selection of the large negative index of refraction can greatly enhance its effectiveness.     

Modeling and simulation of multiple resonant Lorentz dispersive materials for integrated optics

Complex materials are of increasing interest in Finite-Difference Time-Domain (FDTD) modeling. For example, when the particle density becomes large, collisional fluid models of plasmas are an attractive alternative to particle in cell methods. Further, frequency dispersive meta-materials are of increasing interest. Thus, Finite-Difference Time-Domain (FDTD) models are derived for magnetized plasmas and for the Lorentz and Drude material models. This paper models the Lorentz material using the FDTD tool. The modal analysis and various types of simulation results are also presented for these materials.     

一种处理Cole-Cole色散媒质的FDTD改进方案

对柯尔-柯尔(Cole-Cole)色散媒质的时域有限差分(FDTD)方案尝试给出改进:(1)保留了Cole-Cole经验公式中的静态电导率项;(2)扩展到三维情形;(3)人体乳房组织的模型参数源于实验测量.数值算例结果初步证实改进方案的可行性和有效性.改进方案适用于模拟更一般柯尔-柯尔色散媒质中的电波传播.     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

A note on Snell laws for electromagnetic plane waves in lossy media

Based on Maxwell's equations and Ohm's law, we rederived the Snell laws for reflection and transmission of harmonic inhomogeneous plane electromagnetic (EM) waves propagating through planar lossy interfaces. The present results are new, simple and exact and they recover the ordinary Snell laws in the case of lossless media. Besides, these results show that the wave propagation direction strongly depends on the polarization state of the EM wave and the lossy media can behave as a polarizing device. Moreover, we verify that in low frequency regime these traveling waves do not exhibit total internal reflection at interfaces between two adjacent lossy media.     

一种高效计算各向异性磁化等离子体的时域有限差分算法

利用电流密度和电场强度的卷积关系,引入电流密度和电场强度分段线性近似,导出一种适合各向异性磁化等离子体介质的FDTD算法的计算公式.计算各向异性碰撞磁化等离子体平板对平行于磁场传播的电磁波的反射和透射系数,与解析结果比较,验证该算法的高效性和高精度,与电流密度卷积(JEC)算法和分段线性电流密度卷积(PLCDRC)算法相比,计算精度及计算效率显著提高.此外,用该算法验证了法拉第旋转效应.     

A Hybrid Technique Combining FDTD and Series Solution for Near-to-Far-Field Transformation

A new near-to-far-field transformation that combines the FDTD and the series solution is presented. In this method, near-field data is calculated by FDTD and far-field data is also expressed by series solution. The boundary condition is applied to a virtual surface, which encloses arbitrarily shaped scatterers. In order to verify this method, scattering patterns from this transformation are compared with scattering patterns from exact solution. Numerical results show that errors in the far-field data from this method is less than surface equivalence theorem. Thus, this method is valid for obtaining the far-field data.     

Three-Dimensional FDTD Analysis of Chiral Discontinuities in the Waveguide

A simple finite difference time domain (FDTD) scheme for chiral media is developed for modeling the chiral discontinuities in the waveguide with the aid of the second-order backward finite difference. The accuracy of the FDTD scheme is demonstrated by the comparison between our calculated results and those in the literature for achiral and chiral media. Comparisons of characteristics between chiral and achiral discontinuities in the waveguide are presented at last.