Hybrid method for investigation of electromagnetic scattering from conducting target above the randomly rough surface

A current based hybrid method (HM) is proposed which combines the method of moment (MOM) with the Kirchhoff approximation (KA) for the analysis of scattering interaction between a two-dimensional (2D) infinitely long conducting target with arbitrary cross section and a one-dimensional (1D) Gaussian rough surface. The electromagnetic scattering region in the HM is split into KA region and MOM region. The electric field integral equation (EFIE) in MOM region (target) is derived, the computational time of the HM depends mainly on the number of unknowns of the target. The bistatic scattering coefficient for the infinitely long cylinder above the rough surface with Gaussian roughness spectrum is calculated, and the numerical results are compared and verified with those obtained by the conventional MOM, which shows the high efficiency of the HM. Finally, the influence of the size, location of the target, the rms height and correlation length of the rough surface on the bistatic scattering coefficient with different polarizations is discussed in detail.     

任意截面介质波导导模的叠代矩量法分析

本文提出用以计算任意截面介质波导色散关系的叠代矩量法,并给出用叠代矩量法计算椭圆介质波导、三角形波导、平切圆波导、蛋形截面波导等的传播特性的实例.计算结果与有限元法、点匹配法等精确数值计算结果符合很好.本方法简便易行,计算工作量较小.     

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

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

Solution of problem of electromagnetic waves scattering on inhomogeneously layered scatterers using pattern equation method

In this paper, the pattern equation method (PEM) is used to solve the problems of electromagnetic waves scattering by inhomogeneously layered scatterers. The numerical algorithm based on PEM for multilayered scatterers of arbitrary geometry is developed. The numerical results are presented for magneto-dielectric objects coated with magneto-dielectric materials to demonstrate the stability, accuracy, and efficiency of the proposed method. A comparison of the results obtained by the present method and by the other methods are shown and the limits of applicability of method are outlined.     

Investigation of composite electromagnetic scattering from ship-like target on the randomly rough sea surface using FDTD method

Composite electromagnetic scattering from a two-dimensional (2D) ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain (FDTD) method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices. The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique, a 2D infinitely long cylinder over the sea surface is taken into account first. The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally, the influence of the incident angle, the polarization, and the size of the ship-like target on the composite scattering coefficient is discussed in detail.     

Analysis of the Scattering Characteristic of Several Waveguide Components Using Boundary Element Method

Several kinds of waveguide components such as curved waveguide bends, arbitrary angle waveguide bends and T-junctions have been analyzed with boundary element method in this paper. A new discretization method for the boundary element method to solve the waveguide discontinuities has been given. The numerical results obtained agree well with the experimental results and numerical results in other literature. Especially, the scattering characteristics of Forded E-, H-plane T-junctions in 3mm band have been analyzed using boundary element method and the calculation results are presented.     

Improved Analytical Approximations to the Scattering Solutions of the Schrödinger Equation with a Hyperbolical Potential

Using an improved approximate formula to the centrifugal term, we present arbitrary l-state scattering solutions of the hyperbolic potential. The approximate analytical formula of scattering phase shifts and normalized wavefunctions are presented. All data calculated by the above approximate analytical formula are compared with those obtained by using the numerical integration method in the scattering state cases. We find that this improved approximate formula is better than previous one since the calculated results are in good agreement with those exact ones.     

有限元/边界积分方法在海面及其上方弹体目标电磁散射中的应用

本文将有限元/边界积分方法(FE/BIM)结合区域分解方法引入到粗糙海面及其上方目标 的电磁散射问题的研究中. 由于积分边界可以以任意形状设置在距模型表面任意远的距离处, 故本文采用共形人工边界结合区域分解建模方法截断模型的开放计算区域以减少求解未知量, 在截断区域内部采用有限元方法求解, 而计算区域的边界条件通过边界积分方程方法得到. 通过与矩量法获得的数值计算结果进行比较, 证明了该混合算法及模型处理方法的正确性, 进而研究了海面上方弹体目标的电磁散射特性, 并讨论了其双站散射系数随电磁波入射角度、目标高度、海面风速以及弹体尺寸的电磁散射特性变化情况. 本文结果可用于反演复杂背景下的目标信息及目标探测等领域. 关键词： 电磁散射 粗糙海面 目标 有限元/边界积分方法     

Electromagnetic scattering from two-layer rough interfaces in the Kirchhoff approximation

