Next Article in Journal
Normalized Model Reference Adaptive Control Applied to High Starting Torque Scalar Control Scheme for Induction Motors
Next Article in Special Issue
Comparative Research on Ventilation Characteristics of Scattering and Sample Room from Chinese Spallation Neutron Source
Previous Article in Journal
Study on the Uniformity of Secondary Air of a 660 MW Ultra-Supercritical CFB Boiler
Previous Article in Special Issue
Optical Model of Thermal Radiation Loading System for Turbine Vane Leading Edge
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Convective Heat Transfer Characteristics of Supercritical Liquid Hydrogen in a U-Type Tube inside a Moderator

1
Department of Mechanical and Power Engineering, Harbin University of Science and Technology, Harbin 150080, China
2
Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
3
China Spallation Neutron Source Science Center, Dongguan 523803, China
4
School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(10), 3605; https://doi.org/10.3390/en15103605
Submission received: 14 April 2022 / Revised: 8 May 2022 / Accepted: 13 May 2022 / Published: 14 May 2022
(This article belongs to the Special Issue Advanced Thermal Analysis and Management Technology)

Abstract

:
The flow and heat transfer characteristics of supercritical fluid in a U-tube have an important influence on the safe operation of a moderator, and the variation of gravity direction is suitable for special working conditions of the moderator. In this study, the three-dimensional turbulence flow and heat transfers of supercritical liquid hydrogen in a U-tube were investigated at an Re number ranging from 16,425 to 54,750 under constant heat flux (q = 80 kW/m2). The total length of the U-tube was 1725 mm, which had an entrance length L/D of 23, with the inner diameter and wall thickness of D × δ = 10 × 2 mm. The finite volume method was adopted, and the grid independence was verified by the grid convergence index (GCI). The calculation results of three turbulence models (SST k-w, RNG k-ε, Standard k-ε) were compared with the corresponding experimental data to obtain the turbulence model with the smallest error. The convective heat transfer characteristics with different values of heat flux (q = 30 kW/m2~100 kW/m2), mass flow (G = 3 g/s~10 g/s), and gravity (gx, gy, gz) were compared. Meanwhile, the heat transfer characteristics of supercritical liquid and conventional liquid hydrogen were compared. The results show that Nu increased from 5 g/s to 10 g/s by 56.6%, and mass flow rate had a greater impact on the variation of Nu; when gravity direction was consistent with the flow direction of liquid hydrogen (gx direction), the Nu number inside the channel was 4.21% and 5.56% higher than that in gy and gz direction, respectively. Supercritical liquid hydrogen has a stronger heat transfer ability than conventional liquid hydrogen, of which the Nu number is 16.7% higher. This study can provide useful guidance for the design of flow and heat transfer of supercritical liquid hydrogen in a U-tube and its application in moderators. Furthermore, it provides reference technical values for thermal safety and thermal management of the target station to ensure its safe and stable operation.

1. Introduction

In recent years, the rapid development of supercritical fluid technology has been widely used in aerospace, electrical engineering, nuclear technology, low temperatures, etc. The Chinese Spallation Neutron Source (CSNS), as the first spallation neutron source in China, has completed the first phase of its construction. In the process of using supercritical hydrogen as a moderator to improve neutron intensity and slow down neutron energy, it was proposed to use a composite moderator in the upcoming second phase project to provide more energy spectrum options for spectrometers [1]. The moderator is made by a combination of supercritical liquid hydrogen tubes and water. A U-tube can be used in a moderator due to its compact structure and high heat transfer efficiency [2]. Under the condition of high thermal density and near the critical temperature, the physical properties of supercritical fluid will change dramatically, and its flow and heat transfer characteristics will be affected by different channel structures. Meanwhile, the thermodynamic process is very complex [3], so it is necessary to study the flow and heat transfer characteristics of supercritical liquid hydrogen in a U-tube [4]. Different domestic and foreign scholars have carried out relevant research on the flow and heat transfer characteristics of supercritical liquid hydrogen. For instance, Tomokazu et al. [5] carried out simulation experiments under different flow conditions to verify the feasibility of suppressing the local temperature rise to the range of 3 K in the design of a liquid hydrogen-cooled moderator. David [6] carried out experiments to improve the cooling efficiency of liquid hydrogen in a moderator. Hendricks and Simon [7] experimentally studied the flow and heat transfer characteristics of supercritical pressure hydrogen in elbows with different diameters and shapes under a uniform heat flux. Beech et al. [8] found that the flow and heat transfer phenomenon of supercritical hydrogen in elbow was affected by many factors, such as the physical property change and boundary conditions. The physical process was found to be complex, and it was difficult to summarize a universal phenomenon from limited experimental conditions. Zhou and Ji [9] simulated the flow and heat transfer of supercritical hydrogen in a 180° elbow. It was found that heat transfer outside the elbow was strengthened and that inside deteriorated.
In addition to supercritical liquid hydrogen, supercritical water and supercritical nitrogen [10,11] are also hot research topics. Many scholars have conducted numerical simulations based on different types of pipes, such as spiral pipes, loops and vertical circular pipes [12,13,14]. Wang et al. [15] carried out an experimental study on the heat transfer of aviation kerosene in a U-tube under supercritical pressure. Wu [16] discussed the effects of four factors, including inlet temperature, heat flux, inlet pressure and inlet mass flow, on the flow and heat transfer process of supercritical n-decane in a rectangular straight pipe channel. Yu [17] studied the flow and heat transfer characteristics of a supercritical CO2–propane binary mixture. The results showed that the heat transfer coefficient of the CO2–propane mixture decreases with increase of heat flux at same temperature.
There are few studies on the heat transfer characteristics of supercritical liquid hydrogen tubes with a U-shape and constant heat flux on the outer wall. The application value of a supercritical liquid hydrogen tube for a moderator and attention of the neutron source were considered in this study. In the current work, the flow characteristics of supercritical liquid hydrogen in a U-shaped tube were investigated, especially the heat transfer research when gravity direction changes in extreme conditions. Nowadays, the operating conditions of CSNS will be upgraded to G = 30 g/s, q = 500 kW/m2. Therefore, this work is of great significance for the thermal safety and emergency response of the target station.

