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.     

Evaluation of solution methods for radiative heat transfer in gaseous oxy-fuel combustion environments

The exact solution to radiative heat transfer in combusting flows is not possible analytically due to the complex nature of the integro-differential radiative transfer equation (RTE). Many different approximate solution methods for the solution of the RTE in multi-dimensional problems are available. In this paper, two of the principal methods, the spherical harmonics (P₁) and the discrete ordinates method (DOM) are used to calculate radiation. The radiative properties of the gases are calculated using a non-gray gas full spectrum k-distribution method and a gray method. Analysis of the effects of numerical quadrature in the DOM and its effect on computation time is performed. Results of different radiative property methods are compared with benchmark statistical narrow band (SNB) data for both cases that simulate air combustion and oxy-fuel combustion. For both cases, results of the non-gray full spectrum k-distribution method are in good agreement with the SNB data. In the case of oxy-fuel simulations with high partial pressures of carbon dioxide, use of gray method for the radiative properties may cause errors and should be avoided.     

Numerical prediction of heat transfer by natural convection and radiation in an enclosure filled with an isotropic scattering medium

This paper deals with the numerical solution for natural convection and volumetric radiation in an isotropic scattering medium within a heated square cavity using a hybrid thermal lattice Boltzmann method (HTLBM). The multiple relaxation time lattice Boltzmann method (MRT-LBM) has been coupled to the finite difference method (FDM) to solve momentum and energy equations, while the discrete ordinates method (DOM) has been adopted to solve the radiative transfer equation (RTE) using the S8 quadrature. Based on these approaches, the effects of various influencing parameters such as the Rayleigh number (Ra), the wall emissivity (ε_ι), the Planck number (Pl), and the scattering albedo (ω), have been considered. The results presented in terms of isotherms, streamlines and averaged Nusselt number, show that in absence of radiation, the temperature and the flow fields are centro-symmetrics and the cavity core is thermally stratified. However, radiation causes an overall increase in the temperature and velocity gradients along both thermally active walls. The maximum heat transfer rate is obtained when the surfaces of the enclosure walls are regarded as blackbodies. It is also seen that the scattering medium can generate a multicellular flow.     

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.     

A note on the interface functions

The functions B₁(τ), D₂(τ), and D₃(τ), which arise in the study of radiative transfer through a planar medium with a relative refractive index different from unity, have been investigated. These functions reduce to exponential integrals when the refractive index of the medium is equal to that of the surroundings. Their behavior has been examined and is presented graphically and in tabular form for various ratios of the interface refractive index. When the interface refractive index ratio is larger than 1.2, these functions can be accurately approximated by an exponential of the form a e^?bτ.     

Fast approximate technique for the cumulative wavenumber model to modeling radiative transfer in a mixture of real gas media

In the cumulative wavenumber (CW) model, the total range of the absorption cross-section Cη is subdivided into the supplementary absorption cross-section of gray gases C_j, j=1,…,n, where n is the number of gray gases; and the wavenumber region is subdivided into intervals Δ_i=[η_i−1, η_i], i=1, 2,…,p, where p is the number of intervals. The intersection of the two spectral subdivisions is used to define the modeling of the fractional gray gas D_ij. In the CW model, we solve the radiative transfer equation (RTE) in every subinterval D_ij; then it is necessary to solve n x p times the spectral form of the RTE for complete spectral integration. In this work, the CW model is used with a numerical approximation technique based on additive properties of radiative intensity to reduce the solution of RTE to n new fractional gray gas D_j for complete spectral integration. The CW model was first coupled with the discrete ordinates method and the accuracy of the simplified technique and the algorithm was first examined for one-dimensional homogeneous media; results are compared with line-by-line calculations and it is found that the CW model with the simplified technique is exact for the homogeneous media examined. Also, the fast approach is tested in the diffuse reflecting boundaries case. The CW model is implemented in a bi-dimensional enclosure containing real gases in isothermal cases. Afterwards, this approximate technique is extended to non-isothermal and non-homogeneous cases; the results are compared with line-by-line calculations taken from literature and good agreement was found. The results obtained using the acceleration technique for the CW model agree with the results of original CW model. With this acceleration technique the CPU time decreases p times. Spectral database HITRAN and HITEMP are used to obtain the molecular absorption spectrum of the gases.     

The use of the itFN method for radiative transfer problems with reflective boundary conditions

