Analysis of Conduction and Radiation Heat Transfer in a Differentially Heated 2‐D Square Enclosure

Radiative heat transfer with and without conduction in a differentially heated 2‐D square enclosure is analyzed. The enclosure with diffuse gray boundaries contains radiating and/or conducting gray homogeneous medium. Radiatively, the medium is absorbing, emitting and scattering. On the south boundary, four types of discrete heated regions, viz., the full boundary, the left one‐third, left two third and middle one third, are considered. In the absence of conduction, distributions of heat flux along the south boundary are studied for the effect of extinction coefficient. In the presence of conduction, distributions of radiation, conduction and total heat fluxes along the south boundary are analyzed for the effects of extinction coefficient, scattering albedo, conduction–radiation parameter, and south boundary emissivity. Effects of these parameters on centerline temperature distribution are also studied. To assess the performance of three commonly used radiative transfer methods, in all cases, the radiative transfer equation is solved using the discrete ordinate method (DOM), the conventional discrete ordinate method (CDOM) and the finite volume method (FVM). In the combined mode problem, with volumetric radiative information known from one of the three methods, viz., DOM, CDOM, and FVM, the energy equation is solved using the finite difference method (FDM). In all cases, the results from FDM‐DOM, FDM‐CDOM, and FDM‐FVM are in good agreement. Computationally, all three sets of methods are equally efficient.     

Development and comparison of the DTM,the DOM and the FVM formulations for the short-pulse laser transport through a participating medium

The present article deals with the analysis of transient radiative transfer caused by a short-pulse laser irradiation on a participating medium. A general formulation of the governing transient radiative transfer equation applicable to a 3-D Cartesian enclosure has been presented. To solve the transient radiative transfer equation, formulations have been presented for the three commonly used methods in the study of radiative heat transfer, viz., the discrete transfer method, the discrete ordinate method and the finite volume method. To show the uniformity in the formulations in the three methods, the intensity directions and the angular quadrature schemes for computing the incident radiation and heat flux have been taken the same. To validate the formulations and to compare the performance of the three methods, effect of a square short-pulse laser having pulse-width of the order of a femtosecond on transmittance and reflectance signals in case of an absorbing and scattering planar layer has been studied. Effects of the medium properties such as the extinction coefficient, the scattering albedo and the anisotropy factor and the laser properties such as the pulse-width and the angle of incidence on the transmittance and the reflectance signals have been compared. In all the cases, results of the three methods were found to compare very well with each other. Computationally, the discrete ordinate method was found to be the most efficient.     

Chebyshev collocation spectral domain decomposition method for coupled conductive and radiative heat transfer in a 3D L-shaped enclosure

ABSTRACT

This paper presents a Chebyshev collocation spectral domain decomposition method (CSDDM) to study the coupled conductive and radiative heat transfer in a 3D L-shaped enclosure. The partitioned 3D L-shaped enclosure is subdivided into rectangular subdomains based on the concept of domain decomposition. The radiative transfer equation is angularly discretized by the discrete ordinate method with the SRAP_N quadrature scheme and then solved by the CSDDM using the same grid system as in solving the conduction. The effects of the conduction–radiation parameter, the optical thickness, the scattering albedo, and the aspect ratio on thermal behavior of the system are investigated. The results indicate that the 3D CSDDM has a good accuracy and can be considered as a good alternative approach for the solution of the coupled conduction and radiation problems in 3D partitioned domains.     

Combined conductive and radiative heat transfer in an anisotropic scattering participating medium with irregular geometries

This paper deals with the numerical solution for the steady state combined conductive–radiative heat transfer in an anisotropic participating medium within the irregular geometries using the blocked-off method in Cartesian coordinates. The walls of the enclosures were considered to be opaque, diffuse and gray having specified heat flux and temperature boundary conditions. The finite-volume method has been adopted to solve the energy equation and the discrete ordinates method has been employed to solve the radiative transfer equation. The radiative and radiative–conducive models were validated by comparison with the results of specific test cases taken from the literature. The results showed very satisfactory predictions compared with the benchmarked results. As the degree of enclosure complexity (with curved or skewed walls) increased, finer grids were required. Based on this method, the effects of various influencing parameters such as the conduction–radiation parameter, scattering albedo and extinction coefficient have been considered.     

