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.     

Spectral domain analysis of coupled microstrip lines on magnetized ferrite substrates

Dispersion characteristics of coupled microstrip lines on ferrimagnetic substrates are presented and their features discussed. The Hertz vector potentials in the spectral domain method are used to calculate the dyadic Green's functions in an impedance matrix form. The behavior of the propagation constants is obtained by using the spectral Galerkin's technique. Numerical results are found as a function of various geometrical parameters and of the externally applied static magnetic field. Nonreciprocal behavior is mainly observed when the static magnetic field is applied in the transverse direction, parallel to the ground plane.     

Nonlocal transport of passive scalars in turbulent penetrative convection

We present a Green's function approach for quantifying the transport of a passive scalar (tracer) field in three-dimensional simulations of turbulent convection. Nonlocal, nondiffusive behavior is described by a transilient matrix (the discretized Green's function), whose elements contain the fractional tracer concentrations moving from one subvolume to another as a function of time. The approach was originally developed for and applied to geophysical flows, but here we extend the formalism and apply it in an astrophysical context to three-dimensional simulations of turbulent compressible convection with overshoot into convectively stable bounding regions. We introduce a novel technique to compute this matrix in a single simulation by advecting labeled particles rather than solving the passive scalar equation for a large number of different initial conditions. The transilient matrices thus computed are used as a diagnostic tool to quantitatively describe nonlocal transport via matrix moments and transport coefficients in a generalized, multiorder diffusion equation. Results indicate that transport in both the vertical and horizontal directions is strongly influenced by the presence of coherent velocity structures, generally resembling ballistic advection more than diffusion. The transport of a small fraction of tracer particles deep into the underlying stable region is reasonably efficient, a result which has possible implications for the problem of light-element depletion in late-type stars.     

The complex equivalent source method for sound propagation over an impedance plane

The sound field caused by a monopole source above an impedance plane can be calculated by using a superposition of equivalent point sources located along a line in the mirror space below the plane. Originally, such an approach for representing the half-space Green's function was described by Sommerfeld at the beginning of the last century, in order to treat half-space problems of heat conduction. However, the representation converges only for masslike impedances and cannot be used for the more important case of reflecting planes with springlike surface impedances. The singular part of the line integral can be transformed into a Hankel function, which shows that surface waves are contained in the whole solution. Unfortunately, this representation suffers from the lack of validity at certain receiver points and from restrictions on wave number and impedance range to ensure the necessary convergence. The main idea of the present method is to use also a superposition of equivalent point sources, but to allow that these sources can be located at complex source points. The corresponding form of the half-space Green's function is suitable for both masslike and springlike surface impedances, and can be used as a cornerstone for a boundary element method.     

水平层状介质并矢Green函数的递推矩阵方法

采用递推矩阵方法计算任意数目水平层状介质的并矢Green函数.根据层界面处电场和磁场的连续性条件得到3个确定Sommerfeld积分待定系数的矩阵方程组,分别对应于垂向单位电偶极子产生的TM波、水平方向单位电偶极子产生的TE波和TM波,这些方程组均可通过递推方法快速求解.只需改变3个方程组中源项元素的位置,就可以方便地得到当源点和场点在任意层时的并矢Green函数.本文给出的并矢Green函数表达式形式简洁且不含指数增加项,计算时不会出现溢出现象.     

Scattering of a scalar wave from a two-dimensional randomly rough Neumann surface

We present a reciprocity and unitarity preserving formulation of the scattering of a scalar plane wave from a two-dimensional, randomly rough surface on which the Neumann boundary condition is satisfied. The theory is formulated on the basis of the Rayleigh hypothesis in terms of a single-particle Green's function G(q_‖|k_‖) for the surface electromagnetic waves that exist at the surface due to its roughness, where k_‖ and q_‖ are the projections on the mean scattering plane of the wave vectors of the incident and scattered waves, respectively. The specular scattering is expressed in terms of the average of this Green's function over the ensemble of realizations of the surface profile function (G(q_‖|k_‖)). The Dyson equation satisfied by (G(q_‖|k_‖)) is presented, and the properties of the solution are discussed, with particular attention to the proper self-energy in terms of which the averaged Green's function is expressed. The diffuse scattering is expressed in terms of the ensemble average of a two-particle Green's function, which is the product of two single-particle Green's functions. The Bethe-Salpeter equation satisfied by the averaged two-particle Green's function is presented, and properties of its solution are discussed. In the small roughness limit, and with the irreducible vertex function approximated by the sum of the contribution from the maximally-crossed diagrams, which represent the coherent interference between all time-reversed scattering sequences, the solution of the Bethe-Salpeter equation predicts the presence of enhanced backscattering in the angular dependence of the intensity of the waves scattered diffusely.     

Experimental research of the interference and phase structure of the power flux from a local source in shallow water

