时域混合算法在一维海面与舰船目标复合电磁散射中的应用

采用时域积分方程(TDIE)与时域基尔霍夫近似(TDKA)的混合算法研究粗糙海面与舰船目标的复合瞬态电磁散射.该方法将舰船目标及其近邻海面划分为TDIE区域,用TDIE方法精确求解;将剩余电大尺寸的粗糙海面划分为TDKA区域,采用高效的TDKA电流近似求解.通过混合算法和传统TDIE算法结果的对比,表明TDIE-TDKA混合算法能保证计算的精度,同时具有较高的计算效率.最后,讨论了海面上方有无目标、海面上方风速、电磁脉冲入射角、舰船目标尺寸、吃水深度对后向散射磁场的影响.     

新型压缩感知计算模型分析三维电大目标电磁散射特性

为提高基于压缩感知技术的矩量法在三维电大目标双站电磁散射问题中的计算效率和稳定性,提出新的稀疏、测量和重构方法,构建一种新型压缩感知计算模型.不同于基于欠定方程的传统的压缩感知计算模型,新型计算模型首先采用按行均匀抽取阻抗矩阵的方法构造测量矩阵以获得稳定的计算结果;然后,基于Foldy-Lax方程生成多阶特征基函数并作为稀疏基对感应电流进行稀疏转换;再依据少数低阶特征基函数足以近似表征感应电流的先验条件,将恢复算法简化为最小二乘法;最后,将矩阵方程转换为一个超定系统并采用最小二乘法解出电流系数.与传统的计算模型相比,新型计算模型不仅可以获得更加稳定的精确解,还可以显著提高电大目标双站散射问题的求解效率和计算精度.数值仿真结果证明了新方法的可行性和高效性.     

由卷积型变分原理求解瞬态热传导问题的半解析法

主要针对二维问题，基于卷积型变分原理，提出了求解瞬态热传导问题的半解析法，该方法在空间域内作分区域插值离散，在时间域上采用解析函数。     

用变分优化正交的连带Laguerre基函数计算无支撑单层MoS2中激子的能量和波函数

采用基于变分优化正交化连带Laguerre基函数的准确对角化方法计算了无支撑单层MoS2中A型激子的能量和波函数.介电屏蔽效应破坏了SO(3)对称性导致激子能量以反常的轨道角动量顺序出现.该方法利用连带Laguerre多项式构造了既满足束缚态又满足连续态的正交基函数集,并推导了哈密顿量矩阵元的解析表达式.收敛速度不仅与基函数的数量有关,而且与变分参数的大小有关.在变分优化下,收敛速度非常快,表明了该方法的可靠性.采用正交化连带Laguerre基矢可大大减小基函数的数量和计算量.我们计算的激子本征能量即使是在很少的基函数下也与文献中的结果非常吻合.变分优化二维正交归一化连带Laguerre基矢适用于二维材料中激子和原子物理的精确描述.     

改进的二维分形模型在海面电磁散射中的应用

提出了一种改进的二维分形海面模型，其表面谱函数在空间波数小于基波波数及大于基波波数时分别满足正幂率关系和负幂率关系.通过比较可以发现在不同风速时，改进模型的空间自相关函数及表面轮廓谱和有关文献结果有较好的吻合.在满足Kirchhoff近似条件下推导了改进分形模型的散射系数及散射强度系数的计算公式并进行了数值计算，比较了改进模型和经典模型的后向散射强度系数角分布并详细讨论了它们随入射频率、海上风速和风向的变化. 关键词： 改进分形模型 粗糙海面 电磁散射 Kirchhoff近似     

三维目标与粗糙面复合散射的广义稀疏矩阵平面迭代及规范网格算法

提出了快速计算二维导体粗糙面与面上金属目标复合散射的广义稀疏矩阵平面迭代及规范网格法(G-SMFSIA/CAG).推导了二维导体粗糙面与面上目标相互作用的耦合积分方程,用稀疏矩阵平面迭代及规范网格法(SMFSIA/CAG)求解粗糙面部分的表面积分方程,而用基于RWG基函数的矩量法(MOM)计算目标部分的表面积分方程,并通过更新方程的激励项迭代求解目标与粗糙面的相互耦合作用.结合Monte-Carlo方法产生具有PM(Pierson-Moskowitz)海浪谱的随机海洋粗糙面,数值分析了海面上不同形状导体目 关键词： 复合散射 广义稀疏矩阵平面迭代及规范网格法 随机海洋粗糙面 双站散射系数     

