Next Article in Journal
First-Principles Calculations of the Electronic Structure and Optical Properties of Yttrium-Doped ZnO Monolayer with Vacancy
Next Article in Special Issue
Supercritical Phase Inversion: A Powerful Tool for Generating Cellulose Acetate-AgNO3 Antimicrobial Membranes
Previous Article in Journal
On the Microstructure and Isothermal Oxidation at 800 and 1200 °C of the Nb-24Ti-18Si-5Al-5Cr-5Ge-5Sn (at.%) Silicide-Based Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis on Heat Transfer Characteristics of Supercritical CO2 in Heated Vertical Up-Flow Tube

1
The Beijing Key Laboratory of Multiphase Flow and Heat Transfer, North China Electric Power University, Beijing 102206, China
2
Key Laboratory of Power Station Energy Transfer Conversion and System, North China Electric Power University, Ministry of Education, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(3), 723; https://doi.org/10.3390/ma13030723
Submission received: 6 January 2020 / Revised: 29 January 2020 / Accepted: 30 January 2020 / Published: 5 February 2020
(This article belongs to the Special Issue Supercritical Processing in Polymers and Aerogels)

Abstract

:
It is great significance to understand the mechanism of heat transfer deterioration of supercritical CO2 for heat exchanger design and safe operation in the supercritical CO2 Brayton cycle. Three-dimensional steady-state numerical simulation was performed to investigate the behavior of supercritical CO2 heat transfer in heated vertical up-flow tube with inner diameter di = 10 mm and heated length Lh = 2000 mm. Based on the characteristics of inverted-annular film boiling at subcritical pressure, the heat transfer model of supercritical CO2 flowing in the heated vertical tube was established in this paper. The mechanisms of heat transfer deterioration (HTD) and heat transfer recovery (HTR) for supercritical CO2 were discussed. Numerical results demonstrate that HTD is affected by multiple factors, such as the thickness and property of vapor-like film near the wall, the turbulence intensity near the interface between liquid-like and vapor-like, and in the liquid-like core region as well as the distribution of radial velocity vector. Among the above factors, the change of turbulent kinetic energy caused by the buoyancy effect seems to be a more important contributor to HTD and HTR. Furthermore, the influences of heat flux and mass flux on the distribution of wall temperature were analyzed, respectively. The reasons for the difference in wall temperature at different heat fluxes and mass fluxes were explained by capturing detailed thermal physical properties and turbulence fields. The present investigation can provide valuable information for the design optimization and safe operation of a supercritical CO2 heat exchanger.

1. Introduction

Compared with the water–steam Rankine cycle, the supercritical CO2 Brayton cycle has been widely considered for power conversion systems in various energy industries, such as nuclear energy, fossil fuel, solar energy, and geothermal energy due to its high efficiency and system compactness. From the perspective of physical property, CO2 has relatively low critical parameters (critical pressure Pcr = 7.38 MPa, critical temperature Tcr = 304.25 K), and it is difficult to oxidize between CO2 and metal materials, which makes it possible to raise steam parameters including pressure and temperature for improving thermal efficiency [1]. Therefore, supercritical CO2 Brayton cycle has gradually become a research issue all over the world [2,3,4,5,6]. It is no doubt that heat transfer from heat source to supercritical CO2 is carried out through heat exchanger. Consequently, one of the key issues relating supercritical CO2 Brayton cycle is the safety of the heat exchanger. A comprehensive understanding of the heat transfer characteristic of supercritical CO2 flowing in a circular tube is essential to engineering design and the safe operation of heat exchangers.
In general, three heat transfer modes are divided during supercritical fluid flowing in heated tube: normal heat transfer (NHT), heat transfer enhancement (HTE), and heat transfer deterioration (HTD). In order to distinguish from HTD, researchers usually refer to NHT and HTE as non-HTD mode [7]. Compared with non-HTD mode, HTD is characterized by pronounced peaks with sharply increasing wall temperature or an abrupt decrease of heat transfer coefficient h, which will cause an obvious temperature difference between the heated surface and bulk fluid. Further, the heated surface may then overheat and fail. Therefore, HTD is more concerned during heat exchanger design and operation. Currently, many researchers have paid great attention to HTD of supercritical fluids, mainly establishing criterion to judge the onset of HTD for supercritical fluids [8,9,10,11,12] and revealing HTD mechanism.
For the HTD mechanism, numerous studies have mainly focused on the importance of flow acceleration effect and buoyancy effect. Jackson [13] proposed a criterion that was written as Gr/Re2.7 for predicting buoyancy effect on heat transfer. Jackson recommended that buoyancy effect could be ignored during supercritical fluids flowing in vertical tube when Gr/Re2.7 < 10−5, otherwise, it should be considered. Koshizuka et al. [14] numerically analyzed the HTD phenomenon of supercritical water. They considered that HTD was primarily triggered by the buoyancy effect at lower mass flux, and HTD appeared at the boundary between natural and forced convection. Fan et al. [15] presented that the main factor for HTD was the buoyancy effect at low mass flux and high heat flux, while the major reasons were not yet clear at high mass flux. Huang et al. [16] reviewed that both a strong flow acceleration effect and buoyancy effect were induced simultaneously at high heat flux, which caused complex flow and heat transfer behavior. Moreover, Huang et al. also summarized the criteria to judge the influence of buoyancy on heat transfer. They noted that there existed inconsistencies between the criteria and experimental data. They suggested that a new correlation was necessary to reliably predict the buoyancy effect. Liao and Zhao [17] experimentally investigated the convection heat transfer of supercritical CO2 in heated vertical miniature tubes, and found that buoyancy effect was still prominent at high mass flux. In conclusion, both the flow acceleration effect and buoyancy effect cannot thoroughly explain the HTD of supercritical fluids, and thus it is necessary to further investigate the heat transfer mechanism of supercritical CO2.
In the present research, the numerical investigation was conducted to reveal the heat transfer mechanism of supercritical CO2 in a vertical up-flow tube, and the effect of heat flux and mass flux on heat transfer was analyzed. The main content of this paper was organized as follows. Firstly, based on the characteristics of supercritical CO2 physical properties and previous studies, the physical heat transfer model of supercritical CO2 flowing in a vertical tube was established. Secondly, the heat transfer mechanism of supercritical CO2 was revealed according to the heat transfer model. Finally, the difference in supercritical CO2 heat transfer at different heat fluxes and mass fluxes were analyzed. The present investigation can provide the valuable information for optimization design and safe operation of heat exchanger.

