An Effective Multigrid Preconditioned CG Algorithm for Millimeter Wave Scattering by an Infinite Plane Metallic Grating

An effective multigrid based preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using multigrid technique we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the initial guess for the conjugate gradient (CG) method. Numerical results show that our developed method can reach five times improvement of computational complexity.     

Analysis of Millimeter Wave Scattering by an Electrically Large Metallic Grating Using Wavelet-Based Algebraic Multigrid Preconditioned CG Method

An effective wavelet based multigrid preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using wavelet transformation we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the new improvement for the conjugate gradient (CG) method. Numerical results show that our developed wavelet based multigrid preconditioned CG method can reach large improvement of computational complexity. Due to the automaticity of wavelet transformation, this method is potential to be a block box solver without physical background.     

Millimeter Wave Scattering of 2-D Frequency Selective Surface by Efficient Method of Lines with Preconditioned CG Technique

In this paper, both banded and symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) techniques are combined with method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. The electromagnetic wave scattering of 2-D frequency-selective surface is used as the examples to describe its implementation, whose analysis usually needs fast algorithms because of electrically large dimension. For arbitrary incident wave, helmholtz equation and boundary condition are used to calculate the impedance matrix and then to obtain reduced current-voltage linear matrix equation in spatial domain. Both banded and effective symmetric successive overrelaxation preconditioned conjugate gradient iterative method are chosen to solve this matrix equation. Our numerical results show that PCG methods can converge to accurate solution in much fewer iteration steps for analysis of the electromagnetic wave scattering from 2-D frequency-selective surface.     

Application of Preconditioned GMRES Algorithm for Analysis of Millimeter Wave Ferrite Sphere Circulators

The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. In order to speed up the convergence of GMRES, the symmetric successive overrelaxation (SSOR) preconditioning scheme is applied for the analysis of millimeter wave ferrite circulator. Consequently, this preconditioned GMRES (PGMRES) approach can reach convergence ten times faster than GMRES. The reflection and insertion losses of millimeter wave waveguide circulator are compared with those obtained from literature.     

Application of Preconditioned GMRES Algorithm for Analysis of Millimeter Wave Ferrite Sphere Circulators

The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. In order to speed up the convergence of GMRES, the symmetric successive overrelaxation (SSOR) preconditioning scheme is applied for the analysis of millimeter wave ferrite circulator. Consequently, this preconditioned GMRES (PGMRES) approach can reach convergence 19 times faster than GMRES. The isolation and insertion losses of millimeter wave waveguide circulator are compared with those obtained from literature.     

A Fast Finite Element Analysis of Millimeter Wave Waveguide Filters

The multifrontal method is applied for solving a large system of linear equations resulting from the use of edge-based finite element method (FEM). The finite element method combined with the perfectly matched layers (PML) is given for simulation of microwave filters and the algorithm of multifrontal method is described in detail. The reflection and insertion losses of millimeter wave filter are analyzed as the examples and the obtained results are compared with those obtained from literature. In order to demonstrate the efficiency of the multifrontal method, the computational time is compared with that of both symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) and conjugate gradient methods (CG) for the thick-iris waveguide bandpass filer.     

Application of SSOR Preconditioning Technique to Method of Lines for Millimeter Wave Scattering

In this paper, symmetric successive overrelaxation (SSOR) preconditioned CG technique are introduced into method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. Millimeter wave scattering by an infinite plane metallic grating is used as the examples to describe its implementation, whose analysis usually needs fast algorithms because of electrically large dimension. For arbitrary incident wave, Helmholz equation and boundary condition are used to calculate the impedance matrix and then to obtain reduced current-voltage linear matrix equation in spatial domain. An effective symmetric successive overrelaxation preconditioned conjugate gradient iterative method, SSOR-PCG, is chosen to solve this matrix equation. With SSOR as the preconditioner as well as its efficient implementing in CG algorithm, PCG method can converge to accurate solution in much fewer iteration steps.     

Simulation of Millimeter Wave Circulators by FEM with Unifrontal/Multifrontal Technique