The paper studies the interference structure of low-frequency tonal and wideband signals in shallow water, received by four-component vector-scalar modules. The spatial amplitudes and phase characteristics of the scalar field are analyzed, as well as three components of the vibration velocity vector and the power flux vector. A relationship is established between the zone of interference maxima and minima and the phase gradient in the horizontal and vertical plane, the change in direction of the vertical and horizontal components of the vibration velocity vector, and the change in the depression angle of the power flux vector in the horizontal plane.     

The fast multipole method and Fourier convolution for the solution of acoustic scattering on regular volumetric grids

The fast multipole method (FMM) is applied to the solution of large-scale, three-dimensional acoustic scattering problems involving inhomogeneous objects defined on a regular grid. The grid arrangement is especially well suited to applications in which the scattering geometry is not known a priori and is reconstructed on a regular grid using iterative inverse scattering algorithms or other imaging techniques. The regular structure of unknown scattering elements facilitates a dramatic reduction in the amount of storage and computation required for the FMM, both of which scale linearly with the number of scattering elements. In particular, the use of fast Fourier transforms to compute Green’s function convolutions required for neighboring interactions lowers the often-significant cost of finest-level FMM computations and helps mitigate the dependence of FMM cost on finest-level box size. Numerical results demonstrate the efficiency of the composite method as the number of scattering elements in each finest-level box is increased.     

An analytic and numerical solution with spectral Green's function method for transport equation in spherical geometry

Since in many cases curvilinear geometry is more appropriate than cartesian geometry for precise modeling of the complex systems for reactor calculation, we have developed the spectral Green's function (SGF) method which is employed to obtain angular and scalar flux distributions in heterogeneous sphere geometry with isotropic scattering. In this study, we showed that the neutron transport problems of homogeneous spheres could be reduced to the solution of plane geometry equation.Finally, some results are discussed and compared with those already obtained by diamond difference scheme to test the accuracy of the results. The agreement is satisfactory. SGF method is very suitable for the numerical solution of the neutron transport equation with isotropic scattering.     

Iterative approach of high-order scattering solution for vector radiative transfer of inhomogeneous random media

By stratifying a random scatter media into multiple thin layers in the vertical z direction, the first-order scattering solution of each thin layer is employed to derive high-order scattering solution of whole random media. Using the Fourier transform and Mueller matrices in discrete ordinates, an iterative approach to solve high-order scattering solution of vector radiative transfer (VRT) equation is newly developed. Numerical results are well compared with the Mueller matrix solutions of the first order for a single layer medium, second order for a half-space, and the results of the discrete ordinate and eigen analysis method. It demonstrates our approach as feasible, effective and especially applicable to high-order solution of VRT for both bistatic scattering and thermal emission of inhomogeneous non-spherical scatter media.     

基于新型混合场积分方程的半空间三维均匀介质目标电磁散射分析

利用介质半空间格林函数,将一种混合场积分方程(JMCFIE)推广应用到有耗半空间的均匀介质的散射分析中,以减少离散矩阵方程迭代求解时所需的迭代步数,提高矩量法(MoM)对半空间介质目标散射分析的求解效率.数值算例验证了方法的准确性和有效性.     

Closed form solutions for the acoustical impulse response over a masslike or an absorbing plane

The transient sound field caused by a Dirac delta impulse function above an infinite locally reacting plane can be calculated by applying the inverse Fourier transform of the corresponding half-space Green's function in frequency domain. As a starting point, the representation given by Ochmann [J. Acoust. Soc. Am. 116(6), 3304-3311 (2004)] is used, which consists of discrete and continuous superposition of point sources. For a locally reacting plane with masslike character and also with pure absorbing behavior, it is possible to express the resulting impulse response in closed form. Such a result is surprising, because corresponding formulations in the frequency domain are not available yet. Hence, the first main result is the closed form solution Eq. (22) for an impulse response over an infinite plane with a pure imaginary impedance. The second main result is the closed form solution Eq. (53) for an impulse response over an infinite plane with a pure real impedance. As a particular application of both main results, a convolution technique is used for deriving formulas for point sources with a general time dependency. For special signals like an exponentially decaying time signal or a triangular shaped impulse, the resulting sound field can be presented in terms of elementary functions.     

半空间电大涂敷目标散射的高频分析方法

研究了半空间内(海面,地面)电大尺寸涂敷目标散射的高频求解方法.将半空间并矢格林函数引入物理光学方法中,对半空间复杂环境影响进行考虑,并利用阻抗边界条件考虑涂敷目标表面的复杂电磁散射,推导出半空间物理光学分析方法,同时结合图形电磁学,对半空间电大涂敷目标进行消隐判断,提取像素面元法矢量和深度缓存等有效信息,快速有效地计算了半空间电大涂敷目标的雷达散射截面.数值结果证明了方法的有效性和准确性. 关键词： 半空间物理光学方法 半空间并矢格林函数 图形电磁学 雷达散射截面     

