The 3-D multidomain pseudospectral time-domain algorithm for inhomogeneous conductive media

A three-dimensional (3-D) multidomain pseudospectral time-domain (PSTD) method with a strongly well-posed perfectly matched layer (PML) is developed as an accurate and flexible tool for the simulation of electromagnetic wave propagation and scattering in inhomogeneous and conductive media. This approach allows for an accurate treatment of curved geometries by multidomain decomposition and curvilinear coordinate transformation. Numerical experiments show the results agree excellently with analytical solutions and results of other well-known algorithms, and demonstrate a remarkable improvement in accuracy and efficiency over the FDTD method. The 3-D multidomain PSTD algorithm is then applied to calculate radar cross sections (RCS).     

Two dimensional multidomain pseudospectral time-domain algorithm based on alternating-direction implicit method

In order to deal with the stability problem restricted by the finite-difference time-domain (FDTD) and conventional pseudospectral time domain (PSTD), the multidomain PSTD algorithm based on alternating-direction implicit (ADI) technique is proposed in this paper. This algorithm improves the stability and efficiency of conventional PSTD, while it maintains the accuracy and flexibility of conventional PSTD for an accurate treatment of arbitrarily curved objects. A compact matrix form is derived to effectively describe two-dimensional ADI multidomain pseudospectral time domain (ADI-MPSTD) algorithm. Numerical results show an excellent agreement with analytical solutions as well as results obtained by the FDTD algorithm, and fully demonstrate a remarkable improvement in stability and efficiency.     

MPSTD算法中子域分界面匹配条件的比较研究

针对多区域伪谱时域算法(MPSTD)中子域分界面两侧为不同介质的情形,提出了特征变量-物理边界(CV-PB)子域分界面匹配条件,然后在总结现有子域分界面匹配条件的基础上,分别比较了子域分界面两侧为相同和不同介质时各种子域分界面匹配条件的稳定性.数值仿真表明,在子域分界面两侧为相同介质材料时采用特征变量法匹配条件,而在子域分界面两侧为不同介质材料时采用CV-PB匹配条件,可使得MPSTD算法稳定性最好.     

The unconditionally stable pseudospectral time-domain (PSTD) method

This letter presents a new time-domain method for Maxwell's equations, in which the unconditionally stable techniques, the alternating direction implicit (ADI) and the split-step (SS) schemes, are developed for the pseudospectral time-domain (PSTD) algorithm to maintain stability while achieving higher accuracy and efficiency over the FDTD method. The multidomain strategy is employed to allow for a flexible treatment of internal inhomogeneities. Numerical results demonstrate the unconditional stability and the second-order accuracy for both ADI- and SS-PSTD algorithms.     

A pseudospectral method for time-domain computation ofelectromagnetic scattering by bodies of revolution

We present a multidomain pseudospectral method for the accurate and efficient time-domain computation of scattering by body-of-revolution (BOR) perfectly electrically conducting objects. In the BOR formulation of the Maxwell equations, the azimuthal dependence of the fields is accounted for analytically through a Fourier series. The numerical scheme in the (r,z) plane is developed in general curvilinear coordinates and the method of characteristics is applied for patching field values in the individual subdomains to obtain the global solution. A modified matched-layer method is used for terminating the computational domain and special attention is given to proper treatment of the coordinate singularity in the scattered field formulation and correct time-domain boundary conditions along edges. Numerical results for monochromatic plane wave scattering by smooth and nonsmooth axis-symmetric objects, including spheres, cone-spheres, and finite cylinders, is compared with results from the literature, illustrating the accuracy and computational efficiency associated with the use of properly constructed spectral methods. To emphasize the versatility of the presented framework, we compute plane wave scattering by a missile and find satisfactory agreement with method-of-moment (MoM) computations     

A fast higher-order time-domain finite element-boundary integral method for 3-D electromagnetic scattering analysis

