非均匀随机散射介质层的多次散射和热辐射的分层迭代解

为纳入矢量辐射传输方程的多次散射，将非均匀平行分层分布的随机散射介质划分成多个薄 层，利用各薄层的低阶Mueller矩阵解推导出整层的高阶散射迭代解.用这一方法计算了具有 占空比、粒子大小和温度廓线非均匀分布的平行分层随机密集球形粒子积雪层的极化热辐射 .数值结果与一层均匀随机散射层的离散坐标-特征值方法作了比较，讨论了多次散射和热辐 射与散射层各特征参数的关系，与积雪层微波遥感数据进行了对比. 关键词： VRT方程 平行分层 高阶散射解 迭代解     

一层随机非球形粒子高阶散射矢量辐射传输迭代解的参数求逆

根据一层随机非球形粒子矢量辐射传输(VRT)方程的一阶与高阶Mueller矩阵解的差异,提出了在已知一些参数的条件下,用水平与垂直同极化和交叉极化后向散射测量反演非球形粒子复介电常数和单位面积粒子数的迭代方法,并用极化后向散射强度的模拟测量数据及机载合成孔径雷达(Air SAR)实测数据进行了反演试验,并对反演值与各参数的关系以及测量噪声对反演的影响进行了讨论.     

随机非球形粒子全极化散射的时间相关Mueller矩阵解

从与时间相关的矢量辐射传输方程推导一阶Mueller矩阵解,用来模拟Gauss型平面脉冲波入射下,一层随机、非均匀取向非球形粒子的全极化双站散射.数值计算了同极化和去极化脉冲响应,与入射脉冲进行了比较,说明了随机介质的物理参数,如粒子的取向和占空比、入射角、极化以及层厚等对脉冲响应的影响 关键词： 平面脉冲波 非球形粒子 Mueller矩阵     

任意阶球谐函数法求解辐射传输问题的误差分析

针对各向异性散射介质内的红外辐射传输开展数值方法研究.应用球谐函数展开推导了一维吸收、发射、散射性灰介质的辐射传输近似方程组,采用差分方法结合三对角矩阵解法建立了任意阶辐射传输近似方程数值解法.并通过计算算例验证了高阶PN方法的数值无关性问题.计算结果表明,高阶球谐函数数值解与理论解吻合得很好,同时计算验证了不同展开阶...     

粗糙面上粒子层矢量辐射传输方程的二阶迭代解法

为了探究典型粗糙面上随机粒子层中能量传输的多次散射机制,提出了一种基于矢量辐射传输方程的建模二阶计算方法.该建模方法将建模场景(粗糙面上粒子层)在高度维(Z轴)划分为多个传输散射层,基于矢量辐射传输理论中的一阶迭代散射解,利用典型粗糙面的半经验半解析方法,求解出整个场景的二阶迭代散射解.同时,研究粒子层能量在粒子与粒子间的多次散射机制,以及粒子与地表粗糙面间的多次散射机制.数值结果表明,该二阶迭代解法相较于矢量辐射传输方程的一阶迭代散射解,能够更完整地探究互作用的散射特性,且可从能量传输角度解译建模场景中物体间的相干作用,从而可用于植被地物环境下的多次散射机制的解析以及散射系数变化趋势的预估.     

非球形粒子散射光的去偏振特性研究

目标对入射偏振光的散射特性反映了目标的属性信息。在Rayleigh散射理论的基础上，通过单层非球形粒子对入射偏振光的散射数学模型，应用矢量传输方程来计算非球形粒子散射的Mueller矩阵元，求解散射介质的偏振度，讨论在不同介质层厚度、粒子半径和探测角的条件下，非球形粒子散射光去偏振度的变化特性，给出了模拟仿真结果。该方法为研究目标的内部结构、厚度和粗糙度等特征以及目标的探测和识别提供了一种新的途径。     

随机介质中的矢量辐射传输理论

本文讨论了极化电磁波在随机介质中多次散射，传输和热辐射的斯托克斯矢量的辐射传输理论。其中包括随机分布离散的球形和非球形粒子的矢量辐射传输方程，离散坐标——特征分析法，付利叶变换，迭代法等数值解。讨论了非球形粒子的穆勒矩阵。并研究了密集分布的散射粒子介质的辐射传输理论，考虑了密集粒子散射的相干性，计算了有效传播常数。理论及数值结果与实验用了很好的比较。     

一层非球形粒子散射的标量辐射传输迭代解的求逆

