各向异性介质在三维ADI-FDTD中的应用

该文研究一种减小三维交替方向隐式时域有限差分法(ADI-FDTD)数值色散的新方法。通过在三维空间中合理添加各向异性介质,达到调整相速的目的,从而减小数值色散,使计算结果更加精确。首先对添加各向异性介质后的三维ADI-FDTD迭代公式进行变形,并得到新的数值色散关系,从而求解得到各向异性介质的相对介电常数。以空心波导和具有介质不连续性的波导作为数值算例,分析不同的各向异性介质和添加方法对计算精度的影响,并与传统ADI-FDTD得到的结果和计算资源占用情况进行比较。结果表明通过正确选择各向异性介质和添加方法,可以有效地减小三维ADI-FDTD数值色散。     

Genetic Algorithm in Reduction of Numerical Dispersion of 3-D Alternating-Direction-Implicit Finite-Difference Time-Domain Method

A new method to reduce the numerical dispersion of the 3-D alternating-direction-implicit finite-difference time-domain method is proposed. Firstly, the numerical formulations are modified with the artificial anisotropy, and the new numerical dispersion relation is derived analytically. Moreover, theoretical proof of the unconditional stability is shown. Secondly, the relative permittivity tensor of the artificial anisotropy can be obtained by the adaptive genetic algorithm. In order to demonstrate the accuracy and efficiency of this new method, several examples are simulated. The numerical results and the computational requirements of the proposed method are then compared with those of the conventional method and measured data. In addition, the reduction of the numerical dispersion is investigated as the objective function of the genetic algorithm. It is found that this new method is accurate and efficient by choosing a proper objective function     

Improvement on the numerical dispersion of 2-D ADI-FDTD with artificial anisotropy

In this letter, by introducing artificial anisotropy into computational space, a simple and efficient approach to reduce numerical dispersion of the two-dimensional alternating direction implicit finite-difference time-domain (ADI-FDTD) method is proposed. It is shown that performance of the ADI-FDTD method can be improved significantly for both single frequency simulations and relatively wideband problems. Consequently, the usefulness and effectiveness of the ADI-FDTD method can be notably enhanced.     

一种有效减少ADI-FDTD数值色散的方法

ADI—FDTD算法的数值色散效应较为明显，本文的研究表明一种通过添加各向异性媒质来修正相速误差，从而减少FDTD数值色散的方法，同样适用于ADI-FDTD，且收效更为显著。数值运算结果证明该方法能够简单有效地去除较宽频带范围内的色散。     

A study of the numerical dispersion relation for the 2-D ADI-FDTD method

This letter presents a numerical dispersion relation for the two-dimensional (2-D) finite-difference time-domain method based on the alternating-direction implicit time-marching scheme (2-D ADI-FDTD). The proposed analytical relation for 2-D ADI-FDTD is compared with those relations in the previous works. Through numerical tests, the dispersion equation of this work was shown as correct one for 2-D ADI-FDTD.     

基于非零散度关系的交替方向隐式减缩FDTD算法

该文证明了即使在无源区域,交替方向隐式时域有限差分法(ADI-FDTD)所给出的电磁场量不满足零散度关系,同时推导出了该散度关系的具体表达式。基于该非零散度关系,将不受Courant稳定条件限制的ADI-FDTD法和能节约最多达1/3内存的减缩时域有限差分(R-FDTD)法结合,提出了一种新的交替方向隐式减缩FDTD算法。该算法保留了ADI-FDTD能增大时间步长,缩短计算时间的优点,同时与ADI-FDTD相比节约了最多达1/3(三维)或2/5(二维)的内存。与基于零散度关系的ADI/R-FDTD相比,该算法避免了采用长时间步长计算时的发散现象。应用所提出的ADI/R-FDTD算法计算了二维自由空间波的传播及一维频率选择表面垂直入射的问题,计算结果与ADI-FDTD计算结果完全一致,验证了ADI/R-FDTD的正确性和有效性。     

色散媒质中ADI-FDTD的PML

基于交替方向隐式(ADI)技术的时域有限差分法(FDTD)是一种非条件稳定的计算方法,该方法的时间步长不受Courant稳定条件限制,而是由数值色散误差决定。与传统的FDTD相比, ADI-FDTD增大了时间步长, 从而缩短了总的计算时间。该文采用递归卷积(RC)方法导出了二维情况下色散媒质中ADI-FDTD的完全匹配层(PML)公式。应用推导公式计算了色散土壤中目标的散射,并与色散媒质中FDTD结果对比,在大量减少计算时间的情况下,两者结果符合较好。     

