# Study on Thermal-Hydraulic Behavior of an Integral Type Reactor under Heaving Condition

Authors: Beibei Feng; Xingtuan Yang; Shengqiang Li; Yanfei Sun; Jiyuan Tu; Shengyao Jiang
Journal: Science and Technology of Nuclear Installations
 (2014)

---

## Abstract

A self-developed program was used to study the thermal-hydraulic behavior of an integral type reactor under heaving condition. Comparison of calculated results with the data of experiments performed on a natural circulation loop designed with reference to an integral type reactor of Tsinghua University in inclination, heaving, and rolling motions was carried out. Characteristics of natural circulation in heaving motion and effect of motion parameters on natural circulation were investigated. Results indicated that: (1) long-period heaving motion would lead to more significant influence than inclination and rolling motion; (2) it was an alternating force field which consisted of gravity and an additional force that decided the flow temperature and density difference of natural circulation; (3) effect of strengthk and cycle T of heaving motion on flow fluctuation of natural circulation and condensate depression of heating section outlet was performed.


---

## Body

## 1. Introduction

Integral type reactor with its coolant system components in reactor vessel has wide applications due to its passive safety property. 5 MW experimental low temperature nuclear heating reactor (NHR) of Tsinghua University is a low-power integral type reactor with modular and self-stabilization which has been applied to the engineering fields of heat supply and seawater desalination [1, 2]. As part of the studies responding to the application of NHR as ship power, Institute of Nuclear and New Energy Technology (INET) of Tsinghua University built an experimental loop with reference to NHR to investigate the effect of rolling motion and inclination motion (two typical oceanic movements) on natural circulation [3–5]. A program was also developed to provide a tool for investigating the effect of heaving motion (the third typical oceanic movement) on natural circulation which could not realized through the experimental loop.On land, successful operation of NHR has proved the feasibility of its technical solutions. Under oceanic conditions, it needs revalidation through experiments resulting from the effect of oceanic environment on thermal-hydraulic behavior of natural circulation due to the additional movements caused by oceanic motions, such as rolling, heaving, and inclination. So, INET built a rolling experimental apparatus with reference to NHR to simulate oceanic motions, which consisted of a symmetrical two-circuit test loop and auxiliary equipment [6, 7]. Geometric parameters of experimental rolling apparatus was designed completely the same as NHR to investigate the effect of same strength of oceanic motion on natural circulation, and operation parameters were arranged the same as NHR. Due to the restriction of experimental rolling apparatus, effect of heaving motion on thermal-hydraulic behavior of natural circulation could not be investigated experimentally and a self-developed program was accordingly developed to accomplish this work. The aim of the present paper was to describe our comprehensive analysis results concerning natural circulation under heaving motion and effect of motion parameters on natural circulation obtained using the self-developed program. Comparison of effect of inclination motion and rolling motion on natural circulation between calculation and experimental data was conducted to validate the self-developed program. A heaving strength parameter k was proposed to describe heaving motion in alternating force field. Flow rate and condensate depression of natural circulation under heaving motion with various cycle T and strength k were investigated in detail.


## 2. Experimental System

### 2.1. Similarity Analysis

For two similar physical phenomena, parameters describing the law of the phenomena should be proportionable. Similarity analysis is carried out for catching the dynamics of the coupled physics. In order to guarantee a more clear and comparable comprehension to NHR, mathematical equations of experimental system should have the same solution as NHR. Criteria numbers deduced from momentum equation, energy equation, and thermal boundary condition could be obtained as follows:(1)
S
t
=
L
V
·
t
,
F
r
=
V
2
g
·
L
,
E
u
=
P
ρ
·
V
2
,
R
e
=
ρ
V
L
μ
,
P
e
=
ρ
c
V
L
λ
,
N
u
=
∂
L
λ
.For the issue of single-phase flow and heat transfer process, the experimental apparatus and NHR are equivalent when the above criteria numbers are the same. For the heat exchange process of natural circulation, density difference induced by flow temperature and height difference drive the natural circulation. So, Grashof number should be introduced:(2)
G
r
=
g
β
Δ
T
L
3
v
2
.


### 2.2. Experimental Apparatus

