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.     

Curved ray tracing method for one-dimensional radiative transfer in the linear-anisotropic scattering medium with graded index

The curved ray tracing method (CRT) is extended to radiative transfer in the linear-anisotropic scattering medium with graded index from non-scattering medium. In this paper, the CRT is presented to solve one-dimensional radiative transfer in the linear-anisotropic scattering gray medium with a linear refractive index and two black boundaries. The predicted temperature distributions and radiative heat flux at radiative equilibrium are determined by the proposed method, and numerical results are compared with the data in references. The results show that the CRT has a good accuracy for radiative transfer in the linear-anisotropic scattering medium with graded index and the dimensionless emissive power and dimensionless radiative heat flux depend on the dimensionless refractive index gradient. It can also be seen that the dimensionless refractive index gradient has important effects on the temperature discontinuity at the boundaries.     

Two analytic methods applied to radiative transfer in the linear-anisotropic scattering medium with graded index and diffuse gray boundaries

The curved ray-tracing method is extended to radiative transfer in the graded index medium with diffuse gray boundary conditions instead of black boundary conditions and the pseudo-source adding method is extended to the case of the linear-anisotropic scattering medium with graded index from non-scattering medium. Furthermore, the equivalence of the two methods is verified by formulation derivation. As exact analytical solutions, both the methods have high accuracy and fast computational speed. The predicted temperature distributions and dimensionless radiative heat flux at radiative equilibrium are determined by the proposed methods, and the numerical results are compared with the data in references. The results show that the present methods have a good accuracy. Influences of various combinations of refractive index and boundary emissivities on the temperature distributions and dimensionless radiative heat flux are also investigated.     

Thermal shield effect with uniaxial crystals

A model describing radiative transfers inside uniaxial anisotropic media is presented. The transport equations for each electromagnetic mode supported by these media are derived in the limit of geometrical optics and an analytical solution is obtained from a ray tracing method. The temperature field inside such a medium illuminated on both sides by a blackbody radiation is calculated and compared to the temperature field of an isotropic medium submitted to the same conditions. We show that the temperature field in the anisotropic medium is drastically smaller than its counterpart, even for weak anisotropy.     

Radiative equilibrium in a rectangular enclosure bounded by gray walls

Two-dimensional temperature and heat flux distributions are calculated for an absorbing-emitting gray medium at radiative equilibrium in a rectangular enclosure. The bounding walls are gray and diffuse with arbitrary surface temperature distributions, and heat generation may take place inside the medium. As a first approximation, the problem is solved for optically thick systems (differential approximation). These results are subsequently improved by the introduction of a number of geometrical parameters to yield good accuracy for all optical thicknesses. As examples, two cases are discussed in detail: (1) uniform heat generation in a black enclosure and (2) an enclosure with one gray surface at constant temperature. Comparison with some numerical solutions generated by Hottel's zonal method shows excellent agreement.     

The DRESOR method for radiative heat transfer in a one-dimensional medium with variable refractive index

This paper extends the DRESOR (Distribution of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface) method to radiative transfer in a variable refractive index medium. In this method, the intensity is obtained from the source term along the curved integration paths determined only by the variable refractive index, and the DRESOR values are calculated by the Monte Carlo method in which the propagation of the energy bundles are affected by Snell's law. With given temperatures on the black boundaries of a one-dimensional medium, the temperature distribution inside the medium with a variable scattering property is calculated under the condition of radiative equilibrium. It is shown that the DRESOR method has a good accuracy in the cases studied. For an isotropic-scattering medium with the same optical thickness, the scattering albedo has no effect on the temperature distribution, which can be obtained from the general equations and can be seen as an extension of what exists for a constant refractive index; however, the different refractive index causes obvious changes in the temperatures inside the medium. The effect of anisotropic scattering on the temperature distribution cannot be ignored, although it is still weaker than the effect caused by variation in the refractive index.     

多维梯度折射率介质内辐射传递的扩散近似

推导多维梯度折射率介质内稳态辐射传递的扩散近似方程.使用有限元法对扩散近似进行离散和求解,利用两个二维半透明介质的稳态辐射传递问题验证该扩散近似的精度及适用性.算例考虑介质为均匀折射率及梯度折射率两种情况.利用扩散近似分别求解辐射平衡时的边界热流、介质内温度场分布,并与辐射传递方程的求解结果进行对比分析.结果表明:介质折射率变化、散射特性、光学厚度及散射反照率均直接影响扩散近似的精度;在光学厚及强散射条件下,该扩散近似可以作为一种快速算法应用于梯度折射率介质稳态辐射传递的求解.     

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.     

