气溶胶凝聚粒子散射特性的数值计算

利用Cluster-cluster aggregation (CCA)模型,模拟了由相同数目球形原始微粒凝聚而成的四种随机取向气溶胶凝聚粒子．根据物质的电结构,将气溶胶凝聚粒子离散为一系列偶极子,结合离散偶极子近似方法,在获得每一个偶极子的电偶极矩之后,数值计算了气溶胶凝聚粒子散射强度的角分布,并分析了散射强度随入射光入射角度和气溶胶凝聚粒子尺寸参数变化的规律．结果显示：当散射角较小时,气溶胶凝聚粒子取向和入射光的入射角度对散射强度影响不大,当散射角增大时,散射强度则明显依赖于气溶胶凝聚粒子取向和入射光的入射角度;对于不同尺寸参数的气溶胶凝聚粒子,在同一角度入射情况下,随尺寸参数的增加,气溶胶凝聚粒子的散射主要集中于前向散射．     

Sensitivity of depolarized lidar signals to cloud and aerosol particle properties

Measurements from depolarized lidars provide a promising method to retrieve both cloud and aerosol properties and a versatile complement to passive satellite-based sensors. For lidar observations of clouds and aerosols, multiple scattering plays an important role in the scattering process. Monte Carlo simulations are carried out to investigate the sensitivity of lidar backscattering depolarization to cloud and aerosol properties. Lidar parameters are chosen to be similar to those of the upcoming space-based CALIPSO lidar. Cases are considered that consist of a single cloud or aerosol layer, as well as a case in which cirrus clouds overlay different types of aerosols. It is demonstrated that besides thermodynamic cloud phase, the depolarized lidar signal may provide additional information on ice or aerosol particle shapes. However, our results show little sensitivity to ice or aerosol particle sizes. Additionally, for the case of multiple but overlapping layers involving both clouds and aerosols, the depolarized lidar contains information that can help identify the particle properties of each layer.     

Multiple scattering of light by water cloud droplets with external and internal mixing of black carbon aerosols

The mixture of water cloud droplets with black carbon impurities is modeled by external and internal mixing models.The internal mixing model is modeled with a two-layered sphere(water cloud droplets containing black carbon(BC) inclusions),and the single scattering and absorption characteristics are calculated at the visible wavelength of 0.55 μm by using the Lorenz-Mie theory.The external mixing model is developed assuming that the same amount of BC particles are mixed with the water droplets externally.The multiple scattering characteristics are computed by using the Monte Carlo method.The results show that when the size of the BC aerosol is small,the reflection intensity of the internal mixing model is bigger than that of the external mixing model.However,if the size of the BC aerosol is big,the absorption of the internal mixing model will be larger than that of the external mixing model.     

沙尘气溶胶粒子群的散射和偏振特性

根据Mie散射理论,以对数正态分布函数描述沙尘气溶胶粒子群的粒径尺度分布,计算了沙尘气溶胶粒子群在0.2～40μm波段间对太阳短波辐射和地球大气长波辐射的单次散射反照率、散射相矩阵函数,揭示了不同相对湿度时,沙尘粒子群对入射辐射的散射和偏振的特征。结果表明,沙尘粒子群的单次散射反照率随着入射波长的增加有较大起伏,不同相对湿度条件下,变化趋势基本一致;在可见光、近红外波段单次散射反照率随湿度增加而变大,湿度95%时非常接近于1;大于10μm的热红外波段单次散射反照率随相对湿度增加而减小,具有较强的吸收辐射能力。散射辐射强度受湿度影响较小,随散射角的增加呈现先减小后增大的趋势,且增大的趋势随着波长的增加而减弱;不同波段上,线偏振和圆偏振随散射角和相对湿度变化存在差异;在前向和后向仅对入射辐射为圆偏振辐射产生圆偏振散射;散射光的偏振特性及其湿度差异主要表现在后向散射区,多以拱形形式体现。拱顶峰值散射角位置存在差异,且峰值散射角随相对湿度的降低向后向漂移。     

利用DDA方法计算大气气溶胶粒子光学特性

大气气溶胶是大气的重要组成部分,其光学特性是研究大气辐射传输特性的重要参量。本文基于DDA方法,对不同形状气溶胶粒子的光学特性进行计算,得到气溶胶粒子的消光因子、吸收因子随波长变化的数值结果。结合Muller散射矩阵,给出了气溶胶粒子散射强度和极化度的角分布,为研究大气辐射传输提供了有效的方法。     