The combination of unifrontal and multifrontal techniques is applied for solving a large, sparse, and unsymmetric system of linear equations resulting from the use of edge-based finite element method (FEM) for millimeter wave ferrite devices. The formulation of finite element method combined with the perfectly matched layers (PML) is given and the combined algorithm of unifrontal and multifrontal methods is described. The lower data movement of unifrontal method and the lower fill-in of normal multifrontal methods are combined for an algorithm with their features. The performance of typical waveguide junction circulators is analyzed and the calculated results are compared with those obtained from the literature.     

Analysis of Millimeter Wave Scattering by the Infinite Plane Metallic Grating Using PCG-FFT Technique

In this paper, both fast Fourier transformation (FFT) and preconditioned CG technique are introduced into method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. Electromagnetic wave scattering by an infinite plane metallic grating is used as the examples to describe its implementation. For arbitrary incident wave, Helmholz equation and boundary condition are first transformed into new ones so that the impedance matrix elements are calculated by FFT technique. As a result, this Topelitz impedance matrix only requires O(N) memory storage for the conjugate gradient FFT method to solve the current distribution with the computational complexity O(N log N) . Our numerical results show that circulate matrix preconditioner can speed up CG-FFT method to converge in much smaller CPU time than the banded matrix preconditioner.     

Fast Stimulation of Millimeter Wave Waveguide Filters by Boundary Marching Method

The boundary marching method is applied for the analysis of the waveguide bandpass filers. The vector finite element method with the perfectly matched layers (PML) as the matching load is given for parameter extraction of millimeter wave filters. The implementation of the fast boundary marching method is described in detail. The reflection and insertion losses of millimeter wave filters are obtained and compared with those obtained from literature.     

Efficient Analysis of Millimeter Wave Ferrite Circulators by the GMRES Iterative Algorithm

The finite element method (FEM) combined with the perfectly matched layers (PML) is given for simulation of waveguide ferrite circulators. The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. The formulation of FEM and the algorithm of GMRES method are described in detail. The reflection and insertion losses of millimeter wave waveguide circulator are analyzed and the results are compared with those obtained from literature.     

Efficient Analysis of Millimeter Wave Circuit Using Numerical Calibration Technique

In this paper, A numerical technique, called short-open calibration (SOC), in conjunction with edge-based finite element method (FEM) is employed to analyze millimeter wave circuit that can be segmented into two distinct section: static model of feedlines and dynamic model of circuit discontinuity. The derivation of reflection coefficient of 3D discontinuities is arranged in two steps. In the first step, this SOC technique is incorporated into the FEM for mesh truncation of computaional domain. In this way, much faster convergence is achieved for large-sparse linear matrix equations from FEM by this termination than by perfectly matching layers (PML). The field distribution of the dominated mode in uniform feedlines and entire circuit is obtained individually by exciting a pair of even and odd impressed voltages along the struture. In step two, Scattering parameters based on the voltages and current defintion is calculated by integral of electric and magnetic fields. Numerical solutions for a class of planar circuit discontuities are very well compared with those published in the available literatures.     

Discontinuity of NRD-Guide and Its Application in Millimeter Wave Circuits

In this paper, the performances of two typical discontinuous structures of the NRD-guide, the traverse metallic rod to disturb an incident wave and the traverse probe to excite required modes, are analyzed and computed by utilizing the theory of modes orthogonality of the NRD-guide. The moment method is used to evaluate the current distribution on the rod and probe and the affection of some higher-order modes is taken into account simultaneously. On these bases, the multitone exciting steady states of two kinds of important circuits that employ the above discontinuous structures, the NRD-guide crossbar balanced mixer and the LSM ₁₁ ⁰ -mode oscillator, are simulated by the harmonic balance technique based on the multi-dimensional discrete Fourier transformation (MDFT-HB method).     

Coupled FEM/BEM for control of noise radiation and sound transmission using piezoelectric shunt damping

In this paper, we present a coupled finite element/boundary element method (FEM/BEM) for control of noise radiation and sound transmission of vibrating structure by passive piezoelectric techniques. The system consists of an elastic structure (with surface mounted piezoelectric patches) coupled to external/internal acoustic domains. The passive shunt damping strategy is employed for vibration attenuation in the low frequency range. The originality of the present paper lies in evaluating the classically used FEM/BEM methods for structural–acoustics problems when taking account smart systems at the fluid–structure interfaces.     

