Issue 
EPJ Nuclear Sci. Technol.
Volume 2, 2016



Article Number  35  
Number of page(s)  11  
DOI  https://doi.org/10.1051/epjn/2016030  
Published online  09 September 2016 
https://doi.org/10.1051/epjn/2016030
Regular Article
Effect of heat transfer correlations on the fuel temperature prediction of SCWRs
^{1}
Departamento de Sistemas Energéticos, Facultad de Ingeniería, Universidad Nacional Autónoma de México,
C.P. 62550
Jiutepec,
Mor., Mexico
^{2}
Área de Ingeniería en Recursos Energéticos, Universidad Autónoma MetropolitanaIztapalapa,
C.P. 09340
México,
D.F., Mexico
^{3}
Sabbatical leave at the Facultad de Ingeniería of the Universidad Nacional Autónoma de México through the Programa de Estancias Sabáticas del CONACyT,
México,
D.F., Mexico
^{⁎} email: yurihillel@gmail.com
Received:
9
June
2015
Received in final form:
17
May
2016
Accepted:
20
July
2016
Published online: 9 September 2016
In this paper, we present a numerical analysis of the effect of different heat transfer correlations on the prediction of the cladding wall temperature in a supercritical water reactor at nominal operating conditions. The neutronics process with temperature feedback effects, the heat transfer in the fuel rod, and the thermalhydraulics in the core were simulated with a threepass core design.
© E.G. EspinosaMartínez et al., published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The super critical water reactor (SCWR) is one of the most promising and innovative designs selected by the Generation IV International Forum. This is a very highpressure watercooled reactor which will operate at conditions above the thermodynamic critical point. Water enters the reactor core and then exits without change of phase, i.e., no water/steam separation is necessary. There is an increase of thermal efficiency of current nuclear power plants from 30–35% to approximately 45–50%.
Figure 1 shows the difference in the operating conditions of current generation reactor systems in comparison to SCWRs. Compared to existing pressurized water reactors (PWRs), in SCWRs the target is to increase the coolant pressure from 10–16 MPa to about 25 MPa; the inlet temperature to about 350 °C, and the outlet temperature to about 625 °C [1].
In this paper, we presented a numerical analysis of the effect of different heat transfer correlations on the prediction of fuel and wall cladding temperatures in a supercritical water reactor. The neutronics process with temperature feedback effects, the heat transfer in the fuel rod and the thermalhydraulics in the core were simulated. Special attention was given to the thermalhydraulics, which uses a threepass core design with multiple heatup steps, where each step was simulated using an average channel. The first pass called “evaporator” is located in the center of the core. In this region, the moderator water flows downward in gaps between assembly boxes and inside the moderator tubes. The moderator water, heatedup through its path downward to the lower plenum, is mixed with the coolant coming from the downcomer reaching an inlet temperature of around 583 K. The evaporator heats the coolant up to 663 K, flowing upward and around the fuel rods, resulting in an outlet temperature 5 K higher than the pseudocritical temperature of 557.7 K at a pressure of 25 MPa. The second pass, called “superheater”, with downward flow, heats the coolant up to 706 K. After a second mixing in an outer mixing plenum below the core, the coolant will finally be heated up to 803 K with an upward flow in a second superheater (the third pass) located at the core periphery. A transient onedimensional radial conduction model was applied in the fuel rod for each cell in the axial coordinate. Energy balances for the coolant have been implemented using a steady state and a onedimensional model for the axial coordinate. Fuel lattice neutronics calculations were performed with the HELIOS2 code and the reactivity coefficients were used to evaluate the reactivity effects due to changes in the fuel temperature and in the supercritical water density for 177 energy groups. Due to the strong variation of coolant density through the core, five densities were considered. This safety parameter is calculated in order to evaluate the variation of the reactivity due to the Doppler effect, as a function of the fuel temperature, which is related to the resonances broadening when the fuel temperature increases. The coupling of neutronics with the heat transfer in the fuel rod, and the thermalhydraulics is presented, and numerical experiments due to changes in the mass flow rate were accomplished in this study. Effects on fuel temperature predictions with improved heat transfer correlations and classical heat transfer correlations were also compared.
2 Supercritical fluids
The behavior of liquid and gas density with pressure and temperature is illustrated in Figure 2. When the pressure and temperatures are low, there is a significant density difference between the liquid and the gas states. Near the critical point, the density difference between the liquid and gas is small, and above the critical point, the densities of the liquid and the gas have become equal.
The heat transfer process, at critical and supercritical pressures, is influenced by the significant changes in thermophysical properties, as is observed in Figure 3 for specific heat, thermal conductivity, and density obtained from thermal properties taken from [3]. The most significant thermophysical property variations occur near the critical and pseudocritical points. For example, the specific heat of water has a maximum value at the critical point. The exact temperature that corresponds to the specific heat peak at pressures above the critical pressure is known as the pseudocritical temperature [4].
Fig. 3
Behavior of the specific heat (C_{p}), thermal conductivity (k) and density (ρ), as a function of the temperature at 25 MPa. 
3 Supercritical water heattransfer correlation
The practical prediction methods for heat transfer at supercritical pressures are presented in [1,4]. The supercritical water heat transfer correlations applied in this work are shown in Table 1. Dimensionless numbers used in Table 1 are given by: (1) which is the Nussel number. Here H_{∞} is the heat transfer coefficient, k is the thermal conductivity and D_{H} is the hydraulic characteristic length. The subscript b means that bulkfluid temperature is used to calculate the thermophysical properties. These properties can also be calculated with the wall temperature, which will be specified with a subscript w. The Reynolds number is defined by:
(2) where G is the mass flux and μ is the viscosity. The Prandtl number is defined as:
(3) where Cp is the specific heat. The heat transfer coefficient is used in the boundary condition given below in equation (6), and H_{∞} represents the heat transfer from the wall to the coolant. McAdams [6] proposed the use of the DittusBoelter correlation for forcedconvective heat transfer in turbulent flows at subcritical pressures. The only difference between the DittusBoelter and McAdams correlations is that the latter has a larger coefficient. According to Schnurr et al. [10], it agrees with experimental data. However, it was noted that the correlation might produce unrealistic temperature results near the critical and pseudocritical points, due to it being very sensitive to variations in the thermophysical properties.
Bishop et al. [7] conducted experiments in supercritical water flowing upward inside bare tubes and annuli, within the following range of operating parameters: P = 22.8–27.6 MPa, T_{b} = 282–527 °C, m = 651–3662 kg/m^{2} s and q = 0.31–3.46 MW/m^{2}. Their data for heat transfer in tubes were generalized with a fit of ±15%. This correlation uses a crosssectional averaged Prandtl number and the final term in the correlation (1 + 2.4 D/x) accounts for the entranceregion effect. Bishop et al.'s correlation was modified and used without the entranceregion term, because this term depends significantly on the particular design of the inlet of the bare test section.
Swenson et al. [8] have suggested a correlation in which thermophysical properties are mainly based on a wall temperature, as they found that conventional correlations, which use a bulkfluid temperature as a basis for calculating the majority of thermophysical properties, did not work as well.
A dimensional analysis was performed by Mokry et al. [9] in order to obtain a general empirical form of correlation for the heat transfer calculations, and as a result of the experimental data analysis, two correlations for the heat transfer coefficient at supercritical water conditions were obtained.
In the core layout of the SCWR under study, water, as the working fluid, is guided three times through the core (twice up and once down). This design is called the threepass core concept. The first pass, called the evaporator, is situated in the center of the core. In this region, the moderator water flows downward in gaps between assembly boxes and inside the moderator tubes. The moderator water, in its downward path to the lower plenum is heated up, and is mixed with the coolant (1200 kg/s as inlet mass flow) which comes from the downcomer, thereby reaching an inlet temperature of around 583 K. The evaporator heats the coolant up to 663 K, flowing upward around the fuel rods, resulting in an outlet temperature 5 K higher than the pseudocritical temperature of 557.7 K at a pressure of 25 MPa. An inner steam plenum above the core eliminates hot streaks. The second pass, called superheater, with a downward flow, heats the coolant up to 706 K. After a second mixing in an outer mixing plenum below the core, the coolant is finally heated up to 803K with an upward flow in a second superheater located at the core periphery, known as the third pass. Each pass, the evaporator and both superheaters, is built of 52 fuel assembly clusters as shown in Figure 4 [11]. Therefore the complete reactor core is composed of 156 assembly clusters.
The fuel assembly design is taken from the European high performance light water reactor (HPLWR) concept. The 7 × 7 square arrangement design, with 40 fuel rods distributed in dual rows, and a single water tube replacing 9 fuel rods was used [12]. The fuel rods and the water tube are housed within the assembly box and grouped in a cluster of 9 assemblies, in a 3 × 3 arrangement with similar dimensions to a PWR assembly. As found in the PWR, control rods are inserted from the core top into 5 of the 9 water tubes of a cluster (Fig. 4b). The structural material for cladding, assembly boxes and water tubes is stainless steel. The main reactor parameters are presented in Table 2.
Supercritical water heattransfer correlations (HTCs).
Fig. 4
(a) Arrangement of evaporator, superheater 1, superheater 2, and assembly clusters in the core, (b) assembly cluster with water tubes and control rods. 
4 Implementation of the heat transfer correlations
In order to analyze the effect of different heat transfer correlations on the prediction of the wall temperature of the fuel rods, the SCWR numerical code developed by BarragánMartínez [15] was applied using the HTCs shown in Table 1. The numerical model of the heat transfer processes in the fuel element of the HPLWR was obtained using the numerical model of typical reactors [16]. The supercritical water reactor is integrated of cylindrical fuel elements which contain ceramic pellets inside the cladding.
Then, the effect of heat transfer correlations on the fuel temperature prediction of SCWRs was conducted with numerical experiments.
4.1 Fuel heat transfer model
A detailed multinode fuel pin model was developed for this study. The fuel heat transfer formulation is based on the following fundamental assumptions:

