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.     

等离子体的SO-FDTD算法和对电磁波反射系数的计算分析

将一类色散介质的介电常数写成有理分式函数形式，进而导出FDTD中电位移矢量D和电场强度E之间的关系，形成SO-FDTD方法. 应用该方法计算了冷等离子体平板对电磁波的反射系数，通过与解析结果的比较，验证了该算法的高效性和高精度，同时，应用SO-FDTD方法，计算了等离子体层对垂直入射电磁波的反射系数，结果表明：等离子体厚度、电子密度、电子密度的分布形式和入射波频率是影响反射系数的重要因素. 关键词： 等离子体 电磁波 FDTD方法     

Finite-Difference Time-Domain Algorithm for Dispersive Media Based on Runge-Kutta Exponential Time Differencing Method

The electromagnetic propagation in dispersive media is modeled using finite difference time domain (FDTD) method based on the Runge-Kutta exponential time differencing (RKETD) method. The second-order RKETD-FDTD formulation is derived. The high accuracy and efficiency of the presented method is confirmed by computing the transmission and reflection coefficients for a nonmagnetized collision plasma slab in one dimension. The comparison of the numerical results of the RKETD and the exponential time differencing (ETD) algorithm with analytic values indicates that the RKETD is more accurate than the ETD algorithm.     

New FDTD-based method for computing the Lorentz dispersive media

We apply the shift-operator finite-difference time-domain (SO-FDTD) method to the calculation of the Lorentz dispersive media. The Lorentz-medium permittivity is written as a rational fractional function, and the relation between D and E is derived in a time domain. The high efficiency and accuracy of the method is demonstrated by calculating the reflection and transmission of the electromagnetic plane wave impinging on a Lorentz-dispersive-medium slab. The results obtained with the help of the SOFDTD method are compared with analytical results. The SO-FDTD method is proved to have high efficiency and accuracy in computing the Lorentz dispersive media.     

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.     

弱电离尘埃等离子体层反射与透射的SO-FDTD方法分析

将弱电离尘埃等离子体复介电常数表示以jω 为自变量的分式多项式形式, 利用移位算子时域有限差分(SO-FDTD)方法, 给出了弱电离尘埃等离子体频域本构关系的时域形式. 推导了弱电离尘埃等离子体电磁特性的时域递推关系式. 用SO-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.     

Analysis of the temperature influence on the plasma reflectance by a new FDTD method

We apply the shift-operator FDTD (SO-FDTD) method for calculating dispersive media plasmas. We write the dielectric constant of a dispersive medium as a rational fractional function and derive the relation between D and E in the time domain. For a uniform, approximately isothermic, isotropic, warm-collision-plasma slab, the reflection coefficient of electromagnetic waves can be calculated using the SO-FDTD method. We study the influence of temperature on the collisions and absorption of electromagnetic waves in the plasma and discuss the regularity of the relationships obtained.     

各向异性磁化等离子体的SO-FDTD算法

给出了一种新的计算各向异性磁化色散介质的有限差分(FDTD)算法，称为移位算子FDTD(SO-FDTD)算法，它利用算子之间的移位递推关系，将一类色散介质的包含介电常数的表达式写成有理分式函数形式，进而导出FDTD中一系列相关量之间的关系.通过计算各向异性等离子体平板对电磁波的反射系数和透射系数，验证了该算法的高效性和高精度，与JEC算法相比，可使计算效率提高数倍. 关键词： 磁化等离子体 电磁波 FDTD方法 各向异性     

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.     

Finite-difference time-domain studies of low-frequency stop band in superconductor-dielectric superlattice

<正>The transmission coefficients of electromagnetic(EM) waves due to a superconductor-dielectric superlattice are numerically calculated.Shift operator finite difference time domain(SO-FDTD) method is used in the analysis.By using the SO-FDTD method,the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers,from which a stop band starting from zero frequency can be apparently observed.The relation between this low-frequency stop band and relative temperature,and also the London penetration depth at a superconductor temperature of zero degree are discussed,separately.The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.     

Analysis of backscattering RCS of targets coated with parabolic distribution and time-varying plasma media

The backscattering radar cross section (RCS) of targets coated with parabolic distribution and time-varying unmagnetized plasma is analyzed by the finite-difference time-domain (FDTD) method based on the Runge–Kutta exponential time differencing (RKETD) technique. The two-dimensional RKETD-FDTD formulations are derived. The effects of the plasma parameters on the RCS are investigated. The numerical results illustrate that plasma cloaking system can successfully reduce the backscattering RCS, the plasma stealth is effective, and reasonable selection of different parameters of plasma can enhance its effectiveness.     

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

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

