Radiative transfer with internal reflection via the converged discrete ordinates method

In response to the challenge of establishing highly accurate solutions to the plane layer radiative transfer equation with the simplest of methods, the converged discrete ordinates method is presented. With this algorithm of only finite difference, quadrature and acceleration, we show how to obtain highly accurate intensities for radiative transfer in a finite layer with internal surface reflection. The method features angular smoothing and angular interpolation through “faux” quadrature. In addition, a manufactured solution demonstrates the high accuracy of the method for forward peaked scattering. We consider scattering in a heterogeneous medium as a final demonstration.     

Benchmark results in vector atmospheric radiative transfer

In this paper seven vector radiative transfer codes are inter-compared for the case of underlying black surface. They include three techniques based on the discrete ordinate method (DOM), two Monte-Carlo methods, the successive orders scattering method, and a modified doubling-adding technique. It was found that all codes give very similar results. Therefore, we were able to produce benchmark results for the Stokes parameters both for reflected and transmitted light in the cases of molecular, aerosol and cloudy multiply scattering media. It was assumed that the single scattering albedo is equal to one. Benchmark results have been provided by several studies before, including Coulson et al. [22], Garcia and Siewert [7], [8], Wauben and Hovenier [10], and Natraj et al. [11] among others. However, the case of the elongated phase functions such as for a cloud and with a high angular resolution is presented here for the first time. Also in difference with other studies, we make inter-comparisons using several codes for the same input dataset, which enables us to quantify the corresponding errors more accurately.     

Radiative and conductive heat transfer in a nongrey semitransparent medium. Application to fire protection curtains

This paper deals with heat transfer in nongrey media which scatter, absorb and emit radiation. Considering a two dimensional geometry, radiative and conductive phenomena through the medium have been taken into account. The radiative part of the problem was solved using the discrete ordinate method with classical S_n quadratures. The absorption and scattering coefficients involved in the radiative transfer equation (RTE) were obtained from the Mie theory. Conduction inside the medium was linked to the RTE through the energy conservation. Validation of the model has been achieved with several simulation of water spray curtains used as fire protection walls.     

Numerical investigation of the effect of soot aggregation on the radiative properties in the infrared region and radiative heat transfer

The effect of aggregation on soot radiative properties in the infrared region of the spectrum is numerically investigated using Rayleigh-Debye-Gans theory for fractal aggregates (RDG-FA). In order to use the RDG-FA theory for a wide range of aggregate sizes and wavelengths, the predicted phase functions, scattering and absorption coefficients are compared with a more accurate theory, the integral equation formulation for scattering—IEFS. The importance of scattering when compared with absorption is investigated, as well as the effect of aggregation on the phase function shape and on the scattering cross section. It is concluded that in the case of small aggregates formed with small primary particles the scattering coefficient is negligible compared with the absorption coefficient, and scattering and aggregation of primary particles can be ignored. Thus, the Rayleigh approximation can be used leading to isotropic scattering. In the case of large aggregates constituted by large primary particles, aggregation becomes important and the scattering cross section is of the same order of magnitude of the absorption cross section. Moreover, the phase function becomes highly peaked in the forward direction. Therefore, the Rayleigh and the equivalent volume Mie sphere approximations are not valid, and the RDG-FA method emerges as a good compromise between accuracy and simplicity of application. However, radiative transfer calculations between two infinite, parallel, black walls show that scattering may always be neglected in the calculation of total radiative heat source and heat fluxes to the walls. The minor influence of scattering on the accuracy of the predictions is explained by the shift between the spectral region where scattering is important and the region where the spectral radiative heat source is large.     

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.     

Fast yet accurate computation of the complete radiance distribution in the coupled atmosphere-ocean system

The iteration of the source-function technique is used in conjunction with the discrete ordinate radiative transfer (DISORT) method to derive analytic solutions for the radiances at arbitrary polar angles and arbitrary levels in a multiple scattering (random) medium consisting of two strata with a change in the index of refraction across the interface between them. The atmosphere overlying a body of water (the ocean) serves as a prime example of such a system consisting of two adjacent strata with different indices of refraction. The accuracy of the solutions is tested by considering the coupled atmosphere-ocean system where absorption and scattering by molecules and particles are included. The tests show that these analytic solutions are superior to a spline interpolation that has been used previously in connection with the DISORT method applied to a random medium consisting of two strata with different indices of refraction.     

Analysis of collimated radiation in participating media using the discrete transfer method

A general formulation of the discrete transfer method is provided to analyze radiative heat transfer problems in a participating medium subjected to collimated radiation. The formulation is validated by considering 1-D planar absorbing, emitting and anisotropically scattering gray medium in radiative equilibrium. Anisotropy of the medium is approximated by linear anisotropic phase function. For the purpose of comparison, the problem is also solved analytically. Results are obtained for different angles of incidence of the collimated radiation. At a given angle of incidence, results are obtained for forward, isotropic and backward scattering situations. Heat flux results are compared over a wide range of values of the extinction coefficient. Emissive power distributions in the medium are also obtained for some cases. The discrete transfer method results are found to compare very well with the analytic results.     

