Solving the atmospheric scattering optical transfer function using the multi-coupled single scattering method

The atmospheric scattering optical transfer function （OTF） is solved by applying the multi-coupled single scattering （MCSS） method to the three-dimensional radiative transfer equation （RTE） under the periodic ground condition. This approach is a direct hit to the atmospheric scattering OTF using the same original context of modulation transfer function （MTF） measurement, i.e., images of sinusoidal grating at different spatial frequencies. Both the amplitude and phase shift of the OTF at various zenith and azimuth angles can be obtained at an arbitrary spatial frequency.     

气粒混合物非灰辐射特性合并宽窄谱带K分布模型

气粒混合物辐射问题具有全场性、非灰性、耦合性等特点,准确预估高温燃气/粒子非灰辐射特性是非常重要的。本文将合并宽窄谱带K分布模础(CWNBCK)与离散坐标法(DOM)结合,开展了非灰气粒混合物辐射换热问题的模拟工作,分别验证了一维和三维情况下应用该模型的准确性,给出不同工况下的热流源项、壁面热流或辐射热流等。结果表明:该模型能够给出与SNB模型精度基本相同的结果,考虑其计算效率的提高,可以在工程实际中应用该模型计算非灰气粒混合物辐射换热。     

Two-dimensional radiative transfer due to curved Dirac delta line sources

In this article, the two-dimensional radiative transport equation is considered for the curved Dirac delta line source. In order to account this source type, Green’s function of the radiative transport equation for the half-plane is derived in parts of the ballistic and diffuse contribution as well as under consideration of the Fresnel reflection at the boundary. The final results are verified with the Monte Carlo method for different collimated beams and several curved sources such as the elliptic and logarithmic spiral line source.     

Elementary processes involving Rydberg atoms and molecules in an intense laser radiation field

An analytical review of the theory of elementary atomic and molecular processes in the field of intense laser monochromatic radiation is presented. The discussion focuses on near-threshold processes involving Rydberg intermediate complexes, which play an important role in Earth’s upper atmosphere. The possibilities of the stationary method of radiative scattering matrix, based on the formalism of the renormalized Lippmann-Schwinger equations, are demonstrated. The approach is used to describe in a unified way a variety of near-threshold processes, the probability amplitudes of which are elements of the radiative scattering matrix for all possible reaction channels. They also include processes of restructuring of the particles (predissociation, associative ionization, exchange reactions, etc.), which in principle cannot be described using traditional nonstationary theories.     

Post-Hartree-Fock methods and dynamic correlation in atoms and molecules

The methods of calculation of the matrix of the exchange-correlation interaction are considered within the framework of one post-Hartree-Fock one-electron method of investigation of the properties of many-electron systems. Such post-Hartree-Fock methods are based on two-step variational self-consistent calculations of the spin orbitals and superposition coefficients of configurations in the multiconfiguration approximation. The post-Hartree-Fock method used involves an approach related to the extended Koopmans’ theorem, which, in turn, proves to be a high-energy approximation for quantum Green’s functions. Obvious application areas of the calculations of the exchange-correlation interaction within the framework of the method proposed are the multiparticle perturbation theory, the parameterization of the energy representation as a functional of the single-particle density matrix, and the theory of Green’s functions in the multiconfiguration approximation. A relativistic generalization of the method with the aim of calculating the radiative corrections for many-electron atoms and for problems of interaction with an external field in the nonstationary Floquet theory is possible.     

Comparison between direct-sequence and multicarrier spread-spectrum acoustic communications in time-varying channels

Underwater communication experiments have been conducted in the Norwegian Oslofjord. Two modulation schemes are compared in a 7-kHz frequency band on a 14-kHz center frequency. The first scheme is direct-sequence spread spectrum (DSSS), using a 7-chip spreading code to achieve a raw data rate of 1000 bps on a single carrier. The second scheme is multicarrier spread spectrum (MCSS) and accomplishes spreading by using seven subbands. The DSSS receiver equalizes on the chips prior to explicit symbol despreading, whereas MCSS features joint multiband equalization and despreading. Four channels are examined, from nearly static to overspread. In slowly varying channels, MCSS offers the best performance. DSSS has the best tracking potential for rapidly varying channels, where the challenge is to obtain reliable chip decisions before symbol despreading. The tracking potential can be realized to some extent by hypothesis-feedback equalization. It is further shown that adaptive equalizers are capable of code conversion, i.e., the DSSS receiver can demodulate the MCSS waveform, and vice versa. Neither receiver requires knowledge of the spreading code in order to despread the data.     