2. Numerical Simulation

2.1. Physical Model

The U-tube geometry model is shown in Figure 1. According to a previous study of a spallation neutron source for moderators, most pipelines inside moderators can adopt a small diameter. In this study, the spatial direction of the U-tube was vertical, and the fluid flow direction was consistent with the gravity direction (gx). The total length of the U-tube was 1725 mm, withan inner diameter and wall thickness of D × δ = 10 mm × 2 mm. The curvature radius of the elbow was R = 27 mm, and the length of the heating section was L = L1 + L2 + L3 = 525 mm. To weaken the inlet section effect and outlet reflux phenomenon, before and after the heating section of the U-tube, namely, the inlet and outlet were respectively set with a 600 mm adiabatic inlet section and a 600 mm adiabatic outlet section.

2.2. Meshing

The calculation domain of the U-tube was established and meshed. As shown in Figure 2, the upper half grid diagram of the computational domain in the yoz direction and xoz section direction of the horizontal tube inlet were presented. Considering that the physical parameters of fluid in the boundary layer region near the wall will change dramatically, the flow and heat transfer characteristics were relatively complex, and the grid near the wall was encrypted. Through the finer grid, the flow and temperature field around the U-tube were investigated in detail. In this work, the enhanced wall function (y+ < 1) was used to divide the grid, which has certain advantages for studying the heat transfer and flow conditions of the wall.

2.3. Governing Equations

In the current work, it was assumed that fluid flow was in a steady turbulent state and did not have an internal heat source.
The continuity, momentum, and energy conservation equations were described as follows [18]:
The continuity equation:
( ρ u ) = 0
The momentum equation:
( ρ u u ) = p x + τ x x x + τ y x y + τ z x z + ρ f x ( ρ v u ) = p y + τ x y x + τ y y y + τ z y z + ρ f y ( ρ w u ) = p z + τ x z x + τ y z y + τ z z z + ρ f z
The energy equation:
[ ( ρ E + P ) u ] = ( λ e f f T ( τ e f f u ) ) + S h  
where ρ is density; u is the velocity vector; u, v, w are components of velocity u ; τ is the stress tensor; p is the bulk fluid pressure; λ is thermal conductivity; T is bulk fluid temperature; E is total energy; P is static pressure; τeff is shear stress; λeff is effective thermal conductivity; Sh is volumetric heat source.
In this study, the Re value under working conditions was much higher than its critical value, so the flow was in a turbulent state. During the simulation, the standard k-ε turbulence model was selected to describe the turbulent flow [19].
The transport equations are described as follows:
x i ( ρ k u ) = x j [ ( μ + μ t σ k ) k x j ] + G K + G b ρ ε Y M + S k
x i ( ρ ε u ) = x j [ ( μ + μ t σ ε ) ε x j ] + C 1 ε ε k ( G K + C 3 ε G b ) C 2 ε ρ ε 2 k + S ε
where,
μ t = ρ C μ k 2 ε
where Gk represents the generation of turbulence kinetic energy due to the mean velocity gradients, Gb is the generation of turbulence kinetic energy due to buoyancy, YM represents the contribution of the fluctuating dilatation incompressible turbulence to the overall dissipation rate, μ is dynamic viscosity, μt is turbulent viscosity, k is turbulent kinetic energy, ε is the dissipation rate of turbulent kinetic energy and E is total energy. σk and σε are the turbulent Prandtl numbers for k and ε, respectively. Sk and Sε are source terms. The model constants C1ε, C2ε, Cμ, σk, σε had the following default values: C1ε = 1.44, C2ε = 1.92, Cμ = 0.09, σk = 1.0, σε = 1.3 [20].

2.4. Data Reduction Method

The dimensionless independent parameter Re was obtained by Equation,
R e = ρ u D μ
The local Nusselt number Nu was defined as:
N u = h D λ
where D is the hydraulic tube diameter, and h is the local heat transfer coefficient.
The local heat transfer coefficient equation was as follows:
h = q u T w T b ,   T b = f ( H b , P ) ,   H b = A ρ u H D A A ρ u D A
where qu is the local heat flux of wall, Tw is the averaging temperature wall, Tb is the bulk temperature of fluid, u is the axial velocity of the fluid, A is the cross-sectional area of the pipe and H is the bulk liquid enthalpy.

