Discrete ordinates interpolation method for radiative heat transfer problems in three-dimensional enclosures filled with non-gray or scattering medium

The discrete ordinates interpolation method (DOIM) is applied to three groups of problems of radiative heat transfer in three-dimensional rectangular enclosures containing non-gray or scattering medium. The original DOIM is first extended to a gray gas model using a new geometric interpolation scheme. It is applied to participating media for different scattering phase functions and optical thicknesses. For the non-gray gas model, the DOIM coupled with the narrow band-based weighted-sum-of-gray-gases (WSGG) model is developed. A few test problems with real gases such as pure H₂O and a mixture of CO₂, H₂O and N₂ are taken. The wall heat flux is calculated and compared with the exact solutions or reference values. All results of test problems are found to be reliable in this study. The DOIM closely reproduces the Monte Carlo reference solutions for different scattering phase functions and optical thicknesses. The non-gray gas results are compared with reference calculations based on the statistical narrow band model and they also show good agreements. The DOIM shows a remarkable merit in the computation time and the grid compatibility, to prove its usefulness for engineering applications.     

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.     

Network simulation method applied to radiation and viscous dissipation effects on MHD unsteady free convection over vertical porous plate

The effects of thermal radiation and viscous dissipation on magneto-hydrodynamic (MHD) unsteady free-convection flow over a semi-infinite vertical porous plate are analysed. The fluid considered is non-gray (absorption coefficient dependent on wave length). The Network Simulation Method is used to solve the boundary-layer equations based on the finite-difference formulation; only discretization of the spatial co-ordinates is necessary, while time remains as a real continuous variable. This method provides a solution for both transient and steady-state problems at the same time, and programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The velocity, temperature, local skin-friction and local Nusselt number are studied for different parameters, including the radiation parameter, Eckert number, magnetic number and suction (or injection).     

Transient radiative and conductive heat transfer in non-gray semitransparent two-dimensional media with mixed boundary conditions

This paper is devoted to transient heat transfer involving radiation and conduction. Considering a non-gray purely absorbing media, the radiative heat transfer equation (RTE) is solved iteratively with the Discrete Ordinates Method (DOM) using an exponential differencing scheme. The energy balance equation is used to compute temperature at each time step with the Crank–Nicholson technique. Energy equation is coupled to the RTE through the radiative source term. Both equations are discretized with finite differencing schemes. The energy conservation leads to the sparse system of linear equations A× T=B which is solved with a bi-conjugate stabilized gradient technique (BCSG). Validation of the model with different test cases is achieved and application to transient heating of glass is also studied.     

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.     

On the combined heat transfer in the multilayer non-gray porous fibrous insulation

A detailed numerical modeling is performed to investigate coupled heat transfer of natural convection, radiation and conduction in high-temperature multilayer thermal insulation (MTI), which consists of high-porous, non-gray semitransparent fibrous materials and reflective foils. Radiation within fibers, radiation between fibers and the reflective foils, conduction within fibers and convection from the fibers to the surrounding fluid are considered. Macroscopic (porous media) modeling is used to determine velocity, pressure and temperatures fields for fibrous insulation with a random packing geometry under natural convection, whereas the radiative transfer equation (RTE) is used to solve the radiative heat flux for non-gray materials. Key features of the macroscopic model include two-dimensional effects, non-gray radiative exchange, and the relaxation of the local thermodynamic non-equilibrium (LTNE). This model was validated by comparison with experimental data and it was used to investigate natural convection of coupled heat transfer in multilayer insulation, numerical results showed that natural convection is more likely to occur when the heated/cooled rate is low, while natural convection can be ignored in simulating steady-state coupled heat transfer in MTI.     

Numerical simulation of radiative heat transfer from non-gray gases in three-dimensional enclosures

The discrete ordinates and the discrete transfer methods are applied to the numerical simulation of radiative heat transfer from non-gray gases in three-dimensional enclosures. Several gas radiative property models are used, namely the correlated k-distribution (CK), the spectral line-based weighted-sum-of-gray-gases (SLW) and the weighted-sum-of-gray-gases (WSGG) methods. The results are compared with recently published accurate calculations based on the statistical narrow band model. The WSGG model is computationally efficient, but often yields relatively large errors. It should be used only if moderate accuracy is sufficient. The SLW model is the best alternative regarding the compromise between accuracy and numerical efficiency. However, an optimization of the coefficients of the model is essential to reduce the computational requirements, especially in the case of gas mixtures. The CK model is the most accurate of the methods evaluated here, but too time consuming for engineering applications, although recent developments may partly overcome this shortcoming.     

A new inverse approach for the equivalent gray radiative property of a non-gray medium using a modified zonal method and the complex-variable-differentiation method

Non-gray radiative properties of an absorbing, emitting, non-gray participating medium significantly increase the difficulty of solving the radiative transfer equation. This paper presents a new inverse approach for the equivalent gray radiative property of a non-gray medium. In this approach, the unknown equivalent gray radiative properties are treated as the optimization variables, and the errors to be minimized are the differences between the calculated temperatures and the measured ones. The measured data are simulated by solving the direct problem, in which a modified zonal method together with the Edwards exponential wide-band model is employed. In the inverse problem, the sensitivity coefficients are first calculated by the complex-variable-differentiation method, and then the least-square method and the Newton-Raphson iterative method are employed to minimize the target function. The effectiveness and efficiency of the inverse problem are demonstrated in an example, and another case is given to show the accuracy and potential of the proposed algorithm. The effects of the measurement error and the number of measurement points on the accuracy of the inverse analysis are also investigated in detail.     

Unsteady coupling of Navier-Stokes and radiative heat transfer solvers applied to an anisothermal multicomponent turbulent channel flow

Direct numerical simulations (DNS) of an anisothermal reacting turbulent channel flow with and without radiative source terms have been performed to study the influence of the radiative heat transfer on the optically non-homogeneous boundary layer structure. A methodology for the study of the emitting/absorbing turbulent boundary layer (TBL) is presented. Details on the coupling strategy and the parallelization techniques are exposed. An analysis of the first order statistics is then carried out. It is shown that, in the studied configuration, the global structure of the thermal boundary layer is not significantly modified by radiation. However, the radiative transfer mechanism is not negligible and contributes to the heat losses at the walls. The classical law-of-the-wall for temperature can thus be improved for RANS/LES simulations taking into account the radiative contribution.