煤粉燃烧颗粒粒子云辐射特性的研究

分析了在560nm的测量波长处煤粉炉中炭粒和飞灰粒子云的辐射特性。利用Mie理论,分别计算了燃烧器出口和炉膛的上部颗粒粒子云的辐射特性,发现在燃烧器出口区,炭粒对颗粒总体辐射减弱系数的影响较大,占90%以上;而在炉膛的上部,飞灰的影响较大。从数值计算结果中可以看出,飞灰粒子云的散射率高达91.6%,其对辐射的散射作用比吸收更强,而炭粒子的散射率仅为50%。     

模拟辐射传热的离散坐标法的改进

针对离散坐标法模拟求解炉内辐射换热的问题,在对辐射传递方程采用控制容积法离散时,按辐射方向上各界面之间的投影关系划分区域,对辐射传递方程在各子控制体上沿辐射传输方向进行积分,求解各子控制体的辐射强度,再取容积平均值作为最终节点的辐射强度的新离散坐标方法,定义为分段积分离散坐标法,并对三维矩型燃气炉内辐射换热进行了数值模拟,与传统离散坐标方法及区域法进行了比较,比较结果表明新离散坐标法较好的解决了假散射问题,提高了离散坐标法的精度.     

纳米光学辐射传热:从热辐射增强理论到辐射制冷应用

热辐射作为一种无处不在的物理现象,对于科学研究和工程应用都具有重要意义.传统上对热辐射的理解主要是基于普朗克定律,它描述了物体通过辐射交换能量的能力.而近年来的研究表明,由于微纳光学材料在尺寸上远小于热辐射峰值波长,它们的热辐射性质往往很大程度上有别于传统黑体辐射理论所描述的宏观物体.更重要的是,微纳光学材料的热辐射性质可以通过改变它们的几何尺寸和微观构型进行定量的优化设计与精确调控.纳米光学材料与辐射制冷效应的结合,给热辐射效应在能源和环境等相关领域的应用提供了极具前景的应用价值.本文首先从热辐射的基本原理和规律出发,介绍纳米结构热辐射增强的发展进程和最新进展,包括二维材料间的近场热辐射机理以及尺寸效应导致的远场热辐射增强;其次,介绍了近年来纳米光学材料在辐射制冷应用中的重大进展,包括可以实现高效日间辐射制冷的各种纳米光学材料设计;最后,进一步介绍了日间辐射制冷的各种实际应用,包括建筑物制冷、冷凝水收集、舒适衣物与太阳能电池降温等.此外,展望了纳米光学材料的辐射制冷技术在推动荒漠生态环境的治理与改造方面的广阔未来.     

Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere-ocean system

A vector radiative transfer model has been developed for a coupled atmosphere-ocean system. The radiative transfer scheme is based on the discrete ordinate and matrix operator methods. The reflection/transmission matrices and source vectors are obtained for each atmospheric or oceanic layer through the discrete ordinate solution. The vertically inhomogeneous system is constructed using the matrix operator method, which combines the radiative interaction between the layers. This radiative transfer scheme is flexible for a vertically inhomogeneous system including the oceanic layers as well as the ocean surface. Compared with the benchmark results, the computational error attributable to the radiative transfer scheme has been less than 0.1% in the case of eight discrete ordinate directions. Furthermore, increasing the number of discrete ordinate directions has produced computations with higher accuracy. Based on our radiative transfer scheme, simulations of sun glint radiation have been presented for wavelengths of 670 nm and 1.6 μm. Results of simulations have shown reasonable characteristics of the sun glint radiation such as the strongly peaked, but slightly smoothed radiation by the rough ocean surface and depolarization through multiple scattering by the aerosol-loaded atmosphere. The radiative transfer scheme of this paper has been implemented to the numerical model named Pstar as one of the OpenCLASTR/STAR radiative transfer code systems, which are widely applied to many radiative transfer problems, including the polarization effect.     

Verification of numerical simulation method for entropy generation of radiation heat transfer in semitransparent medium

