Application of the finite element method to Monte Carlo simulations of scattering of waves by random rough surfaces: penetrable case

The finite element method (FEM) of Monte Carlo simulations of random rough surface scattering is extended to penetrable rough surface scattering. The attraction of the method is the banded nature of the resulting matrix equation. The method yields a system of linear algebraic equations which is solved by a direct sparse symmetric matrix inversion. Convergence and accuracy of the method is demonstrated and established by varying various input parameters such as the number of evanescent waves, the number of sampling points and the surface lengths. Results with incident plane wave TE polarization are presented for both the mean reflected scattered intensity and the mean transmitted scattered intensity as a function of surface parameters such as RMS surface heights and correlation lengths. The numerical results are compared against the tapered wave integral equation (TWIE) method. The results of a tapered wave solution of the integral equation averaging over many realizations are in good numerical agreement with FEM if large surface lengths are used in the integral equation method. It is found that a large surface length is required in the TWIE method to have a narrow incident angular spectlum to accurateiy predict the transmitted scattered intensity, whereas a relatively small surface length is sufficient in the FEM. The total CPU time and memory storage requirements for the FEM are much less than that of the TWIE method for eases when the number of horizontal sampling points is much larger than the number of vertical sampling points in the region of discretization. The percentage error in conservation of energy for the FEM is shown to be less than 0.4% for all the examples presented. The total CPU time, memory storage requirements and the percentage error comparisons between the FEM and the TWIE are presented.     

Acoustic scattering from a rough sea surface: the mean field by the integral equation method

The scattering of an acoustic signal incident from below at low angles on a rough sea surface is treated by the integral equation method in the parabolic approximation. Equations are obtained allowing the mean scattered field to be calculated even when the surface causes a large phase modulation in the incident wave. Solutions are found using the method of Laplace transforms and some results are presented for a specific type of rough surface.     

Acoustic scattering from a rough sea surface: the mean field by the integral equation method

Abstract

The scattering of an acoustic signal incident from below at low angles on a rough sea surface is treated by the integral equation method in the parabolic approximation. Equations are obtained allowing the mean scattered field to be calculated even when the surface causes a large phase modulation in the incident wave. Solutions are found using the method of Laplace transforms and some results are presented for a specific type of rough surface.     

Rigorous and approximate methods for modeling wave scattering from a locally perturbed perfectly conducting surface

Rigorous and approximate methods are considered for solving the problem of harmonic plane wave scattering from a plane surface arbitrarily perturbed along one dimension on a finite interval. This problem is treated using the Fredholm integral equations of the second kind and the Kirchhoff and Rayleigh approximations. The estimates of the computational efficiency of the integral equation method and the Rayleigh approximation are compared by calculating fields scattered from random rough surfaces in the resonance region (i.e., when the roughness height is comparable to or smaller than the incident wavelength) for an arbitrary incidence of a plane wave. Scattering patterns calculated using the integral equations and the Kirchhoff approximation are discussed in the case of large-scale random rough surface scattering. Particular attention is paid to scattering at near-grazing incidence.     

Specular propagation over rough surfaces: numerical assessment of Uscinski and Stanek's mean Green's function technique

The purpose of this paper is to numerically evaluate the effectiveness and accuracy of Uscinski and Stanek's mean Green's function technique for computing the mean field of a wave scattered by a rough surface. We present here a direct comparison of this technique with a rigorous numerical method, the forward scattering integral equation method, and another analytical method, the first-order smoothing approximation. Furthermore, we compare the roughness generated equivalent admittance using the three methods. Numerical computations reveal that the scattered field calculated by this technique is not accurate particularly for the equivalent admittance at low grazing angles, even though the mean surface current density is recovered when the wave has traversed several correlation lengths on the surface.     

雪层覆盖地面电磁散射的矩量法

