应用改进的物理光学法和图形计算电磁学近似算法快速计算导体目标电磁散射特性

结合改进的物理光学法和图形计算电磁学法实现了考虑边缘绕射情况下复杂目标的高频电磁波散射的高效且精确求解.传统的考虑边缘绕射的物理光学算法不能直接计算出目标的雷达截面,它需要先计算绕射贡献,然后加上物理光学的散射贡献,最终才能得到目标的雷达截面.通过运用改进的物理光学法对图形计算电磁学法进行修正,直接修正表面法向量,从而修正了表面电流,这样就考虑了边缘处的绕射,提高了算法的效率.这不但充分利用了计算机硬件优势,借助于计算机显示技术,由图形加速卡完成最困难、最费时的消隐工作,而且利用图形计算电磁学的积分公式,将三维空间的积分转化为屏幕像素的二维空间积分,使得计算大幅简化.数值结果表明了所提出方法的精确性和高效性.     

任意形状凸粗糙物体高斯光束相干散射研究

研究了任意形状凸粗糙物体对高斯光束的相干散射特性。由平面波谱展开法推导出粗糙面高斯波束散射场表达式．并根据物理光学近似和稳相法原理得到相干散射截面的理论公式。与平面波解不同的是,在公式中引入了波束因子这一表征波束作用的重要参量,该参量与束腰半径、束腰中心与物体中心距离和物体照射面积以及入射、散射方向有关。最后数值计算了1．06μm激光波束对粗糙球的散射,分析了波束因子．介电常量和极化状态对红外激光相干散射截面的影响,重点讨论了波束因子的变化规律。分析表明．当波束半径远大于物体尺寸时,波束入射的结果可退化到平面波情况。     

Geometrical optics approximation for light scattering by absorbing spherical particles

By means of geometrical optics, an approximation method is presented to compute the light scattering intensity of absorbing spherical particles illuminated by a plane wave. For absorbing particles, the effective refractive index and the effective refractive angle are related to the complex refractive index and incident angle. The formulas for calculation of the break of phases of reflection and refraction, which are different from the case of transparent particles, are exactly derived. Verification of the geometrical optics approximation (GOA) was performed by case studies and comparison of the present results with the Mie scattering. It is found that agreement between the GOA and the Mie theory is excellent in forward directions for weakly/moderately absorbing particles. Differently, for strongly absorbing particles, good agreement between the calculation methods is in the forward directions and large scattering angles. The agreement between the GOA and the Mie theory is better for larger particles.     

水中气泡上的体散射函数的模拟与计算

基于几何光学的基本原理,推导出一种可以计算水中大尺度气泡上光散射角度与强度的关系式。推导中避免引入衰减因子G,较Davis模型更为简单。该模型可应用于光在水中单个气泡上散射的数值计算。最后,模拟计算了平行光束入射水中气泡的体散射函数曲线,发现水中气泡的前向散射远大于后向散射;当气泡半径在远大于光波波长的前提下变化时,气泡上散射光强分布规律与气泡半径无关;而介质相对折射率的增大会削弱前向散射而增强后向散射光强。     

折射率正负梯度交替表面的研究

利用不均匀材料提出了一种折射率正负梯度交替表面, 并利用几何光学法和电磁场数值仿真方法对其机理进行了研究与验证, 发现其可以操控电磁波的传播: 在一定的参数下, 其可以将空间中的入射电磁波一直束缚在介质中, 或在介质中传播一定距离后再实现电磁波的释放; 同时该结构具有宽带、极化无关等特性. 利用该结构可引导能流, 减小后向散射截面, 可用作隐身表面.     

Bubble acoustical scattering cross section under multi-frequency acoustic excitation

The acoustical scattering cross section is usually employed to evaluate the scattering ability of the bubbles when they are excited by the incident acoustic waves. This parameter is strongly related to many important applications of performance prediction for search sonar or underwater telemetry, acoustical oceanography, acoustic cavitation, volcanology, and medical and industrial ultrasound. In the present paper, both the analytical and numerical analysis results of the acoustical scattering cross section of a single bubble under multi-frequency excitation are obtained. The nonlinear characteristics(e.g.,harmonics, subharmonics, and ultraharmonics) of the scattering cross section curve under multi-frequency excitation are investigated compared with single-frequency excitation. The influence of several paramount parameters(e.g., bubble equilibrium radius, acoustic pressure amplitude, and acoustic frequencies) in the multi-frequency system on the predictions of scattering cross section is discussed. It is shown that the combination resonances become significant in the multi-frequency system when the acoustic power is big enough, and the acoustical scattering cross section is promoted significantly within a much broader range of bubble sizes and acoustic frequencies due to the generation of more resonances.     