金属纳米柱天线的快速电磁分析

基于分段正弦基函数和矩量法,通过求解离散电流节点格林函数的封闭解得到金属纳米柱天线激发表面等离子体的阻抗矩阵.与使用其它基函数矩量法相比,该方法可以减少矩阵方程的维数.仿真结果表明:使用此方法只需求解很小的矩阵方程就可以求解出纳米天线表面极化电流,从而实现对纳米天线的散射特征及谐振模式的快速分析;其结果与时域有限差分仿真结果吻合良好且速度具有显著的优势,尤其在计算斜入射问题时计算优势更加明显.本文的方法对文中计算的模型有效,同样为其他形状纳米柱天线和碳纳米管器件散射特性仿真提供了快速有效的电磁分析方法.     

基于数值方法的激光海面漫反射通信能量分布特性

 针对激光点对点通信方式的不足,提出了利用海面作为激光漫反射媒介进行组网通信的设想,采用前后向迭代的数值方法结合Green函数谱积分加速算法(FBM/SAA)对激光海面漫反射通信的能量分布特性进行了研究。通过对激光光束入射海面后产生的散射场的分析、计算及实验验证,得出了较为准确的2维激光海面双站散射系数,并对激光光束入射海面后的散射场进行了分析。研究结果表明：入射激光束经过粗糙海面散射后,能量大部分集中在前向散射区域上,而后向散射强度很弱,且在散射场的边缘处能量迅速衰减,说明了激光海面漫反射组网通信方法的可行性。通过与前后向迭代法和Kirchhoff近似方法的计算结果的比较,说明了FBM/SAA是一种高效、准确的计算方法。     

二维粗糙海面的光散射及其红外成像

首先根据JONSWAP海面功率谱模型数值模拟出二维粗糙海面，采用几何光学近拟与基尔霍夫（Kirchhoff）标量近似计算了二维海面的光散射，计算中将每一面元看成一具有微粗糙度的粗糙面而不是近似地当作平面，并利用投影法与射线追踪法数值计算了一定入射角和散射角下的遮挡函数，有效地提高了海面光散射计算的精确性。最后利用太阳光的光谱辐照度数值模拟了海面的3μm-5μm红外散射图像，对于红外探测器抑制海面反射太阳光造成的亮带干扰具有一定的参考价值。     

基于压缩感知的一维海面与二维舰船复合后向电磁散射快速算法研究

矩量法作为数值方法中积分方程方法的代表, 具有计算精度高、所用格林函数自动满足辐射条件、无须额外设置边界条件等优点. 但是在舰船目标与海面复合后向电磁散射仿真中, 传统矩量法需针对每个入射角反复求解矩阵方程组, 导致其在处理后向散射问题时计算量大, 耗时长, 仿真效率低下. 为解决上述问题, 本文提出了一种基于压缩感知技术的矩量法的改进算法. 该算法在求解复合后向散射问题时, 首先利用观测矩阵与传统矩量法中的电压矩阵相乘, 得到一组新的低维度的电压矩阵; 其次通过求解新电压矩阵下的矩阵方程组, 获得电流矩阵的观测值; 最后利用恢复算法(本文采用正交匹配追踪算法)重构出所需的原始入射源照射下的电流系数. 通过与传统矩量法的计算结果对比, 表明本文所提算法能够在保证计算精度的前提下, 明显减少计算时间, 提高计算效率.     

Modelling of the electromagnetic scattering by sea surfaces at grazing incidence. Application to HF surface wave radars

The object of the study is to model the electromagnetic scattering by time-evolving sea surfaces illuminated by high frequency surface wave radars. A common way to simulate this case is to consider the scattering at grazing incidence. Recent studies have focused on so-called exact methods. As such methods are very time and memory consuming, it becomes important to implement techniques to reduce the calculation costs. Another important issue is to better understand the interactions between the electromagnetic waves and moving sea surfaces. We present in this article our simulator main characteristics and a comparison between simulated data and measurements.     

Theoretical and computational aspects of scattering from periodic surfaces: one-dimensional transmission interface

