Fast iterative approach to the difference scattering from a dielectric target above a rough surface

1 Introduction Electromagnetic scattering from the target above or beneath a rough surface has at- tracted much interest during recent years, because of extensive applications to radar surveillance, target detection, ground radar probing, and so on[1―3]. In order to numeri- cally simulate scattering from composite model of the target and underlying rough sur- face, some fast numerical methods, such as general forward backward method and spec- trum acceleration algorithm (GFBM/SAA)[4,5], fin…     

随机Gauss粗糙面上三维导体目标散射差场的随机泛函解析计算方法

提出一种解析的随机泛函方法(SFA),计算导体Gauss粗糙面上三维导体目标的复合电磁散射.推导粗糙面的随机Green函数,用一种新的四路径模型描述面体复合散射机理,用SFA求解双站差场雷达散射截面.以导体球目标为算例,与其他数值计算方法比较后验证了SFA的有效性与准确性,同时讨论了粗糙度、体目标尺寸以及距离粗糙面高度等参量变化对结果的影响,给出复杂形状体目标的双站差场雷达散射截面的空间角分布. 关键词： 随机泛函方法 粗糙面随机Green函数 差场雷达散射截面 面体复合散射     

三维随机粗糙面与目标复合电磁散射的FDTD方法

用时域有限差分方法(FDTD)研究三维周期性延拓的随机粗糙面与上方目标复合电磁散射.用周期性延拓消除数值计算中截取有限大小粗糙面产生的边缘效应,讨论一个周期单元粗糙面的边长与其相关长度之间的关系.给出在FDTD方法中向粗糙面加载入射波的方法,建立了粗糙面上单个三维目标的复合散射FDTD计算模型.数值结果给出粗糙面与目标散射的近场分布,应用近远场变换得到全方位散射角的双站散射系数.比较了三维与二维散射模型的区别.结果显示当粗糙面上放置目标时,其后向散射显著增强.     

三维随机粗糙面上导体目标散射的解析-数值混合算法

提出三维导体目标与导体粗糙面复合散射的解析-数值混合迭代算法,推导出三维目标与粗糙面的耦合积分方程,以及粗糙面散射的Kirchhoff近似(KA)计算式.粗糙面的KA解析计算大大降低了粗糙面求解的复杂度,与目标矩量法的混合迭代保证了计算结果的精度,使得三维体-面目标复合散射计算变得可行.由于体-面两者的高阶耦合作用明显减小,保证了该混合迭代算法的收敛性.与镜像Green函数方法的比较表明该混合算法的有效性,并讨论了粗糙面长度选择对计算结果的影响.结合Monte-Carlo方法,数值分析了理想导体Gauss 关键词： 复合散射 Kirchhoff近似 共轭梯度法 互耦迭代     

分形粗糙面双站散射的快速前后向迭代法数值模拟

为模拟复杂分形表面特别是在低掠角入射条件下的双站散射,发展了一种结合前后向迭代方法(FBM)与谱加速算法(SAA)快速求解散射场的Monte Carlo数值方法,计算了在TE,TM锥形波入射在一维分形导体粗糙面的双站散射以及有规则异物存在时的双站散射,讨论了分形粗糙面双站散射的角度性分布与其分数维的关系. 关键词： GFBM/SAA 分形粗糙面 双站散射     

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

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

FDTD Investigation on Electromagnetic Scattering from Two-Layered Rough Surfaces under UPML Absorbing Condition

Electromagnetic scattering from one-dimensional two-layered rough surfaces is investigated by using finite-difference time-domain algorithm （FDTD）. The uniaxial perfectly matched layer （UPML） medium is adopted for truncation of FDTD lattices, in which the finite-difference equations can be used for the total computation domain by properly choosing the uniaxial parameters. The rough surfaces are characterized with Gaussian statistics for the height and the autocorrelation function. The angular distribution of bistatic scattering coefficient from single-layered perfect electric conducting and dielectric rough surface is calculated and it is in good agreement with the numerical result with the conventional method of moments. The influence of the relative permittivity, the incident angle, and the correlative length of two-layered rough surfaces on the bistatic scattering coefficient with different polarizations are presented and discussed in detail.     

Application of a multiregion model to the EM scattering from a rough surface with or without a target above it