The experimental setup is made up of three loops (as shown in Figure1(a)). Test section in experiments is a symmetrical two-circuit loop (primary loop) which is designed exactly based on NHR’s structure provided with the same height and width. Primary loop equals to a cut from NHR in the diameter direction, as shown in Figure 1(b). Experimental condition parameters are arranged the same as NHR to simulate the actual environment and achieve high comparative accuracy.Experimental apparatus [6].
(a)
(b)
(c)Primary loop consisted of heating section, upper plenum, steam generator, downcomer, and lower plenum. A rolling test bench is built to fix primary loop and oceanic motions, such as rolling and inclination, could be achieved through the movement of the test bench driven by two motors. Three electric heaters (EH) located at the shared part of the two circuits supply the heat source with outlet temperature 220°C. To meet the requirement of heating power, heater tubes are divided into two groups of eight tubes each with 75 V DC operating voltage and 3.7 Kw single-robust power would be achieved. The arrangement of heater tubes and electrifying mode is shown in Figure1(c). Two double-pipe heat interchangers located on the upper part of each circuit were used, but the amount of heat exchange tubes was reduced proportionality to ensure same heat load on heat exchange area as NHR. The primary loop was mounted on a ship motions simulation platform and could achieve rolling motion with maximum angles of 45° and 20° in x direction and in y direction, respectively. The working fluid of the primary loop was demineralized water with a pressure of 4.0 MPa.The secondary loop consisted of a pump, a heat exchanger, two mixers, and a pressurizer. The coolant for the secondary loop was water driven by a pump rather than natural circulation. The mixer and control valves were used to equally set the coolant temperature and flow rates of the two circuits in order to avoid an unbalanced load. The secondary loop was kept at high pressure (4.0 MPa) to prevent the occurrence of two-phase flow due to the simplicity of the equipment. The pressurizer was filled with high-pressure nitrogen gas in order to adjust and maintain the pressure. The main equipment for the third loop was the cooling tower, through which heat input from the secondary loop was finally emitted into the environment with the coolant circulation.Main shape and operation parameters are listed in Table1. Compared with public experimental studies, testing section in this work is designed with following characteristics.(a)
A symmetrical two-circuit loop which simulate natural circulation of an actual reactor under oceanic motions is more reasonable.
(b)
Same height and width of primary loop as NHR will reflect the real strength of rolling motion.
(c)
Operation parameters of primary loop are arranged the same as NHR to simulate NHR’s operation.Table 1
Comparison of operation parameters of experimental apparatus and NHR.
Power (Kw)
Pressure bearing capacity(MPa)
Inlet/outlet temperature of heating section    (°C)
Experimental apparatus
Operation value
163.054
5.0
179.8/210.1
Design value
164
5.0
180/211
NHR
3 × 104
5.0
180/211
Temperature difference between inlet and outlet of heating section (°C)
Average temperature of heating section   (°C)
Outlet velocity of heating section  (m·s−1)
Experimental apparatus
Operation value
30.4
194.95
0.3619
Design value
31
195.5
0.37321
NHR
31
195.5
/
Outlet velocity of heat exchanger  (m·s−1)
Natural circulation flow(kg·s−1)
Amount of fuel rod
Experimental apparatus
Operation value
0.4292
1.2062
48
Design value
0.4318
1.2404
/
NHR
/
227
8784
Amount of heat exchangers
Type of heat exchanger
Experimental apparatus/NHR ratio
Experimental apparatus
Operation value
2
Telescopic
1/183
Design value
/
/
/
NHR
18
Telescopic
/


## 2.1. Similarity Analysis

For two similar physical phenomena, parameters describing the law of the phenomena should be proportionable. Similarity analysis is carried out for catching the dynamics of the coupled physics. In order to guarantee a more clear and comparable comprehension to NHR, mathematical equations of experimental system should have the same solution as NHR. Criteria numbers deduced from momentum equation, energy equation, and thermal boundary condition could be obtained as follows:(1)
S
t
=
L
V
·
t
,
F
r
=
V
2
g
·
L
,
E
u
=
P
ρ
·
V
2
,
R
e
=
ρ
V
L
μ
,
P
e
=
ρ
c
V
L
λ
,
N
u
=
∂
L
λ
.For the issue of single-phase flow and heat transfer process, the experimental apparatus and NHR are equivalent when the above criteria numbers are the same. For the heat exchange process of natural circulation, density difference induced by flow temperature and height difference drive the natural circulation. So, Grashof number should be introduced:(2)
G
r
=
g
β
Δ
T
L
3
v
2
.


## 2.2. Experimental Apparatus