axissymmetric radial heat transfer,

the heat conduction in the axial direction is negligible,

the volumetric heat rate generation in the fuel is uniform in each radial node, and

storage of heat in the fuel cladding and gap is negligible.
Under these assumptions, the transient temperature distribution in the fuel pin, and the initial and boundary conditions are given in the following conditions: (4) (5) (6) (7)
In equation (4) q‴(t) = 0, for r_{f} ≤ r_{cl}. In these equations, r is the cylindrical radial coordinate, r_{0}, r_{f} and r_{cl} are the centroid, fuel and clad radius, respectively, q‴(t) = P(t)/V_{f} at each axial node, where P is the neutronics power, T_{m} is the moderator temperature, and H_{∞} is the convective heat transfer coefficient.
The differential equations described previously are transformed into discrete equations using the control volume formulation technique in an implicit form [17]. The control volume formulation enables the equations for fuel, gap, and cladding to be written as a single set of algebraic equations for the sweep in the radial direction: (8) where , and are unknowns, a_{j}, b_{j}, c_{j} and d_{j} are coefficients, which are computed at the time t. When these equations are put into a matrix form, the coefficient matrix is tridiagonal. The solution procedure for the tridiagonal system is the Thomas algorithm, which is the most efficient algorithm for this type of matrices. The coefficients a_{j}, b_{j}, and c_{j} are dependent on thermophysical properties, i.e., thermal conductivity, density and specific heat; and since they are function of , at least one iteration is needed.
4.2 Thermalhydraulics model
The basic equations for describing the thermalhydraulics behavior in the three representative heated channels (one channel for each pass core) assuming the supercritical fluid is a single phase fluid, are presented as following. Incompressible flow was also considered in this study, i.e., the mass flux (G) is a constant. Under this consideration, the energy equation at steady state is shown as follows: (9) where T is the temperature, f is the friction factor, P_{H} is the heated perimeter, A_{f} is the flow area. The heat transfer from the wall to the coolant is obtained with Newton's law of cooling:
(10) The temperature in each node of the channel is obtained numerically as:
(11) where Δz is the node length and i is the node number.
4.3 Reactor power model
The reactor power is given by (12) where F(z) is the axial power factor, P_{0} is nominal power and n(t) is the normalized neutron flux, which is calculated by using a point reactor kinetics model with six groups of delayed neutrons:
(13) (14) where C_{i} is a delayed neutron concentration of the ith precursor group normalized with the steadystate neutron density, ρ is the net reactivity, β is the neutron delay fraction, Λ is the neutron generation time and β_{i} is the portion of neutrons generated by the ith group. The initial conditions are given by n(0) = n_{0} and c_{i}(0) = β_{i}n_{0}/Λλ_{i}. The parameters of the kinetics model are presented in Table 3.
The net reactivity in this work includes three main components: Doppler effects due to fuel temperature, coolant density, and reactor control rods.
The kinetics point equations are stiff in the coefficients because they differ in several orders of magnitude. The implicit variable integration method was used to solve equation (13), and the Euler method in an explicit form was used to solve the delayed precursor concentration given by equation (14).
The reactivity coefficient due to variations in fuel temperature was studied for the square fuel assembly design proposed by [18]. The calculations were done for the fuel assembly model along the active core height. Due to the strong variation of coolant density through the axial direction of the core, five densities: 0.74, 0.45, 0.31, 0.17 and 0.09 g/cm^{3} were considered. This safety parameter is calculated in order to evaluate the variation of the reactivity due to the Doppler Effect, as a function of the fuel temperature, which is related to the resonances broadening when the temperature increases. The values of the reactivity as a function of the coolant density and fuel temperature are presented in Figure 5. The values of the infinite multiplication factor obtained with HELIOS2 for 177 energy groups were used to determine the reactivity [19].
Fig. 5
Reactivity coefficients obtained with HELIOS2 for 177 energy groups at different densities. 
4.4 Representative SCWR nodalization
The fuel rod temperature distribution was obtained for the radial nodes at each of the twenty one thermalhydraulics axial nodes in the core. The arrangement of the computational nodes of the thermalhydraulics model is illustrated in Figure 6.
Figure 7 shows the grid used in calculations. Half control volume near the boundary, radial nodes 1, 2, 3, 4, and 5 for the fuel; radial node 6 was used for the gap; radial nodes 7 and 8 for the clad. Radial nodes 1 and 8 were used for the boundary condition.
Fig. 6
Arrangement of the computational nodes in the thermalhydraulics core model of the SCWR. 
Fig. 7
Arrangement of the computational cells of fuel, gap, and clad. 
5 Numerical experiments
Each channel in the core was based on an hydraulic unit cell whose parameters are: P_{H} = 0.025 m, D_{H} = 0.054 m, and A_{f} = 0.34 m^{2}. The parameters of the fuel element are: r_{f} = 5.207 × 10^{−3} m for the fuel, r_{g} = 5.321 × 10^{−3} m for the gap, and r_{cl} = 6.134 × 10^{−3} m for the clad. The active height of the fuel cell (4.2 m) was divided into 21 equidistant axial nodes (Δz = 0.2 m). The distribution axial of power for each channel was imposed with the idea that the heat flux is not uniform. The thermal physical properties used were taken from Wagner and Kretzschmar [3]. 73, 48 and 35 assembly clusters for Channel 1, Channel 2 and Channel 3, respectively, were used in the simulation, in order to reach a better power distribution within the core.
Figure 8 presents the results for Channel 1, showing the wall temperature behavior for different correlations presented in Table 1. It should be noted that the last node temperature (at 4 m) is practically the same, and the trend is very similar for all the correlations, except for a short zone where the Swenson correlation yields a lower temperature while Mokry's correlations (both preliminary and final) yield a higher temperature, the same was noted for the Bishop's correlations (with and without ERE).
In Figure 9 the results for Channel 2 are presented, showing the wall temperature behavior for the correlations in Table 1. Similar results were obtained, however contrary to what was observed in Channel 1, the Swenson's correlation yields slightly higher temperatures along the entire channel meanwhile the Bishop's (with and without ERE) and Mokry's correlations yield slightly lower temperatures along the entire channel.
Figure 10 presents the results for Channel 3, showing the wall temperature behavior for the correlations presented in Table 1. In this case, the trends that most resemble each other are presented. Again, the Swenson's correlation deviates the most, yielding slightly higher temperatures than other correlations.
In Figure 11 the results along the three channels are presented. It should be noted that Swenson's correlation is the one with greater deviation from DittusBoelter's correlation, with a difference of 10 K in Channel 1.
There is a wall temperature reduction at the end of each channel; especially for Channel 2, and this is due to the axial distribution of thermal power which has a minimum in this bottom core zone. This is an undesired result of the three pass core concept.
Another finding in this numerical analysis was that Swenson's correlation gave the most conservative predictions, in terms of safety, because higher temperatures are calculated due to the use of the wall temperature for the Re and Pr calculations, while the other correlations use the bulk temperature.
Fig. 8
Simulation results for Channel 1 showing the wall temperature behavior for different HTCs. 
Fig. 9
Simulation results for Channel 2 showing the wall temperature behavior for different HTCs. 
Fig. 10
Simulation results for Channel 3 showing the Wall Temperature behavior for different HTCs. 
Fig. 11
Simulation results showing the wall temperature behavior across the three channels for different HTCs. 
6 Conclusions
The correlation, which agrees most with DittusBoelter, is McAdams. The only difference in the equation is the value of the coefficient. Bishop's correlations, with and without EntranceRegion Effect (ERE) have little differences among them in the prediction of the wall temperatures, meaning that, for this simulation the ERE is not important; predictions compared to the DittusBoelter correlation are a little higher in the first channel and slightly lower in Channels 2 and 3. With preliminary and final Mokrýs correlations, higher temperature predictions were found in Channel 1, but were very similar to DittusBoelter in Channels 2 and 3. Swenson's correlation showed the most deviated results, yielding lower temperatures in the first channel and higher in Channels 2 and 3.
Swenson's correlation uses the wall temperature for calculating the Re and Pr numbers, while the others used the bulk temperature and we found the greatest differences compared to other HTCs. For this reason Swenson's correlation could be very useful in order to find the most conservative results for Channel 3, where high wall temperatures could affect the fuel rod integrity.
Acknowledgments
Special thanks to the National Council for Sciences and Technology (CONACYT) for the scholarship provided to the Master Student Erick G. EspinosaMartinez, and to the National Autonomous University of Mexico for the PAPIIT IN113213 project funds.
References
 I. Pioro, R. Duffey, Heat Transfer and Hydraulic Resistance at Supercritical Pressures in Power Engineering Applications (ASME Press, New York, 2007) [CrossRef] [Google Scholar]
 H. Thind, Heattransfer analysis of doublepipe heat exchangers for indirectcycle SCW NPP, Master Thesis, University of Ontario Institute of Technology, Ontario, 2012 [Google Scholar]
 W. Wagner, H.J. Kretzschmar, International Steam Tables. Properties of Water and Steam Based on the Industrial Formulation IAPWSIF97 (Springer, Berlin, 2008), 2nd ed. [Google Scholar]
 I.L. Pioro, H.F. Khartabil, R.B. Duffey, Heat transfer to supercritical fluids flowing in channels – empirical correlations (survey), Nucl. Eng. Des. 230, 69 (2004) [CrossRef] [Google Scholar]
 F.W. Dittus, L.M. Boelter, Heat transfer in automobile radiators of the tubular type, Int. Commun. Heat Mass Transf. 12, 3 (1930) [CrossRef] [Google Scholar]
 W. McAdams, Heat Transmission (McGrawHill, New York, 1942), 2nd ed. [Google Scholar]
 A.A. Bishop, R.O. Sandberg, L.S. Tong, High Temperature Supercritical Pressure Water Loop: Part IV, Forced Convection Heat Transfer to Water at NearCritical Temperatures and SuperCritical Pressures (Westinghouse Electric Corporation, Pittsburgh, Pennsylvania, 1964) [Google Scholar]
 H.S. Swenson, J.R. Carver, C.R. Kakarala, Heat transfer to supercritical water in smoothbore tubes, J. Heat Transf. Trans. ASME Series C 87, 477 (1965) [CrossRef] [Google Scholar]
 S. Mokry, A. Farah, K. King, S. Gupta, I. Pioro, P. Kirillov, Development of supercritical water heattransfer correlation for vertical bare tubes, Nucl. Eng. Des. 241, 1126 (2011) [CrossRef] [Google Scholar]
 N.M. Schnurr, V.S. Sastry, A.B. Shapiro, A numerical analysis of heat transfer to fluids near the thermodynamic critical point including the thermal entrance region, J. Heat Transf. Trans. ASME 98, 609 (1976) [CrossRef] [Google Scholar]
 T. Schulenberg, J. Starflinger, Core design concepts for high performance light water reactors, Nucl. Eng. Technol. 39, 249 (2007) [CrossRef] [Google Scholar]
 J. Hofmeister, C. Waata, J. Starflinger, T. Schulenberg, E. Laurien, Fuel assembly design study for a reactor with supercritical water, Nucl. Eng. Des. 237, 1513 (2007) [CrossRef] [Google Scholar]
 T. Reiss, S. Fehér, S. Czifrus, Coupled neutronics and thermohydraulics calculations with burnup for HPLWRs, Prog. Nucl. Energy 50, 52 (2008) [CrossRef] [Google Scholar]
 T. Schulenberg, J. Starflinger, High Performance Light Water Reactor. Design and Analyses (KIT Scientific Publishing, Germany, 2012) [Google Scholar]
 A.M. BarragánMartínez, Diseño neutrónico y termohidráulico de un reactor nuclear enfriado con agua supercrítica, PhD Thesis, Universidad Nacional Autónoma de México, Mexico City, 2013 [Google Scholar]
 G. EspinosaParedes, E.G. EspinosaMartínez, Fuel rod model based on NonFourier heat conduction equation, Ann. Nucl. Energy 36, 680 (2009) [CrossRef] [Google Scholar]
 S.V. Patankar, Numerical Heat Transfer and Fluid Flow (McGrawHill, New York, 1980) [Google Scholar]
 A.M. BarragánMartínez, C. MartindelCampo, J.L. François, G. EspinosaParedes, MCNPX and HELIOS2 comparison for the neutronics calculations of a Supercritical Water Reactor HPLWR, Ann. Nucl. Energy 51, 181 (2013) [CrossRef] [Google Scholar]
 J. Duderstadt, L. Hamilton, Nuclear Reactor Analysis (John Wiley & Sons, United States of America, 1976) [Google Scholar]