A novel hybrid time-domain finite element-boundary integral method for analyzing three-dimensional (3-D) electromagnetic scattering phenomena is presented. The method couples finite element and boundary integral field representations in a way that results in a sparse system matrix and solutions that are devoid of spurious modes. To accurately represent the unknown fields, the scheme employs higher-order vector basis functions defined on curvilinear tetrahedral elements. To handle problems involving electrically large objects, the multilevel plane-wave time-domain algorithm is used to accelerate the evaluation of the boundary integrals. Numerical results demonstrate the accuracy and versatility of the proposed scheme.     

A fast-decoupled method for time-domain simulation of powerconverters

Fast long-term dynamic simulation of power converters is discussed. Because of switching actions, the converter can operate in different configurations and, therefore, it can be represented by a state-space time-varying linear system. This paper presents a new method for calculating switching time from an accurate and explicit form of the switching condition, hence the step-by-step search of the switching instant can be greatly simplified. In the proposed algorithm, a converter is partitioned into subsystems, which are simulated in a decoupled manner and coordinated through data exchange     

Analysis of modes in a freestanding microsquare resonator by 3-D finite-difference time-domain simulation

Modes in a microsquare resonator slab with strong vertical waveguide consisting of air/semiconductor/air are analyzed by three-dimensional (3-D) finite-difference time-domain simulation, and compared with that of two-dimensional (2-D) simulation under effective index approximation. Mode frequencies and field distributions inside the resonator obtained by the 3-D simulation are in good agreement with those of the 2-D approximation. However, field distributions at the boundary of the resonator obtained by 3-D simulation are different from that of the 2-D simulation, especially the vertical field distribution near the boundary is great different from that of the slab waveguide, which is used in the effective index approximation. Furthermore the quality factors obtained by 3-D simulation are much larger than that by 2-D simulation for the square resonator slab with the strong vertical waveguide.     

Piecewise-exponential approximation for fast time-domain simulation of 2-D cellular neural networks

Cellular neural networks (CNNs) were introduced as promising image processing systems. However, since analytical design techniques are rarely available, extensive simulation is the main practical tool for developing significant applications. This paper presents a new algorithm for fast simulation of large-scale CNNs. It is based on the discretization of the sigmoid generating the output from the state of each cell. This discretization leads to a piecewise exponential approximation of the time-domain solution. Computation is only required when the output of a cell jumps to a different discrete level and involves only this cell and its neighbors. The algorithm is spatially adaptive since the computational effort is concentrated on the most rapidly evolving portions of the array.     

A method for fast time-domain simulation of networks with switches

A method for fast time-domain simulation of piecewise-linear networks with switches is described in this paper. The method is based on a discrete-time switch model that consists of a constant conductance in parallel with a current source. In each simulation step, the value of the current source is updated as a function of known network signals. The function takes one of two forms, depending on the state (on or off) of the switch. Since the system matrix is constant, regardless of the states of the switches, simulation time is essentially the same as for a linear, time-invariant network of the same complexity. The paper discusses selection of the model and simulation parameters. The simulation algorithm is described and an example is included. It is shown that the method is not only efficient but also quite general and void of convergence problems. Its primary application is for long-term transient simulation of power electronic systems such as switching power converters     

求解三维麦克斯韦方程的时域无网格算法研究

发展了用于求解三维麦克斯韦方程的时域无网格算法.算法基于生成的无网格点云,通过泰勒级数展开结合加权最小二乘逼近计算点云中心点上的空间导数,并构造近似黎曼解处理空间离散涉及的通量运算;空间离散后的半离散方程则采用四步Runge-Kutta格式推进求解.结合求解三维麦克斯韦方程,给出了时域无网格算法的具体实施过程,并基于发展的算法,成功地模拟出金属球、立方体及进气道模型等三维散射目标的电磁散射场,获得的雷达散射截面能与理论解、矩量法或精确控制法等结果吻合.     

Discontinuous galerkin time-domain method for GPR simulation in dispersive media