2.5. Calculation Conditions

In this study, ANSYS CFX 2021 software was used to conduct a numerical simulation of heat transfer characteristics of supercritical liquid hydrogen in U-shaped tubes. The Semi-Implicit Method for Pressure-Linked Equations of pressure coupled equations was selected. The finite volume method was used to discrete the control equation and other equations in which the second-order upwind was selected for the discretization scheme. The room mean square (RMS) value of the residual value was set to 10−6, and the number of iterations was set to 10,000 times, until a stable independent convergence solution was obtained. The relaxation factors were all based on the default value. In the residual setting, except that the energy equation term was less than 10−6, the variables were less than 10−4.
The boundary conditions of the model were set to be simulated, including its fluid domain and solid domain.
  • Boundary condition setting: Firstly, the reference working condition was established, As shown in Figure 1, the gravity direction was the positive direction of the x-axis, and constant heat flow q was used in wall of the heating section. The boundary conditions are shown in the following.
    q = λ T x = 80 kW m 2 ,   G = 5   g / s ,   g x = 9.8   m / s
    Tin = 15 K, Pout = 1.5 MPa
  • Parameter settings of physical properties: Thermophysical properties of liquid hydrogen under supercritical conditions can be obtained from data provided by the National Institutes of Standards and Technology (NIST) [21]. The physical properties of liquid hydrogen can be regarded as a single value function of temperature. The corresponding thermal properties were inputted through the piecewise-linear interpolation function of CFX 2021. In this work, GASPAK software was used to calculate the thermophysical parameters of supercritical fluid hydrogen. When the temperature was inputted, the thermophysical parameters of the corresponding fluid were obtained, such as λ, cp, μ, ρ. The temperature range selected in this study was 15 K~40 K. The thermophysical parameters of supercritical liquid hydrogen were plotted by the obtained thermophysical parameters, as shown in Figure 3.
  • As shown in Figure 3, the thermophysical properties of liquid hydrogen change with temperature at P = 1.5 MPa. For supercritical fluids, they were near the pseudo-critical temperature (the temperature corresponding to the maximum specific heat capacity of fluid at supercritical constant temperature). The physical properties of fluid will change sharply, which will affect the heat transfer of fluid. In general, the critical pressure and temperature of liquid hydrogen are 1.29 MPa and 33.15 K, respectively. It can be seen from Figure 3 that density, ρ, and dynamic viscosity, μ, both showed a downward trend with the increase in temperature and, finally, tended to converge. The changing trend of specific heat, cp, and thermal conductivity, λ, was similar. Still, the changing trend of specific heat, cp, was more significant, which is not conducive to the stability of liquid hydrogen in heat transfer.

3. Results and Discussion

3.1. Verification of Grid Independence

The quality and quantity of grids have a significant influence on calculation results, so it is necessary to verify grid independence. In order to ensure the reliability of the near-wall calculation, the enhanced wall function was selected, and the y+ value should have been less than 1. In this study, five groups of grids were used to test the independence of grids. The maximum temperature of liquid hydrogen, the average flow rate and the maximum temperature of the outlet were monitored. As shown in Figure 4, simulation results of different grids are presented, among which five groups of grids were 0.56, 1.02, 1.53, 1.98 and 2.41 million, respectively. As the number of grids increased from 1.98 million, the trend of the three groups of variables tended to be gentle, and a convergence solution could be obtained. Considering the calculation accuracy and efficiency, the grid number of 1.98 million was selected for this study to ensure that the calculation results had grid independence.
The grid convergence index (GCI) [22] was defined as:
G C I = F s | ε | r p 1   ε = f 1 f 2 f
r k , k + 1 = N k + 1 N k 3
where Fs is the safety factor value of 1.25, f is the maximum values of f1 and f2, f1 and f2 are the outlet temperature values, Nk is the number of grids and p is the convergence accuracy, which is 1.97. GCI12 was 0.75%, GCI23 was 2.76%, GCI34 was 0.81% and GCI45 was 1.23%, and the values were less than 3%, indicating that the numerical simulation results were independent of number of grids. Considering the calculation accuracy and efficiency, the fourth data group was sufficient to ensure grid independence of the calculation results.

3.2. Sensitivity Analysis of Turbulence Model

To prove that the selected turbulence model well-described the flow of supercritical liquid hydrogen in the tube, a variety of turbulence models were chosen to compare with the experimental data as references at the same boundary conditions [23]. In the experiment, the thermal thermocouple and automatic voltage digitizing system were used to measure the pipe wall temperature and fluid temperature, respectively. In the condition of G = 4.5 kg/m2−45 kg/m2, P = 1 MPa−5 MPa, 12 groups of data were collected to compare the results of liquid hydrogen in the supercritical state and pressure critical state. The specific model was a circular tube with a diameter of 8.51 mm and a length of 0.9144 mm. The inlet flow rate of supercritical liquid hydrogen was 4.5 kg/m2, and the inlet temperature and pressure were 25.33 K and 4.947 MPa, respectively.
It can be seen from Figure 5 that the partition trends of the supercritical liquid hydrogen temperature and tube wall temperature calculated by different turbulence models were in good agreement with the experimental results. However, there were still some differences between them. By observing and analyzing the distribution of wall temperature and liquid hydrogen temperature, it was found that the standard k-ε model had the best consistency. The error of the liquid hydrogen temperature distribution was not more than 5%, and the maximum error between the experimental value and simulation value were within 20%. The overall temperature distribution of other models was lower than the experimental data, so the standard k-ε model was finally used as the turbulence model in the current work.

