# Research and Evaluation for Passive Safety System in Low Pressure Reactor

Authors: Peng Chuanxin; Zhuo Wenbin; Chen Bingde; Nie Changhua; Huang Yanping
Journal: Science and Technology of Nuclear Installations
 (2015)

---

## Abstract

Low pressure reactor is a small size advanced reactor with power of 180 MWt, which is under development at Nuclear Power Institute of China. In order to assess the ability and feasibility of passive safety system, several tests have been implemented on the passive safety system (PSS) test facility. During the LOCA and SBO accident, the adequate core cooling is provided by the performance of passive safety system. In addition the best-estimate thermal hydraulic code, CATHARE V2.1, has been assessed against cold leg LOCA test. The calculation results show that CATHARE is in a satisfactory agreement with the test for the steady state and transient test.


---

## Body

## 1. Introduction

The use of passive safety systems eliminates the costs associated with the installation, maintenance, and operation of active safety systems that require multiple pumps with independent and redundant electric power supplies [1]. As a result, passive safety systems are being considered for numerous reactor concepts, such as AP600/AP1000 [2], APR1400 [3], APWR+ [4], SMART [5], IRIS [6], and MASLWR [7].The AP600 was the first passively safe nuclear plant to be certified in the United States. It employs passive safety systems that rely on gravity (core make-up tank, in-containment refueling water storage tank), compressed gas (accumulator), natural circulation (passive residual heat removal heat exchanger), and evaporation to provide for long-term cooling in the event of an accident.The certification of AP600/AP1000 was based on comprehensive integral system and separate effects testing conducted by Westinghouse and the US Department of Energy at the SPES test facility in Italy and at the APEX test facility at Oregon State University. The purpose of SPES-2 tests was to simulate AP600 response to postulated accident events and to validate safety analysis codes [8]. The second integral facility APEX emphasized the depressurization, IRWST injection, and long-term cooling during small break LOCA [9]. Besides, US NRC conducted independent integral test programs at ROSA in Tokai-mura, Japan, focusing on the high pressure and depressurization phase, the initiation of IRWST injection [10], and the passive residual heat removal [11].Low pressure reactor is a small size advanced reactor with power of 180 MWt, which is under development at Nuclear Power institute of China. The primary system of reactor consists of a reactor pressure vessel, four loops, and the tank. Each loop is made up of heat exchanger and reactor coolant pump, hot leg and cold leg, and so forth. The reactor produces 200,000 m3/day of hot water (90°C) for low temperature seawater desalination. The low pressure reactor operates at low temperature (temperature at the outlet is 98°C) and atmospheric pressure (pressure at the top of tank is 0.1 MPa). It follows that the consequences of accidents are less severe than typical PWR. In addition, it adopts passive engineered safety systems such as emergency core cooling system and passive residual heat removal system to mitigate consequences of accidents. Figure 1 shows the schematic diagram of the low pressure reactor. The tank is a significant component of low pressure reactor, which is a very large pool filled with cold borated water with a volume of 1310 m3. During the normal operation, the tank is used to maintain the primary system pressure. And during the shutdown or non-LOCA, tank is treated as heat sink for the core decay heat. Another important role of the tank is to provide water source for core injection if LOCA occurs. As the operational pressure is about 0.4 MPa, the head pressure of the water in the tank is adequate for core injection. Therefore the design of Accumulator and Automatic Depressurization System is not necessary for the low pressure reactor.Figure 1
Schematic diagram of low pressure reactor.In order to assess the ability and feasibility of low pressure reactor passive safety system, several tests have been implemented on the passive safety system (PSS) test facility [12]. In this paper the experimental investigations on the characteristics of the passive safety systems under LOCA and SBO are reported in detail. In addition the experimental data are used to assess CATHARE simulating the cold leg LOCA test.


## 2. Passive Safety System Test Facility