This paper presents a newly developed high-order discontinuous Galerkin time-domain (DGTD) method for solving Maxwell's equations in linear dispersive media with UPML boundary treatment. A unified formulation is derived for linear dispersive media of Debye type and the artificial material in the UPML regions with the help of auxiliary differential equations. The DGTD employs finite-element-type meshes, and uses piecewise high-order polynomials for spatial discretization and Runge-Kutta method for time integrations. Arbitrary high-order accuracy can be obtained for scattering of various objects in dispersive media. After validating the numerical convergence of the DGTD method together with the second-order Yee's scheme, we apply this new method to the ground-penetrating radar for the detection of buried objects in a lossy half space.     

Efficient compact 2-D time-domain method with weighted Laguerre polynomials

An efficient time-domain method based on a compact two-dimensional (2-D) finite-difference time-domain (FDTD) method combined with weighted Laguerre polynomials has been proposed to analyze the propagation properties of uniform transmission lines. Starting from Maxwell's differential equations corresponding to the compact 2-D FDTD method, we use the orthonormality of weighted Laguerre polynomials and Galerkin's testing procedure to eliminate the time variable. Thus, an implicit relation, which results in a marching-on-in-degree scheme, can be obtained. To verify the accuracy and efficiency of the hybrid method, we compare the results with those from the conventional compact 2-D FDTD and compact 2-D alternating-direction-implicit (ADI) FDTD methods. The hybrid method improves the computational efficiency notably, especially for complex problems with fine structure details that are restricted by stability constrains in the FDTD method.     

计算电大尺寸建筑物内电波场强的PSTD方法

本文采用一种新的时域数值方法－伪谱时域（ＰＳＴＤ）法来计算电大尺寸建筑物内电波场强。提出了由初始条件技术和一维ＰＳＴＤ方程对入射平面波脉冲进行模拟,并利用纯散射场法和线性插值把平面波引入求解问题空间,有效地解决了ＰＳＴＤ方法中入射波设置问题。数值结果表明这种新方法用于模拟电大尺寸建筑物内电波场强的精度和有效性。     

The transfinite-element time-domain method

This paper presents an efficient time-domain method for computing the propagation of electromagnetic waves in microwave structures. The procedure uses high-order vector bases to achieve high-order accuracy in space, Newmark's method to provide unconditional stability in time, and the transfinite-element method to truncate the waveguide ports. The resulting system matrix is real, symmetric, positive-definite, and can be solved by using the highly efficient multilevel preconditioned conjugate gradient algorithm. Since the method allows large time steps and nonuniform grids, the computational complexity for problems with irregular geometries is superior to that of the finite-difference time-domain method.     

Parametric oscillations in photonic crystal slabs 3-D time-domain analysis

Slab waveguides with a two-dimensional periodic distribution of the refractive index are proposed and investigated as optical microcavities for efficient parametric oscillations through four-wave mixing in isotropic materials. We carry out a case study based on the complete solution of three-dimensional Maxwell equations, including material dispersion and cubic nonlinear response.     

一种基于宽带MIMO雷达时域成像的阵列布阵模型

多输入多输出(MIMO)霄达足近几年发展起来的一种新慨念雷达体制.为了进一步降低宽带MIMO雷达的阵列规模和硬件复杂度,使MIMO雷达成像技术实用化成为可能,结合改进型后向投影(BP)成像算法,提出了一种非均匀线阵成像系统阵列布阵模型.提出的阵列布阵模型通过成像系统方位向分辨率婴求和目标方位向成像场景大小确定阵元间距,使得系统的阵列规模和硬件复杂度显著降低.同时,时域成像算法的使用避免了频域算法存在的采样定理限制,有效增强了阵列设计的灵活性,并使得成像阵列规模和硬件复杂度得到进一步降低.最后仿真实验验证了该成像系统布阵模型的正确性和有效性.     

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.     

A semivectorial finite-difference time-domain method (optical guided structure simulation)

A semivectorial finite-difference time-domain method (FDTD) that solves the vector wave equations for the transverse electric fields is presented and validated. By taking into consideration the boundary conditions for the transverse, electric fields in the finite-difference scheme, the polarization effect of the electromagnetic waves can be modeled. In comparison with the full vector FDTD, the present approach requires less memory and is more computational efficient. The method is validated by a comparison with the exact analytical solutions as well as the full vector FDTD results and is shown to be very accurate.<>