2. Thermal-Physical Properties and Idealized Heat Transfer Model

Figure 1a shows the thermal-physical properties (such as specific heat cp, density ρ, thermal conductivity λ and viscosity μ) of supercritical CO2 at P = 15.0 MPa. Different from the step change in properties at critical temperature under subcritical pressures, the thermal-physical properties of supercritical CO2 continuously change when temperature increases. As shown in Figure 1a, specific heat shows a peak at specified supercritical pressure. The temperature corresponding to maximum specific heat is defined as the pseudocritical temperature Tpc. Figure 1b shows the P/Pcr-T/Tcr curves at subcritical and supercritical pressures. According to classical thermodynamics, there exists no transition from liquid to vapor at supercritical pressures. However, the thermal-physical properties vary drastically and non-linearly within a narrow range of temperature near the pseudocritical temperature. Under supercritical pressures, the fluid behaves like vapor when its temperature is higher than Tpc, and like liquid when its temperature is lower than Tpc. The phenomenon is similar to the transition from liquid to vapor at subcritical pressures. Thus, for supercritical CO2, the fluid state is defined as liquid-like at T < Tpc, and is defined as vapor-like at T > Tpc. The pseudocritical temperature is regarded as the “pseudo-phase transition” temperature between liquid-like and vapor-like.
Holman et al. [18] visually analyzed forced convection heat transfer to supercritical Freon 12 in a vertical annulus and observed vapor trails. Holman et al. proposed that it seemed reasonable to name the phenomenon after “pseudo-boiling” even if there existed no liquid-vapor interface at supercritical temperature, because boiling was also caused by sharp density gradient at subcritical pressure. Wang et al. [19] concluded that the HTD of supercritical pressure water was mainly ascribed to low-density layer covering the heated surface from the macroscopic point of view. Zhang et al. [20] experimentally investigated the characteristics of supercritical CO2 heat transfer in heated vertical-flow tube, and mentioned that liquid-like near the tube wall transited to vapor-like when supercritical CO2 was gradually heated. They considered that the laminarization of vapor-like layer near the heated surface was a significant HTD mechanism.
Based on the above viewpoint, the physical heat transfer model of supercritical CO2 flowing in a heated vertical circular tube is established according to the characteristics of inverted-annular film boiling at subcritical pressure [21], which considers the effect of the low density fluid layer near the wall on heat transfer, as displayed in Figure 2. Under Tw,i > Tp,c > Tf, low density fluid that is treated as vapor-like appears near the wall, and high density fluid that is treated as liquid-like flows in the center of the tube and is separated from the heated surface by vapor-like fluid, where Tw,i represents inner wall temperature and Tf represents local fluid temperature in the core region. The heat transfer process in the physical model mainly involves three steps: heat transfer from tube wall to vapor-like film, subsequently from vapor-like to interface between vapor-like and liquid-like, and then from the interface to liquid-like core.
In contrast to phase transition and the existence of a vapor-liquid interface at subcritical pressures, it should be noted that “pseudo-phase transition” process of supercritical CO2 near pseudocritical temperature is continuous, and the surface tension of supercritical CO2 is zero [22]. For the traditional numerical simulation of two-phase heat transfer at subcritical pressure, the surface tension of phase interface should be considered in the multiphase flow model [23,24]. Obviously, the multiphase flow model to investigate the boiling process at subcritical pressures is not suitable for supercritical fluids. Therefore, the numerical investigation of the heat transfer of supercritical CO2 still falls within the scope of single-phase flow and heat transfer.

3. Numerical Model and Method

3.1. The Governing Equations

The finite volume method is adopted to solve the transport equations of mass, momentum, and energy. The above governing equations are solved in a Cartesian coordinate system and can be expressed as below:
Mass
ρ ¯ u ˜ i x i = 0
Momentum
ρ ¯ u ˜ i u ˜ j x j = p ¯ x i + x j μ u i x j + u j x i 2 3 δ i j u k x k ρ u i u j ¯ + ρ g i
Energy
ρ ¯ u ˜ i h ˜ x i = x i μ 1 Pr + μ t / μ Pr t h ˜ x i
where superscript ‒ represents the time average scalar, superscript ~ represents the Favre average scalar, while ρ, u, λ and h are density, velocity, thermal conductivity, and fluid enthalpy, respectively.
For supercritical fluids, the turbulence model is great significant to the accuracy of calculation results. Liu et al. [25] and Wang et al. [26] discussed the abilities of different turbulence models to predict the experimental data of supercritical fluids, and found that shear stress transport (SST) kω low Reynolds model gave a more accurate prediction than other turbulence models. Consequently, the SST kω turbulence model is used in the current simulation.
The transport equations for SST kω model can be described as below:
Turbulent kinetic energy
ρ ¯ u ˜ i k x i = x j μ + μ t σ k k x j + G k Y k
Specific dissipation rate
ρ ¯ u ˜ i ω x i = x j μ + μ t σ ω ω x j + G ω Y ω + D ω
μt in Equations (4) and (5) is calculated by
μ t = ρ k ω max 1 α * , F 2 2 S i j S i j α 1 ω 1
where
α * = α * α 0 * + Re t / R k 1 + Re t / R k , S i j = u i x j + u j x i
Gk, Yk, Gw, Yw, Dw in Equations (4) and (5) are calculated by
G k = u t 2 S i j S i j
Y k = ρ β * k ω
G ω = α α * α 0 + Re t / R ω 1 + Re t / R ω ρ μ t G k
Y ω = ρ β ω 2
D ω = 2 1 F 1 ρ σ ω , 2 1 ω k x j ω x j
In the above equations, the relevant constants can be referred in Ref. [27].