Radiative transfer theory with time delay for the effect of pulse entrapping in a resonant random medium: General transfer equation and point-like scatterer model

Abstract

A pulse propagation of a vector electromagnetic wave field in a discrete random medium under the condition of Mie resonant scattering is considered on the basis of the Bethe–Salpeter equation in the two-frequency domain in the form of an exact kinetic equation which takes into account the energy accumulation inside scatterers. The kinetic equation is simplified using the transverse field and far wave zone approximations which give a new general tensor radiative transfer equation with strong time delay by resonant scattering. This new general radiative transfer equation, being specified in terms of the low-density limit and the resonant point-like scatterer model, takes the form of a new tensor radiative transfer equation with three Lorentzian time-delay kernels by resonant scattering. In contrast to the known phenomenological scalar Sobolev equation with one Lorentzian time-delay kernel, the derived radiative transfer equation does take into account effects of (i) the radiation polarization, (ii) the energy accumulation inside scatterers, (iii) the time delay in three terms, namely in terms with the Rayleigh phase tensor, the extinction coefficient and a coefficient of the energy accumulation inside scatterers, respectively (i.e. not only in a term with the Rayleigh phase tensor). It is worth noting that the derived radiative transfer equation is coordinated with Poynting's theorem for non-stationary radiation, unlike the Sobolev equation. The derived radiative transfer equation is applied to study the Compton–Milne effect of a pulse entrapping by its diffuse reflection from the semi-infinite random medium when the pulse, while propagating in the medium, spends most of its time inside scatterers. This specific albedo problem for the derived radiative transfer equation is resolved in scalar approximation using a version of the time-dependent invariance principle. In fact, the scattering function of the diffusely reflected pulse is expressed in terms of a generalized time-dependent Chandrasekhar H-function which satisfies a governing nonlinear integral equation. Simple analytic asymptotics are obtained for the scattering function of the front and the back parts of the diffusely reflected Dirac delta function incident pulse, depending on time, the angle of reflection, the mean free time, the microscopic time delay and a parameter of the energy accumulation inside scatterers. These asymptotics show quantitatively how the rate of increase of the front part and the rate of decrease of the rear part of the diffusely reflected pulse become slower with transition from the regime of conventional radiative transfer to that of pulse entrapping in the resonant random medium.     

Analysis of radiative transport in a cylindrical enclosure—An application of the modified discrete ordinate method

Application of the modified discrete ordinate method (MDOM) proposed by Mishra et al. [Mishra SC, Roy HK, Misra N. Discrete ordinate method with a new and simple quadrature scheme. J Quant Spectrosc Radiat Transfer 2006;101:249-262.] has been extended for calculation of volumetric radiative information in a cylindrical enclosure. Radiatively, the medium inside a diffuse gray 1-D concentric cylinder is absorbing, emitting and scattering. Three types of problems, viz., an isothermal medium representing non-radiative equilibrium case, a non-isothermal medium representing radiative equilibrium situation and the case of a combined mode conduction and radiation heat transfer have been used to test the robustness of the MDOM. Temperature/emissive power and heat flux/energy flow rate distributions in the medium have been found for the effects of various parameters like the extinction coefficient, the scattering albedo, the boundary emissivity and the conduction-radiation parameter. To check the accuracy of the results of the MDOM, results have been compared with those available in the literature and also by obtaining the radiative information using the finite volume method. MDOM has been found to provide accurate results.     

MHD oscillatory flow on free convection–radiation through a porous medium with constant suction velocity

