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

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

色散媒质中ADI-FDTD的PML

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

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.     

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.     

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.     

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

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

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

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

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

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.     

色散媒质中采用Z变换的ADI-FDTD方法

基于Z变换方法将ADI-FDTD推广应用于色散媒质,得到了二维情况下色散媒质中的迭代差分公式,同时给出了一种采用ⅡR滤波器结构来减少存储量的方法.本文方法适合于任意M阶色散媒质,最后,给出了一个算例,数据仿真结果表明,本文的算法与传统色散媒质中的FDTD相比,在计算结果吻合的情况下,存储量相当,计算效率却更高,时间步长仅仅由计算精度来决定.     

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.     

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.     

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.      

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.     

A comprehensive study of numerical anisotropy and dispersion in 3-DTLM meshes

This paper presents a comprehensive analysis of the numerical anisotropy and dispersion of 3-D TLM meshes constructed using several generalized symmetrical condensed TLM nodes. The dispersion analysis is performed in isotropic lossless, isotropic lossy and anisotropic lossless media and yields a comparison of the simulation accuracy for the different TLM nodes. The effect of mesh grading on the numerical dispersion is also determined. The results compare meshes constructed with Johns' symmetrical condensed node (SCN), two hybrid symmetrical condensed nodes (HSCN) and two frequency domain symmetrical condensed nodes (FDSCN). It has been found that under certain circumstances, the time domain nodes may introduce numerical anisotropy when modelling isotropic media     

等离子体中的PLCDRC-ADI-FDTD方法

采用分段线性电流密度递归卷积(Piecewise Linear Current Density Recursive Convolution)方法将交替方向隐式时域有限差分方法(ADI-FDTD)推广应用于色散介质—等离子体中,得到了二维情况下等离子体中的迭代差分公式,为了验证该方法的有效性和可靠性,计算了等离子体涂敷导体圆柱的RCS和非均匀等离子体平板的反射系数,数据仿真结果表明,此算法与传统的FDTD相比,在计算结果吻合的情况下,存储量相当,计算效率更高,时间步长仅仅由计算精度来决定.     

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.     

Analysis of the numerical One-Step method for the study applied on bio electromagnetics

The development of wireless technologies arises important questions about the effects of the wave propagation in the human body. To study accurately these effects, we have to use rigorous numerical methods. In this paper, we present and analyze the One-Step time domain method. This method, which was proposed by De Raedt [Phys Rev E 67(056706):1–12, 2003] for lossless media, is known to be unconditionally stable and so it can be used for applications for which the Courant–Friedrich–Levy (CFL) stability condition can be a limiting factor, e.g., for bioelectromagnetic applications. The numerical dispersion and the insertion of lossy media in the One-Step method are evaluated. The perfectly matched layer (PML) absorbing conditions are also introduced in our study.     

Conformal method to eliminate the ADI-FDTD staircasing errors

The alternating-direction-implicit finite-difference time-domain (ADI-FDTD) method is an unconditionally stable method and allows the time step to be increased beyond the Courant-Friedrich-Levy (CFL) stability condition. This method is potentially very useful for modeling electrically small but complex features often encountered in applications. As the regular FDTD method, however, the spatial discretization in the ADI-FDTD method is only first-order accurate for discontinuous media; several researchers have shown that the errors can be very high when the regular ADI-FDTD method is applied to such discontinuous media. On the other hand, the conformal FDTD method has recently emerged as an efficient FDTD method with higher order accuracy. In this work, a second-order accurate ADI-FDTD method using the conformal approximation of spatial derivatives is proposed. This new scheme, called the ADI-CFDTD method, retains the second-order accuracy in both temporal and spatial discretizations even for discontinuous media with metallic structures, and is unconditionally stable. 2D and 3D examples demonstrate the efficacy of this method and its application in EMC problems.     

色散媒质中ADI-FDTD

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