The F_N method is used to compute the net radiative heat flux relevant to radiative transfer in an anisotropically scattering, plane-parallel medium with specularly and diffusely reflecting boundaries.     

Elliptic PDE formulation and boundary conditions of the spherical harmonics method of arbitrary order for general three-dimensional geometries

The inherent complexity of the radiative transfer equation makes the exact treatment of radiative heat transfer impossible even for idealized situations and simple boundary conditions. Therefore, a wide variety of efficient solution methods have been developed for the RTE. Among these solution methods the spherical harmonics method, the moment method, and the discrete ordinates method provide means to obtain higher-order approximate solutions to the equation of radiative transfer. Although the assembly of the governing equations for the spherical harmonics method requires tedious algebra, their final form promises great accuracy for any given order, since it is a spectral method (rather than finite difference/finite volume in the case of discrete ordinates). In this study, a new methodology outlined in a previous paper on the spherical harmonics method (P_N) is further developed. The new methodology employs successive elimination of spherical harmonic tensors, thus reducing the number of first-order partial differential equations needed to be solved simultaneously by previous P_N approximations (=(N+1)²). The result is a relatively small set (=N(N+1)/2) of second-order, elliptic partial differential equations, which can be solved with standard PDE solution packages. General boundary conditions and supplementary conditions using rotation of spherical harmonics in terms of local coordinates are formulated for the general P_N approximation for arbitrary three-dimensional geometries. Accuracy of the P_N approximation can be further improved by applying the “modified differential approximation” approach first developed for the P₁-approximation. Numerical computations are carried out with the P₃ approximation for several new two-dimensional problems with emitting, absorbing, and scattering media. Results are compared to Monte Carlo solutions and discrete ordinates simulations and a discussion of ray effects and false scattering is provided.     

Comparison of radiative transfer approximations for a highly forward scattering planar medium

We examine critically the accuracy of the two-flux, spherical harmonics and discrete ordinates methods for predicting radiative transfer in a planar, highly-forward scattering and absorbing medium. Numerical results for the radiative fluxes show that the two-flux and P₃-approximations yield accurate results compared to solutions based on the F_N-method. Indeed, these approximate methods are relatively simple and have potential for generalization to predict radiative transfer in multidimensional systems, as long as an appropriate simplification of the phase function is utilized.     

Application of the modified differential approximation for radiative transfer to arbitrary geometry

The first-order spherical harmonics method (or P₁ approximation) has found prolific usage for approximate solution of the radiative transfer equation (RTE) in participating media. However, the accuracy of the P₁ approximation deteriorates as the optical thickness of the medium is decreased. The modified differential approximation (MDA) was originally proposed to remove the shortcomings of the P₁ approximation in optically thin situations. This article presents algorithms to apply the MDA to arbitrary geometry—in particular, geometry with obstructions, and inhomogeneous media. The wall-emitted component of the intensity was computed using a combined view-factor and ray-tracing approach. The Helmholtz equation, arising out of the medium-emitted component, was solved using an unstructured finite-volume procedure. The general procedure was validated for both two-dimensional (2D) and three-dimensional (3D) geometries against benchmark Monte Carlo results. The accuracy of MDA was found to be superior to the P₁ approximation for all optical thicknesses. Its accuracy, when compared with the discrete ordinates method (both S₆ and S₈), was found to be clearly superior in optically thin situations, but problem dependent in optically intermediate and thick situations. For 3D geometries, calculation and storage of the view-factor matrix was found to be a major shortcoming of the MDA. In addition, for inhomogeneous media, calculation of optical distances requires a ray-tracing procedure, which was found to be a bottleneck from a computational efficiency standpoint. Several strategies to reduce both memory and computational time are discussed and demonstrated.     

Time-dependent radiative transfer using Chapman-Enskog maximum entropy method

The time-dependent problems of radiative transfer involve a coupling between radiation and material energy fields and are nonlinear because of proposed temperature dependence of the medium characteristics in semi-infinite medium with Rayleigh anisotropic scattering. By means of the limited flux, Chapman-Enskog and maximum entropy technique the time-dependent radiative transfer equation has been solved explicitly. The maximum entropy method is used to solve the resulting differential equation for radiative energy density. The calculations are carried out for temperature (normalized dimensionless) Θ(x,τ), radiative energy density and net flux with Rayleigh and anisotropic scattering for different space at different times.     