Lattice Boltzmann Method and Modified Discrete Ordinate Method Applied to Radiative Transport in a Spherical Medium with and without Conduction

This article deals with the application of the modified discrete ordinate method (MDOM) to calculate volumetric radiative information with and without conduction in a concentric spherical enclosure containing a participating medium. With radiative information known from the MDOM, the energy equation of the combined mode transient conduction and radiation heat transfer is formulated and solved using the lattice Boltzmann method (LBM). Without conduction, for pure radiation case, two benchmark problems, representing nonradiative and radiative equilibrium situations are taken up. In the case of non-radiative equilibrium, an isothermal medium is bounded by cold walls and medium is the source of radiation, while in the case of radiative equilibrium, nonisothermal medium is confined between a hot and a cold wall, and the hot (inner sphere) wall is the radiation source. Depending upon the problem, heat flux, energy flow rate, emissive power, and temperature distributions in the medium are calculated for different values of parameters such as the extinction coefficient, the scattering albedo, the conduction-radiation parameter, the boundary emissivity, and the radius ratio. To validate the MDOM and the LBM-MDOM formulations, problems are also solved using the finite volume method (FVM) and the finite-difference method (FDM)–FVM approach, in which the FVM is used to calculate the volumetric radiation and the energy equation is also solved using the FDM. Results of the MDOM, LBM–MDOM, FVM and FDM–FVM are also benchmarked against those available in the literature. MDOM and LBM–MDOM have been found to provide accurate results.     

Analysis of Conduction and Radiation Heat Transfer in a 2-D Cylindrical Medium Using the Modified Discrete Ordinate Method and the Lattice Boltzmann Method

This article deals with the analysis of radiative transport with and without conduction in a finite concentric cylindrical enclosure containing absorbing, emitting, and scattering medium. Isothermal medium as the radiation source confined between the cold cylinders and a nonisothermal medium with the inner cylinder as the radiation source are the two nonradiative and radiative equilibrium problems. They involve only calculation of radiative information. In the third problem, a conducting-radiating medium is thermally perturbed by raising the temperature of the inner cylinder. In all problems, radiative information is computed using the modified discrete ordinate method (MDOM), and in the third problem, the lattice Boltzmann method (LBM) is used to formulate and solve the energy equation. Depending on the problems, effects of various parameters such as the extinction coefficient, the scattering albedo, the boundary emissivity, the conduction-radiation parameter, and the radius ratio are studied on temperature and heat flux distributions. The MDOM and the LBM-MDOM results are compared with those available in the literature. To further establish the accuracy of the MDOM and the LBM-MDOM results, in all problems, comparisons are made with the results obtained from the finite volume method (FVM) and the finite difference method-FVM approach, in which FVM provides the radiative information. The selection of the FDM-FVM for the third problem is also with the objective that for this problem, not much work is reported in which the FVM is used to calculate the radiative information. MDOM and LBM-MDOM results are found to compare well with those available in the literature, and in all cases they are in excellent agreement with FVM and FDM-FVM approaches.     

Interaction of Participating Medium Radiation with Thermosolutal Convection—A Critical Appraisal

The present study addresses the interaction effect of participating media radiation with the onset of double‐diffusive convection in a square enclosure. Vertical walls are imposed with constant temperature and concentration, and the horizontal walls are impermeable and adiabatic. The boundaries of the enclosure are diffuse‐gray, and the enclosed fluid isotropically scatters, emits, and absorbs thermal radiation. Numerical simulations have been performed for both aiding and opposing buoyancy conditions. The buoyancy ratio has been varied to simulate the effect of buoyancy driven flow and compositionally driven flow, along with transition of flow between the above. Optical properties like opacity of medium, scattering albedo, Planck number, and wall emissivity have been varied to depict their influence on flow and heat transfer. The modified MAC method is used for the solution of convective transport equations. Gradient dependent consistent hybrid upwind scheme of second order (GDCHUSSO) is used for the discretization of the convective terms. The discrete ordinate method, with S8 approximation, is used to solve the radiative transport equation. The parametric results are provided in graphical and tabular form. Flow lines, isotherms, and isoconcentration contour maps are provided to bring clarity in the understanding of the momentum, heat, and solute transport phenomenon. The stabilization effect of thermal radiation at critical buoyancy ratio for buoyancy opposed flow is observed. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21108     

