# Analysis on Steady-State Operation and Heat Loss of Chinese Integrated Pressurized Water Reactor

Authors: Zhang Fan; Lu Dao-Gang; Sui Dan-Ting; Guo Chao; Yuan Bo
Journal: Science and Technology of Nuclear Installations
 (2014)

---

## Abstract

Chinese integrated pressurized water reactor (CIPWR) has compact configuration and high inherent safety, which is appropriate for nuclear power plants of small and medium scale. Heat balance model has been adopted widely in thermal power calibration of PWRs because of its advantage of accuracy. In this paper, a package based on FORTRAN language is developed and added into RELAP5 to calculate the heat loss value needed in heat balance analysis. The steady-state operation of CIPWR is modelled correctly by RELAP5. The heat loss of CIPWR is calculated by the package, and the comparison of the main values of parameters needed in the heat loss calculation between RELAP5 and the package has been done. It shows that the package has high calculation accuracy and can be applied in reactor design and monitoring.


---

## Body

## 1. Introduction

Reactor Excursion and Leak Analysis Program (RELAP5) is a reactor analysis program developed by Idaho National Engineering Laboratory (INEL), approved by Nuclear Regulatory Commission (NRC) for engineering review of thermal hydraulic calculation of PWRs [1]. RELAP5 has been generally applied to reactor design and safety analysis of PWRs, which is convenient for users and has high calculation efficiency. With the increasing study requirements, lots of routine packages are added to RELAP5 by users to meet their own study needs.To make reactor operate more safely and economically, the accurate calibration of real thermal power is important for design, analysis, and monitoring of the reactor, which has been applied in PWRs currently [2]. The heat balance model has been adopted widely in thermal power calibration with advantage of accuracy. The heat loss (or increase) from other heat sources (except reactor core) in primary loop is needed in the heat balance calculation, which has to be obtained from operating experience sometimes. A package based on FORTRAN language is developed and added into RELAP5 in this study, which makes it have the ability to calculate the heat loss in heat balance calculation.Integral pressurized water reactor (IPWR) is one of the hot issues in the field of small-size nuclear systems, which has high inherent safety. The structure of IPWR is more compact for the reason of the less use of pipes. The characteristic of small size makes IPWR more easily used in different conditions than other reactors, such as ships and remote places. There are various IPWR designs in the world with the same characteristic [3, 4]; that is, all primary-loop components are arranged inside the pressure vessel to eliminate or minimize the possibility of loss of coolant accident (LOCA). IRIS [5] reactor developed by USA, System-Integrated Modular Advanced Reactor (SMART) [6] developed by the Korea Atomic Energy Research Institute (KAERI), and CAREM-25 [7] developed by Argentina are the typical IPWR designs with coolant pumps. The multiapplication small light water reactor (MASLWR) [8] without any pumps uses natural circulation in both steady and transient operations. There are several different concepts of IPWR proposed by Nuclear Power Institute of China [9, 10]. One of these designs of IPWR with coolant pumps is our study object, which is called CIPWR in this paper.Our study aims at evaluating the performance of CIPWR in stable and transient operations. In this paper, RELAP5/MODE3.4 is used to model the stable condition of CIPWR and calculate the heat loss needed in heat balance analysis with the package we developed. The series of transient accidents of CIPWR will be studied in the future.


## 2. Theoretical Modelling

### 2.1. Heat Loss Calculation Package

The package added to RELAP5 is used to calculate the heat loss in heat balance analysis. All the input values in the package can be obtained from the calculation result of steady condition by RELAP5, such as the value of temperature, pressure, and vapour fraction. Figure1 shows the flow chart of the package.Figure 1
Structure of the heat loss package.


#### 2.1.1. The Principle of Heat Balance

