Scattering of electromagnetic waves from two-dimensional perfectly conducting random rough surfaces – study with the curvilinear coordinate method

We present a method giving the bi-static scattering coefficient of two-dimensional (2-D) perfectly conducting random rough surface illuminated by a plane wave. The theory is based on Maxwell's equations written in a nonorthogonal coordinate system. This method leads to an eigenvalue system. The scattered field is expanded as a linear combination of eigensolutions satisfying the outgoing wave condition. The boundary conditions allow the scattering amplitudes to be determined. The Monte Carlo technique is applied and the bi-static scattering coefficient is estimated by averaging the scattering amplitudes over several realizations. The random surface is represented by a Gaussian stochastic process. Results are compared to published numerical and experimental data. Comparisons are conclusive.     

Scattering of a plane wave by one-dimensional dielectric random rough surfaces—study with the curvilinear coordinate method

We present a method giving the bi-static scattering coefficient of a one-dimensional dielectric random rough surface illuminated by a plane wave. The theory is based on Maxwell's equations written in a nonorthogonal coordinate system. For each medium, this method leads to an eigenvalue system. The scattered field is expanded as a linear combination of eigensolutions satisfying the outgoing wave condition. The boundary conditions allow the diffraction amplitudes to be determined. The Monte Carlo technique is applied and the bi-static scattering coefficient is estimated by averaging the scattering amplitudes over several realizations. The results of a Gaussian random process with a Gaussian roughness spectrum are compared to published experimental and numerical data. Comparisons are conclusive.     

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.     

Implementation and validation of the curvilinear coordinate method for the scattering of electromagnetic waves from two-dimensional dielectric random rough surfaces

The curvilinear coordinate method is applied for analysing 2-D dielectric random rough surfaces. The theory is based on Maxwell's equations written in a non-orthogonal coordinate system. For each medium, this method leads to an eigenvalue system. The scattered fields within two media are expanded as linear combinations of eigensolutions satisfying the outgoing wave condition. The boundary conditions allow the scattering amplitudes to be determined. The coherent and incoherent intensities are estimated by averaging the scattering amplitudes over several realizations. The theory is verified by comparison with results obtained by other exact method. A discussion on the C-method and the Sparse-Matrix CAnonical Grid method is proposed in terms of accuracy and computation time.     

Scattering of a plane wave from a periodic random surface: a probabilistic approach

This paper deals with a probabilistic formulation of the wave scattering from a periodic random surface. When a plane wave is incident on a random surface described by a periodic stationary stochastic process, it is shown by a group-theoretic consideration that the scattered wave may have a stochastic Floquet form, i.e. a product of a periodic stationary random function and an exponential phase factor. Such a periodic stationary random function is then written by a harmonic series representation similar to a Fourier series, where Fourier coefficients are mutually correlated stationary processes instead of constants. The mutually correlated stationary processes are represented by Wiener - Hermite functional series with unknown coefficient functions called Wiener kernels. In case of a slightly rough surface and TE wave incidence, low-order Wiener kernels are determined from the boundary condition. Several statistical properties of the scattering are calculated and illustrated in figures.     

两个相邻目标对平面波、高斯波束的光散射

基于等效原理和互易性定理研究了两个靠近目标对平面波、高斯波束的光散射问题,给出了这一复合光散射模型的二阶散射结果。通常一阶散射结果容易求解,但由于耦合效应的复杂性,很难给出二阶散射结果的解析形式。为了解决这一问题,应用互易性定理给出了求解任意相邻介质目标二阶散射场的公式,同时借助等效原理将求解散射场公式中的体积分简化为面积分的形式,从而降低了求解难度。求解了两相邻球形粒子的复合散射场,并将求解结果与应用时域积分方程法求得的结果进行了比较。同时,还讨论了束腰半径、目标位置对散射截面及偏振度的影响。     

Effect of the illumination length on the statistical distribution of the field scattered from one-dimensional random rough surfaces: analytical formulae derived from the small perturbation method

We consider a dielectric plane surface with a local cylindrical perturbation illuminated by a monochromatic plane wave. The perturbation is represented by a random function assuming values with a Gaussian probability density with zero mean value. Outside the perturbation zone, the scattered field can be represented by a superposition of a continuous spectrum of outgoing plane waves. The stationary phase method leads to the asymptotic field, the angular dependence of which is given by the scattering amplitudes of the propagating plane waves. The small perturbation method applied to the Rayleigh integral and the boundary conditions gives a first-order approximation of the scattering amplitudes. We show that the real part and the imaginary part of the scattering amplitudes are Gaussian stochastic variables with zero mean values and unequal variances. The values of variances depend on the length of the perturbation zone. In most cases, the probability density function for the amplitude is a Hoyt distribution and the phase is not uniformly distributed between -π and π. The standard Rayleigh and uniform distributions are obtained for special values of the length and in the case of an infinite illumination length.     

