三维FDTD中TF/SF边界的时域群速补偿法

TF／SF（Total Field／Scattered Field）法是FDTD（Finite-Difference Time—Domain）数值算法中引入平面波源的有效方法之一。本文对TF／SF法中的一维辅助渡源法进行改进，首先推导出三维FDTD网格内的电磁波传播群速和相速的近似解析解。为预测离散网格内电磁波传播提供理论基础；然后以此为依据，提出对一维辅助波源的时域群速补偿法。经过实际三维FDTD问题的数值模拟，证实了本文方法在保持一维辅助波源法高效简洁的同时，进一步减少了泄露误差，给FDTD计算散射问题提供了有力工具。     

基于Split-Field FDTD法的近埋地无限长散射体二维算法

提出了一种基于周期结构split-field FDTD法的近埋地无限长散射体二维算法.该方法根据散射体轴向均匀性将三维split-field FDTD法转化为二维算法,减少了内存和计算量,可分析斜入射脉冲波照射下近地、埋地无限长散射体散射问题.为了进一步减小计算量,连接边界上的入射波(地上为原始入射波和反射波的叠加,地下为透射波)采用一维FDTD法引入.吸收边界采用了UPML匹配层,导出了适用于split-field FDTD算法与有耗介质匹配的UPML方程.通过数值算例,验证了此二维算法应用于近埋地无限长散射体问题的有效性和实用性.     

三维FDTD方法中高效平面波源的引入

应用分裂平面波时域有限差分（Splitting Plane wave Finite Difference Time Domain,SP-FDTD）方法到三维时域有限差分（Finite Difference Time Domain,DFDTD）中引入高效平面波源.该方法基于分裂场思想,在一维FDTD上构造了新的迭代公式,使得一维FDTD和三维FDTD离散网格之间的数值相速度一致,消除了由于相速不一致而在总场区引起的泄漏误差以及插值带来的数值误差.通过数值算例验证了SP-FDTD方法对不同波源在任意角度（斜入射）下的平面波入射都是有效的,且泄露误差均在-300dB水平.     

垂直人射HEMP近地面电磁环境特性研究

高空核爆电磁脉冲(HEMP)具有覆盖范围广、峰值场强高等特点,对电子设备构成严重威胁.针对有关HEMP标准中主要涉及自由空间HEMP波形描述的情况,研究了双指数HEMP平面波垂直入射地面时,地面附近的电磁脉冲环境特性.基于时域有限差分(FDTD)法提出使用一维合成平面波解决总场空间激励源的引入问题.通过设置不同的入射波...     

一维光子晶体慢光波导的FDTD研究

设计了一种新型一维光子晶体慢光波导结构。在常规波导一侧进行了特殊的设计,使波导具有周期性结构,从而具有特殊的色散关系,获得慢光效应。基于麦克斯韦方程利用平面波展开法对光子晶体慢光波导的色散关系进行了分析,获得了波导模以及相应的慢光频率。并利用时域有限差分法（FDTD）对脉冲在波导的传播进行了时域上的模拟,对慢光效应进行验证。     

PML-FDTD及总场-散射场区连接边界条件在三维柱坐标系下的实现及应用

基于直接形式的麦克斯韦方程，在三维扩展柱坐标系下推导了完全匹配层吸收边界条件的时域有限差分表达式(UPML—FDTD)，将FDTD方法推广到三维柱坐标系中，并对其反射系数进行了测试。给出了柱坐标系下完整的总场一散射场区的连接边界条件以及平面波源的设置方法。最后利用这种方法模拟了平面波在自由空间中的传播过程，计算了理想导电圆柱体在TEM波照射下的散射方向图。     

太赫兹波在二维光子晶体波导中的传输特性研究

首先用平面波展开法(PWM)计算了二维光子晶体的能带结构,然后提出了适合THz 波传输的光子晶体波导模型,并采用时域有限差分法(FDTD)研究了THz波在这种波导中的传输特性.在波导的输入输出口采样场值经过傅立叶变换以后进行比较,结果很合理.分析结果表明,位于光子晶体禁带内的THz波在这种波导中的传输是几乎没有损耗的,这为开发性能优良的THz 器件提供了理论依据.     

平面波在三维完全匹配层中的传播特性

完全匹配层(PML)良好的吸波特性已成功地用于FDTD问题。Katz给出了三维问题分量场所满足的12个微分方程,但未对波在PML中的传播特性进行讨论。本文基于Katz的微分方程组,推导出了平面波解、波数及其阻抗。     

激励源对二维光子晶体传输特性的影响

用时域有限差分法（FDTD）分析了二维光子晶体的传输特性，研究了纯平面波源、高斯波源对传输特性的影响，指出在同一种激励源情况，二维光子晶体的传输特性与其入射角有关。     

FDTD计算中非均匀网格网络的分析及综合

本文依据微波网络理论及网格波阻抗的概念,对时域有限差分(FDTD)计算中的非均匀网格反射波进行分析,提出了网格网络分析和综合的概念,分析了网格网络的反射特性,给出了网格网络的综合过程.基于上述理论和概念,在一维(TEM波)和三维微带结构以及方波导中实现了有关的非均匀网格的网络特性,其特性可用于沿传播方向的网格非均匀性的匹配,减小网格反射波及改进波的传播特性等.     

FDTD Simulation of Electromagnetic Wave Transformation in a Dynamic Magnetized Plasma

Finite-Difference Time-Domain (FDTD) algorithms are applied to study the transformation of a pre-existing electromagnetic plane wave by prescribed time variation of a cold magnetoplasma with a magnetic field along the propagation direction. A one-dimensional FDTD code is used to verify the results obtained earlier using analytical approximations based on (a) WKB method for slow switching of the plasma medium and (b) Green's function technique for rapid switching of the plasma medium. A novel successive reduction method has been developed and applied to obtain the amplitudes and the frequencies of the new modes generated by the switching of the medium.     