Heat balance method is based on the steam generator approximation enthalpy balance [11]; that is, the enthalpy increase of feedwater is the main part of the thermal power. The thermal power can be calculated through the enthalpy difference between the steam and feedwater in the second loop, counting the internal heat loss at the same time.As is shown in Figure2, there are several ways for CIPWR primary system to obtain or lose heat (for other PWRs, the pressurizer electrical heating should be taken into consideration). The heat increase consists of core heat production, the heat brought in by charging stream, and heat production during the pump operation. The heat decrease consists of heat dissipation of components and the heat taken by letdown stream. The feedwater in the second loop exchanges heat with primary loop by U tubes. The heat balance principle of the CIPWR can be described by
(1)WR=WSG-WΔPr⁡,
where the WR (MW) is the core thermal power, which can be obtained from the RELAP5 calculation result and WSG (MW) is the steam generator power obtained from the primary loop, which can be obtained by the package. WΔPr⁡ (MW) is the heat loss, which is the sum of thermal power coolant obtained from other heat sources except core and lost because of the heat dissipation of components and pipes. The heat loss WΔPr⁡ (MW) is the calculation result of the package.Figure 2
Heat balance principle diagram of CIPWR.In the second loop, unsaturated feedwater with the mass flow rateQe (kg/s) and specific enthalpy He flows through the heat transfer tube and convects with the first side of steam generator and turns to the steam with the specific enthalpy Hv (kJ/kg). Equation (2) can be used to calculate the steam generator power:
(2)WSG=(Hv-He)Qe,
where Hv is the steam enthalpy of steam in the generator export, He (kJ/kg) is the feedwater enthalpy, and Qe (kg/s) is the flow rate of feedwater.


#### 2.1.2. Thermodynamic Properties Subroutine

The enthalpy of water and steam included in the heat balance calculation can be computed by the thermodynamic properties subroutine according to IAPWS-IF97 [12] (IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam). IAPWS-IF97 replaced the previous industrial formulation IFC-67 [13] for industrial use for its advantage of calculation with speed and accuracy.The range of IAPWS-IF97 has five subareas, each of which has a basic equation that can be used to calculate the water and steam properties by derivation. As one of the thermophysical properties, enthalpy can be obtained directly from the pressure and temperature values by the basic equations except region 3. The feedwater in the entrance of steam generator is unsaturated water belonging to region 1, whose thermophysical properties can be computed by (3). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(3)g(p,T)RT=γ(π,τ)=∑i=134ni(7.1-π)Ii(τ-1.222)Ji,
whereπ=p/p* and τ=T*/T with p*=16.53 MPa, T*=1386 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients ni and exponents Ii and Ji are given in the IAPWS-IF97.The enthalpy of feedwater can be computed by(4)He=g-T(∂g∂T)p=τRT×∑i=134ni(7.1-π)IiJi(τ-1.222)Ji-1.The steam in the export of steam generator loop is wet saturated steam consisting of saturated steam and saturated water with different enthalpy. The total enthalpy of steam can be calculated according to the vapor fractionx as
(5)Hv=xHvs+(1-x)Hls,
where Hls (kJ/kg) is the saturated water enthalpy, which can be calculated by (4). Hvs is the saturated steam enthalpy belonging to region 2, which can be calculated by (6). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(6)g(p,T)RT=γ(π,τ)=ln⁡π+∑i=19nioτJio+∑i=143niπIi(τ-0.5)Ji.The enthalpy of saturated steamHvs (kJ/kg) can be computed by
(7)Hvs=g-T(∂g∂T)p=τRT×{1π+∑i=143niIiπIi-1(τ-0.5)Ji},
where π=p/p* and τ=T*/T with p*=1 MPa, T*=540 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients nio, Jio, ni, Ii, and Ji are given in the IAPWS-IF97.As shown in Table1, the validation of this subroutine has been done by comparing it with the value given in IAPWS-IF97.Table 1
Comparison of calculation results between package subroutine and IAPWS-IF97.
Region
Parameters
Given by
h (kJ kg−1)
Region 1
T = 300 K
IF97
115.331273
P = 3 MPa
Package
115.33127302
T = 300 K
IF97
184.142828
P = 80 MPa
Package
184.14282773
Region 2
T = 300 K
IF97
2549.91145
P = 0.0035 MPa
Package
2549.9114508
T = 2000 K
IF97
6571.22604
P = 30 MPa
Package
6571.2260386


### 2.2. Modelling of Steady Operation of CIPWR