SLA (Second-law analysis) of transient radiative transfer processes

This paper concerns a SLA (second-law analysis) of transient radiative heat transfer in an absorbing, emitting and scattering medium. Based on Planck’s definition of radiative entropy, transient radiative entropy transfer equation and local radiative entropy generation in semitransparent media with uniform refractive index are derived. Transient radiative exergy transfer equation and local radiative exergy destruction are also derived based on Candau’s definition of radiative exergy. The analytical results are consistent with the Gouy–Stodola theorem of classical thermodynamics. As an application concerning transient radiative transfer, exergy destruction of diffuse pulse radiation in a semitransparent slab is studied. The transient radiative transfer equation is solved using the discontinuous finite element based discrete ordinates equation. Transient radiative exergy destruction is calculated by a post-processing procedure.     

RADIATIVE HEAT TRANSFER IN FINITE CYLINDRICAL HOMOGENEOUS AND NONHOMOGENEOUS SCATTERING MEDIA EXPOSED TO COLLIMATED RADIATION

Radiative heat transfer is studied in a finite axisymmetrical cylindrical enclosure exposed to collimated radiation. The integral equations for radiative transfer are solved by the YIX method and the quadrature method for comparison. Integrated intensity and radiative heat flux are presented in homogeneous and nonhomogeneous scattering media exposed to both uniform and Gaussian distributions of normal collimated incident radiation. The effects of aspect ratio, different incident radiation, and anisotropic scattering phase function as well as nonhomogeneous property distribution are discussed. Ray effects appear in the YIX solution for the case of a nonhomogeneous step change in the extinction coefficient. In order to eliminate the ray effect, an adaptive angular quadrature scheme is described and applied.     

Numerical simulation of coupled radiation‐convection heat transfer of high‐temperature developing flow in tube

By combining the discrete ordinate method with the control volume method, the coupled radiation‐convection heat transfer of high‐temperature developing laminar flow in a tube is investigated numerically. The radiative transfer is solved by the discrete ordinate method and its contribution to thermal balance is dealt with as a source term in the energy equation, which is solved, as well as the momentum equation, by the control volume method. The effects of medium optical thickness and tube wall temperature on the temperature distribution in medium as well as the heat flux and local Nusselt number on wall are analyzed. The results show that the radiation heat transfer of high‐temperature medium influences the temperature distribution and convection heat transfer greatly, and plays an important role in the heat transfer of developing laminar flow in a tube. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 33(1): 53–63, 2004; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10135     

Coupled radiation‐convection heat transfer of high‐temperature participating medium in heated/cooled tubes

The coupled radiation‐convection heat transfer of high‐temperature participating medium in heated/cooled tubes is investigated numerically. The medium flows in a laminar and fully developed state with a Poiseuille velocity distribution, but the thermal status is developing. By the discrete ordinate method, the nonlinear integrodifferential radiative transfer equation in a cylindrical coordinate form is solved to give the radiative source term in the energy equation of coupled heat transfer. The energy equation is solved by the control volume method. The local Nusselt number and wall heat flux of convection as well as the total wall heat flux are employed to evaluate the influence of radiation heat transfer on convection. The analysis shows that the radiation heat transfer weakens the convection effect, promotes the temperature development, and significantly shortens the tube length with obvious heated/cooled effect. There is an obvious difference between the coupled heat transfer in a heated tube and that in a cooled tube, even though the medium properties are kept constant. The wall emissivity, the medium thermal conductivity and scattering albedo have significant influences on the coupled heat transfer, but the effect of medium scattering phase function is small. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 33(1): 64–72, 2004; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10137     