为实现分层介质粗糙面电磁散射的矩量法研究,给出一种分层介质粗糙面电磁积分方程的区域分解方法.将格林定理应用于粗糙面所分的各子空间,结合波动方程和格林函数推导分层粗糙面的电磁积分方程,利用矩量法对其进行离散,数值计算得到雪层覆盖地面散射系数的角分布曲线,其中,粗糙表面由一维带限Weierstrass分形谱和Monte Carlo方法模拟.通过与时域有限差分法数值结果的比对,验证该方法的准确性,并分析散射系数随雪和地面参数、介质参数以及入射波参数的变化,获得了较完整的散射特征.     

Specular propagation over rough surfaces: numerical assessment of Uscinski and Stanek's mean Green's function technique

The purpose of this paper is to numerically evaluate the effectiveness and accuracy of Uscinski and Stanek's mean Green's function technique for computing the mean field of a wave scattered by a rough surface. We present here a direct comparison of this technique with a rigorous numerical method, the forward scattering integral equation method, and another analytical method, the first-order smoothing approximation. Furthermore, we compare the roughness generated equivalent admittance using the three methods. Numerical computations reveal that the scattered field calculated by this technique is not accurate particularly for the equivalent admittance at low grazing angles, even though the mean surface current density is recovered when the wave has traversed several correlation lengths on the surface.     

Time-Domain Scattering Modelling of Two-Dimensional Cylinder Located on a Rough Surface

Numerical modelling on the transient electromagnetic scattering by a two-dimensional （21）） cylinder located on a time-evolving rough surface is presented by using time-domain integral equations. The proposed special choice of a tapered Gauss pulse incident wave removes the truncation error from the rough surface. Additionally, a two-level averaging technique is utilized to overcome the instability from the time marching procedure of solving integral equations. Excellent correspondences between the surface current distributions, as well as the far-zone fields, computed by the proposed method and that obtained by the traditional method of moments associated with the inverse discrete Fourier transformation scheme demonstrate the accuracy of the modelling.     

Direct Solution of the Inverse Problem for Rough Surface Scattering

We consider the inverse scattering problem for a scalar wave field incident on a perfectly conducting one-dimensional rough surface. The Dirichlet Green function for the upper half-plane is introduced, in place of the free-space Green function, as the fundamental solution to the Helmholtz equation. Based on this half-plane Green function, two reasonable approximate operations are performed, and an integral equation is formulated to approximate the total field in the two-dimensional space, then to determine the profile of the rough surface as a minimum of the total field. Reconstructions of sinusoidal, non-sinusoidal and random rough surface are performed using numerical techniques. Good agreement of these results demonstrates that the inverse scattering method is reliable.     

Inversion of synthetic and experimental acoustical scattering data for the comparison of two reconstruction methods employing the Born approximation

This work is concerned with the reconstruction, from measured (synthetic or real) data, of a 2D penetrable fluid-like object of arbitrary cross-section embedded in a fluid of infinite extent and insonified by a plane acoustic wave. Green's theorem is used to provide a domain integral representation of the scattered field. The introduction therein of the Born approximation gives rise to a linearized form of the inverse problem. The actual inversion is carried out by two methods. The first diffraction tomography (DT), exhibits the contrast function very conveniently and explicitly in the form of a wave number/incident angle Fourier transform of the far backscattered field and thus requires measurements of this field for incident waves all around the object and at all frequencies. The second discretized domain integral equation with Born approximation method, is numerically more intensive, but enables a wider choice of configurations and requires less measurements (one or several frequencies, one or several incident waves, choice of measurement points) than the DT method. A comparison of the two methods is carried out by inversion of both simulated and experimental scattered field data.     

SH-wave propagation and resonance phenomena in a periodically layered composite structure with a crack

