利用高阶矢量基函数求解时域磁场积分方程

本文利用一种新的高阶矢量基函数求解了三维时域磁场积分方程,该基函数定义在一个曲边三角形贴片上并用拉格朗日插值多项式来表示每一个贴片内的未知电流密度.该基函数的实质就是将拉格朗日插值多项式的插值点选为高斯积分结点,极大地简化和加快了时域积分方程矩量法的繁琐的时间和空间积分运算;另外,该基函数不要求网格为规范网格,给复杂目标的网格剖分带来很大方便.在空间上利用点匹配方法求解了时域磁场积分方程,数值计算结果表明了该方法求解时域积分方程的精确性和高效性.     

介质散射的双线性高阶叠层基函数矩量法分析

提出了一种基于双线性思想的高阶基函数降低矩量法分析介质目标电磁散射问题时的未知量.论文给出了双线性高阶叠层基函数的构造过程, 并将其应用于介质的面积分方程分析中.数值算例比较了不同阶数时所需未知量以及计算精度, 表明该高阶叠层基函数在满足相同积分方程的计算精度时能比低阶基函数节省未知量.     

三维磁各向异性目标电磁散射的MOM-CGM-FFT方法

本文给出了研究任意形状三维磁各向异性目标电磁散射问题的一种混合计算方法。该方法以磁场作为未知函数建立频域体积分 微分方程 ,使用脉冲基函数和点匹配函数的矩量法 (MOM)将之转化为线性代数方程组。在求解过程中应用共轭梯度法 (CGM)和快速傅立叶变换 (FFT)相结合的方法降低所需计算机内存和CPU时间。对均匀磁各向同性球的数值计算结果与相应Mie理论结果吻合较好。在此基础上给出了磁各向异性球和立方体的雷达散射截面。计算结果表明该计算方法兼容性强 ,是一种求解三维电磁散射问题的有效途径。     

粗糙海面上三维金属目标的电磁散射特性分析

该文采用矩量法(MoM)计算粗糙海面上三维金属目标的电磁散射特性。计算了位于半空间媒质中的电偶极子和磁偶极子的矢量位并矢格林函数和标量位格林函数,并将其应用于矩量法中。把海水视为下半空间媒质,粗糙海面为位于上半空间中的介质表面。通过建立介质和金属混合目标的积分方程,并采用迭代方法求解矩阵方程以得到该模型的散射特性,数值结果验证了本文方法的有效性。     

E-P-AMCBFM分析介质与PEC混合体的电磁散射

使用基于表面积分方程的矩量法来分析介质与理想导体混合体的电磁散射是计算电磁学的一大热点。对理想导体目标体表面建立电场积分方程,在介质目标体表面建立PMCHW方程组,与基于矩阵分块技术的自适应修正特征基函数法结合,对介质涂敷理想导体目标体的电磁散射进行分析,将其称之为EFIE-PMCHW-AMCBFM(E-P-AMCBFM)。并讨论不同参数如基函数阶数,矩阵块间重叠区域等对计算效率的影响,数值结果表明E-P-AMCBFM对于处理介质-理想导体混合体的电磁散射问题具有较高的精度和效率。     

导体目标瞬态电磁散射的时空矩量法研究

时域积分方程的矩量法是求解瞬态电磁散射的方法之一。研究了基于加权Lagurre函数和RWG基函数分别作为时间、空间基函数的时空矩量法,给出了时域磁场积分方程时空矩量法的全部计算公式,编制了相应串行和并行计算程序。计算结果表明：该方法具有很好的时域稳定特性,为宽带电磁散射分析提供了可能,同时也指出了其应用的局限性,为改进其方法提供了参考。     

并行高阶矩量法分析舰队RCS和其它电磁特性

随着计算机多核技术和计算机群技术的发展和MPI并行技术的应用,使得并行高阶矩量法可以解决复杂大型目标的电磁特性。文中介绍了并行高阶矩量法的基本理论：高阶基函数、并行矩阵填充和并行矩阵求解,并在这些理论的基础上应用其分析舰队的RCS和其它电磁特性来说明这种数值算法的可行性和快速性。     

高阶MoM及其与PO的混合法计算电磁散射问题

