Modification of Classical SPM for Slightly Rough Surface Scattering with Low Grazing Angle Incidence

Based on the impedance/admittance rough boundaries, the reflection coefficients and the scattering cross section with low grazing angle incidence are obtained for both VV and HH polarizations. The error of the classical perturbation method at grazing angle is overcome for the vertical polarization at a rough Neumann boundary of infinite extent. The derivation of the formulae and the numerical results show that the backscattering cross section depends on the grazing angle to the fourth power for both Neumann and Dirichlet boundary conditions with low grazing angle incidence. Our results can reduce to that of the classical small perturbation method by neglecting the Neumann and Dirichlet boundary conditions.     

Rough surface Green's function based on the first-order modified perturbation and smoothed diagram methods

This paper presents an analytical theory of rough surface Green's functions based on the extension of the diagram method of Bass, Fuks, and Ito with the smoothing approximation used by Watson and Keller. The method is a modification of the perturbation method and is applicable to rough surfaces with small RMS height. But the range of validity is considerably greater than for the conventional perturbation solutions. We consider one-dimensional rough surfaces with a Dirichlet boundary condition. The coherent Green's function is obtained from the smoothed Dyson's equation using a spatial Fourier transform. The mutual coherence function for the Green's function is obtained by first-order iteration of the smoothing approximation applied to the Bethe-Salpeter equation in terms of a quadruple Fourier transform. These integrals are evaluated by the saddle-point technique. The equivalent bistatic cross section per unit length of the surface is compared with that for the conventional perturbation method and the Watson-Keller result. With respect to the Watson-Keller result, it should be noted that our result is reciprocal, while the Watson-Keller result is non-reciprocal. Included in this paper is a discussion of the specific intensity at a given observation point. The theory developed will be useful for RCS signature related problems and low grazing angle scattering when both the transmitter and the object are close to the surface.     

Universal behaviour of scattering amplitudes for scattering from a plane in an average rough surface for small grazing angles

Abstract

It is shown that for scattering from a plane in an average rough surface, the scattering cross section of the range of small grazing angles of the scattered wave demonstrates a universal behaviour. If the angle of incidence is fixed (in general it should not be small), the diffusive component of the scattering cross section for the Dirichlet problem is proportional to θ² where θ is the (small) angle of elevation, and for the Neumann problem it does not depend on θ. For the backscattering case these dependences correspondingly become θ⁴ and θ°. The result is obtained from the structure of the equations that determine the scattering problem rather than by use of an approximation.     

Electromagnetic Scattering from One-Dimensional Time-Varying Sea Surface with Gaussian-Like Beam Incidence

A widely used iterative technique named the Method of Ordered Multiple Interactions is given for calculat-ing the Gaussian-like beam scattering from a time-varying sea surface with the Pierson-Moskowitz (P-M) spectrum. Thisis done by solving the magnetic field integral equation for the current induced on an infinite rough surface. Following thediscretization of the integral equation, the unknown currents can be determined more rapidly with the LU decomposition.Numerical results are presented with emphasis on the electromagnetic backscattering at low grazing angle incidence. Itis shown that the backscattering cross section is proportional to the nearly fourth power of the grazing angle for thebackscattering cross section with a different beam waist, surface length, and velocity of the wind are discussed     

Electromagnetic Scattering from One-Dimensional Time-Varying Sea Surface with Gaussian-Like Beam Incidence

A widely used iterative technique named the Method of Ordered Multiple Interactions is given for calculating the Gaussian-like beam scattering from a time-varying sea surface with the Pierson-Moskowitz (P-M) spectrum. This is done by solving the magnetic field integral equation for the current induced on an infinite rough surface. Following the discretization of the integral equation, the unknown currents can be determined more rapidly with the LU decomposition. Numerical results are presented with emphasis on the electromagnetic backscattering at low grazing angle incidence. It is shown that the backscattering cross section is proportional to the nearly fourth power of the grazing angle for the plane and beam incidence. This is consistent with the result given in some references. The angular distributions of the backscattering cross section with a different beam waist, surface length, and velocity of the wind are discussed.     

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.     

The intensity of the field generated by a point source and reflected from a rough surface