高阶ADI-FDTD算法的数值色散分析

该文给出高阶交替方向隐时域优先差分(ADI-FDTD)算法,即在ADI-FDTD迭代公式的基础上对时间的差分仍然采用二阶中心差分格式,而对空间的差分则采用四阶中心差分格式,并解析地证明了所给出的高阶ADI-FDTD算法仍然满足无条件稳定方程,同时对增长因子相位的分析,得到数值色散关系,最后对其数值色散误差进行了分析,研究表明与普通ADI-FDTD相比,其色散误差较小。     

An Arbitrary-Order LOD-FDTD Method and its Stability and Numerical Dispersion

An arbitrary-order unconditionally stable three-dimensional (3-D) locally-one- dimensional finite-difference time-method (FDTD) (LOD-FDTD) method is proposed. Theoretical proof and numerical verification of the unconditional stability are shown and numerical dispersion is derived analytically. Effects of discretization parameters on the numerical dispersion errors are studied comprehensively. It is found that the second-order LOD-FDTD has the same level of numerical dispersion error as that of the unconditionally stable alternating direction implicit finite-difference time-domain (ADI-FDTD) method and other LOD-FDTD methods but with higher computational efficiency. To reduce the dispersion errors, either a higher-order LOD-FDTD method or a denser grid can be applied, but the choice has to be carefully made in order to achieve best trade-off between the accuracy and computational efficiency. The work presented in this paper lays the foundations and guidelines for practical uses of the LOD method including the potential mixed-order LOD-FDTD methods.      

ADI-FDTD+GRT在波导电路分析中的应用

本文研究时域有限差分法(FDTD)的一种新的时空压缩技术,并应用于波导电路的分析.首先分析了软激励条件下的改进的几何重置技术(GRT),研究了合理选择源面与参考面的放置位置,使GRT不仅减小了吸收边界对计算结果的影响,而且节省了计算空间,还可以精确得到全部散射参量.另外阐述了与交替方向隐式时域有限差分法(ADI-FDTD)相结合,使计算空间和时间同时被压缩,达到节省计算资源的目的.为了衡量ADI-FDTD+GRT算法的计算精度和效率,分析了包含不连续结构的波导作为算例,将其数值计算结果分别与传统FDTD和HFSS作比较,并将端面和参考面不同间距的ADI-FDTD+GRT与传统ADI-FDTD在仿真结果和资源占用方面进行对比,结果表明本文算法是精确和高效的.     

UPML媒质中的包络ADI-FDTD

由于交替方向隐式时域有限差分法(Alternating-Direction Implicit Finite-Difference Time Domain,ADI-FDTD)的数值色散会随着时间步长的增加而增加,文中讨论了单轴各向异性完全匹配层(uniaxial perfectly matched layer,UPML)媒质中包络交替方向隐式时域有限差分法(Envelope ADI-FDTD),推导了二维Envelope ADI-FDTD UPML的迭代公式,并提出一种新的离散方法。与ADI-FDTD UPML相比,改进后的Envelope ADI-FDTD UPML的时间步长可以取得更大,且能有效地修正相速误差,从而减少数值色散,提高计算精度。     

Approximate Crank-Nicolson schemes for the 2-D finite-difference time-domain method for TE/sub z/ waves

Two implicit finite-difference time-domain (FDTD) methods are presented in this paper for a two-dimensional TE/sub z/ wave, which are based on the unconditionally-stable Crank-Nicolson scheme. To treat PEC boundaries efficiently, the methods deal with the electric field components rather than the magnetic field. The "approximate-decoupling method" solves two tridiagonal matrices and computes only one explicit equation for a full update cycle. It has the same numerical dispersion relation as the ADI-FDTD method. The "cycle-sweep method" solves two tridiagonal matrices, and computes two equations explicitly for a full update cycle. It has the same numerical dispersion relation as the previously-reported Crank-Nicolson-Douglas-Gunn algorithm, which solves for the magnetic field. The cycle-sweep method has much smaller numerical anisotropy than the approximate-decoupling method, though the dispersion error is the same along the axes as, and larger along the 45/spl deg/ diagonal than ADI-FDTD. With different formulations, two algorithms for the approximate-decoupling method and four algorithms for the cycle-sweep method are presented. All the six algorithms are strictly nondissipative, unconditionally stable, and are tested by numerical computation in this paper. The numerical dispersion relations are validated by numerical experiments, and very good agreement between the experiments and the theoretical predication is obtained.     

An efficient method to reduce the numerical dispersion in the ADI-FDTD