The passive safety system test facility is employed to simulate various accidental scenarios (such as loss-of-coolant accident and station black-out accident) that might occur in the low pressure reactor. To preserve the similar thermal hydraulic transients, the PSS facility has a geometrical scaling ratio of 1/1 in length and height and 1/45 in volume, and it operates at full pressure and temperature conditions of prototype [13]. Table 1 shows the scaling ratios of the facility to the low pressure reactor for some major geometrical and thermal hydraulic parameters. The schematic diagram of PSS facility is shown in Figure 2. And its main features are summarized as follows.(1)
The reactor vessel is simulated by the core simulation vessel and downcomer annulus simulation vessel. A total number of 96 electrically heated rods with a maximum total power of 600 kW are installed in the core simulation vessel. The break simulating system consists of a quick open valve, a venturi flow meter, and the core external cavity simulation vessel. The break flow rate was directly measured with a venturi flow meter. The venturi flow meter was installed at the upstream of the break nozzle for the break flow to be single-phase water flow at the measurement location. The accumulated water in the core external cavity simulation vessel shows good agreement with the break flow rate measured by venturi flow meter.
(2)
The PSS test facility is designed to be operated by a combination of manual and automatic operations. The controlled components include the electrically heated rod, the primary coolant pump, the primary loop flow control valve, and the heat exchanger cooling water flow control valve. All the safety related accidents should be initiated and controlled by automatic control logics as they require automatic reactor trip logics to initiate the passive core cooling systems and the passive residual heat removal system. The corresponding trip logics, tables, and set points are programmed to control a sequence of events.
(3)
The data acquisition system involves 55 instruments to measure system pressure, differential pressure, fluid temperature, water level, and flow rate. The precision of nickel-chrome and nickel-silicon thermocouples is within ±1.5°C. The pressure of the system is measured by the pressure transducer, whose precision is within ±0.5%. The flow rate is measured by wide range venturi flow meters, the precision of which is within ±1%.Table 1
Major scaling ratios of PSS.
Parameters
Scaling ratio
PSS
Parameters
Scaling ratio
PSS
Length (height)
l
oR
1 : 1
Diameter
d
oR
1 : 6.71
Area
d
oR
2
1 : 45
Volume
l
oR
d
oR
2
1 : 45
Core temperature rise
T
oR
1 : 1
Velocity
l
oR
1
/
2
1 : 1
Time
l
oR
1
/
2
1 : 1
Power/volume
l
oR
1
/
2
1 : 1
Heat flux
l
oR
1
/
2
1 : 1
Core power
l
oR
1
/
2
d
oR
2
1 : 1
Diameter of rod
1
1 : 1
No. of rods
d
oR
2
1 : 45
Flow rate
l
oR
1
/
2
d
oR
2
1 : 45
Pressure drop
l
oR
1 : 1Figure 2
Schematic diagram of PSS test facility.


## 3. Description of the Tests

A series of tests devoted to validate the efficiency of low pressure reactor passive safety system were performed in the recent years. The experimental program includes two series with altogether eighteen experiments. The experimental series and conditions are listed in Table2. The first series focus on validation and research of passive core cooling behaviour, such as cold leg LOCA, hot leg LOCA, and the influence of sensitive parameters. The second series study SBO accident and the influence of sensitive parameters on the passive residual heat removal behaviour [5]. In this paper the cold leg LOCA and SBO test results are presented.Table 2
Typical test matrix.
Experiment
Condition(Break location, break size, and water level and temperature in tank)
Objectives
First series
PSS-11
Cold leg, 22 mm, 23 m, 32°C
Validation of passive core injection
PSS-12
Hot leg, 22 mm, 23 m, 32°C
Validation of passive core injection
PSS-13
Cold leg, 18 mm, 23 m, 32°C
Break size effect
PSS-14
Cold leg, 14 mm, 23 m, 32°C
Break size effect
PSS-15
Cold leg, 10 mm, 23 m, 32°C
Break size effect
PSS-16
Hot leg, 18 mm, 23 m, 32°C
Break size effect
PSS-17
Hot leg, 14 mm, 23 m, 32°C
Break size effect
PSS-18
Hot leg, 10 mm, 23 m, 32°C
Break size effect
PSS-19
Cold leg, 22 mm, 18 m, 32°C
Water level effect
PSS-110
Cold leg, 22 mm, 13 m, 32°C
Water level effect
PSS-111
Cold leg, 22 mm, 23 m, 22°C
Water temperature effect
PSS-112
Cold leg, 22 mm, 23 m, 42°C
Water temperature effect
Experiment
Condition(Resistance*, startup mode*, and water temperature in tank)
Objectives
Second series
PSS-21
100% of prototype, cold startup, 32°C
Validation of passive residual heat removal
PSS-22
120% of prototype, cold startup, 32°C
Resistance effect
PSS-23
150% of prototype, cold startup, 32°C
Resistance effect
PSS-24
100% of prototype, hot startup, 32°C
Startup mode effect
PSS-25
100% of prototype, cold startup, 22°C
Water temperature effect
PSS-26
100% of prototype, cold startup, 42°C
Water temperature effect
*Resistance includes.
Startup mode includes cold startup and hot startup.
Cold startup: the fluid in the PL is cold before the natural circulation is established.
Hot startup: the fluid in the PL is warm before the natural circulation is established.


### 3.1. Cold Leg LOCA Test

The passive system is designed to cool the reactor in more than 48 h, and the fluid temperature in the tank should be less than 80°C, when the LOCA accident occurs. The motivation of cold leg LOCA test (PSS-11) is to evaluate the safety features of low pressure reactor.The postulated cold leg LOCA is the one initiated by opening the break simulation valve, causing the system pressure to decrease. Thereafter, the check valve of injection line is actuated automatically. Then the coolant in the tank was injected into reactor through the surge line. The discharge of coolant in the tank causes liquid level to decrease. A scram signal trips, when the pressure at the core outlet is less than 350 kPa. Following the S-signal, the reactor is tripped; the reactor coolant pumps coasted down; the feedwater supply of heat exchanger is terminated; the passive residual heat removal line and injection line isolation valves are open. According to the characteristics of the test, the overall transients of the cold leg LOCA test can be roughly characterized by typical two stages: injection stage and natural circulation stage.The typical results of cold leg LOCA test are shown in Figures3~6. Figure 3 shows the core pressure. As the break occurs at the cold leg, the system pressure decreases. The top of tank is open to atmosphere, and the initial liquid level is 23 m (the bottom of tank is 0.15 m higher than the top of reactor vessel). Therefore the water in the tank can be injected into reactor vessel by gravity.Figure 3
The system pressure.Figure4 shows the flow rate of some lines. As the liquid level in the tank descends, the mass flow rate of break, IL, SL, and PL decreases. When the level in the tank is equal to that in the core external cavity simulation vessel, the break flow stops and the natural circulation stage begins.Figure 4
The mass flow rate of the main lines.The temperature of fluid in the core is illustrated in Figure5. At the initial period, the fluid temperature drops quickly, and the core outlet temperature is lower than that in the inlet. Because most of injection of SL and IL bypasses the core simulation vessel, however the cold water form PL can be injected into the upper plenum. As a result, the reverse flow happens in the core simulation vessel. As the flow rate injected into the core decreases, the temperature begins to rise after 2.8 h. When the fluid in the core is stagnant, the temperature in the core inlet and outlet rises to 52.4°C and 69.4°C, respectively. The fluid temperature in the core does not rise until the stable natural circulation is established.Figure 5
Temperature at the core inlet and outlet.Figure 6
Fluid temperatures in tank.In the natural circulation stage, cold water in the tank is injected into the core through SL and IL. Hot water from the core outlet rises along PL and returns to the tank, which is driven by buoyancy. The core residual heat is removed into the tank by natural circulation. Figure6 shows the fluid temperatures in tank. When the natural circulation begins (about 4.1 h after accident happens), the fluid temperature at the top of tank rises. After 48 h, the fluid temperature in the tank is less than 60°C. Therefore, the adequate core cooling is provided by passive safety system dominantly, during the cold leg LOCA.


### 3.2. Station Black-Out Accident