This paper studies the intensity of the acoustic field generated by a point source above a rough surface with the Dirichlet or Neumann boundary conditions. The derived equations are valid for arbitrary distances between the source, receiver and rough surface, including the case when these distances are smaller than the correlation radius of the surface roughness. It is believed that the proposed method is an improvement of the more conventional approach, which is based on integration over individual areas of the rough surface and that is valid when the source, receiver, and surface are at large distances from each other. The main limitation in deriving the expressions for the acoustic field intensity is the condition that the mutual shadowing of the surface points is small, which is close to the small slope approximation for the rough surface profile. The derivation includes the limiting cases which lead to the traditional small perturbation method and Kirchhoff approximation.     

Behaviour of scattering from a rough surface at small grazing angles

The particular problem of wave scattering at low grazing angles is of great interest because of its importance for the long-distance propagation of radio waves along the Earth's surface, radar observation of near surface objects, as well as solving many other fundamental and applied problems of remote sensing. One of the main questions is: how do the scattering amplitude and specific cross section behave for extremely small grazing angles? We consider the process of wave scattering by a statistically rough surface with the Neumann boundary condition. This model corresponds to sound scattering from a perfectly 'hard' surface (for example, the interface between air and the sea surface) or 'vertically' polarized electromagnetic waves scattered by a perfectly conducting one-dimensional (i.e. cylindrical) surface when the magnetic field vector is directed along the generating line of this cylindrical surface. We assume that the surface roughness is sufficiently small (in the sense of the Rayleigh parameter) and the surface is rigorously statistically homogeneous and therefore, infinite. We confine ourselves only to the first-order approximation of small perturbation theory and therefore consider every act of wave scattering in the Born approximation when the Bragg scattering process takes place. Only one resonant Fourier component of surface roughness is responsible for the scattering in a given direction. However, we take into account the attenuation of incident and scattered waves due to the multiple scattering processes on the path 'before' and 'after' a scattering event in a given direction. Also we consider every one of these multiple scattering events only in the Born approximation. The main result we have obtained is that for small grazing angles the scattering cross section of the diffuse component decreases as the second power of the grazing angles with respect to the incident and scattered directions, and as the fourth power of the grazing angle for the backscattering (radar) situation. Generalizing our results from plane-wave scattering to finite beams allows us to obtain the criterion on the beamwidth. For sufficiently narrow beams the multiple scattering processes do not play any role because of a short 'interaction path', and only single Bragg scattering determines the scattering amplitude (which does not tend to zero for small grazing angles). However, for sufficiently wide beams the result obtained for infinite plane waves becomes valid: due to the above-mentioned multiple scattering processes, the scattering amplitude tends to zero for small grazing angles. Consequently, the behaviour of the scattering cross section for small grazing angles depends on the radiation pattern width of the transmitting and receiving antennae: for sufficiently wide beams the scattering cross section decreases to zero at small grazing angles, but for narrow beams it tends to the finite non-zero value.     

A practical cross-section for scattering from rough surfaces at very low grazing angles in both the forward and backward directions

In this paper we develop an extension of the small slope approximation (SSA) for scattering from randomly rough Dirichlet surfaces, which includes some multiple scattering. This extension is designated by SSA+. We focus on scattering at very low grazing angles where multiple scattering of both the incident and scattered fields is of importance. Numerical results for the SSA+ bistatic scattering cross-section for very low forward grazing angles are presented using the Gaussian roughness spectrum and for both very low forward and very low backward grazing angles using the Pierson–Moskowitz and modified power law spectra. The results are restricted to an angle of incidence of 80°. It is shown that when the lowest-order SSA gives reasonably accurate results, the SSA+ increases the accuracy up to at least the final 0.2° of grazing in the forward direction. In the backward direction, the SSA+ gives good results for the Pierson–Moskowitz spectrum, but the results are less dramatic.     

Low-frequency scattering of acoustic waves by a bounded rough surface in a half-plane