We consider the scattering from and transmission through a one-dimensional periodic surface. For this problem, the electromagnetic cases of TE and TM polarization reduce to the scalar acoustic examples. Three different theoretical and computational methods are described, all involving the solution of integral equations and their resulting discrete matrix system of equations for the boundary unknowns. They are characterized by two sample spaces for their discrete solution, coordinate space and spectral space, and labelled by the sampling of the rows and columns of the discretized matrices. They are coordinate-coordinate (CC), the usual coordinate-space method, spectral-coordinate (SC) where the matrix rows are discretized or sampled in spectral space and spectral-spectral (SS) where both rows and columns are sampled in spectral space. The SS method uses a new topological basis expansion for the boundary unknowns.

Equations are derived for infinite surfaces, then specialized and solved for periodic surfaces. Computational results are presented for the transmission problem as a function of roughness, near-grazing incidence as well as many other angles, density and wavenumber ratios. Matrix condition numbers and different sampling methods are considered. An error criterion is used to gauge the validity of the results.

The computational results indicated that the SC method was by far the fastest (by several orders of magnitude), but that it became ill-conditioned for very rough surfaces. The CC method was most reliable, but often required very large matrices and was consequently extremely slow. It is shown that the SS method is computationally efficient and accurate at near-grazing incidence and can be used to fill a gap in the literature. Extensive computational results indicate that both SC and SS are highly robust computational methods. Spectral-based methods thus provide viable computational schemes to study periodic surface scattering.     

Theoretical and computational aspects of scattering from periodic surfaces: one-dimensional transmission interface

Abstract

We consider the scattering from and transmission through a one-dimensional periodic surface. For this problem, the electromagnetic cases of TE and TM polarization reduce to the scalar acoustic examples. Three different theoretical and computational methods are described, all involving the solution of integral equations and their resulting discrete matrix system of equations for the boundary unknowns. They are characterized by two sample spaces for their discrete solution, coordinate space and spectral space, and labelled by the sampling of the rows and columns of the discretized matrices. They are coordinate-coordinate (CC), the usual coordinate-space method, spectral-coordinate (SC) where the matrix rows are discretized or sampled in spectral space and spectral-spectral (SS) where both rows and columns are sampled in spectral space. The SS method uses a new topological basis expansion for the boundary unknowns.

Equations are derived for infinite surfaces, then specialized and solved for periodic surfaces. Computational results are presented for the transmission problem as a function of roughness, near-grazing incidence as well as many other angles, density and wavenumber ratios. Matrix condition numbers and different sampling methods are considered. An error criterion is used to gauge the validity of the results.

The computational results indicated that the SC method was by far the fastest (by several orders of magnitude), but that it became ill-conditioned for very rough surfaces. The CC method was most reliable, but often required very large matrices and was consequently extremely slow. It is shown that the SS method is computationally efficient and accurate at near-grazing incidence and can be used to fill a gap in the literature. Extensive computational results indicate that both SC and SS are highly robust computational methods. Spectral-based methods thus provide viable computational schemes to study periodic surface scattering.     

Electromagnetic volume scattering in the half-space: a moment-method analysis without Sommerfeld integrals or their approximations

Traditionally, in moment-method analyses of electromagnetic scattering, the elements of the impedance matrix are calculated as convolutions of the basis elements with the appropriate dyadic Green's function. However, for scattering in the half-space, the vertical and azimuthal copolar terms of the Green's function require evaluation of Sommerfeld integrals which are computationally burdensome. In this paper, it is shown that, in populating the impedance matrix for the half-space problem, evaluation of Sommerfeld integrals is, in fact, not necessary. For monochromatic excitation, the plane-wave expansion of the scattered field constitutes a Fourier transform, in the horizontal plane, of a vector spectral function. This vector function results from the convolution, in the vertical dimension, of the respective angular spectra of the Green's function and the equivalent current. On application of the moment method, through the Weyl identity, the impedance-matrix elements corresponding to the singular terms of the Green's function are convolutions in the horizontal plane of spherical potentials, and Fourier transforms of scalar spectral functions. These scalar functions are derived from the basis elements and, with a judicious choice of basis, they are well behaved and of compact support, and consequently their Fourier transforms can be computed as FFTs.     

一种基于压缩感知的三维导体目标电磁散射问题的快速求解方法

