Least-squares finite element method for radiation heat transfer in graded index medium

To avoid the complicated and time-consuming computation of curved ray trajectories, a least-squares finite element method based on discrete ordinate equation is extended to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Four cases of radiative heat transfer are examined to verify this least-squares finite element method. Linear and nonlinear graded index are considered. The predicted dimensionless net radiative heat fluxes are determined by the least-squares finite element method and compared with the results obtained by other methods. The results show that the least-squares finite element method is stable and has a good accuracy in solving the multi-dimensional radiative transfer problem in a semitransparent graded index medium, while the Galerkin finite element method sometimes suffers from nonphysical oscillations.     

Discontinuous finite element method for radiative heat transfer in semitransparent graded index medium

To avoid the complicated and time-consuming computation of curved ray trajectories, a discontinuous finite element method based on discrete ordinate equation is extended to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Two cases of radiative heat transfer in two-dimensional rectangular gray semitransparent graded index medium enclosed by opaque boundary are examined to verify this discontinuous finite element method. Special layered and radial graded index distributions are considered. The predicted dimensionless net radiative heat fluxes and dimensionless temperature distributions are determined by the discontinuous finite element method and compared with the results obtained by the curved Monte Carlo method in references. The results show that the discontinuous finite element method has a good accuracy in solving the multi-dimensional radiative transfer problem in a semitransparent graded index medium.     

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.     

Finite element approach for radiative transfer in multi-layer graded index cylindrical medium with Fresnel surfaces

Because the optical plane defined by the incidence and reflection direction at a cylindrical surface has a complicated relation with the local azimuthal angle and zenith angle in the traditional cylindrical coordinate system, it is difficult to deal with the specular reflective boundary condition in the solution of the traditional radiative transfer equation for cylindrical system. In this paper, a new radiative transfer equation for graded index medium in cylindrical system (RTEGCN) is derived based on a newly defined cylindrical coordinate system. In this new cylindrical coordinate system, the optical plane defined by the incidence and reflection direction is just the isometric plane of the local azimuthal angle, which facilitates the RTEGCN in dealing with cylindrical specular reflective boundaries. A least squares finite element method (LSFEM) is developed for solving radiative transfer in single and multi-layer cylindrical medium based on the discrete ordinates form of the RTEGCN. For multi-layer cylindrical medium, a radial basis function interpolation method is proposed to couple the radiative intensity at the interface between two adjacent layers. Various radiative transfer problems in both single and multi-layer cylindrical medium are tested. The results show that the present finite element approach has good accuracy to predict the radiative heat transfer in multi-layer cylindrical medium with Fresnel surfaces.     

Benchmark numerical solutions for radiative heat transfer in two-dimensional medium with graded index distribution

In graded index media, the ray goes along a curved path determined by Fermat principle. Generally, the curved ray trajectory in graded index media is a complex implicit function, and the curved ray tracing is very difficult and complex. Only for some special refractive index distributions, the curved ray trajectory can be expressed as a simple explicit function. Two important examples are the layered and the radial graded index distributions. In this paper, the radiative heat transfer problems in two-dimensional square semitransparent with layered and radial graded index distributions are analyzed. After deduction of the ray trajectory, the radiative heat transfer problems are solved by using the Monte Carlo curved ray-tracing method. Some numerical solutions of dimensionless net radiative heat flux and medium temperature are tabulated as the benchmark solutions for the future development of approximation techniques for multi-dimensional radiative heat transfer in graded index media.     

Meshless approach for coupled radiative and conductive heat transfer in one-dimensional graded index medium

A meshless local Petrov-Galerkin (MLPG) approach is employed for solving the coupled radiative and conductive heat transfer in a one-dimensional slab with graded index media. The angular distribution term in discrete ordinate equation of radiative transfer within a one-dimensional graded index slab is discretized by a step scheme, and the meshless approach for radiative transfer is based on the discrete ordinate equation. A moving least-squares approximation is used to construct the shape function. Two particular test cases for coupled radiative and conductive heat transfer within a one-dimensional graded index slab are examined to verify this new approximate method. The temperatures and the radiative heat fluxes are obtained. The results are compared with the other benchmark approximate solutions. By comparison, the results show that the MLPG approach has a good accuracy in solving the coupled radiative and conductive heat transfer in one-dimensional graded index media.     

Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems.     

Verification of numerical simulation method for entropy generation of radiation heat transfer in semitransparent medium

The numerical simulation method of radiative entropy generation in participating media presented by Caldas and Semiao [Entropy generation through radiative transfer in participating media: analysis and numerical computation. JQSRT 2005;96:423-37] is extended to analyze the radiative entropy generation in the enclosures filled with semitransparent media. A discrete ordinates method is used to solve radiative transfer equation and radiative entropy generation. Two different examples are employed to verify the numerical simulation method of radiative entropy generation in the enclosure. Numerical results of dimensionless radiative entropy generation of enclosure are identical to that of entire thermodynamics analysis for the enclosure system. This numerical simulation method can be used in the entropy generation analysis of high-temperature systems such as boilers and furnaces, in which radiation is the dominant mode of heat transfer.     

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.     

Modified finite volume method applied to radiative transfer in 2D complex geometries and graded index media

A modified finite volume method with unstructured triangular meshes is proposed to solve the RTE in 2D complex geometries and for graded index media. In such media, the RTE has an additional term corresponding to “angular redistribution”. This term is due to the change in the orientation of the direction of propagation for the radiation along curved optical paths. Some benchmark cases applied to a slab (1D) and a square cavity (2D) with linear and nonlinear refractive graded index are used to validate the new method. New results are presented for a disk with radial graded index.     

Spectral collocation method with a flexible angular discretization scheme for radiative transfer in multi-layer graded index medium