Spectral collocation method with a flexible angular discretization scheme for radiative transfer in multi-layer graded index medium

The spectral collocation method (SCM) is employed to solve the radiative transfer in multi-layer semitransparent medium with graded index. A new flexible angular discretization scheme is employed to discretize the solid angle domain freely to overcome the limit of the number of discrete radiative direction when adopting traditional S_N discrete ordinate scheme. Three radial basis function interpolation approaches, named as multi-quadric (MQ), inverse multi-quadric (IMQ) and inverse quadratic (IQ) interpolation, are employed to couple the radiative intensity at the interface between two adjacent layers and numerical experiments show that MQ interpolation has the highest accuracy and best stability. Variable radiative transfer problems in double-layer semitransparent media with different thermophysical properties are investigated and the influence of these thermophysical properties on the radiative transfer procedure in double-layer semitransparent media is also analyzed. All the simulated results show that the present SCM with the new angular discretization scheme can predict the radiative transfer in multi-layer semitransparent medium with graded index efficiently and accurately.     

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.     

A numerical model for multigroup radiation hydrodynamics

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M₁) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well-established kinetic model through Marshak wave tests with frequency-dependent opacities. Very good agreement between the multigroup M₁ and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.     

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.     

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.     

Differential approximations for transient radiative transfer in refractive planar media with pulse irradiation

Transient radiative transfer in an anisotropically scattering refractive planar medium with pulse irradiation is solved by various approximation methods, such as P?1, P?1 parabolic, P_1/3 and two-flux. The time-resolved transmittance and reflectance are calculated for various radiative parameters, and are compared with those obtained by the discrete ordinate method (DOM). Among the approximation methods considered, the P_1/3 approximation is the better one, because its results are in overall good agreement with those obtained by the more rigorous DOM, except the transmittance around the peak for neither thin nor very thick slabs. It is found that the curved paths of radiation and the internal reflection of the back scattered radiation enhance the effect of scattering.     

Application of Synthetic Kernel (SKN) method to participating linearly anisotropically scattering solid spherical medium

The Synthetic Kernel (SK_N) method is applied to a solid spherical absorbing, emitting and linearly anisotropically scattering homogeneous and inhomogeneous medium. The SK_N method relies on approximating the integral transfer kernels by Synthetic Kernels. The radiative integral transfer equation is then reducible to a set of coupled second-order differential equations. The SK_N method, which uses Gauss quadratures, is tested against integral equation and the discrete-ordinates S₈ solutions for various optical radius and scattering albedo variations.     

Determination of optical properties of slices of turbid media by diffuse CW laser light scattering profilometry

In this work we present a method for determining the optical parameters of turbid media, namely its absorption coefficient (μ_a) and its reduced scattering coefficient . It is based on the measurement of CW transmittance profiles and analysis of the experimental data by a theoretical model based on the diffusion approximation (DA) of the radiative transfer equation (RTE). The method developed has been investigated with solid polymer probes but it could be applied for liquid materials as well. Experimental results are presented and compared to those of other authors together with a discussion about the accuracy of measurements. In addition, measurements using integrating spheres as well as Monte Carlo simulations are also presented to validate these results.     

Simultaneous solution of Kompaneets equation and radiative transfer equation in the photon energy range 1-125 keV

Radiative transfer equation in plane parallel geometry and Kompaneets equation is solved simultaneously to obtain theoretical spectrum of 1-125 keV photon energy range. Diffuse radiation field are calculated using time-independent radiative transfer equation in plane parallel geometry, which is developed using discrete space theory (DST) of radiative transfer in a homogeneous medium for different optical depths. We assumed free-free emission and absorption and emission due to electron gas to be operating in the medium. The three terms n, n² and (∂n/∂x_k) where n is photon phase density and x_k=(hν/kT_e), in Kompaneets equation and those due to free-free emission are utilized to calculate the change in the photon phase density in a hot electron gas. Two types of incident radiation are considered: (1) isotropic radiation with the modified black body radiation I^MB[1] and (2) anisotropic radiation which is angle dependent. The emergent radiation at τ=0 and reflected radiation τ=τ_max are calculated by using the diffuse radiation from the medium. The emergent and reflected radiation contain the free-free emission and emission from the hot electron gas. Kompaneets equation gives the changes in photon phase densities in different types of media. Although the initial spectrum is angle dependent, the Kompaneets equation gives a spectrum which is angle independent after several Compton scattering times.