3.2. Physical Model and Gird Generation

Figure 3a,b shows the physical model of a vertical tube with an inner diameter of 10 mm. The computational domain includes three parts: inlet adiabatic section, heated section and outlet adiabatic section. The heated length is 2000 mm. The length of inlet adiabatic section is 600 mm to ensure the fluid being fully developed before supercritical CO2 enters the heated section. The length of outlet adiabatic section is also 600 mm to ensure that supercritical CO2 flows steadily at outlet of heated section. The wall thickness of the simulated tube is ignored.
The structured grid technology is employed to generate hexahedron grid in the computational domain by ANSYS ICEM 15.0 software (ANSYS ICEM CFD 15.0,). Figure 3c,d show the grid generation in axial and radial direction, respectively. The axial grids are uniformly distributed. The O-type grid generation method is employed for a cross section of the tube. The grids near the wall-fluid interface are refined enough to accurately capture the dramatic variation of thermal-physical properties, and the non-dimensional wall distance y+ of the first node near the wall is always far less than 1 to meet the SST kω low Reynolds turbulence model requirement [28]. In the present study, the distance from the first node to the inner wall is set as 0.0015 mm and the growth ratio is set as 1.1.

3.3. The Boundary Conditions

The ANSYS Fluent 15.0 is adopted to perform numerical simulation for the flow and heat transfer of supercritical CO2 in vertical tube at three-dimensional steady state. The SIMPLEC algorithm is employed to solve the equations coupling pressure and velocity. The second-order upwind scheme is used to discretize the governing equations in Section 3.1, and the under-relaxation factors are set as default to achieve faster convergence. The maximum residual values are set as 1.0 × 10−6 for the mass and momentum equations and 1.0 × 10−7 for energy equation. The computation is assumed to converge when the iteration stops and the fluid temperature at outlet remains constant, and (|qm,inlet-qm,outlet|/qm,inlet) < 0.1%, where qm represents mass flow rate at the same time. Wang et al. [29] analyzed the sensitivity of heat transfer to thermal-physical properties and considered that heat transfer coefficient was obviously affected by the thermal-physical properties. Thus, the NIST real gas model is employed and a dynamic connection will be formed between the REFPROP database and the Fluent solver. With each iteration, the REFPROP database will provide the corresponding thermal-physical properties to the Fluent solver, which can timely reflect the influence of drastic changes in thermal-physical properties on the characteristics of flow and heat transfer. Compared with the system pressure, the pressure drop along the tube can be neglected. Therefore, it is assumed that the system pressure is constant during flow and heat transfer [30]. The inlet and outlet boundary conditions are set as the mass-flow-inlet and pressure-outlet, respectively. The inner-wall surface is considered to be heated with a uniform and constant heat flux.

3.4. Grid Independence

In general, grid independence analysis should be performed to ensure the validity and accuracy of the simulated results. Table 1 shows seven grid systems for grid independence analysis. The influence of the number of nodes in axial and radial directions on numerical simulation is thoroughly considered in the current simulation.
Figure 4 shows the results of grid independence analysis using seven grid systems under P = 8.0 MPa, G = 1000 kg/m2s, qw = 200 kW/m2, Tin = 290 K. In Figure 4a, the Grid No. 4 is selected as the comparative item. It can be observed that there is no significant deviation in the calculated inner wall temperature Tw,i between Grid No. 3 and Grid No. 4. As shown in Figure 4b, the calculated inner wall temperature Tw,i between Grid No. 6 and Grid No. 7 are almost coincident when the result calculated by Grid No. 7 is selected as the standard. In addition, according to the Richardson extrapolation [31,32], the grid convergence index (GCI) is employed to further evaluate the grid quality. In general, it can be observed in Figure 4 that the calculation results are more sensitive to distributions of radial grids than that in axial direction. As a result, the coarse-grid, medium-grid, and fine-grid with grid numbers 1.13 × 106 (Grid No. 3), 2.03 × 106 (Grid No. 6), and 3.35 × 106 (Grid No. 7) are selected for the GCI calculation.
GCI is defined as
GCI = F s · ε r p 1
where Fs represents the factor of safety and Fs = 1.25 [31], ε represents the relative error. r represents the ratio of grids number for the fine-grid and coarse-grid, p represents the formal order of algorithm accuracy, and the ε, r, and p can be expressed as [33]
ε = f 1 f 2 f 1
r = N coarse N fine 1 D
p = ln f 3 f 2 / f 2 f 1 ln r
where f1, f2, and f3 are the numerical solutions for fine-grid, medium-grid, and coarse-grid, respectively. D represents the geometry dimension. N is the number of control volumes. The value of GCI for the inner wall temperature at outlet is 0.037%. Therefore, Grid No. 6 with 2.03 × 106 nodes is regarded as the optimal grid system in the subsequent numerical investigation in order to save computer resources and to simultaneously ensure the quality of the grid.

3.5. Model Validation

The experimental data from Ref. [12] are selected to validate the numerical model to ensure the accuracy of numerical simulation. The experimental data in Figure 5 are measured in Beijing Key Laboratory of Multiphase Flow and Heat Transfer. The experimental system, parameter measurements and uncertainties, as well as the date reduction method are detailed in Ref. [12], which are not repeated here. During the selection of experimental runs, the effects of operation pressure, heat flux, and mass flux on heat transfer are considered. The boundary conditions of experimental runs involved in this section are listed in Table 2, where inlet pressure P, mass flux G, heat flux qw, inlet temperature Tin and inlet Reynolds number Rein are included. Figure 5a–f shows the comparisons between predicted heat transfer coefficient and experimental results. The variation trends of the predicted heat transfer coefficient show good consistency with the experimental data. Consequently, the numerical model is capable of simulating the flow and heat transfer of the supercritical CO2 in the present study.
In Figure 5, the mean relative error (errorave) and normalized mean square error (errorNMS) between the predicted heat transfer coefficient hpre and experimental hexp are expressed as below:
error ave = 1 n i = 1 n e i × 100 %
error NMS = 1 n i = 1 n e i e ¯ × 100 %
where ei represents the error for a single data point, e i = h pre h exp / h exp . e ¯ represents the mean relative error.