本文提出了一种新的以高阶矩量法(MoM)与物理光学法相结合的混合法(MoM-PO).该方法采用曲面参数化的离散方法,保证了建模的精确性.计算过程中将散射表面灵活划分为MoM区和PO区,在各自区域可以灵活确定离散单元的大小和密度.MoM区域的高阶矩量法,采用基于Lagrange插值的高阶矢量基函数,结合点匹配技术,比传统的高阶法简单,易于实现.计算结果表明,本文的高阶矩量法及其与物理光学法结合的混合方法能准确有效的计算目标的电磁散射特性.     

金属-各向异性介质体组合目标频域体表积分方程矩量法

利用体表积分方程矩量法求解了具有任意的介电常数张量和磁导率张量的各向异性介质与金属的组合目标的电磁散射问题.给出了基于RWG面基函数和SWG体基函数的体表积分方程阻抗矩阵元素表达式并详细推导了阻抗矩阵元素所涉及的各种积分运算的计算方法;通过数值计算实例与解析解或其它数值方法的详细对比分析,证明了计算公式的正确性.     

埋地三维复杂介质目标电磁散射特性分析

采用混合位积分方程（MPIE）和基于四面体元基函数的矩量法分析计算了埋地三维介质目标电磁散射问题,利用二级离散复镜像（DCIM）和广义函数束（GPOF）相结合的方法求解近场Sommerfeld积分,很好地解决了多层媒质中复杂目标电磁散射计算中的棘手问题,其方法简练、精确、高效,数值分析结果与有关文献吻合很好,证实了该方法的正确性和通用性。此外,该文还通过计算比较了不同观察点、不同目标埋地深度及介质参数的电磁散射特性。     

Integral equation based analysis of scattering from 3-D inhomogeneous anisotropic bodies

This paper presents an integral equation based scheme to analyze scattering from inhomogeneous bodies with anisotropic electromagnetic properties. Both the permittivity and permeability are assumed to be generalized tensors. Requisite integral equations are derived using volume equivalence theorem with the electric and magnetic flux densities being the unknown quantities. Matrix equations are derived by discretizing these unknowns using three dimensional Rao-Wilton-Glisson basis functions. Reduction of the integral equation to a corresponding matrix equation is considerably more involved due to the presence of anisotropy and the use of vector basis function; methods for evaluation of the integrals involved in the construction of this matrix is elucidated in detail. The method of moments technique is augmented with the fast multipole method and a compression scheme. The latter two enable large scale analysis. Finally, several numerical results are presented and compared against analytical solutions to validate the proposed scheme. An appendix provides analytical derivations for the formulae that are used to validate numerical method, and the necessary formulae that extends the approach presented herein to the analysis of scattering bianisotropic bodies.     

A higher order parallelized multilevel fast multipole algorithm for3-D scattering

A higher order multilevel fast multipole algorithm (MLFMA) is presented for solving integral equations of electromagnetic wave scattering by three-dimensional (3-D) conducting objects. This method employs higher order parametric elements to provide accurate modeling of the scatterer's geometry and higher order interpolatory vector basis functions for an accurate representation of the electric current density on the scatterer's surface. This higher order scheme leads to a significant reduction in the mesh density, thus the number of unknowns, without compromising the accuracy of geometry modeling. It is applied to the electric field integral equation (EFIE), the magnetic field integral equation (MFIE), and the combined field integral equation (CFIE), using Galerkin's testing approach. The resultant numerical system of equations is then solved using the MLFMA. Appropriate preconditioning techniques are employed to speedup the MLFMA solution. The proposed method is further implemented on distributed-memory parallel computers to harness the maximum power from presently available machines. Numerical examples are given to demonstrate the accuracy and efficiency of the method as well as the convergence of the higher order scheme     

A potential integral equations method for electromagneticscattering by penetrable bodies

A method for computing the electromagnetic scattering by general inhomogeneous penetrable bodies is presented. The method is based on the volume equivalence principle and it uses the electromagnetic potentials as unknowns. The resulting coupled integral equations system is solved by the method of moments in combination with cubical and curvilinear meshes in the special case of purely dielectric scatterers. To show the accuracy of the method, numerical results of the transmitted and of the scattered fields are compared with existing analytical and experimental results     

粗糙面与目标电磁散射统计特性分析