This paper presents a dynamic analysis of time-harmonic plane SH-waves propagating in periodically multilayered elastic composites with a strip-like crack. The total wave field in the multilayered elastic structure is described as a sum of incident wave field modeled by the transfer matrix method and the scattered wave field governed by an integral representation containing the crack-opening-displacement. The integral equation derived from the boundary conditions on the crack-faces is solved numerically by a Galerkin method. The paper focuses on resonant and non-resonant regimes of anti-plane wave motion in a stack of elastic layers weakened by a single strip-like crack and wave localization in the vicinity of the crack. The scattered extra displacement induced by the presence of the crack is investigated in detail for both situations of high and low contrast in material properties. Numerical results for the average crack-opening-displacement, the transmission coefficient, the stress intensity factor and the average energy flow are presented and discussed to reveal wave resonance and localization phenomena within the band-gaps and the pass-bands.     

粗糙海面高功率微波传输特性矩量法分析

应用蒙特卡罗方法实现了粗糙海面的仿真与模拟,建立了基于双积分方程的高功率微波（HPM）近海面传输特性矩量法计算模型。模型采用光滑窗函数对均匀平面波进行调制,把均匀平面波入射调制为锥形波,消除了粗糙海面突然被截断而引起的边缘效应的影响;重新推导了锥形波入射下的基尔霍夫近似公式,并在满足基尔霍夫近似的条件下,通过对比分析,验证了模型的正确性;采用模型计算分析了不同海面几何参数和海水媒质参数对HPM近海面传输系数的影响。结果表明:粗糙海面的均方根高度对HPM传输系数影响明显,均方根高度越大,传输系数越小,能量分布越均匀;另外随着海水介电常数实部和虚部的增加,传输系数均有所增加,并且实部的影响更明显。     

基于前后向迭代法的3维标量激光漫反射特性

 针对激光点对点通信方式的不足,提出了利用海面作为激光漫反射媒介进行组网通信的设想,并对激光入射理想导体表面后产生的散射场进行了的研究。首先利用3维标量锥形波对入射激光束进行了模拟;然后列出了积分方程,并根据矩量法进行了离散;最后通过前后向迭代法求解矩阵方程,得出了3维散射系数,准确地表示了基于理想导体界面的3维标量激光漫反射特性。     

Unification of the Kirchhoff approximation and the method of moment for optical wave scattering from the lossy dielectric Gaussian random rough surface

The optical wave scattering from one-dimensional （1D） lossy dielectric Gaussian random rough surface is studied. The tapered incident wave is introduced into the classical Kirchhoff approximation （KA）, and the shadowing effect is also taken into account to make the KA results have a high accuracy. The definition of the bistatic scattering coefficient of the modified KA and the method of moment （MOM） are unified. The characteristics of the optical wave scattering from the lossy dielectric Gaussian random rough surface of different parameters are analyzed by implementing MOM.     

不同土壤类型的粗糙地面与其下方埋藏目标复合电磁散射研究

利用MOM研究了不同土壤类型的一维有耗介质粗糙地面与其下方二维埋藏目标的复合电磁散射问题,在数值计算中,引入了锥形入射波以减少粗糙面模拟中由于人为截断所引起的计算误差.当给定入射波频率、土壤成分及土壤温度和土壤湿度后,可以根据相关模型计算出土壤的相对介电常数.数值结果以二维无限长导体圆柱为例,利用 Monte Carlo方法统计计算了双站复合散射截面随土壤类型、土壤湿度、入射波频率及均方根高度等的变化情况,并做了详细分析及讨论.该研究结果对于探地雷达在民用、军事、工业等多领域的应用具有一定的理论指导价值. 关键词： 粗糙地面 土壤类型 相对介电常数 四成分模型     

Numerical analysis of immersed finite cylindrical shells using a coupled BEM/FEM and spatial spectrum approach