An efficient multiregion model is introduced to calculate the electromagnetic scattering from a perfectly electrical conducting(PEC) rough surface with or without a PEC target above it.In the multiregion model,the rough surface is split into multiple regions depending on their position along the rough surface.Two intermediate regions are chosen as the dominant region.If a target is located above the rough surface,the target will also be included in the dominant region.The method of moments(MOM) is only adopted on the dominant region to ensure validity.Hence,the new model can greatly reduce the number of unknowns associated with full MOM analysis.The induced electric currents on the other regions are obtained by approximately considering the mutual coupling between different regions along the rough surface.Compared with the published hybrid method,this new model is not only suitable for EM scattering from a target above a rough surface but also applicable for just rough surfaces.Several numerical simulations are presented to show the validity and efficiency of the multiregion model.     

Bistatic scattering and depolarization by randomly rough surfaces: application to the natural rough surfaces in X-band

Abstract

The scattering of waves by random rough surfaces has important applications in the remote sensing of oceans and land. The problem of developing a model for rough surfaces is very difficult since, at best, the scattering coefficient σ⁰ is dependent upon (at least) the radar frequency, geometrical and physical parameters, incident and observation angles, and polarization. The problem of electromagnetic scattering from a randomly rough surface is analysed using the Kirchhoff approximation (stationary phase, scalar approximation), the small-perturbation model and the two-scale models. A first major new consideration in this paper is the polarimetric signature calculations as a function of the transmitter location and receiver location for a bistatic radio-link. We calculate the like- and cross-polarized received power directly using the scattering coefficients, without calculating the Mueller matrix. Next, a study of the regions of validity of the Kirchhoff and small-perturbation rough surface scattering models (in the bistatic case) is presented. Comparisons between the numerical calculations and the models are made for various surface rms heights and correlation lengths both normalized to the incident wavenumber (denoted by σ and L, respectively). By using these two parameters to form a two-dimensional space, the approximate regions of validity are then established. The second major new consideration is the development of a theoretical two-scale model describing bistatic reflectivity as well as the numerical results computed for the bistatic radar cross section from rough surfaces especially from the sea and snow-covered surfaces. The results are used to show the Brewster angle effect on near-grazing angle scattering.     

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.     

Modulated Phase of a Potts Model with Competing Binary Interactions on a Cayley Tree

We study the phase diagram for Potts model on a Cayley tree with competing nearest-neighbor interactions J ₁, prolonged next-nearest-neighbor interactions J _p and one-level next-nearest-neighbor interactions J _o . Vannimenus proved that the phase diagram of Ising model with J _o =0 contains a modulated phase, as found for similar models on periodic lattices, but the multicritical Lifshitz point is at zero temperature. Later Mariz et al. generalized this result for Ising model with J _o ≠0 and recently Ganikhodjaev et al. proved similar result for the three-state Potts model with J _o =0. We consider Potts model with J _o ≠0 and show that for some values of J _o the multicritical Lifshitz point be at non-zero temperature. We also prove that as soon as the same-level interactionJ _o is nonzero, the paramagnetic phase found at high temperatures for J _o =0 disappears, while Ising model does not obtain such property. To perform this study, an iterative scheme similar to that appearing in real space renormalization group frameworks is established; it recovers, as particular case, previous work by Ganikhodjaev et al. for J _o =0. At vanishing temperature, the phase diagram is fully determined for all values and signs of J ₁,J _p and J _o . At finite temperatures several interesting features are exhibited for typical values of J _o /J ₁.     

Bistatic scattering from forest components. Part I: coherent polarimetric modelling and analysis of simulated results

This article deals with the study of bistatic scattering by forested areas using a coherent scattering model. The contribution of the scatterers is considered (trunks, branches), as well as their interaction with the ground. This code has been extended from monostatic to bistatic configurations and we present all the tests we have performed with this tool in order to validate it as far as possible. The transmitter is at a fixed location and the receiver moves either in the incidence plane (only θ _s varies) or in all positions (θ _s and φ _s both vary). Our objective is to determine the main trends characterizing bistatic scattering by a forest. We study the variation of the scattering coefficient, when evolving from the backscattering to the specular configuration (only θ _s varies). Both radiometric and polarimetric aspects are analysed, first when considering a simple forest of trunks and then a forest with trunks and branches. Finally, we analyse the impact of both θ _s and φ _s for two cases: a vertical dielectric cylinder in free space and a set of vertical dielectric cylinders, arranged in different ways. Here again, the radiometric and polarimetric aspects are investigated to check the electromagnetic behaviour.     

Bethe-ansatz solution of a model for a mixed valent impurity with two magnetic configurations: II. Numerical solution of the thermodynamic equations

A mixed valence impurity with two magnetic configurations of total angular momentumJ ₂ andJ ₁=J ₂+1/2, respectively, coupled by conduction electrons with total angular momentum 1/2 via a hybridization matrix element is considered. The thermodynamic Bethe-ansatz equations derived previously are solved numerically for various values ofJ ₂. Thef-level occupation, the entropy, the magnetic susceptibility and the specific heat are obtained as a function of temperature for variousf-level positions. The magnetic field dependence is also discussed in the limit of integer valence (exchange model).Supported by the CONICET, Argentina     