The experimental setup is made up of three loops (as shown in Figure1(a)). Test section in experiments is a symmetrical two-circuit loop (primary loop) which is designed exactly based on NHR’s structure provided with the same height and width. Primary loop equals to a cut from NHR in the diameter direction, as shown in Figure 1(b). Experimental condition parameters are arranged the same as NHR to simulate the actual environment and achieve high comparative accuracy.Experimental apparatus [6].
(a)
(b)
(c)Primary loop consisted of heating section, upper plenum, steam generator, downcomer, and lower plenum. A rolling test bench is built to fix primary loop and oceanic motions, such as rolling and inclination, could be achieved through the movement of the test bench driven by two motors. Three electric heaters (EH) located at the shared part of the two circuits supply the heat source with outlet temperature 220°C. To meet the requirement of heating power, heater tubes are divided into two groups of eight tubes each with 75 V DC operating voltage and 3.7 Kw single-robust power would be achieved. The arrangement of heater tubes and electrifying mode is shown in Figure1(c). Two double-pipe heat interchangers located on the upper part of each circuit were used, but the amount of heat exchange tubes was reduced proportionality to ensure same heat load on heat exchange area as NHR. The primary loop was mounted on a ship motions simulation platform and could achieve rolling motion with maximum angles of 45° and 20° in x direction and in y direction, respectively. The working fluid of the primary loop was demineralized water with a pressure of 4.0 MPa.The secondary loop consisted of a pump, a heat exchanger, two mixers, and a pressurizer. The coolant for the secondary loop was water driven by a pump rather than natural circulation. The mixer and control valves were used to equally set the coolant temperature and flow rates of the two circuits in order to avoid an unbalanced load. The secondary loop was kept at high pressure (4.0 MPa) to prevent the occurrence of two-phase flow due to the simplicity of the equipment. The pressurizer was filled with high-pressure nitrogen gas in order to adjust and maintain the pressure. The main equipment for the third loop was the cooling tower, through which heat input from the secondary loop was finally emitted into the environment with the coolant circulation.Main shape and operation parameters are listed in Table1. Compared with public experimental studies, testing section in this work is designed with following characteristics.(a)
A symmetrical two-circuit loop which simulate natural circulation of an actual reactor under oceanic motions is more reasonable.
(b)
Same height and width of primary loop as NHR will reflect the real strength of rolling motion.
(c)
Operation parameters of primary loop are arranged the same as NHR to simulate NHR’s operation.Table 1
Comparison of operation parameters of experimental apparatus and NHR.
Power (Kw)
Pressure bearing capacity(MPa)
Inlet/outlet temperature of heating section    (°C)
Experimental apparatus
Operation value
163.054
5.0
179.8/210.1
Design value
164
5.0
180/211
NHR
3 × 104
5.0
180/211
Temperature difference between inlet and outlet of heating section (°C)
Average temperature of heating section   (°C)
Outlet velocity of heating section  (m·s−1)
Experimental apparatus
Operation value
30.4
194.95
0.3619
Design value
31
195.5
0.37321
NHR
31
195.5
/
Outlet velocity of heat exchanger  (m·s−1)
Natural circulation flow(kg·s−1)
Amount of fuel rod
Experimental apparatus
Operation value
0.4292
1.2062
48
Design value
0.4318
1.2404
/
NHR
/
227
8784
Amount of heat exchangers
Type of heat exchanger
Experimental apparatus/NHR ratio
Experimental apparatus
Operation value
2
Telescopic
1/183
Design value
/
/
/
NHR
18
Telescopic
/


## 3. Calculation Model and Validation

### 3.1. Mathematical Model

Mathematical model based on the analysis of experimental system was built first [8–10]. To simplify handling natural circulation in heaving condition, coordinate system was established on the experimental apparatus. It was mathematically an issue in noninertial reference frame of natural circulation under heaving condition. An inertia force was brought into the fluid particle momentum equation in heaving condition resulting from the additional force of oceanic motions.Continuity equation is as follows:(3)
∂
ρ
∂
t
+
∇
·
(
ρ
V
→
)
=
0
.Momentum equation is as follows:(4)
D
V
→
D
t
=
f
→
-
1
ρ
∇
P
+
1
ρ
∇
·
(
2
μ
S
)
.Energy equation is as follows:(5)
ρ
C
D
T
D
t
=
∇
·
(
λ
∇
T
)
,
where ρ, f, μ, and λ were the density, mass force, viscosity, and thermal conductivity of the fluid and S was rate of strain tensor. 								(6)
f
→
=
g
→
-
a
→
o
-
w
˙
→
×
r
→
-
w
→
×
(
w
→
×
r
→
)
-
2
w
→
×
V
r
→
,
where g was acceleration of gravity and a
→
o
=
d
2
r
o
→
/
d
t
2 as an additional force induced by the translational motion of experimental apparatus in x, y, and z direction. a
→
t
=
w
˙
→
×
r
→ was tangential force; a
→
n
=
w
˙
→
×
(
w
˙
→
×
r
→
) was normal force; and a
→
c
=
2
w
˙
→
×
V
→
r was Coriolis force.In the case of one-dimensional flow, the Coriolis force was negligible. Assume thata
o
z
=
k
g
·
sin
(
w
d
·
t
); thus g
-
a
o
z
=
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
·
g.Heaving motion was the vertical translational motion which was in the same direction of gravity. So,(
g
-
a
o
z
) could be treated as alternating gravitational field which was viewed as the product of acceleration of gravity g and a time-varying coefficient (
1
-
k
sin
(
w
d
·
t
)
).Relative velocity of natural circulation to experimental apparatus was used to describe natural circulation in the above formulas. Symmetric double circulation of primary loop that consisted of two inner circuits, one outer circuit, and three mediate channels was observed, as shown in Figure2. F101, F102, and F103 were used to mark the right heating section, the mediate heating section and the left heating section of primary loop, respectively, and F104, F105 were used to mark the two-side heat exchangers.Figure 2
Structure of primary loop consisted of inner circuit and outer circuit.


### 3.2. Program Verification

Verification of the developed program was carried out through comparison between calculation results and the data of experiments on land and in rolling, inclination conditions. For one-dimensional natural circulation, the key factor was whether the calculated flow resistance was consistent with the experimental results. So, calculated flow resistance would be modified through experiments to make sure of the correction of program. Then, use this program to calculate natural circulation in heaving condition.Figures3 and 4 plotted the comparison of experimental results and calculated results of rolling motion with swing angle ψ
m
=
22.5
° and rolling cycle T
=
13 s in zero-power condition and in hot state, respectively. It can be inferred that curves of calculated results were almost the same as that of experimental results of rolling motion with swing angle ψ
m
=
22.5
° and rolling cycle T
=
13 s in zero-power condition and of rolling motion with swing angle ψ
m
=
45
° and rolling cycle T
=
18 s in full-power condition.Comparison of experimental results and calculated results of rolling motion with swing angleψ
m = 22.5° and rolling cycle T
=
13 s in zero-power condition: (a) calculated result; (b) experimental result.
(a)
(b)Comparison of experimental results and calculated results of rolling motion with swing angleψ
m = 45° and rolling cycle T
=
18 s in full power condition: (a) calculated result; (b) experimental results.
(a)
(b)