首先给出用于描述和产生粗糙海面的模型，然后在锥形平面波入射条件下，用矩量法求解了粗糙面及与圆柱形近地小目标复合散射所满足的积分方程，最后给出了不同条件下粗糙面、粗糙面与目标复合散射特性的统计结果及相互间的比较。     

Electromagnetic Pulse Penetration into Reinforced-Concrete Buildings

The purpose of this study is to define the protection against electromagnetic pulses (EMP's) provided by reinforced-concrete buildings. We want to define the electromagnetic field inside a building modeled by parallel wires because attenuation due to concrete itself can generally be neglected. The problem resolution is achieved by using integral equations in the time domain. These equations are solved numerically by applying the method of moments associated with the utilization of Lagrange interpolation polynomials. Applications in the field of building protection are shown.     

The use of the transfinite interpolation in the method of momentsapplied to electromagnetic scattering by dielectric cylinders

The method of moments (MoM) solution of electromagnetic scattering presents two major numerical difficulties: the number of unknowns and the computation time necessary to calculate the matrix elements. To circumvent these problems, a MoM using the transfinite interpolation and a reduced integration scheme is presented here. The so-called h and p versions of the new method are applied to the scattering of an electromagnetic wave by an infinite dielectric cylinder (TM case) in the Richmond formulation. The transfinite and classical methods are compared in terms of the convergence rates of the radar cross section and of the total electric field inside the dielectric. The results confirm the superiority of the new schemes as predicted by the theory     

基于特征基函数法的一维理想导体粗糙海面电磁散射快速算法研究

针对传统矩量法在处理具有较多未知量的理想导体粗糙海面电磁散射问题时对计算机内存的需求过大, 耗时过长的缺陷, 文中引入了特征基函数法, 并根据Foldy-Lax多径散射方程构造特征基函数, 首先只考虑离散子域本身的自相互作用构造主要特征基函数, 然后考虑各离散子域间的互耦效应构造次要特征基函数, 最后由主要特征基函数和次要特征基函数的加权叠加构造特征基函数.通过与传统矩量法仿真结果的对比, 讨论了不同次要特征基函数的阶数或不同离散子域的个数对计算精度和计算效率的影响.仿真结果表明了本文所采用的算法能够在保证计算精度的前提下, 减少计算时间, 并能够通过离散子域尺寸的选取控制实际操作矩阵的维数.     

A coupled surface-volume integral equation approach for thecalculation of electromagnetic scattering from composite metallic andmaterial targets

A coupled surface-volume integral equation approach is presented fur the calculation of electromagnetic scattering from conducting objects coated with materials. Free-space Green's function is used in the formulation of both integral equations. In the method of moments (MoM) solution to the integral equations, the target is discretized using triangular patches for conducting surfaces and tetrahedral cells for dielectric volume. General roof-top basis functions are used to expand the surface and volume currents, respectively. This approach is applicable to inhomogeneous material coating, and, because of the use of free-space Green's function, it can be easily accelerated using fast solvers such as the multilevel fast multipole algorithm     

Numerical Analysis of Quasioptical Multireflector Antennas in 2-D With the Method of Discrete Singularities: E-Wave Case

Numerical method for modeling the E-polarized wave scattering by electrically large quasioptical two-dimensional (2-D) reflectors is presented. Reflectors are assumed zero-thickness and perfectly electrically conducting. Efficient numerical solution is obtained from the coupled singular integral equations discretized using new quadrature formulas of interpolation type. It has controlled accuracy and deals with small-size matrices. To simulate a small-horn feeding, the incident field is taken as a beam generated by a complex-source-point (CSP) current. Presented numerical results validate empirical rule of -10 dB edge illumination needed to provide the best electromagnetic performance of reflector     

二维有耗介质目标重建的Newtoin迭代方法

本文给出了一种由已知的散射场数据重建二线非均匀有耗目标的复介电常数的迭代算法。由积分方程出发，利用点匹配技术导出了依赖于未知参数的解析逆散射公式。由此可以以解析的形式计算场量对未知参数的导数（Jacobian和Hessian矩阵）。本文采用Newton优化方法迭代末解道散射问题，具有二次收敛特性。为了克服逆散射中解的不适定性，连续采用多个方向的TM波照射目标，并采集目标区域外的散射场数据，以及采用共轭梯度法（CGM）求解逆问题．数值结果表明了本文所提方法的可行性和灵活性。