Unsteady two-dimensional hydromagnetic free convection and thermal radiation flow of an electrically conducting viscous-incompressible fluid, through a highly porous medium bounded by a vertical plane surface of constant temperature are presented. The Rosseland diffusion approximation is used to describe the radiative heat flux in the energy equation. Expressions for the velocity and temperature are obtained. The free-stream velocity of the fluid vibrates about a mean constant value and the surface absorbs the fluid with constant velocity. Effects of varying R (radiative parameter), G (Grashof number), k′ (permeability of the porous medium) and M (magnetic parameter upon the velocity field and the effect of varying R and Pr (Prandtl number) on the temperature are discussed.     

Temperature field of radiative equilibrium in a semitransparent slab with a linear refractive index and gray walls

The temperature field in a semitransparent slab of absorbing-emitting gray medium at radiative equilibrium is solved in this paper. The medium has a linear refractive index and the two boundaries are diffuse gray walls. A curved ray tracing technique is combined with a pseudo-source adding method to deduce the radiative intensities on the gray walls. And on the basis of the previous work done by Ben Abdallah and Le Dez, the discrete temperature field in the slab is deduced. The influences of refractive index distribution, boundary wall emissivities and optical thickness on the radiative equilibrium temperature field are examined. The results display the significant influences of the refractive index distribution and the boundary wall emissivities.     

Application de la méthode des ordonnées discrètes au transfert radiatif dans un milieu bidimensionnel gris á géométrie complexe

Application of the discrete method to the radiative heat transfer in a two-dimensional grey medium of complex geometry. This paper describes a new approach in determining the radiative intensity and the temperature fields in a semi-transparent medium enclosed in a two-dimensional cavity the boundary surfaces of which are uniformly grey and purely isotropic diffuse reflectors, with the help of a new combination of ray tracing, finite volumes and discrete ordinates method. Since the grid used can be unstructured, the technique is applicable to the calculation of radiative transfer in enclosures of complex geometry. The basic equations are given, followed by results for cases of simple geometry compared with the exact solutions and the treatment of other cases of more complex geometry. The method eliminates oscillations in the intensity field and yields accurate results.     

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.     

Pairs of opposite decaying evanescent waves in random media and tunneling radiative transfer

A tunneling mechanism of radiative transfer through a dielectric random medium is revealed applying technique of Dyson and Bethe-Salpeter equations for electromagnetic wave multiple scattering by medium inhomogeneities (scatterers) with near fields effects in scattered fields. The mechanism consists in existing inside of a random inhomogeneity a pair of virtually opposite decaying evanescent waves whose interference results in energy flux.     

MOL solution of DOM for transient radiative transfer in 3-D scattering media

A methodology based on the method of lines solution of discrete ordinates method for solution of the 3-D transient radiative transfer equation is introduced. The method is applied to the prediction of transient and steady state transmittances in a cubical enclosure containing purely scattering medium and validated against Monte Carlo solutions from the literature. The flexibility of the method for implementation of linear spatial differencing schemes, flux limiters and weighted essentially non-oscillatory methods is demonstrated. Van Leer flux limiter is found to provide stable, accurate and efficient solutions.     

A modification of the cumulative wavenumber method to compute the radiative heat flux in non-uniform media

This paper presents a modification of the cumulative wavenumber (CW) method to determine the radiative heat flux in non-uniform participating gases. Previous works in the literature have shown that the CW method renders accurate estimates of the radiative volumetric heat source in the medium. However, as will be shown in this work, the radiative heat flux can present considerable deviation of the correct solution, which results from the radiative energy balance not being satisfied by the CW method. A modification of the method is devised in this work to satisfy the radiative energy balance while keeping the same value of the radiative volumetric heat source. The proposed methodology is applied together with the discrete ordinates method to solve the radiation heat transfer in a one-dimensional slab containing a non-isothermal layer of carbon dioxide. The results are compared to the benchmark line-by-line (LBL) integration, and show that the modified CW method can satisfy the radiative energy balance, improving the estimation of the radiative heat flux in the medium.     

Expansions for two-dimensional radiative transfer in an isotropically scattering semi-infinite medium

Small spatial frequency expansions for the source function and radiative flux are obtained for a purely scattering, semi-infinite, two-dimensional medium. Both collimated and diffuse boundary conditions are analyzed. With these expansions, other expansions are obtained which are valid at large optical distances away from the incident radiation. Expansions are presented for a finite strip, circular disk and a Gaussian distribution.     

Transient conductive and radiative transfer in a planar layer with Arrhenius heat generation

Transient conductive and radiative energy transfer in a gray absorbing-emitting planar medium bounded by black walls is studied. The temperature distribution is uniform when the medium starts releasing energy according to the Arrhenius equation. The Crank-Nicolson method is used in the radiative-dominant case with good accuracy. Time-dependent temperature and heat flux distributions are calculated until steady-state is reached. Confirmation of the solution technique is made by comparison with previous investigations. The dimensionless activation energy required for ignition varies from 4.1 for pure conduction to 7.5 for pure radiation.