Electromagnetic wave scattering from multilayers consisting of two two-layer Gaussian rough surfaces with lossless media is investigated in the Kirchhoff approximation (KA), with consideration of the shadowing effects. The tapered incident wave is introduced into the classic KA, and the bistatic scattering coefficient is redetermined. The advantage of this method is that it is faster in computation than the exact numerical methods. The numerical results show that the bistatic scattering coefficient calculated in the KA is in good agreement with that obtained by using the method of moment (MOM) over a most angular range, which indicates the validity of the KA proposed in our paper. Finally, the effects of the relative permittivity, the root-mean-square (RMS) height, the correlative length, and the average height between the two interfaces on the bistatic scattering coefficient are discussed in detail.     

Higher-order quadrature-based moment methods for kinetic equations

Kinetic equations containing terms for spatial transport, body forces, and particle–particle collisions occur in many applications (e.g., rarefied gases, dilute granular gases, fluid-particle flows). The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent variables. A useful alternative is to reformulate the problem in terms of the moments of the velocity distribution function. Closure of the moment equations is challenging for flows sufficiently far away from the Maxwellian limit. In previous work, a quadrature-based third-order moment closure was derived for approximating solutions to the kinetic equation for arbitrary Knudsen number. A key component of quadrature-based closures is the moment-inversion algorithm used to find the non-negative weights and velocity abscissas. Here, a robust inversion procedure is proposed for three-component velocity moments up to ninth order. By reconstructing the velocity distribution function, the spatial fluxes in the moment equations are treated using a kinetic-based finite-volume solver. Because the quadrature-based moment method employs the moment transport equations directly instead of a discretized form of the kinetic equation, the mass, momentum and energy are conserved for arbitrary Knudsen and Mach numbers. The computational algorithm is tested for the Riemann shock problem and, for increasing Knudsen numbers (i.e. larger deviations from the Maxwellian limit), the accuracy of the moment closure is shown to be determined by the discrete representation of the spatial fluxes.     

Electromagnetic modeling of inhomogeneous and anisotropic structures by volume integral equation methods

Given the application of inhomogeneous and anisotropic structures in different application areas, it is of critical importance to develop accurate and efficient modeling methods. Among various methods, volume integral equations (VIEs) using moment method are efficient solutions for electromagnetic modeling of inhomogeneous and anisotropic structures. In this paper, we investigate the solutions of the VIE method and augmented volume integral equation (A-VIE) method for solving inhomogeneous and anisotropic structures with arbitrary shapes. The moment method solutions are presented and curl-conforming bases are used to discretize the electric and magnetic field distributions inside the structures. When the structures contain inhomogeneous magnetic materials, A-VIE method is applied. Various numerical examples are shown to demonstrate the accuracy and efficiency of the algorithms.     

粗糙海面及其上方导体目标复合电磁散射的混合算法研究

利用基于电流计算的矩量法结合高频算法基尔霍夫近似的混合算法,分析了一维PM谱粗糙海面及其上方二维无限长任意截面导体目标的双站复合电磁散射特性.混合算法将粗糙面和目标分别划分到KA区域和MOM区域,由于无需数值求解粗糙海面区域的表面极化电流,该算法的运算时间和对计算机内存的需求主要取决于粗糙面上方目标的网格划分情况.数值结果以无限长导体圆柱为例计算了其与一维下垫PM谱粗糙海面的复合双站散射截面,并将计算结果与经典MOM结果进行了比对和验证,结果表明混合方法具有较高的计算效率.最后应用混合方法讨论了不同极化状态、海上不同风速以及目标不同尺寸和位置对复合散射截面的影响. 关键词： 粗糙海面 电磁散射 混合算法 矩量法     

Improved Analytical Approximations to the Scattering Solutions of the Schrdinger Equation with a Hyperbolical Potential

Using an improved approximate formula to the centrifugal term, we present arbitrary l-state scattering solutions of the hyperbolic potential. The approximate analytical formula of scattering phase shifts and normalized wavefunctions are presented. All data calculated by the above approximate analytical formula are compared with those obtained by using the numerical integration method in the scattering state cases. We find that this improved approximate formula is better than previous one since the calculated results are in good agreement with those exact ones.     

A review of the numerical investigation on the scattering of Gaussian beam by complex particles

The study of light scattering by various particles is an active and important subject of research with myriad practical applications. During the years the scattering of plane wave by various particles has been investigated extensively. In recent years, with the development of laser sources and the tremendous expansion of their application, there has been a growing interest in the study of light scattering by various particles illuminated by a focused Gaussian beam. Since the analytical methods are only suitable for the analysis of Gaussian beam scattering by some regular particles, for complex particles with arbitrary shape and structure, one has to resort to the numerical methods. In this article, we review the recent numerical investigation on the scattering of Gaussian beam by systems of complex particles, including arbitrarily shaped conducting particles, dielectric particles, composite particles with inclusions, as well as random discrete particles and fractal soot aggregates. The essential formulations of the proposed numerical methods are outlined and the numerical results for some complex particles are also presented. This review is expected to provide useful help for the study of the interaction between the laser beams and the complex particles.     