提出了一种将压缩感知和特征基函数结合的方法来计算三维导体目标的雷达散射截面.利用压缩感知理论,将随机选择的矩量法阻抗矩阵作为测量矩阵,将激励电压视为测量值,然后再用恢复算法可实现二维或二维半目标感应电流的求解.对于三维导体目标,使用Rao-Wilton-Glisson基函数表示的感应电流在常用的离散余弦变换基、小波基等稀疏基上不稀疏.为此,本文将计算出的目标特征基函数作为稀疏基,用广义正交匹配追踪算法作为恢复算法来加速恢复过程,并应用到三维导体目标的雷达散射截面计算中.数值结果证明了本文方法的准确性与高效性.     

Theoretical and computational aspects of scattering from periodic surfaces: two-dimensional perfectly reflecting surfaces using the spectral-coordinate method

We consider the scattering from a two-dimensional periodic surface. From our previous work on scattering from one-dimensional surfaces (1998 Waves Random Media 8 385) we have learned that the spectral-coordinate (SC) method was the fastest method we have available. Most computational studies of scattering from two-dimensional surfaces require a large memory and a long calculation time unless some approximations are used in the theoretical development. By using the SC method here we are able to solve exact theoretical equations with a minimum of calculation time.

We first derive in detail (part I) the SC equations for scattering from two-dimensional infinite surfaces. Equations for the boundary unknowns (surface field and/or its normal derivative) result as well as an equation to evaluate the scattered field once we have solved for the boundary unknowns. Special cases for the perfectly reflecting Dirichlet and Neumann boundary value problems are presented as is the flux-conservation relation.

The equations are reduced to those for a two-dimensional periodic surface in part II and we discuss the numerical methods for their solution. The two-dimensional coordinate and spectral samples are arranged in one-dimensional strings in order to define the matrix system to be solved.

The SC equations for the two-dimensional periodic surfaces are solved in part III. Computations are performed for both Dirichlet and Neumann problems for various periodic sinusoidal surface examples. The surfaces vary in roughness as well as period and are investigated when the incident field is far from grazing incidence ('no grazing') and when it is near-grazing. Extensive computations are included in terms of the maximum roughness slope which can be computed using the method with a fixed maximum error as a function of the azimuthal angle of incidence, the polar angle of incidence and the wavelength-to-period ratio.

The results show that the SC method is highly robust. This is demonstrated with extensive computations. Furthermore, the SC method is found to be computationally efficient and accurate for near-grazing incidence. Computations are presented for grazing angles as low as 0.01°. In general, we conclude that the SC method is a very fast, reliable and robust computational method to describe scattering from two-dimensional periodic surfaces. Its major limiting factor is high slopes and we quantify this limitation.     

Electromagnetic volume scattering in the half-space: a moment-method analysis without Sommerfeld integrals or their approximations

Abstract

Traditionally, in moment-method analyses of electromagnetic scattering, the elements of the impedance matrix are calculated as convolutions of the basis elements with the appropriate dyadic Green's function. However, for scattering in the half-space, the vertical and azimuthal copolar terms of the Green's function require evaluation of Sommerfeld integrals which are computationally burdensome. In this paper, it is shown that, in populating the impedance matrix for the half-space problem, evaluation of Sommerfeld integrals is, in fact, not necessary. For monochromatic excitation, the plane-wave expansion of the scattered field constitutes a Fourier transform, in the horizontal plane, of a vector spectral function. This vector function results from the convolution, in the vertical dimension, of the respective angular spectra of the Green's function and the equivalent current. On application of the moment method, through the Weyl identity, the impedance-matrix elements corresponding to the singular terms of the Green's function are convolutions in the horizontal plane of spherical potentials, and Fourier transforms of scalar spectral functions. These scalar functions are derived from the basis elements and, with a judicious choice of basis, they are well behaved and of compact support, and consequently their Fourier transforms can be computed as FFTs.     

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.     

Transfer-matrix for the multichannel scattering problem

A transfer-matrix for the multichannel scattering problem is obtained. The elements of this matrix are expressed in terms of transmission and reflection amplitudes. On the basis of the matrix for a system of N localized and nonoverlapped scattering centers the recurrent equations for the transfermatrix elements are derived and the initial conditions are defined.