The passive residual heat removal system is designed to remove the decay heat in more than 72 h, and the fluid temperature in the tank should be less than 80°C, when the station black-out accident occurs. The motivation of station black-out accident (PSS-21) is to validate the ability of the passive residual heat removal system.The postulated station black-out accident is the one initiated by stopping the reactor coolant pumps. Following the accident, the reactor is tripped; only the passive safety systems are operational to remove the core decay heat or to mitigate the consequence of the accident.The typical results of station black-out accident test are shown in Figures7~9. Figure 7 shows the flow rate in the two parallel natural circulation paths. When the station black-out accident is initiated, the natural circulation begins at the time of reactor coolant pumps coastdown. The two parallel natural circulations “core—residual heat removal line—tank—safety injection line—core” and “core—residual heat removal line—tank—surge line—core” remove the core decay heat to the tank.Figure 7
Natural circulation flow rate.Figure8 shows the temperature of fluid at the core inlet and outlet. At the beginning of the accident, the natural circulation is unable to remove the core decay heat effectively. And a small portion of coolant at the core outlet evaporates.  The two-phase mixture flows into the residual heat removal line, resulting in the increase of natural circulation. The measurements of the flow in the PL line, considered just for single-phase flow, are corrected by flow rate of IL and SL. When the natural circulation comes to being stable, the temperature in the core drops rapidly.Figure 8
Core inlet and outlet fluid temperatures.Figure 9
Fluid temperatures in tank.Figure9 shows the temperature of coolant in the tank. The residual heat removal line connects to the top of tank; therefore, the water in the top becomes hot firstly. As the core decay heat is taken to the tank constantly, the cold water in the tank is replaced by hot water gradually. After 72 h, the coolant in the core is still lower than 60°C. Therefore the experimental results show that the single-phase natural circulation flows steadily into the two parallel loops that it effectively removes the decay heat from the core to tank in more than 72 hours.


## 3.1. Cold Leg LOCA Test

The passive system is designed to cool the reactor in more than 48 h, and the fluid temperature in the tank should be less than 80°C, when the LOCA accident occurs. The motivation of cold leg LOCA test (PSS-11) is to evaluate the safety features of low pressure reactor.The postulated cold leg LOCA is the one initiated by opening the break simulation valve, causing the system pressure to decrease. Thereafter, the check valve of injection line is actuated automatically. Then the coolant in the tank was injected into reactor through the surge line. The discharge of coolant in the tank causes liquid level to decrease. A scram signal trips, when the pressure at the core outlet is less than 350 kPa. Following the S-signal, the reactor is tripped; the reactor coolant pumps coasted down; the feedwater supply of heat exchanger is terminated; the passive residual heat removal line and injection line isolation valves are open. According to the characteristics of the test, the overall transients of the cold leg LOCA test can be roughly characterized by typical two stages: injection stage and natural circulation stage.The typical results of cold leg LOCA test are shown in Figures3~6. Figure 3 shows the core pressure. As the break occurs at the cold leg, the system pressure decreases. The top of tank is open to atmosphere, and the initial liquid level is 23 m (the bottom of tank is 0.15 m higher than the top of reactor vessel). Therefore the water in the tank can be injected into reactor vessel by gravity.Figure 3
The system pressure.Figure4 shows the flow rate of some lines. As the liquid level in the tank descends, the mass flow rate of break, IL, SL, and PL decreases. When the level in the tank is equal to that in the core external cavity simulation vessel, the break flow stops and the natural circulation stage begins.Figure 4
The mass flow rate of the main lines.The temperature of fluid in the core is illustrated in Figure5. At the initial period, the fluid temperature drops quickly, and the core outlet temperature is lower than that in the inlet. Because most of injection of SL and IL bypasses the core simulation vessel, however the cold water form PL can be injected into the upper plenum. As a result, the reverse flow happens in the core simulation vessel. As the flow rate injected into the core decreases, the temperature begins to rise after 2.8 h. When the fluid in the core is stagnant, the temperature in the core inlet and outlet rises to 52.4°C and 69.4°C, respectively. The fluid temperature in the core does not rise until the stable natural circulation is established.Figure 5
Temperature at the core inlet and outlet.Figure 6
Fluid temperatures in tank.In the natural circulation stage, cold water in the tank is injected into the core through SL and IL. Hot water from the core outlet rises along PL and returns to the tank, which is driven by buoyancy. The core residual heat is removed into the tank by natural circulation. Figure6 shows the fluid temperatures in tank. When the natural circulation begins (about 4.1 h after accident happens), the fluid temperature at the top of tank rises. After 48 h, the fluid temperature in the tank is less than 60°C. Therefore, the adequate core cooling is provided by passive safety system dominantly, during the cold leg LOCA.