3.3. Analysis of Heat Transfer Characteristics under Benchmark Conditions

3.3.1. Analysis of Longitudinal Temperature, Density and Velocity under the Reference Condition

As mentioned above, the reference operating parameters were q = 80 kW/m2, G = 5 g/s and gx = −9.8 m/s2 to analyze the longitudinal and transverse velocity, temperature and density of the U-tube. Figure 6 is the distribution of velocity (a), temperature (b) and density (c) of xoy plane; the velocity streamline (positive along x-direction) (d) and the temperature iso-surface (e) in the computational domain. It can be seen from the diagram that liquid hydrogen flowed evenly before entering bending channel and separated after entering the bending pipe. A low temperature and high-density region were formed on the outer wall of the tube, and a high temperature and low-density region were formed on the inner wall of tube. Along the flow direction, the low-temperature zone was in front, and the high-temperature zone was in the back. Finally, in the back of the elbow, two regions produced a mixing phenomenon. The velocity distribution in the flow channel was different. In the first half of the bend, the velocity on the inner side was abruptly increased from 0.9 m/s to 1.2 m/s, while velocity on the outer side was decreased from 0.9 m/s to 0.8 m/s. In the latter half of the elbow, the inner flow velocity decreased gradually, while the outer flow velocity increased and was higher than the inner. The main reason was that liquid hydrogen was just entering elbow. Due to the sudden emergence of the centrifugal force, the fluid inside the tube had a lateral displacement from inside to outside and produced fluid extrusion. The flow velocity of the outside fluid was reduced under extrusion, while space inside became larger, and flow velocity increased significantly. Then, with the flow of fluid, the distribution was gradually uniform. Under the action of centrifugal force, the velocity of the outer fluid was gradually accelerated, and the velocity of the inner fluid was slower. Figure 7 illustrates the distribution of the temperature, density, and secondary flow trace on the different cross-sections of the bent tube. When the fluid flowed into the elbow, the lateral flow became obvious, and the Dean vortex was formed. The secondary flow on the cross-section took the cold fluid from the core to the outside of the U-tube elbow, and the heated fluid on the outside (near the wall) flowed along the wall to the inner-side and then reached the center of the pipe. The traces near the wall generated by the buoyancy force moved up gradually. In contrast, the Dean traces generated by the centrifugal force moved down gradually, and all of them could be seen in the α = 120° section. In addition, the thermal boundary layer became thinner and disappeared during circulation, and there was a sudden temperature difference between the wall and the fluid near the wall. Because under the influence of the enhanced secondary flow, a thin layer of fluid nearby may absorb more heat, and the bulk temperature of supercritical hydrogen in this layer may have increased rapidly due to the significant decrease in density. To summarize, the secondary flow phenomenon could conduce thinning of the boundary layer of supercritical hydrogen, thus enhancing heat transfer.

3.3.2. Analysis of Section Temperature, Density and Velocity at the Bend Pipe under the Reference Condition

Figure 7 shows the temperature, density and trace contour of the calculated elbow on the cross-section of α (0°, 60°, 90°, 120°, 180°) at different angles. In general, all contour images exhibited left–right symmetry. From Figure 7a, it can be seen that temperature and density distribution inside and outside the flow channel were the same near the entrance section of the elbow at α = 0°, which was consistent with the theoretical analysis. The temperature and density contour showed a concave shape in α = 90°−120°. It could be seen from the temperature distribution contour that after entering the elbow section, the temperature distribution was closely related to the velocity distribution, and the velocity field determined the temperature field. As mentioned above, the flow rate of liquid hydrogen inside the pipeline was lower than that of the outside, which led to the heat transfer capacity of the inside pipeline being weaker than that of the outside. Therefore, the temperature of liquid hydrogen inside the pipeline was higher than that outside. It can be seen from the trace diagram that the secondary flow double vortex generated by the combined action of the centrifugal force and buoyancy appeared at α = 60°, which was located inside, and the vortex area as well as the velocity increased with the increase of the α angle. From α = 90–180°, it can be seen that with the increase in angle, the vortex near the wall and internal vortex gradually move from the inside to the outside. At α = 180°, due to the disappearance of centrifugal force, only the secondary flow vortex generated by buoyancy was left, and the vortex at this time accounted for most of the region. The temperature, density and velocity distribution characteristics of liquid hydrogen in the elbow can be ultimately observed in Figure 6 and Figure 7.

3.3.3. Analysis of Tw, h and Nu in the Tube

