PML吸收边界的时域有限差分法的数值色散研究

PML吸收边界以其优异的吸收能力与效果而倍受人们的关注。本文针对运用PML吸收边界的时域有限差分法的数值色散问题进行了研究,并得到了较为满意的结果,PML吸收边界在有效减少电磁 波在边界上的反射的情况下,并没有带来对数值色 散的不良影响。     

时域有限差分法的Matlab仿真

介绍了时域有限差分法的基本原理，并利用Matlab仿真，对矩形波导谐振腔中的电磁场作了模拟和分析。     

二维时域有限差分法对接地板宽度有限的共面波导的特性参数分析

本文利用二维时域有限差分（FDTD）法分析接地板宽度有限的共面波导（FGCPW）在高达1THZ频段内的频变特性。将计算结果与三维FDTD法分析结果相比较发现，其一致性较好，且大大提高了仿真的效率，通过比较拉地板宽度对FGCPW的色散和衰减特性的影响，说明设计共面波导的关键在于对几何尺寸参数的正确选择。文中还引出各参数间的一般约束关系。     

时域有限差分法模拟量子阱红外探测器光栅的光耦合

由于量子选择定则的限制,对于量子阱红外探测器(QWIP),必须利用衍射光栅增强其光学偶合效率,本文给出了一种基于时域有限差分法(FDTD)的数值方法,计算制备在QWIP器件上的金属光栅的衍射效应,模拟计算的结罘表明,FDTD方法是解析这种复杂结构内电磁场问题的有效手段．可以计算QWIP器件内各点电磁场所有分量的详细分布,进而可以估算衍射光栅的偶合效率,以及优化QWIP结构设计。     

时域有限差分法分析混晌室中场的均匀性

针对目前仿真分析混响室内的场分布过程中计算量大、耗费时间长和对计算机配置要求高的现状,用在一定的约束条件下随机分布于球面上的一组等效源辐射的平面电磁波叠加而成的一维推进的重叠平面渡模拟搅拌器搅拌过程中混响室内的电磁波,并采用时域有限差分法(FDTD)法对基于上述模型的混响室内的场分布进行了仿真计算,得出了使混响室内场分...     

一种非条件稳定的隐式时域有限差分法

介绍一种基于交替方向隐式(ADI)技术的时域有限差分法(FDTD).该方法是非条件稳定的,时间步长不再受到Courant稳定条件的限制,而是由数值色散误差来确定.与传统的FDTD相比,ADI-FDTD增大了时间步长,从而缩短了总的计算时间,特别是当空间网格远小于波长时,优点更加突出.首次把完全匹配层(PML)边界条件应用到ADI-FDTD计算中,采用幂指数形式的时间步进算法,推导了相应的迭代公式.进行了实例计算,并与传统FDTD的结果对比,验证了ADI-FDTD的有效性与优越性.     

时域有限差分法的收敛性

本文通过计算实例讨论了时域有限差分法的收敛性问题。计算了在平面波照射下无限长介质柱中的电磁场分布。把不同空间和时间步长时的计算结果,与用解析法所获得的结果进行比较,证明时域有限差分法有较好的收敛特性。对具有不同空间分辨率的人体二维非均匀块状模型中吸收电磁能量分布的计算表明,在高度非均匀介质中的场分布,仍然保持较好的收敛性。这些讨论从一个侧面提供了选取网格步长的依据。     

新型分裂步长时域有限差分法

提出一种新型的分裂步长时域有限差分(NSS-FDTD)法,并对其数值色散进行分析。该方法基于Split-Step方案和Crank-Nicolson方案,采用新的矩阵分解形式,与传统的FDTD算法、SS-FDTD算法相比,减少了计算复杂度。新型算法的推导程序简单,且具有良好的数值色散特性,还加入了一阶Mur吸收边界条件,给出一阶Mur吸收边界差分方程。将数值实验的结果和传统FDTD方法及理论值进行比较,数值结果一致性较好。     

二维金属型光子晶体的有限差分法分析

运用空域有限差分法对二维金属型光子晶体的带隙特性进行了分析,介绍了在周期边界以及金属柱面边界条件下,光子晶体单元内差分矩阵方程的建立过程,通过对所得差分矩阵行列式等于零的特征方程的求解,得出不同周期边界条件下相应的本征频点,从而求得光子晶体的能带结构图。对不同结构的二维金属型光予晶体进行了计算,计算结果经时域有限差分法（FDTD）给出了有效验证。     

色散金属的时域有限差分方法

表面等离子体激元具有使光场局域化和局域电磁场增强等特性,在纳米光子学和微观检测等诸多领域显示出广泛的应用潜力.时域有限差分（FDTD）数值计算方法能仿真激光与亚波长金属微结构相互作用的表面等离子体效应.金属具有色散性质,其相对介电常数模型有Drude模型和Lorentzs模型及它们相结合的Drude-Lorentzs模型,能拟合金属在可见光和近红外部分或全部波段的色散特性.FDTD数值计算要采用增加辅助变量和相应的辅助差分方程的方法使FDTD迭代计算稳定.     

A three-dimensional angle-optimized finite-difference time-domain algorithm

We present a three-dimensional finite-difference time-domain (FDTD) algorithm to minimize the numerical dispersion error at preassigned angles. Filtering schemes are used to further optimize its frequency response for broad-band simulations. A stability analysis of the resulting FDTD algorithm is also provided. Numerical results show that the dispersion error around any preassigned angle can be reduced significantly in a broad range of frequencies with small computational overhead.     

A novel finite-difference time-domain wave propagator