Effect of the illumination length on the statistical distribution of the field scattered from one-dimensional random rough surfaces: analytical formulae derived from the small perturbation method

We consider a dielectric plane surface with a local cylindrical perturbation illuminated by a monochromatic plane wave. The perturbation is represented by a random function assuming values with a Gaussian probability density with zero mean value. Outside the perturbation zone, the scattered field can be represented by a superposition of a continuous spectrum of outgoing plane waves. The stationary phase method leads to the asymptotic field, the angular dependence of which is given by the scattering amplitudes of the propagating plane waves. The small perturbation method applied to the Rayleigh integral and the boundary conditions gives a first-order approximation of the scattering amplitudes. We show that the real part and the imaginary part of the scattering amplitudes are Gaussian stochastic variables with zero mean values and unequal variances. The values of variances depend on the length of the perturbation zone. In most cases, the probability density function for the amplitude is a Hoyt distribution and the phase is not uniformly distributed between –π and π. The standard Rayleigh and uniform distributions are obtained for special values of the length and in the case of an infinite illumination length.     

Effective dynamic properties of 3D composite materials containing rigid penny-shaped inclusions

The propagation of time-harmonic plane elastic waves in infinite elastic composite materials consisting of linear elastic matrix and rigid penny-shaped inclusions is investigated in this paper. The inclusions are allowed to translate and rotate in the matrix. First, the three-dimensional (3D) wave scattering problem by a single inclusion is reduced to a system of boundary integral equations for the stress jumps across the inclusion surfaces. A boundary element method (BEM) is developed for solving the boundary integral equations numerically. Far-field scattering amplitudes and complex wavenumbers are computed by using the stress jumps. Then the solution of the single scattering problem is applied to estimate the effective dynamic parameters of the composite materials containing randomly distributed inclusions of dilute concentration. Numerical results for the attenuation coefficient and the effective velocity of longitudinal and transverse waves in infinite elastic composites containing parallel and randomly oriented rigid penny-shaped inclusions of equal size and equal mass are presented and discussed. The effects of the wave frequency, the inclusion mass, the inclusion density, and the inclusion orientation or the direction of the wave incidence on the attenuation coefficient and the effective wave velocities are analysed. The results presented in this paper are compared with the available analytical results in the low-frequency range.     

Investigation on electromagnetic scattering from rough soil surface of layered medium using the small perturbation method

Electromagnetic scattering from a rough surface of layered medium is investigated, and the formulae of the scattering coefficients for different polarizations are derived using the small perturbation method. A rough surface with exponential correlation function is presented for describing a rough soil surface of layered medium, the formula of its scattering coefficient is derived by considering the spectrum of the rough surface with exponential correlation function; the curves of the bistatic scattering coefficient of HH polarization with variation of the scattering angle are obtained by numerical calculation. The influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the roughness surface parameters and the frequency of the incident wave on the blstatic scattering coefficient is discussed. Numerical results show that the influence of the permittivity of layered medium, the mean layer thickness of intermediate medium, the rms and the correlation length of the rough surface, and the frequency of the incident wave on the bistatic scattering coefficient is very complex.     

具有二维fBm特征的分层介质粗糙面电磁散射的特性研究

运用微扰法研究了平面波入射分层介质粗糙面的电磁散射,推出了不同极化状态下的双站散射系数公式.采用二维fBm分形粗糙面来模拟实际的分层介质粗糙面,结合二维fBm分形粗糙面的功率谱导出了平面波入射二维fBm分形分层介质粗糙面的散射系数计算公式.通过数值计算得到了HH极化下双站散射系数随散射角的变化曲线,讨论了分维、底层介质介电常数、中间介质介电常数和厚度及入射波频率对双站散射系数的影响,得到了二维fBm分形分层介质粗糙面散射系数的分维特征、基本特征、分区特征和随频率变化的特征. 关键词： 电磁散射 二维fBm分形粗糙面 分层介质 微扰法     

On the “quantum rainbow” in surface scattering