近红外波段气溶胶的衰减特性研究

在已知大气气溶胶折射率和气溶胶谱分布的基础上, 对近红外波段的气溶胶衰减特性进行了研究。利用Mie散射理论计算并讨论了气溶胶的消光、散射、吸收效率因子随尺度参数的变化和消光系数随半径和波长的变化, 并且在MATLAB中对各种变化情况进行了仿真。结果表明, 三种气溶胶粒子的消光和散射能力依次为沙尘性粒子, 水溶性粒子, 烟煤。消光系数在粒子半径和入射波长相近时达到最大, 并且粒子半径对消光、散射、吸收系数的影响比入射波长更明显。这些结论可以为红外辐射在大气中的衰减计算和分析提供依据。     

Scattering of elastic waves by multiple elastic circular cylinders with imperfect interface

In the current paper a general method is presented for the rigorous solution for the scattering of elastic waves by a cluster of elastic circular cylinders of infinite length. The interface separating the cylinder from the surrounding media is considered to be homogeneous imperfect. Specifically, the tractions are continuous but the displacements are discontinuous and proportional in terms of interface stiffness parameters to their respective traction components. Using the exact theory of multipole expansion, analytic solutions for the scattered and internal fields excited by an incident plane P-wave, an incident cylindrical P-wave and an incident plane SV-wave are derived.

Numerical results for directivity patterns and scattering cross-sections are presented for a finite hexagonal array of elastic circular inclusions with imperfect interface. The results show that the sequence of maxima and minima in the curves of scattered cross-sections becomes more undistinguishable as the interface becomes more imperfect. Also, the results reveal that large low-frequency peaks of the scattered cross-sections, which correspond to resonance scattering, can be observed for both the low-velocity and high-velocity elastic cylinders with extremely imperfect interface while the small high-frequency peaks of the scattered cross-sections can appear for low-velocity elastic cylinders with relatively perfect interface. Furthermore, the results clearly show that the interaction effects between cylinders cannot be ignored for an incident plane SV-wave as compared to an incident plane P-wave. More importantly is the fact that the reciprocity relations, which hold for elastic wave scattering by a single cylinder, no longer apply for elastic wave scattering by multiple cylinders.     

Scattering of elastic waves by multiple elastic circular cylinders with imperfect interface

In the current paper a general method is presented for the rigorous solution for the scattering of elastic waves by a cluster of elastic circular cylinders of infinite length. The interface separating the cylinder from the surrounding media is considered to be homogeneous imperfect. Specifically, the tractions are continuous but the displacements are discontinuous and proportional in terms of interface stiffness parameters to their respective traction components. Using the exact theory of multipole expansion, analytic solutions for the scattered and internal fields excited by an incident plane P-wave, an incident cylindrical P-wave and an incident plane SV-wave are derived.

Numerical results for directivity patterns and scattering cross-sections are presented for a finite hexagonal array of elastic circular inclusions with imperfect interface. The results show that the sequence of maxima and minima in the curves of scattered cross-sections becomes more undistinguishable as the interface becomes more imperfect. Also, the results reveal that large low-frequency peaks of the scattered cross-sections, which correspond to resonance scattering, can be observed for both the low-velocity and high-velocity elastic cylinders with extremely imperfect interface while the small high-frequency peaks of the scattered cross-sections can appear for low-velocity elastic cylinders with relatively perfect interface. Furthermore, the results clearly show that the interaction effects between cylinders cannot be ignored for an incident plane SV-wave as compared to an incident plane P-wave. More importantly is the fact that the reciprocity relations, which hold for elastic wave scattering by a single cylinder, no longer apply for elastic wave scattering by multiple cylinders.     

A New Way of Looking at Nuclear Resonant Scattering