4. Results and Discussion

In the current simulation, Tb represents the bulk fluid temperature and is defined as [34]
T b = A ρ u c p T d A / A ρ u c p d A
where A is the cross-section area of the circular tube, m2, ρ represents the local fluid density, kg/m3. u represents the local fluid velocity, m/s, and cp represents the local fluid specific heat, kJ/kg K.
The bulk fluid enthalpy ib can be defined as
i b = A ρ u i d A / A ρ u d A
where i represents the local fluid enthalpy, kJ/kg.

4.1. Mechanism of Heat Transfer Deterioration

In this section, Case 1 in Table 2 is selected to analyze the heat transfer mechanism of supercritical CO2. Figure 6 shows the variations of predicted inner wall temperature Tw,i against bulk enthalphy ib under P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2. In Figure 5a and Figure 6, it can be found that the characteristic of wall temperature Tw,i and heat transfer coefficient h along the tube is consistent with HTD definition mentioned in the Introduction section. The characteristic locations A, B, and C are designated to reveal the mechanism of supercritical CO2 heat transfer, as shown in Figure 6. For convenient analysis, the HTD stage is defined from location A to B, and the HTR stage is defined from location B to C.
Detailed information on the thermal-physical properties and turbulence fields at locations A, B, and C marked in Figure 6 is shown in Figure 7. It should be noted that Tpc is used to characterize the interface between liquid-like and vapor-like at steady state to evaluate the influence of vapor-like film on heat transfer. In other words, the distance from the location of Tpc to the tube wall indicates the vapor-like film thickness, as shown in Figure 7a. It is well known that the thicker the vapor-like film is, the farther large-property-variation region will deviate from tube wall. Further, the smaller the vapor-like film density is, the worse the film property characterized by thermal conductivity λ and specific heat cp is. Figure 7a,b shows that the thickness of vapor-like film gradually increases along the tube, and the thermal conductivity of vapor-like film generally decreases, which results in a large thermal resistance near the wall to prevent heat from the heated surface transferring to colder bulk fluid. The large specific heat region moves away from the tube wall at the same time, the ability of vapor-like film near the wall to absorb heat is weakened. The above phenomenon is very disadvantageous to the heat transfer of supercritical CO2. In other words, the thicker vapor-like film and the worse vapor-like property are the important reasons for HTD.
Figure 7c clearly depicts the radial distribution of axial velocity u at characteristic locations. Owing to the flow acceleration effect, the axial velocity of bulk fluid increases continuously. Especially, due to buoyancy effect, the axial velocity of vapor-like near the wall increases obviously and is deformed into M-shape, which has a momentous influence on the turbulence production. As shown in Figure 7d, because shear stress caused by the velocity gradient is gradually increased, the turbulent kinetic energy in vapor-like film increases generally along the tube. The characteristic can enhance heat transfer in vapor-like film near the wall. It should be noted that the change of turbulent kinetic energy in radial direction is mainly caused by local buoyancy effect near the wall.
Further, Figure 7d illustrates that the level of turbulent kinetic energy near the interface between liquid-like and vapor-like and in the liquid-like core region is gradually decreased at HTD stage, which means that the fluid near the vapor-liquid interface and in the core region approaches to laminarization. Owing to the inadequate momentum and heat diffusion, the heat passing through vapor-like film cannot be effectively transferred to the core region. “Heat transfer stagnation” almost occurs near vapor-liquid interface, and heat transfer is further seriously suppressed. As a result, the wall temperature rises sharply. At HTR stage, compared with the location B, although the vapor-like film is thicker at the location C, the turbulent kinetic energy near the vapor-liquid interface and in the core region is larger, which indicates the sufficient momentum and heat diffusion. Consequently, the heat from the heated surface can be effectively transferred to the liquid-like region, and then the wall temperature is decreased. Therefore, the ability of turbulence diffusion near the interface between the liquid-like and vapor-like, and in the liquid-like core region seems to be a more important reason contributing to HTD. Obviously, the mechanism of HTD analyzed by numerical simulation in the present investigation is different from that proposed in Ref. [20], where Zhang et al. considered that the laminarization of the vapor-like layer near the wall was an important HTD mechanism.
It is interesting to notice that the HTD and HTR seem to be also related to the distribution of the radial velocity vector, as shown in Figure 8. At the HTD stage, the radial velocity near the wall is initially toward the center of tube cross section, which implies that the fluid in cross section has a tendency to flow toward the center of tube, as illustrated in Figure 8a. As a result, the tube wall cannot be effectively cooled in radial direction, the occurrence of HTD is accelerated. However, near the highest wall temperature at HTD stage, the radial velocity near the wall begins to turn, and then at HTR stage, the radial velocity vector towards the tube wall extends gradually to central region, as clearly shown in Figure 8b,c. This trend enhances mixing between the hot and cold fluid in cross section, and the tube wall can be cooled by the fluid in radial direction, which contributes to enhancing heat transfer and wall temperature recovery.
Therefore, the occurrence of HTD is caused by multi-factors including the thickness and property of vapor-like film near the wall, turbulent kinetic energy distribution affected by buoyancy as well as radial velocity vector. In other words, it is not only the buoyancy effect that leads to HTD.

4.2. The Effect of Heat Flux