A new approach to reduce the numerical dispersion in the finite-difference time-domain (FDTD) method with alternating-direction implicit (ADI) is studied. By adding anisotropic parameters into the ADI-FDTD formulas, the error of the numerical phase velocity can be controlled, causing the numerical dispersion to decrease significantly. The numerical stability and dispersion relation are discussed in this paper. Numerical experiments are given to substantiate the proposed method.     

Envelope ADI–FDTD Method and Its Application in Three-Dimensional Nonuniform Meshes

The envelope alternating-direction-implicit finite difference time domain (ADI-FDTD) method in 3-D nonuniform meshes was proposed and studied. The phase velocity error for the envelope ADI-FDTD and ADI-FDTD methods in uniform and nonuniform meshes and different temporal increments were studied. A cavity problem was studied using the envelope ADI-FDTD and ADI-FDTD methods in graded meshes and the conventional FDTD method in a uniform mesh. The simulation results show that the envelope ADI-FDTD performs better than the ADI-FDTD in numerical accuracy     

Theoretical Proof of Unconditional Stability of the 3-D ADI-FDTD Method

In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.     

Dispersion analysis of the three-dimensional complex envelope ADI-FDTD method

In this paper, numerical dispersion properties of the three-dimensional complex envelope (CE) alternate-direction implicit finite-difference time-domain (ADI-FDTD) method are studied. The variations of dispersion errors with propagation direction, ratio of carrier to envelope frequencies, and spatial and temporal steps are presented. It is found that the CE ADI-FDTD scheme have much better accuracy and efficiency over the ADI-FDTD, especially with a higher ratio of carrier to envelope frequencies. Therefore, the CE ADI-FDTD is recommended for use in efficient narrow bandwidth electromagnetic modeling.     

Comments on "An efficient PML implementation for the ADI-FDTD method"

For original paper see Wang and Teixeira (IEEE Microwave Wireless Comp. Lett., vol.13, p.72-4, 2003 February). In this paper, a more precise way to evaluate the actual performance of the perfectly matched layer (PML) used for the alternating direction implicit finite-difference time-domain (ADI-FDTD) method is presented. It is shown that the intrinsic numerical dispersion error of the ADI-FDTD method must be taken into account when the actual performance of the ADI-PML (as well as the ADI-FDTD method) is evaluated. Most importantly, it is demonstrated that the ADI-PMLs implemented with either the traditional manner or the way proposed in have almost the same level of accuracy when the performance of the ADI-PML is correctly evaluated.     

Efficient High-Order Absorbing Boundary Condition for the ADI-FDTD Method

In this letter, an efficient high-order absorbing boundary condition for the alternating-direction-implicit-finite-difference time-domain (ADI-FDTD) method is developed. The proposed high-order absorbing boundary condition (ABC) has very high numerical efficiency because it can keep the tri-diagonal matrix form during the ADI iteration. The absorbing performance of the proposed ABC is analyzed theoretically. To verify the performance of the proposed high-order ABC, a rectangular waveguide problem which has very strong dispersion characteristic is analyzed with the ADI-FDTD method. Numerical results show that the proposed high-order ABC has excellent absorbing performance.     

3-D ADI-FDTD method-unconditionally stable time-domain algorithmfor solving full vector Maxwell's equations

We previously introduced the alternating direction implicit finite-difference time domain (ADI-FDTD) method for a two-dimensional TE wave. We analytically and numerically verified that the algorithm of the method is unconditionally stable and free from the Courant-Friedrich-Levy condition restraint. In this paper, we extend this approach to a full three-dimensional (3-D) wave. Numerical formulations of the 3-D ADI-FDTD method are presented and simulation results are compared to those using the conventional 3-D finite-difference time-domain (FDTD) method. We numerically verify that the 3-D ADI-FDTD method is also unconditionally stable and it is more efficient than the conventional 3-D FDTD method in terms of the central processing unit time if the size of the local minimum cell in the computational domain is much smaller than the other cells and the wavelength     

An improved ADI-FDTD method and its application to photonicsimulations

When the alternating direction implicit-finite difference time domain method (ADI-FDTD) is applied to simulating photonic devices, full efficiency can not be achieved if reasonable accuracy is to be kept, due to numerical errors such as numerical dispersion. A simple modification to ADI-FDTD is proposed by calculating the envelope rather than the fast-varying field, so that errors are minimized. A factor of two-five in speed can usually be gained while retaining the same level of accuracy compared with conventional FDTD. The efficiency and the accuracy of this improved approach is demonstrated on several problems, from simple waveguide structures to complex photonic crystal structures