Figure 8 shows the change of Tw, h and Nu with the direction of L/D when supercritical liquid hydrogen and conventional liquid hydrogen are heated in a U-tube. Point lines at L/D = 22, 26 and 30 indicate the upstream below the inlet, mid-point and downstream elbow outlet of the U-type. As shown in the figure, the overall trends of the two working fluids were similar. It could be found that Tw at the straight pipe section increased continuously along the flow direction of liquid hydrogen inside the pipe. After liquid hydrogen entered the elbow section (a = 0°), the wall temperature, h and Nu curves all showed inflection points. The wall temperature gradually decreased from 79° to 63.8°, and the temperature difference of convective heat transfer increased. The variation trends of convective heat transfer coefficients h and Nu were opposite to the wall temperature. This was mainly due to the centrifugal force, resulting in a secondary flow, so that the velocity of liquid hydrogen in the outer tube was higher than that of the inner tube. It can be seen from Figure 8 that at the elbow (α = 0°~180°), h and Nu between liquid hydrogen and the pipe wall increased significantly. More specifically, the heat transfer was strengthened. Finally, at the position flowing out of the elbow (α = 180°, L/D = 30), the centrifugal force inside the pipe disappeared, and the curve showed a turning point. Nu and h gradually decreased with the increase of the boundary layer. The temperature near the wall of the fluid slowly increased, and heat transfer deteriorated. During the flow of the fluid in the upstream and downstream straight pipes (L/D = 0~22, 30~40), the thermophysical properties of supercritical liquid hydrogen changed with temperature. Since Nu is the function of h and λ, the trends of Nu and h were not completely consistent. Because the heat transfer intensity determined the number of Tw, Tw presents a completely opposite trend with h and Nu, which directly reflects the heat transfer characteristics of the U-tube. The secondary flow generated by the centrifugal force at the elbow was dominant, and both mutations occurred. At the L/D = 23 to 30 position, the numerical mutation phenomenon appeared and reached the peak, respectively. The main reason is that this position was located in the inlet and outlet area of the elbow. At this time, the centrifugal force and the physical parameters of the supercritical liquid hydrogen were particularly violent, so the numerical value changed drastically. Due to the fact that the thermophysical properties of conventional liquid hydrogen do not change with temperature, the numerical changes at the elbow were not obvious. The overall trend of the Nu curve was similar to h, and the peak point of Nu and h appeared in the same position. Nu could accurately reflect the heat transfer of liquid hydrogen in the U-tube, which is the same as that of supercritical carbon dioxide [24]. For conventional liquid hydrogen, its thermophysical properties do not change with temperature. Therefore, compared with supercritical liquid hydrogen, the change of Tw, h and Nu in the conventional liquid hydrogen flow in a U-tube was relatively flat, and the heat transfer ability was also weaker.

3.4. Analysis of the Average Nu and Minimum Nu under Different Working Conditions

It was found that based on the above research, when the inlet mass flow rate G of liquid hydrogen in a single change tube was 3~10 (g/s), the heat flux q of the tube wall was 30~100 (kW/m2). The variation of Tw, h and Nu with L/D was similar to the above variables. The difference was as follows: After entering the elbow section, with the increase of heat flux, the valley point of Tw shifted to the left, and the peak point of Nu gradually shifted to left, which corresponded to the valley point of the wall temperature. The average and minimum Nu curves for above G range of 3~10 (g/s) and the tube wall heat flux q range of 30~100 (kW/m2) are shown in Figure 9a,b. At the same time, it was compared with conventional liquid hydrogen, and the final results of the two fluids were compared and analyzed.
As shown in Figure 9, under different Re and heat flux densities, the variation trends of average Nu and minimum Nu in two different fluids were the same. The inlet flow G = 5 g/s and heat flux q = 80 kW/m2 were taken as reference conditions for analysis. It can be seen from Figure 9a that heat flux was 80 kW/m2, as the reference condition for the analysis. The average Nu and the minimum Nu decreased from 5 g/s to 3 g/s by 25.9% and 26.8%, respectively, and increased from 5 g/s to 10 g/s by 56.6% and 58.9%, respectively. According to Figure 9b, when the mass flow rate was 5 g/s, the average Nu and minimum Nu decreased by 13.4% and 14.6%, respectively, from 80 kW/m2 to 30 kW/m2 and increase by 3.5% and 3.7%, respectively, from 80 kW/m2 to 100 kW/m2. It can be found from Figure 9a,b that the average Nu and minimum Nu of conventional liquid hydrogen were smaller than those of supercritical liquid hydrogen. Based on the above conclusions, it can be inferred that the heat transfer ability of supercritical liquid hydrogen is stronger than conventional liquid hydrogen. For a straight pipe, when Re reaches 1 × 105, it tends to be almost constant, but for an elbow pipe (α = 0°~180°), Nu increases with the increase of the Reynolds number and heat flux, showing an approximate linear upward trend. To sum up, the Reynolds number has a greater influence on the change of Nu.

3.5. Effect of Gravity Direction on Flow and Heat Transfer of Supercritical Liquid Hydrogen

