Nonaxisymmetric radiative transfer in inhomogeneous cylindrical media with anisotropic scattering

In this paper, the control volume finite element method (CVFEM) is applied for the first time to solve nonaxisymmetric radiative transfer in inhomogeneous, emitting, absorbing and anisotropic scattering cylindrical media. Mathematical formulations as well as numerical implementation are given and the final discretized equations are based on similar meshes used for convective and conductive heat transfer in computational fluid dynamic analysis. In order to test the efficiency of the developed method, four nonaxisymmetric problems have been examined. Also, the grid dependence and the false scattering of the CVFEM are investigated and compared with the finite volume method and the discrete ordinates interpolation method.     

Three-dimensional radiative transfer modeling using the control volume finite element method

In the present study, a three-dimensional algorithm for the treatment of radiative heat transfer in emitting, absorbing and scattering media is developed. The approach is based on the utilization of control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied at the first time to 3D radiative heat transfer in participating media. The accuracy of the present algorithm is tested by comparing its predictions to other published works. Comparisons show that CVFEM produces good results. Moreover, this approach permits compatibility with other numerical methods used for computational fluids mechanics problems.     

Parallel Jacobian-free Newton Krylov solution of the discrete ordinates method with flux limiters for 3D radiative transfer

The present study introduces a parallel Jacobian-free Newton Krylov (JFNK) general minimal residual (GMRES) solution for the discretized radiative transfer equation (RTE) in 3D, absorbing, emitting and scattering media. For the angular and spatial discretization of the RTE, the discrete ordinates method (DOM) and the finite volume method (FVM) including flux limiters are employed, respectively. Instead of forming and storing a large Jacobian matrix, JFNK methods allow for large memory savings as the required Jacobian-vector products are rather approximated by semiexact and numerical formulations, for which convergence and computational times are presented. Parallelization of the GMRES solution is introduced in a combined memory-shared/memory-distributed formulation that takes advantage of the fact that only large vector arrays remain in the JFNK process. Results are presented for 3D test cases including a simple homogeneous, isotropic medium and a more complex non-homogeneous, non-isothermal, absorbing–emitting and anisotropic scattering medium with collimated intensities. Additionally, convergence and stability of Gram–Schmidt and Householder orthogonalizations for the Arnoldi process in the parallel GMRES algorithms are discussed and analyzed. Overall, the introduction of JFNK methods results in a parallel, yet scalable to the tested 2048 processors, and memory affordable solution to 3D radiative transfer problems without compromising the accuracy and convergence of a Newton-like solution.     

Numerical analyses of radiative heat transfer in any arbitrarily-shaped axisymmetric enclosures

A numerical approach for the treatment of radiative heat transfer in any irregularly-shaped axisymmetric enclosure filled with absorbing, emitting and scattering gray media is developed. Radiative transfer equation (RTE) is formulated for a general axisymmetric geometrical configurations, and the discretized equation is conducted using an unstructured meshes, generated by an appropriate computer algorithm, and the control volume finite element method which frequently adopted in CFD problems. A computer procedure has been done to solve the discretized RTE and to examine the accuracy and the computational efficiency of the proposed numerical approach. By using this computer algorithm, five test cases, a cylindrical enclosure with absorbing and emitting medium, a diffuser shaped axisymmetric enclosure, a finite axisymmetric cylindrical enclosure with a curved wall, a furnace with axially varying medium temperature and a rocket nozzle, are treated and the obtained results agree very well with other published works. Furthermore, the developed computer procedure has an accurate CPU time and it can be coupled easily with CFD codes.     

Prediction of radiative heat transfer in 2D irregular geometries using the collocation spectral method based on body-fitted coordinates

Recently, an efficient numerical method, which is called the collocation spectral method (CSM), for radiative heat transfer problems, has been proposed by the present authors. In this numerical method there exists the exponential convergence rate, which can obtain a very high accuracy even using a small number of grids. In this article, the CSM based on body-fitted coordinates (BFC) is extended to simulate radiative heat transfer problems in participating medium confined in 2D complex geometries. This numerical method makes simultaneously the use of the merits of both the CSM and BFC. In this numerical approach, the radiative transfer equation (RTE) in orthogonal Cartesian coordinates should be transformed into the equation in body-fitted nonorthogonal curvilinear coordinates. In order to test the efficiency of the developed method, several 2D complex irregular enclosures with curved boundaries and containing an absorbing, emitting and scattering medium are examined. The results obtained by the CSM are assessed by comparing the predictions with those in references. These comparisons indicate that the CSM based on BFC can be recommended as a good option to solve radiative heat transfer problems in complex geometries.     

Prediction of radiative heat transfer in 3D complex geometries using the unstructured control volume finite element method

In this paper, a 3D algorithm for the treatment of radiative heat transfer in emitting, absorbing, and scattering media is developed. The numerical approach is based on the utilization of the unstructured control volume finite element method (CVFEM) which, to the knowledge of the authors, is applied for the first time to simulate radiative heat transfer in participated media confined in 3D complex geometries. This simulation makes simultaneously the use of the merits of both the finite element method and the control volume method. Unstructured 3D triangular element grids are employed in the spatial discretization and azimuthal discretization strategy is employed in the angular discretization. The general discretization equation is presented and solved by the conditioned conjugate gradient squared method (CCGS). In order to test the efficiency of the developed method, several 3D complex geometries including a hexahedral enclosure, a 3D equilateral triangular enclosure, a 3D L-shaped enclosure and 3D elliptical enclosure are examined. The results are compared with the exact solutions or published references and the accuracy obtained in each case is shown to be highly satisfactory. Moreover, this approach required a less CPU time and iterations compared with those of even parity formulation of the discrete ordinates method.     