## 3.2. Station Black-Out Accident

The passive residual heat removal system is designed to remove the decay heat in more than 72 h, and the fluid temperature in the tank should be less than 80°C, when the station black-out accident occurs. The motivation of station black-out accident (PSS-21) is to validate the ability of the passive residual heat removal system.The postulated station black-out accident is the one initiated by stopping the reactor coolant pumps. Following the accident, the reactor is tripped; only the passive safety systems are operational to remove the core decay heat or to mitigate the consequence of the accident.The typical results of station black-out accident test are shown in Figures7~9. Figure 7 shows the flow rate in the two parallel natural circulation paths. When the station black-out accident is initiated, the natural circulation begins at the time of reactor coolant pumps coastdown. The two parallel natural circulations “core—residual heat removal line—tank—safety injection line—core” and “core—residual heat removal line—tank—surge line—core” remove the core decay heat to the tank.Figure 7
Natural circulation flow rate.Figure8 shows the temperature of fluid at the core inlet and outlet. At the beginning of the accident, the natural circulation is unable to remove the core decay heat effectively. And a small portion of coolant at the core outlet evaporates.  The two-phase mixture flows into the residual heat removal line, resulting in the increase of natural circulation. The measurements of the flow in the PL line, considered just for single-phase flow, are corrected by flow rate of IL and SL. When the natural circulation comes to being stable, the temperature in the core drops rapidly.Figure 8
Core inlet and outlet fluid temperatures.Figure 9
Fluid temperatures in tank.Figure9 shows the temperature of coolant in the tank. The residual heat removal line connects to the top of tank; therefore, the water in the top becomes hot firstly. As the core decay heat is taken to the tank constantly, the cold water in the tank is replaced by hot water gradually. After 72 h, the coolant in the core is still lower than 60°C. Therefore the experimental results show that the single-phase natural circulation flows steadily into the two parallel loops that it effectively removes the decay heat from the core to tank in more than 72 hours.


## 4. Computer Code Simulations

### 4.1. Code Description

CATHARE is a system code developed by CEA, EDF, FRAMATOME-ANP, and IRSN for pressurized water reactor (PWR) safety analysis. CATHARE is designed to perform best-estimate analyses of LOCA, operational transients, and other accident scenarios in pressurized water reactor. It can also model phenomena occurring in experimental facility designed to simulate transients in reactor systems. CATHARE modules are based on a six-equation two-fluid model (mass, energy, and momentum equations for each phase), with additional optional equations for noncondensable gases and radiochemical components. A specific treatment of the residual phases exists in order to manage their appearance and disappearance while minimizing convergence problems and with a quasiperfect mass and energy conservation [14]. The version of the code used in this paper is CATHARE V2.1.


### 4.2. System Model