The problem of the scattering of harmonic plane waves by a rough half-plane is studied here. The surface roughness is finite. The slope of the irregularity is taken as arbitrary. Two boundary conditions are considered, those of Dirichlet and Neumann. An asymptotic solution is obtained, when the wavelength lambda of the incident wave is much larger than the characteristic length of the roughness iota, by means of the method of matched asymptotic expansions in terms of the small parameter epsilon= 2piiota/lambda. For the Dirichlet problem, the solution of the near and far fields is obtained up to O(epsilon2). The far field solution is given in terms of a coefficient that have a simple explicit expression, which also appears in the corresponding solution to the Neumann problem, already solved. Also the scattering cross section is given by simple formulas to O(epsilon3). It is noted that, for the Dirichlet problem, the leading term is of order epsilon3 which, by contrast, is different from that of the circular cylinder in full space, that is, of order epsilon(-1) (log epsilon)(-2). Some examples display the simplicity of the general results based on conformal mapping, which involve arcs of circle, polygonal lines, surface cracks and the like.     

High-frequency wave diffraction by an impedance segment at oblique incidence

The plane problem of high-frequency acoustic wave diffraction by a segment with impedance boundary conditions is considered. The angle of incidence of waves is assumed to be small (oblique). The paper generalizes the method previously developed by the authors for an ideal segment (with Dirichlet or Neumann boundary conditions). An expression for the directional pattern of the scattered field is derived. The optical theorem is proved for the case of the parabolic equation. The surface wave amplitude is calculated, and the results are numerically verified by the integral equation method.     

求解非线性Poisson方程的格子演化算法

非线性Poisson方程在化学、化工及生物等领域有着广泛的应用。本文发展了一种基于格子演化的新算法-格子Poisson方法(LPM),并且给出了Dirichlet边界条件和Neumann边界条件的实现方法。本方法不需要对方程进行线化处理,直接求解非线性方程,适用范围广泛。Dirichlet边界与Neumann边界的数值模拟结果与多重网格法等结果符合很好,验证了该方法在求解非线性Poisson方程的正确性与有效性。本方法非常适合并行计算,并方便扩展到三维情况。     

Region of validity of perturbation theory for dielectrics and finite conductors

We present a study of region of validity of first-order perturbation theory applied to rough surface scattering. The scattering problem is solved numerically for the case of periodic surface or gratings varying in one dimension. Scattering of electromagnetic waves from an ensemble of gratings of sufficiently long period will give a good approximation to the case of an infinite rough surface. We use this to test the validity of the first-order perturbation theory. Use of an infinite periodic surface allows us to give results for a range of angle of incidence covering those representing a low grazing angle, near 90° from the mean surface normal. We consider the case for perfect dielectrics and finite conductors. The real and imaginary parts of the refractive index used were limited to less than three due to the numerical instability of the numerical calculation method involved. We find that for perfect dielectrics the first-order small perturbation theory remains for TE polarization valid for all incidence angles, while for TM polarization it seems to fail if the incidence angle approaches the Brewster angle.     

Studies of scattering theory using numerical methods

Abstract

Numerical simulations, using both exact and approximate methods, are used to study rough surface scattering in both the smd and large roughness regimes. This study is limited lo scattcring lrom rough one-dimensional surfaces that obey the Dirichlet boundary condition and have a Gaussian roughness spectrum. For surfdces with small roughness (kh?1, where k is the radiation wavenumber and h is the root-mean-square (RMS) Surface height), perturbation theory is known to be valid. However, it is shown numerically that when kh?1 and kl?6 (where I is the surface correlation length) the Kirchhoffapprorimation is valid except at low grazing angles, and one must sum the first three orders of perturbation theory obtain the correct result. For kh?1 and kl?1, first-order perturbation theory is accurate. In this region, the accuracy of the first two terms of the iterative series solution of the exact integral equation is examined; the first term a1 this series is the Kirchhoff approximation, It is shown numerically that lor very small kh these first two terms reduce to first-order perturbation theory. However, lor this reduction to occur, kh must be made smaller than necessdry lor first-order perturbation theory to be accurate. In the regime of large roughness (kh?1) backscattering enhancement occurs when the RMS slope is on the order of unity. Several investigators have recently shown that the second term of the iterative series solution (the double-scattering term) replicates the properties of backscattering enhancement reasonably well. However, the double-scattering term has a lundamental flaw: predictions lor the scattering cross section per unit length based on the double-scattering term increase as the surfdce length is increased. This is shown here with numerical simulations and with an approximate analytical result based on the high frequency limit. The physical significance of this finding is also discussed. The final topic is the use of the double-scattering approximation to study the mechanism for backscattering enhancement with the Dirichlet boundary condition. This mechanism is usually assumed to be interference between reciprocal scattering paths. When the interlerence between reciprocal scattering paths is removed, the enhancement is eliminated. This shows that interference between reciprocal paths is almost certainly the dominant mechanism for backscattering enhancement in the scattering regime studied.     