瞬变等离子体中电磁波频率漂移特性研究

推导一维瞬变磁化等离子体的时域有限差分(FDTD)递推式,并以一维金属谐振腔作为计算模型,分析了瞬变等离子体中电磁波频率漂移特性。从理论分析出发,得到了一维瞬变等离子体对电磁波作用的解析解。通过选取相同的模型和参数,将FDTD数值解与解析解进行对比,验证了所用FDTD方法的准确性,在此基础上研究了瞬变等离子体中电磁波的频率漂移规律。     

一种新的FDTD入射场设置方法

本文提出了一种新的入射场设置方法.采用一维FDTD算法来模拟入射波.所得结果与通常的入射场解析表达式的结果比较表明,采用一维FDTD算法引入的入射波在总场区分布均匀,在散射场区泄漏小.     

FDTD入射波解析设置的约束条件

用ＦＤＴＤ方法求解电磁散射问题中，入射波的设置起着关键性作用。本文从ＦＤＴＤ零场初始条件出发，导出了入射波解析设置时的约束条件，该约束条件表明，在ＦＤＴＤ总场边界上要求入射波逐点延时设置，通过与自由空间及半空间散射数值例子对比，验证了该约束条件下入射波解析设置方法的正确性。     

网格长度对FDTD计算的影响

基于平面波传播的物理性质，提出一种FDTD改进方法——在研究稳态问题时，为保证数字波前与入射方向垂直，网格划分不能是固定的，其长度之比应随人射方向变化。研究了改进方法对数值色散的影响，并与传统方法作了比较。在三维Yee网格中对平面波的模拟结果表明，改进方法的收敛性能大为提高，对远场散射的误差影响小。     

Spectral FDTD: a novel technique for the analysis of oblique incident plane wave on periodic structures

This paper introduces a new technique which calculates the reflection coefficient for the plane wave incident on planar periodic structures. The method referred to as spectral finite-difference time-domain (SFDTD) replaces the conventional single-angle incident wave, with a constant transverse wavenumber (CTW) wave. Because the transverse wavenumbers are constant, the fields have no delay in the transverse plane (x-y plane), and PBC (periodic boundary condition) can be directly implemented in the time domain for both oblique and normal incident waves. The stability criterion for this new FDTD technique is angle-independent and therefore this method works well for incident angles close to grazing (/spl theta/=90/spl deg/) as well as normal incident (/spl theta/=0/spl deg/). This shows the efficiency of the method compared to other available FDTD techniques for the same purpose that force a more restricted stability criterion as angles turns to grazing. The validity of this method is verified by comparing the reflection coefficient calculated by this method with the analytical results of a grounded slab. The results of this technique are also compared with method of moments for a periodic array of metallic patches and a good agreement is observed. A periodic array of metallic patches above a PEC plate is analyzed and the reflection coefficient is calculated over a wide frequency band for angles varying from 0/spl deg/ to close to 90/spl deg/.     

Validity and Accuracy of Equivalent Circuit Models of Passive Inductive Meshes. Definition of a Novel Model for 2D Grids

Accuracy of equivalent circuit models of periodic grids is investigated in amplitude and phase in the visible region. The grids studied here are one-dimensional (1D) and two-dimensional (2D) inductive thin metal meshes. They are located in free space and are illuminated by a plane wave under normal incidence. The range of validity and the accuracy of conventional circuit models are defined by comparison with rigorous results obtained with the Finite-Difference Time-Domain (FDTD) method. In particular, it is shown that electrical models of 1D grids are accurate, whereas equivalent circuits of 2D grids should be used very cautiously. Then, a new formulation is proposed to overcome this major drawback. In the non-diffraction region, the agreement between our model and the FDTD results is within 2% for the power reflectivity and 1° for the phase over a very wide range of strip widths.     

FDTD Simulation of Electromagnetic Reflection of Conductive Plane Covered with Imhomogeneous Time-Varying Plasma

A finite-difference time-domain (FDTD) algorithm is applied to study the electromagnetic reflection of conduction plane covered with inhomogeneous, collision, warm, time-varying plasma. The collision frequency of plasma is a function of electron density and plasma temperature. Under the one-dimensional case, transient electromagnetic propagation through various plasmas have been obtained, and the reflection coefficient of EM wave through inhomogeneous time-varying plasma (ITVP), homogeneous time-varying plasma (HTVP) and inhomogeneous plasma (IP) are calculated under different conditions. The results illustrate that a plasma cloaking system can successfully absorb the incident EM wave.     

High-accuracy FDTD solution of the absorbing wave equation, and conducting Maxwell's equations based on a nonstandard finite-difference model

We previously introduced high-accuracy finite-difference time-domain (FDTD) algorithms based on nonstandard finite differences (NSFD) to solve the nonabsorbing wave equation and the nonconducting Maxwell equations. We now extend our methodology to the absorbing wave equation and the conducting Maxwell equations. We first derive an exact NSFD model of the one-dimensional wave equation, and extend it to construct high-accuracy FDTD algorithms to solve the absorbing wave equation, and the conducting Maxwell's Equations in two and three dimensions. For grid spacing h, and wavelength /spl lambda/, the NSFD solution error is /spl epsiv//spl sim/(h//spl lambda/)/sup 6/ compared with (h//spl lambda/)/sup 2/ for ordinary FDTD algorithms using second-order central finite-differences. This high accuracy is achieved not by using higher-order finite differences but by exploiting the analytical properties of the decaying-harmonic solution basis functions. Besides higher accuracy, in the NSFD algorithms the maximum time step can be somewhat longer than for the ordinary second-order FDTD algorithms.