Radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and particles using the modified discrete-ordinates method

The radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and a cloud of particles is investigated by using the combined finite-volume and discrete-ordinates method, named modified discrete-ordinates method (MDOM), which integrates the radiative transfer equation (RTE) over a control volume and a control angle simultaneously like in the finite-volume method (FVM) and treats the angular derivative terms due to spherical geometry as the conventional discrete-ordinates method (DOM). The radiative properties involving non-gray gas and particle behavior are modeled by using the extended weighted sum of gray gases model (WSGGM) with particles. Mathematical formulation and final discretization equations for the RTE are introduced by considering the behavior of a two-phase mixture of non-gray gas and particles in a spherically symmetric concentric enclosure. The present approach is validated by comparing with the results of previous works including gray and non-gray radiative heat transfer. Finally, a detailed investigation of the radiative heat transfer with non-gray gases and/or a two-phase mixture is conducted to examine the dependence of the radiative heat transfer upon temperature ratio between inner and outer spherical enclosure, particle concentration, and particle temperature.     

Light Reflection and Transmission by layers with oriented spheroidal particles

On the basis of the numerical solution of the radiative transfer equation (RTE), the influence of non-sphericity of oriented spheroidal particles on radiation intensity formation in a layer with a multiple scattering regime is considered. For solving the RTE with characteristics depending on the propagation direction of the radiation beam, an algorithm is developed based on the method of layer doubling. Calculations are presented of angle correlations for intensities of radiation in a layer with Fresnel reflection on the boundaries, and reflection and transmission coefficients.     

涂层的热红外发射率计算模型

根据Kubelka-Munk 模型和Mie散射理论,解出了涂覆在高发射率衬底上的红外涂层里的辐射传输方程（RTE）,建立了涂层在热红外（8～14 μm）波段的发射率的计算模型.在运用铝粉粒子作为颜料粒子的情况下,根据计算结果,得出了涂层发射率与粒子半径、涂层厚度的关系,以及最优的粒子半径范围.     

A theoretical study of radiation between two concentric spheres using a modified discrete ordinates method associated with Legendre transform

A modified discrete ordinates method (DOM) is used in spherical participating media. The radiative intensity is broken up into two components. One component is traced back to the enclosure's source. It is called direct intensity. The other component is rather traced back to the contribution of the medium itself. It is called diffuse intensity. Thus, the radiative transfer equation (RTE) is transformed into two simultaneous equations: a direct RTE and a diffuse RTE. The direct RTE is solved analytically. The diffuse RTE is solved numerically using the DOM. The streaming angular derivative term appearing in spherical geometry is modeled by making use of the Finite Legendre Transform. We study a pure radiation transfer problem between two concentric spheres. The medium is assumed to be gray and isotropically scattering. The limiting spheres are considered to be opaque, gray, diffusely emitting and diffusely reflecting with uniform emissivity over each surface. The obtained results are compared with available cases reported in the literature. In particular, relative importance of the direct radiation in optically thin media is studied.     

利用T-matrix计算非球形粒子散射特性的研究

本文用T-matrix方法计算了非球形气溶胶粒子的光学特性,得到了气溶胶粒子的消光截面、散射截面、吸收截面与气溶胶粒子形状的关系,不同形状气溶胶粒子的有着相同的散射相函数和不同的偏振度,非球形气溶胶粒子的散射相函数对其复折射指数的实部和虚部都不太敏感,而偏振相函数对其实部和虚部都比较敏感.此结论为研究大气辐射传输提供了较好的方法,尤其是偏振度与偏振相函数的提出为用偏振的方法进一步的反演气溶胶的光学参数提供了理论基础.     

Radiative properties of densely packed spheres in semitransparent media: A new geometric optics approach