给出了一层有下垫反射面的非球形粒子散射的标量辐射传输方程一、二阶迭代解,推导了相函数各阶Legendre展开系数与随机小椭球粒子相函数的对应关系,提出了一层随机小椭球粒子介电常数和单位面积粒子数的迭代反演方法.通过两次各方位角上双站散射测量,反演随机小椭球粒子的介电常数和单位面积粒子数.本方法比现有文献的一些反演方法易于实现,可应用于颗粒性复合材料介电与结构特性以及地表背景参数的反演研究 关键词： 辐射传输方程 迭代解 反演 相函数 介电常数     

基于新型混合场积分方程的半空间三维均匀介质目标电磁散射分析

利用介质半空间格林函数,将一种混合场积分方程(JMCFIE)推广应用到有耗半空间的均匀介质的散射分析中,以减少离散矩阵方程迭代求解时所需的迭代步数,提高矩量法(MoM)对半空间介质目标散射分析的求解效率.数值算例验证了方法的准确性和有效性.     

随机介质中的矢量辐射传输理论

本文讨论了极化电磁波在随机介质中多次散射,传输和热辐射的斯托克斯矢量的辐射传输理论。其中包括随机分布离散的球形和非球形粒子的矢量辐射传输方程,离散坐标-特征分析法,付利叶变换,迭代法等数值解。讨论了非球形粒子的穆勒矩阵。并研究了密集分布的散射粒子介质的辐射传输理论,考虑了密集粒子散射的相干性,计算了有效传播常数。理论及数值结果与实验作了很好的比较。     

Iterative approach of high-order scattering solution for vector radiative transfer of inhomogeneous random media

By stratifying a random scatter media into multiple thin layers in the vertical z direction, the first-order scattering solution of each thin layer is employed to derive high-order scattering solution of whole random media. Using the Fourier transform and Mueller matrices in discrete ordinates, an iterative approach to solve high-order scattering solution of vector radiative transfer (VRT) equation is newly developed. Numerical results are well compared with the Mueller matrix solutions of the first order for a single layer medium, second order for a half-space, and the results of the discrete ordinate and eigen analysis method. It demonstrates our approach as feasible, effective and especially applicable to high-order solution of VRT for both bistatic scattering and thermal emission of inhomogeneous non-spherical scatter media.     

Exact solution of Maxwell's equations for optical interactions with a macroscopic random medium

We report what we believe to be the first rigorous numerical solution of the two-dimensional Maxwell equations for optical propagation within, and scattering by, a random medium of macroscopic dimensions. Our solution is based on the pseudospectral time-domain technique, which provides essentially exact results for electromagnetic field spatial modes sampled at the Nyquist rate or better. The results point toward the emerging feasibility of direct, exact Maxwell equations modeling of light propagation through many millimeters of biological tissues. More generally, our results have a wider implication: Namely, the study of electromagnetic wave propagation within random media is moving toward exact rather than approximate solutions of Maxwell's equations.     

Coherent effects of multiple scattering for scalar and electromagnetic fields: Monte-Carlo simulation and Milne-like solutions

Based on an expansion of the Bethe-Salpeter equation in scattering orders a semi-analytic approach for simulation of coherent phenomena of multiple scattering in random media has been developed. We found that for scalar field the manifestation of these phenomena, observed as temporal field correlation function and coherent backscattering, are universal and well agreed with the results predicted by diffusion approximation. For the electromagnetic field the temporal correlation function and coherent backscattering are noticeably differ from those found for the scalar field, depending strongly on the scattering anisotropy. The obtained numerical results, for the first time to our knowledge, are compared directly with the known generalizations of the Milne solution.     

Reflection of waves from a strong scatterer buried in a random medium

Reflection of waves from a mirror covered by a random layer of isotropic, absorbing scatterers is studied and the angular distribution of the scattered intensity is calculated both analytically and numerically. It is shown that backscattering enhancement as well as an enhancement of the incoherent signal in the specular direction take place even in the singly scattered random field. The dependence of the retroreflected intensity is shown to be a non-monotonic function of the depth of the mirror, with a maximum at a depth of the order of the scattering mean free path. Possibilities for employing the results obtained to detect buried strong scatterers and to retrieve parameters of the random media are discussed. In particular, it is shown that in the case of strong absorption the reflecting plane manifests itself by the presence of a peak in the retroreflected intensity which is missing from the scattering diagram of a free-standing or an infinitely thick random layer.     