Adjoint problem of radiative transfer for a pseudo-spherical atmosphere and general viewing geometries

We formulate the adjoint radiative transfer for a pseudo-spherical atmosphere and various retrieval scenarios. The single scattering radiance is computed in a spherical atmosphere by using the source integration technique, while for the multiple scattering radiance we formulate an one-dimensional adjoint radiative transfer equation in a plane-parallel atmosphere. The adjoint solution of the radiative transfer equation is obtained by employing the discrete ordinate method with matrix exponential. We provide an abbreviated derivation of our formalism as well as a discussion of the numerical implementation of the theory.     

紫外波段多分散系气溶胶散射相函数随机抽样方法研究

基于电磁散射与辐射传输中的基本理论,对紫外波段霾尺度范围内满足特定分布的多种气溶胶粒子的散射相函数进行了研究.提出了一种直接随机抽样拟合散射相函数的方法.比较了H-G相函数、改进的H-G相函数及随机抽样拟合的相函数与多分散系Mie相函数的偏离程度.数值计算了不同相函数拟合方法对应气溶胶的传输特性.计算结果表明,相函数的准确模拟计算对于蒙特卡罗方法等辐射传输问题的解决具有十分重要的意义.     

Temperature field of radiative equilibrium in a two-dimensional graded index medium with gray boundaries

In this paper, a numerical method is presented for the study of the radiative transfer in a two-dimensional graded index semitransparent medium with diffuse gray boundaries. The numerical method is a combination of the linear refractive index bar model, the discrete curved ray-tracing technique and the pseudo source adding method (LRIB-CRTP). In the traditional ray-tracing technique, it is difficult to deal with the diffuse gray boundary while solving the radiative transfer. Using the pseudo source adding method, the diffuse gray boundary of the medium can be treated as a black boundary. We have also studied the radiative equilibrium temperature field of the medium and analyzed the influence of some parameters involved. The results show that the directional discrete number is important for the medium having a large absorption coefficient. The results also show that the refractive index distribution greatly influences the temperature field, whereas the linear absorption coefficient distribution has little influence on the temperature field.     

Probe into the reflection from GaP nanoparticles via different solutions of radiative transfer equation

The reflection and absorption spectra of gallium phosphide (GaP) nanoparticles were measured. The radiative transfer equation (RTE) for the medium with scattering and absorption is solved by three different solutions. The ratio of the absorption and scattering coefficients (E _a/E _s) of the GaP nanoparticles layer is calculated from the reflection spectrum via the three solutions, respectively, and the result derived with the three-flux model is closest to the exact solution given by Giovanelli. The E _a/E _s curves all exhibit the energy band gaps of GaP nanoparticles, which are consistent with the absorption spectrum measurement. The shape of the reflection spectrum is mainly determined by the absorption, and the scattering only influences its intensity. The energy band structure of the powder sample plays an important role in the reflection phenomenon, and the reflectance data can be used for quantitative analyses.     

三维燃烧介质和壁面温度的非接触联合重建研究

温度分布在线实时测量对于燃烧过程优化和污染物控制具有重要意义, 针对以往非接触三维温度分布重建过程的耗时性问题和忽略壁面辐射的不足, 本文提出了一种新的离散重建模型, 用于三维吸收、 发射和散射性高温燃烧介质以及壁面温度的快速联合非接触测量. 该模型以四个CCD(Charge Coupled Device) 为测量传感器, 通过构建辐射逆问题求解方程, 从CCD输出的辐射投影图像重建温度分布. 介质中不同投影方向内的辐射传递过程通过离散传递法来描述, 介质的散射和壁面反射则通过离散坐标法来近似. 离散后计算局部辐射强度的病态方程通过最小二乘余量法来求解, 论文对其计算速度进行了优化. 通过非对称温度分布测量算例分析了该模型的有效性, 讨论了测量噪音、 介质和壁面辐射特性对重建精度的影响, 并与其他方法对比分析了模型的重建速度. 计算结果表明本文提出的离散模型可以有效地用于大型高温燃烧介质和壁面温度分布的联合非接触测量. 即使在有噪声的情况下, 该模型也能获得准确的测量结果, 与其他计算方法相比, 采用改进的最小二乘余量法, 能有效地提高温度分布的重建计算速度.     

Generalized form of the direct and adjoint radiative transfer equations

The main goal of this paper is to give a rigorous derivation of the generalized form of the direct (also referenced as forward) and adjoint radiative transfer equations. The obtained expressions coincide with expressions derived by Ustinov [Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211]. However, in contrast to [Ustinov EA. Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211] we formulate the generalized form of the direct radiative transfer operator fully independent from its adjoint. To illustrate the application of the derived adjoint radiative transfer operator we consider the angular interpolation problem in the framework of the discrete ordinate method widely used to solve the radiative transfer equation. It is shown that under certain conditions the usage of the solution of the adjoint radiative transfer equation for the angular interpolation of the intensity can be computationally more efficient than the commonly used source function integration technique.     