A method to separate conservative and magnetically-induced electric fields in calculations for MRI and MRS in electrically-small samples

This work presents a method to separately analyze the conservative electric fields (E_c, primarily originating with the scalar electric potential in the coil winding), and the magnetically-induced electric fields (E_i, caused by the time-varying magnetic field B1) within samples that are much smaller than one wavelength at the frequency of interest. The method consists of first using a numerical simulation method to calculate the total electric field (E_t) and conduction currents (J), then calculating E_i based on J, and finally calculating E_c by subtracting E_i from E_t. The method was applied to calculate electric fields for a small cylindrical sample in a solenoid at 600 MHz. When a non-conductive sample was modeled, calculated values of E_i and E_c were at least in rough agreement with very simple analytical approximations. When the sample was given dielectric and/or conductive properties, E_c was seen to decrease, but still remained much larger than E_i. When a recently-published approach to reduce heating by placing a passive conductor in the shape of a slotted cylinder between the coil and sample was modeled, reduced E_c and improved B1 homogeneity within the sample resulted, in agreement with the published results.     

微粗糙面上方球形粒子的光散射及其散射截面的计算

基于互易定理研究了光波入射时微粗糙面与其上方球形粒子复合模型的光散射。根据粗糙表面电流积分方程并利用表面微扰展开，得到了微粗糙面表面极化电流的迭代解，给出了耦合电场的计算方法。结合散射耦合场散射矩阵和已有的微粗糙面及球形粒子的散射矩阵，给出了复合模型散射截面的计算公式，数值计算了复合模型的后向散射截面并进行了详细讨论。     

雪层覆盖的粗糙地面与上方矩形截面柱复合电磁散射的时域有限差分法研究

采用土壤介电常数的四成分模型和雪的介电常数模型分别表示实际的地面和雪层的介电特性, 应用指数型分布粗糙面模型和Monte Carlo方法模拟实际的粗糙地面, 运用时域有限差分方法研究了雪层覆盖的粗糙地面与上方矩形截面柱的复合电磁散射问题. 得出了复合电磁散射系数的角分布曲线, 计算了双站复合散射系数随土壤与雪层粗糙度参数、介电参数、 矩形截面柱几何参数、介电参数等的变化情况, 并做了详细分析与讨论. 得到了雪层覆盖的粗糙地面与上方矩形截面柱复合电磁散射特性.     

Analysis of a penny-shaped crack in a magneto-electro-elastic medium

The solution for an ellipsoidal cavity in an infinite transversely isotropic magneto-electro-elastic medium under remotely applied axisymmetric, combined mechanical–electric–magnetic loading is derived. The exact solution for a penny-shaped crack is obtained by letting the minor axis of the ellipsoidal cavity approach zero. The results demonstrate that the stress intensity factor depends only on the applied mechanical loading, but the electric displacement intensity factor and the magnetic induction intensity factor depend on the applied combined loading, as well as the two ratios of β/α and γ/α. Parameter α is the ratio of the minor axis to the major axis of the ellipsoidal cavity, β is the ratio of the dielectric constant of the material in the cavity to the effective dielectric constant of the magneto-electro-elastic medium, and γ is the ratio of the magnetic permeability of the material in the cavity to the effective magnetic permeability of the magneto-electro-elastic medium. The two ratios characterize the permeability of the crack to electric and magnetic fields. Several limiting cases for β/α and γ/α are studied. A self-consistent method is adopted to determine the real crack opening α under the combined mechanical–electric–magnetic loading. The stress, electric displacement and magnetic induction intensity factors of a penny-shaped crack in BaTiO₃–CoFe₂O₄ composites are calculated for different volume fractions and different applied combined loadings.     

Surface optical waves at the boundary with a nonlinear Kerr medium

A rigorous solution consistent with a plane wave approximation is given to the boundary problem for Maxwell’s equations for surface optical waves at the boundary with a nonlinear Kerr medium. Exact formulas for the flux intensity (J ₀) and energy density (W ₀) of these waves are derived depending on the parameters of the adjacent media and the propagation constant (ξ). It is shown that these variables as functions of ξ have minima. Thus, J ₀ and W ₀ increase sharply as the propagation constant deviates from the minimum value ξ_min. Their values are greater, the larger the difference between the dielectric constants of the linear and nonlinear media is. An expression for the propagation velocity of a nonlinear surface wave is also obtained.     

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.