A successive scattering methodology with application to two cylinders in the Fresnel region of each other

We present an efficient approach to compute the second-order scattering of an electromagnetic wave by two discrete scatterers in proximity to each other. Such a two-body system represents the simplest canonical arrangement to address near-field volume scattering phenomena in microwave remote sensing models of vegetation. Using an analytical wave-based approach, a successive scattering methodology is employed to derive the first interaction term in multiple scattering by two arbitrary scatterers in terms of their transition operators. The general formulation is applied to find the second-order bistatic scattering amplitude for a pair of finite length thin cylinders at arbitrary interaction distances using the exact Green's function. To improve computational efficiency, the solution is then specialized to the Fresnel region. These second-order bistatic scattering amplitude results are in agreement with the exact Green's function model when the scatterers are in the Fresnel region of each other. Additionally, it is demonstrated that using the far field approximation in the Fresnel region can yield significant deviations from the exact results. The Fresnel model, unlike the far field approximation, accurately predicts the scattering amplitude peak values and null locations, and is suited to fast solutions in realistic canopy simulations.     

Modeling of high-frequency propagation in inhomogeneous background random media

When a high-frequency electromagnetic wave propagates in a complicated scattering environment, the contribution at the observer is usually composed of a number of field species arriving along different ray trajectories. In order to describe each contribution separately the parabolic extension along an isolated ray trajectory in an inhomogeneous background medium was performed. This leads to the parabolic wave equation along a deterministic ray trajectory in a randomly perturbed medium with the possibility of presenting the solution of the high-frequency field and the higher-order coherence functions in the functional path-integral form. It is shown that uncertainty considerations play an important role in relating the path-integral solutions to the approximate asymptotic solutions. The solutions for the high-frequency propagators derived in this work preserve the random information accumulated along the propagation path and therefore can be applied to the analysis of double-passage effects where the correlation between the forward-backward propagating fields has to be accounted for. This results in double-passage algorithms, which have been applied to analyze the resolution of two point scatterers. Under strong scattering conditions, the backscattering effects cannot be neglected and the ray trajectories cannot be treated separately. The final part is devoted to the generalized parabolic extension method applied to the scalar Helmholtz's equation, and possible approximations for obtaining numerically manageable solutions in the presence of random media.     

Diffraction tomography: a geometrical distortion free procedure

Ductal echography (DE) is a recent anatomically-led method of investigation of internal mammary structures that provides direct observation of ductolobular structures in mammary lobes. Indeed, breast cancer initiates from the epithelium (hypoechogenic) and develops first in ductolobular structures. Our research aims at developing an ultrasonic scanner (hemispherical antenna) that will allow three-dimensional measurement of the field diffracted by the gland. The present reconstruction procedure is an alternative solution to the standard elliptic back-projection (EBP) technique based on the Born approximation (weak scattering assumption). The latter implicitly assumes that the sound speed is (almost) constant within the breast. This unrealistic a priori is inconsistent with our heterogeneous media characterization purpose, and leads both to geometrical distortion in the reconstruction and to poor focusing of the backpropagated waves (low contrast imaging). However, the EBP technique retains a two-fold advantage: firstly it does not make any other assumption concerning the distribution of mechanical parameters apart from low gradients. Secondly, it shows great robustness, a high resolving power, and is easy to implement. Thus, in order to account for strong wavefront distortions, we integrate temporal compensation of the scattered signals acquired into the scattering EBP technique. The adjustment relies on a time of flight estimation based on a "layer stripping" approach. Numerical tests based on finite difference time domain simulations of data scattered by a random tissue-like phantom are proposed.     

气动外形表面雷达吸波涂层问题电磁散射场的数值模拟

研究了将无覆盖区域分裂与精确控制法相结合、用于气动外形表面雷达吸波涂层问题电磁散射场的数值求解方法.无覆盖区域分裂是按空气和涂层两种不同介质物理分区要求进行的,子域交界处离散网格彼此之间可以是非匹配的,以便处理不同介质的不同网格尺度要求.电磁波通过子域交界面的约束条件是通过引入Lagrange乘子而弱满足.整个求解过程包含基于精确控制的外迭代和基于Lagrange乘子控制的区域分裂内迭代.最后计算给出了二维翼型气动外形涂层问题的电磁散射场,并与对应的无涂层外形的计算结果进行了比较,反映出表面涂层对电磁散射场的影响.