On the finite volume method and the discrete ordinates method regarding radiative heat transfer in acute forward anisotropic scattering media

The ability of the finite volume method (FVM) and the discrete ordinates method (DOM) to model radiative heat transfer in acute forward anisotropic scattering media has been investigated. The test case involves a purely scattering medium in a cubic enclosure, irradiated by one boundary with diffuse emission. Four phase functions have been considered: three of the Henyey-Greenstein type with respective asymmetry factors of 0.2, 0.8 and 0.93, and a Mie phase function with a strong forward scattering peak (computed for a size parameter of 245 and corresponding to an asymmetry factor of 0.93). Results obtained with the FVM are in good agreement with Monte Carlo reference solutions, whatever the level of acute anisotropic scattering (for asymmetry factors up to 0.93). The DOM combined with the renormalization procedures of the phase function proposed by Kim and Lee (Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular enclosures. Int J Heat Mass Transfer 1988;31:1711-21. [1]) and Wiscombe (On initialization error and flux conservation in the doubling method. JQSRT 1976;18:637-58. [2]) provides accurate results only for the smallest asymmetry factor. As the asymmetry factor increases, the renormalization procedures induce strong modifications in the values of the discretized phase function resulting in an underestimation of the effective attenuation by scattering. This error has been found to increase with optical thickness. In fact, when using the DOM, results would be more accurate combining this method with a Delta-Eddington approximation of the phase function, instead of using the actual phase function which is altered too much by renormalization.     

Influence of complex component and particle polydispersity on radiative properties of soot aggregate in atmosphere

The effects of morphological structure, water coating, dust mixing and primary particle size distribution on the radiative properties of soot fractal aggregates in atmosphere are investigated using T-matrix method. These fractal aggregates are numerically generated using a combination of the particle-cluster and cluster-cluster aggregation algorithms with fractal parameters representing soot aggregate in atmosphere. The radiative properties of compact aggregate notably deviate from that of the branched one, and the effect of morphology changes on the radiative properties in wet air cannot be neglected. However it is reasonable to use realization-averaged radiative properties to represent that of the aggregates with certain morphology. In wet air, the scattering, absorption and extinction cross-section and symmetry parameter of soot aggregates coated with water notably increase with water shell thickness. The mixing structures of dust have little effect on radiative properties of aggregates, but the volume fraction of dust has an obvious effect on extinction, scattering and absorption cross-section of aggregates when the size parameters are above the Rayleigh limit. Although the primary particle size distribution of soot aggregate has mild effect on the scattering albedo and asymmetry parameter, the deviations of the extinction, scattering, absorption cross-section among the three size distributions are significant in this study. The size distribution has a significant effect on forward scattering of phase function, while the effect can be neglected as the size parameter approaches to the Rayleigh limit.     

Diffuse reflectance of TiO2 pigmented paints: Spectral dependence of the average pathlength parameter and the forward scattering ratio

Diffuse reflectance spectra of paint coatings with different pigment concentrations, normally illuminated with unpolarized radiation, have been measured. A four-flux radiative transfer approach is used to model the diffuse reflectance of TiO₂ (rutile) pigmented coatings through the solar spectral range. The spectral dependence of the average pathlength parameter and of the forward scattering ratio for diffuse radiation, are explicitly incorporated into this four-flux model from two novel approximations. The size distribution of the pigments has been taken into account to obtain the averages of the four-flux parameters: scattering and absorption cross sections, forward scattering ratios for collimated and isotropic diffuse radiation, and coefficients involved in the expansion of the single particle phase function in terms of Legendre polynomials.     

求解辐射传递方程的离散坐标法

阐述了含吸收散射性介质三维空腔内辐射传递方程的离散坐标解法。讨论了入射散射项积分格式的选取,以及假散射和射线效应对解精度的影响。对三维矩形炉膛内辐射传递过程进行了数值模拟,并与区域法和离散传递法进行比较。比较结果表明离散坐标法具有较好的精度,是目前燃烧室内辐射传热过程数值模拟的一种较好的方法。     

Approximate analytical solution to normal emittance of semi-transparent layer of an absorbing, scattering, and refracting medium

A two-step approximate analytical solution for the normal emittance of a plane layer of an absorbing, scattering and refracting medium is derived analytically. The analysis is based on the transport approximation and the two-step solution method for radiative transfer. The high accuracy of the approximate solution, examined by comparing its results to those obtained independently by the discrete ordinates and Monte Carlo methods, makes it suitable for application in combined experimental-analytical studies to identify selected spectral radiative properties of dispersed media in the range of semi-transparency.