The whole reactor system can be modelled by different modules of RELAP5, which contains hydrodynamic components, branches, heat structures, reactor kinetics, control system, and so on. Hydrodynamic components consist of the parts where the fluid flows through, such as pipes, annulus, valves, pumps, and time-dependent volumes. Heat structures simulate the heat transfer in the system by setting the thermal conductivity and heat capacity of materials. Control system can set the change with time of physical quantities and trips and also can be used in self-initialization. Reactor kinetics consist of multigroup neutron transport model and point reactor model, which calculate the power of reactor by setting the initial power and reactivity feedback coefficient of various influent factors.Figure3 shows the sketch of CIPWR with two loops with coolant pumps. Once-through steam generators (OTSG) are adopted for the advantage that the configuration is compact enough to enable their location in the pressure vessel. OTSGs are arranged above the reactor core for the better natural circulation in case of SBO accident. The upper pressure vessel space acts as pressurizer in CIPWR, which eliminates lots of pipes effectively. These unique designs of CIPWR make it a compact PWR and reduce the risk of some accidents such as LOCA effectively.Figure 3
Sketch of CIPWR.Figure4 shows the node map of CIPWR based on RELAP5/MODE3.4. Reactor kinetics chose point reactor model, and the reactor power is 100 MW. Average channel model is adopted in modelling reactor core, which can meet the requirements of steady-state calculation. Although the real arrangement in core should be modelled to compute more accurately, lots of fuel elements nodes tremendously increase computation time, which is not necessary in this study. As is shown in Figure 5, the core channel is modelled by a pipe and a heat structure with three materials, which are uranium dioxide for core fuel, helium for the gap, and zirconium alloy for fuel cladding from the inside out.Figure 4
Node map of CIPWR.Figure 5
Sketch of core materials.OTSG has high power density and can be adjusted rapidly, which makes it conventionally applied in IPWR. Figure6 shows the node map of bushing OTSG CIPWR used. The feedwater in second side obtains heat by convection with first side of OTSG and turns to the steam along the flow channel. The design of two pipes in the primary loop makes heat transfer effective and structure compact. The material of heat structure between two sides of OTSG is Inconel 600 MA, which is widely used for steam generator tubes.Figure 6
Node map of OTSG.


## 2.1. Heat Loss Calculation Package

The package added to RELAP5 is used to calculate the heat loss in heat balance analysis. All the input values in the package can be obtained from the calculation result of steady condition by RELAP5, such as the value of temperature, pressure, and vapour fraction. Figure1 shows the flow chart of the package.Figure 1
Structure of the heat loss package.


### 2.1.1. The Principle of Heat Balance

Heat balance method is based on the steam generator approximation enthalpy balance [11]; that is, the enthalpy increase of feedwater is the main part of the thermal power. The thermal power can be calculated through the enthalpy difference between the steam and feedwater in the second loop, counting the internal heat loss at the same time.As is shown in Figure2, there are several ways for CIPWR primary system to obtain or lose heat (for other PWRs, the pressurizer electrical heating should be taken into consideration). The heat increase consists of core heat production, the heat brought in by charging stream, and heat production during the pump operation. The heat decrease consists of heat dissipation of components and the heat taken by letdown stream. The feedwater in the second loop exchanges heat with primary loop by U tubes. The heat balance principle of the CIPWR can be described by
(1)WR=WSG-WΔPr⁡,
where the WR (MW) is the core thermal power, which can be obtained from the RELAP5 calculation result and WSG (MW) is the steam generator power obtained from the primary loop, which can be obtained by the package. WΔPr⁡ (MW) is the heat loss, which is the sum of thermal power coolant obtained from other heat sources except core and lost because of the heat dissipation of components and pipes. The heat loss WΔPr⁡ (MW) is the calculation result of the package.Figure 2
Heat balance principle diagram of CIPWR.In the second loop, unsaturated feedwater with the mass flow rateQe (kg/s) and specific enthalpy He flows through the heat transfer tube and convects with the first side of steam generator and turns to the steam with the specific enthalpy Hv (kJ/kg). Equation (2) can be used to calculate the steam generator power:
(2)WSG=(Hv-He)Qe,
where Hv is the steam enthalpy of steam in the generator export, He (kJ/kg) is the feedwater enthalpy, and Qe (kg/s) is the flow rate of feedwater.


### 2.1.2. Thermodynamic Properties Subroutine