Directional hemispherical radiative properties of random dielectric rough surfaces

In this paper the directional hemispherical reflectivity and transmissivity of one-dimensional, randomly rough, dielectric surfaces are determined by the use of the integral method. This method is derived from electromagnetic theory without any restrictive hypotheses. Since this exact approach is computationally very intensive, a geometric optics approximation method is also developed. Curves displaying radiative properties versus the correlation length for a constant mean square deviation of the surface from flatness are presented. In this respect, the influence on the validity of the approximate method of multiple scattering, the shadowing effect and the real index of refraction of the dielectric have been investigated. Transverse electric and transverse magnetic polarized incident plane waves are considered. For the latter, our interest is focused on the influence of roughness on the reflected and transmitted intensities for an angle of incidence close to the Brewster angle.     

Directional hemispherical radiative properties of random dielectric rough surfaces

In this paper the directional hemispherical reflectivity and transmissivity of one-dimensional, randomly rough, dielectric surfaces are determined by the use of the integral method. This method is derived from electromagnetic theory without any restrictive hypotheses. Since this exact approach is computationally very intensive, a geometric optics approximation method is also developed. Curves displaying radiative properties versus the correlation length for a constant mean square deviation of the surface from flatness are presented. In this respect, the influence on the validity of the approximate method of multiple scattering, the shadowing effect and the real index of refraction of the dielectric have been investigated. Transverse electric and transverse magnetic polarized incident plane waves are considered. For the latter, our interest is focused on the influence of roughness on the reflected and transmitted intensities for an angle of incidence close to the Brewster angle.     

Scattering characteristics of nonmagnetic invisibility cloak with minimized scattering

A kind of transformation functions is proposed to realize the nonmagnetic invisibility cloak with minimized scattering on the basis of generalized transformation. By matching the impedance at the outer surface of the cloak, the transformations with two parameters are determined. The good performance of the cloak is indicated by the full wave simulation based on the finite element method. Furthermore, based on the calculation of total scattering cross section, it is shown that the scattering cross section is very sensitive to the different parameters even though the impedance at the exterior boundary matches perfectly with the free space. In addition, from the effective media theory, an alternating layered system composed of two isotropic materials is proposed to realize experimentally the cloak.     

Scattering and absorption of light by ice particles: Solution by a new physical-geometric optics hybrid method

A new physical-geometric optics hybrid (PGOH) method is developed to compute the scattering and absorption properties of ice particles. This method is suitable for studying the optical properties of ice particles with arbitrary orientations, complex refractive indices (i.e., particles with significant absorption), and size parameters (proportional to the ratio of particle size to incident wavelength) larger than ∼20, and includes consideration of the edge effects necessary for accurate determination of the extinction and absorption efficiencies. Light beams with polygon-shaped cross sections propagate within a particle and are traced by using a beam-splitting technique. The electric field associated with a beam is calculated using a beam-tracing process in which the amplitude and phase variations over the wavefront of the localized wave associated with the beam are considered analytically. The geometric-optics near field for each ray is obtained, and the single-scattering properties of particles are calculated from electromagnetic integral equations. The present method does not assume additional physical simplifications and approximations, except for geometric optics principles, and may be regarded as a “benchmark” within the framework of the geometric optics approach. The computational time is on the order of seconds for a single-orientation simulation and is essentially independent of the size parameter. The single-scattering properties of oriented hexagonal ice particles (ice plates and hexagons) are presented. The numerical results are compared with those computed from the discrete-dipole-approximation (DDA) method.     

The effect of a layer of varying density on spatial correlation of sea-bottom backscattering signal

The research is about the effect of a layer of varying density of sea-bottom sediments on spatial correlation of sea-bottom backscattering.The relationship between scattering cross section and spatial correlation is that backscattering cross section decreases quickly and the spatial correlation becomes stronger as the incident angle increases.Therefore,the densitydepth profile is introduced into sea-bottom high-frequency backscattering echo model,which is used to simulate sea-bottom backscattering and calculate the function of spatial correlation.The influence of the density gradient on spatial correlation of sea-bottom backscattering is investigated by analyzing the relations between vertical gradient of density and the scattering cross section.As can be seen from the simulation results,the impact of the density gradient on the spatial correlation is found more significant.While the density gradient increases,the scattering cross-section and the radius of the spatial correlation broaden,the spatial correlation becomes stronger.At the same time,the scattering cross-section decreases more quickly as the incident angle increases.     