In this letter we propose a new method to compute the scattering amplitudes of atoms from a corrugated wall model of a crystal surface. The method method requires reasonable computing time and converges rapidly to the solution of the scattering amplitudes (within 0.01% in each beam and satisfying unitarity) for any corrugation such that kζ₀ < 4.5 and ζ₀ < 0.35 Å, where k is the magnitude of the incident wave vector and ζ₀ the amplitude of the corrugation. These conditions for kζ₀ and ζ₀ are satisfied for any crystal surface and light scattering atoms (H and He) at thermal energies which makes the method appropriate for the system He, LiF for which results are presented comparing them with recent experimental data.     

An investigation of transmission coefficients for finite and semi-infinite coupled plate structures

This paper introduces a method for determining the transmission coefficient for finite coupled plates using an analytical waveguide model combined with a scattering matrix. In the scattering matrix method, the amplitudes of the structural waves impinging on a junction are separated into incident, reflected, and transmitted components. The energy flow due to each of these waves is obtained using a wave impedance method, which is subsequently used to determine the transmission coefficient. Transmission coefficients for semi-infinite and finite L-shaped plates are investigated for single and multiple point force excitations, and for controlled incident wave sources. It is shown that the transmission coefficients can also be calculated from details of the modal transmission coefficients and the modal composition of the energy incident on the junction. Results show that the modal transmission coefficients are largely independent of whether the plates have finite or semi-infinite boundary conditions, and are only dependent on the details of the coupling. Finally, frequency averaged transmission coefficients are compared for semi-infinite and finite structures. In the cases considered, it is found that the semi-infinite system is a good approximation for finite systems after frequency averaging, especially if the system is excited with multiple point force excitation.     

The use of the virtual source technique in computing scattering from periodic ocean surfaces

In this paper the virtual source technique is used to compute scattering of a plane wave from a periodic ocean surface. The virtual source technique is a method of imposing boundary conditions using virtual sources, with initially unknown complex amplitudes. These amplitudes are then determined by applying the boundary conditions. The fields due to these virtual sources are given by the environment Green's function. In principle, satisfying boundary conditions on an infinite surface requires an infinite number of sources. In this paper, the periodic nature of the surface is employed to populate a single period of the surface with virtual sources and m surface periods are added to obtain scattering from the entire surface. The use of an accelerated sum formula makes it possible to obtain a convergent sum with relatively small number of terms (～40). The accuracy of the technique is verified by comparing its results with those obtained using the integral equation technique.     

海面冰层对声波的反射和散射特性

北极海面冰层复杂多变,其对声波的反射和散射严重影响冰下水声信道的传输特性,建立海面冰层的声波反射和散射模型对冰下水声通信研究具有重要意义.假设海面冰层为多层固体弹性介质且冰-水界面粗糙,满足微扰边界条件,导出声波从海水介质入射到海面冰层时相干反射系数满足的线性方程组.对相干反射系数随声波频率、掠射角、冰层厚度的变化进行数值分析.进一步引入根据散射声场功率谱密度计算散射系数的方法,改变掠射角,对冰层厚度、散射掠角对散射系数的影响进行研究.     

起伏海面环境下水声信道特性及估计方法*

为了研究随机起伏海面的散射对水声信道特性的影响,以及不同信道估计算法在该信道条件下的性能问题,该文从利用Monte Carlo方法模拟的基于PM海浪谱的二维随机起伏海面出发,提取了与风向不同夹角的一维随机起伏海面,并利用Kirchhoff近似法计算了其散射强度。将得到的声散射特性和海面环境信息融入到信道模型中,建立了基于PM海浪谱的随机信道模型。通过仿真分析了不同风速、传播距离条件下,声波的传播损失和信道冲激响应变化,并利用l_0-最小均方误差法、匹配追踪算法和正交匹配追踪算法进行了信道估计。给出了不同环境条件下,三种算法估计性能的综合比较结果,验证了算法针对该信道模型的有效性。     

分形分层粗糙面电磁波透射特性的混合算法

运用基于矩量法(method of moment,MoM)结合基尔霍夫近似(Kirchhoff approximation,KA)的混合算法研究一维带限Weierstrass分形分层介质粗糙面的电磁波透射问题.在混合算法中将上层粗糙面和下层粗糙面分别划分到MoM区域和KA区域,数值计算得到透射系数随透射波的散射角的变化,讨论粗糙面分维数、高度起伏均方根、底层介质介电常数、中间层介质介电常数和厚度、入射波频率对透射系数的影响,结果表明上层粗糙面参数对透射系数有显著影响,而下层粗糙面参数对透射系数影响较小.     

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

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

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.     

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

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