The enthalpy of water and steam included in the heat balance calculation can be computed by the thermodynamic properties subroutine according to IAPWS-IF97 [12] (IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam). IAPWS-IF97 replaced the previous industrial formulation IFC-67 [13] for industrial use for its advantage of calculation with speed and accuracy.The range of IAPWS-IF97 has five subareas, each of which has a basic equation that can be used to calculate the water and steam properties by derivation. As one of the thermophysical properties, enthalpy can be obtained directly from the pressure and temperature values by the basic equations except region 3. The feedwater in the entrance of steam generator is unsaturated water belonging to region 1, whose thermophysical properties can be computed by (3). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(3)g(p,T)RT=γ(π,τ)=∑i=134ni(7.1-π)Ii(τ-1.222)Ji,
whereπ=p/p* and τ=T*/T with p*=16.53 MPa, T*=1386 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients ni and exponents Ii and Ji are given in the IAPWS-IF97.The enthalpy of feedwater can be computed by(4)He=g-T(∂g∂T)p=τRT×∑i=134ni(7.1-π)IiJi(τ-1.222)Ji-1.The steam in the export of steam generator loop is wet saturated steam consisting of saturated steam and saturated water with different enthalpy. The total enthalpy of steam can be calculated according to the vapor fractionx as
(5)Hv=xHvs+(1-x)Hls,
where Hls (kJ/kg) is the saturated water enthalpy, which can be calculated by (4). Hvs is the saturated steam enthalpy belonging to region 2, which can be calculated by (6). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(6)g(p,T)RT=γ(π,τ)=ln⁡π+∑i=19nioτJio+∑i=143niπIi(τ-0.5)Ji.The enthalpy of saturated steamHvs (kJ/kg) can be computed by
(7)Hvs=g-T(∂g∂T)p=τRT×{1π+∑i=143niIiπIi-1(τ-0.5)Ji},
where π=p/p* and τ=T*/T with p*=1 MPa, T*=540 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients nio, Jio, ni, Ii, and Ji are given in the IAPWS-IF97.As shown in Table1, the validation of this subroutine has been done by comparing it with the value given in IAPWS-IF97.Table 1
Comparison of calculation results between package subroutine and IAPWS-IF97.
Region
Parameters
Given by
h (kJ kg−1)
Region 1
T = 300 K
IF97
115.331273
P = 3 MPa
Package
115.33127302
T = 300 K
IF97
184.142828
P = 80 MPa
Package
184.14282773
Region 2
T = 300 K
IF97
2549.91145
P = 0.0035 MPa
Package
2549.9114508
T = 2000 K
IF97
6571.22604
P = 30 MPa
Package
6571.2260386


## 2.1.1. The Principle of Heat Balance

Heat balance method is based on the steam generator approximation enthalpy balance [11]; that is, the enthalpy increase of feedwater is the main part of the thermal power. The thermal power can be calculated through the enthalpy difference between the steam and feedwater in the second loop, counting the internal heat loss at the same time.As is shown in Figure2, there are several ways for CIPWR primary system to obtain or lose heat (for other PWRs, the pressurizer electrical heating should be taken into consideration). The heat increase consists of core heat production, the heat brought in by charging stream, and heat production during the pump operation. The heat decrease consists of heat dissipation of components and the heat taken by letdown stream. The feedwater in the second loop exchanges heat with primary loop by U tubes. The heat balance principle of the CIPWR can be described by
(1)WR=WSG-WΔPr⁡,
where the WR (MW) is the core thermal power, which can be obtained from the RELAP5 calculation result and WSG (MW) is the steam generator power obtained from the primary loop, which can be obtained by the package. WΔPr⁡ (MW) is the heat loss, which is the sum of thermal power coolant obtained from other heat sources except core and lost because of the heat dissipation of components and pipes. The heat loss WΔPr⁡ (MW) is the calculation result of the package.Figure 2
Heat balance principle diagram of CIPWR.In the second loop, unsaturated feedwater with the mass flow rateQe (kg/s) and specific enthalpy He flows through the heat transfer tube and convects with the first side of steam generator and turns to the steam with the specific enthalpy Hv (kJ/kg). Equation (2) can be used to calculate the steam generator power:
(2)WSG=(Hv-He)Qe,
where Hv is the steam enthalpy of steam in the generator export, He (kJ/kg) is the feedwater enthalpy, and Qe (kg/s) is the flow rate of feedwater.