## 3.1. Mathematical Model

Mathematical model based on the analysis of experimental system was built first [8–10]. To simplify handling natural circulation in heaving condition, coordinate system was established on the experimental apparatus. It was mathematically an issue in noninertial reference frame of natural circulation under heaving condition. An inertia force was brought into the fluid particle momentum equation in heaving condition resulting from the additional force of oceanic motions.Continuity equation is as follows:(3)
∂
ρ
∂
t
+
∇
·
(
ρ
V
→
)
=
0
.Momentum equation is as follows:(4)
D
V
→
D
t
=
f
→
-
1
ρ
∇
P
+
1
ρ
∇
·
(
2
μ
S
)
.Energy equation is as follows:(5)
ρ
C
D
T
D
t
=
∇
·
(
λ
∇
T
)
,
where ρ, f, μ, and λ were the density, mass force, viscosity, and thermal conductivity of the fluid and S was rate of strain tensor. 								(6)
f
→
=
g
→
-
a
→
o
-
w
˙
→
×
r
→
-
w
→
×
(
w
→
×
r
→
)
-
2
w
→
×
V
r
→
,
where g was acceleration of gravity and a
→
o
=
d
2
r
o
→
/
d
t
2 as an additional force induced by the translational motion of experimental apparatus in x, y, and z direction. a
→
t
=
w
˙
→
×
r
→ was tangential force; a
→
n
=
w
˙
→
×
(
w
˙
→
×
r
→
) was normal force; and a
→
c
=
2
w
˙
→
×
V
→
r was Coriolis force.In the case of one-dimensional flow, the Coriolis force was negligible. Assume thata
o
z
=
k
g
·
sin
(
w
d
·
t
); thus g
-
a
o
z
=
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
·
g.Heaving motion was the vertical translational motion which was in the same direction of gravity. So,(
g
-
a
o
z
) could be treated as alternating gravitational field which was viewed as the product of acceleration of gravity g and a time-varying coefficient (
1
-
k
sin
(
w
d
·
t
)
).Relative velocity of natural circulation to experimental apparatus was used to describe natural circulation in the above formulas. Symmetric double circulation of primary loop that consisted of two inner circuits, one outer circuit, and three mediate channels was observed, as shown in Figure2. F101, F102, and F103 were used to mark the right heating section, the mediate heating section and the left heating section of primary loop, respectively, and F104, F105 were used to mark the two-side heat exchangers.Figure 2
Structure of primary loop consisted of inner circuit and outer circuit.


## 3.2. Program Verification

Verification of the developed program was carried out through comparison between calculation results and the data of experiments on land and in rolling, inclination conditions. For one-dimensional natural circulation, the key factor was whether the calculated flow resistance was consistent with the experimental results. So, calculated flow resistance would be modified through experiments to make sure of the correction of program. Then, use this program to calculate natural circulation in heaving condition.Figures3 and 4 plotted the comparison of experimental results and calculated results of rolling motion with swing angle ψ
m
=
22.5
° and rolling cycle T
=
13 s in zero-power condition and in hot state, respectively. It can be inferred that curves of calculated results were almost the same as that of experimental results of rolling motion with swing angle ψ
m
=
22.5
° and rolling cycle T
=
13 s in zero-power condition and of rolling motion with swing angle ψ
m
=
45
° and rolling cycle T
=
18 s in full-power condition.Comparison of experimental results and calculated results of rolling motion with swing angleψ
m = 22.5° and rolling cycle T
=
13 s in zero-power condition: (a) calculated result; (b) experimental result.
(a)
(b)Comparison of experimental results and calculated results of rolling motion with swing angleψ
m = 45° and rolling cycle T
=
18 s in full power condition: (a) calculated result; (b) experimental results.
(a)
(b)


## 4. Results and Discussion

### 4.1. Characteristics of Natural Circulation in Heaving Motion