Case 3 and Case 4 in Table 2 are selected to evaluate the effect of heat fluxes on the flow and heat transfer of supercritical CO2. The distribution curves of predicted inner wall temperature Tw,i against bulk enthalphy ib at different heat fluxes are shown in Figure 9. At the relatively high heat flux (qw = 351.22 kW/m2), the wall temperature is obviously increased, even HTD is more serious. Heat transfer is also deteriorated at the relatively low heat flux (qw = 294.50 kW/m2) in present simulation, but Tw,i is gently changed. As shown in Figure 9, the characteristic locations A and A’ (ib,A = ib,A’ = 276.0 kJ/kg), B and B’ (ib,B = ib,B’ = 327.8 kJ/kg) are marked to analyze the reason for the difference in heat transfer of supercritical CO2 at different heat fluxes.
The distributions of thermal-physical properties and turbulent flow at locations A and A’, B and B’ are presented in Figure 10. At high heat flux, the vapor-like film is thicker, and the property of vapor film characterized by specific heat and thermal conductivity is worse, as shown in Figure 10a,b. As a result, the ability of fluid near the wall to absorb and transfer heat is weakened. In other words, the thermal resistance caused by vapor-like film is larger at high heat flux. In addition, as can be found in Figure 10c,d, there is no obvious difference in the turbulent flow fields between location A and A’. However, compared with the location B’, due to a more significant buoyancy effect, the axial velocity u near the wall increases more obviously at location B. Therefore, the shear stress near the wall is increased, which is ultimately reflected by the distribution of turbulent kinetic energy, as displayed in Figure 10d. This implies that a change of turbulent flow near the wall can promote heat from the heated surface to pass through the vapor-like film at high heat flux. Comparing the locations B and B’, it should be noted that there is little difference in turbulent kinetic energy in the core region.
Consequently, due to the thicker vapor-like film and the worse film property near the wall, the more heat cannot be removed in a timely manner from the heated surface at relatively high heat flux, which contributes to higher wall temperature and even more severe HTD.
In addition, Figure 11 presents the distributions of radial velocity vector at characteristic locations at different heat fluxes. The magnitude order of velocity between high and low heat fluxes is the same in present simulated case, which illustrates that the distribution of radial velocity vector has little effect on the difference in wall temperature at different heat fluxes. It is noteworthy that the distribution characteristics of radial velocity vector at HTD stages shown in Figure 11a,b are similar to those of Case 1 in Section 4.1. Thus, it is further demonstrated that the distribution of wall temperature along the tube is also related to the radial velocity vector at specified operation parameters during supercritical CO2 heat transfer.

4.3. The Effect of Mass Flux

Case 5 and Case 6 in Table 2 are selected to evaluate the effect of mass fluxes on heat transfer of supercritical CO2. The distribution curves of predicted inner wall temperature Tw,i against bulk enthalphy ib at different mass fluxes are shown in Figure 12. It is observed that Tw,i increases as the mass flux decreases. In the present analysis, although the HTD occurs at both high and low mass fluxes, the amplitude of wall temperature overshoot is smaller at the relatively high mass flux. It illustrates that high mass flux can suppress HTD occurrence. The characteristic locations A and A’ (ib,A = ib,A’ = 264.6 kJ/kg), B and B’ (ib,B = ib,B’ = 280.3 kJ/kg) are marked to analyze the reason for difference in heat transfer of supercritical CO2 at different mass fluxes, as shown in Figure 12.
Figure 13 shows the distributions of thermal-physical properties and turbulent flow at the locations A and A’, B and B’ marked in Figure 12. At low mass flux, the vapor-like film is thicker, and its property is worse, as shown in Figure 13a,b. This illustrates that the thermal resistance of vapor-like film covering the heated surface is larger at low mass flux. Simultaneously, at low mass flux, the lower axial velocity cannot timely remove heat from the heated surface, and the turbulent kinetic energy is also generally smaller, especially in the liquid-like core region, which weakens heat transfer in the radial direction, as shown in Figure 13c,d.
Figure 14 shows the distributions of radial velocity vector at characteristic locations under different mass fluxes. It is of interest to see that, in contrast to the distribution of radial velocity vector at high mass flux, the velocity vector direction of supercritical CO2 in the core region is always toward the tube wall at low mass flux, as shown in Figure 14a,b. However, at different mass fluxes, there is a common ground that the velocity vector near the wall is always toward the center of tube at the beginning of the HTD stage. At low mass flux, it should be noted that velocity vector also begins to turn near the highest wall temperature at HTD stage, as shown in Figure 14b. This conclusion that the radial velocity vector affects the heat transfer of supercritical CO2 is consistent with Section 4.1 and Section 4.2.

5. Conclusions

In this paper, the numerical analysis was performed to investigate the flow and heat transfer characteristics of supercritical CO2 in heated vertical up-flow tube. The mechanisms of HTD and HTR for supercritical CO2 were revealed. The reasons for the difference in wall temperature distribution at different heat fluxes and mass fluxes were analyzed. Based on the present study, the main conclusions are presented as follows:
  • The occurrence of HTD is ascribed to the comprehensive effects of multiple factors. The vapor-like film with a large thermal resistance covering the heated surface prevents heat from the tube wall from transferring to the cooler bulk fluid. Simultaneously, it is another important reason for HTD that the fluid near the interface between vapor-like and liquid-like, and in the core region gradually approaches laminarization. Due to inadequate momentum and heat diffusion, heat transfer in the core fluid is seriously suppressed. Further, the radial velocity vector of fluid near the wall is initially toward the center of the tube, which results in the heated surface failing to be cooled effectively by the fluid in radial direction.
  • The formation of HTR is mainly attributed to two factors. The increased turbulent kinetic energy promotes heat transfer in radial direction, especially near the interface between liquid-like and vapor-like and in the liquid-like core region. The radial velocity vector of supercritical CO2 is towards the tube wall at HTR stage, thus the tube wall can be cooled by the fluid in radial direction.
  • For the relatively high heat flux, owing to the thicker vapor-like film and the worse film property near the wall, the more heat cannot be removed timely from the heated wall, which results in a higher wall temperature and even more serious HTD. For the relatively low mass flux, the vapor-like film is thicker, and its property is worse, and the lower axial velocity and the level of turbulent kinetic energy weaken heat transfer in axial and radial directions, which result in HTD being more likely to occur.