In unique working conditions, disasters such as earthquakes may have occurred, and the gravity direction of the U-tube where the supercritical liquid hydrogen is located will change accordingly, which has a particular impact on the research of a spallation neutron source. Therefore, to explore the influence of different gravity directions on low and heat transfer of supercritical liquid hydrogen in U-tube, the benchmark conditions were established: G = 5 g/s, q = 80 kW.
As shown in Figure 10a,b, point lines at L/D = 22, 26 and 30 indicate the position of upstream elbow inlet, mid-point and downstream elbow outlet of the U-type. Figure 10c shows the direction of gravity. When gravity directions were (gx = −9.8 m/s2, 0, 0), (0, gy = −9.8 m/s2, 0, and (0, 0, gz = −9.8 m/s2), respectively, Tw and Nu changed with the direction of L/D, where (gx = −9.8 m/s2, 0, 0) were taken as condition parameters of the above study. It can be seen from Figure 10a that wall temperature was the highest when gravity was along the z-direction and the smallest when gravity was along the x-direction. When fluid flowed in the flow channel, the wall temperature decreased with the increase of gravity, and wall temperature decreased gradually after fluid entered the elbow section. The wall temperature rose as the fluid left the elbow. The numerical mutation at the elbow was due to the thermal properties of supercritical hydrogen and the centrifugal force, which is mentioned above.
It can be seen from Figure 10b that as gravity changed from the z-axis direction to the x-axis direction, Nu also increased, and the peak point of the Nu curve corresponded to the valley point of the wall temperature change curve. The main reason was that with the flow process of supercritical liquid hydrogen in the channel, there is a temperature difference in the channel to form a density difference. Under the action of gravity, buoyancy is formed in the U-tube and can produce a secondary flow to strengthen the heat transfer in the channel. In the flow channel, the direction of buoyancy is directly related to heat transfer ability. When buoyancy is in the radial direction, it is conducive to convective heat transfer. When fluid is in the tangential direction, it will lead to the deterioration of heat transfer [25]. When gravity is along x-axis, the same as the flow direction of liquid hydrogen, the heat transfer capacity was the strongest, and Nu was the largest. When along the z-axis direction, the velocity and temperature distribution in the flow channel were unevenly distributed, which deteriorated the heat transfer and minimized Nu.

4. Conclusions

The working conditions considered in this study, G = 3 g/s~10 g/s (Re = 16,425~54,750), q = 30 kW/m2~100 kW/m2 and gravity (gx, gy, gz), have a certain theoretical guiding significance for the construction of the CSNS. The flow and heat transfer technology of supercritical liquid hydrogen in the U-tube inside the moderator determine the normal operation of the moderator. Based on the importance of flow and heat transfer technology of supercritical liquid hydrogen in the U-tube, the following conclusions could be obtained in the elbow:
  • During the flow of supercritical liquid hydrogen in U-tube, the heat transfer in the elbow part of the U-tube (α = 0°~180°) was strengthened, and heat transfer deterioration occurred at the elbow outlet (α = 180° L/D = 30).
  • For the straight tube, when Re reached 1 × 105, Nu approached to be constant. However, for the vertical elbow, when Re and q were from 27,375 to 54,750 and 80 kW/m2 and 100 kW/m2, Nu increased by 56.6% and 3.7%, and Re had a more significant impact on Nu variation.
  • When the gravity direction was consistent with the flow direction of liquid hydrogen (gx direction), the Nu number inside the channel was 4.21%, which is 5.56% higher than that in the gy and gz direction, respectively. The heat transfer capacity of fluid in flow channel along x-axis direction was the strongest, followed by the y-axis direction, and the z-axis direction was the weakest.

Author Contributions

W.F. and Y.L. (Yiping Lu) are the main authors of this manuscript. W.F., Y.L. (Yiping Lu), F.S., L.M., S.W., Y.L. (Youlian Lu), L.Z., S.F. and J.T. contributed to this manuscript. S.W. and J.T. proposed the research topic and gave the idea of the scheme. J.T. and Y.L. (Youlian Lu) analyzed the data and contributed analysis tools. L.Z. and J.T. designed the overall framework and provided guiding technical support. J.T. and L.M. checked, reviewed, and revised the paper. J.T. performed the final proofreading and supervised this research. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Program for Guangdong Introducing Innovative and Enterpre-neurial Teams (project no. 2017ZT07S225), the National Natural Science Foundation of China (2017ZT07S225), the National Natural Science Foundation of China (Grant Nos. U1932219), the National Natural Science Foundation of China (No. 11805033) and State Key Laboratory of Nuclear Power Safety Monitoring Technology and Equipment (No. K-A2019.418).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Across-sectional area: m2
Dtube hydraulic diameter, mm
δwall thickness, mm
τstress tensor
σk, σεturbulent Prandtl numbers for k and ε
Tbulk fluid temperature, K
Etotal energy, kJ/kg
Pstatic pressure, Pa
Gk, Gbturbulence kinetic energy, kg/m2s
Lcharacteristic length, mm
Rcurvature radius of the bend, mm
Hbulk liquid enthalpy, kJ/kg
NuNussel number
R e Reynolds number
Twaveraging temperature of wall, K
Tbbulk temperature of fluid, K
rbend radius, mm
α bending angle
kturbulence kinetic energy, m2/s2
pbulk liquid pressure, Pa
uvelocity vector, m/s
hheat transfer coefficient, W/m2·K
ρ density, kg/m3
λ thermal conductivity of the fluid, W/m·K−1
εdissipation rate of turbulent kinetic energy, m2/s3
μ dynamic viscosity, Pa·S
μ tturbulent viscosity, Pa·S
qwlocal heat flux of the wall, kW/m2
Sk, Sεsource terms
ωaxial velocity, m/s
c p specific heat capacity, J/kg·K−1
νkinematic viscosity, m2/s
λeffeffective thermal conductivity, W/m·K−1
τeffshear stress, N/m2
Shvolumetric heat source

References

  1. Chen, Z.J. China Spallation Neutron Source (CSNS). J. Chin. Acad. Sci. 2011, 66, 726–728. [Google Scholar]
  2. Yin, W.; Liang, T.; Yu, Q. Neutronics design for the coupled para-hydrogen moderator for CSNS. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2011, 631, 105–110. [Google Scholar] [CrossRef]
  3. Freels, J.D. COMSOL Validation Progress on Supercritical Hydrogen Heat Transfer. Source 2008, 6, 22–24. [Google Scholar]
  4. Hu, J.; Li, Z.; Gao, Y.; Yang, J.; Liu, L. Investigation into Abnormal Heat Transfer Mechanism of Supercritical Carbon Dioxide in Vertically Upward Tube. Chin. J. Ship Res. 2021, 16, 170–179. [Google Scholar]
  5. Aso, T.; Kaminaga, M.; Terada, A.; Hino, R. Thermal-Hydraulic Experiments and Analyses on Cold Moderator. In Proceedings of the ICANS-XV, Tsukuba, Japan, 11 November 2000; pp. 6–9. [Google Scholar]
  6. Aso, T.; Monde, M.; Sato, H.; Hino, R.; Tatsumoto, H.; Kato, T. Experimental Investigation of Characteristics of Impinging Jet Heat Transfer and Application to JSNS Moderator Design. Trans. At. Energy Soc. Jpn. 2006, 5, 179–189. [Google Scholar] [CrossRef]
  7. Hendricks, R.C.; Simon, F.F. Heat transfer to hydrogen flowing in a curved tube. Chem. Prop. 1963, 7, 7–9. [Google Scholar]
  8. Beech, J.C. Forced Convection Heat Transfer to Supercritical Cryogenic Hydrogen: Part 1. Lit. Surv. 1969, 12, 47–50. [Google Scholar]
  9. Zhou, B.; Ji, Y.; Sun, J.; Sun, Y.; Shi, L. Numerical Investigations of Supercritical Hydrogen Flow and Heat Transfer in Curved Tube. Hedongli Gongcheng/Nucl. Power Eng. 2019, 40, 197–201. [Google Scholar]
  10. Ortiz, F.J.G.; López-Guirao, F.; Jiménez-Espadafor, F.J.; Benjumea, J.M. Heat Transfer Limitations in Supercritical Water Gasification. Energies 2022, 15, 177. [Google Scholar] [CrossRef]
  11. Toor, S.S.; Shah, A.A.; Sharma, K.; Seehar, T.H.; Pedersen, T.H.; Rosendahl, L.A. Bio-Crude Production from Protein-Extracted Grass Residue through Hydrothermal Liquefaction. Energies 2022, 15, 364. [Google Scholar] [CrossRef]
  12. Diao, L. Numerical Simulation of Flow and Heat Transfer Characteristics of Supercritical CO2 in Straight and Spiral Tubes. Master’s Thesis, Guangdong University of Technology, Guangzhou, China, 2020. [Google Scholar]
  13. Zhang, K.; Han, C.; Ren, J.; Zhou, Y.; Bi, M. Numerical Simulation on Heat Transfer of Supercritical LNG in Coil Tubes of Submerged Combustion Vaporizer. Huagong Xuebao/CIESC J. 2015, 66, 4788–4795. [Google Scholar]
  14. Yu, B.W.; Fan, Y.Y.; Wang, Y.S.; Hong, W.P.; Niu, X.J. Numerical study on heat transfer characteristics of supercritical CO2-DME mixture in vertical circular tube. Eng. Thermophys. 2021, 42, 481–488. [Google Scholar]
  15. Wang, Y.; Lu, Y.; Li, S.; Dong, M. Numerical Study on Heat Transfer of Supercritical-Pressure RP-3 Aviation Kerosene in U-Turn Circular Tubes. Huagong Xuebao/CIESC J. 2021, 72, 4639–4648. [Google Scholar]
  16. Wu, M.T. Numerical Study on Flow and Heat Transfer of n-Decane in Supercritical State; Northeast Electric Power University: Jilin City, China, 2021. [Google Scholar]
  17. Wu, X.M.; Zhao, R.; Wei, Z.F. Study on flow boiling heat transfer characteristics of CO2/propane mixture in horizontal tube. J. Eng. Thermophys. 2013, 34, 4. [Google Scholar]
  18. Keshavarzzadeh, A.H.; Ahmadi, P.; Safaei, M.R. Assessment and Optimization of an Integrated Energy System with Electrolysis and Fuel Cells for Electricity, Cooling and Hydrogen Production Using Various Optimization Techniques. Int. J. Hydrogen Energy 2019, 44, 21379–21396. [Google Scholar] [CrossRef]
  19. Han, C.L.; Ren, J.J.; Dong, W.P.; Bi, M.S. Numerical Investigation of Supercritical LNG Convective Heat Transfer in a Horizontal Serpentine Tube. Cryogenics (Guildf) 2016, 78, 1–13. [Google Scholar] [CrossRef]
  20. Matsson, J.E. An Intrduction to ANSYS Fluent 2021; SDS Publications: Mission, KS, USA, 2021; pp. 75–83. [Google Scholar]
  21. Lemmon, E.; McLinden, M.O.; Huber, M.L. NIST Reference Fluid Thermodynamic and Transport Properties-REFPROP. Chem. Prop. 2002, 23, v7. [Google Scholar]
  22. An, E.K.; Zhang, R.; Han, Y.F. Grid independence analysis of numerical simulation of multiphase turbulent combustion. Boil. Technol. 2018, 49, 5. [Google Scholar]
  23. Friedman, R.; Graham, R.W.; Hendricks, R.C.; Hsu, Y.Y. Experimental Heat-Transfer Results for Cryogenic Hydrogen Flowing in Tubes at Subcritical and Supercritical Pressures to 800 Pounds per Square Inch Absolute. Natl. Aeronaut. Space Admin. 2016, 16, 1–23. [Google Scholar]
  24. Zhu, L.; Lu, Y.; Tong, J.; Liang, T.; Lu, Y.; Fu, W.; Wang, B.; Zhang, Y. Sensitivity Analysis of Influencing Factors of Supercritical Methane Flow and Heat Transfer in a U-Tube. Energies 2021, 14, 5714. [Google Scholar] [CrossRef]
  25. Mao, S.; Zhou, T.; Wei, D.; Liu, W.; Zhang, Y. Heat transfer characteristics of supercritical water in channels: A systematic literature review of 20 years of research. Appl. Therm. Eng. 2021, 197, 117403. [Google Scholar] [CrossRef]
Figure 1. U-type tube model.
Figure 1. U-type tube model.
Energies 15 03605 g001
Figure 2. (a) Cross yoz section grid; (b) xoz section grid.
Figure 2. (a) Cross yoz section grid; (b) xoz section grid.
Energies 15 03605 g002
Figure 3. Curve of thermophysical properties of liquid hydrogen with temperature.
Figure 3. Curve of thermophysical properties of liquid hydrogen with temperature.
Energies 15 03605 g003
Figure 4. Grid independence validation of CHM.
Figure 4. Grid independence validation of CHM.
Energies 15 03605 g004
Figure 5. Sensitivity analysis of the turbulence model.
Figure 5. Sensitivity analysis of the turbulence model.
Energies 15 03605 g005
Figure 6. Velocity, temperature and density contours of xoy plane in U-type tube heating section.
Figure 6. Velocity, temperature and density contours of xoy plane in U-type tube heating section.
Energies 15 03605 g006aEnergies 15 03605 g006b
Figure 7. Temperature, density and velocity contours on the α section of the bend pipe with different bending angles.
Figure 7. Temperature, density and velocity contours on the α section of the bend pipe with different bending angles.
Energies 15 03605 g007
Figure 8. Variation of Tw, h and Nu with L/D.
Figure 8. Variation of Tw, h and Nu with L/D.
Energies 15 03605 g008
Figure 9. Comparison of average Nu and minimum Nu of supercritical and traditional liquid hedrogen under different conditions of G and q (a,b).
Figure 9. Comparison of average Nu and minimum Nu of supercritical and traditional liquid hedrogen under different conditions of G and q (a,b).
Energies 15 03605 g009
Figure 10. Variation of Tw, Nu with L/D under different gravity directions and gravity direction diagram.
Figure 10. Variation of Tw, Nu with L/D under different gravity directions and gravity direction diagram.
Energies 15 03605 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, W.; Lu, Y.; Shen, F.; Mei, L.; Wang, S.; Lu, Y.; Zhu, L.; Fu, S.; Tong, J. Study on Convective Heat Transfer Characteristics of Supercritical Liquid Hydrogen in a U-Type Tube inside a Moderator. Energies 2022, 15, 3605. https://doi.org/10.3390/en15103605

AMA Style

Fu W, Lu Y, Shen F, Mei L, Wang S, Lu Y, Zhu L, Fu S, Tong J. Study on Convective Heat Transfer Characteristics of Supercritical Liquid Hydrogen in a U-Type Tube inside a Moderator. Energies. 2022; 15(10):3605. https://doi.org/10.3390/en15103605

Chicago/Turabian Style

Fu, Weida, Yiping Lu, Fei Shen, Longwei Mei, Songlin Wang, Youlian Lu, Lingbo Zhu, Shinian Fu, and Jianfei Tong. 2022. "Study on Convective Heat Transfer Characteristics of Supercritical Liquid Hydrogen in a U-Type Tube inside a Moderator" Energies 15, no. 10: 3605. https://doi.org/10.3390/en15103605

APA Style

Fu, W., Lu, Y., Shen, F., Mei, L., Wang, S., Lu, Y., Zhu, L., Fu, S., & Tong, J. (2022). Study on Convective Heat Transfer Characteristics of Supercritical Liquid Hydrogen in a U-Type Tube inside a Moderator. Energies, 15(10), 3605. https://doi.org/10.3390/en15103605

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop