Adaptability analysis of radiative transport diffusion approximation in planar-graded-index media

Transcription

1 Research Article Adaptability analysis of radiative transport diffusion approximation in planar-graded-index media Advances in Mechanical Engineering 2018, ol. 10(11) 1 6 Ó The Author(s) 2018 DOI: / journals.sagepub.com/home/ade Lijun Liu, Hong Shi and Xiaoyan Liu Abstract To study the applicability of diffusion approximation method, the influence of the optical thickness, scattering albedo, and wall emissivity of the examined medium on the accuracy of the diffusion approximation is analyzed. The accuracy is determined by comparison with results obtained by using the Monte Carlo method. By calculating the dimensionless radiation heat flux distributions at the bottom of the medium for various parameter combinations, the influence of each optical parameter is analyzed. In addition, the adaptability of the diffusion approximation method to the case of a gradient-index medium is examined. The results show that in the homogeneous-refractive-index medium and the gradient-index medium, the accuracy of diffusion approximation method is improved by an increase in the optical thickness or the scattering albedo and by a decrease in the wall emissivity. Moreover, the maximum relative errors are primarily distributed at the medium boundary and center. In the case of increased optical thickness or decreased wall emissivity, the error of the diffusion approximation method becomes stable and remains at a small value. Keywords Radiative transfer, diffusion approximation, gradient index, numerical analysis, relative error Date received: 16 March 2018; accepted: 1 October 2018 Handling Editor: Oronzio Manca Introduction Radiation transfer, as a basic transfer mode, occupies an important position in the field of industrial applications. The research is not only limited to the field of solar radiation transmission 1 and high-temperature gas radiation 2 but also the micro-scale radiative heat transfer, 3 spacecraft thermal control design, 4 and radiation inverse problem method to identify the temperature field of target 5 and so on. Obviously, the study of radiative transfer is still of great significance. At present, the calculation of the radiative transfer equation mainly depends on different numerical methods, such as Monte Carlo method (MCM), 6 discrete transfer method, 7 discontinuous finite element method, 8 natural element method, 9 and so on. However, these methods have different applicable conditions. For instances, the MCM can deal with all kinds of complex problems, the complexity of the simulation calculation increases proportionately with the complexity of the problem and the amount of calculation increases greatly when dealing with the radiation of non-uniform medium. The discrete transfer method is considered to be very convenient for dealing with homogeneous and inhomogeneous medium radiation, but it is not easy to converge when solving the energy equation. The discontinuous finite element method has strong adaptive School of Civil Engineering, Northeast Petroleum University, Daqing, People s Republic of China Corresponding author: Lijun Liu, School of Civil Engineering, Northeast Petroleum University, 99 Xuefu Street, Daqing , People s Republic of China. lljdqpi@126.com Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 2 Advances in Mechanical Engineering ability and high accuracy, but the computation is still slightly larger. The P1 diffusion approximation (DA) method, which is based on spherical harmonics, eliminates the angle coordinates and also reduces the dimension processing, which greatly simplifies this problem. In recent years, the DA method is mainly applied to the field of medical imaging because of its fast computation. Ruan 10 used diffusion equation to solve the forward problem in the shape diffusion optical tomography technology. Furthermore, Xu et al. 11 investigated total time-resolve diffusion reflectance from a semi-infinite turbid media based on the DA method. For the research on the gradient-index medium, to our knowledge, Liu and colleagues 12,13 have derived the DA equation of steady-state radiative transfer in multidimensional gradient-index medium, and the application conditions of the DA method are briefly described, so that the degree of influence of each optical parameter has not been further revealed. In this work, the impact of t L, e, and v on the accuracy of the DA method, for media having a homogeneous refractive index and a gradient index distribution, is analyzed and compared. By calculating the distribution of the radiation heat flux at the bottom of the medium, the degree of influence of each optical parameter is determined. Note that, in order to simplify this work, the two-dimensional (2D) steady-state radiative transfer equation is considered. Method The radiative transfer equation for participating media can be expressed as n 2 d I(r, O) ds n 2 = n 2 ki b (r) (k + ss) I(r, O) + s 1Þ s I(r, O 0 )F(O 0, O)dO 0 4p 4p where I(r, O) is the radiation intensity, n is the refractive index of the medium, s is the geometric path length, k is the absorbing coefficient, s s is the scattering coefficient,i b (r) is the blackbody radiation intensity, and F(O 0, O) is the scattering phase function. Based on P1 DA method, I can be rewritten in terms of the incident radiation G, and radiation heat flux q can be written as I(r, O)= 1 ½G(r)+3q(r) ^sš 4p 2Þ Combined with formulas (1) and (2), the steady-state radiative transfer equation can be obtained " # " # 1 1 r b 1 A rg(r) + r 1v 3 b 1 A rn2 1v n b(1 v)½g(r) 4pn 2 I b (r)š = 0 The boundary condition are expressed as 3Þ G = 2 e 2q n w pn 2 I bw 4Þ e where n w is the outer normal vector of the wall. Using the conventional Galerkin finite element method, finally we can obtain 8 X< rg i (r) rg j (r)d + (rlnn 2 (r) rg i (r))g j (r)d : i + (k½b(3 A 1 v)š+r 2 lnn 2 (r))g i (r)g j (r)d 9 rn 2 n 2 n w + 3 e = w b 0 G i (r)g j (r)da 22 e w ; G i = 4p½b(3 A 1 v)škn 2 I b (r)g j (r)d 3 e w b 0 4pn 2 I bw G j (r)da 5Þ 22 e w where G i is the shape function, and A represent the computing domains and the area, respectively, and b is the extinction coefficient. The discrete formula (5) can be written in a matrix form K ji G i = H j 6Þ Each matrix element is defined as 8 K ji = X < rg i (r) rg j (r)d : i + (r ln n 2 (r) rg i (r))g j (r)d + (k½b(3 A 1 v)š + r 2 ln n 2 (r))g i (r)g j (r)d 9 rn 2 n 2 n w + 3 e = w b 0 G i (r)g j (r)da 2 2 e w ; H j = 4p½b(3 A 1 v)škn 2 I b (r)g j (r)d 3 e w b 0 4pn 2 I bw G j (r)da 2 2 e w 7Þ 8Þ