RADIATION ELEMENT METHOD FOR TRANSIENT HYPERBOLIC RADIATIVE TRANSFER IN PLANE-PARALLEL INHOMOGENEOUS MEDIA

In this study the radiation element method is formulated to solve transient radiative transfer with light radiation propagation effect in scattering, absorbing, and emitting media with inhomogeneous property. The accuracy of the method is verified by good agreement between the present calculations and Monte Carlo simulations. The sensitivity of the method against element size, ray emission number, and time increment size is examined. The transient effect of radiation propagation is essential in short-pulse laser radiation transport when the input pulse width is not considerably larger than the system radiation propagation time. The transient characteristics of radiative transfer are investigated in the media subject to collimated laser irradiation and/or diffuse irradiation withtemporal Gaussian and/or square profiles. The inhomogeneous profile of extinction coefficient of the medium affects strongly the transient radiative flux divergence inside the medium.     

Discrete ordinate solutions associated with the finite Legendre transform for radiative transfer in a slab with sinusoidal refractive index

In this work, we develop an alternative discrete ordinate approximation for radiative transfer in a refractive slab. The present method treats the angular derivative term of the radiative transfer equation for a planar medium with varying refractive index (VRI) by using a finite Legendre transform which gives a simple expression of the angular derivative term. Thus, the solution procedure does not march along direction, and so is not restricted to a monotonic variation of refractive index. We apply this method to study radiative heat transfer in a cold slab with anisotropic scattering, diffuse boundaries and sinusoidal VRI. We also solve the problems by the discrete curved ray tracing (DCRT). The hemispherical reflectance and transmittance of slabs with irradiation from the upper surroundings obtained by the present method and those obtained by the DCRT are in excellent agreement. For a slab of a sinusoidal refractive index with the minimum at the center plane, the gradient of refractive index causes the internal reflection of a part of irradiation, which reduces the transmittance of the slab. Other effects of the VRI, the optical thickness, the scattering albedo, the anisotropically scattering coefficient and the boundary reflection are also investigated.     

Multilayer differential discrete ordinate method for inhomogeneous participating media

This paper presents a multilayer differential discrete ordinate method to solve the radiative transfer equation for an absorbing, emitting and scattering inhomogeneous plane parallel medium. This method reduces the integro-differential equation into a set of coupled first order ordinary differential equations with two point boundary conditions on using a suitable quadrature scheme. These equations are then solved numerically. Numerical validation of the method for gray medium is done by comparing the results obtained with benchmark cases available in the literature. Validation for a non-gray medium is done by considering a problem concerning radiative transfer from the atmosphere. The brightness temperature at the top of the atmosphere is calculated at various frequencies and validated with those obtained by several other numerical methods.     

The finite volume method approach to the collapsed dimension method in analyzing steady/transient radiative transfer problems in participating media

With the finite volume formulation (FVM) approach applied to the collapsed dimension method (CDM), this article deals with the application of the CDM to analyze radiative heat transfer problems in a participating medium subjected to a continuous diffuse or a continuous/short-pulse collimated boundary radiative loading. The planar medium contained between diffuse gray boundaries is absorbing, emitting and anisotropically scattering. With three categories of thermal boundary radiative loadings, for the four types of problems considered, the CDM results are compared for a wide range of radiative parameters with that of the FVM.     

Application of the lattice Boltzmann method for solving the energy equation of a 2-D transient conduction-radiation problem

A lattice Boltzmann method (LBM) is used to solve the energy equation in a test problem involving thermal radiation and to thus investigate the suitability of scalar diffusion LBM for a new class of problems. The problem chosen is transient conductive and radiative heat transfer in a 2-D rectangular enclosure filled with an optically absorbing, emitting and scattering medium. The energy equation of the problem is solved alternatively with a previously used finite volume method (FVM) and with the LBM, while the radiative transfer equation is solved in both cases using the collapsed dimension method. In a parametric study on the effects of the conduction-radiation parameter, extinction coefficient, scattering albedo, and enclosure aspect ratio, FVM and LBM are compared in each case. It is found that, for given level of accuracy, LBM converges in fewer iterations to the steady-state solution, independent of the influence of radiation. On the other hand, the computational cost per iteration is higher for LBM than for the FVM for a simple grid. For coupled radiation-diffusion, the LBM is faster than the FVM because the radiative transfer computation is more time-consuming than that of diffusion.     