The spectral collocation method (SCM) is employed to solve the radiative transfer in multi-layer semitransparent medium with graded index. A new flexible angular discretization scheme is employed to discretize the solid angle domain freely to overcome the limit of the number of discrete radiative direction when adopting traditional S_N discrete ordinate scheme. Three radial basis function interpolation approaches, named as multi-quadric (MQ), inverse multi-quadric (IMQ) and inverse quadratic (IQ) interpolation, are employed to couple the radiative intensity at the interface between two adjacent layers and numerical experiments show that MQ interpolation has the highest accuracy and best stability. Variable radiative transfer problems in double-layer semitransparent media with different thermophysical properties are investigated and the influence of these thermophysical properties on the radiative transfer procedure in double-layer semitransparent media is also analyzed. All the simulated results show that the present SCM with the new angular discretization scheme can predict the radiative transfer in multi-layer semitransparent medium with graded index efficiently and accurately.     

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.     

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 formulation of the first- and second-order discrete ordinates equations for radiative heat transfer calculation in three-dimensional participating media

Two finite element methods (FEMs), FEDOM1 and FEDOM2 (standing for the first and the second finite element discrete ordinates methods, respectively), are formulated and numerically tested. The reference second-order discrete equation is modified in its scattering terms and is applied to the problems of absorbing/emitting and anisotropically scattering media by using the FEM. Numerical features of the developed FEMs are compared with one of the discrete ordinates interpolation method (DOIM), which uses a finite difference scheme. Prediction results of radiative heat transfer by these two FEMs are compared with reference solutions and verified in three-dimensional enclosures containing participating media. The results of FEDOM1 and FEDOM2 agree well with exact solutions for the problem of absorbing/emitting medium with various range of optical thickness. Generally, the two FEMs show more accurate results than DOIM. And FEDOM1 shows more accurate results than FEDOM2 in most of the test problems. Both of the developed FEMs show reasonable results compared with published Monte Carlo solutions for the tested absorbing/emitting and anisotropically scattering media. Although the FEDOM2 is not as accurate as the FEDOM1, it shows its own advantages that it reduces CPU time and memory space of dependent variable to half.     

Least-squares collocation meshless approach for radiative heat transfer in absorbing and scattering media

A least-squares collocation meshless method is employed for solving the radiative heat transfer in absorbing, emitting and scattering media. The least-squares collocation meshless method for radiative transfer is based on the discrete ordinates equation. A moving least-squares approximation is applied to construct the trial functions. Except for the collocation points which are used to construct the trial functions, a number of auxiliary points are also adopted to form the total residuals of the problem. The least-squares technique is used to obtain the solution of the problem by minimizing the summation of residuals of all collocation and auxiliary points. Three numerical examples are studied to illustrate the performance of this new solution method. The numerical results are compared with the other benchmark approximate solutions. By comparison, the results show that the least-squares collocation meshless method is efficient, accurate and stable, and can be used for solving the radiative heat transfer in absorbing, emitting and scattering media.     

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.     

Discrete transfer method applied to radiative transfer in a variable refractive index semitransparent medium

Application of the discrete transfer method (DTM) has been extended to the analysis of radiative heat transfer in a variable refractive index participating medium. To validate the DTM formulation, radiative heat transfer in an absorbing, emitting and isotropically scattering planar medium was considered. The participating medium was assumed to be in radiative equilibrium. For both constant and variable refractive indices of the medium, the DTM results were compared with those available in the literature. The DTM was found to provide accurate results.     

Discontinuous spectral element method for solving radiative heat transfer in multidimensional semitransparent media

A discontinuous spectral element method (DSEM) is presented to solve radiative heat transfer in multidimensional semitransparent media. This method is based on the general discontinuous Galerkin formulation. Chebyshev polynomial is used to build basis function on each element and both structured and unstructured elements are considered. The DSEM has properties such as hp-convergence, local conservation and its solutions are allowed to be discontinuous across interelement boundaries. The influences of different schemes for treatment of the interelement numerical flux on the performance of the DSEM are compared. The p-convergence characteristics of the DSEM are studied. Four various test problems are taken as examples to verify the performance of the DSEM, especially the performance to solve the problems with discontinuity in the angular distribution of radiative intensity. The predicted results by the DSEM agree well with the benchmark solutions. Numerical results show that the p-convergence rate of the DSEM follows exponential law, and the DSEM is stable, accurate and effective to solve multidimensional radiative transfer in semitransparent media.     

Finite element method for the two-dimensional atmospheric radiative transfer

The finite element method is applied to the solution of the two-dimensional atmospheric radiative transfer. The analysis is mainly focussed on the derivation of the cell or element equation. The Galerkin method and several hybrid methods using the integral and finite difference form of the radiative transfer equation are employed to obtain the cell equation. The assembled system of equations relating the radiances at the lower and upper boundary of the domain is solved by a direct method.     

Radiative intensity solution and thermal emission analysis of a semitransparent medium layer with a sinusoidal refractive index

The radiative intensity in a sinusoidal refractive index semitransparent medium layer is solved by the curved ray-tracing method in combination with the pseudo-source adding method. One boundary of the medium layer is an opaque diffuse substrate wall. The other boundary is a semitransparent specular or diffuse surface, from which the medium thermal emission emerges. With considering a linear temperature distribution, the radiative intensity formulae are, respectively, deduced under the two boundary conditions. On the basis of the radiative intensity solutions, the directional and hemispherical emission of the medium layer with a specular surface as well as the hemispherical emission of that with a diffuse surface are calculated. The influences of the optical thickness, sinusoidal refractive index distribution and linear temperature distribution on the thermal emission are investigated. The results show that the effects of refractive index and temperature distribution are significant and are different under the two reflecting modes of the surface.