
Accurate prediction of the angular and spatial distributions of radiative intensity is a very important and challenging issue for the coupled radiation and conduction problem with nonlinear anisotropic scattering medium. Different with the traditional hybrid spectral methods, spectral collocation method associated with discrete ordinate method (SCM-DOM), the spectral collocation method is extended to discretized both angular and spatial domains of governing equations in concentric cylinders. The angular and spatial derivative terms of governing equations in the cylindrical coordinate system are approximated by high order Chebyshev polynomials instead of the low order finite difference schemes. The performance of SCM is evaluated by comparing with available data in literature. Numerical results show that convergence rates of angular and spatial nodes approximately follow the exponential decaying law. In addition, for nonlinear anisotropic scattering medium, the SCM provides smoother results and mitigates the ray effect. The SCM is a successful and efficient method to deal with coupled radiative and conductive heat transfer in concentric cylinders. Furthermore, the effects of various geometric and thermal physical parameters on dimensionless temperature and heat flux are comprehensively investigated.
Citation: Yasong Sun, Jiazi Zhao, Xinyu Li, Sida Li, Jing Ma, Xin Jing. Prediction of coupled radiative and conductive heat transfer in concentric cylinders with nonlinear anisotropic scattering medium by spectral collocation method[J]. AIMS Energy, 2021, 9(3): 581-602. doi: 10.3934/energy.2021028
[1] | Xiaofeng Zhang, Qing Ai, Kuilong Song, Heping Tan . Bidirectional monte carlo method for thermal radiation transfer in participating medium. AIMS Energy, 2021, 9(3): 603-622. doi: 10.3934/energy.2021029 |
[2] | Caliot Cyril, Flamant Gilles . Pressurized Carbon Dioxide as Heat Transfer Fluid: In uence of Radiation on Turbulent Flow Characteristics in Pipe. AIMS Energy, 2014, 1(2): 172-182. doi: 10.3934/energy.2014.2.172 |
[3] | Hong-Yu Pan, Chuang Sun, Xue Chen . Transient thermal characteristics of infrared window coupled radiative transfer subjected to high heat flux. AIMS Energy, 2021, 9(5): 882-898. doi: 10.3934/energy.2021041 |
[4] | Chao-Jen Li, Peiwen Li, Kai Wang, Edgar Emir Molina . Survey of Properties of Key Single and Mixture Halide Salts for Potential Application as High Temperature Heat Transfer Fluids for Concentrated Solar Thermal Power Systems. AIMS Energy, 2014, 2(2): 133-157. doi: 10.3934/energy.2014.2.133 |
[5] | Cheng Ziming, Lin Bo, Shi Xuhang, Wang Fuqiang, Liang Huaxu, Shuai Yong . Influences of atmospheric water vapor on spectral effective emissivity of a single-layer radiative cooling coating. AIMS Energy, 2021, 9(1): 96-116. doi: 10.3934/energy.2021006 |
[6] | Vladimir A. Kaminsky, Nina Yu. Obvintseva, Svetlana A. Epshtein . The estimation of the kinetic parameters of low-temperature coal oxidation. AIMS Energy, 2017, 5(2): 163-172. doi: 10.3934/energy.2017.2.163 |
[7] | Hong-Wei Chen, Fu-Qiang Wang, Yang Li, Chang-Hua Lin, Xin-Lin Xia, He-Ping Tan . Numerical design of dual-scale foams to enhance radiation absorption. AIMS Energy, 2021, 9(4): 842-853. doi: 10.3934/energy.2021039 |
[8] | Simona Silvia Merola, Luca Marchitto, Cinzia Tornatore, Gerardo Valentino . Spray-combustion process characterization in a common rail diesel engine fuelled with butanol-diesel blends by conventional methods and optical diagnostics. AIMS Energy, 2014, 2(2): 116-132. doi: 10.3934/energy.2014.2.116 |
[9] | Wei An, Yifan Zhang, Bo Pang, Jun Wu . Synergistic design of an integrated pv/distillation solar system based on nanofluid spectral splitting technique. AIMS Energy, 2021, 9(3): 534-557. doi: 10.3934/energy.2021026 |
[10] | Thomas Wetzel, Julio Pacio, Luca Marocco, Alfons Weisenburger, Annette Heinzel, Wolfgang Hering, Carsten Schroer, Georg Muller, Jurgen Konys, Robert Stieglitz, Joachim Fuchs, Joachim Knebel, Concetta Fazio, Markus Daubner, Frank Fellmoser . Liquid metal technology for concentrated solar power systems: Contributions by the German research program. AIMS Energy, 2014, 2(1): 89-98. doi: 10.3934/energy.2014.1.89 |
Accurate prediction of the angular and spatial distributions of radiative intensity is a very important and challenging issue for the coupled radiation and conduction problem with nonlinear anisotropic scattering medium. Different with the traditional hybrid spectral methods, spectral collocation method associated with discrete ordinate method (SCM-DOM), the spectral collocation method is extended to discretized both angular and spatial domains of governing equations in concentric cylinders. The angular and spatial derivative terms of governing equations in the cylindrical coordinate system are approximated by high order Chebyshev polynomials instead of the low order finite difference schemes. The performance of SCM is evaluated by comparing with available data in literature. Numerical results show that convergence rates of angular and spatial nodes approximately follow the exponential decaying law. In addition, for nonlinear anisotropic scattering medium, the SCM provides smoother results and mitigates the ray effect. The SCM is a successful and efficient method to deal with coupled radiative and conductive heat transfer in concentric cylinders. Furthermore, the effects of various geometric and thermal physical parameters on dimensionless temperature and heat flux are comprehensively investigated.
Nomenclature: An: the coefficient of Legendre expansion of order n; B: backward scattering phase function; Di,j: the element of the first order derivative matrix; D(2)i,j: the element of the second order derivative matrix; er,eψ,ez: the spatial vectors; Er: integral averaged relative error; F: forward scattering phase function; G: incident radiative energy; h: Lagrange interpolation polynomials; I: radiative intensity; k: thermal conductivity; n: unit outward normal vector; Ncr: conduction-radiation parameter; Nr∗,Nφ,Nθ: number of collocation points in dimensionless radius, azimuthal angle and polar angle, respectively; Pn: Legendre polynomials of order n; qc: dimensionless conductive heat flux; qr: dimensionless radiative heat flux; qt: dimensionless total heat flux; r: radial coordinate of cylindrical coordinate system; r∗: dimensionless radius; Rref: available data from references; RSCM: numerical solution by SCM; S: source term; T: temperature; wθ,wφ: quadrature weight in polar angle and azimuthal angle;
Greek Symbols: αr∗,αθ,αφ: standard computational domain in dimensionless radius, polar angle and azimuthal angle, respectively; β: extinction coefficient; δ: standard deviation; ε: emissivity of boundary surface; η: direction cosine in eψ direction; θ, θ′: polar angle; Θ: dimensionless temperature; κa: absorption coefficient; κs: scattering coefficient; κt: transmission coefficient; μ: direction cosine in er direction; ξ: direction cosine in ez direction; σ: Stefan Boltzmann constant; τ: optical thickness; φ, φ′: azimuthal angle; Φ(Ω,Ω′): scattering phase function from the incident direction Ω′ to the scattering direction Ω; ψ: dimensionless radiative intensity; Ψ: scattering angle; ω: scattering albedo; Ω, Ω′: the direction of radiative intensity;
Subscripts: B3: a typical backward scattering; b: black body radiative intensity; F3: a typical forward scattering; i,j,k,m,m′: solution node indexes; in: inner wall of concentric cylinders; isotropic scattering: isotropic scattering; out: outer wall of concentric cylinders; r∗in: value at inner wall; r∗out: value at outer wall; w: value at wall; αr∗: value in standard radius computational domain; αφ: value in standard azimuthal angle computational domain;
Superscripts: CG: Chebyshev-Guass points; CGL: Chebyshev-Guass-Lobatto points; m,m′;n,n′: angular direction of radiation;
The combined radiative-conductive heat transfer in participating medium [1] plays a dominant role in high temperature equipment, such as aeroengine combustor, nuclear reactor and industry furnaces, etc. Accurate prediction of temperature and heat flux requires solving radiative transfer equation (RTE) and energy equation simultaneously [2,3,4,5,6]. Different with RTE in Cartesian coordinate, RTE in cylindrical coordinate exists the angular derivative term which would increase the mathematical complexity.
In recent years, the coupled radiative-conductive heat transfer in participating medium of cylindrical geometry has evoked wide interests of many researchers. As early as 1982, Fernandes and Francis [7] gave the rigorous formulations of combined conduction and radiation in concentric cylinders and numerically solved by Galerkin finite element method. Pandey [8] employed undetermined parameters method to solve this coupled problem for gray and nongray gases contained between infinitely long concentric cylinders with black surfaces. Krishnaprakas [9] used the hybrid strategy to analyze combined conduction and radiation in cylindrical geometries. In this paper, energy equation was solved by finite different method, and RTE was solved by discrete ordinates method in conjunction with Crank-Nicolson scheme. The effects of thermal-physical parameters, namely emittance, scattering albedo, scattering phase function, conduction-radiation parameters on heat fluxes were investigated. Dlala et al. [10] investigated coupled radiative-conductive heat transfer in gray hollow spheres and cylinders. They used finite Chebyshev transform (FCT) to improve the performance of discrete ordinates method, and adopted Chebyshev polynomials to approximate the angular derivative term instead of finite difference scheme. The FCT was more accurate than traditional discrete ordinate method. Mishra et al. [11,12] developed the modified discrete ordinate method and lattice Boltzmann method to analyze coupled radiative-conductive heat transfer in infinite and finite concentric cylinders with absorbing, emitting, and scattering medium. Authors claimed this modified discrete ordinate method was not require complicated and intensive calculation to determine the discrete directions and directional weights, allowed freedom of direction selection. Zhou et al. [13] extended the MDOM, which based on superposition technique and considering the contributions of the walls and medium, to the cylindrical medium. And the results showed that the cost computational time was comparable to DOM and the ray effect can be mitigated effectively.
Different with above numerical algorithms, spectral methods [14,15,16] are useful tools to solve ordinary differential equations or partial differential equations with high order accuracy, and usually the best choice for solving problem with smooth solution [17]. Benefiting from high accuracy, simple implementation, and exponential convergence characteristic, they have been widely used to solve problems in many fields, such as computational fluid dynamics [18,19], magnetohydrodynamics [20], and optics [21]. Using spectral methods to solve thermal radiative heat transfer already has a history of two decades [22,23,24,25,26]. Among these researches, several hybrid spectral methods such as spectral collocation method associated with discrete ordinate method (SCM-DOM) [22,23,24] and spectral element method combined with discrete ordinate method (SEM-DOM) [25] were developed, in which the spatial domain was discretized by spectral methods, and the angular domain was discretized by DOM. Recently, Wang et al. [26] further extended the SCM-DOM for solving polarized radiative transfer problems in multi-layered participating media.
Different with the above hybrid spectral methods, Zhou et al. [27,28] taken advantage of SCM to discretize the entire spatial and angular domain rather than only the spatial domain. Their numerical test showed that SCM can achieve the high accuracy both in spatial and angular directions. They further developed SCM to solve radiative integro-differential transfer equation [29] in one-dimensional medium which only contains unknown radiative flux and eliminate the effect of angular derivative term. However, in common situations, radiative heat transfer is coupled to other models of heat transfer. To the best of authors' knowledge, no research to date has aimed to using SCM to solving coupled radiative and conductive heat transfer so far.
The objective of this work is to extend the SCM to solve coupled radiative and conductive heat transfer in concentric cylindrical medium. This paper is organized as follow: In section 2, the physical and mathematical models of coupled radiative and conductive heat transfer in concentric cylinders are presented. In section 3, the SCM formulations of RTE and energy equation are deduced. In section 4, the performance of SCM is verified from the available data in the literature. In section 5, the effects of various geometric and thermo-physical parameters are comprehensively investigated. Finally, the conclusions are summarized in section 6.
As shown in Figure 1, the present study considers the coupled radiative and conductive heat transfer processing in cylindrical coordinate system. The absorbing, emitting and anisotropic scattering medium is filled in concentric cylinders.
In concentric infinite cylinders, the non-conservation form of RTE for a gray medium is [30,31,32]
μ∂I(r,Ω)∂r−ηr∂I(r,Ω)∂φ+(κa+κs)I(r,Ω)=κaIb(r)+κs4π∫4πI(r,Ω′)Φ(Ω′,Ω)dΩ′ | (1) |
where I(r,Ω) is the radiative intensity at spatial position r along angular direction Ω; The direction Ω can be expressed by the direction cosines μ=sinθcosφ, η=sinθsinφ and ξ=cosθ, where θ is polar angle and φ is azimuthal angle; κa and κs are absorption coefficient and scattering coefficient, respectively. The anisotropic scattering phase function Φ(Ω′,Ω) represents the probability that a radiative beam along angular direction Ω′ is scattered to angular direction Ω, and is approximated by a finite series of Legendre polynomials as
Φ(Ω′,Ω)=1+N∑n=1AnPn(cosΨ) | (2) |
where An is the coefficient of Legendre expansion of order n which are listed in Table 1 [33]; Pn is the Legendre polynomials; Ψ is the included angle between the incident direction Ω′ and the scattering direction Ω.
coefficients | scattering phase functions | |||||
F1 | F2 | F3 | B1 | B2 | B3 | |
A0 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 |
A1 | 2.53602 | 2.00917 | 1.00000 | −0.56524 | −1.20000 | −1.00000 |
A2 | 3.56549 | 1.56339 | 0.29783 | 0.50000 | ||
A3 | 3.97976 | 0.67407 | 0.08571 | |||
A4 | 4.00292 | 0.22215 | 0.01003 | |||
A5 | 3.66401 | 0.04725 | 0.00063 | |||
A6 | 3.01601 | 0.00671 | ||||
A7 | 2.23304 | 0.00068 | ||||
A8 | 1.30251 | 0.00005 | ||||
A9 | 0.53463 | |||||
A10 | 0.20136 | |||||
A11 | 0.05480 | |||||
A12 | 0.01099 |
For the gray, opaque and diffuse boundary, the boundary conditions of RTE can be expressed as
I(rw,Ω)=εwIb,w+1−εwπ∫nw⋅Ω′<0I(rw,Ω′)|nw⋅Ω′|dΩ′,nw⋅Ω′⩾0 | (3) |
where εw is the wall emissivity, Ib,w is the blackbody radiative intensity at the wall, nw is unit wall normal.
As shown in Figure 2, the following symmetric condition is satisfied
I(r,θ,φ)=I(r,π−θ,φ)=I(r,θ,2π−φ)=I(r,π−θ,2π−φ) | (4) |
Thus, Eq (4) is solved on the three-dimensional domain (r,θ,φ)∈(rin,rout)×(0,π/2)×(0,π).
For the infinite concentric cylinders, the steady-state energy equation can be written as
dqrdr+1rqr=4πκak{Ib[T(r)]−14πG(r)} | (5) |
with the boundary conditions
{T(rin)=TinT(rout)=Tout | (6) |
where G is the incident radiative energy,
G(r)=4∫π0∫π20Isinθdθdφ | (7) |
For convenience of analysis, the following dimensionless parameters [10] are introduced
Θ=TTin,r∗=rrout,ψ=πIσT4in,Ncr=kβ4σT3in,ω=κsβ,τout=βrout,qr=∫4πψμdΩ,qc=−dΘdr∗ | (8) |
Then, RTE, energy equation and the corresponding boundary conditions can be transformed into the dimensionless forms as
μτout∂ψ∂r∗−ητoutr∗∂ψ∂φ=−ψ+(1−ω)Θ4(r∗)+ω4π∫4πψ(r∗,Ω′)Φ(Ω,Ω′)dΩ′ | (9) |
{ψ(r∗in)=εinΘ4in+1−εinπ∫nin⋅Ω′>0ψ(r∗in,Ω′)|nin⋅Ω′|dΩ′,nin⋅Ω′<0ψ(r∗out)=εoutΘ4out+1−εoutπ∫nout⋅Ω′>0ψ(r∗out,Ω′)|nout⋅Ω′|dΩ′,nout⋅Ω′<0 | (10) |
d2Θdr∗2+1r∗dΘdr∗=(1−ω)Ncrτ2out(Θ4−14π∫4πψdΩ) | (11) |
{Θ(r∗in)=1Θ(r∗out)=Θout | (12) |
As shown in Figure 3, the spatial-angular domain [rin,rout]×[0,π/2)×[0,π) is discretized into {r∗1,r∗2,⋯,r∗Nr∗}, {θ1,θ2,⋯,θNθ}, and {φ1,φ2,⋯,φNφ} along r∗, θ and φ directions, respectively. According the theory of SCM, the discretized spatial-angular domain should be transferred to the standard Chebyshev domain by the following relationship
{r∗i=1−r∗in2αr∗,i+1+r∗in2θj=π4(αθ,j+1)φk=π2(αφ,k+1) | (13) |
where αi is the Gauss-Lobatto points.
Then, dimensionless radiative intensity and dimensionless temperature can be approximated by Lagrange interpolation polynomials and collocation points [16]
ψ(r∗,θ,φ)≈Nr∗∑i=1Nθ∑j=1Nφ∑k=1ψi,j,khi(r∗)hj(θ)hk(φ) | (14) |
Θ(s)≈N∑i=1Θihi(s) | (15) |
where hi are Lagrange interpolation polynomials.
Substituting Eq (14) into Eq (9), weighting by weight function and integrating over the computational domain, Eq (9) can be discretized as
2μm,nτout(1−r∗in)Nr∗∑j=0DCGLαr∗,i,jψm,nj−2ηm,nπτoutr∗Nφ∑m′=0DCGαφ,m,m′ψm′,ni+ψm,ni=(1−ω)Θ4(r∗i)+ωπ8Nφ∑m′=0Nθ∑n′=0ψm′,n′iΦ(Ω,Ω′)sinθn′wn′θwm′φi=0,1,...,Nr∗ | (16) |
The corresponding boundary conditions are discretized as
{ψm,nr∗in=εinψb,in+π(1−εin)2Nφ∑m′=0,μm′,n′<0Nθ∑n′=0ψm′,n′r∗in|μm′,n′|sinθn′wn′θwm′φ,μm′,n′>0ψm,nr∗out=εoutψb,out+π(1−εout)2Nφ∑m′=0,μm′,n′>0Nθ∑n′=0∑ψm′,n′r∗out|μm′,n′|sinθn′wn′θwm′φ,μm′,n′<0 | (17) |
Similarly, the energy equation is discretized as
(21−r∗in)2Nr∗∑j=0D(2)CGLαr∗,i,jΘj+1r∗21−r∗inNr∗∑j=0DCGLαr∗,i,jΘj=(1−ω)Ncrτ2out(Θ4i−π8Nφ∑m′=0Nθ∑n′=0ψm′,n′isinθn′wn′θwm′φ) | (18) |
where Di,j and D(2)i,j [14] are the first and second order derivative matrix, respectively.
Figure 4 shows the flow chart of the SCM for the coupled radiation-conduction problem in concentric cylinders.
Based on the above described SCM model for coupled radiative and conductive heat transfer in concentric cylinders with participating medium. In the following, several test cases are adopted to verify the performance of SCM model. Compared with available data in references, the accuracy and efficiency of SCM for coupled radiative-conductive heat transfer in concentric cylinders are validated.
In order to quantitatively evaluate the accuracy of SCM, the integral averaged relative error is defined as
Er=∫|RSCM(r∗)−Rref(r∗)|dr∗∫|Rref(r∗)|dr∗×100% | (19) |
where Rref is available data from references.
A piratical case of coupled radiative-conductive heat transfer in concentric cylinders is considered with the conduction-radiation parameter Ncr=0.01, the ratio of inner and outer cylinders r∗in/r∗out=0.5 and dimensionless temperature at outer surface Θout=0.1. There are blackbody surfaces and non-scattering medium. The extinction coefficient is β=1. This case has also been adopted by Mishra et al. [34] for lattice Boltzmann method associated with finite volume method (LBM-FVM).
The distribution of dimensionless temperature within concentric cylinders by SCM is plotted in Figure 5, and compared to LBM-FVM results. The SCM results is very close to those of the LBM-FVM results, and the integral average relative error is 0.875%.
Furthermore, Table 2 lists the dimensionless total heat fluxes at inner and outer surfaces for different conduction-radiation parameters Ncr and scattering albedo ω. In this table, the results in first and second column are copied from Ref. [10], and obtained by DOM and FCT, respectively. Different with the first case, dimensionless temperature at outer surface is Θout=0.5. It can be seen that for both conduction dominated (Ncr=1.00) and radiation dominated (Ncr=0.01) situations, the results of SCM and FCT are in good agreement with each other. The maximum relative error between SCM and FCT is 1.023%.
Ncr | ω | qt(r∗in) | qt(r∗out) | ||||
DOM [10] | FCT [10] | SCM | DOM [10] | FCT [10] | SCM | ||
1 | 0.9 | 1.6460 | 1.6436 | 1.6421 | 0.8230 | 0.8218 | 0.8210 |
0.5 | 1.6502 | 1.6488 | 1.6468 | 0.8251 | 0.8244 | 0.8234 | |
0.1 | 1.6542 | 1.6537 | 1.6512 | 0.8271 | 0.8268 | 0.8256 | |
0.1 | 0.9 | 3.4764 | 3.4523 | 3.4363 | 1.7382 | 1.7261 | 1.7183 |
0.5 | 3.5181 | 3.5045 | 3.4840 | 1.7592 | 1.7522 | 1.7422 | |
0.1 | 3.5578 | 3.5529 | 3.5271 | 1.7789 | 1.7763 | 1.7638 | |
0.01 | 0.9 | 21.7839 | 21.5403 | 21.3807 | 10.8919 | 10.7700 | 10.6921 |
0.5 | 22.1593 | 21.9907 | 21.7937 | 11.0796 | 10.9953 | 10.8988 | |
0.1 | 22.4352 | 22.3172 | 22.0889 | 11.2176 | 11.1586 | 11.0465 |
Then, we would verify the accuracy in the whole domains by testing whether the condition r∗q∗t=constant is satisfied. This is because as long as the equilibrium is satisfied, the energy equation always writes div[q∗t(r)]=0 where q∗t is the dimensionless total heat flux and r is spatial position. For a one-dimensional system, spatial position r is simplified as radial distance r and the divergence equation in cylindrical coordinates simply writes ddr(r∗q∗t)=0 which means r∗q∗t=constant. In Table 3, the results obtained with different collocation point numbers ¯r∗q∗r±δ are given, where the standard deviation δ=√1N∑Ni=1(¯r∗q∗t−r∗q∗t)2, the ¯r∗q∗t is the mean value of the r∗q∗t, i.e, ¯r∗q∗t=1N∑Ni=1r∗q∗t.
Nr | Nφ | Nθ | ¯r∗q∗t±δ |
12 | 28 | 28 | 7.342084±7.99×10−4 |
20 | 7.341663±8.49×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341664±8.56×10−5 | ||
44 | 7.341664±8.57×10−5 | ||
28 | 12 | 28 | 7.382330±1.71×10−3 |
20 | 7.351596±2.74×10−4 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.336948±3.89×10−5 | ||
44 | 7.334239±2.21×10−5 | ||
28 | 28 | 12 | 7.341663±8.53×10−5 |
20 | 7.341663±8.53×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341663±8.53×10−5 | ||
44 | 7.341663±8.53×10−5 |
It can be seen that in all case of different number of collocation points, r∗q∗t remains almost constant everywhere in the medium, with a very small deviation insuring four or five correct digits. Therefore, the accuracy of SCM is verified everywhere in the domain.
In order to the efficiency of this SCM model, Figure 6 depicts the effect of collocation point numbers on the integral averaged relative error. As shown in Figure 6, the horizontal axes are the number of radial points, azimuthal angle points and polar angle points, respectively. The vertical axis is the integral averaged relative error for the case of Ncr=0.03, r∗in/r∗out=0.5, εin=εout=1, Θout=0.1, β=2 and ω=0.5. When the number of radial points is near to 30, the integral averaged relative error is less than 1e-6 (Er⩽10−6). For Nr∗<30, the integral averaged relative error decreases very fast and approximately follows the exponential law with the increasing of the number of radial points. The similar trends are also found for the numbers of azimuthal angle points and polar angle points.
In order to comprehensively analyze this coupled heat transfer, the effects of geometric and thermo-physical parameters on dimensionless temperature and heat flux are investigated. In the subsection, the effects of different kinds of scattering phase functions are firstly studied, and then investigate the effect of various geometric and thermal physical parameters based on nonlinear anisotropic F1 scattering phase function.
Figure 7 presents the effect of scattering phase function on dimensionless temperature distribution. In these scattering phase functions, F3 scattering phase function ΦF3(Ω,Ω′)=1+cosΨ is a typical forward scattering phase function, B3 scattering phase function ΦB3(Ω,Ω′)=1−cosΨ is a typical backward scattering phase function. Figure 8 shows the dimensionless temperature distribution for three kinds of typical cases, F3 phase function, B3 phase function and isotropic scattering phase function Φ(Ω,Ω′)=1. Other parameters are fixed as Ncr=0.01, r∗in/r∗out=0.5, εin=εout=1, Θout=0.1, β=2 and ω=0.5. When r∗ is small, namely the region near the inner wall with high temperature. The dimensionless temperature ΘB3>Θisotropicscattering>ΘF3. While in the region near the outer wall with low temperature, the dimensionless temperature of the three kinds of scattering phase function is almost the same.
As shown in Figures 9 and 10, these phenomena can be explained by studying the dimensionless temperature for pure radiation problem (Ncr=0) and the ratio of radiative heat flux to total heat flux qr/(qr+qc) for this coupled radiation-conduction problem.
Firstly, the region near the high temperature inner wall is investigated. As shown in Figure 9, the dimensionless temperature distribution ΘB3>Θisotropicscattering>ΘF3. The reason is that the forward scattering phase function means more forward transmit power while the backward scattering phase function means more backward scattering energy. At the same time, Figure 10 shows that the ratio of radiative heat flux to total heat flux in the region near the inner wall is more than 0.9. Therefore, for the case of forward scattering phase function F3, the temperature near the inner wall is higher than that of the case of isotropic scattering. Meanwhile, compared with the case of isotropic scattering, the case of backward scattering phase function B3 leads to a lower temperature.
Secondly, the region near the outer wall with the low temperature is analyzed. As shown in Figure 9, the higher temperature distribution near the inner wall leads to the lower temperature distribution near the outer wall, namely ΘF3>Θisotropicscattering>ΘB3. However, as shown in Figure 10, the ratio of radiative heat flux to total heat flux near the outer wall is smaller than that of other regions. This means the influence of radiation is reduced near the outer wall. Consequently, the increase in temperature caused by the forward scattering phase function and the decrease in temperature caused by the backward scattering phase function are both weakened. Therefore, the dimensionless temperature distributions of the three kinds of scattering phase functions are almost the same in the regions near the outer wall.
As shown in Figure 11, the dimensionless temperature distribution tends to be linear with the increasing of conduction-radiation parameter Ncr from 0.01 to 10. The conduction-radiation parameter is defined as the ratio of conductive heat transfer and radiation heat transfer. It is obvious that, for the large value of Ncr, conduction plays a dominant role in this coupled heat transfer problem. On the contrary, the small value of Ncr means radiation becomes much more pronounced. For the conduction-dominated problem, the variation of dimensionless temperature tends to linear. Meanwhile, the variation of dimensionless temperature tends to nonlinear for the radiation-dominated problem. Thus, with the increasing of conduction-radiation parameter, the distribution of dimensionless temperature tends to linear.
For the cases of Ncr=1 (Figure 12a) and Ncr=0.1 (Figure 12b), the dimensionless temperature distributions are rising all the time with the increasing of radius ratio r∗in/r∗out. However, for the case of Ncr=0.01, there are two different trends of dimensionless temperature distributions in Figure 12c and d. In Figure 12c, the dimensionless temperature increases as the radius ratio changes from 0.1 to 0.5. But, in Figure 12d, the dimensionless temperature shows a downward trend instead of continuing to rise, and the dimensionless temperature profiles tends to be linear.
These phenomena can be explained by studying the radiative heat flux and the ratio of radiative heat flux to total heat flux, which are depicted in Figures 13 and 14, respectively.
In the three cases of Ncr=1,0.1,0.01, the radiative heat flux would augment steady with the increasing of the radius ratio. But the ratio of radiative heat flux to total heat flux firstly increases with radius ratio changing from 0.1 to 0.5, and then decreases with radius ratio increasing from 0.5 to 0.9.
It's worth to note that, for the case of Ncr=1 and Ncr=0.1, the change trend of dimensionless temperature is similar with that of radiative heat flux, namely, rise steady with the increasing of the radius ratio. However, for the case of Ncr=0.01, the change trend of dimensionless temperature is more like that of the ratio of radiative heat flux to total heat flux (increase first and then decrease).
In the cases of Ncr=1 and Ncr=0.1, conductive heat transfer plays a dominant role compared with radiative heat transfer. Therefore, the profile of the dimensionless temperature distribution tends to be linear. And the increase of radiative heat flux would raise the dimensionless temperature. However, the lower conduction-radiation parameter (Ncr=0.01) means the radiative heat transfer plays a dominant part. When the radius ratio changes from 0.1 to 0.5, the radiative heat flux is increasing. Consequently, the dimensionless temperature would rise and the profile presents the nonlinear characteristics. When the radius ratio keeps increasing from 0.5 to 0.9, the radiative heat flux continues to increase, but the ratio of radiative heat flux to total heat flux decrease. Considering the dominant role of radiative heat transfer, the decrease of ratio would definitely change the profile of temperature distribution, namely the profile tends to be linear which shows the increasing influence of conductive heat transfer. These are reasons why the dimensionless temperature rises firstly, then decrease and tends to be linear in the case of Ncr=0.01.
As shown in Figure 15, the dimensionless temperature decreases with the increasing of scattering albedo, and the distribution of dimensionless temperature also tends to be linear. The scattering albedo means the relative magnitude of absorption coefficient and scattering coefficient. ω=0 means the no-scattering medium, and ω=1 indicates the pure scattering medium. When scattering albedo approaches to 1, the scattering becomes stronger. This means the less radiant energy is absorbed. Consequently, the distribution of dimensionless temperature shows a downward trend. Meanwhile, for the case of a larger scattering albedo, the less radiant energy is absorbed means that conduction dominates a larger part in this coupled problem. Therefore, the distribution of dimensionless temperature tends to be linear.
SCM is developed to solve the coupled radiation-conduction problem in concentric cylinders with absorbing, emitting and nonlinear anisotropic scattering medium. Both the two involved RTE and steady-state energy equation are solved by SCM. In the solving process, the spatial and angular domains of RTE, and the spatial domain of energy equation are discretized by high order Chebyshev polynomials and Chebyshev collocation points. Compared with available data in references, accuracy and efficiency of the SCM for the coupled radiative-conductive heat transfer are validated. The high order accuracy can be obtained in a few nodes, and the exponential convergence characteristic of SCM exists in both spatial and angular domains. Considering that the SCM can obtain the high order accuracy and exponential convergence rate, the SCM model is an efficient model to solve the coupled radiative-conductive heat transfer in concentric cylinders with nonlinear anisotropic scattering medium. Besides, the effects of scattering phase function, conduction-radiation parameter, radius ratio and scattering albedo on dimensionless temperature and heat flux are comprehensively investigated.
This work was supported by the National Natural Science Foundation of China (Nos. 51976173, 51976014), Key Research and Development Program of Shaanxi (No. 2018SF-387), and the Fundamental Research Funds for the Central Universities (Nos. 3102019OMS701, 3102020OSC702).
We wish to confirm that there are no known conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.
[1] |
Howell JR, Menguc MP, Siegel R (2015) Thermal Radiation Heat Transfer, sixth ed., New York: CRC Press. doi: 10.1201/b18835
![]() |
[2] |
Kim TY, Baek SW (1991) Analysis of combined conductive and radiative heat-transfer in a 2-dimensional rectangular enclosure using the discrete ordinates method. Int J Heat Mass Transfer 34: 2265-2273. doi: 10.1016/0017-9310(91)90052-G
![]() |
[3] | Li ZH, Li XL, Xia XL, et al. (2019) A hybrid strategy for solving radiation-conduction in irregular geometries filled with gray semitransparent medium using Monte Carlo method combined with blocked-off and embedded boundary treatments. Numer Heat Transfer B 77: 22-41. |
[4] |
Razzaque MM, Howell JR, Klein DE (1984) Coupled radiative and conductive heat transfer in a two-dimensional rectangular enclosure with gray participating media using finite elements. J Heat Transfer 106: 613-619. doi: 10.1115/1.3246723
![]() |
[5] |
Sun YJ, Zhang XB (2018) A hybrid strategy of lattice Boltzmann method and finite volume method for combined conduction and radiation in irregular geometry. Int J Heat Mass Transfer 121: 1039-1054. doi: 10.1016/j.ijheatmasstransfer.2018.01.067
![]() |
[6] |
Bouzgarrou F, Askri F, Ali HB, et al. (2017) Analyses of unsteady conduction-radiation heat transfer using unstructured Lattice Boltzmann method. Int J Therm Sci 116: 287-309. doi: 10.1016/j.ijthermalsci.2017.03.002
![]() |
[7] |
Fernandes R, Francis J (1982) Combined conductive and radiative heat transfer in and absorbing, emitting, and scattering cylindrical medium. J Heat Transfer 104: 594-601. doi: 10.1115/1.3245173
![]() |
[8] |
Pandey DK (1989) Combined conduction and radiation heat transfer in concentric cylindrical media. J Thermophys 3: 75-82. doi: 10.2514/3.128
![]() |
[9] |
Krishnaprakas CK (1998) Combined conduction and radiation heat transfer in a cylindrical medium. J Thermophys 12: 605-608. doi: 10.2514/2.6385
![]() |
[10] |
Dlala NA, Sghaier T, Seddiki E (2007) Numerical solution of radiative and conductive heat transfer in concentric spherical and cylindrical media. J Quant Spectrosc Radiat Transfer 107: 443-457. doi: 10.1016/j.jqsrt.2007.02.012
![]() |
[11] |
Mishra SC, Krishna CH (2011) Analysis of radiative transport in a cylindrical enclosure--an application of the modified discrete ordinate method. J Quant Spectrosc Radiat Transfer 112: 1065-1081. doi: 10.1016/j.jqsrt.2010.11.011
![]() |
[12] |
Mishra SC, Krishna CH, Kim MY (2011) Analysis of conduction and radiation heat transfer in a 2D cylindrical medium using the modified discrete ordinate method and the lattice Boltzmann method. Numer Heat Transfer A-Appl 60: 254-287. doi: 10.1080/10407782.2011.588581
![]() |
[13] |
Zhou RR, Li BW (2019) The modified discrete ordinates method for radiative heat transfer in two-dimensional cylindrical medium. Int J Heat Mass Transfer 139: 1018-1030. doi: 10.1016/j.ijheatmasstransfer.2019.05.071
![]() |
[14] | Trefethen LN (2000) Spectral Methods in MATLAB. Philadelphia: Society for Industrial and Applied Mathematics. |
[15] | Canuto C, Hussaini MY, Quarteroni A, et al. (2006) Spectral Methods: Fundamentals in Single Domains, Berlin: Springer. |
[16] | Shen J, Tang T, Wang LL (2011) Spectral Methods: Algorithms, Analysis and Applications, Berlin: Springer. |
[17] | Shen J, Tang T (2006) Spectral and High-Order Methods with Applications, Beijing: Science Press. |
[18] | Peyret R (2002) Spectral Methods for Incompressible Viscous Flow, New York: Springer. |
[19] |
Canuto C, Hussaini MY, Quarteroni A, et al. (2007) Spectral Methods: Evolution to Complex Geometries and Applications to Fluid Dynamics, Berlin: Springer. doi: 10.1007/978-3-540-30728-0
![]() |
[20] |
Kornet K, Potherat A (2015) A method for spectral DNS of low Rm channel flows based on the least dissipative modes. J Comput Phys 298: 266-279. doi: 10.1016/j.jcp.2015.05.018
![]() |
[21] |
Abdrabou A, Heikal AM, Obayya SSA (2016) Efficient rational Chebyshev pseudo-spectral method with domain decomposition for optical waveguides modal analysis. Opt Express 24: 10495-10511. doi: 10.1364/OE.24.010495
![]() |
[22] |
Li BW, Sun YS, Yu Y (2008) Iterative and direct Chebyshev collocation spectral methods for one-dimensional radiative heat transfer. Int J Heat Mass Transfer 51: 5887-5894. doi: 10.1016/j.ijheatmasstransfer.2008.04.048
![]() |
[23] |
Sun YS, Li BW (2009) Chebyshev collocation spectral method for one-dimensional radiative heat transfer in graded index media. Int J Thermal Sci 48: 691-698. doi: 10.1016/j.ijthermalsci.2008.07.003
![]() |
[24] |
Li GJ, Ma J, Li BW (2015) Collocation spectral method for the transient conduction-radiation heat transfer with variable thermal conductivity in two-dimensional rectangular enclosure. J Heat Transfer 137: 032701. doi: 10.1115/1.4029237
![]() |
[25] |
Zhao JM, Liu LH (2007) Spectral element approach for coupled radiative and conductive heat transfer in semitransparent medium. J Heat Transfer 129: 1417-1424. doi: 10.1115/1.2755061
![]() |
[26] |
Wang CH, Feng YY, Yang YH, et al. (2020) Chebyshev collocation spectral method for vector radiative transfer equation and its applications in two-layered media. J Quant Spectrosc Radiat Transfer 243: 106822. doi: 10.1016/j.jqsrt.2019.106822
![]() |
[27] |
Zhou RR, Li BW (2017) Chebyshev collocation spectral method to solve radiative transfer equation in one-dimensional cylindrical medium. Int J Heat Mass Transfer 111: 1206-1217. doi: 10.1016/j.ijheatmasstransfer.2017.04.094
![]() |
[28] | Zhou RR, Li BW, Sun YS (2020) Predicting radiative heat transfer in axisymmetric cylindrical enclosures using the collocation spectral method. Int Commun Heat Mass Transfer 243: 106822. |
[29] |
Zhou RR, Li BW (2017) Chebyshev collocation spectral method for one-dimensional radiative heat transfer in linearly anisotropic scattering cylindrical medium. J Quant Spectrosc Radiat Transfer 189: 206-220. doi: 10.1016/j.jqsrt.2016.11.021
![]() |
[30] | Modest MF (2013) Radiative Heat Transfer, New York: Academic Press. |
[31] |
Chui EH, Raithby GD, Hughes P (1992) Prediction of radiative transfer in cylindrical enclosures with the finite volume method. J Thermophys Heat Transfer 6: 605-611. doi: 10.2514/3.11540
![]() |
[32] | Zhao JM, Liu LH (2018) Radiative transfer equation and solutions, Berlin: Springer, 933-978. |
[33] |
Kim TK, Lee H (1988) Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular enclosures. Int J Heat Mass Transfer 31: 1711-1721. doi: 10.1016/0017-9310(88)90283-9
![]() |
[34] |
Mishra SC, Kim MY, Das R, et al. (2009) Lattice Boltzmann method applied to the analysis of transient conduction radiation problems in a cylindrical medium. Numer Heat Transfer A-Appl 56: 42-59. doi: 10.1080/10407780903107162
![]() |
1. | Jiazi Zhao, Yasong Sun, Yifan Li, Changhao Liu, Investigation of coupled radiation-conduction heat transfer in cylindrical systems by discontinuous spectral element method, 2022, 6, 2515-3080, 354, 10.33737/jgpps/156350 | |
2. | Shreen El-Sapa, Noura S. Alsedais, Effect of slippage on a translational motion of two interacting non-concentric spheres squeezed by couple stress fluid, 2024, 0019-5588, 10.1007/s13226-024-00683-4 |
coefficients | scattering phase functions | |||||
F1 | F2 | F3 | B1 | B2 | B3 | |
A0 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 |
A1 | 2.53602 | 2.00917 | 1.00000 | −0.56524 | −1.20000 | −1.00000 |
A2 | 3.56549 | 1.56339 | 0.29783 | 0.50000 | ||
A3 | 3.97976 | 0.67407 | 0.08571 | |||
A4 | 4.00292 | 0.22215 | 0.01003 | |||
A5 | 3.66401 | 0.04725 | 0.00063 | |||
A6 | 3.01601 | 0.00671 | ||||
A7 | 2.23304 | 0.00068 | ||||
A8 | 1.30251 | 0.00005 | ||||
A9 | 0.53463 | |||||
A10 | 0.20136 | |||||
A11 | 0.05480 | |||||
A12 | 0.01099 |
Ncr | ω | qt(r∗in) | qt(r∗out) | ||||
DOM [10] | FCT [10] | SCM | DOM [10] | FCT [10] | SCM | ||
1 | 0.9 | 1.6460 | 1.6436 | 1.6421 | 0.8230 | 0.8218 | 0.8210 |
0.5 | 1.6502 | 1.6488 | 1.6468 | 0.8251 | 0.8244 | 0.8234 | |
0.1 | 1.6542 | 1.6537 | 1.6512 | 0.8271 | 0.8268 | 0.8256 | |
0.1 | 0.9 | 3.4764 | 3.4523 | 3.4363 | 1.7382 | 1.7261 | 1.7183 |
0.5 | 3.5181 | 3.5045 | 3.4840 | 1.7592 | 1.7522 | 1.7422 | |
0.1 | 3.5578 | 3.5529 | 3.5271 | 1.7789 | 1.7763 | 1.7638 | |
0.01 | 0.9 | 21.7839 | 21.5403 | 21.3807 | 10.8919 | 10.7700 | 10.6921 |
0.5 | 22.1593 | 21.9907 | 21.7937 | 11.0796 | 10.9953 | 10.8988 | |
0.1 | 22.4352 | 22.3172 | 22.0889 | 11.2176 | 11.1586 | 11.0465 |
Nr | Nφ | Nθ | ¯r∗q∗t±δ |
12 | 28 | 28 | 7.342084±7.99×10−4 |
20 | 7.341663±8.49×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341664±8.56×10−5 | ||
44 | 7.341664±8.57×10−5 | ||
28 | 12 | 28 | 7.382330±1.71×10−3 |
20 | 7.351596±2.74×10−4 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.336948±3.89×10−5 | ||
44 | 7.334239±2.21×10−5 | ||
28 | 28 | 12 | 7.341663±8.53×10−5 |
20 | 7.341663±8.53×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341663±8.53×10−5 | ||
44 | 7.341663±8.53×10−5 |
coefficients | scattering phase functions | |||||
F1 | F2 | F3 | B1 | B2 | B3 | |
A0 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 | 1.00000 |
A1 | 2.53602 | 2.00917 | 1.00000 | −0.56524 | −1.20000 | −1.00000 |
A2 | 3.56549 | 1.56339 | 0.29783 | 0.50000 | ||
A3 | 3.97976 | 0.67407 | 0.08571 | |||
A4 | 4.00292 | 0.22215 | 0.01003 | |||
A5 | 3.66401 | 0.04725 | 0.00063 | |||
A6 | 3.01601 | 0.00671 | ||||
A7 | 2.23304 | 0.00068 | ||||
A8 | 1.30251 | 0.00005 | ||||
A9 | 0.53463 | |||||
A10 | 0.20136 | |||||
A11 | 0.05480 | |||||
A12 | 0.01099 |
Ncr | ω | qt(r∗in) | qt(r∗out) | ||||
DOM [10] | FCT [10] | SCM | DOM [10] | FCT [10] | SCM | ||
1 | 0.9 | 1.6460 | 1.6436 | 1.6421 | 0.8230 | 0.8218 | 0.8210 |
0.5 | 1.6502 | 1.6488 | 1.6468 | 0.8251 | 0.8244 | 0.8234 | |
0.1 | 1.6542 | 1.6537 | 1.6512 | 0.8271 | 0.8268 | 0.8256 | |
0.1 | 0.9 | 3.4764 | 3.4523 | 3.4363 | 1.7382 | 1.7261 | 1.7183 |
0.5 | 3.5181 | 3.5045 | 3.4840 | 1.7592 | 1.7522 | 1.7422 | |
0.1 | 3.5578 | 3.5529 | 3.5271 | 1.7789 | 1.7763 | 1.7638 | |
0.01 | 0.9 | 21.7839 | 21.5403 | 21.3807 | 10.8919 | 10.7700 | 10.6921 |
0.5 | 22.1593 | 21.9907 | 21.7937 | 11.0796 | 10.9953 | 10.8988 | |
0.1 | 22.4352 | 22.3172 | 22.0889 | 11.2176 | 11.1586 | 11.0465 |
Nr | Nφ | Nθ | ¯r∗q∗t±δ |
12 | 28 | 28 | 7.342084±7.99×10−4 |
20 | 7.341663±8.49×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341664±8.56×10−5 | ||
44 | 7.341664±8.57×10−5 | ||
28 | 12 | 28 | 7.382330±1.71×10−3 |
20 | 7.351596±2.74×10−4 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.336948±3.89×10−5 | ||
44 | 7.334239±2.21×10−5 | ||
28 | 28 | 12 | 7.341663±8.53×10−5 |
20 | 7.341663±8.53×10−5 | ||
28 | 7.341663±8.53×10−5 | ||
36 | 7.341663±8.53×10−5 | ||
44 | 7.341663±8.53×10−5 |