Application of the lattice Boltzmann method to transient conduction and radiation heat transfer in cylindrical media

In this paper, the lattice Boltzmann method (LBM) is applied to solve the energy equation of a transient conduction-radiation heat transfer problem in a two-dimensional cylindrical enclosure filled with an emitting, absorbing and scattering media. The control volume finite element method (CVFEM) is used to obtain the radiative information. To demonstrate the workability of the LBM in conjunction with the CVFEM to conduction-radiation problems in cylindrical media, the energy equation of the same problem is also solved using the finite difference method (FDM). The effects of different parameters, such as the grid size, the scattering albedo, the extinction coefficient and the conduction-radiation parameter on temperature distribution within the medium are studied. Results of the present work are compared with those available in the literature. LBM-CVFEM results are also compared with those given by the FDM-CVFEM. In all cases, good agreement has been obtained.     

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.     

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.     

二维稳态辐射传输方程的有限差分求解法

针对扩散光学层析在小动物成像中的应用问题并基于混浊介质空间光子三维散射的实际物理效应,提出的二维稳态辐射传输方程的有限差分数值求解新方法.在此基础上,研究了不同的空间剖分网格和角度离散密度对求解准确度的影响,并通过将所提方法与蒙特卡洛模拟进行比对,验证方法的正确性.研究表明：在均匀组织体内,当离散角度达到一定数量时,由辐射传输方程的有限差分解获得的透射面和侧面的光子密度对空间网格大小并不敏感,而在反射面上光子密度计算则需要较密的空间网格才能够达到一定准确度.本研究为发展基于辐射传输方程的扩散光学层析理论奠定了基础.     

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.     

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.     

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.     

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.     

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.     

非规则形状介质内辐射-导热耦合传热的间断有限元求解

采用间断有限元法(discontinuous finite element method,DFEM)求解非规则形状介质内的辐射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热问题的高精度数值结果.和传统连续型有限元方法不同,DFEM将计算区域划分成相互独立的离散单元,形函数的构造、未知量的加权近似以及控制方程的求解均在每一个离散单元上进行.通过在单元之间施加迎风格式的数值通量,DFEM保证了整个计算区域的连续性,因此这种方法兼具良好的几何灵活性和局部守恒性.推导了辐射传输方程和能量扩散方程的射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热的高精度数值结果.     

Evaluation method for radiative heat transfer in polydisperse water droplets

Simplifications of the model for nongray radiative heat transfer analysis in participating media comprised of polydisperse water droplets are presented. Databases of the radiative properties for a water droplet over a wide range of wavelengths and diameters are constructed using rigorous Mie theory. The accuracy of the radiative properties obtained from the database interpolation is validated by comparing them with those obtained from the Mie calculations. The radiative properties of polydisperse water droplets are compared with those of monodisperse water droplets with equivalent mean diameters. Nongray radiative heat transfer in the anisotropic scattering fog layer, including direct and diffuse solar irradiations and infrared sky flux, is analyzed using REM². The radiative heat fluxes within the fog layer containing polydisperse water droplets are compared with those in the layer containing monodisperse water droplets. Through numerical simulation of the radiative heat transfer, polydisperse water droplets can be approximated by using the Sauter diameter, a technique that can be useful in several research fields, such as engineering and atmospheric science. Although this approximation is valid in the case of pure radiative transfer problems, the Sauter diameter is reconfirmed to be the appropriate diameter for approximating problems in radiative heat transfer, although volume-length mean diameter shows better accordance in some cases. The CPU time for nongray radiative heat transfer analysis with a fog model is evaluated. It is proved that the CPU time is decreased by using the databases and the approximation method for polydisperse particulate media.     

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.     

Finite element method for modeling radiative transfer in semitransparent graded index cylindrical medium

Both Galerkin finite element method (GFEM) and least squares finite element method (LSFEM) are developed and their performances are compared for solving the radiative transfer equation of graded index medium in cylindrical coordinate system (RTEGC). The angular redistribution term of the RTEGC is discretized by finite difference approach and after angular discretization the RTEGC is formulated into a discrete-ordinates form, which is then discretized based on Galerkin or least squares finite element approach. To overcome the RTEGC-led numerical singularity at the origin of cylindrical coordinate system, a pole condition is proposed as a special mathematical boundary condition. Compared with the GFEM, the LSFEM has very good numerical properties and can effectively mitigate the nonphysical oscillation appeared in the GFEM solutions. Various problems of both axisymmetry and nonaxisymmetry, and with medium of uniform refractive index distribution or graded refractive index distribution are tested. The results show that both the finite element approaches have good accuracy to predict the radiative heat transfer in semitransparent graded index cylindrical medium, while the LSFEM has better numerical stability.     

Non-gray radiative transfer in plane-parallel media with reflecting boundaries

A generalized equation of radiative transfer in the two-group picket-fence model is analyzed for a plane parallel, emitting, absorbing and isotropically scattering medium containing uniform heat sources and having boundary surfaces which are diffuse emitters and diffuse reflectors and are maintained at uniform but arbitrary temperatures. The solution of the general problem is expressed by the superposition of simpler problems which are solved by the application of the normal-mode-expansion technique. Highly accurate numerical results are presented for the temperature distribution and the radiative heat flux in the medium.