浅海波导中目标散射的边界元方法

提出了一种计算三维散射体在声速剖面随深度变化、距离无关浅海波导中散射声场的数值方法波导边界元方法。当散射体不十分靠近波导界面因而边界多次散射可以忽略时,在边界元计算中可以用自由场格林函数近似波导格林函数。应用镜像法和球波函数加法定理推导了理想波导中球体散射声场的解析解,用来验证波导边界元方法的计算精度,证明该数值方法是准确的。对浅海波导中水下潜器散射声场数值模拟的结果表明,浅海波导海面、海底界面反射、声速剖面等对目标散射声场的幅值和方向性都有很大的影响。     

Microwave Doppler spectra of sea return at small incidence angles: specular point scattering contribution

It is well known that the sea return echo contains contributions from at least two scattering mechanisms. In addition to the resonant Bragg scattering, the specular point scattering plays an important role as the incidence angle becomes smaller (≤20o). Here, in combination with the Kirchhoff integral equation of scattering field and the stationary phase approximation, analytical expressions for Doppler shift and spectral bandwidth of specular point scattering, which are insensitive to the polarization state, are derived theoretically. For comparison, the simulated results related to the two-scale method (TSM) and the method of moment (MOM) are also presented. It is found that the Doppler shift and the spectral bandwidth given by TSM are insufficient at small incidence angles. However, a comparison between the analytical results and the numerical simulations by MOM in the backscatter configuration shows that our proposed formulas are valid for the specular point scattering case. In this work, the dependences of the predicted results on incidence angle, radar frequency, and wind speed are also discussed. The obtained conclusions seem promising for a better understanding of the Doppler spectra of the specular point scattering fields from time-varying sea surfaces.     

Study of femtosecond soliton dynamics in photonic crystal fiber using the moment method

We investigate the dynamics of femtosecond solitons in photonic crystal fibers (PCFs) by including high-order dispersion terms until to sixth-order in the generalized nonlinear Schrödinger equation, in addition to the nonlinear effects of the self phase modulation, self steepening and Raman scattering. We calculate theoretically the pulse parameters using the moment method. In the case of the fundamental soliton, our computed equations are coupled and difficult to solve analytically. However, we use the finite difference method to calculate numerically pulse parameters using an initially hyperbolic secant pulse at 1550-nm with different peak powers along 10m-PCF. Our numerical results show that the nonlinear regimes allow obtaining pulse compressions and initial pulse amplitudes. Furthermore, we remark a pulse broadening, and weak shifts of the peak power positions and frequencies in the critical and dispersive regimes. The use of an initial chirp provides a better pulse compressions and especially for low input powers. Also, the initial positive chirp reduces the optimal compression position lengths, while the negative one increases them. Therefore, we conclude that our theoretical calculations and numerical simulation results show that the moment method associated with the finite differences method is effective for the study of femtosecond pulse dynamics in PCFs.     

Scattering states of modified PSschl-Teller potential in D-dimension

We present a new approximation scheme for the centrifugal term,and apply this new approach to the Schrdinger equation with the modified Pschl-Teller potential in the D-dimension for arbitrary angular momentum states.The approximate analytical solutions of the scattering states are derived.The normalized wave functions expressed in terms of the hypergeometric functions of the scattering states on the k/2π scale and the calculation formula of the phase shifts are given.The numerical results show that our results are in good agreement with those obtained by using the amplitude-phase method(APM).     

Transient response in flexure to general uni-directional loads of variable cross-section beam with concentrated tip inertias immersed in a fluid

This paper presents a method for solving problems of transient response in flexure due to general unidirectional dynamic loads of beams of variable cross section with tip inertias. An elastodynamic theory which includes effects of continuous mass and rigidity of the beam has been applied. In the analysis the general dynamic load is expanded into a Fourier series and the beam is divided into many small uniform thickness segments. The equation of motion of each segment is mapped onto the complex domain by use of the Laplace transform method. The solutions of each set of adjoining segments are related to each other at the boundaries by the use of the transfer matrix method. The displacement, the bending slope, the bending moment and the shearing force at each boundary and at arbitrary time are obtained from the Laplace transform inversion integral by using the residue theorem. The theoretical results given in this paper are applicable to problems of dynamic response due to arbitrary loads varying with time of beams of arbitrary shape with concentrated tip inertias. As applications of the present theoretical results, numerical calculations have been carried out for two cases: a uniform beam with a tip inertia and a non-uniform beam (a truncated cone) with a tip inertia. Both are immersed in a fluid and subjected to large waves such as cnoidal waves.