3 Liu et al. 3 Figure 1. Dimensionless radiation heat flux distributions along the bottom wall for e = 1:0, v = 0:5 and (c)t L = 1:0, 5:0, and 10:0, respectively. (c) Results and discussion To determine the performance and accuracy of the presented DA method, we examine radiative heat transfer problems in a rectangle with opaque and diffuse surfaces. The boundaries are cold (0 K), and the medium temperature is maintained at 1000 K. For the following numerical study, the medium is uniformly divided into units. Depending on the t L, e, v, and n of the material, different results can be obtained. Several sample problems are tested to explore the accuracy and computation efficiency of P1 DA for different t L, e, and v. In subsection Square enclosure filled with homogeneous-refractive-index medium, three examples regarding homogeneous-refractive-index media are presented. In subsection Square enclosure filled with graded-index medium, examples concerning gradient-index medium are considered. Square enclosure filled with homogeneous-refractiveindex medium The DA algorithm is tested through comparison with the MCM. The examined medium has a fixed n(x)=5. For the case of e = 1:0 and v = 0:5, the dimensionless heat flux distributions at t L =1:0, 5:0, and 10:0 are shown in Figure 1 (c), respectively. It is apparent that the DA results near the boundaries are larger than those obtained by using the MCM when t L is small and that the relative error is approximately 12%. Furthermore, the deviation near the center is small, with a maximum relative error of only 3.4%. With increasedt L (t L = 5:0), the error is primarily distributed near the center of the medium, and the maximum relative error is approximately 6%. Finally, for t L = 10:0, the maximum relative error is less than 5.5%. Overall, the DA method produces large errors in optically thin media, and these errors decrease as t L increases. Aftert L reaches a certain value, the relative error is no longer significantly reduced with further increases to the thickness. For e = 0:5,t L = 10:0, and different v, the resultant dimensionless heat flux distribution are presented in Figure 2. This figure demonstrates that the maximum relative error appears at the wall in this case. Under the condition of v = 0, the maximum relative error of the DA method is 6.3%, which is higher than that obtained for v = 0:5. The reason is that high scattering albedo

4 4 Advances in Mechanical Engineering Figure 2. Dimensionless radiation heat flux distributions along the bottom wall for e = 0:5, t L = 10:0, and v = 0:5 and v =0. Figure 3. Dimensionless radiation heat flux distribution on homogeneous-refractive-index medium bottom for t L = 10:0, e = 0:1, and v = 0:5. makes wall radiation intensity with angle distribution more uniform. In Figure 3, the distributions of the dimensionless heat flux as obtained using both the DA method and the MCM are plotted for the case of t L = 10:0, e = 0:1, and v = 0:5. It can be observed from Figures 1(c), 2, and 3 that, for the same t L and v and with e = 1:0, 0:5, and 0:1, respectively, the maximum relative errors are 5.5%, 3%, and 0.53%, respectively. Therefore, the DA method exhibits higher accuracy for small e. This is because the lower the emissivity, the lower the absorption rate, and the more uniform the radiation intensity distribution. Square enclosure filled with graded-index medium Next, we consider a special layered graded-index medium, where the refractive index is modeled by n(x)=5 1 0:95 2 x 2 0:5 9Þ Figure 4. Dimensionless radiation heat flux distribution on graded-index medium bottom for t L = 1:0, e = 1:0, and v = 0:5. Thus, this medium has a variable refractive index n(x). We set t L = 1:0, e = 1:0, and v = 0:5, and the resultant dimensionless heat flux distributions are shown in Figure 4. Hence, it is apparent that the dimensionless heat flux density is asymmetric and the peak moves to the left in the gradient-index medium, and the uniform radiation intensity distribution cause the error at the left wall to increase, while that at the right wall decrease. Figure 5 and shows the dimensionless heat flux distributions for t L =5:0 and 10:0, respectively. The other conditions are as above. Comparison of these results with Figure 4, indicates that, for t L = 1:0, the DA method is inaccurate at the left wall, and the maximum relative error is 23%. However, with increased t L, the left-wall error decreases gradually. The maximum relative error of t L = 5:0 is approximately 5%. However, with further increases in t L, for example, in the case of t L = 10:0, the maximum relative error remains close to 5%. Thus, a threshold optical thickness t L0 exists, in the case of t L <t L0, the DA method exhibits low accuracy; otherwise, for higher accuracy is achieved.

