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.