二维单步交替方向隐式时域有限差分法吸收边界性能分析

给出了两种适用于二维单步交替方向隐式时域有限差分(2-D Leapfrog ADI-FDTD)方法的吸收边界:Mur边界和卷积完全匹配层(CPML)边界。单步交替方向隐式时域有限差分(Leapfrog ADI-FDTD)方法是一种无条件稳定的全隐式差分算法,由于二维空间Leapfrog ADI-FDTD的迭代同时存在显式和隐式方程,故而不同电磁分量的边界条件也存在差异。从原理出发,推导了适用于2-D Leapfrog ADI-FDTD方法的CPML边界条件,并与一阶Mur边界进行比较,利用自由空间的反射误差来表征两种边界的吸收性能,简要总结了两种吸收边界的优缺点。     

一种新时域交替隐式差分算法在散射问题中的应用

提出一种新时域交替隐式有限差分(ADI-FDTD)算法格式. 传统ADI-FDTD算法的 场量步进方程涉及周围若干网格的较多场量，导致两个区域的步进方程处理较困难：一个是 邻近完全匹配层(PML)和散射场交界区，另一个是邻近连接边界区. 特别是后者，考虑入射 波影响需对场量所在区域判断，根据不同情况对原有方程进行修正，一维和二维散射问题 相对简单，可三维问题修正式有数十种之多而几乎无法完成. 本方法基于分裂场形式的ADI- FDTD技术，使得散射场区和PML吸收层区的表达形式完全一致，从而忽略两者差别.另 关键词： 时域交替隐式有限差分算法 电磁散射     

平面光波导放大器中光场分布的一种计算方法

用交变隐式差分方向的时域有限差分法(ADI-FDTD)计算了掺Er³⁺玻璃平面光波导放大器中光场的分布.其算法简单,结果无条件稳定,能完全模拟波导中的光场特性.     

Er与Yb共掺波导放大器的ADI-FDTD数值模拟

提出了一种新的简单的数值计算方法模拟Er与Yb共掺磷酸盐玻璃波导放大器(EYCDWA)的增益特性.该方法是基于交变隐式差分方向的时域有限差分法(ADI-FDTD)的迭代方法,用ADI-FDTD计算波导中的光场强度分布,用迭代法研究光波的传输方程.数值结果表明,增加光波导长度和提高Er³⁺浓度是提高EYCD WA增益的两种重要途径,Yb³⁺起能量传递的作用,Er/Yb的比率取1~4较好.     

二维圆柱坐标系中五角形共形网格算法

 时域有限差分（FDTD）算法求解电磁问题时，如果目标结构不能和传统网格体系共形，往往采用局部共形网格或非正交网格算法。详细讨论了二维圆柱坐标系中五角形共形网格算法，该算法能方便地对倾斜薄层媒质或电导率不均匀的倾斜媒质等复杂结构建模,并以锥形电磁脉冲天线为例，验证了该算法的有效性。     

三维周期结构弱无条件稳定时域有限差分算法

从麦克斯韦方程出发,导出了适用于周期结构的三维弱无条件稳定时域有限差分（FDTD）算法的迭代式,在理论上证明了其稳定性条件,其稳定性条件比普通FDTD的稳定性条件要宽松,并将周期边界条件应用到迭代式中,得到了可直接编程迭代的方程,给出了对其特有的一类非三对角阵的解决办法。最后通过算例将计算结果与常规FDTD及ADI-FDTD计算结果比较,验证了算法的精确性和有效性。     

三维时变等离子体目标的电磁散射特性研究

基于电流密度拉普拉斯变换方法改进的时域有限差分(LTJEC-FDTD)算法, 研究时变等离子体目标的电磁散射特性.由Maxwell方程和等离子体本构方程出发, 利用拉普拉斯变换和拉普拉斯逆变换, 推导出计算三维时变问题的时域有限差分(FDTD)算法的迭代式. 采用模式匹配方法验证了FDTD迭代式的正确性, 并通过计算等离子体球的雷达散射截面(RCS)验证了算法相关边界的正确性. 最后用LTJEC-FDTD算法分析了涂覆时变等离子体的战斧式巡航导弹的RCS. 关键词： 时变等离子体 雷达散射截面 模式匹配方法 时域有限差分方法     

基于辛Runge-Kutta-Nystrom方法的雷达散射截面计算

从基本的差分概念和Maxwell方程出发, 引入电磁场方程的Hamilton函数.提出一种基于Runge-Kutta-Nystrom辛算法的高阶时域有限差分方法，该方法保持了系统的相空间体积不变和总能量不变，并导出了迭代公式.在此基础上计算了一种金属圆柱的雷达散射截面.计算结果表明该方法的正确性及快速、精确的特性. 关键词： 雷达散射截面 高阶算法 辛Runge-Kutta-Nystrom方法 时域有限差分     

