logo

Nusselt number correlation for turbulent heat transfer of helium–xenon gas mixtures

NUCLEAR ENERGY SCIENCE AND ENGINEERING

Nusselt number correlation for turbulent heat transfer of helium–xenon gas mixtures

Biao Zhou
Yu Ji
Jun Sun
Yu-Liang Sun
Nuclear Science and TechniquesVol.32, No.11Article number 128Published in print 01 Nov 2021Available online 20 Nov 2021
42800

A gas-cooled nuclear reactor combined with a Brayton cycle shows promise as a technology for high-power space nuclear power systems. Generally, a helium–xenon gas mixture with a molecular weight of 14.5–40.0 g/mol is adopted as the working fluid to reduce the mass and volume of the turbomachinery. The Prandtl number for helium–xenon mixtures with this recommended mixing ratio may be as low as 0.2. As the convective heat transfer is closely related to the Prandtl number, different heat transfer correlations are often needed for fluids with various Prandtl numbers. Previous studies have established heat transfer correlations for fluids with medium-high Prandtl numbers (such as air and water) and extremely low-Prandtl fluids (such as liquid metals); however, these correlations cannot be directly recommended for such helium–xenon mixtures without verification. This study initially assessed the applicability of existing Nusselt number correlations, finding that the selected correlations are unsuitable for helium–xenon mixtures. To establish a more general heat transfer correlation, a theoretical derivation was conducted using the turbulent boundary layer theory. Numerical simulations of turbulent heat transfer for helium–xenon mixtures were carried out using Ansys Fluent. Based on simulated results, the parameters in the derived heat transfer correlation are determined. It is found that calculations using the new correlation were in good agreement with the experimental data, verifying its applicability to the turbulent heat transfer for helium–xenon mixtures. The effect of variable gas properties on turbulent heat transfer was also analyzed, and a modified heat transfer correlation with the temperature ratio was established. Based on the working conditions adopted in this study, the numerical error of the property-variable heat transfer correlation was almost within 10%.

Gas-cooled nuclear reactorSpace nuclear powerHelium–xenon mixturesConvective heat transferNusselt number
1

Introduction

Deep-space exploration, a landing on Mars, and the construction of planetary bases will become established activities in the near future. The development of advanced and reliable space power will become hugely important in order to achieve these ambitious milestones. Compared to typical power sources used in space, such as chemical fuel cells or solar photovoltaic arrays, the space nuclear reactor power is characterized by wider power coverage and longer operating duration; as such, it is regarded as a promising technology for human space exploration[1-3]. Among existing technology roadmaps of space nuclear reactors, the gas-cooled nuclear reactor combined with a closed Brayton cycle is perceived as the most rational scheme for the system power exceeding 200 kW[4-6]; this is because it can achieve a higher energy conversion efficiency and a smaller specific mass. Typically, pure helium is adopted as the working fluid for the terrestrial Brayton cycle because of its excellent thermal and transport properties[7]. However, its small molecular weight inevitably results in a larger aerodynamic loading of the turbomachinery, which in turn increases the mass and volume of the system[8]. To reduce the aerodynamic loads and obtain an acceptable heat transfer capacity, a helium–xenon mixture with a molecular weight of 14.5–40.0 g/mol is generally considered as the working fluid for Brayton cycles[9-13]. The Prandtl number (Pr) of helium–xenon mixtures in this molecular weight range is approximately 0.21–0.30, which is significantly lower than that of conventional air or water (Pr > 0.70). For such low-Pr fluids, the thermal boundary layer is thicker, and the similarity of velocity and temperature will be disrupted. Hence, the convective heat transfer characteristics of low-Pr helium–xenon mixtures differ from those of fluids with medium-high Prandtl numbers. Previous studies have also shown that the Dittus-Bolter (DB)[14] and Colburn equations[15] obtained from experimental data of those conventional fluids overpredicted the Nusselt number of low-Pr gas mixtures[16, 17]. Therefore, to better describe the convective heat transfer of helium–xenon mixtures, it is necessary to establish a suitable heat transfer correlation.

The particularity of flow and heat transfer for low-Pr fluids has been widely concerned[18-20]. However, most studies on this issue have mainly focused on liquid metals (Pr ≅ 0.01), where many heat transfer correlations for liquid metals were proposed[21-23]. Lyon et al.[24] first derived a semi-empirical heat transfer correlation for a boundary condition with constant heat flux. Subsequently, Azer et al.[25] and Sleicher et al.[26] proposed their heat transfer correlations for liquid metals using theoretical hypotheses. Based on experimental data for liquid metals, Stomquist et al.[27], Skupinski et al.[28], and Churchill[29] presented their empirical heat transfer correlations. The Pr of liquid metals was found to be one or two orders of magnitude lower than that of the helium–xenon mixture; importantly, the flow and heat transfer characteristics of these two fluids significantly differed. To better illustrate the difference in convective heat transfer, Fig. 1 presents the temperature profiles of the air, helium–xenon mixture, and liquid mercury in a heated tube[30]. The temperature profile of the liquid metal was flat; however, the temperature distribution of the helium–xenon mixture was similar to that of air and relatively steep near the wall. These results show that the thermal resistance of liquid metal is evenly distributed on the radial section, indicating that heat transfer was dominated by molecular heat conduction. The thermal resistance of the helium–xenon mixture was concentrated in the limited near-wall region, and turbulent heat diffusion continued to be dominant. Due to the difference in flow and heat transfer, the heat transfer correlations from liquid metals could not directly be used to calculate flow and heat transfer for helium–xenon mixtures without verification. As such, there is a need to conduct in-depth research on the heat transfer correlation for helium–xenon mixtures within the Pr range.

Fig. 1
(Color online) Temperature profiles of different fluids
pic

This study evaluates the applicability of existing Nusselt number correlations to helium–xenon mixtures. Subsequently, a theoretical derivation of the Nusselt number correlation was conducted to establish a more general heat transfer correlation. Numerical simulations of the turbulent heat transfer for helium–xenon mixtures were also carried out; the simulated results were used to determine related parameters in the derived correlation. Finally, the effects of varying gas properties on turbulent heat transfer were considered, and a new property-variable heat transfer correlation was presented.