The numerical simulation method of radiative entropy generation in participating media presented by Caldas and Semiao [Entropy generation through radiative transfer in participating media: analysis and numerical computation. JQSRT 2005;96:423-37] is extended to analyze the radiative entropy generation in the enclosures filled with semitransparent media. A discrete ordinates method is used to solve radiative transfer equation and radiative entropy generation. Two different examples are employed to verify the numerical simulation method of radiative entropy generation in the enclosure. Numerical results of dimensionless radiative entropy generation of enclosure are identical to that of entire thermodynamics analysis for the enclosure system. This numerical simulation method can be used in the entropy generation analysis of high-temperature systems such as boilers and furnaces, in which radiation is the dominant mode of heat transfer.     

Linear and nonlinear optical properties of excitons in semiconductor quantum wells and microcavities

Linear and nonlinear light propagation in single and multiple quantum wells and in semiconductor microresonators are studied on the basis of Maxwell’s equations. The treatment includes radiative broadening of quantum-confined excitons, radiative coupling between quantum wells as well as coupling of quantum wells to the cavity field of a microresonator for steady state or ultrashort pulse excitation. The dynamical evolution of the coherent quantum-well polarization under the influence of many-body effects is studied within a microscopic model. The theory is used to investigate the influence of exciton saturation and dephasing on pulse propagation and excitonic normal-mode coupling.     

一维瞬态辐射传输的有限元法模拟

随着超短脉冲激光的快速发展,吸收散射性介质内的瞬态辐射传输引起了人们的广泛关注.本文基于离散坐标法和最小二乘有限元法(LSFEM),提出了模拟多维吸收散射性介质内瞬态辐射传输的数值模型.该模型有效地克服了在标准Galerkin有限元法(GFEM)中发生的伪振荡现象,在时间步长较大的情况下仍然可以得到光滑无振荡的解.而且,最小二乘法产生的求解系数矩阵是对称正定的,与GFEM中的系数矩阵相比,仅需要存储一半的非零系数,可以应用许多高效的迭代求解方法进行求解.为了检验模型,本文研究了一维吸收散射性介质内瞬态辐射传输问题,其结果与蒙特卡洛法(MCM)和积分模型法(IE)的结果进行了比较,结果证实:本文的方法可以精确、高效地模拟参与性介质内的瞬态辐射传输.     

An improved separability approximation for line radiative transport in nonhomogeneous media

A simple modification to the constant half-width approximation of Wilson and Greif, which is an extension of the well-known Curtis-Godson method to the treatment of temperature variations along the integration path, is introduced which permits a more accurate evaluation of line radiative transport in non- homogeneous gases. To demonstrate the method's accuracy, comparisons are made with Wilson and Greif and numerical frequency integrated results for the line equivalent width and radiative flux in a planar slab with prescribed Planck function and line half-width spatial variations are chosen to represent typical shock layer conditions. It is found that the modified procedure reduces the inaccuracies inherent in Wilson and Greif's approximation by factors ranging from 5 to 10, while retaining the latter method's ease of application.     

Probabilities of radiative transitions between Stark states in orthohelium

The wavefunctions, matrix elements, and probabilities of radiative transitions between Stark sub-levels of atomic multiplets are calculated as a function of the strength of a dc electric field. The general expressions for the wavefunction of a multiplet state in a field obtained by perturbation theory for close-lying levels with the use of completely reduced Green’s function allow one to determine the field dependence of both the dipole-allowed and dipole-forbidden radiation amplitudes. A decomposition of the second-order amplitude for the transition between fine-structure sublevels of two levels of equal parity into irreducible components is obtained. Numerical calculations of the probabilities of radiative transitions between triplet states of helium show the possibility of experimental observation of the emergence and vanishing of Stark lines of radiative transitions in the vicinity of anticrossing fields.     

Transfert de rayonnement : méthodes itératives

We review several numerical approaches which have been developed during the last 10 years in order to perform more realistic radiative modelling of solar and stellar atmospheres. However, we shall restrict ourselves to the same ‘family’ of method, the ones which are derived from the so-called Accelerated Lambda-Iteration technique.     

基于能量密度分布的辐射源粒子空间抽样方法研究