Three-dimensional vector radiative transfer in a semi-infinite, Rayleigh scattering medium exposed to spatially varying polarized radiation: generalized reflection matrix

Three-dimensional vector radiative transfer in a semi-infinite medium exposed to spatially varying, polarized radiation is studied. The problem is to determine the generalized reflection matrix for a multiple scattering medium characterized by a 4×4 scattering matrix. A double integral transform is used to convert the three-dimensional vector radiative transfer equation to a one-dimensional form, and a modified Ambarzumian's method is then applied to derive a nonlinear integral equation for the generalized reflection matrix. The spatially varying backscattered radiation for an arbitrarily polarized incident beam can be found from the generalized reflection matrix. For Rayleigh scattering and normal incidence and emergence, the generalized reflection matrix is shown to have five non-zero elements. Benchmark results for these five elements are presented and compared to asymptotic results. When the incident radiation is polarized, the vector approach used in this study correctly predicts three-dimensional behavior, while the scalar approach does not. When the incident radiation is unpolarized, both the vector and scalar approaches predict a two-dimensional distribution of the intensity, but the error in the scalar prediction can be as high as 20%.     

基于半空间FDTD的近远场外推方法: TE情形

研究基于半空间时域有限差分(FDTD)计算中TE情形的远区散射场计算方法.先将半空间Green函数代入势函数,推导远区纵向磁场的计算公式,对各空间内电磁流和对应的相位进行详细分析.再计算两个半空间模型的远区散射场,结果表明:该方法是正确和有效的.针对日本96式装甲车的计算发现,单站雷达探测系统位于战车正上方时具有最好的探测效果.     

Resonance analysis of a 2D alluvial valley subjected to seismic waves

The T-matrix formalism and an ultrasonic experiment are developed to study the scattering of in-plane waves for an alluvial valley embedded in a two-dimensional half-space. The solution of the in-plane scattering problem can be determined by the T-matrix method, where the basis functions are defined by the singular solutions of Lamb's problems with surface loading in both horizontal and vertical directions. In the experiment, a thin steel plate with a semicircular aluminum plate attached on the edge is used to simulate the two-dimensional alluvial valley in the state of plane stress. Based on the spectra of displacement signals measured at the free edge of the scatterer, the resonance frequencies where the peaks appear can be identified. It can be shown that the nondimensional resonance frequency is one of the characteristic properties of the scattering system. Furthermore, it is noted that the nondimensional resonance frequencies measured experimentally are in good agreement with those calculated theoretically.     

并行时域电场积分方程方法在动态海面二维电磁散射中的应用

采用并行时域电场积分方程方法对动态海面的二维瞬态散射特性进行研究。为了保证该方法的后期稳定性,时间基函数和空间基函数采用二阶B样条基函数和三角基函数,矩阵元素采用时间维度精确解析、空间奇异部分精确解析进行计算;为了减少对无限海面进行截断带来的边缘效应,入射波采用锥形调制高斯脉冲;结合信息传递接口(MPI)技术和稀疏矩阵压缩存储技术,对不同时刻的海面进行瞬态散射分析。大量的数值算例证明了该方法在计算动态海面的二维瞬态散射问题时的正确性,还可以保证后期的稳定性,提高计算效率,减少对计算机内存需求。     

Relations between chirp transform and Fresnel diffraction, Wigner distribution function and a fast algorithm for chirp transform

Two physical interpretations of chirp transform related to Fresnel diffraction and Wigner distribution function are given.The chirp transform can be regarded as a Fresnel diffraction observed on a spherical tangent to the diffraction plane,or a rotation and stretching transformation of the Wigner distribution function space.A general fast algorithm for the numerical calculation of chirp transform is developed by employing two fast Fourier transform algorithms.The algorithm,by which a good evaluation can be achieved,unifies the calculations of Fresnel diffraction,arbitrary fractionalorder Fourier transforms and other scalar diffraction systems.The algorithm is used to calculate the Fourier transform of a Gaussian function and the Fourier transform,the Fresnel transform,the Fractional-order Fourier transforms of a rectangle function to evaluate the performance of this algorithm.The calculated results are in good agreement with the analytical results,both in the amplitude and phase.     

跨界面目标电磁散射Sommerfeld积分的双重广义函数束拟合离散复镜像方法

混合势积分方程（mixed-potential electric field integral equation, MPIE）由于其Green函数的低阶奇异性是矩量法（MoM）等数值方法求解散射问题的主要方程.在跨界面目标的散射问题中,MPIE中的矢量势和标量势Green函数包含Sommerfeld类型的谱域积分,利用离散复镜像法（discrete complex image method, DCIM）和Sommerfeld恒等式,将其转化为有限项复镜像Green函数的求和运算,避免了烦杂的谱域积分运算. 关键词： Sommerfeld积分 离散复镜像法 广义函数束