## 2.1.2. Thermodynamic Properties Subroutine

The enthalpy of water and steam included in the heat balance calculation can be computed by the thermodynamic properties subroutine according to IAPWS-IF97 [12] (IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam). IAPWS-IF97 replaced the previous industrial formulation IFC-67 [13] for industrial use for its advantage of calculation with speed and accuracy.The range of IAPWS-IF97 has five subareas, each of which has a basic equation that can be used to calculate the water and steam properties by derivation. As one of the thermophysical properties, enthalpy can be obtained directly from the pressure and temperature values by the basic equations except region 3. The feedwater in the entrance of steam generator is unsaturated water belonging to region 1, whose thermophysical properties can be computed by (3). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(3)g(p,T)RT=γ(π,τ)=∑i=134ni(7.1-π)Ii(τ-1.222)Ji,
whereπ=p/p* and τ=T*/T with p*=16.53 MPa, T*=1386 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients ni and exponents Ii and Ji are given in the IAPWS-IF97.The enthalpy of feedwater can be computed by(4)He=g-T(∂g∂T)p=τRT×∑i=134ni(7.1-π)IiJi(τ-1.222)Ji-1.The steam in the export of steam generator loop is wet saturated steam consisting of saturated steam and saturated water with different enthalpy. The total enthalpy of steam can be calculated according to the vapor fractionx as
(5)Hv=xHvs+(1-x)Hls,
where Hls (kJ/kg) is the saturated water enthalpy, which can be calculated by (4). Hvs is the saturated steam enthalpy belonging to region 2, which can be calculated by (6). This equation is a fundamental equation for the specific Gibbs free energy g expressed in dimensionless form:
(6)g(p,T)RT=γ(π,τ)=ln⁡π+∑i=19nioτJio+∑i=143niπIi(τ-0.5)Ji.The enthalpy of saturated steamHvs (kJ/kg) can be computed by
(7)Hvs=g-T(∂g∂T)p=τRT×{1π+∑i=143niIiπIi-1(τ-0.5)Ji},
where π=p/p* and τ=T*/T with p*=1 MPa, T*=540 K, and R=0.461526 kJ kg−1 K−1 and p(MPa) and T (K) are pressure and temperature value of the unsaturated water. The coefficients nio, Jio, ni, Ii, and Ji are given in the IAPWS-IF97.As shown in Table1, the validation of this subroutine has been done by comparing it with the value given in IAPWS-IF97.Table 1
Comparison of calculation results between package subroutine and IAPWS-IF97.
Region
Parameters
Given by
h (kJ kg−1)
Region 1
T = 300 K
IF97
115.331273
P = 3 MPa
Package
115.33127302
T = 300 K
IF97
184.142828
P = 80 MPa
Package
184.14282773
Region 2
T = 300 K
IF97
2549.91145
P = 0.0035 MPa
Package
2549.9114508
T = 2000 K
IF97
6571.22604
P = 30 MPa
Package
6571.2260386


## 2.2. Modelling of Steady Operation of CIPWR