扩展性微动目标回波模拟与特征参数提取研究

研究了微动目标的多普勒回波模拟及特征参数提取技术. 提出了一种基于物理光学法和等效电磁流法的扩展性微动目标回波模拟方法. 将在目标坐标系下计算得到的后向散射场通过坐标转换, 成为雷达坐标系下的目标回波, 通过与解析信号模型对比验证方法的正确性. 分析了圆锥与带翼弹头的进动特性, 为获得较好的时频聚集性同时避免交叉项采用S-method方法对获取的回波信号进行时频分析, 分析了不同雷达波入射角度, 不同运动状态及不同几何外形的时频分布特点. 对时频分布图进行逆Radon变换, 将正弦曲线映射到参数空间, 从而获取目标的微动参数. 该研究结合电磁散射与信号处理技术, 通过对典型弹道目标的仿真, 获得一些不同于传统微动模型的结果, 结合电磁散射理论, 对这些现象进行了解释分析. 该研究成果在弹道目标的探测识别领域具有重要的理论与应用价值.     

Coupling to the elastic channel in electron impact spectroscopy: A simple second born model and its application to excitation of the 61P1 state of mercury and to the excitation of molecules

The matrix element in the infinite channel close coupling approximation responsible for coupling to the elastic channel in electron impact inelastic encounters is investigated. The contribution from the imaginary part of the energy denominator in the elastic coupling matrix element for dipole allowed transitions is shown to yield large angle differential cross sections in good agreement with experiment. This coupling mechanism predicts that the shape of the inelastic differential cross section will be dominated by the shape of the elastic cross section in the large angle high energy limit. In fact the coupling matrix element exhibits a dependence on incident energy, k², and momentum transfer, K, of the form 1/kK² which is in agreement with the theoretical predictions of Huo and means that in the limit of high incident energy the non-first-Born elastic coupling will dominate the angular dependence of the inelastic differential cross section at large scattering angles. In the case of molecular electron impact spectra it is shown that the analog of the Massey—Moore coherence features depending on the symmetry of the states involved in the excitation process will also occur in the coupling contribution. It is suggested that this new mechanism for producing coherent features in inelastic differential cross sections may be the explanation of the coherent features observed experimentally by Karle and Swick.It can be concluded on the basis of the results obtained here that the coupling to the elastic cross section provided by the imaginary contribution from the second Born energy denominator is sufficient to explain presently available experimental data on the large angle differential cross section and spin polarization. The simple coupling model was found to be inadequate to explain the small angle differential cross section in the range 10° < θ < 30° even at incident energies as high as 400 eV. The calculations also showed significant differences between first and second Born calculations at zero scattering angle. No conclusion can be drawn about this observation as all the omitted terms should make significant contributions in the small angle region. It is important to again emphasize that the large angle scattering even in the limit of high incident electron energy will be completely dominated by the coupling to the elastic channel^{7, 11}. On the basis of this work it appears that the coherent structure in the large angle inelastic differential cross section observed by Swick and Karle^{12, 13} at incident electron energies in the keV region may well be due to coupling to the elastic channel.     

海底密度梯度对底回波信号空间相关特性的影响研究

研究了海底表层沉积物的垂直密度梯度对底回波空间相关特性强度的影响。底回波散射截面与空间相关特性间的关系为:随入射角的增大,回波散射截面下降得越快,底回波空间相关特性越强。因此,将密度-深度模型引入海底高频回波模型进行回波仿真,计算回波的空间相关函数,通过分析密度梯度对底回波散射截面大小的影响,考察其对底回波空间相关特性强度的影响。从仿真分析结果可以看出,海底表层沉积物密度梯度对空间相关特性的影响较为显著,随着密度梯度的增大,回波散射截面变大且随入射角的增大下降趋势变快,空间相关半径变长,空间相关特性变强。      

Difference scattering field properties between periodic defect particles and three-dimensional slightly rough optical surface

Based on the practical situation of nondestructive examination, the calculation model of the composite scattering is established by using a three-dimensional half-space finite difference time domain, and the Monte Carlo method is used to solve the problem of the optical surface with roughness in the proposed scheme. Moreover, the defect particles are observed as periodic particles for a more complex situation. In order to obtain the scattering contribution of defects inside the optical surface, a difference radar cross section is added into the model to analyze the selected calculations on the effects of numbers, separation distances, different depths and different materials of defects. The effects of different incident angles are also discussed. The numerical results are analyzed in detail to demonstrate the best position to find the defects in the optical surface by detecting in steps of a fixed degree for the incident angle.     