Scattering of bulk acoustic waves with different polarizations on a statistically rough free surface of a solid at oblique incidence

In the first Born approximation of the perturbation theory by a Green's function method developed by Maradudin, Mills [7] and Kosachev, Lokhov, Chukov [8,9] the problem of scattering bulk acoustic waves with different polarizations at oblique incidence on a statistically rough free boundary of an isotropic solid was solved. When the correlation function of the surface roughness is of a Gaussian form, the expressions for the transformation energy factor of the incident wave in the scattered volume and surface Rayleigh waves with respect to polarization, frequency and grazing angle of the incident wave as well as the roughness parameters and the Poisson coefficient of the medium were obtained. These results are helpful in accounting for the experiments on residual losses [15–17].     

相位微扰法在粗糙面光散射中的应用

根据粗糙面散射理论，用相位微扰法推导出了从随机粗糙介质表面散射的激光雷达散射截面的理论计算公式，计算了几种粗糙表面样品在１．０６μｍ下的单位面积激光雷达散射截面，数值结果与实验数据基本吻合。     

Path integral solution of linear second order partial differential equations I: the general construction

A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schrödinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette.     

Tracking gauge symmetry factorizability on intervals

We track the gauge symmetry breaking pattern by boundary conditions on fifth and higher-dimensional intervals. It is found that, with Dirichlet–Neumann boundary conditions, the Kaluza–Klein decomposition in five-dimension for arbitrary gauge group can always be factorized into that for separate subsets of at most two gauge symmetries, and so is completely solvable. Accordingly, we present a simple and systematic geometric method to unambiguously identify the gauge breaking/mixing content by general set of Dirichlet–Neumann boundary conditions. We then formulate a limit theorem on gauge symmetry factorizability to recapitulate this interesting feature. Albeit the breaking/mixing, a particularly simple check of orthogonality and normalization of fields' modes in effective 4-dim picture is explicitly obtained. An interesting chained-mixing of gauge symmetries in higher dimensions by Dirichlet–Neumann boundary conditions is also explicitly constructed. This study has direct applications to higgsless/GUT model building.     

Reliable approach for bistatic scattering of three-dimensional targets from underlying rough surface based on parabolic equation

A parabolic equation(PE) based method for analyzing composite scattering under an electromagnetic wave incidence at low grazing angle, which composes of three-dimensional(3-D) electrically large targets and rough surface, is presented and discussed. A superior high-order PE version is used to improve the accuracy at wider paraxial angles, and along with the alternating direction implicit(ADI) differential technique, the computational efficiency is further improved. The formula of bistatic normalized radar cross section is derived by definition and near-far field transformation. Numerical examples are given to show the validity and accuracy of the proposed approach, in which the results are compared with those of Kirchhoff approximation(KA) and moment of method(Mo M). Furthermore, the bistatic scattering properties of composite model in which the 3-D PEC targets on or above the two-dimensional Gaussian rough surfaces under the tapered wave incidence are analyzed.     

小攻角下三维细长体定常空化形态研究

基于边界元方法, 使空泡表面和细长体表面分别满足Dirichlet 边界条件和Neumann边界条件, 数值迭代获得小攻角下三维细长体的定常空化形态. 采用线性三角形单元, 将控制点布置在网格节点上, 应用局部正交坐标系并采用迭代方法获得空泡表面的速度势, 进而通过边界积分方程确定空泡厚度的分布. 采用拉格朗日插值方法得到空泡末端的厚度, 避免了迭代过程中网格的重新划分. 数值结果与实验值符合良好, 验证了该方法的合理性. 系统研究了空化数、攻角以及锥角对于三维细长体空化形态的影响规律. 数值结果表明: 攻角使得细长体的空化形态呈现不对称性, 出现空泡向背流面“堆积”的现象; 而空化数越小或锥角越大, 空泡形态的不对称性将越严重. 关键词： 边界元方法 三维细长体 局部空化 攻角