辛时域有限差分算法研究等离子体光子晶体透射系数

相较于传统的时域有限差分算法,辛时域有限差分算法具有高准确度性和低色散性.传统的时域有限差分算法的计算准确度较低,数值色散误差较大,并且破坏了麦克斯韦方程的辛结构,从而导致其稳定性较差.然而辛时域有限差分算法可以克服这些缺点,从而保证了整个仿真计算的准确性和稳定性.本文基于辛时域有限差分算法,对等离子体光子晶体的带隙特性,透射系数等进行了研究,并与传统的时域有限差分算法进行了对比,验证了辛时域有限差分算法的优势和可行性.     

谱方法用于非定常流动计算的隐式求解

本文对谱方法用于周期性非定常流动的隐式求解方法进行了探讨,分析了影响计算稳定性和收敛速度的因素.提出了结合多重网格的隐式求解方法并对算法进行了验证,初步计算表明本文算法具有良好的稳定性和收敛速度.对于周期性非定常流动,结合本文提出的隐式求解的时域谱方法可以达到很高的精度且具有良好的计算效率.     

2D FDTD Algorithm for the Analysis of Lossy and Dispersive Millimeter-Wave Coplanar Waveguide

The dispersion and loss characteristics of millimeter-wave coplanar waveguides (CPWs) are analyzed by using an efficient two-dimensional finite-difference time-domain (2D FDTD) algorithm. Combined with graded mesh technique, alternating direction implicit technique and curve-fitting technique, this proposed algorithm verifies its accuracy and efficiency through an example of infinite-ground CPW. Moreover, the influence of ground-plane width on CPW's loss characteristics is analyzed. Bulk micromachined CPWs with low loss are also analyzed by the proposed algorithm.     

One-step alternating direction implicit FDTD algorithm

In this paper, a novel unconditionally stable alternating direction implicit finite-different time-domain method(ADI-FDTD) called the one-step ADI-FDTD method is presented, where the calculation for one discrete time step is performed using only one procedure, but not the original two sub-updating procedures. Consequently, the proposed one-step ADI-FDTD methods have consumed less computer memory and computation resources and have been faster than the conventional ADI-FDTD methods. We analytically and numerically verified that the new algorithm is unconditionally stable and free from the Courant condition.     

基于交替隐式有限差分法的快速早期乳腺癌时域微波断层成像

采用时域微波断层成像技术进行早期乳腺癌检测能够准确地获得乳房的电参数分布,具有明确的物理解释和医学诊断价值.临床应用讲究即时性,为了提高检测的速度,本文将交替隐式有限差分法应用到乳腺癌检测中,基于正反演时间步进成像算法进行成像分析,结果显示在保证精度的前提下,采用交替隐式有限差分法的成像时间可缩短为传统时域有限差分法的23%,提高了微波断层成像技术的临床可应用性.     

An Efficient Numerical Model for Planar Spiral Inductor On-Chip

In this paper, a very useful numerical technique has been developed for analyzing the transient characteristics of a planar-spiral inductor on-chip. A locally conformal technique and an alternating-direction implicit scheme are applied to the finite-difference time-domain method. A formulation for solving three dimensional Maxwell’s equations is proposed. Using the proposed method, various parameters of the planar-spiral inductors have been analyzed and an equivalent circuit, which includes frequency-independent circuit elements, has been introduced. Highly computational efficiency is implemented. Numerical results show excellent agreement with the measured data over a wide frequency range.     

NUMERICAL MODELING OF OPTICALLY CONTROLLED DIELECTRIC REASONATORS

To analyze the optically controlled dielectric resonators accurately, an efficient numerical modelling technique is proposed in this paper. By using alternating-direction implicit finite-difference time-domain method and conformal technique, the resonant frequency of a dielectric resonator is calculated. In addition, an optical generation of plasma is used as a possible means of controlling the resonate frequency, and the effect that solid state plasmas have on the resonator’s frequency is described. The numerical results agree very well with measurements in estimation of the optically induced resonator’s frequency-shift.     

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.     

A Study of Infrared Photonic Crystal Devices Using New ADI-FDTD Algorithm

To investigate the infrared photonic crystal devices numerically, the traditional finite-difference time-domain (FDTD) method has been modified by combining with a new alternating direction implicit (ADI) algorithm. An improvement of two-five in speed over previous FDTD methods can be obtained by calculating the envelope rather than the fast-varying field, and the numerical errors are minimized. Consider the isolated localized coupled-cavity modes, the phenomenon of eigenmode splitting has been observed when the coupled-cavity structures in two dimension triangular dielectric photonic crystals are simulated. The results are in good agreement with experiments.