Investigation on transient electromagnetic scattering from a randomly rough surface and the perfect electric conductor target with an arbitrary cross section above it

The Time-Domain-Integral-Equation (TDIE) method is proposed to analyze transient scattering interaction between a two-dimensional infinitely long conducting target with an arbitrary cross section and a one-dimensional rough surface. Based on the electric-field-integral-equation in time domain, the explicit and implicit solutions of MOT (Marching-on-time) are derived and presented. The current response at the center of the rough surface and the far electric field response with time in the composite model are calculated and analyzed. The numerical results are compared and verified with those obtained by conventional MOM-IDFT (Method of Moment-inverse discrete Fourier transform). Finally, the influence of the size, the location of the target and the incident angle on the current response and the far electric fields response are discussed in detail. Supported by the National Natural Science Foundation of China (Grant No. 60571058) and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20070701010)     

基片及其上方回转椭球粒子极化光散射

基于BV理论建立基片及其上方回转椭球粒子的复合散射模型,通过矢量球谐函数展开,对散射过程进行了分析,对散射场及微分散射截面详细求解,并给出了数值计算结果,与离散源方法做了比较,同时退化为球粒子与扩展Mie理论做了比较,说明了此方法的有效性。并详细讨论分析了微分散射截面随不同入射角,散射角,回转椭球粒子的尺寸、长短轴比例,距基片的距离,介电常数,粒子取向角的变化关系。结果表明:同一散射角下入射角越大,其微分散射截面越大;粒子尺寸越大,相互作用越大,其微分散射截面越大;长短轴比例越大,其微分散射截面越小;距离基片的距离越大,微分散射截面越大;微分散射截面的变化主要依赖于相对介电常数实部、虚部数值较大的一方,并且随粒子取向角的增大而增大。     

Light atom scattering from adsorbates at low coverage

A multiple scattering theory is developed for the scattering of light atoms from a disordered adsorbate on a smooth surface. Using gas phase potentials for adatom potentials, excellent agreement is found with data for He scattering from adsorbed Xe and surprisingly good agreement with the extensive data for He scattering from adsorbed CO. The attractive adatom potential must be included in the calculation if quantitative comparison with experimental data is to be made. An important contribution to the total cross section of an adatom are quantum mechanical oscillations similar to glory oscillations in gas phase scattering. The dependence of the total cross section on the incident angle is strongly influenced by the varying number of adatoms seen by the probe.     

Factors influencing electromagnetic scattering from the dielectric periodic surface

The scattering characteristics of the periodic surface of infinite and finite media are investigated in detail.The Fourier expression of the scattering field of the periodic surface is obtained in terms of Huygens' s principle and Floquet's theorem.Using the extended boundary condition method(EBCM) and T-matrix method, the scattering amplitude factor is solved,and the correctness of the algorithm is verified by use of the law of conservation of energy.The scattering cross section of the periodic surface in the infinitely long region is derived by improving the scattering cross section of the finite period surface.Furthermore, the effects of the incident wave parameters and the geometric structure parameters on the scattering of the periodic surface are analyzed and discussed.By reasonable approximation, the scattering calculation methods of infinite and finite long surfaces are unified.Besides, numerical results show that the dielectric constant of the periodic dielectric surface has a significant effect on the scattering rate and transmittance.The period and amplitude of the surface determine the number of scattering intensity peaks, and, together with the incident angle, influence the scattering intensity distribution.     

Analysis of multiple scattering for (e, 2e) experiments on thin films

The coincidence rate for (e, 2e) scattering from thin films has significant contributions from multiple collisions in contrast to (e, 2e) experiments on gas targets where the multiple scattering rate is negligible. The most likely multiple scattering events involve one or more small angle collisions in addition to the wide angle (e, 2e) scattering. The total coincidence rate can be written as the convolution of a “smearing function” with the (e, 2e) cross section. The smearing function is an infinite series of multiple order, small-angle scattering events and can be determined from the measured, small-angle cross section of the particular material. The (e, 2e) cross section can be deconvoluted from the smearing function by standard techniques. A general expression for the smearing function is derived in this paper and is evaluated for the specific case of amorphous carbon. As an example in applying this analysis, (e, 2e) data on amorphous carbon are deconvoluted from the smearing function by two different methods.