Author Contributions

Conceptualization, C.Y. and J.X.; methodology, C.Y. and J.X.; validation, C.Y., J.X., B.Z. and G.L.; investigation, C.Y., J.X., B.Z. and G.L.; resources, J.X.; data curation, C.Y. and J.X.; writing—original draft preparation, C.Y., J.X., B.Z. and G.L.; writing—review and editing, C.Y., J.X., B.Z. and G.L.; visualization, C.Y., J.X., B.Z. and G.L.; supervision, C.Y. and J.X.; funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China, grant number 2017YFB0601801 and the Fundamental Research Funds for the Central Universities, grant number 2018ZD02.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Holcomba, G.R.; Carney, C.; Doğan, Ö.N.F. Oxidation of alloys or energy applications in supercritical CO2 and H2O. Corros. Sci. 2016, 109, 22–35. [Google Scholar] [CrossRef]
  2. Dostal, V. A Supercritical Carbon Dioxide Cycle for Next Generation nNuclear Reactors. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2004. [Google Scholar]
  3. Al-Sulaiman, F.A.; Atif, M. Performance comparison of different supercritical carbon dioxide Brayton cycles integrated with a solar power tower. Energy 2015, 82, 61–71. [Google Scholar] [CrossRef]
  4. Mecheri, M.; Moullec, Y.L. Supercritical CO2 Brayton cycles for coal-fired power plants. Energy 2016, 103, 758–771. [Google Scholar] [CrossRef]
  5. Yang, Y.; Bai, W.; Wang, Y.; Zhang, Y.F.; Li, H.Z.; Yao, M.Y.; Wang, H.J. Coupled simulation of the combustion and fluid heating of a 300MW supercritical CO2 boiler. Appl. Therm. Eng. 2017, 113, 259–267. [Google Scholar] [CrossRef]
  6. Xu, J.L.; Sun, E.H.; Li, M.J.; Liu, H.; Zhu, B.G. Key issues and solution strategies for supercritical carbon dioxide coal fired power plant. Energy 2018, 157, 227–246. [Google Scholar] [CrossRef]
  7. Zhang, Q.; Li, H.X.; Lei, X.L.; Zhang, J.; Kong, X.F. Study on identification method of heat transfer deterioration of supercritical fluids in vertically heated tubes. Int. J. Heat Mass Transf. 2018, 127, 674–686. [Google Scholar] [CrossRef]
  8. Yamagata, K.; Nishikawa, K.; Hasegawa, S.; Fujii, T.; Yoshida, S. Forced convective heat transfer to supercritical water flowing in tubes. Int. J. Heat Mass Transf. 1972, 15, 2575–2593. [Google Scholar] [CrossRef]
  9. Mokry, S.; Pioro, I.; Farah, A.; King, K.; Gupta, S.; Peiman, W.; Kirillov, P. Development of supercritical water heat-transfer correlation for vertical bare tubes. Nucl. Eng. Des. 2011, 241, 1126–1136. [Google Scholar] [CrossRef]
  10. Chen, W.W.; Fang, X.D. A new heat transfer correlation for supercritical water flowing in vertical tubes. Int. J. Heat Mass Transf. 2014, 78, 156–160. [Google Scholar] [CrossRef]
  11. Kline, N.; Feuerstein, F.; Tavoularis, S. Onset of heat transfer deterioration in vertical pipe flows of CO2 at supercritical pressures. Int. J. Heat Mass Transf. 2018, 118, 1056–1068. [Google Scholar] [CrossRef]
  12. Zhu, B.G.; Xu, J.L.; Wu, X.M.; Xie, J.; Li, M.J. Supercritical “boiling” number, a new parameter to distinguish two regimes of carbon dioxide heat transfer in tubes. Int. J. Therm. Sci. 2019, 136, 254–266. [Google Scholar] [CrossRef]
  13. Jackson, J.D. Fluid flow and convective heat transfer to fluids at supercritical pressure. Nucl. Eng. Des. 2013, 264, 24–40. [Google Scholar] [CrossRef]
  14. Koshizuka, S.; Takano, N.; Oka, Y. Numerical analysis of deterioration phenomena in heat transfer to supercritical water. Int. J. Heat Mass Transf. 1995, 38, 3077–3084. [Google Scholar] [CrossRef]
  15. Fan, Y.H.; Tang, G.H.; Li, X.L.; Yang, D.L.; Wang, S.Q. Correlation evaluation on circumferentially average heat transfer for supercritical carbon dioxide in non-uniform heating vertical tubes. Energy 2019, 170, 480–496. [Google Scholar] [CrossRef]
  16. Huang, D.; Wu, Z.; Sunden, B.; Li, W. A brief review on convection heat transfer of fluids at supercritical pressures in tubes and the recent progress. Appl. Energy 2016, 162, 494–505. [Google Scholar] [CrossRef]
  17. Liao, S.M.; Zhao, T.S. An experimental investigation of convection heat transfer to supercritical carbon dioxide in miniature tubes. Int. J. Heat Mass Transf. 2002, 45, 5025–5034. [Google Scholar] [CrossRef]
  18. Holman, J.P.; Rea, S.N.; Howard, C.E. Forced convection heat transfer to Freon12 near the critical state in a vertical annulus. Int. J. Heat Mass Transf. 1965, 8, 1095–1100. [Google Scholar] [CrossRef]
  19. Wang, J.G.; Li, H.X.; Yu, S.Q.; Chen, T.K. Investigation on the characteristics and mechanisms of unusual heat transfer of supercritical pressure water in vertically-upward tubes. Int. J. Heat Mass Transf. 2011, 54, 1950–1958. [Google Scholar] [CrossRef]
  20. Zhang, S.J.; Xu, X.X.; Liu, C.; Liu, X.X.; Dang, C.B. Experimental investigation on the heat transfer characteristics of supercritical CO2 at various mass flow rates in heated vertical-flow tube. Appl. Therm. Eng. 2019, 157, 113687. [Google Scholar] [CrossRef]
  21. Elias, E.; Chambre, P. Inverted-annular film boiling heat transfer from vertical surfaces. Nucl. Eng. Des. 1981, 64, 249–257. [Google Scholar] [CrossRef]
  22. Jasper, J.J. The Surface Tension of Pure Liquid Compounds. J. Phys. Chem. Ref. Data 1972, 1, 841. [Google Scholar] [CrossRef] [Green Version]
  23. Yang, Z.; Peng, X.F.; Ye, P. Numerical and experimental investigation of two phase flow during boiling in a coiled tube. Int. J. Heat Mass Transf. 2008, 51, 1003–1016. [Google Scholar] [CrossRef]
  24. Da, R.E.; Del, C.D. Numerical Simulation of Laminar Liquid Film Condensation in a Horizontal Circular Minichannel. J. Heat Transf. 2012, 134, 051019. [Google Scholar]
  25. Liu, L.; Xiao, Z.J.; Yan, X.; Zeng, X.K.; Huang, Y.P. Heat transfer deterioration to supercritical water in circular tube and annular channel. Nucl. Eng. Des. 2013, 255, 97–104. [Google Scholar] [CrossRef]
  26. Wang, K.Z.; Xu, X.X.; Wu, Y.Y.; Liu, C.; Dang, C.B. Numerical investigation on heat transfer of supercritical CO2 in heated helically coiled tubes. J. Supercrit. Fluids 2015, 99, 112–120. [Google Scholar] [CrossRef]
  27. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  28. Li, Z.H.; Wu, Y.X.; Tang, G.L.; Zhang, D.L.; Lu, J.F. Comparison between heat transfer to supercritical water in a smooth tube and in an internally ribbed tube. Int. J. Heat Mass Transf. 2015, 84, 529–541. [Google Scholar] [CrossRef]
  29. Wang, X.C.; Xiang, M.J.; Huo, H.J.; Liu, Q. Numerical study on nonuniform heat transfer of supercritical pressure carbon dioxide during cooling in horizontal circular tube. Appl. Therm. Eng. 2018, 141, 775–787. [Google Scholar] [CrossRef]
  30. Shen, Z.; Yang, D.; Wang, S.Y.; Wang, W.Y.; Li, Y.D. Experimental and numerical analysis of heat transfer to water at supercritical pressures. Int. J. Heat Mass Transf. 2017, 108, 1676–1688. [Google Scholar] [CrossRef]
  31. Roache, P.J. Quantification of uncertainty in computational fluid dynamics. Annu. Rev. Fluid Mech. 1997, 29, 123–160. [Google Scholar] [CrossRef] [Green Version]
  32. Mansour, A.; Laurien, E. Numerical error analysis for three-dimensional CFD simulations in the two-room model containment THAI+: Grid convergence index, wall treatment error and scalability tests. Nucl. Eng. Des. 2018, 326, 220–233. [Google Scholar] [CrossRef]
  33. Luo, C.; Wu, S.; Song, K.W.; Hua, L.; Wang, L.B. Thermo-hydraulic performance optimization of wavy fin heat exchanger by combining delta winglet vortex generators. Appl. Therm. Eng. 2019, 163, 114343. [Google Scholar] [CrossRef]
  34. Wang, L.L.; Chen, Z.J.; Meng, H. Numerical study of conjugate heat transfer of cryogenic methane in rectangular engine cooling channels at supercritical pressures. Appl. Therm. Eng. 2013, 54, 237–246. [Google Scholar] [CrossRef]