This contribution presents a new Ray-tracing method for calculating effective radiative properties of densely packed spheres in non-absorbing or semitransparent host medium. The method is restricted to the geometric optic objects and neglects the wave effects. The effective radiative properties such as the absorption and scattering coefficients, and phase function are retrieved from the calculation of mean-free paths of scattering and absorption, and the angular scattering probability of radiation propagating in the dispersed medium. The model accounts for the two geometric effects called here as non-point scattering and ray transportation effects. The successful comparison of the current model with data of radiative properties and transmittances of particle beds in a non-absorbing medium reported in the literature confirm its suitability. It is shown that: (i) for opaque or absorbing particles (not systematically opaque), the non-point scattering is the dominant geometric effects whereas both non-point scattering and ray transportation effects occur for weakly absorbing and transparent particles. In the later cases, these two geometric effects oppose and may cancel out. This may explain why the Independent scattering theory works well for packed of quasi-transparent particles; (ii) the non-point scattering and ray transportation effects can be captured through the scattering and absorption coefficients while using the classical form of phase function. This enables using the standard radiative transfer equation (RTE); (iii) the surrounding medium absorption can be accounted for without any homogenization rule. It contributes to increasing the effective absorption coefficient of the composite medium as expected but, at the same time, it reduces the particle extinction; and (iv) the current transfer calculation predicts remarkably the results of direct Monte Carlo (MC) simulation. This study tends therefore to confirm that the RTE can be applied to densely packed media by using effective radiative properties.     

Comparison of the radiative transfer equations for constant and spatially varying refraction

The radiative transfer equation for scattering media with constant refraction index (RTE) and the radiative transfer equation for scattering media with spatially varying refraction index (RTEvri) are compared by using the principle of conservation of energy. It is shown that the RTEvri, not only accounts for the spatial variations of refraction index, but also contains a term that accounts for the divergence of the rays. The latter term is missing in the RTE. A corrected RTE is proposed.     

Directional radiometry and radiative transfer: A new paradigm

Measurements with directional radiometers and calculations based on the radiative transfer equation (RTE) have been at the very heart of weather and climate modeling and terrestrial remote sensing. The quantification of the energy budget of the Earth's climate system requires exquisite measurements and computations of the incoming and outgoing electromagnetic energy, while global characterization of climate system's components relies heavily on theoretical inversions of observational data obtained with various passive and active instruments. The same basic problems involving electromagnetic energy transport and its use for diagnostic and characterization purposes are encountered in numerous other areas of science, biomedicine, and engineering. Yet both the discipline of directional radiometry and the radiative transfer theory (RTT) have traditionally been based on phenomenological concepts many of which turn out to be profound misconceptions. Contrary to the widespread belief, a collimated radiometer does not, in general, measure the flow of electromagnetic energy along its optical axis, while the specific intensity does not quantify the amount of electromagnetic energy transported in a given direction.The recently developed microphysical approach to radiative transfer and directional radiometry is explicitly based on the Maxwell equations and clarifies the physical nature of measurements with collimated radiometers and the actual content of the RTE. It reveals that the specific intensity has no fundamental physical meaning besides being a mathematical solution of the RTE, while the RTE itself is nothing more than an intermediate auxiliary equation. Only under special circumstances detailed in this review can the solution of the RTE be used to compute the time-averaged local Poynting vector as well as be measured by a collimated radiometer. These firmly established facts make the combination of the RTE and a collimated radiometer useful in a well-defined range of applications. However, outside the domain of validity of the RTT the practical usefulness of measurements with collimated radiometers remains uncertain, while the theoretical modeling of these measurements and the solution of the energy-budget problem require a more sophisticated approach than solving the RTE.     

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.     

Accurate and efficient Monte Carlo solutions to the radiative transport equation in the spatial frequency domain

We present an approach to solving the radiative transport equation (RTE) for layered media in the spatial frequency domain (SFD) using Monte Carlo (MC) simulations. This is done by obtaining a complex photon weight from analysis of the Fourier transform of the RTE. We also develop a modified shortcut method that enables a single MC simulation to efficiently provide RTE solutions in the SFD for any number of spatial frequencies. We provide comparisons between the modified shortcut method and conventional discrete transform methods for SFD reflectance. Further results for oblique illumination illustrate the potential diagnostic utility of the SFD phase-shifts for analysis of layered media.     

Dynamic control of the terahertz rainbow trapping effect based on a silicon-filled graded grating

We theoretically propose a scheme to realize the dynamic control of the properties of the terahertz(THz) rainbow trapping effect(RTE) based on a silicon-filled graded grating(SFGG) in a relatively broad band via optical pumping.Through the theoretical analysis and finite-element method simulations, it is conceptually demonstrated that the band of the RTE can be dynamically tuned in a range of ~0.06 THz. Furthermore, the SFGG can also be optically switched between a device for the RTE and a waveguide for releasing the trapped waves. The results obtained here may imply applications for the tunable THz plasmonic devices, such as on-chip optical buffers, broad band slow-light systems, and integrated optical filters.     

A quasi-analytic solution of the radiative transfer equation for three-dimensional atmospheres

In this study, we present a new solution of the three-dimensional (3-D) radiation transfer equation (RTE). The solution employs a discretization technique to separate the independent variables involved in the 3-D RTE, and the doubling-adding method to solve the RTE explicitly and quasi-analytically. The remarkable feature of the present solution is the application of scaling-function expansion to those terms that are dependent on horizontal coordinates. Scaling-function expansion is suitable for representing irregular horizontal inhomogeneities with small-scale variations. By applying scaling-function expansion, the 3-D RTE can be formulated in the form of a vector-matrix differential equation; matrices involved in the equation are generally sparse and dominantly diagonal matrices, and this considerably reduces the labor involved in matrix calculations. We tested the performance of the present solution via radiative transfer calculations of solar radiation in horizontally inhomogeneous two-dimensional cloud models. The calculated results indicate that even if the resolution of the scaling-function expansion is too coarse in regions around small-scale variations, the influence does not spread problematically to other regions far from the variations; this illustrates the advantage of the scaling-function expansion. The present solution can be used to investigate quantitatively and to estimate the effects of cloud spatial inhomogeneity on the corresponding radiation field.     

Refractive index imaging from radiative transfer equation-based reconstruction algorithm: Fundamentals

We present the first reconstruction algorithm for refractive index imaging, which is based on the radiative transfer equation (RTE). An objective function is iteratively minimized to find a solution to the problem of inversion of the refractive index field. The function describes the discrepancies of the emerging light measurements on the surface of the sample to be probed with predicted data from the corresponding numerical model. The unknown refractive index field is updated within each reconstruction iteration according to a search direction on the index distribution given by the adjoint model to the RTE. In this paper, emphasis is placed on the theoretical aspects. Preliminary tests are demonstrated on generic phantoms.     

Dielectronic recombination and resonant transfer excitation processes for helium-like krypton

The relativistic configuration interaction method is employed to calculate the dielectronic recombination(DR) cross sections of helium-like krypton via the 1s2lnl ’(n = 2,3,...,15) resonances.Then,the resonant transfer excitation(RTE) processes of Kr 34+ colliding with H,He,H 2,and CH x(x = 0-4) targets are investigated under the impulse approximation.The needed Compton profiles of targets are obtained from the Hartree-Fock wave functions.The RTE cross sections are strongly dependent on DR resonant energies and strengths,and the electron momentum distributions of the target.For H 2 and H targets,the ratio of their RTE cross sections changes from 1.85 for the 1s2l2l ’ to 1.88 for other resonances,which demonstrates the weak molecular effects on the Compton profiles of H 2.For CH x(x = 0-4) targets,the main contribution to the RTE cross section comes from the carbon atom since carbon carries 6 electrons;as the number of hydrogen increases in CH x,the RTE cross section almost increases by the same value,displaying the strong separate atom character for the hydrogen.However,further comparison of the individual orbital contributions of C(2p,2s,1s) and CH 4(1t 2,2a 1,1a 1) to the RTE cross sections shows that the molecular effects induce differences of about 25.1%,19.9%,and 0.2% between 2p-1t 2,2s-2a 1,and 1s-1a 1 orbitals,respectively.     

Multilayer modeling of radiative transfer by SLW and CW methods in non-isothermal gaseous medium

