GENETIC ALGORITHM IN REDUCTION OF NUMERICAL DISPERSION OF 3-D ADI-FDTD METHOD

A new method to reduce the numerical dispersion of the three-dimensional Alternating Direction Implicit Finite-Difference Time-Domain （3-D ADI-FDTD） method is proposed. Firstly, the numerical formulations of the 3-D ADI-FDTD method are modified with the artificial anisotropy, and the new numerical dispersion relation is derived. Secondly, the relative permittivity tensor of the artificial anisotropy can be obtained by the Adaptive Genetic Algorithm （AGA）. In order to demonstrate the accuracy and efficiency of this new method, a monopole antenna is simulated as an example. And the numerical results and the computational requirements of the proposed method are compared with those of the conventional ADI-FDTD method and the measured data. In addition the reduction of the numerical dispersion is investigated as the objective function of the AGA. It is found that this new method is accurate and efficient by choosing proper objective function.     

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

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

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

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

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.     

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.     

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.     

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

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

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     

Investigation on the characteristics of the envelope FDTD based on the alternating direction implicit scheme

This letter presents numerical characteristics of recently developed the envelope FDTD based on the alternating direction implicit scheme (envelope ADI-FDTD). Through numerical simulations, it is shown that the envelope ADI-FDTD is unconditionally stable and we can get better dispersion accuracy than the traditional ADI-FDTD by analyzing the envelope of the signal. This fact gives the opportunity to extend the temporal step size to the Nyquist limit in certain cases. Numerical results show that the envelope ADI-FDTD can be used as an efficient electromagnetic analysis tool especially in the single frequency or band limited systems.     

Stability and Dispersion Analysis for ADI-FDTD Method in Lossy Media

The stability and dispersion analysis for the alternating-direction-implicit finite-difference time-domain (ADI- FDTD) method in lossy media is presented. Although the stability and numerical dispersion have been analyzed for the ADI-FDTD method, most of the analysis is dedicated to the cases of lossless media. Here, the stability and dispersion analysis is performed for the method in lossy media. The stability analysis theoretically proves the unconditional stability of the ADI-FDTD method in lossy media. Meanwhile, the dispersion analysis reveals the numerical loss and dispersion characteristics of this method. This will be meaningful for the evaluation and further development of the ADI-FDTD method in lossy media     

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的时间步长可以取得更大,且能有效地修正相速误差,从而减少数值色散,提高计算精度。     

基于非零散度关系的交替方向隐式减缩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的正确性和有效性。     

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.      

Alternating-direction implicit (ADI) formulation of the finite-difference time-domain (FDTD) method: algorithm and material dispersion implementation

The alternating-direction implicit finite-difference time-domain (ADI-FDTD) technique is an unconditionally stable time-domain numerical scheme, allowing the /spl Delta/t time step to be increased beyond the Courant-Friedrichs-Lewy limit. Execution time of a simulation is inversely proportional to /spl Delta/t, and as such, increasing /spl Delta/t results in a decrease of execution time. The ADI-FDTD technique greatly increases the utility of the FDTD technique for electromagnetic compatibility problems. Once the basics of the ADI-FDTD technique are presented and the differences of the relative accuracy of ADI-FDTD and standard FDTD are discussed, the problems that benefit greatly from ADI-FDTD are described. A discussion is given on the true time savings of applying the ADI-FDTD technique. The feasibility of using higher order spatial and temporal techniques with ADI-FDTD is presented. The incorporation of frequency dependent material properties (material dispersion) into ADI-FDTD is also presented. The material dispersion scheme is implemented into a one-dimensional and three-dimensional problem space. The scheme is shown to be both accurate and unconditionally stable.     

On the dispersion relation of ADI-FDTD

In this letter, we analyze the alternating direction implicit finite-difference time-domain (ADI-FDTD) dispersion relation and find the numerical plane-wave relationship between the magnetic and electric fields, showing that the scheme is not divergence-free. We also show that negative-group-velocity modes with positive phase velocities may appear for high Courant numbers. A parallel comparison is made with the behavior of the classical Yee FDTD and the Crank-Nicolson schemes.     

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.     

Higher-order alternative direction implicit FDTD method

A second-order in time and fourth-order in space alternative direction implicit (ADI) finite difference time domain (FD-TD) method has been developed. The dispersion relation is derived and compared with the ADI-FDTD method. Numerical results demonstrated that the higher-order ADI-FDTD has a better accuracy compared to the ADI-FDTD method.     

ADI-FDTD Method With Fourth Order Accuracy in Time

This letter presents an unconditionally stable alternating direction implicit finite-difference time-domain (ADI-FDTD) method with fourth order accuracy in time. Analytical proof of unconditional stability and detailed analysis of numerical dispersion are presented. Compared to second order ADI-FDTD and six-steps SS-FDTD, the fourth order ADI-FDTD generally achieves lower phase velocity error for sufficiently fine mesh. Using finer mesh gridding also reduces the phase velocity error floor, which dictates the accuracy limit due to spatial discretization errors when the time step size is reduced further.     

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.