The whole reactor system can be modelled by different modules of RELAP5, which contains hydrodynamic components, branches, heat structures, reactor kinetics, control system, and so on. Hydrodynamic components consist of the parts where the fluid flows through, such as pipes, annulus, valves, pumps, and time-dependent volumes. Heat structures simulate the heat transfer in the system by setting the thermal conductivity and heat capacity of materials. Control system can set the change with time of physical quantities and trips and also can be used in self-initialization. Reactor kinetics consist of multigroup neutron transport model and point reactor model, which calculate the power of reactor by setting the initial power and reactivity feedback coefficient of various influent factors.Figure3 shows the sketch of CIPWR with two loops with coolant pumps. Once-through steam generators (OTSG) are adopted for the advantage that the configuration is compact enough to enable their location in the pressure vessel. OTSGs are arranged above the reactor core for the better natural circulation in case of SBO accident. The upper pressure vessel space acts as pressurizer in CIPWR, which eliminates lots of pipes effectively. These unique designs of CIPWR make it a compact PWR and reduce the risk of some accidents such as LOCA effectively.Figure 3
Sketch of CIPWR.Figure4 shows the node map of CIPWR based on RELAP5/MODE3.4. Reactor kinetics chose point reactor model, and the reactor power is 100 MW. Average channel model is adopted in modelling reactor core, which can meet the requirements of steady-state calculation. Although the real arrangement in core should be modelled to compute more accurately, lots of fuel elements nodes tremendously increase computation time, which is not necessary in this study. As is shown in Figure 5, the core channel is modelled by a pipe and a heat structure with three materials, which are uranium dioxide for core fuel, helium for the gap, and zirconium alloy for fuel cladding from the inside out.Figure 4
Node map of CIPWR.Figure 5
Sketch of core materials.OTSG has high power density and can be adjusted rapidly, which makes it conventionally applied in IPWR. Figure6 shows the node map of bushing OTSG CIPWR used. The feedwater in second side obtains heat by convection with first side of OTSG and turns to the steam along the flow channel. The design of two pipes in the primary loop makes heat transfer effective and structure compact. The material of heat structure between two sides of OTSG is Inconel 600 MA, which is widely used for steam generator tubes.Figure 6
Node map of OTSG.


## 3. Calculation Result

### 3.1. Steady Operation of CIPWR Modelled by RELAP5

Steady state of CIPWR has been modelled by RELAP5. Figure8 shows the pressure values of reactor core, feedwater in the entrance of OTSG, steam in the export of OTSG, and the equivalent pressurizer in the upper pressure vessel space. As shown in Figure 7, the pressure values are nearly stable during the operation; that is, the pressure value of core remains 7.1127 MPa, the pressure value of upper plenum remains 7.0989 MPa, the pressure value of low plenum remains 7.1216 MPa, and the pressure value of the equivalent pressurizer remains 7.07127 MPa, which shows that the upper pressure vessel space can act as the pressurizer. The distribution of pressure is in agreement with the design pressure value (7.063 MPa) of pressurizer.Figure 7
Pressure of reactor core, upper plenum, low plenum, steam, and pressurizer in steady state.Figure 8
Temperature of reactor core, upper plenum, low plenum, steam, and second side of SG U tube in steady state.Figure8 shows the temperature values of reactor core, feedwater in the entrance of OTSG, steam in the export of OTSG, and the equivalent pressurizer in the upper pressure vessel space. As shown in the figure, the temperature values reach steady state after about 10 seconds and are consistent with the trend of heat transfer.Table2 shows calculation values of pressure and temperature in different positions modelled by RELAP5.Table 2
RELAP5 calculation values of pressure and temperature in different positions.
Parameter
Position
Values
Pressure (MPa)
Core
7.1127
Upper plenum
7.0989
Low plenum
7.1216
Pressurizer
7.07127
SG feedwater
3.0 MPa
SG steam
2.943 MPa
Temperature (K)
Core
0.9667
Upper plenum
557.59
Low plenum
515.17
SG primary side
547.92
SG feedwater
423.15
SG steam
505.93


### 3.2. The Calculation of Heat Loss

When steady-state operation finished, the heat loss calculation package started, using the value in the RELAP calculation result. The calculation result is shown in Table3; the heat loss of CIPWR is 1.3454 MW.Table 3
Values of heat loss calculation.
Belong to
Parameters
Values
RELAP5
Feedwater temperature
423.15 K
Feedwater pressure
3.0 MPa
Steam temperature
506.05 K
Steam pressure
2.943 MPa
Vapor fraction
0.9667
Reactor power
96.58 MW
Feedwater flow rate
43.9 kg/s
Heat loss calculation package
Feedwater enthalpy
633.813 kJ/kg
Steam enthalpy
2802.974 kJ/kg
OTSG power
95.226 MW
Heat loss
1.354 MWAs is shown in Table4, all of the errors of the main values of parameters needed in the heat loss calculation between RELAP5 and the package are less than 0.3%, which means that the reasonable calculation result can be obtained by the heat loss package.Table 4
Difference and error of calculation results between RELAP5 and package.
Parameter
RELAP5 calculation results
Package calculation results
Difference
Error (%)
Feedwater enthalpy (kJ/kg)
633.602
633.813
0.211
0.0333
Steam enthalpy (kJ/kg)
2798.343
2802.974
4.631
0.1655
OTSG Power (MW)
95.032
95.226
0.194
0.2041
Heat Loss (MW)
1.548
1.354
0.194
—