Non-isothermal gaseous medium is modeled using the multilayer approach, breaking the one-dimensional system into a series of isothermal layers. Spectral integration of the radiative transfer equation (RTE) is performed for the spectral line weighted-sum-of-gray-gases and cumulative wavenumber approaches to modeling the spectral nature of the gas radiation. An exact analytical solution of RTE is obtained for the layers with both black and gray walls. Predictions show high accuracy, even with surprisingly few layers.     

Effects of the approximations of light propagation on quantitative photoacoustic tomography using two-dimensional photon diffusion equation and linearization

Quantitative photoacoustic tomography (QPAT) employing a light propagation model will play an important role in medical diagnoses by quantifying the concentration of hemoglobin or a contrast agent. However, QPAT by the light propagation model with the three-dimensional (3D) radiative transfer equation (RTE) requires a huge computational load in the iterative forward calculations involved in the updating process to reconstruct the absorption coefficient. The approximations of the light propagation improve the efficiency of the image reconstruction for the QPAT. In this study, we compared the 3D/two-dimensional (2D) photon diffusion equation (PDE) approximating 3D RTE with the Monte Carlo simulation based on 3D RTE. Then, the errors in a 2D PDE-based linearized image reconstruction caused by the approximations were quantitatively demonstrated and discussed in the numerical simulations. It was clearly observed that the approximations affected the reconstructed absorption coefficient. The 2D PDE-based linearized algorithm succeeded in the image reconstruction of the region with a large absorption coefficient in the 3D phantom. The value reconstructed in the phantom experiment agreed with that in the numerical simulation, so that it was validated that the numerical simulation of the image reconstruction predicted the relationship between the true absorption coefficient of the target in the 3D medium and the reconstructed value with the 2D PDE-based linearized algorithm. Moreover, the the true absorption coefficient in 3D medium was estimated from the 2D reconstructed image on the basis of the prediction by the numerical simulation. The estimation was successful in the phantom experiment, although some limitations were revealed.

     

Multiple scattering of light transmission in a smoke layer

The visibility in a fire scene decreases because of the existence of smoke produced by the flammable materials. With the growth of smoke concentration, the relationship between light and smoke becomes complicated due to the multiple scattering. In this paper, the radiative transfer equation (RTE) that considers the multiple scattering was applied to calculate the light transmission in a smoke layer. As input parameters of RTE, the single scattering albedo, asymmetry parameter and extinction cross section of single smoke agglomerate were calculated by the discrete dipole approximation (DDA) method. The effects of smoke agglomerate diameter, number density of smoke layer, and the incident light wavelength were considered. The results show that the light transmitted flux decreases with the growth of smoke diameter and number density, and increases with the growth of wavelength. The smoke diameter is dominant among the three parameters, and the light transmitted flux tends to be stable when the wavelength reaches a certain value.     

Coupled radiative and conductive heat transfer in a non-grey absorbing and emitting semitransparent media under collimated radiation

This paper deals with heat transfer in non-grey semitransparent two-dimensional sample. Considering an homogeneous purely absorbing medium, we calculated the temperature field and heat fluxes of a material irradiated under a specific direction. Coupled radiative and conductive heat transfer were considered. The radiative heat transfer equation (RTE) was solved using a S₈ quadrature and a discrete ordinate method. Reflection and absorption coefficients of the medium were calculated with the silica optical properties. The conduction inside the medium was linked to the RTE through the energy conservation. Validation of the model and two original cases are also presented.     

Modified finite volume method applied to radiative transfer in 2D complex geometries and graded index media

A modified finite volume method with unstructured triangular meshes is proposed to solve the RTE in 2D complex geometries and for graded index media. In such media, the RTE has an additional term corresponding to “angular redistribution”. This term is due to the change in the orientation of the direction of propagation for the radiation along curved optical paths. Some benchmark cases applied to a slab (1D) and a square cavity (2D) with linear and nonlinear refractive graded index are used to validate the new method. New results are presented for a disk with radial graded index.