A new model for nuclear-resonant scattering of gamma radiation from resonant matter has been developed and is summarized here. This “coherent-path” model has lead to closed-form, finite-sum solutions for radiation scattered in the forward direction. The solution provides a unified microscopic picture of nuclear-resonant scattering processes. The resonant absorber or scatterer is modeled as a one-dimensional chain of “effective” nuclei or “effective” planes. The solution is interpreted as showing that the resonant radiation undergoes sequential scattering from one absorber “nucleus” or “plane” to another before reaching the detector. For recoil-free processes the various “paths” to the detector contribute coherently. The solution for this case gives calculated results that are indistinguishable from those using the classical optical model approach, although the forms of the solutions are completely different. The coherent-path model shows that the “speed-up” and “dynamical beating” effects are primarily a consequence of the fact that the single “effective” nuclear scattering processes are 180° out of phase with the incident radiation while the double nuclear scattering processes are in phase with the incident radiation. All multiple scattering paths are, and must be, included. The model can also treat the incoherent processes, i.e., processes involving gamma emission with recoil or conversion-electron emission. The source of the resonant gamma radiation can be from a radioactive source or from synchrotron radiation: both cases are treated. The model is used to explain and understand the results when each of the following experimental procedures is applied: time-differential Mössbauer spectroscopy, time-differential synchrotron radiation spectroscopy, enhanced-resolution resonant-detector Mössbauer spectroscopy, and the “gamma echo”.     

Refocusing of a high-power femtosecond laser pulse after a layer of atmospheric aerosol

The refocusing of a high-power femtosecond laser pulse after filamentation in a layer of atmospheric aerosol has been numerically investigated. It is shown that the presence of an on-path cloud of aerosol particles in the atmosphere results in the decrement of the filament energy, and the process of pulse refocusing becomes stochastic after the aerosol layer due to the competition between the influence of perturbations induced by aerosol scattering and filtration of radiation spatial modes.     

热红外对地遥感中的大气散射效应

从辐射传输方程出发,严格地分析了以往工作中被忽略的大气层中悬浮微粒散射对热红外辐射能量传输的影响,特别研究了对能量平衡的贡献和对传感器信噪比的影响,对以遥感器高度和能见度作为参量的仿真计算结果说明,与传统的不考虑大气中悬浮微粒的散射情况相比差1~2个数量级。这说明,悬浮微粒的散射不仅对红外辐射有消光作用,而且是研究热红外遥感系统中新的辐射通量出现的重要因素。     

Infrared radiative transfer modelling in a 3D scattering cloudy atmosphere: Application to limb sounding measurements of cirrus

The Monte Carlo cloud scattering forward model (McClouds_FM) has been developed to simulate limb radiative transfer in the presence of cirrus clouds, for the purposes of simulating cloud contaminated measurements made by an infrared limb sounding instrument, e.g. the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). A reverse method three-dimensional Monte Carlo transfer model is combined with a line-by-line model for radiative transfer through the non-cloudy atmosphere to explicitly account for the effects of multiple scattering by the clouds. The ice cloud microphysics are characterised by a size distribution of randomly oriented ice crystals, with the single scattering properties of the distribution determined by accurate calculations accounting for non-spherical habit.A comparison of McClouds_FM simulations and real MIPAS spectra of cirrus shows good agreement. Of particular interest are several noticeable spectral features (i.e. H₂O absorption lines) in the data that are replicated in the simulations: these can only be explained by upwelling tropospheric radiation scattered into the line-of-sight by the cloud ice particles.     

USING THE MULTIPLE SCATTERING THEORY FOR CALCULATION OF THE RADIATION FLUXES FROM EXPERIMENTAL AEROSOL DATA

The ground-based measurements of the spectral direct solar radiation fluxes and spectral sky radiances were basis for obtaining the aerosol characteristics such as columnar distribution function and effective refractive index, respectively. These characteristics were retrieved by a numerical solution of the inverse problem. The particle columnar size distribution function was chosen in the form of modified gamma distribution. The obtained aerosol parameters were utilized as input data for the multiple radiation scattering model (MRSM) to calculate the radiation fluxes. The dependence of calculation error of individual radiation characteristics (sky radiance, global and diffuse radiation) on a given number of terms NT (scattering orders) in the MRSM is studied. The NT necessary for successful calculation of radiation characteristics (for chosen accuracy) for actual atmospheric conditions is shown (aerosol optical properties were characterized, e.g., by spectral optical thickness). The approximate formulae for calculating the total sky radiances and total radiation fluxes, are find. Applying these formulae the value NT′ smaller than NT can be used for calculating radiation characteristics with the same accuracy. The scattering orders higher than NT′ are calculated on the basis of a known structure of previous scattering orders (1→NT′).     

空间垂直方向气溶胶消光特性的研究