5 Liu et al. 5 Figure 5. Dimensionless radiation heat flux distributions along the bottom wall for e = 1:0, v = 0:5, and t L = 5:0 and t L = 10:0. Figure 6. Dimensionless radiation heat flux distributions along the bottom wall for t L = 10:0, e = 0:5 and v = 0:5 and v = 0. In Figure 6, the dimensionless heat flux distributions are plotted for t L = 10:0, e = 0:5, and different v values of 0.5 and 0. It can observed that an increased value of v can improve the DA method accuracy. That is, the maximum relative error for v = 0:5 is 2.79%, whereas that for v = 0 is 3.61%. Because t L = 10:0 in this case, the DA method has higher accuracy, and the relative error is not large. Thus, for an optically thick media, increasing v can improve the accuracy of DA method; however, the effect of the increased v is not significant. Figure 7 shows the distributions of dimensionless heat flux under the conditions of t L = 10:0, e = 0:1, and v = 0:5. Comparison of this figure with Figures 5 and 6 demonstrate that, for the samet L and v, the maximum relative errors of DA method are approximately 5%, 2.8%, and 2.8% for e = 1:0, 0:5, and 0:1, respectively. Although the relative errors obtained from Figures 6 and 7 are essentially identical, some error is apparent near the medium center in Figure 6. In contrast, the DA method has a significant error in the Figure 7. Dimensionless radiation heat flux distributions on graded-index medium bottom for t L = 10:0, e = 0:1, and v = 0:5. small range near the left wall surface of the medium only in Figure (7), this is in good agreement with the benchmark solution in other ranges.

6 6 Advances in Mechanical Engineering Conclusion The DA method, which is a rapid algorithm, was extended to solve radiative transfer problems in a 2D rectangle containing a participating medium. Then, the applicability and precision of this method for a variety of optical parameters were examined. Under appropriate boundary conditions, the DA method was found to be consistent with the MCM, and the errors were primarily distributed at the boundaries and the medium center. Comparison of the influence of the optical parameters on the DA method in case of a homogeneousmedium-refractive-index and a graded-index medium indicated the accuracy of the DA. That is, the accuracy of the DA method increase with increased optical thickness or scattering albedo and with decreased wall emissivity. When the optical thickness increases or the wall emissivity decreases to a certain value, the error of the DA method becomes stable and keeps a small value. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) received no financial support for the research, authorship, and/or publication of this article. ORCID ids Lijun Liu Xiaoyan Liu References 1. Avezov RR, Avezova NR and Samiev KA. Transmission coefficient of the transparent insulation of flat diffuse solar radiation heating systems. Appl Sol Energ 2007; 43: iskanta R. Overview of convection and radiation in high temperature gas flows. Int J Eng Sci 1998; 36: Jing L. Micron/nano scale heat transfer. Beijing, China: Science Press, Hu J. Progress in thermal control technology for Chinese spacecraft. Spacecr Eng 2001; 1: Liu L, Tan H and Yu Q. Inverse radiation problem of temperature field in three-dimensional rectangular furnaces. Int Commun Heat Mass 1999; 26: Ruan L-M, Tan J-Y, Dong S-K, et al. The precision analysis about radiative transfer Monte-Carlo method and numerical experiment research. J Eng Thermophys 2003; 24: Nirgudkar H, Kumar S and Srivastava A. Solution of radiative transfer equation using discrete transfer method for two-dimensional participating medium. Int Commun Heat Mass 2015; 61: Liu LH and Liu LJ. Discontinuous finite element method for radiative heat transfer in semitransparent graded index medium. J Quant Spectrosc Ra 2007; 105: Yong Z. Numerical research on the radiative transfer in participant media by natural element method. Harbin, China: Harbin Institute of Technology, Ruan PQ. Investigations on shape-based diffuse optical tomography: theory and methodology. Tianjin, China: Tianjin University, Xu T, Zhang C, Wang X, et al. Accuracy of the diffusion approximation for total time resolved reflectance from a semi-infinite turbid medium. Optik Int J Light Electron Optics 2003; 114: Liu LJ and Zhao JM. Diffusion approximation of steady radiative transport in graded index media. Comput Phys 2013; 30: Liu LJ and Liu LH. P1 approximation of radiative transfer equation in gradient index medium. In: Proceedings of the Chinese Engineering Thermophysics Society Heat And Mass Transfer academic conference, Qingdao, China, October 2009.