## 3.1. Steady Operation of CIPWR Modelled by RELAP5

Steady state of CIPWR has been modelled by RELAP5. Figure8 shows the pressure values of reactor core, feedwater in the entrance of OTSG, steam in the export of OTSG, and the equivalent pressurizer in the upper pressure vessel space. As shown in Figure 7, the pressure values are nearly stable during the operation; that is, the pressure value of core remains 7.1127 MPa, the pressure value of upper plenum remains 7.0989 MPa, the pressure value of low plenum remains 7.1216 MPa, and the pressure value of the equivalent pressurizer remains 7.07127 MPa, which shows that the upper pressure vessel space can act as the pressurizer. The distribution of pressure is in agreement with the design pressure value (7.063 MPa) of pressurizer.Figure 7
Pressure of reactor core, upper plenum, low plenum, steam, and pressurizer in steady state.Figure 8
Temperature of reactor core, upper plenum, low plenum, steam, and second side of SG U tube in steady state.Figure8 shows the temperature values of reactor core, feedwater in the entrance of OTSG, steam in the export of OTSG, and the equivalent pressurizer in the upper pressure vessel space. As shown in the figure, the temperature values reach steady state after about 10 seconds and are consistent with the trend of heat transfer.Table2 shows calculation values of pressure and temperature in different positions modelled by RELAP5.Table 2
RELAP5 calculation values of pressure and temperature in different positions.
Parameter
Position
Values
Pressure (MPa)
Core
7.1127
Upper plenum
7.0989
Low plenum
7.1216
Pressurizer
7.07127
SG feedwater
3.0 MPa
SG steam
2.943 MPa
Temperature (K)
Core
0.9667
Upper plenum
557.59
Low plenum
515.17
SG primary side
547.92
SG feedwater
423.15
SG steam
505.93


## 3.2. The Calculation of Heat Loss

When steady-state operation finished, the heat loss calculation package started, using the value in the RELAP calculation result. The calculation result is shown in Table3; the heat loss of CIPWR is 1.3454 MW.Table 3
Values of heat loss calculation.
Belong to
Parameters
Values
RELAP5
Feedwater temperature
423.15 K
Feedwater pressure
3.0 MPa
Steam temperature
506.05 K
Steam pressure
2.943 MPa
Vapor fraction
0.9667
Reactor power
96.58 MW
Feedwater flow rate
43.9 kg/s
Heat loss calculation package
Feedwater enthalpy
633.813 kJ/kg
Steam enthalpy
2802.974 kJ/kg
OTSG power
95.226 MW
Heat loss
1.354 MWAs is shown in Table4, all of the errors of the main values of parameters needed in the heat loss calculation between RELAP5 and the package are less than 0.3%, which means that the reasonable calculation result can be obtained by the heat loss package.Table 4
Difference and error of calculation results between RELAP5 and package.
Parameter
RELAP5 calculation results
Package calculation results
Difference
Error (%)
Feedwater enthalpy (kJ/kg)
633.602
633.813
0.211
0.0333
Steam enthalpy (kJ/kg)
2798.343
2802.974
4.631
0.1655
OTSG Power (MW)
95.032
95.226
0.194
0.2041
Heat Loss (MW)
1.548
1.354
0.194
—


## 4. Conclusions

The comparison between the simulated value and design value shows that the steady state of CIPWR has been modelled by RELAP5 correctly. According to the distribution of simulated values of pressure in primary loop, the upper pressure vessel space can act as pressurizer practically with the advantage of eliminating lots of pipes which reduces the risk of some accidents effectively.The heat loss calculation package which is added to RELAP5 was developed based on heat balance method. The comparison of the main values needed in the heat loss calculation between RELAP5 and the package shows that the heat loss package can be applied in reactor design and monitoring. Further study will aim at calculating the series of transient accidents of CIPWR based on the steady state by RELAP5.



---

*Source: 707384-2014-09-01.xml*