A novel time-domain wave propagator is introduced. A two-dimensional (2-D) finite-difference time-domain (FDTD) algorithm is used to analyze ground wave propagation characteristics. Assuming an azimuthal symmetry, surface, and/or elevated ducts are represented via transverse and/or longitudinal refractivity and boundary perturbations in 2-D space. The 2-D FDTD space extends from x=0 (bottom) to x→∞ (top), vertically and from z→-∞ (left) to z→∞ (right), horizontally. Perfectly matched layer (PML) blocks on the left, right, and top terminate the FDTD computation space to simulate a semi-open propagation region. The ground at the bottom is simulated either as a perfectly electrical conductor (PEC) or as a lossy second medium. A desired, initial vertical field profile, which has a pulse character in time, is injected into the FDTD computation space. The PML blocks absorb field components that propagate towards left and top. The ground wave components (i.e., the direct, ground-reflected and surface waves) are traced longitudinally toward the right. The longitudinal propagation region is covered by a finite-sized FDTD computation space as if the space slides from left to right until the pulse propagates to a desired range. Transverse or longitudinal field profiles are obtained by accumulating the time-domain response at each altitude of range and by applying the discrete Fourier transformation (DFT) at various frequencies     

Numerical analysis of traveling-wave photodetectors' bandwidth using the finite-difference time-domain method

We present full-wave analysis of traveling-wave photodetectors (TWPDs) using the finite-difference time-domain (FDTD) method. Impulse response in the frequency domain is obtained after time-domain data are calculated by the FDTD method. The impulse response includes the optical field profile, carrier transit time, microwave loss, microwave dispersion, and velocity mismatch all together. Three-decibel bandwidth is analyzed with the thickness of an i-layer and waveguide width as the design parameters. It is shown how transit time and microwave characteristics affect the bandwidth according to the TWPD's length. Three-decibel bandwidth is dominated by carrier transit time in case the device length is shorter than 300-500 /spl mu/m under the conditions given in this paper. However, if the device length gets longer, microwave characteristics affect the bandwidth.     

An isotropy-improved nonstandard finite-difference time-domain method

By introducing extra degree of freedom into the space-domain stencils, an isotropy-improved nonstandard finite-difference (NSFD) (2,2) method is presented in this paper. The technique is also extended to the higher-order (2,4) stencils. Both analytical and numerical results are presented and compared to demonstrate the accuracy and the strengths of these methods.     

一种高阶辛时域有限差分法的研究

从电磁场方程的Hamilton函数出发,提出了一种基于辛时域积分技术的高阶时域有限差分方法.该方法对时域的离散采用了能够保证系统的相空间体积不变和总能量不变的辛格式,对于空间的离散采用中心差分格式.计算结果表明与传统的时域高阶差分方法--Runge-Kutta法比较,该方法计算速度和计算精度都有较大的提高.     

Fast finite-difference time-domain analysis of resonators usingdigital filtering and spectrum estimation techniques

The use of digital filtering and spectrum estimation techniques for improving the efficiency of the FD-TD algorithm in solving eigenvalue problems is discussed. The great improvement of the efficiency of the method is demonstrated by means of both numerical and measurement results. In addition, several improvements to the present FD-TD method for eigenvalue analysis are presented. These include the analysis of open dielectric resonators and the extraction of the resonant frequencies from the FD-TD results. The result for the open dielectric resonator analysis is validated using measured data     

Stability condition for body-of-revolution finite-difference time-domain method

Although the body-of-revolution finite-difference time-domain method has been successfully applied to solve Maxwell's equations for rotationally symmetric structures, the stability condition of this method is inferred empirically. A new stability condition is given, which is derived analytically, and it demonstrates that the empirical stability condition is inaccurate. It is shown that the new stability restriction is weaker than the empirical one.     

非嵌入式多项式混沌法在爆轰产物JWL参数评估中的应用

多项式混沌展开(polynomial chaos expansion, PCE)方法对于分析随机等离子体电磁传播不确定性已经展示出了相当大的潜力。然而,由于构建多项式混沌模型的计算量随着不确定性输入维数的增加呈指数增长,数值模拟耗时长,导致“维数灾难”问题。因此,本文基于非侵入式多项式混沌(non-intrusive polynomial chaos,NIPC)方法,采用混合蒙特卡洛(Monte Carlo, MC)/PCE方法研究了多层等离子体平板电子密度不确定性对透射系数的影响,并验证了所提出方法的实用性。与MC方法相比,本文方法可以大大加快仿真的速度,有效缓解了多项式展开项的数量随着随机变量维数的增加而快速增加的缺点,同时大大降低了MC方法的仿真次数,有利于电磁模型的多参数不确定性分析。

     

Hybrid finite-volume finite-difference time-domain modelling ofcurved surfaces

A computationally efficient hybrid finite-volume finite-difference time-domain (FV-FDTD) method is proposed to model the scattering from curved objects. The hybrid FV-FDTD method uses a local conforming mesh consisting of only two FVTD cells, one either side of the scatterer's contour, giving a better approximation to the surface. The improved definition of the scatterer using the hybrid FV-FDTD method is shown to improve the solution in comparison to the FDTD method     

A finite-difference time-domain method applied to anisotropicmaterial

The popularity of the finite-difference time-domain (FDTD) method stems from the fact that it is not limited to a specific geometry and it does not restrict the constitutive parameters of a scatterer. Furthermore, it provides a direct solution to problems with transient illumination, but can also be used for harmonic analysis. However, researchers have limited their investigation to materials that are either isotropic or that have diagonal permittivity, conductivity, and permeability tensors. The authors derive the necessary extension to the FDTD equations to accommodate nondiagonal tensors. Excellent agreement between FDTD and exact analytic results is obtained for a one-dimensional anisotropic scatterer