A FDTD formulation for high order dispersive media of muscle-equivalent

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 muscle slab. For human tissues where multiterm Debye relaxation equations must generally be used. We describe a new differential equation approach, which can be used for general dispersive media. In this method D(t) is expressed in terms of E(t) by means of a differential equation involving D, E, and their time derivatives. The method is illustrated by means of example of media for which relative permittivity is given by a multiterm Debye equation, and for an approximate two-thirds muscle-equivalent model of the human body.     

时空非均匀等离子体鞘套中太赫兹波的传播特性

高超声速飞行器再入地面的过程中,其周围等离子体的电子密度是非均匀且随时间变化的.对于不同的再入高度,飞行器周围的温度和压强也会发生改变.因此,研究电磁波在时空非均匀等离子体鞘套中的传播特性意义重大.首先建立了时变非均匀的等离子体鞘套模型,然后通过经验公式得到温度、压强与碰撞频率三者的关系.采用时域有限差分方法计算了太赫兹波段中不同电子密度弛豫时间、温度、压强时的反射系数、透射系数和吸收率.研究结果表明:在太赫兹波段中,电子密度的弛豫时间越长,温度越高,压强越大,电磁波越容易穿透等离子体;弛豫时间越短,温度越低,压强越小,等离子体对电磁波吸收率的变化越明显.这些结果为解决"黑障"问题提供了理论依据.     

一种基于数字信号处理技术的改进通用色散介质移位算子时域有限差分方法

借鉴数字信号处理技术中的无限脉冲响应滤波器实现的思想,提出一种适用于通用色散介质的改进移位算子时域有限差分（shift operator finite-difference time-domain,SO-FDTD）方法.与原SO-FDTD方法相比,改进的SO-FDTD方法内存占用减少33%以上,同时计算时间也有所减少.最后通过一阶、二阶色散介质雷达散射截面计算验证了方法的通用性和有效性. 关键词： 时域有限差分 色散介质 移位算子     

耦合势有限体积法高效模拟各向异性地层中海洋可控源的三维电磁响应

本文基于电场矢势与标势分解的耦合势有限体积法研究建立一套各向异性地层中海洋可控源电磁法的三维响应的高效数值模拟技术.首先引入电场的矢势和标势,将电场分解为无散场和无旋场之和,Maxwell方程转换为关于矢势与标势的混合Helmholtz方程,克服低感应数问题.在此基础上,借助Yee氏交错网格和有限体积法以及非均质单元中等效电导率公式,建立混合Helmholtz方程的离散方程.并采用直接法求解器PARDISO求解离散方程,有效保证在大的求解空间中仍然能够获得电磁场稳定可靠的数值解.此外,在数值模拟中利用差异场技术,克服源的奇异性问题,尽可能提高近场的计算精度.与解析解的对比证明了该算法的有效性.数值模拟结果表明,海洋可控源电磁法沿测线方向的电场,对油气藏的纵向电阻率敏感,对横向电阻率不敏感;对油气藏上方的覆盖层的纵向电阻率和横向电阻率都敏感.     

磁化铁氧体覆盖二维导体目标FDTD分析的移位算子方法

采用移位算子方法把时域有限差分法推广应用于二维磁各向异性色散介质—磁化铁氧体中.证明了电磁波横向入射二维轴向磁化铁氧体目标情形下,电磁波可按目标的轴向分解为横电波(TE波)和横磁波(TM波),且TE波的散射特性与铁氧体介质无关,而TM波的散射特性与介质电磁参量密切相关,同时对其物理原因进行了分析.通过采用移位算子方法处理磁化铁氧体频域本构关系,得到该情形下目标电磁散射的移位算子时域有限差的迭代计算公式,同时解决了电磁波在各向异性和频率色散介质中传播的问题.计算了轴向磁化铁氧体涂敷VonKarman型导体柱的TM波双站雷达散射截面,分析了铁氧体参量对目标双站雷达散射截面的影响.结果表明:恰当地选择铁氧体参量能有效地减少目标的雷达散射截面,本文时谐因子取exp〔jωt〕.     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

磁化铁氧体电磁散射的移位算子FDTD分析

将移位算子（shift operator,SO）-时域有限差分（fintie difference time domain,FDTD）方法推广到外磁场沿任意方向时磁化铁氧体的情形.从磁化铁氧体的相对磁导系数张量和各向异性情形下Maxwell旋度方程出发,推导出外磁场沿任意方向情形下磁化铁氧体电磁散射的SO-FDTD迭代公式.应用该方法计算了磁化铁氧体层、磁化铁氧体球及含磁化铁氧体涂层的Von Karman形雷达罩电磁散射,数值计算结果表明了该算法的正确有效性. 关键词： 磁化铁氧体 FDTD 移位算子法