This study numerically analyzes submerged cylindrical shells using a coupled boundary element method (BEM) with finite element method (FEM) in conjunction with the wave number theory, in which the spatial Fourier transform of surface velocity for cylinders is directly related to pressure in a far field. The acoustic loading is formulated using a symmetric complex matrix derived from a boundary integral equation where the symmetry is based on an acoustic reciprocal principle for surface acoustics. In this formulation the acoustic loading matrix is a large acoustic element whose degree of freedom is connected to the normal displacement of the vibrating structures. The coupled BEM/FEM equation is a banded, symmetric matrix, and thus its bandwidth can be minimized using a proper algorithm. This formulation significantly increases numerical efficiency. The computed normal velocity is thus transformed to wave number representation to examine acoustic radiation. A finite plane cylindrical shell, without attached stiffeners, and a shell with internal ring stiffeners are chosen to demonstrate the present analysis procedure. The far field pressure computed directly from the integral equation and predicted by wave number theory correlates closely with increasing vibrating frequency. Meanwhile, the influences of the internal ring structures on acoustic radiation are examined using the wave number theory, which helps in understanding how internal structures influence radiated noise.     

Scattering and diffraction of a plane wave by a randomly rough half-plane

The scattering and diffraction of a TE (transverse electric) plane wave by a randomly rough half-plane are studied by a combination of three techniques: the Wiener-Hopf technique, the small perturbation method and a probabilistic method based on the shift-invariance of a homogeneous random function. By use of the D^a-Fourier transformation based on the shift-invariance, it is shown that the scattered wave is written by an inverse Fourier transformation of a homogeneous random function with a complex parameter. For a small rough case, such a random function with a complex parameter is expanded in a perturbation series and then the first-order solution is obtained exactly in an integral form. The first-order solution involves two physical processes such that the edge-diffracted wave is scattered by the randomly rough plane and the scattered wave, due to roughness, is diffracted by the half-plane. The solution is transformed into a sum of the Fresnel integrals with complex arguments, an integral along the steepest descent path and a branch-cut integral, which are evaluated numerically. Then, intensities of the coherently scattered wave and incoherent wave are calculated in the region near the edge and illustrated in figures.     

Scattering and diffraction of a plane wave by a randomly rough half-plane*

Abstract

The scattering and diffraction of a TE (transverse electric) plane wave by a randomly rough half-plane are studied by a combination of three techniques: the Wiener-Hopf technique, the small perturbation method and a probabilistic method based on the shift-invariance of a homogeneous random function. By use of the D^a-Fourier transformation based on the shift-invariance, it is shown that the scattered wave is written by an inverse Fourier transformation of a homogeneous random function with a complex parameter. For a small rough case, such a random function with a complex parameter is expanded in a perturbation series and then the first-order solution is obtained exactly in an integral form. The first-order solution involves two physical processes such that the edge-diffracted wave is scattered by the randomly rough plane and the scattered wave, due to roughness, is diffracted by the half-plane. The solution is transformed into a sum of the Fresnel integrals with complex arguments, an integral along the steepest descent path and a branch-cut integral, which are evaluated numerically. Then, intensities of the coherently scattered wave and incoherent wave are calculated in the region near the edge and illustrated in figures.     

三维粗糙面电磁双站散射的直接型区域分解计算

提出三维粗糙面双站电磁散射的直接型有限元-区域分解方法.首先建立含有迭代Robin边界条件(IRBC)的区域分解法耦合模型,再用内视法导出高度稀疏分块的分区耦合矩阵,之后给出缩减耦合矩阵带宽的子区域排序方法和IRBC的FFT加速算法.用有限元-完全匹配层和未分区的有限元-IRBC方法验证数值结果.     

高斯分层介质粗糙面光波透射问题的微扰法研究

运用微扰法研究平面波入射分层介质粗糙面的光波透射问题,推出了不同极化状态的透射光波散射系数公式。采用高斯粗糙面来模拟实际的分层介质粗糙面,结合高斯粗糙面的功率谱导出了平面波入射高斯分层介质粗糙面的透射系数计算公式。通过数值计算得到HH极化透射系数随透射光波散射角变化的曲线,讨论底层介质介电常数、中间介质介电常数和厚度、粗糙面参数及入射光波长对透射系数的影响。数值计算结果表明：底层介质介电常数、中间介质介电常数和厚度、粗糙面参数及入射光波长对透射系数的影响是非常复杂的。