Combined conduction and radiation heat transfer in a gray anisotropically scattering planar medium with diffuse-specular boundaries

Coupled conduction and radiation heat transfer in a gray planar nonlinearly anisotropic scattering medium bounded between two plane parallel surfaces reflecting both diffusely and specularly is analyzed. The governing integrodifferential equations are solved by a numerical iterative method consisting of Numerov's method to solve the energy equation and Chandarsekhar's discrete ordinates method in conjunction with the Crank-Nicolson method to solve the radiative transfer equation. Convergence of the solution is enhanced by Ng-acceleration. The numerical algorithm described is found to be fast and reliable. Numerical results based on S₃₂ method indicate that anisotropy plays an important role, and difference between the diffuse and specular reflections is found to be insignificant.     

A Conduction–Radiation Mixture Model for Laser-Assisted Phase Change of Semitransparent Material

A one-dimensional transient coupled conduction-radiation numerical model is developed to investigate the laser melting of semitransparent material under a continuous collimated laser pulse in a convective cooling environment. The medium is considered absorbing, emitting, and scattering. The thermophysical properties are taken to be different for different phase fields. Volumetric radiation is incorporated in the proposed model. The radiation information is obtained by solving the equation of transfer. The temperature field is obtained by solving the energy equation with internal radiation source. The finite-volume method is used to discretize both the equation of transfer and the energy equation. The enthalpy formulation is adopted to capture the continuously evolving solid–liquid interface during the phase change. The laser source is approximated with the collimated radiation source. Collimated intensity is captured directly (without splitting the total intensity into two parts: diffuse and collimated) by adjusting the control angles. The present model is first validated with the existing phase-change model in the literature. Then the effects of different parameters such as optical thickness, scattering albedo, and the conduction–radiation parameter on the liquid fractions and temperature distribution in the medium are studied. It is observed that when the radiation is dominant, the temperature in the medium is high and hence the liquid fraction is more, in contrast to conduction-dominated phase change.     

NATURAL CONVECTION AND RADIATION IN A SQUARE ENCLOSURE

Natural convection of a radiating fluid in a square enclosure is studied numerically. The coupled momentum, energy, and radiative transfer equations are solved by an iterative procedure. The solutions to the equation of radiative transfer are obtained by the discrete ordinates method using S4 and S8 quadratures. The method is based on control volume formulation and is fully compatible with the SIMPLER algorithm used to solve the momentum and energy equations. The effects of optical thickness and scattering on the flow and temperature fields and heat transfer rates are analyzed. The changes in the buoyant flow patterns and temperature distributions due to the presence of radiation in inclined or heat generating enclosures are also studied. Comparative results obtained by the P-I differential approximation are presented.     

Simultaneous estimation of extinction coefficient distribution, scattering albedo and phase function of a two-dimensional medium

In this work, we develop a scheme based on solving the general integro-differential equation of radiation transport to estimate simultaneously the distribution of the extinction coefficient, the scattering albedo and the phase function of a two-dimensional inhomogeneneous medium with less diffusive radiation. The forward problem for a cylindrical medium subjected to collimated incident radiation is solved by the discrete-ordinate method. The inverse radiation problem is formulated as a least square problem that minimizes the discrepancy between the measured and the calculated leaving radiative fluxes. The Levenberg-Marquardt algorithm is applied to the least square problems for a variety of cases. The results obtained show that this scheme can reconstruct accurate enough results for most of the cases considered. Comparisons of the results show that the accuracy of the estimated results decreases with the increase of the scattering albedo and we need more discrete ordinates to generate accurate enough estimated results for an optically thin case. The estimated results obtained from the measurement data with moderate errors are still acceptable.