Preliminary Development of Thermal Power Calculation Code H-Power for a Supercritical Water Reactor

SCWR (Supercritical Water Reactor) is one of the promising Generation IV nuclear systems, which has higher thermal power efficiency than current pressurized water reactor. It is necessary to perform the thermal equilibrium and thermal power calculation for the conceptual design and furthermonitoring and calibration of the SCWR.One visual software namedH-Power was developed to calculate thermal power and its uncertainty of SCWR, in which the advanced IAPWS-IF97 industrial formulation was used to calculate the thermodynamic properties of water and steam. The ISO-5167-4: 2003 standard was incorporated in the code as the basis of orifice plate to compute the flow rate. New heat balance model and uncertainty estimate have also been included in the code. In order to validate H-Power, an assessment was carried out by using data published by US and Qinshan Phase II. The results showed that H-Power was able to estimate the thermal power of SCWR.


Introduction
SCWR (Supercritical Water Reactor) is one of the promising Generation IV nuclear systems, which has higher thermal power efficiency than current pressurized water reactor.Compared to the currently running water coolant reactors, SCWR has the higher thermal efficiency and safety, which makes it a more promising advanced nuclear energy system.The concept of SCWR was originally put forward by Westinghouse and GE (General Electric) in the 1950s and preliminarily studied by the United States and the former Soviet Union from 1950s to 1960s.In the 1990s, Dobashi et al. [1] proposed the concept of SCWR again and made a further development in this area.
For reactor operating more safely, stably, and economically, the accurate calibration for real thermal power has important significance in reactor design and analysis.It is necessary to develop a special thermal power calculation code for SCWR, while there is no available program published in the world.This paper preliminarily developed visual software named H-Power to calculate SCWR thermal power and its uncertainty based on heat balance method, which has been applied to current PWRs [2].
According to the various conceptual SCWR designs [3][4][5], H-Power divided SCWR into two groups: one is the single-loop type and the other is the double-loop type.The validation of H-Power consists of two parts: for the singleloop SCWR, the data published by Jacopo [6] was used, while for the double-loop SCWR, the values of Qinshan Phase II [7] were used; though it is not the SCWR, we just use it for preliminary validation since the exact detail design data was rarely published.
IAPWS is an international nonprofit association concerned with the properties of water and steam, particularly thermodynamic properties.In 1997, IAPWS adopted the "IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" for industrial use [8], called IAPWS-IF97 for short, which replaced the previous industrial formulation IFC-67 [9] published in 1967.
We have done lots of researches and work in comparing IAPWS-IF97 with IFC-67 in several aspects.The result showed that IAPWS-IF97 improves significantly in boundary consistency and accuracy and calculation speed [10], which makes it more extensive use in power industry nowadays.
H-Power adopted IAPWS-IF97 to calculate thermodynamic properties (including specific volume, specific

Theoretical Modeling
The modules of H-Power are given in Figure 1, where READDT is the reading data subroutine, NUSOLV is numerical calculation subroutine, THPROP is the subroutine which calculates the thermodynamic properties of water and steam, FLOWRT is the subroutine for flow rate calculation, ERCALC subroutine is for the uncertainty estimation, and OUTPUT is the output data subroutine.
H-Power provides a visual interface for users and has high module independence, which makes it a practical code in SCWR analysis.

Subroutine Flow Chart.
The range of IAPWS-IF97 is divided into five subareas; each of them has a basic equation.To make the program more concise and readable, this subroutine used the compiling method of modularization; that is, separate districts of thermodynamics equation were written into separate subroutine.
We can directly obtain the thermophysical properties from pressure and temperature values using the basic equations except region 3. Subregion 3 does not have the formula which can calculate other various thermal parameters directly through the temperature and pressure; it can only rely on iteration calculation.In this subroutine, the first pressure calculation uses temperature and a small specific volume which was chosen as initial value.If the difference between the calculated value of pressure and user input values is positive, the value of specific volume increases incremental quantity ratio by half and continues the iteration until the pressure difference is less than condition of convergence.The iteration is over and other properties can be calculated.
Through this kind of iteration and close proximity to the input value by increasing incremental quantity in halving bisection method, not only the disadvantages of the increased specific volume increment are too small to calculate quickly was avoided, but also the accuracy of calculation can be ensured.
The subroutine of water and steam thermodynamic properties has been totally developed according to the flow chart shown in Figure 2.

Subroutine Calculation Verification.
In order to do the verification, we list the calculation results together with those given by IAPWS-IF97 in Table 1.We can see that the present results are in good agreement with those given by IAPWS-IF97, so that the thermodynamic properties calculated by H-Power are valid.

Flow Rate Calculation Subroutine.
The main feedwater flow rate of the steam generator can be measured by orifice plate, which is widely used in measurement of fluid flow.This part is based on the ISO 5167 standard in [11].
As the schematic diagram shown in Figure 3, upstream fluid pressure  1 is obtained by (1) taking pressure loss Δ caused by the existence of orifice plate into account: H-Power adopted the following equation for flow rate measurement: where   is mass flow rate,  is ratio of  (diameter of orifice at operating conditions) and  (diameter of internal pipe at operating conditions),  is expansibility factor (for compressible fluid,  = 1), Δ is differential pressure,   is the density of upstream, and  is discharge coefficient determined by Reader-Harris/Gallagher equation as follows: Science and Technology of Nuclear Installations 3

Error and quit
Error and quit If  < 71.12 mm, (4) should add the following term: ) . ( As we can see from the above equation, discharge coefficient  depends on the Reynolds number Re  , which is dependent on   in turn, so the calculation of mass flow rate   has to be iterated. Diameters in the formula for calculating should be corrected due to the difference of temperature between working condition and measurement.If there is drain hole (diameter is   ) on the plate, the  should be corrected by the following equation: where the   have already been corrected by temperature as follows: where  0 is the measurement value under the standard temperature  0 ,  is the coefficient of linear expansion, and  is the operating temperature.Re  is the Reynolds number calculated by  (diameter of internal pipe at operating conditions), dimensionless parameters presenting the ratio of the inertia force and viscous force of upstream, calculated by the following equation: where  is viscosity of the fluid at the working conditions, which is obtained by the thermodynamic properties subroutine.
1 (=  1 /) is the ratio of the distance from the upstream tapping to the upstream face of the orifice plate and the pipe diameter, and the value selection of its conditions is as follows: 2 (=   2 /) is the ratio of the distance from the downstream tapping to the downstream face of the orifice plate and the pipe diameter, and the value selection of its conditions is as follows: The uncertainty of the feedwater flow measured and calculated by orifice plate and relevant equations will be given in Section 3.2.The computer verification of this subroutine will be integrated in Section 4.2.

Single-Loop SCWR Program.
In single-loop type SCWR, the coolant is supercritical water which will be operated above the critical point of water (  = 647.096K,   = 22.064 MPa).The coolant remains single-phase when operated above the critical point, which is a stupendous advantage for it will eliminate coolant boiling throughout the system.The no need for recirculation results in the elimination of pressurizer, jet pumps, steam generators, and steam separators.This direct cycle type is greatly simplified as shown in Figure 4.
The thermal power in this type can be calculated with thermal balance using measured values, including inlet and outlet temperatures, working pressure, flow rate of coolant, and the practical speed of pump.The thermal power () was obtained by the following equation: where  (kg/s) is the flow rate of the coolant in the loop, Ω (tr/min) is the practical speed of pump, Ω 0 (tr/min) is the ideal speed of pump, and  ℎ −   (kj/kg) is the difference between the enthalpy at the hot leg and the cold leg.The enthalpies are calculated using THPROP subroutine.

Double-Loop SCWR Program
. Some conceptual designs of double-loop SCWR have been put forward, like CANDU-SCWR with a steam generator proposed by Dr. William Fatoux [12] and a pressurized water reactor cooled with supercritical water in the primary loop proposed by Vogt et al. [13].Double-loop SCWR drives a turbine by the steams from the steam generator in the second loop indirectly as shown in Figure 5. Heat balance method in this type of SCWR is based on the steam generator approximation enthalpy balance.H-Power calculate the thermal power of reactor by the theory of heat balance, which obtains thermal power (of generator) in second loop from primary loop through the measurement values of temperature, pressure, and flow rate in the second loop, counting the internal heat loss in at the same time.The advantages of the working conditions of second loop for easy measurement makes it an important method of obtaining thermal power of the double-loop SCWR.
Figure 6 shows the heat in and out the system (expressed by arrows), which can be described by the following heat balance equation: where the   is the core thermal power (MW),  SG  is the thermal power obtained from the primary loop in the second loop in th steam generator, and  ΔPr is the thermal power coolant system obtained from other heat sources, (MW), which can be estimated from the design of reactor.
The main physical process occurring in the steam generator is as follows: unsaturated water with the mass flow rate   and specific enthalpy   flows into the steam generator through the  tube, exchanges heat with the coolant in primary loop, and mostly becomes wet steam with the mass flow  V and specific enthalpy  V ; after the increase of endothermic enthalpy, the others are the saturated blowdown water with mass flow rate   and specific enthalpy   .The thermal power in ith steam generator can be calculated with following equation: where  V is the wet steam enthalpy in steam generator export,   is the blowdown enthalpy,   is the flow rate of blowdown in second loop,   is the feedwater enthalpy in second loop, and   is the flow rate of feedwater in second loop.The feed water satisfies the law of conservation of mass as follows: Simplify ( 13) and ( 14): where feedwater enthalpy   and blowdown enthalpy   are calculated by the thermodynamic properties subroutine, flow rate of feedwater   can be obtained from flow rate calculation subroutine, and the flow rate of blowdown   is a measured value.The wet steam always consists of saturated steam and saturated water, which is determined by vapour fraction , so the enthalpy of wet steam is characterized by the following equation: where the  V is the enthalpy of saturated steam and   is the enthalpy of saturated water; both are calculated by the thermodynamic properties subroutine.

Uncertainty Analysis of Double-Loop SCWR
The relative uncertainty of thermal power given in ( 17) is obtained by the uncertainty transmission formula: where the Δ SG / SG is the uncertainty of thermal power in steam generation in each second loop and Δ( ΔPr )/ ΔPr is the uncertainty of thermal power of coolant system in primary loop: where [Δ SG /Δ]  is the max uncertainty of each steam generation.
Assuming that We can obtain the following equation: 3.1.The Calculation of Δ  /  .The following equation can be obtained by error propagation formula: Assume that the uncertainty of saturated water in the wet steam is "1" as in ( 22) and inaccuracy of measurement should be in this range: Equation ( 22) can be deduced by the above formula: Δ V / V is the uncertainty of saturated steam enthalpy, mainly caused by the uncertainty of pressure measurement and water and steam thermodynamic properties as follows: and the uncertainty caused by pressure measurement can be calculated by following equation: where Δ V is uncertainty caused by pressure measurement, root mean square of the uncertainty of pressure transmitter, data acquisition system, and calibration, determined by the design of pressure measurement: where   ,   ,   is given in IAPWS-IF97, and  = 0.1 V / *  =  * /,  * = 1MPa,  * = 540 K, and is the uncertainty caused by the calculation of water and steam thermodynamic properties, where Δ V can be selected in the corresponding section of IAPWS-IF97.

The Calculation of Δ𝐻
/  .The uncertainty of saturated water enthalpy Δ  /  and feedwater enthalpy Δ  /  can be estimated by the method the same as Δ V / V mentioned above.
The relative uncertainty of mass flow rate   /  can be described by the following equation: The uncertainty of discharge coefficient is given by If  > 0.5 and Re  < 10000, the above values should add the following relative uncertainty: The uncertainty of expansibility factor  is 3.5(Δ/ 1 )%, for compressible fluid,  = 0.
The uncertainty of internal diameter of the pipe line can be estimated by the following equation: where 0.05%×  is the error of the measuring instrument; "2" is an extension coefficient of the confidence interval of 95%; √ 3 is the constant of rectangular distribution.The volume elasticity coefficient of water is quite large determining that its compressibility is very small, so the influence of pressure on the density can be ignored.The uncertainty of feed water density / is determined by the uncertainty of temperature measurement and the calculation of water properties as follows:  The calculation procedures and method of the parameters are the same as the Δ V / V mentioned above.
The relative uncertainty of Δ( ΔPr )/ ΔPr is estimated by experience values.

Validation of H-Power Code
4.1.Single-Loop Type Program.We have done the verification of single-loop program in H-Power using the main parameters design by Jacopo Buongiorno as in Table 2 and assuming the ideal and practical pump speed according to the level of pump speed at present.
The calculation result of thermal power 3564.3396MW given in Table 3 shows a high consistency with the value 3575MW designed by US in Table 3.

Double-Loop SCWR Program.
Due to the lack of design data of SCWR, the validation of this part has been done by using the values of Qinshan Phase II [13] in Table 4 though it is not SCWR.Table 5 shows the comparison between the results of H-Power and the values of Qinshan Phase II.
Although the main calculation result trend of SCWR is same as CNP650 of Qinshan Phase II, there is many detail differences in thermodynamic properties of water and steam and uncertainties calculation.The inconsistencies in Table 5 also come from the different input parameters in measuring instrument.
It is a valid way to verify the H-Power while the main calculation process of both PWR and SCWR is the same and there is no particular data of SCWR fitting for the code.The result in Table 5 shows that H-Power runs well and has a valid trend in thermal power and its relevance computation.So it can be used in SCWR design and calculation.

Conclusions
The thermal power calculation code H-Power for SCWR was developed based on heat balance method, which can be easily applied in designs and operation.In this paper thermodynamic properties calculation subroutine based on IAPWS-IF97 has been verified, and the uncertainty analysis has been proposed.The validation of H-Power has been carried out in two parts: single-loop SCWR by the data given by Jacopo Buongiorno and double-loop SCWR by the values of Qinshan Phase II.Although Qinshan Phase II is a PWR, the result can show that H-Power is valid because the main calculation process of both PWR and SCWR is the same.In the future, the validation will be carried out using the values of SCWR when published.

Figure 4 :
Figure 4: Sketch of SCWR with one loop.

Table 1 :
Comparison of calculation results between H-Power routine and IAPWS-IF97.

Table 2 :
The main parameters reference design of the US SCWR.

Table 3 :
The input and output values calculated by H-Power.

Table 4 :
Input parameters of H-Power using the values of Qinshan Phase II.

Table 5 :
Comparison between the results of the values of Qinshan Phase II and H-Power.