2

Application of existing correlations

Table 1 shows that this study selected common correlations and explored their applicability to helium–xenon mixtures. Among these correlations, DB and Colburn’s correlations are empirical correlations obtained from experimental data relating to medium-high Pr fluids. Churchill’s correlation was determined by numerically fitting experimental data for different fluids with a wide range of Prandtl numbers. The uncertainty associated with the above three correlations is provided in Table 1. The Lyon’s and Stomquist’s equations were frequently used for comparison in many studies on the convective heat transfer of liquid metals; however, clear uncertainties have rarely been provided for these equations in published literature. Additionally, it should be noted that the applicable range of most correlations in Table 1 was intentionally extrapolated to evaluate their suitability for helium–xenon mixtures.

Table 1
Existing correlations of Nusselt number
Name Correlations Range of application Max deviation
Dittus-Bolter et al. [14] Nu=0.023Re0.8Pr0.4 0.7≤Pr≤120 ±20%
Colburn [15] Nu=0.023Re0.8Pr1/3 0.5≤Pr≤100 ±25%
Churchill [29] Nu=6.3+0.079RefPr(1+Pr0.8)5/61/f=2.21ln(Re/7) 0.001≤Pr≤200 ±30%
Stomquist et al. [27] Nu=3.6+0.025(PrRe)0.8 Pr<0.1 -
Lyon [24] Nu=7.0+0.025(PrRe/2)0.8 Pr<0.1 -
Show more

Figure 2 shows that calculations were compared with experimental data[17]; the correlations from DB and Colburn had clearly overpredicted the Nusselt number of helium–xenon mixtures, while other calculations were underestimates compared with experimental data. Although predictions by Churchill’s correlation were relatively closer to the experimental data when the Reynolds number (Re) equaled 84 000, it still exceeded the error limit of 10% at Re = 34 000. This indicates that the selected correlations from liquid metals are unsuitable for helium–xenon mixtures.

Fig. 2
(Color online) Predictions of Nusselt number by various correlations
pic
3

Theoretical derivation

To propose a more universal Nusselt number correlation for helium–xenon mixtures, a theoretical derivation based on the turbulent boundary layer theory and eddy viscosity hypothesis was conducted. Obtaining the temperature profile was key to deriving the theoretical heat transfer correlation [31]. First, the derivation processes for the general logarithmic law of temperature was briefly reviewed. Ignoring the varying properties of gas mixtures, the wall shear stress and wall heat flux in the fully developed region may be expressed as follows[32]:

τw=ρ(ν+εm)u¯y, (1) qw=ρCp(α+εh)d(TwT¯)dy. (2)

A two-layer model was used to simplify the boundary layer[33], in which the flow boundary layer was divided into linear and logarithmic regions. Similarly, the thermal boundary layer was divided into molecular conduction and turbulent conduction regions. The linear and molecular conduction regions were dominated by molecular viscosity (ν) and molecular heat diffusivity (α), respectively, neglecting the effect of turbulent momentum diffusivity (εm) and turbulent heat diffusivity (εh). Conversely, turbulent diffusion predominates in the logarithmic and turbulent conduction regions, thus ignoring molecular viscosity and molecular heat diffusivity, respectively. Using the mixing-length theory[34], the profile of dimensionless velocity may be derived by integration as follows:

{U+=y+, y+ycirt+U+=ln(y+/ycirt+)κ+ycirt+, y+>ycirt+ (3)

whereycrit+ denotes the demarcation point between the linear and logarithmic regions, and κ denotes the Karman constant. Based on the definition of turbulent Prandtl number (Prt), Eq. (2) may be normalized, and the profile of dimensionless temperature can be acquired[35], as follows:

{T+=Pry+, y+ycirt,T+T+=Prtln(y+/ycirt,T+)κ+Prycirt,T+, y+>ycirt,T+ (4)

where ycirt,T+ is the demarcation point between the molecular conduction and turbulent conduction regions, the above derivation is reported in several studies[36, 37]. Many empirical logarithmic laws of temperature have been proposed using regression analysis of experimental data. However, these empirical correlations lack a theoretical basis, and as such, their scope of application is limited. There is also an absence of relevant radial experimental data for helium–xenon gas mixtures in the current literature.

Therefore, in present investigation, the derivation of Eq. (4) was sequentially conducted to obtain the dimensionless mass-average temperature (Tb+) and the theoretical Nusselt number correlation. To facilitate this derivation, it was assumed that the thickness of the linear and molecular conduction regions was small and therefore negligible. One study reported that the demarcation point of the flow boundary layer was almost fixed[32]; that is, ycrit+=11, where y+ was defined as per Eq. (5):

y+=uτyν=τw/ρνy. (5)

In which the wall shear (τw) may be expressed as follows:

τw=fρub28. (6)

By adopting the Blasius formula to calculate the friction factor (f), the expression of y+ can be simplified to Eq. (7):

y+=ubf/8νy=0.2Re0.125ubνy. (7)

In which Re=ubνbD. The bulk viscosity (νb) can be approximated to local viscosity (ν), and as there is no significant difference in magnitude, then:

y+=0.2Re0.875yD. (8)

Using this, the ratio of the height of the demarcation point (r) from the wall to the diameter of the tube may be estimated. For ycrit+=11, r/D = 0.035 when Re = 10 000 and r/D = 0.0086 when Re = 50 000; this indicates that the linear region in the tube is a very thin layer. As such, it is reasonable to dismiss the thickness of the linear region when carrying out integration derivation. Additionally, Zhao et al.[38] proposed a correlation to describe the relationship between ycrit,T+  and ycrit+ as:

ycrit,T+ =PrtPrycrit+. (9)

Although this method may not be optimal for helium–xenon mixtures with 0.21 ≤ Pr ≤ 0.30 (as per this study), it may be used as a rough estimate. When Prt = 0.9, then ycrit,T+  is 33.0 < ycrit,T+  < 47.1; when combined with Eqs. (8) and (9), it is possible to calculate the ratio of the height of the demarcation point for the thermal boundary layer (rt) from the wall to the diameter of the tube. Specifically, rt/D = 0.052~0.075 when Re = 10 000, and rt/D = 0.013~0.018 when Re = 50 000, which indicates that the molecular heat conduction region takes up a small proportion of area in the tube; as such, it is reasonable to ignore the thickness of the molecular heat conduction region.

At the wall surface (i.e., y = 0), the axial velocity is reduced to zero. At the center of tube (i.e., y = r0), the axial velocity is the highest and can be described as follows:

Umax+=1κlnr0++ycrit+1κlnycrit+. (10)

By integrating the area under Eq. (3), the expression of dimensionless area-average velocity, Ua+, can be obtained:

Ua+=0r0U+2πrdrπr02. (11)

By subtracting Eq. (10) from Eq. (11), the following relationship can be obtained:

Ua+Umax+=2κr020r0[rln(r0r)rlnr0]dr=-32κ. (12)

Similarly, expressions for dimensionless temperature at the center (Tmax+) can be derived as follows:

Tmax+=Prtκlnr0+ycrit, T++Pr ycrit, T+. (13)

Furthermore, the relationship between Tmax+ and the dimensionless area-average temperature, Ta+, can be obtained using Eq. (14):

Ta+Tmax+=3Prt2κ. (14)

Combining Eqs. (12)–(14), the expression for Ta+ can be derived using Eq. (15):

Ta+=PrtUa+Prtκlnycrit,T++Prycrit,T++Prtκlnycrit+Prtlnycrit+. (15)

However, to derive the theoretical correlation of the Nusselt number, the Tb+ still needs to be determined. Several previous studies[39] have assumed that Tmax+ as Tb+, because differences between the two parameters were <10% for turbulent heat transfer as per Kays et al.[32]. In the present study, it was assumed that Tb+=Ta+ to further reduce deviation. Specifically, a correlation to calculate Tb+ was proposed based on theoretical derivation[39], as:

Tb+Tmax+=3Prt2κ2f83Prt2κ. (16)

Combine Eqs. (14) and (16), then,

Ta+Tb+=13PrtNu2κ2RePr. (17)

Based on a working condition with Re = 30 000 and Pr = 0.30, where the Churchill correlation was used to calculate the Nusselt number, the second term on the right of Eq. (17) equals 0.030 and decreases with an increase in Re; as such, it was reasonable to assume Tb+=Ta+.

Using the definition of Nusselt number and T+, the Nusselt number correlation may be transformed as follows:

{Nu=qwTwTbDλT+=ρCpuτ(TwT)qwf=8τw/ρub2 Nu RePrf/8Tb+, (18)

where f is the friction factor. Combining Eqs. (15) and (18), the following theoretical expression for the Nusselt number was obtained:

Nu=RePrf/8Prt8/fPrtκlnycrit,T++Prycrit,T++Prtκlnycrit+Prtycrit+. (19)

Eq. (19) shows that the appropriate f, ycrit+, ycirt,T+, and average Prt still need to be determined in order to obtain the new Nusselt number correlation. To do this, numerical simulations of turbulent heat transfer for helium–xenon mixtures were carried out, and these parameters were determined using the simulation results.

4

Numerical simulations

4.1
Thermophysical and transport properties

Accurate properties of helium–xenon mixtures are essential to numerical simulations. Typically, existing designs of space gas-cooled nuclear reactor systems have an operating pressure of 2 MPa and a temperature range of approximately 400–1300 K. Under these conditions, the theory of diluted gases based on the Chapman-Enskog approach was not found to be applicable. Tournier et al.[40] reviewed the experimental data for the properties of dense noble gas mixtures across a pressure range of 0.1–20 MPa and temperature up to 1400 K, summarizing semi-empirical property correlations. Based on the study of Tournier et al.[40], a property-calculation code for helium–xenon mixtures was developed, and its accuracy was validated in a previous study[30]. Figure 3 presents the calculation of variations in gas mixture properties with temperature.

Fig. 3
(Color online) Variations in the properties of helium–xenon mixtures with temperature at P = 2 MPa
pic

These results indicate that variations in thermal conductivity (λ) and dynamic viscosity (μ) are proportional to temperature, while the specific heat capacity (Cp) and density (ρ) decrease as temperature increases. As the value of the specific heat capacity varies very slightly with temperature, it can be set as a constant for a particular gas mixture in numerical simulations. Other property correlations were implemented in Ansys Fluent using user-defined functions.

4.2
Numerical approach

The turbulent heat transfer of the helium–xenon mixtures in a tube was simulated using Ansys Fluent. The calculation model refers to the test section of the experiment in Taylor et al.[17], as shown in Fig. 4.

Fig. 4
(Color online) Test section of experiment for helium–xenon mixtures by Taylor et al. (a) Schematic diagram of experimental device; (b) Heated tube of test section
pic

The experimental loop was a closed circuit driven by a gas booster pump; the flow and heat transfer of the gas mixture was investigated in the test section. Figure 4(b) shows the geometric parameters of the test section, which consisted of a straight tube with an inner diameter (D) of 5.87 mm. The total length (L) was 680.92 mm, of which the first part (L1) was adiabatic, while the remainder was heated by a resistance wire to obtain uniform heat flux. A two-dimensional calculation domain was adopted, as the gravitational effect was neglected.

With respect to the simulated boundary conditions, the inlet was set to ‘mass-flow-inlet’, the outlet was set to ‘pressure-out’, and the wall surface adopted a non-slip boundary. The simulation parameters were set according to experimental conditions, listed in Table 2. A previous paper[41] evaluated the applicability of different turbulence models to helium–xenon mixtures; based on the findings from this paper, the SST κ–ω turbulence model was selected for this study. To improve the accuracy of the numerical simulation for turbulent heat transfer of low-Pr helium-xenon mixtures, a local Prt model (shown in Eqs. (20) and (21)), was proposed in Zhou et al.[30]. Based on comparison with more experimental data, in this study the first constant in the expression of Prt was adjusted slightly from the original 0.80 to 0.86.

Table 2
Simulation parameters of helium–xenon mixtures
M (g/mol) No. Pr Rei q (W/m2) Tw/Tb (max) Pout (Pa)
14.5 Run711 0.30 40869 45712 1.15 993150
  Run712 0.30 40554 97522 1.32 982023
  Run713 0.30 35374 175055 1.61 828120
  Run714 0.30 38962 242158 1.84 810084
  Run715 0.30 34042 296622 2.04 806581
28.3 Run707 0.23 53390 157094 1.67 928019
  Run709 0.23 48411 240978 2.06 835497
  Run710 0.23 49001 33809 1.16 779391
40.0 Run696 0.21 36183 136770 1.99 563474
  Run698 0.21 34280 78313 1.59 502892
  Run700 0.21 34865 18151 1.15 465129
  Run704 0.21 59851 31593 1.17 771105
  Run705 0.21 61772 159072 1.77 990763
83.8 Run687 0.25 85802 43637 1.42 471502
  Run689 0.25 87373 96326 1.89 480141
  Run693 0.25 33513 37260 1.62 262358
Show more
Prt=1/{12Prt,+CPet1Prt,(CPet)2×[1exp(1CPetPrt,)]} (20) Prt=0.86+30Relocal0.888Pr, Relocal=uDν (21)

Grid independence verification was investigated to ensure the validity of the numerical calculation. Four mesh files were tested: mesh1: 20*3000, mesh2: 30*4000, mesh3: 40*6000, and mesh4: 60*8000. The bulk temperature (Tb) of helium–xenon mixtures and the wall temperature (Tw) at x = 0.45, 0.55, 0.60, and 0.65 m were calculated. Figure 5 shows that when the grid number exceeded 120 000, the difference in the bulk and wall temperatures as calculated by various meshes was negligible, and the maximum error was <0.8%. Based on the calculation accuracy and efficiency, mesh3 was adopted in this study.

Fig. 5
(Color online) Simulated bulk temperature (Tb) and wall temperature (Tw) with various meshes
pic

Figure 6 presents the grid of mesh3; the grid refinement was applied in the near-wall region, where the spacing of the first grid was 0.004 mm, spacing ratio was 1.08, and maximal aspect ratio was 28.4. Using the Run696 as an example, the corresponding y+ of the first grid was y+ = 0.22. In fact, the wall-resolved method was enabled for the SST κ-ω model when the mesh in the near-wall region was sufficient. At that time, the y+ of the first grid was generally required to be <1.0; as such, mesh3 was able to meet the calculation requirements. In order to further verify the simulation models adopted in the present study, additional experimental data on turbulent heat transfer for helium-xenon gas mixtures were introduced for comparison. Figure 7 shows that the calculated deviation of the wall temperature under various working conditions was essentially 5%. Therefore, the modified turbulence model may be used to study the turbulent heat transfer characteristics of helium–xenon mixtures.

Fig. 6
(Color online) Grid of mesh3
pic
Fig. 7
(Color online) Comparison of simulated wall temperature (Tw) by the modified turbulent model and experimental data
pic
4.3
Determination of parameters

There is a need to analyze the momentum diffusion characteristics of helium-xenon mixtures to determine the ycrit+ of the flow boundary layer. At the position of ycrit+, the turbulent momentum diffusion should be equal to the molecular momentum diffusion. Figure 8 shows that the molecular viscosity (ν) and turbulent momentum diffusivity (εm) of helium–xenon mixtures with two different mixing ratios in the fully developed region (x = 600 mm) were calculated. The results show that the intersections of the curves (ν = εm) were always close to y+ = 11; as such, the demarcation points of the flow boundary layer for the helium–xenon mixtures were determined to be ycrit+=11. Notably, the demarcation points of the flow boundary layer for water and air were also at approximately y+ = 11 in Schlichting et al.[33]. These results indicate that Pr has little effect on the turbulent momentum diffusion of fluids.

Fig. 8
(Color online) Variations in molecular viscosity (ν) and turbulent momentum diffusivity (εm) with y+ for helium–xenon mixtures
pic

To investigate the appropriate equation for the friction factor (f), four typical equations from conventional fluids[42-44] (see Table 3), were selected for testing. One study[16] suggested that the experimental friction factor of noble binary gases may be calculated using Eqs. (22)–(24). The term f is the overall average friction factor, calculated using the distance from the start of the heated section to the second pressure tap (‘pt2’ in Fig. 4):

Table 3
Four equations used to determine the friction factor
Name Formula Range of application
Blasius [42] f=0.316/Re0.25 5.0×103Re105
Colebrook [43] 1/f=1.8lg(Re/6.9) 5.0×103Re5.0×107
Drewr [32] f=0.0056+0.5×Re0.32 3.0×103Re3.0×106
Taitel and Dukler [44] f=0.184/Re0.2 Re3.0×103
Show more
ΔPf=ΔPΔPm, (22) ΔPm=G2Rggc(Tb2P2Tb1P1), (23) f=2ΔPfρ¯DlG2, (24)

where ΔP denotes the total pressure drop; ΔPm denotes the acceleration pressure drop; and ΔPf denotes the friction pressure drop. Figure 9 shows that the calculated values were relatively close to each other and calculation errors were generally within 5%. Based on the good agreement between the experimental data and the concise expression, the Blasius equation was used to calculate the friction factor for helium–xenon mixtures in Eq. (19).

Fig. 9
(Color online) Comparison of calculations using different friction factor equations with experimental data
pic

Notably, the Karman constant (κ) exists in the velocity distribution (see Eq. (3)); as such, the simulated velocity distribution may be numerically fitted to determine κ for helium–xenon mixtures. Figure 10(a) shows that the fitting logarithmic law of velocity for Pr = 0.30 was initially obtained by setting κ = 0.42. Then, a comparison between the calculated logarithmic law with κ = 0.42 and the simulation results of Pr = 0.21 was conducted as shown in Fig. 10(b). The calculated logarithmic law with κ = 0.38 was still in good agreement with the simulation results; as such, κ = 0.42 was adopted in Eq. (19).

Fig. 10
(Color online) Fitted velocity profiles of helium–xenon mixtures 
pic

To obtain the appropriate demarcation point of the thermal boundary layer for helium–xenon mixtures, variations in the molecular heat diffusivity (α) and the turbulent heat diffusivity (εh) were calculated. In this study, the criterion for defining ycrit,T+ is that α = εh. Figure 11 shows that it is evident that the demarcation points were not fixed for different Prs; the smaller the Pr, the farther the demarcation point from the wall. Therefore, it was necessary to establish a model of turbulent thermal diffusivity for helium–xenon mixtures to fix the value of ycrit,T+.

Fig. 11
(Color online)Variations in molecular heat diffusivity (α) and turbulent heat diffusivity (εh) with the y+ of helium–xenon mixtures
pic

The Reynolds-averaged Navier-Stokes equations (RANS) model was developed based on the Boussinesq eddy viscosity hypothesis[45], in which eddy viscosity was calculated using mixing-length theory; as such, εm/ν was a function related only to y+. By analyzing the Prt model (Eqs. (20) and (21)), the Prt was also found to be a function only related to y+ for a specific helium–xenon mixture, as Pet and Relocal were expressed as functions of y+. Therefore, the following relationship can be obtained:

εhν=εmν1Prt =f(y+) (25)

It is feasible to obtain Eq. (25) by numerical fitting. Since the thermal boundary layer of low-Pr helium–xenon mixture was thicker than the flow boundary layer, the fitting range of y+ should ensure that y+ > 11. First, the expression was initially determined by fitting the simulation results of Pr = 0.30, as shown in Fig. 12(a); this expression could be described by Eq. (26). Thereafter, the cases with the Prandtl number to 0.23 and 0.21 were simulated. The results are shown in Fig. 12(b) and 12(c), where it is shown that the calculation of Eq. (26) was in good agreement with simulated data points. Therefore, for 0.21 ≤ Pr ≤ 0.30 and y+ < 50, Eq. (26) was able to reasonably describe the turbulent heat diffusivity of helium-xenon mixtures.

Fig. 12
Numerical fitting of the expression for εh/ν
pic
εhν = 0.0112(y+)1.82 (26)

As α = ν/Pr, the ycrit,T+ may be determined as follows:

ycrit,T+ = 11.83Pr-0.55. (27)

The influence of the Pr on turbulent heat transfer may be quantitatively analyzed using the numerical model of the turbulent heat diffusivity. For air and water, the relationship between turbulent heat diffusivity and y+ was εh/ν(y+)3 [36, 37]. Compared with Eq. (26), it can be seen that the lower Prandtl number will lead to a smaller turbulent thermal diffusivity at the same radial position, resulting in weak convective heat-transfer intensity. This conclusion may be used to understand the results in Fig. 2. The Prandtl number of liquid metals was lower than that of the binary gas mixtures. As the turbulent heat diffusion was weaker for the same Re, the Nusselt number correlations obtained from liquid metals underestimated the experimental data of binary gas mixtures. Conversely, the Pr of conventional fluids was higher than that of helium–xenon mixtures, resulting in more intense turbulent heat diffusion. Therefore, predictions from the DB and Colburn correlations represented overestimations of experimental data.

For working fluids with Pr≫1.0, the Prt tends to be constant of 0.85. Differently, for low-Pr fluids it tends to be larger[46]. However, it is difficult to obtain the section average Prt directly from the local function expressions. In this study, a constant value was used to roughly represent the section average Prt. To determine this parameter, the radial Prt distribution of Run711 was first calculated using the Zhou Prt model (Eqs. (20) and (21)). As shown in Fig. 13, by comparing with the prediction of the Zhou model, a Prt of 0.9 may be a candidate value; as such, this value was initially adopted as the section average Prt in Eq. (19) for the helium–xenon mixtures.

Fig. 13
Radial predictions of Prt at x = 0.60 m for Run711
pic

Once all parameters were determined, the semi-theoretical Nusselt number correlation was quantified using Eq. (28):

Nu=0.20PrRe0.8754.53Re0.125+11.83Pr0.45+1.18lnPr10.05. (28)
5

Verification and analysis

Combining Eqs. (4) and (27), the radial temperature distribution of the corresponding helium–xenon mixtures with 0.21 ≤ Pr ≤ 0.30, was determined as shown in Eq. (29):

{T+=Pry+, y+11.83Pr0.55T+=3.02(lny++0.55lnPr2.47)+11.83Pr0.45, y+>11.83Pr0.55. (29)

To verify the correctness of the derivation, the results using Eq. (29) where the Pr was 0.20, were compared with predictions of other temperature distribution correlations (shown in Table 4). Additionally, the direct numerical simulation results of Kawamura et al.[47] were introduced as benchmarks. As shown in Fig. 14, the calculations using Eq. (29) showed better agreement with the DNS results than those of other correlations, preliminarily proving the validity of Eq. (29) and the correctness of this derivation.

Table 4
Existing correlations for temperature distribution
Authors Correlations Scope of application
Duponcheel et al. [19] T+=5ln(1+0.2Pry+) 0.01 < Pr < 0.025
Kader [31] T+=2.12lny++(3.85Pr1/31.3)22.12lnPr 0.025 < Pr < 170
Gowen and Smith [35] T+=2.5lny++5ln[(5Pr+1)/30]+8.55+5Pr 0.7 < Pr < 14.3
Kays et al. [32] T+=2.2lny++13.39Pr2/35.66 0.7 < Pr < 5.9
Show more
Fig. 14
Comparison of prediction by Eq. (29) with other temperature correlations
pic

To test the accuracy of the new Nusselt number correlation, predictions using Eq. (28) were compared with experimental data, as shown in Fig. 15. The results show that the calculations are in good agreement with the experimental data of gas mixtures under both working conditions. This further verifies the accuracy of the derivation and its applicability to the turbulent heat transfer for helium–xenon mixtures where 0.21 ≤ Pr ≤ 0.30.

Fig. 15
Comparison of predictions using Eq. (28) with experimental data
pic

A sensitivity analysis of each parameter was carried out to evaluate the influence of the parameters that were determined on the derived Nusselt number correlation. As the empirical range of the Karman constant was 0.43 ≥ κ ≥ 0.35, five different constants were tested. Figure 16 shows that even when the value of κ fluctuates in the entire empirical range, there is only a very small difference between the calculated Nusselt numbers. This proves that the calculation of the derived Nusselt number correlation is insensitive to κ. Thereafter, three additional constant Prt numbers were used for sensitivity analysis. Figure 17 demonstrates that the derived Nusselt number correlations vary with different Prt numbers. The predictions where Prt was 0.9 were in better agreement with experimental data, in which a fluctuation of 0.9±0.1 appeared to be valid. Additionally, an error propagation analysis of the demarcation point of the flow boundary layer (ycrit+) was also conducted. Figure 18 shows that the calculations where ycrit+=11 were closer to experimental data, in which a fluctuation of 11±1.0 may also be considered reasonable.

Fig. 16
Sensitivity analysis on the Karman constant for the derived Nusselt number correlation
pic
Fig. 17
Sensitivity analysis on the Prt number for derived Nusselt number correlation
pic
Fig. 18
(Color online) Sensitivity analysis on the ycrit+ for derived Nusselt number correlation
pic

Importantly, helium–xenon mixtures flow through a strongly heated tube in space gas-cooled reactors; under these conditions, the wall temperature is generally higher than the bulk temperature. The physical properties of helium–xenon mixtures at different radial positions will be significantly different. Therefore, the Nusselt number correlation deduced from constant physical properties may need to be modified. To examine the influence of variable properties on turbulent heat transfer, the Nusselt number calculated using Eq. (28) (Nucp) was compared with measured experimental data (Ref. [13]) (NuExp) at various temperature ratios, as shown in Fig. 19.

Fig. 19
Comparison of calculations using Eq. (28) with measured experimental data
pic

The ranges of the tested conditions are 18000 < Re < 60000, 0.21 ≤ Pr ≤ 0.30, and Tw/Tb < 2. Figure 19 shows that the calculation deviation becomes obvious as Tw/Tb increases. When Tw/Tb = 1.8, the numerical error exceeded 30%. A seemingly exponential relationship between NuExp /NuCp and Tw/Tb was observed when NuExp/NuCp were presented in logarithmic coordinates. Moreover, as Tw/Tb approaches 1.0, the constant-property correlation would become valid, and NuExp/NuCp would approach 1.0. Therefore, by numerically fitting data points, the modified, property-variable Nusselt number correlation for the helium–xenon mixture was obtained as follows:

Nuvp=0.20PrRe0.8754.53Re0.125+11.83Pr0.45+1.18lnPr10.05(TwTb)0.63. (30)

The results in Fig. 19 show that the error in predictions using Eq. (30) was almost within 5% and was consistently within 10%. Notably, the modified Nusselt correlation (Eq. (30)) was only determined using existing and limited experimental data; to fully verify its applicability, additional experimental data are required.

6

Conclusion

This study evaluates existing Nusselt number correlations and derives a theoretical Nusselt number correlation for helium–xenon mixtures using turbulent boundary layer theory. The turbulent heat transfer in a tube for helium–xenon mixtures was numerically simulated using a modified Prt model. Based on numerical simulations, a new turbulent heat diffusivity model in the thermal boundary layer was developed, and related parameters were determined in the theoretical Nusselt number expression. A preliminary modified Nusselt number correlation was presented based on variations in gas properties. The main conclusions from this study were:

(1) Compared with the extrapolated experimental data of gas mixtures, the selected Nusselt number correlations from conventional fluids and liquid metals were unsuitable for helium–xenon mixtures with 0.21 ≤ Pr ≤ 0.30;

(2) The calculation error of the friction factor for helium–xenon mixtures using the Blasius equation was within 5%, and the demarcation point of U+ for helium–xenon mixtures and conventional fluids was approximately equal;

(3) For 0.21 ≤ Pr ≤ 0.30 and y+ < 40, εh/ν = 0.0112(y+)1.82 can reasonably describe the turbulent heat diffusivity of helium–xenon mixtures. A lower Pr led to a reduced turbulent thermal diffusivity, resulting in weak convective heat-transfer intensity;

(4) Predictions by the derived semi-theoretical Nusselt number correlation were in good agreement with experimental data, verifying the accuracy of the derivation and its applicability to turbulent heat transfer for helium–xenon mixtures.

Future research should focus on obtaining additional experimental research on low-Pr helium–xenon mixtures. This additional experimental data may be used to verify correlations proposed in this study. Furthermore, the new Nusselt number correlation should be used to develop a thermal-hydraulic transient analysis code for helium–xenon cooled reactor systems.

Reference
1 El-Genk M.S., Schriener T.M.,

Long operation life reactor for lunar surface power

. Nucl. Eng. Des. 241, 2339-2352 (2011). doi: 10.1016/j.nucengdes.2011.02.024
Baidu ScholarGoogle Scholar
2 Howe T., Howe S., Miller J.,

Novel deep space nuclear electric propulsion spacecraft

. Nucl. Technol. 207, 866-875 (2020). doi: 10.1080/00295450.2020.1832814
Baidu ScholarGoogle Scholar
3 Wang C., Liu T., Tang S. et al.,

Thermal-hydraulic analysis of space nuclear reactor TOPAZ-II with modified RELAP5

. Nucl. Sci. Tech. 30, 12 (2019). doi: 10.1007/s41365-018-0537-3
Baidu ScholarGoogle Scholar
4 El-Genk M.S., Tournier J.P.,

Noble gas binary mixtures for closed-Brayton-cycle space reactor power systems

. J. Propul. Power. 23, 863-873 (2007). doi: 10.2514/1.27664
Baidu ScholarGoogle Scholar
5 Tournier J.P., El-Genk M.S.,

Properties of noble gases and binary mixtures for closed Brayton Cycle applications

. Energy Conv. Manag. 49, 469-492 (2008). doi: 10.1016/j.enconman.2007.06.050
Baidu ScholarGoogle Scholar
6 Mason L.S..

A comparison of Brayton and stirling space nuclear power systems for power levels from 1 kilowatt to 10 megawatts

, in Proceedings of American Institute of Physics Conference (2001). doi: 10.1063/1.1358045
Baidu ScholarGoogle Scholar
7 Zhang X., She D., Chen F. et al.,

Temperature analysis of the HTR-10 after full load rejection

. Nucl. Sci. Tech. 30, 11 (2019).
Baidu ScholarGoogle Scholar
8 Gallo B.M., El-Genk M.S.,

Brayton rotating units for space reactor power systems

. Energy Conv. Manag. 50, 2210-2232 (2009). doi: 10.1016/j.enconman.2009.04.035
Baidu ScholarGoogle Scholar
9 El-Genk M.S., Tournier J.P.,

Noble gas binary mixtures for gas-cooled reactor power plants

. Nucl. Eng. Des. 238, 1353-1372 (2008). doi: 10.1016/j.nucengdes.2007.10.021
Baidu ScholarGoogle Scholar
10 Meng T., Cheng K., Zhao F. et al.

Computational flow and heat transfer design and analysis for 1/12 gas-cooled space nuclear reactor

. Ann. Nucl. Energy. 135, 106986, (2020). doi: 10.1016/j.anucene.2019.106986
Baidu ScholarGoogle Scholar
11 Li Z.G., Sun J., Lang M.G. et al.,

Design of a hundred-kilowatt level integrated gas-cooled space nuclear reactor for deep space application

. Nucl. Eng. Des. 361, 110569, (2020). doi: 10.1016/j.nucengdes.2020.110569
Baidu ScholarGoogle Scholar
12 Schriener T.M., El-Genk M.S.,

Effects of decreasing fuel enrichment on the design of the pellet bed reactor (PeBR) for lunar outpost

. Prog. Nucl. Energy. 104, 288-297 (2018). doi: 10.1016/j.pnucene.2017.10.010
Baidu ScholarGoogle Scholar
13 Solonin V.I., Satin A.A., Dunaitsev A.A. et al.,

Hydrodynamic characteristics of coolant tracts in high-temperature gas-cooled nuclear reactor

. Atom. Energy. 124, 302-308 (2018). doi: 10.1007/s10512-018-0414-5
Baidu ScholarGoogle Scholar
14 Dittus F.W., Boelter L.M.K.,

Heat transfer in automobile radiators of the tubular type

. Int. Commun. Heat Mass Transf. 12, 3-22 (1985). doi: 10.1016/0735-1933(85)90003-X
Baidu ScholarGoogle Scholar
15 Colburn A.P.,

A method of correlating forced convection heat-transfer data and a comparison with fluid friction

. Int. J. Heat Mass Transf. 7, 1359-1384 (1964). doi: 10.1016/0017-9310(64)90125-5
Baidu ScholarGoogle Scholar
16 Pickett P.E., Taylor M.F., McEligot D.M.,

Heated turbulent flow of helium-argon mixtures in Tubes

. Int. J. Heat Mass Transf. 22, 705-719 (1979). doi: 10.1016/0017-9310(79)90118-2
Baidu ScholarGoogle Scholar
17 Taylor M.F., Bauer K.E., McEligot D.M.,

Internal forced convection to low-Prandtl-number gas mixtures

. Int. J. Heat Mass Transf. 31, 13-25 (1988). doi: 10.1016/0017-9310(88)90218-9
Baidu ScholarGoogle Scholar
18 Shams A., Santis A.D., Koloszar L.K. et al.,

Status and perspectives of turbulent heat transfer modelling in low-Prandtl number fluids

. Nucl. Eng. Des. 353, UNSP 110220 (2019). doi: 10.1016/j.nucengdes.2019.110220
Baidu ScholarGoogle Scholar
19 Duponcheel M., Bricteux L., Manconi M. et al.

Assessment of RANS and improved near-wall modeling for forced convection at low Prandtl numbers based on LES up to Reτ = 2000

. Int. J. Heat Mass Transf. 75, 470-482 (2014). doi: 10.1016/j.ijheatmasstransfer.2014.03.080
Baidu ScholarGoogle Scholar
20 Vitovsky O.V., Makarov M.S., Nakoryakov V.E. et al.,

Heat transfer in a small diameter tube at high Reynolds numbers

. Int. J. Heat Mass Transf. 109, 997-1003 (2017). doi: 10.1016/j.ijheatmasstransfer.2017.02.041
Baidu ScholarGoogle Scholar
21 Shams A., Santis A.D.,

Towards the accurate prediction of the turbulent flow and heat transfer in low-Prandtl fluids

. Int. J. Heat Mass Transf. 130, 290-303 (2019). doi: 10.1016/j.ijheatmasstransfer.2018.10.096
Baidu ScholarGoogle Scholar
22 Bricteus L., Duponcheel M., Winckelmans G. et al.,

Direct and large eddy simulation of turbulent heat transfer at very low Prandtl number: Application to lead-bismuth flows

. Nucl. Eng. Des. 246, 91-97 (2012). doi: 10.1016/j.nucengdes.2011.07.010
Baidu ScholarGoogle Scholar
23 Zier O., Zimmermann W., Pesch W.,

On low-Prandtl-number convection in an inclined layer of liquid mercury

. J. Fluid Mech. 874, 76-101 (2019). doi: 10.1017/jfm.2019.432
Baidu ScholarGoogle Scholar
24 Lyon R.N.,

Liquid metal heat-transfer coefficients

. Chem. Eng. Prog. 47, 75-79 (1951).
Baidu ScholarGoogle Scholar
25 Azer N.Z., Chao B.T.,

A mechanism of turbulent heat transfer in liquid metals

. Int. J. Heat Mass Transf. 1, 121 (1960). doi: 10.1016/0017-9310(60)90016-8
Baidu ScholarGoogle Scholar
26 Notter R.H., Sleicher C.A.,

Solution to turbulent Graetz problem.3. fully developed and entry region heat-transfer rates

. Chem. Eng. Sci. 27, 2073 (1972). doi: 10.1016/0009-2509(72)87065-9
Baidu ScholarGoogle Scholar
27 Chen F., Huai X.L., Cai J. et al.,

Investigation on the applicability of turbulent-Prandtl-number models for liquid lead-bismuth eutectic

. Nucl. Eng. Des. 257, 128-133 (2013). doi: 10.1016/j.nucengdes.2013.01.005
Baidu ScholarGoogle Scholar
28 Cheng X., Tak N.,

Investigation on turbulent heat transfer to lead-bismuth eutectic flows in circular tubes for nuclear applications

. Nucl. Eng. Des. 236, 385-393 (2006). doi: 10.1016/j.nucengdes.2005.09.006
Baidu ScholarGoogle Scholar
29 Churchill S.W.,

Comprehensive correlating equations for heat, mass and momentum transfer in fully developed flow in smooth tubes

. Ind. Eng. Chem. Fundam. 16, 109-116 (1977). doi: 10.1021/i160061a021
Baidu ScholarGoogle Scholar
30 Zhou B., Ji Y., Sun J. et al.,

Modified turbulent Prandtl number model for helium–xenon gas mixture with low Prandtl number

. Nucl. Eng. Des. 366, 110738, (2020). doi: 10.1016/j.nucengdes.2020.110738
Baidu ScholarGoogle Scholar
31 Kader B.A.,

Temperature and concentration profiles in fully turbulent boundary layers

. Int. J. Heat Mass Transf. 24, 1541-1544 (1981). doi: 10.1016/0017-9310(81)90220-9
Baidu ScholarGoogle Scholar
32 Kays W.M., Crawford M.E., Weigand B., Convective heat and mass transfer, Vol. 4 McGraw-Hill Science (2005).
33 Schlichting H., Gersten K., Krause E., Boundary-layer theory, Springer (2017).
34 Weigand B., Ferguson J.R., Crawford M.E.,

An extended Kays and Crawford turbulent Prandtl number model

. Int. J. Heat Mass Transf. 40, 4191-4196 (1997). doi: 10.1016/S0017-9310(97)00084-7
Baidu ScholarGoogle Scholar
35 Gowen R.A., Smith J.W.,

The effect of the Prandtl number on temperature profiles for heat transfer in turbulent pipe flow Chem

. Eng. Sci. 22, 1701 (1967). doi: 10.1016/0009-2509(67)80205-7
Baidu ScholarGoogle Scholar
36 Aravinth S.,

Prediction of heat and mass transfer for fully developed turbulent fluid flow through tubes

. Int. J. Heat Mass Transf. 43, 1399-1408 (2000). doi: 10.1016/S0017-9310(99)00218-5
Baidu ScholarGoogle Scholar
37 Slaiman Q.J., Abu-Khader M.M., Hasan B.O.,

Prediction of heat transfer coefficient based on eddy diffusivity concept Prediction of Heat Transfer Coefficient Based on Eddy Diffusivity Concept

,” Eng. Res. Des. 85, 455-464 (2007). doi: 10.1205/cherd06002
Baidu ScholarGoogle Scholar
38 Zhao H., Li X., Wu X.,

New friction factor and Nusselt number equations for turbulent convection of liquids with variable properties in circular tubes

. Int. J. Heat Mass Transf. 124, 454-462 (2018). doi: 10.1016/j.ijheatmasstransfer.2018.03.082
Baidu ScholarGoogle Scholar
39 Kader B.A., Yaglom A.M.,

Heat and mass-transfer laws for fully turbulent wall flows

. Int. J. Heat Mass Transf. 15, 2329-3000 (1972). doi: 10.1016/0017-9310(72)90131-7
Baidu ScholarGoogle Scholar
40 Tournier J.P., El-Genk M.S., Gallo B. M.,

Best estimates of binary gas mixtures properties for closed Brayton cycle space applications

, in Proceedings of 4th International Energy Conversion Engineering Conference and Exhibit(2006).
Baidu ScholarGoogle Scholar
41 Zhou B., Ji Y., Sun J. et al.,

CFD simulations of the helium–xenon mixture flow and heat transfer

, in Proceedings of International Conference on Nuclear Engineering (2019).
Baidu ScholarGoogle Scholar
42 Rott N.,

Note on the history of the Reynolds number

. Rev. Fluid Mech. 22, 1-11 (1990). doi: 10.1146/annurev.fl.22.010190.000245
Baidu ScholarGoogle Scholar
43 Colebrook C.F., White C.M.,

Experiments with fluid friction in roughened pipes

, in Proceedings of the Royal Society of London Series A-Mathematical and Physical Sciences(1937).
Baidu ScholarGoogle Scholar
44 Taitel Y., Dukler A.E.,

A model for predicting flow regime transitions in horizontal and near horizontal gas-liquid flow

. AICHE J. 22, 47-55 (1976). doi: 10.1002/aic.690220105
Baidu ScholarGoogle Scholar
45 Grifoll J., Giralt F.,

The near wall mixing length formulation revisited

. Int. J. Heat Mass Transf. 43, 19, 3743-3746 (2000). doi: 10.1016/S0017-9310(00)00009-0
Baidu ScholarGoogle Scholar
46 Kays W.M.,

Turbulent Prandtl number-where are we

, J. Heat Transf.-Trans. ASME. 116, 284-295 (1994). doi: 10.1115/1.2911398
Baidu ScholarGoogle Scholar
47 Kawamura H., Ohsaka K., Abe H. et al.,

DNS of turbulent heat transfer in channel flow with low to medium-high Prandtl number fluid

. Int. J. Heat Mass Transf. 19, 482-491 (1998). doi: 10.1016/S0142-727X(98)10026-7
Baidu ScholarGoogle Scholar