CATHARE has a flexible modular structure for the thermal-hydraulic modeling in applications ranging from simple experimental test facilities to large and complex installations like nuclear power plants. The main hydraulic components or elements are pipes (1D), volumes (0D), 3D vessel, and boundary conditions, connected to each other by junctions. These models except 3D vessel are based on the 6-equation model, which are two mass balance equations, two momentum balance equations, and two energy balance equations. In CATHARE, the critical flow rate is also determined by the some characteristic of the 6-equation model. The physical closure laws play a dominant role in the calculation of critical flow rate. The critical flow rate for nonhomogeneous nonequilibrium two-phase flow is obtained by solving the Jacobian matrix at each time step. The liquid and vapour velocity, obtained by solving the 6-equation system, are then compared with the choking criterion (known as the GROS D’AILLON correlation) [15]. This method is qualified using the results of experimental facilities, such as MOBY DICK facility, BETHSY BREAK NOZZLE characterization tests, and MARVIKEN facility [16]. These facilities are high pressure experiments. In this paper, the break flow rate is validated by experimental data from low pressure:(1)
G
CR
=
S
×
20
×
0.45
6
2
×
ρ
m
×
(
P
-
P
sat
⁡
(
T
L
)
)
α
+
10
×
0.45
6
2
×
1
-
α
×
1
+
0.012
×
L
/
D
+
k
Δ
P
s
,
ρ
m
=
α
×
ρ
g
+
1
-
α
ρ
L
,
P
-
P
sat
⁡
T
L
=
P
-
0.9
×
1
0
5
×
T
Leq
+
17.95
117.8
1
/
0.223
,
if
T
L
<
T
sat
⁡
P
,
P
10
,
if
T
L
>
T
sat
⁡
P
,
T
Leq
=
T
L
+
H
V
-
H
L
×
α
×
ρ
V
ρ
m
×
CP
L
,where G
CR is the critical flow rate choking criterion, S is the break area, ρ
m is the mean density, ρ
g is the gas density, ρ
L is the liquid density, T
Leq is the equilibrium temperature (after steam condensation), T
L is liquid temperature, T
g is the gas temperature, P is the pressure, H
L is the liquid enthalpy, H
g is the gas enthalpy, α is void fraction, and CP
L is the liquid specific heat.In this paper the core and bypass are modeled with 1D pipes, and the WALL is used to model the power produced by the electrically heated rods and the heat loss of test facility. The upper plenum, lower plenum, and tank are modeled by 0D volume. TEE + AXIAL+ BCONDIT is used to model the break. Figure10 illustrates a schematic nodalization of the major elements of the PSS test facility for CATHARE simulation.Figure 10
CATHARE nodalization diagram of the PSS facility.


### 4.3. Calculation Results

#### 4.3.1. Calculation of the Initial Steady State

At first, a steady state calculation with constant boundary conditions was performed for each test. With help of these calculations, the temperature and pressure distribution, the heat losses, the pressure losses, and the water level were adjusted. When the stable conditions are reached, the steady state calculation is stopped. The main parameters obtained at the end of steady state calculation are compared with the experimental data shown in Table3. The small difference between the calculation result and the experimental data remains below the sensors uncertainties. Therefore, we can consider that there is a good agreement between CATHARE code predictions for the initial steady state.Table 3
Measured and calculated initial conditions.
Parameter
Unit
Test
Calculation
Error %
Core power
kW
600
600.5
0.02
Pressure
kPa
398.0
396.3
−0.43
Loop 1 flow rate
kg/s
4.339
4.396
1.3
Loop 2 flow rate
kg/s
4.354
4.401
1.1
Core inlet fluid temperature
°C
91.2
91.0
−0.22
Core outlet fluid temperature
°C
98.0
98.6
0.61
Water level in tank
m
23.0
22.8
−0.87


#### 4.3.2. Calculation of Transient State

The transient calculation results of cold leg LOCA test are compared with experimental data in Figures11~13. The comparison of break flow rate is shown in Figure 11. We can see that the calculated break flow rate is in acceptable agreement with experimental data. For low pressure and highly subcooled inlet conditions, the break flow rate is obtained by solving the Bernoulli equation, and CATHARE can predict it accurately. There is a little discrepancy, which occurs at the end of injection stage. As for the instrumental error, the small break flow rate cannot be monitored by the wide range venturi flow meter.Figure 11
The break flow rate.Figure12 illustrates the calculated system pressure, which shows reasonable agreement with the experimental data, owing to the good predicted break mass flow rate. During the transition between injection stage and natural circulation stage, the main phenomenon is that the flow reverses in passive residual heat removal line. It is also predicted well by CATHARE as shown in Figure 13. In addition, the other calculated system parameters, such as the temperature of coolant at the core inlet and outlet, the IL mass flow rate, SL mass flow rate, and water level in tank, are in accordance with experimental data.Figure 12
The system pressure.Figure 13
The PL flow rate.


## 4.1. Code Description

CATHARE is a system code developed by CEA, EDF, FRAMATOME-ANP, and IRSN for pressurized water reactor (PWR) safety analysis. CATHARE is designed to perform best-estimate analyses of LOCA, operational transients, and other accident scenarios in pressurized water reactor. It can also model phenomena occurring in experimental facility designed to simulate transients in reactor systems. CATHARE modules are based on a six-equation two-fluid model (mass, energy, and momentum equations for each phase), with additional optional equations for noncondensable gases and radiochemical components. A specific treatment of the residual phases exists in order to manage their appearance and disappearance while minimizing convergence problems and with a quasiperfect mass and energy conservation [14]. The version of the code used in this paper is CATHARE V2.1.