Reference [8] has reported the experimental verification of natural circulation of an actual reactor on land. The program has been proved to be correct through the comparison with experimental data. Comparison of flow characteristics of natural circulation on land with heaving motion was conducted to study the effect of heaving motion on natural circulation.The additional force induced by heaving motion was acting on natural circulation in a parallel field of force. Effect of heaving motion exhibited a symmetry property for a symmetrical natural circulation loop. Figure5 plotted the flow and condensate depression of channel outlet of the left and the intermediate heating section in heaving motion with strength (
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s. Dotted line in Figure 5 was the dynamic time-average value in steady heaving motion. Figure 6 plotted the flux of primary side of heat exchanger and condensate depression of outlet of secondary side.Flow and condensate depression of channel outlet of the left and the intermediate heating section in heaving motion with strength(
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s: (a) normalized flow; (b) condensate depression.
(a)
(b)Flux of primary side of heat exchanger and condensate depression of outlet of secondary side: (a) normalized flow; (b) condensate depression.
(a)
(b)Experimental results indicated the following.(1)
Flux of natural circulation, outlet temperature of heating section, and temperature of secondary side of the heat exchanger will rise and fall in the same period as heaving motion. The value of fluctuate cycle was about 13 s as shown in Figures5 and 6.
(2)
The flux and flow temperature of left heating section were the same as that of the intermediate heating section which indicated that the uniform distribution of flow was realized on the bottom of reactor core. Slight difference of flux and temperature were generated from the resistance of channel bend.
(3)
In heaving motion, time-average value of relative flux of heating section and the part of heat exchangers was both less than 1 which indicated a decreased capability of natural circulation of the primary loop. The amount of reduction was about 5% in heaving motion with strength(
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s.


### 4.2. Mechanism of Natural Circulation in Heaving Motion

Heaving motion of the experimental apparatus would result in an additional force acting on the natural circulation which was parallel to gravity. Fixing the coordinate system on the experimental apparatus, natural circulation in heaving motion could be the issue in the noninertia coordinate system and treated as natural circulation under alternating vertical force field which could be calculated through the following formula:(7)
g
-
a
o
z
=
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
·
g
.The driving force of natural circulation could be expressed as(8)
Δ
P
d
=
Δ
ρ
(
t
)
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
,
where Δ
ρ
(
t
) was the time-varying density difference between the hot fluid and the cold fluid. The driving force of natural circulation in heaving motion included the gravity and an additional force resulting from the heaving motion. Slight flow density difference in the primary loop could be achieved in the case of low level of temperature fluctuation caused by heaving motion. If the above factor was ignored, flow rate of natural circulation in heaving motion could be equal to the superposition of the steady flow and the additional flow resulting from the heaving motion. Gravity and the induced force were independent and the effect on natural circulation could be superposed of the two above force directly. In the case of stronger flow temperature difference resulting from heaving motion, the former superposition was not correct in this case any more. Compared with the steady natural circulation on land, flow fluctuation in heaving motion would result in flow density difference. Simultaneously, flow density difference affects the driving force of natural circulation. Consequently, natural circulation in strong heaving motion was a result of alternating force field and flow density field.Figure7 plotted the relationship chain of factors affecting natural circulation in heaving motion: the alternating field which consisted of gravity and additional force would drive natural circulation. Flow fluctuation resulting from the alternating field would cause the flow temperature fluctuation which resulted in flow density difference. Further, density difference simultaneously affected the driving force of natural circulation. Accordingly, flux of natural circulation changed. Relationship of the factors affecting natural circulation was illustrated in Figure 7.Figure 7
Relationship chain of factors affecting natural circulation in heaving motion.


### 4.3. Effect of Heaving Cycle and Heaving Strength on Natural Circulation

Considering the effect of heaving motion on natural circulation, a heaving itema
o
z
=
k
·
g
·
sin
(
w
d
·
t
) was imported into the momentum equation, where  w
d
=
2
π
/
T. So, flow velocity and displacement would be calculated through the following equations:
(9)
V
o
z
=
-
k
g
w
¯
d
·
k
·
g
·
cos
⁡
(
w
¯
d
·
t
)
,
A
o
z
=
-
k
g
w
¯
d
·
k
·
g
·
s
in
(
w
¯
d
·
t
)
.It can be inferred that phase of displacement of heaving motion was inverted to acceleration and the amplitude was proportional to frequency squared. The parameterk was proposed to describe the strength of heaving motion.Figure8 plotted the flux of natural circulation, condensate depression of heating section’s outlet, and the secondary side of heat exchanger under heaving motion with cycle T
=
3 s and k
=
0.2
,
0.4
,
0.6, respectively. Results indicated the following.(1)
Fluctuation of flux and temperature would be stronger with the increasing heaving strength. When heaving strength(
g
-
a
o
z
)
=
0
.
6 g and cycle T
=
3 s, fluctuation of flux was up to 17%.
(2)
Changes of the condensate depression of heating section’s outlet exhibited regular fluctuation along with the heaving motion experiencing several heaving motion cycles with the start of steady state on land, as shown in Figure8(b). In the initial stage, the inlet temperature of heating section kept constant while about 17% flow fluctuation was observed in the process of fluid flowing into heating section channel from the down comer. Accordingly, condensate depression of heating section’s outlet exhibited serious fluctuation in the first heaving cycle. Weaker temperature fluctuation of heating section’s outlet and the secondary side of heat exchanger was achieved in the later transitory stage. This was because the temperature of heating section’s outlet was decided by inlet temperature and flux fluctuation which would go weaker with the end of transitory stage.Flux of natural circulation, condensate depression of heating section’s outlet, and the secondary side of heat exchanger under heaving motion with cycleT
=
3 s and k
=
0.2, 0.4, and 0.6: (a) flux of natural circulation; (b) condensate depression of heating section’s outlet; (c) secondary side of heat exchanger.
(a)
(b)
(c)Whilek was constant, the flow and temperature of natural circulation would exhibit similar fluctuation with the variational cycles of heaving motion. In Figure 9, the peak and valley of natural circulation were given. Flow fluctuation would rise nonlinearly with the heaving cycle T changing from 3 s to 8 s when heaving strength k was constant. Longer heaving cycle meant that the additional force generating from heaving motion would act on natural circulation in a longer period of time which led to bigger flow change and thermal saturation, or even boiling. Figure 10 plotted the calculated results of condensate depression of heating section outlet in heaving motion with strength k
=
0.6 and cycle T
=
3 s, 8s, 13s, 18s, and 23s, respectively. Results indicated an obvious fluctuation of temperature with heaving cycle changing from 3s to 8s. When heaving cycle T
=
8 s, the temperature difference from peak to valley was about 16°C. The critical point of heaving cycle was T
=
13 s when heaving strength k
=
0.6: there was no thermal saturation when T
<
13 s, but thermal saturation when T
>
13 s.Figure 9
Peak and valley of flow changes of natural circulation with variable heaving cycleT when heaving strength k
=
0.6.Figure 10
Calculated results of condensate depression of heating section outlet in heaving motion with strengthk
=
0.6 and cycle T
=
3 s, 8 s, 13 s, 18 s, and 23 s.Little effect of short‐cycle heaving motion (for instance T
=
3 s) on natural circulation was observed, while significant effect, or even serious accident, would occur under heaving motion with long cycle.


### 4.4. Comparison of Characteristics of Natural Circulation in Heaving, Inclination, and Rolling Motions

Effect of basic oceanic movements, such as heaving motion, inclination, and rolling motion on natural circulation was investigated through the program and experiments. Comparison of experimental results and simulation results of natural circulation under inclination was carried out, as shown in Figure11. Total amount of natural circulation decreased slightly, while the flow going through one heat exchanger would increase, but another would decrease under inclination condition, which was related to the resistance characteristics, geometric parameters, and facility arrangement of the experimental apparatus. Under inclination with tilt angle 45°, 12% reduction of the total amount of natural circulation was achieved. From Figure 12, little effect of inclination on natural circulation was observed through comparison of calculated results with experimental data.Figure 11
Comparison of experimental results and simulation results of natural circulation in inclination condition with different angles.Effect of heaving motion with strength 0.6 g and cycle 8 s, vertical rolling with swing angle 22.5° and cycle 8 s, and compound motions consisted of the above two movements on the normalized flow of natural circulation.
(a)
(b)Effect of heaving motion with strength 0.6 g and cycle 8 s, vertical rolling with swing angle 22.5° and cycle 8 s, compound motions which consisted of the above two movements, heaving motion with strength 0.6 g and cycle 8 s, transverse rolling with swing angle 22.5° and cycle 13 s, and compound motion which consisted of the above two movements on the normalized flow of natural circulation was shown in Figure12. Rolling motion along the axis parallel to this page was defined as vertical rolling motion, and rolling motion along the axis perpendicular to this page was defined as transverse rolling motion. An increase of flow fluctuation was observed compared heaving motion 0.6 g, 8 s with vertical rolling 13 s, 22.5°. And, stronger effect of the compound motion consisted of the mentioned heaving motion and rolling motion on natural circulation was confirmed. Due to the identical cycle of vertical rolling motion and heaving motion, cycle of natural circulation in the combination of the above two movements was synchronous. Considering the effect of heaving motion and transverse rolling, much stronger fluctuation of natural circulation in heaving motion with strength 0.6 g and cycle 8 s was observed than in transverse rolling motion with swing angle 22.5° and cycle 13 s. Stronger fluctuation could be predicted in heaving motion with cycle 13 s. A complex flow waveform resulting from the motion of heaving together with transverse rolling was observed and the heaving component was considered the main driving factor. Same flow fluctuation cycle was observed as the cycle of heaving component, and magnitude of flow fluctuation would be affected by the transverse component.For large-scale ships, oceanic low frequency motion would lead to significant effect on natural circulation. Long-period rolling motion took slight effect on natural circulation due to small amount of the induced additional force which inclination could be treated as a special type of rolling motion with infinite cycle. But long-period heaving motion would lead to significant effect on natural circulation. Through the comparison, we could find out that long-period heaving motion may lead to the worst effect on natural circulation than inclination and rolling motion.


## 4.1. Characteristics of Natural Circulation in Heaving Motion

Reference [8] has reported the experimental verification of natural circulation of an actual reactor on land. The program has been proved to be correct through the comparison with experimental data. Comparison of flow characteristics of natural circulation on land with heaving motion was conducted to study the effect of heaving motion on natural circulation.The additional force induced by heaving motion was acting on natural circulation in a parallel field of force. Effect of heaving motion exhibited a symmetry property for a symmetrical natural circulation loop. Figure5 plotted the flow and condensate depression of channel outlet of the left and the intermediate heating section in heaving motion with strength (
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s. Dotted line in Figure 5 was the dynamic time-average value in steady heaving motion. Figure 6 plotted the flux of primary side of heat exchanger and condensate depression of outlet of secondary side.Flow and condensate depression of channel outlet of the left and the intermediate heating section in heaving motion with strength(
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s: (a) normalized flow; (b) condensate depression.
(a)
(b)Flux of primary side of heat exchanger and condensate depression of outlet of secondary side: (a) normalized flow; (b) condensate depression.
(a)
(b)Experimental results indicated the following.(1)
Flux of natural circulation, outlet temperature of heating section, and temperature of secondary side of the heat exchanger will rise and fall in the same period as heaving motion. The value of fluctuate cycle was about 13 s as shown in Figures5 and 6.
(2)
The flux and flow temperature of left heating section were the same as that of the intermediate heating section which indicated that the uniform distribution of flow was realized on the bottom of reactor core. Slight difference of flux and temperature were generated from the resistance of channel bend.
(3)
In heaving motion, time-average value of relative flux of heating section and the part of heat exchangers was both less than 1 which indicated a decreased capability of natural circulation of the primary loop. The amount of reduction was about 5% in heaving motion with strength(
g
-
a
o
z
)
=
0.4 g and cycle T
=
13 s.


## 4.2. Mechanism of Natural Circulation in Heaving Motion

Heaving motion of the experimental apparatus would result in an additional force acting on the natural circulation which was parallel to gravity. Fixing the coordinate system on the experimental apparatus, natural circulation in heaving motion could be the issue in the noninertia coordinate system and treated as natural circulation under alternating vertical force field which could be calculated through the following formula:(7)
g
-
a
o
z
=
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
·
g
.The driving force of natural circulation could be expressed as(8)
Δ
P
d
=
Δ
ρ
(
t
)
(
1
-
k
sin
⁡
(
w
d
·
t
)
)
,
where Δ
ρ
(
t
) was the time-varying density difference between the hot fluid and the cold fluid. The driving force of natural circulation in heaving motion included the gravity and an additional force resulting from the heaving motion. Slight flow density difference in the primary loop could be achieved in the case of low level of temperature fluctuation caused by heaving motion. If the above factor was ignored, flow rate of natural circulation in heaving motion could be equal to the superposition of the steady flow and the additional flow resulting from the heaving motion. Gravity and the induced force were independent and the effect on natural circulation could be superposed of the two above force directly. In the case of stronger flow temperature difference resulting from heaving motion, the former superposition was not correct in this case any more. Compared with the steady natural circulation on land, flow fluctuation in heaving motion would result in flow density difference. Simultaneously, flow density difference affects the driving force of natural circulation. Consequently, natural circulation in strong heaving motion was a result of alternating force field and flow density field.Figure7 plotted the relationship chain of factors affecting natural circulation in heaving motion: the alternating field which consisted of gravity and additional force would drive natural circulation. Flow fluctuation resulting from the alternating field would cause the flow temperature fluctuation which resulted in flow density difference. Further, density difference simultaneously affected the driving force of natural circulation. Accordingly, flux of natural circulation changed. Relationship of the factors affecting natural circulation was illustrated in Figure 7.Figure 7
Relationship chain of factors affecting natural circulation in heaving motion.


## 4.3. Effect of Heaving Cycle and Heaving Strength on Natural Circulation

Considering the effect of heaving motion on natural circulation, a heaving itema
o
z
=
k
·
g
·
sin
(
w
d
·
t
) was imported into the momentum equation, where  w
d
=
2
π
/
T. So, flow velocity and displacement would be calculated through the following equations:
(9)
V
o
z
=
-
k
g
w
¯
d
·
k
·
g
·
cos
⁡
(
w
¯
d
·
t
)
,
A
o
z
=
-
k
g
w
¯
d
·
k
·
g
·
s
in
(
w
¯
d
·
t
)
.It can be inferred that phase of displacement of heaving motion was inverted to acceleration and the amplitude was proportional to frequency squared. The parameterk was proposed to describe the strength of heaving motion.Figure8 plotted the flux of natural circulation, condensate depression of heating section’s outlet, and the secondary side of heat exchanger under heaving motion with cycle T
=
3 s and k
=
0.2
,
0.4
,
0.6, respectively. Results indicated the following.(1)
Fluctuation of flux and temperature would be stronger with the increasing heaving strength. When heaving strength(
g
-
a
o
z
)
=
0
.
6 g and cycle T
=
3 s, fluctuation of flux was up to 17%.
(2)
Changes of the condensate depression of heating section’s outlet exhibited regular fluctuation along with the heaving motion experiencing several heaving motion cycles with the start of steady state on land, as shown in Figure8(b). In the initial stage, the inlet temperature of heating section kept constant while about 17% flow fluctuation was observed in the process of fluid flowing into heating section channel from the down comer. Accordingly, condensate depression of heating section’s outlet exhibited serious fluctuation in the first heaving cycle. Weaker temperature fluctuation of heating section’s outlet and the secondary side of heat exchanger was achieved in the later transitory stage. This was because the temperature of heating section’s outlet was decided by inlet temperature and flux fluctuation which would go weaker with the end of transitory stage.Flux of natural circulation, condensate depression of heating section’s outlet, and the secondary side of heat exchanger under heaving motion with cycleT
=
3 s and k
=
0.2, 0.4, and 0.6: (a) flux of natural circulation; (b) condensate depression of heating section’s outlet; (c) secondary side of heat exchanger.
(a)
(b)
(c)Whilek was constant, the flow and temperature of natural circulation would exhibit similar fluctuation with the variational cycles of heaving motion. In Figure 9, the peak and valley of natural circulation were given. Flow fluctuation would rise nonlinearly with the heaving cycle T changing from 3 s to 8 s when heaving strength k was constant. Longer heaving cycle meant that the additional force generating from heaving motion would act on natural circulation in a longer period of time which led to bigger flow change and thermal saturation, or even boiling. Figure 10 plotted the calculated results of condensate depression of heating section outlet in heaving motion with strength k
=
0.6 and cycle T
=
3 s, 8s, 13s, 18s, and 23s, respectively. Results indicated an obvious fluctuation of temperature with heaving cycle changing from 3s to 8s. When heaving cycle T
=
8 s, the temperature difference from peak to valley was about 16°C. The critical point of heaving cycle was T
=
13 s when heaving strength k
=
0.6: there was no thermal saturation when T
<
13 s, but thermal saturation when T
>
13 s.Figure 9
Peak and valley of flow changes of natural circulation with variable heaving cycleT when heaving strength k
=
0.6.Figure 10
Calculated results of condensate depression of heating section outlet in heaving motion with strengthk
=
0.6 and cycle T
=
3 s, 8 s, 13 s, 18 s, and 23 s.Little effect of short‐cycle heaving motion (for instance T
=
3 s) on natural circulation was observed, while significant effect, or even serious accident, would occur under heaving motion with long cycle.


## 4.4. Comparison of Characteristics of Natural Circulation in Heaving, Inclination, and Rolling Motions

Effect of basic oceanic movements, such as heaving motion, inclination, and rolling motion on natural circulation was investigated through the program and experiments. Comparison of experimental results and simulation results of natural circulation under inclination was carried out, as shown in Figure11. Total amount of natural circulation decreased slightly, while the flow going through one heat exchanger would increase, but another would decrease under inclination condition, which was related to the resistance characteristics, geometric parameters, and facility arrangement of the experimental apparatus. Under inclination with tilt angle 45°, 12% reduction of the total amount of natural circulation was achieved. From Figure 12, little effect of inclination on natural circulation was observed through comparison of calculated results with experimental data.Figure 11
Comparison of experimental results and simulation results of natural circulation in inclination condition with different angles.Effect of heaving motion with strength 0.6 g and cycle 8 s, vertical rolling with swing angle 22.5° and cycle 8 s, and compound motions consisted of the above two movements on the normalized flow of natural circulation.
(a)
(b)Effect of heaving motion with strength 0.6 g and cycle 8 s, vertical rolling with swing angle 22.5° and cycle 8 s, compound motions which consisted of the above two movements, heaving motion with strength 0.6 g and cycle 8 s, transverse rolling with swing angle 22.5° and cycle 13 s, and compound motion which consisted of the above two movements on the normalized flow of natural circulation was shown in Figure12. Rolling motion along the axis parallel to this page was defined as vertical rolling motion, and rolling motion along the axis perpendicular to this page was defined as transverse rolling motion. An increase of flow fluctuation was observed compared heaving motion 0.6 g, 8 s with vertical rolling 13 s, 22.5°. And, stronger effect of the compound motion consisted of the mentioned heaving motion and rolling motion on natural circulation was confirmed. Due to the identical cycle of vertical rolling motion and heaving motion, cycle of natural circulation in the combination of the above two movements was synchronous. Considering the effect of heaving motion and transverse rolling, much stronger fluctuation of natural circulation in heaving motion with strength 0.6 g and cycle 8 s was observed than in transverse rolling motion with swing angle 22.5° and cycle 13 s. Stronger fluctuation could be predicted in heaving motion with cycle 13 s. A complex flow waveform resulting from the motion of heaving together with transverse rolling was observed and the heaving component was considered the main driving factor. Same flow fluctuation cycle was observed as the cycle of heaving component, and magnitude of flow fluctuation would be affected by the transverse component.For large-scale ships, oceanic low frequency motion would lead to significant effect on natural circulation. Long-period rolling motion took slight effect on natural circulation due to small amount of the induced additional force which inclination could be treated as a special type of rolling motion with infinite cycle. But long-period heaving motion would lead to significant effect on natural circulation. Through the comparison, we could find out that long-period heaving motion may lead to the worst effect on natural circulation than inclination and rolling motion.


## 5. Conclusion

The comparison of calculated results with experimental data on a natural circulation loop designed with reference to 5 MW experimental low temperature nuclear heating reactor of Tsinghua University in inclination, heaving, and rolling motions was carried out. Flow characteristics and effect of motion parameters on natural circulation in heaving motion were investigated and mechanism of natural circulation in heaving motion was stated. Results indicated the following.(1)
Significant effect of heaving motion on natural circulation was observed. Flow fluctuation would become stronger in heaving motion with long period and high strength.
(2)
Slight effect of short-period heaving motion on natural circulation was verified.
(3)
Relationship chain of factors affecting natural circulation in heaving motion was given, the alternating field consisted of gravity and additional force, and flow density difference would result in the flow fluctuation of natural circulation.
(4)
Long-period heaving motion would lead to more significant influence than inclination and rolling motion and may cause serious accident.



---

*Source: 857507-2014-04-08.xml*