利用隐式蒙特卡罗方法模拟热辐射光子在物质中的输运过程时,物质辐射源粒子是需要细致处理的物理量.传统的物质辐射源粒子抽样方法是体平均抽样方法,对于大多数问题,这样处理不会带来大的偏差.但是对于一些辐射吸收截面大、单一网格内温差显著的问题,体平均抽样方法的计算结果偏差较大.分析了产生偏差原因,提出一种基于辐射能量密度分布的辐射源粒子空间位置抽样方法,并推导了相应的抽样公式以解决此类问题.数值实验表明,新方法计算结果明显优于原方法且与解析结果基本一致.     

Scaling algorithms for the calculation of solar radiative fluxes

We derived new scaling formulae based on the method of successive orders of scattering to calculate solar radiative flux. In this report, we demonstrate a multiple scaling method, in which we introduce scaling factors for each scattering order independently. The formula of radiative transfer by the method of successive orders of scattering cannot be solved rapidly except in the case of optically thin atmospheres. Then we further derived a double scaling method, which scales the ordinary radiative transfer equation by two scaling factors. We applied the double scaling method to two-stream and four-stream approximations of the discrete ordinates method. Comparing the results of the double scaling method with those of the delta-M method, we found that the double scaling method improved the accuracy of radiative fluxes at large solar zenith angles, especially in the optically thin region, and that in the region where multiple scattering dominates, its accuracy was comparable to that of the delta-M method. Once we determined the scaling factors appropriately, the double scaling method calculated radiative fluxes as rapidly as the delta-M method in the two-stream and four-stream approximations. This method, therefore, is useful for accurate computation of solar radiative fluxes in general circulation models.     

Development of a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI)

Satellite-observed radiance is a nonlinear functional of surface properties and atmospheric temperature and absorbing gas profiles as described by the radiative transfer equation (RTE). In the era of hyperspectral sounders with thousands of high-resolution channels, the computation of the radiative transfer model becomes more time-consuming. The radiative transfer model performance in operational numerical weather prediction systems still limits the number of channels we can use in hyperspectral sounders to only a few hundreds. To take the full advantage of such high-resolution infrared observations, a computationally efficient radiative transfer model is needed to facilitate satellite data assimilation. In recent years the programmable commodity graphics processing unit (GPU) has evolved into a highly parallel, multi-threaded, many-core processor with tremendous computational speed and very high memory bandwidth. The radiative transfer model is very suitable for the GPU implementation to take advantage of the hardware’s efficiency and parallelism where radiances of many channels can be calculated in parallel in GPUs.In this paper, we develop a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI) launched in 2006 onboard the first European meteorological polar-orbiting satellites, METOP-A. Each IASI spectrum has 8461 spectral channels. The IASI radiative transfer model consists of three modules. The first module for computing the regression predictors takes less than 0.004% of CPU time, while the second module for transmittance computation and the third module for radiance computation take approximately 92.5% and 7.5%, respectively. Our GPU-based IASI radiative transfer model is developed to run on a low-cost personal supercomputer with four GPUs with total 960 compute cores, delivering near 4 TFlops theoretical peak performance. By massively parallelizing the second and third modules, we reached 364× speedup for 1 GPU and 1455× speedup for all 4 GPUs, both with respect to the original CPU-based single-threaded Fortran code with the –O₂ compiling optimization. The significant 1455× speedup using a computer with four GPUs means that the proposed GPU-based high-performance forward model is able to compute one day’s amount of 1,296,000 IASI spectra within nearly 10 min, whereas the original single CPU-based version will impractically take more than 10 days. This model runs over 80% of the theoretical memory bandwidth with asynchronous data transfer. A novel CPU–GPU pipeline implementation of the IASI radiative transfer model is proposed. The GPU-based high-performance IASI radiative transfer model is suitable for the assimilation of the IASI radiance observations into the operational numerical weather forecast model.     

Hara’s theorem and W—exchange in hyperon weak radiative decays

We reconsider Hara’s theorem and relate it to the well-known properties of β-decay. All assumptions necessary for the theorem to be true are explicitly formulated. Further, we study the W—exchange contribution to weak radiative decays and show that it does not violate Hara’s theorem. However, this contribution reveals the essential role of particle mixing in symmetry considerations and some peculiar features of gauge-invariant amplitudes under perturbative expansion. Together they explain an effect, which was treated as contradicting Hara’s theorem, without any violation. The properties of W—exchange we describe here may have more general importance and should be taken into account in further detailed calculations of weak processes.