## 4.2. System Model

CATHARE has a flexible modular structure for the thermal-hydraulic modeling in applications ranging from simple experimental test facilities to large and complex installations like nuclear power plants. The main hydraulic components or elements are pipes (1D), volumes (0D), 3D vessel, and boundary conditions, connected to each other by junctions. These models except 3D vessel are based on the 6-equation model, which are two mass balance equations, two momentum balance equations, and two energy balance equations. In CATHARE, the critical flow rate is also determined by the some characteristic of the 6-equation model. The physical closure laws play a dominant role in the calculation of critical flow rate. The critical flow rate for nonhomogeneous nonequilibrium two-phase flow is obtained by solving the Jacobian matrix at each time step. The liquid and vapour velocity, obtained by solving the 6-equation system, are then compared with the choking criterion (known as the GROS D’AILLON correlation) [15]. This method is qualified using the results of experimental facilities, such as MOBY DICK facility, BETHSY BREAK NOZZLE characterization tests, and MARVIKEN facility [16]. These facilities are high pressure experiments. In this paper, the break flow rate is validated by experimental data from low pressure:(1)
G
CR
=
S
×
20
×
0.45
6
2
×
ρ
m
×
(
P
-
P
sat
⁡
(
T
L
)
)
α
+
10
×
0.45
6
2
×
1
-
α
×
1
+
0.012
×
L
/
D
+
k
Δ
P
s
,
ρ
m
=
α
×
ρ
g
+
1
-
α
ρ
L
,
P
-
P
sat
⁡
T
L
=
P
-
0.9
×
1
0
5
×
T
Leq
+
17.95
117.8
1
/
0.223
,
if
T
L
<
T
sat
⁡
P
,
P
10
,
if
T
L
>
T
sat
⁡
P
,
T
Leq
=
T
L
+
H
V
-
H
L
×
α
×
ρ
V
ρ
m
×
CP
L
,where G
CR is the critical flow rate choking criterion, S is the break area, ρ
m is the mean density, ρ
g is the gas density, ρ
L is the liquid density, T
Leq is the equilibrium temperature (after steam condensation), T
L is liquid temperature, T
g is the gas temperature, P is the pressure, H
L is the liquid enthalpy, H
g is the gas enthalpy, α is void fraction, and CP
L is the liquid specific heat.In this paper the core and bypass are modeled with 1D pipes, and the WALL is used to model the power produced by the electrically heated rods and the heat loss of test facility. The upper plenum, lower plenum, and tank are modeled by 0D volume. TEE + AXIAL+ BCONDIT is used to model the break. Figure10 illustrates a schematic nodalization of the major elements of the PSS test facility for CATHARE simulation.Figure 10
CATHARE nodalization diagram of the PSS facility.


## 4.3. Calculation Results

### 4.3.1. Calculation of the Initial Steady State

At first, a steady state calculation with constant boundary conditions was performed for each test. With help of these calculations, the temperature and pressure distribution, the heat losses, the pressure losses, and the water level were adjusted. When the stable conditions are reached, the steady state calculation is stopped. The main parameters obtained at the end of steady state calculation are compared with the experimental data shown in Table3. The small difference between the calculation result and the experimental data remains below the sensors uncertainties. Therefore, we can consider that there is a good agreement between CATHARE code predictions for the initial steady state.Table 3
Measured and calculated initial conditions.
Parameter
Unit
Test
Calculation
Error %
Core power
kW
600
600.5
0.02
Pressure
kPa
398.0
396.3
−0.43
Loop 1 flow rate
kg/s
4.339
4.396
1.3
Loop 2 flow rate
kg/s
4.354
4.401
1.1
Core inlet fluid temperature
°C
91.2
91.0
−0.22
Core outlet fluid temperature
°C
98.0
98.6
0.61
Water level in tank
m
23.0
22.8
−0.87


### 4.3.2. Calculation of Transient State