Parameters Extraction of Millimeter Wave Circuits by Hybrid Vector Finite Element and Moment Method Combined with SOC Technique

This paper presents a hybrid method, which couples the vector finite element method (FEM) and method of moment (MOM) for analyzing the field and current distribution of the millimeter wave circuits. The FEM is applied to handle the interior region of dielectric bodies and MOM is used to solve surface integral equations. Then, These integral expressions are coupled into the FEM equations through the continuity of the tangential fields across the connection boundaries. Simultaneously, the short-open calibration (SOC) technique is used for predicting accurately the scattering parameters of the circuits. Numerical results are well compared with those published in the previous literatures.     

Vector-FEM Analysis of Millimeter Wave Circulator with SOC Truncation Technique

In this paper, the application of the edge-based vector finite element method combined with the short-open calibration (SOC) technique to waveguide junction circulators was presented. By applying permittivity and permeability tensors, the functional formula for tetrahedral vector elements are analytically derived in terms of the electric field strength. And the SOC technique is directly accommodated in the FEM algorithm and used to truncate the computational domain. The SOC technique removes or separates unwanted parasitics brought by the approximation of the impressed voltage source and also the problem of resulting inconsistency between different simulations. Truncation of computational domain by SOC technique makes the iterative solvers for large-sparse linear matrix equations from FEM converge much faster than by perfectly matching layers (PML). The analysis of three 2-port networks is sufficient to form the admittance matrix of the corresponding three-port waveguide junction circulators. To validate the theory, typical three-dimensional numerical results of isotropic dielectric-filled discontinuity are first presented and compared with other hybrid numerical techniques. Further, the ferrite filled waveguide junction circulators are analysed and compared with the measured and available publications. The comparison shows that good agreement is obtained.     

激光激发声表面波在缺陷板材中散射过程的有限元分析

利用有限元法模拟了金属板材中激光激发的声表面波经过缺陷位置时发生散射的瞬态过程,采用线状激光源作为超声导波的激发源.针对三种不同深度的表面缺陷以及三种亚表面缺陷的模型进行了对比计算,结果显示缺陷的深度及位置对声表面波的时域特征存在显著的影响.表面缺陷深度越深将产生较大幅度的表面反射回波,亚表面缺陷的影响将取决于缺陷上顶面距离板材上表面的距离.因此,数值模拟结果表明通过分析激光产生的表面波形可以判定近表面缺陷的尺寸和所处的位置.     

全矢量有限元模型及其在光波导中的应用

为了研究光波导和光子晶体光纤的模式特性和传输特性,从矢量波动方程出发,推导出了各向异性介质中场微分方程复数泛函表达式,利用棱边/节点混合元离散了该泛函,加入了各向异性介质匹配层边界条件,得到关于传播常量的广义特征值方程.以矩形波导为例,对各向异性介质匹配层边界条件的吸收特性进行了研究,得到了基模以及几个高阶模的场分布、色散曲线和损耗曲线.结果表明该方法可靠有效.对正六边形晶格光子晶体光纤进行了分析.数据表明：光纤有效折射率随空气孔直径或波长的增大而减小,但与空气孔圈数无关;光纤限制损耗(confinement loss)随波长增大近似成指数增大,而增加空气孔直径或者空气孔圈数则可使之显著降低.     

Some Weighted Averaging Methods for Gradient Recovery

We propose some new weighted averaging methods for gradient recovery, and present analytical and numerical investigation on the performance of these weighted averaging methods. It is shown analytically that the harmonic averaging yields a superconvergent gradient for any mesh in one-dimension and the rectangular mesh in two-dimension. Numerical results indicate that these new weighted averaging methods are better recovered gradient approaches than the simple averaging and geometry averaging methods under triangular mesh.     

太赫兹波光子晶体光纤传输特性分析

利用矢量有限元法分析了太赫兹波光子晶体光纤单模截止频率和波导色散随光纤结构的变化特性．结果表明，太赫兹波光子晶体光纤的单模截止频率随着光纤空气孔占空比的变大而降低，零波导色散点频率随着空气孔占空比变大而增加，约束损耗随着空气孔的圈数增加而降低．