Cite this article as: ErickGilberto EspinosaMartínez, Cecilia MartindelCampo, JuanLuis François, Gilberto EspinosaParedes, Effect of heat transfer correlations on the fuel temperature prediction of SCWRs, EPJ Nuclear Sci. Technol. 2, 35 (2016)
All Tables
All Figures
Fig. 1
Operating conditions of current nuclear reactors and SCWRs [1]. 

In the text 
Fig. 2
Schematic behavior of liquid and gas density with pressure and temperature [2]. 

In the text 
Fig. 3
Behavior of the specific heat (C_{p}), thermal conductivity (k) and density (ρ), as a function of the temperature at 25 MPa. 

In the text 
Fig. 4
(a) Arrangement of evaporator, superheater 1, superheater 2, and assembly clusters in the core, (b) assembly cluster with water tubes and control rods. 

In the text 
Fig. 5
Reactivity coefficients obtained with HELIOS2 for 177 energy groups at different densities. 

In the text 
Fig. 6
Arrangement of the computational nodes in the thermalhydraulics core model of the SCWR. 

In the text 
Fig. 7
Arrangement of the computational cells of fuel, gap, and clad. 

In the text 
Fig. 8
Simulation results for Channel 1 showing the wall temperature behavior for different HTCs. 

In the text 
Fig. 9
Simulation results for Channel 2 showing the wall temperature behavior for different HTCs. 

In the text 
Fig. 10
Simulation results for Channel 3 showing the Wall Temperature behavior for different HTCs. 

In the text 
Fig. 11
Simulation results showing the wall temperature behavior across the three channels for different HTCs. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.