Figure 1. Thermal-physical properties and thermodynamic state plane of supercritical CO2: (a) The thermal-physical properties curves at P = 15.0 MPa; (b) The P/PcrT/Tcr curves characterizing the transition between liquid and vapor at subcritical pressures and between liquid-like and vapor-like at supercritical pressures.
Figure 1. Thermal-physical properties and thermodynamic state plane of supercritical CO2: (a) The thermal-physical properties curves at P = 15.0 MPa; (b) The P/PcrT/Tcr curves characterizing the transition between liquid and vapor at subcritical pressures and between liquid-like and vapor-like at supercritical pressures.
Materials 13 00723 g001
Figure 2. Schematic physical model of heat transfer for supercritical CO2 in present study.
Figure 2. Schematic physical model of heat transfer for supercritical CO2 in present study.
Materials 13 00723 g002
Figure 3. Physical model and computational grids adopted in the present simulation: (a) and (b) Physical model and coordinate system; (c) and (d) Grid generation in computational domain.
Figure 3. Physical model and computational grids adopted in the present simulation: (a) and (b) Physical model and coordinate system; (c) and (d) Grid generation in computational domain.
Materials 13 00723 g003
Figure 4. Validation of grid independence: (a) The influence of the number of axial nodes on the calculation results; (b) The influence of the number of radial nodes on the calculation results.
Figure 4. Validation of grid independence: (a) The influence of the number of axial nodes on the calculation results; (b) The influence of the number of radial nodes on the calculation results.
Materials 13 00723 g004
Figure 5. Comparison between predicted heat transfer coefficient h and experimental data. (a) P=8.221 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (b) P=20.821 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (c) P=15.565 MPa, G=1001.5 kg/m2s, qw=351.22 kW/m2; (d) P=15.565 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (e) P=8.221 MPa, G=744.5 kg/m2s, qw=244.33 kW/m2; (f) P=8.221 MPa, G=1001.5 kg/m2s, qw=244.33 kW/m2.
Figure 5. Comparison between predicted heat transfer coefficient h and experimental data. (a) P=8.221 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (b) P=20.821 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (c) P=15.565 MPa, G=1001.5 kg/m2s, qw=351.22 kW/m2; (d) P=15.565 MPa, G=1001.5 kg/m2s, qw=294.50 kW/m2; (e) P=8.221 MPa, G=744.5 kg/m2s, qw=244.33 kW/m2; (f) P=8.221 MPa, G=1001.5 kg/m2s, qw=244.33 kW/m2.
Materials 13 00723 g005
Figure 6. Variations of predicted inner wall temperature against bulk enthalphy at P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2.
Figure 6. Variations of predicted inner wall temperature against bulk enthalphy at P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2.
Materials 13 00723 g006
Figure 7. Thermal-physical properties and turbulence fields at characteristic locations A, B and C (P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Figure 7. Thermal-physical properties and turbulence fields at characteristic locations A, B and C (P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Materials 13 00723 g007
Figure 8. Distributions of radial velocity vector at characteristic locations A, B and C (P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2). (a) location A; (b) location B; (c) location C.
Figure 8. Distributions of radial velocity vector at characteristic locations A, B and C (P = 8.221 MPa, G = 1001.5 kg/m2s and qw = 294.50 kW/m2). (a) location A; (b) location B; (c) location C.
Materials 13 00723 g008
Figure 9. Variations of inner wall temperature against bulk enthalphy at different heat fluxes (P = 15.565 MPa, G = 1001.5 kg/ms).
Figure 9. Variations of inner wall temperature against bulk enthalphy at different heat fluxes (P = 15.565 MPa, G = 1001.5 kg/ms).
Materials 13 00723 g009
Figure 10. Thermal-physical properties and turbulent flow fields at characteristic locations (P = 15.565 MPa, G = 1001.5 kg/m2s): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Figure 10. Thermal-physical properties and turbulent flow fields at characteristic locations (P = 15.565 MPa, G = 1001.5 kg/m2s): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Materials 13 00723 g010
Figure 11. Distributions of radial velocity vector at characteristic locations at different heat fluxes (P = 15.565 MPa, G = 1001.5 kg/m2s): (a) and (b) qw = 351.22 kW/m2; (c) and (d) qw = 294.50 kW/m2.
Figure 11. Distributions of radial velocity vector at characteristic locations at different heat fluxes (P = 15.565 MPa, G = 1001.5 kg/m2s): (a) and (b) qw = 351.22 kW/m2; (c) and (d) qw = 294.50 kW/m2.
Materials 13 00723 g011
Figure 12. Variations of inner wall temperature against bulk enthalphy at different mass fluxes (P = 8.221 MPa, qw = 244.33 kg/m2s).
Figure 12. Variations of inner wall temperature against bulk enthalphy at different mass fluxes (P = 8.221 MPa, qw = 244.33 kg/m2s).
Materials 13 00723 g012
Figure 13. Thermal-physical properties and turbulent flow fields at characteristic locations (P = 8.221 MPa, qw = 244.33 kg/m2s): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Figure 13. Thermal-physical properties and turbulent flow fields at characteristic locations (P = 8.221 MPa, qw = 244.33 kg/m2s): (a) Thermal conductivity; (b) Specific heat; (c) Axial velocity; (d) Turbulent kinetic energy.
Materials 13 00723 g013
Figure 14. Distributions of radial velocity vector at characteristic locations at different mass fluxes (P = 8.221 MPa, qw = 244.33 kg/m2s): (a) and (b) G = 744.3 kg/m2s; (c) and (d) G = 1001.5 kg/m2s.
Figure 14. Distributions of radial velocity vector at characteristic locations at different mass fluxes (P = 8.221 MPa, qw = 244.33 kg/m2s): (a) and (b) G = 744.3 kg/m2s; (c) and (d) G = 1001.5 kg/m2s.
Materials 13 00723 g014
Table 1. Cases for grid independence analysis.
Table 1. Cases for grid independence analysis.
Grid No.Number of NodesTotal Number of Nodes
RadialAxial
1441000.29 × 106
2442300.67 × 106
3443901.13 × 106
4444801.39 × 106
5503901.44 × 106
6653902.03 × 106
7753903.35 × 106
Table 2. Boundary conditions of experiment runs for model validation [12].
Table 2. Boundary conditions of experiment runs for model validation [12].
Case No.P (MPa)G (kg/m2s)qw (kW/m2)Tin (K)Rein
18.2211001.5294.50298.01.473 × 105
220.8211001.5294.50299.11.051 × 105
315.5651001.5351.22298.01.148 × 105
415.5651001.5294.50295.81.112 × 105
58.221744.5244.33296.51.053 × 105
68.2211001.5244.33296.51.417 × 105

Share and Cite

MDPI and ACS Style

Yan, C.; Xu, J.; Zhu, B.; Liu, G. Numerical Analysis on Heat Transfer Characteristics of Supercritical CO2 in Heated Vertical Up-Flow Tube. Materials 2020, 13, 723. https://doi.org/10.3390/ma13030723

AMA Style

Yan C, Xu J, Zhu B, Liu G. Numerical Analysis on Heat Transfer Characteristics of Supercritical CO2 in Heated Vertical Up-Flow Tube. Materials. 2020; 13(3):723. https://doi.org/10.3390/ma13030723

Chicago/Turabian Style

Yan, Chenshuai, Jinliang Xu, Bingguo Zhu, and Guanglin Liu. 2020. "Numerical Analysis on Heat Transfer Characteristics of Supercritical CO2 in Heated Vertical Up-Flow Tube" Materials 13, no. 3: 723. https://doi.org/10.3390/ma13030723

APA Style

Yan, C., Xu, J., Zhu, B., & Liu, G. (2020). Numerical Analysis on Heat Transfer Characteristics of Supercritical CO2 in Heated Vertical Up-Flow Tube. Materials, 13(3), 723. https://doi.org/10.3390/ma13030723

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