在已知大气气溶胶折射率和气溶胶谱分布的基础上,对空间垂直方向的气溶胶消光特性进行了研究。文章利用Mie散射理论计算讨论了气溶胶的消光、散射、吸收系数随各参数的变化,并在MATLAB中对变化情况进行了仿真。结果表明,粒子半径与入射波长大小相近时气溶胶的消光最强,并且随着能见度的增大,气溶胶衰减系数减小。这些结论可以为红外辐射的大气衰减计算和分析提供依据。     

随机取向气溶胶凝聚粒子光学特性研究

利用Cluster-cluster aggregation (CCA)模型,模拟了由64个球形原始微粒凝聚而成的四种随机取向气溶胶凝聚粒子。结合离散偶极子近似方法,分别得到了不同入射角和不同尺寸参数情况下气溶胶凝聚粒子的吸收、散射、消光效率因子和不对称因子的数值结果,并比较分析了四种形状气溶胶凝聚粒子吸收、散射、消光效率因子和不对称因子的差异。结果显示：对于相同数目原始微粒的气溶胶凝聚粒子,其光学特性明显依赖于入射光的入射角度和气溶胶凝聚粒子的形状;对于不同尺寸参数的气溶胶凝聚粒子,当入射光波长给定时,吸收、散射和消光效率因子最初随尺寸参数增大而快速增大,随后又随尺寸参数增大缓慢变小,但在尺寸参数变化过程中存在一个极大值;不对称因子则随尺寸参数的增大一直呈增大趋势,且随尺寸参数的增大趋近于1。     

Numerical calculation of apparent IR radiation of cloud

The infrared detection has been one of important approaches for aerial target detection, but the existence of clouds in the sky makes target detection difficult, so it is of great significance to research the features of infrared radiation of clouds. Combining Mie scattering theory, the calculation model of infrared radiation on the cloud appearance was put forward and this model mainly considered two parts including direct and scattering radiation. The calculation methods of direct radiation, primary scattering and multiple scattering were discussed in detail to put forward the specific methods for numerical computation. Based on the calculation model, cumulus was taken as the example to calculate its radiation of wave band of long-wave infrared and analyze its radiation characters. The calculation results showed that: both perpendicular incidence and scattering could make great contributions to apparent radiation of cloud and could not be ignored; as for the distribution of direction angles, the radiation of clouds decreased with increase of direction angles as a whole; in the aspect of spectrum distribution, spectral radiance feature of the clouds was complicated function of wave length. Finally, the calculation results were verified using clouds imaging experiment of which the result also effectively proved the correctness of the calculation model.     

Incoherent properties of induced radiation

The evolution of a quantized electromagnetic field in a thermally excited dispersion medium is determined by two scattering channels. The coherent channel is formed exclusively by the elastic scattering of quanta. The incoherent channel, along with elastic scattering processes, necessarily contains inelastic scattering processes, including induced radiation. Interference between the channels is absent because of the orthogonality of the wave functions of the medium in its final states, which correspond to different scattering channels. Therefore, in an excited medium, interference processes that are not described by its refractive index may arise. An interference pattern of this kind can be formed, in particular, as a result of the superposition of the resonance radiation incident on an excited medium and the radiation reflected from this medium. In this case, the conventional perturbation theory proves to be inadequate.     

Time-domain numerical simulations of multiple scattering to extract elastic effective wavenumbers

Elastic wave propagation is studied in a heterogeneous two-dimensional medium consisting of an elastic matrix containing randomly distributed circular elastic inclusions. The aim of this study is to determine the effective wavenumbers when the incident wavelength is similar to the radius of the inclusions. A purely numerical methodology is presented, with which the limitations usually associated with low scatterer concentrations can be avoided. The elastodynamic equations are integrated by a fourth-order time-domain numerical scheme. An immersed interface method is used to accurately discretize the interfaces on a Cartesian grid. The effective field is extracted from the simulated data, and signal-processing tools are used to obtain the complex effective wavenumbers. The numerical reference solution thus obtained can be used to check the validity of multiple scattering analytical models. The method is applied to the case of concrete. A parametric study is performed on longitudinal and transverse incident plane waves at various scatterer concentrations. The phase velocities and attenuations determined numerically are compared with predictions obtained with multiple scattering models, such as the Independent Scattering Approximation model, the Waterman–Truell model, and the more recent Conoir–Norris model.