The transient calculation results of cold leg LOCA test are compared with experimental data in Figures11~13. The comparison of break flow rate is shown in Figure 11. We can see that the calculated break flow rate is in acceptable agreement with experimental data. For low pressure and highly subcooled inlet conditions, the break flow rate is obtained by solving the Bernoulli equation, and CATHARE can predict it accurately. There is a little discrepancy, which occurs at the end of injection stage. As for the instrumental error, the small break flow rate cannot be monitored by the wide range venturi flow meter.Figure 11
The break flow rate.Figure12 illustrates the calculated system pressure, which shows reasonable agreement with the experimental data, owing to the good predicted break mass flow rate. During the transition between injection stage and natural circulation stage, the main phenomenon is that the flow reverses in passive residual heat removal line. It is also predicted well by CATHARE as shown in Figure 13. In addition, the other calculated system parameters, such as the temperature of coolant at the core inlet and outlet, the IL mass flow rate, SL mass flow rate, and water level in tank, are in accordance with experimental data.Figure 12
The system pressure.Figure 13
The PL flow rate.


## 4.3.1. Calculation of the Initial Steady State

At first, a steady state calculation with constant boundary conditions was performed for each test. With help of these calculations, the temperature and pressure distribution, the heat losses, the pressure losses, and the water level were adjusted. When the stable conditions are reached, the steady state calculation is stopped. The main parameters obtained at the end of steady state calculation are compared with the experimental data shown in Table3. The small difference between the calculation result and the experimental data remains below the sensors uncertainties. Therefore, we can consider that there is a good agreement between CATHARE code predictions for the initial steady state.Table 3
Measured and calculated initial conditions.
Parameter
Unit
Test
Calculation
Error %
Core power
kW
600
600.5
0.02
Pressure
kPa
398.0
396.3
−0.43
Loop 1 flow rate
kg/s
4.339
4.396
1.3
Loop 2 flow rate
kg/s
4.354
4.401
1.1
Core inlet fluid temperature
°C
91.2
91.0
−0.22
Core outlet fluid temperature
°C
98.0
98.6
0.61
Water level in tank
m
23.0
22.8
−0.87


## 4.3.2. Calculation of Transient State

The transient calculation results of cold leg LOCA test are compared with experimental data in Figures11~13. The comparison of break flow rate is shown in Figure 11. We can see that the calculated break flow rate is in acceptable agreement with experimental data. For low pressure and highly subcooled inlet conditions, the break flow rate is obtained by solving the Bernoulli equation, and CATHARE can predict it accurately. There is a little discrepancy, which occurs at the end of injection stage. As for the instrumental error, the small break flow rate cannot be monitored by the wide range venturi flow meter.Figure 11
The break flow rate.Figure12 illustrates the calculated system pressure, which shows reasonable agreement with the experimental data, owing to the good predicted break mass flow rate. During the transition between injection stage and natural circulation stage, the main phenomenon is that the flow reverses in passive residual heat removal line. It is also predicted well by CATHARE as shown in Figure 13. In addition, the other calculated system parameters, such as the temperature of coolant at the core inlet and outlet, the IL mass flow rate, SL mass flow rate, and water level in tank, are in accordance with experimental data.Figure 12
The system pressure.Figure 13
The PL flow rate.


## 5. Conclusions

The designed passive engineered safety systems are an important part of the low pressure reactor. The passive safety characteristics of the reactor following several safety related events have been experimentally investigated in the PSS test facility. In case of LOCA and SBO accident, the transient trend of the reactor primary loop, emergency core cooling system, and passive residual heat removal system are studied. In the cold leg LOCA test, the water level is enough for passive injection driven by gravity. After the injection stage stops, the natural circulation is established to remove the decay heat to tank continuously. During the station black-out accident, the natural circulation effectively removes the decay heat from the core to tank in more than 72 hours. In addition, the best-estimate thermal hydraulic code, CATHARE V2.1, has been assessed against cold leg LOCA test. The calculation results show that CATHARE predicts the break flow rate accurately for low pressure and highly subcooled inlet conditions. And CATHARE is fully able to calculate the overall scenario observed in the cold leg LOCA test. The calculated parameters are in acceptable agreement with the experimental data.



---

*Source: 179235-2015-01-26.xml*