# Effect of Heaving Movement on Flow Instability in U-Tubes of Marine Steam Generator under Natural Circulation

Authors: Jianli Hao; Wenzhen Chen; De Zhang
Journal: Science and Technology of Nuclear Installations
 (2013)

---

## Abstract

Under heaving movement conditions, the single phase flow instability in U-tubes is affected by the additional force, which will influence the marine reactor operation. In the present work, one-dimensional thermal-hydraulic model in U-tubes under heaving movement conditions is established, and the critical pressure drop (CPD) and critical mass flow rate (CMFR) which relate to the occurrence of reverse flow in U-tubes are proposed and analyzed. The effects of the heaving period and heaving acceleration amplitude on the flow instability in U-tubes with the different length are discussed. It is shown that (1) the CPD and CMFR are obviously affected by the heaving movement, which means that the reverse flow characteristic in U-tubes will be changed; (2) the fluctuation periods of the CPD and CMFR are the same as the heaving period, but the fluctuation magnitude of them is little affected by the heaving period; (3) the relative changes of CPD and CMFR are the linear function of heaving acceleration amplitude; and (4) the U-tube length has little influence on the relative changes of CPD and CMFR compared with the heaving acceleration amplitude, which means that the heaving movement has little influence on the space distribution of reverse flow in the U-tubes of marine steam generator.


---

## Body

## 1. Introduction

Natural circulation operation is significant for the marine nuclear power plant in terms of passive safety, efficiency, and noise decrease. The operation performance of steam generator (SG) can deeply influence the nuclear power plant safety. Under natural circulation conditions, it is shown that the single phase flow in the parallel U-tubes of SG may be unstable, and reverse flow occurs within some U-tubes [1, 2]. Because of the occurrence of reverse flow, the effective heat transfer area of SG primary side is reduced, and the flow resistance coefficient under the natural circulation is obviously larger than that under the forced circulation. So the actual nature circulation capability of primary loop is lower than the needed value due to the reverse flow in U-tubes, which has negative influence on the operation of the marine nuclear power plant [3].The flow instability in U-tubes is considered to be a typical Ledinegg-type single phase flow instability [4]. Yang et al. [5] developed a lumped-distribution model to calculate the reverse flow in the inverted U-tubes. Walter and Linzer [6] discussed the influence of the operating pressure on the reverse flow in natural circulation system. They thought that the design procedure for natural circulation systems with variable operating pressure should include the analysis of flow stability. Sanders [7] and Hao et al. [8] derived the relationships between the pressure drop and velocity and mass flow rate. They found that the relationship curves had a negative slope zone, where the flow in U-tubes was unstable and an excursion would occur. The pressure drop and velocity or mass flow rate at the inflexion of the curves are named the critical pressure drop (CPD) and critical velocity (CV) or critical mass flow rate (CMFR) by Hao et al. [9, 10].Under the ocean conditions, the marine nuclear power plant has different characteristic from the land-based nuclear power plant due to the heaving and rolling movement [3], and the phenomena of reverse flow in U-tubes for marine steam generator are more complicated. Pendyala et al. [11] experimentally studied the mass flow rate and pressure drop in a vertical tube under the heaving movement; the results show that the magnitude of the mass flow rate fluctuation decreases with the increase in Reynolds number. While the Reynolds number is greater than 5000, the effect of heaving movement on the friction factor in the tube can be ignored. The heaving movement brings on the additional pressure drop, which has large influence on the total pressure drop in the tubes. So the heaving movement essentially affects the critical mass flow rate and critical pressure drop in U-tubes. In this paper, the relationships between the pressure drop and mass flow rate under heaving movement conditions are analyzed, and the flow instability in the U-tubes is discussed.


## 2. Theoretical Model

Under low frequency conditions, the heaving movement has little influence on the turbulent flow [11]. So the flow and heat transfer in the U-tubes can be simplified as follows: (1) it is acceptable that the one-dimensional approach is sufficient for modeling the flow in the U-tubes; (2) while the inlet mass flow rate is determined, the fluid velocity in the U-tubes is not affected by the heaving movement; (3) the water density in the U-tubes can be expressed as a linear function of temperature; (4) the water in the secondary side of the SG is in the saturated state. The schematic of U-tube under heaving movement conditions is shown in Figure 1.Figure 1
Schematic of U-tube under heaving movement conditions.The total pressure drop in U-tubes under heaving movement conditions can be written as follows:(1)Δp=Δpf+Δpg+Δpa,
where Δpf is the flow resistance pressure drop, Δpg is the gravitational pressure drop, and Δpa is the additional pressure drop caused by the heaving movement. Δpf, Δpg, and Δpa can be written as follows [12]:
(2)Δpf=(fld0+ζ)m˙22A2ρ-,(3)Δpg=-ΔρgH,(4)Δpa=-ΔρHa,
where m˙ is the mass flow rate, f is the friction resistance coefficient, f=0.3164Re-0.25, Re is the Reynolds number, l is the U-tube length, d0 is the U-tube inside diameter, ζ is the local resistance coefficient, ζ=0.262+0.326(d0/ru)3.5, ru is the U-tube bending radius, A is the U-tube flow area, ρ- is the average density, Δρ is the density difference between the cold side and hot side of the U-tube, g is the gravitational acceleration, H is the U-tube height, and a is the heaving acceleration.The heaving acceleration can be written as follows [3]:
(5)a=-a0sin(2πTt),
where a0 is the heaving acceleration amplitude, T is the heaving period, and t is the flow time.The additional pressure drop caused by the heaving movement can be written as follows:(6)Δpa=-Δρa0Hsin(2πTt).The temperature fluctuation caused by the heaving movement is little; the density change along the U-tube caused by the heaving movement conditions can be ignored [13]. Based on Boussinesq approximation, the density can be written as follows:
(7)ρ=ρ0(1-β(Tw-T0)),
where ρ0 is the reference density, β is the thermal expansion coefficient, Tw and T0 are the flow water temperature and reference temperature, respectively.The fluid energy conservation equation can be obtained as follows [3]:
(8)∂Tw∂s=-hspPm˙cp(Tw-T0),
where s is the coordinate in the normal flow direction along the U-tube. Substituting (7) into (8) results in
(9)∂ρ∂s=hspPm˙cp(ρ0-ρ).
While the inlet mass flow rate and fluid temperature are given, integrating (9) along the U-tube length leads to
(10)ρ=ρ0-(ρ0-ρin)e-hspPs/m˙cp,
where ρin is the inlet fluid density, cp is the specific heat capacity, P is the wetted perimeter of U-tube, hsp is the overall heat transfer coefficient, hsp=1/(1/h0+(d0/2λwall)ln(d1/d0)+d0/h1d1), h0 is the surface heat transfer coefficient of U-tube inside wall, λwall is the U-tube thermal conductivity, d1 is the U-tube outside diameter, and h1 is the surface heat transfer coefficient of U-tube outside wall.From (10), the expressions of ρ- and Δρ can be derived as follows:
(11)ρ-=ρ0-(ρ0-ρin)m˙cphspPl(1-e-hspPl/m˙cp),(12)Δρ=(ρ0-ρin)2m˙cphspPl(1-e-hspPl/2m˙cp)2.Inserting (11) into (2), Δpf can be obtained as follows:
(13)Δpf=(m˙2(fl/d0+ζ)2A2)×(ρ0-(ρ0-ρin)cpm˙hspPl(1-e-hspPl/m˙cp))-1.Inserting (12) into (3) and (4), respectively, Δpg and Δpa can be obtained as follows:
(14)Δpg=-(ρ0-ρin)cpm˙hspPl/2(1-e-hspPl/2m˙cp)2Hg,Δpa=-2(ρ0-ρin)cpm˙hspPl(1-e-hspPl/2m˙cp)2Ha0sin(2πTt).SubstitutingΔpf, Δpg, and Δpa into (1) results in
(15)Δp=(m˙2(fl/d0+ζ)2A2)×(ρ0-(ρ0-ρin)cpm˙hspPl(1-e-hspPl/m˙cp))-1-(ρ0-ρin)cpm˙hspPl/2×(1-e-hspPl/2m˙cp)2H(g+a0sin(2πTt)).While∂Δp/∂m˙<0, the flow in U-tubes is unstable [8]. The criterion equation of flow instability in U-tubes under the heaving movement is ∂Δp/∂m˙=0; then the following equation can be obtained:
(16)0=m˙cA2ρ0(fld0+ζ)-ρ0-ρinhspPcp×(1-e-hspPl/2m˙ccp)2H(g+a0sin2πTt)-2ρ0-ρinm˙c(g+a0sin2πTt)×H(-e-hspPl/2m˙ccp+e-hspPl/m˙ccp),
where m˙c is the critical mass flow rate (CMFR). Equation (16) is the implicit expression of m˙c, which can be solved by numerical iteration method.Substituting CMFR obtained by (16) into (15), the critical pressure drop (CPD) can be given as follows:
(17)Δpc=-m˙c22A2ρ0(fld0+ζ)=-2H(ρ0-ρin)e-hspPl/2m˙ccp[g+a0sin(2πTt)].While the mass flow rate (or total pressure drop) in U-tubes is lower thanm˙c (or Δpc), the flow in the U-tubes is unstable, and the reverse flow will occur.


## 3. Numerical Experiments

In order to validate the effect of heaving movement on the flow instability in U-tubes, a certain type of marine steam generator is selected [14], the parameters of the U-tubes are shown in Table 1, and l0 is the longest tube length of the marine steam generator. The inner and outside diameter and straight tube length of the U-tubes are the same, but the bending radiuses of the U-tubes are different.Table 1
Parameters of marine steam generator U-tubes.
Name of parameter
Value of parameter
Inner diameter of U-tubes, m
d
0
Length of U-tubes, m
0.7l0~l0
Secondary side operating pressure of steam generator, MPa
3.0
Inlet fluid temperature, K
T
in
Outside diameter of U-tubes, m
d
1
Primary side operating pressure of steam generator, MPa
14.0
Mass flow rate, kg/s
0.006~0.06
Gravitational acceleration, m/s2
9.8


### 3.1. Relationships between Total Pressure Drop and Mass Flow Rate under Heaving Movement

The heaving movement period is set to be 10 s, and the acceleration amplitude is 0.5 m/s2  [3]. The variation of the heaving movement acceleration with time is shown in Figure 2. Based on (15), the variations of total pressure drop with mass flow rate for the different heaving movement time are shown in Figure 3.Figure 2
Variation of heaving movement acceleration versus time.Figure 3
Variation of total pressure drop with mass flow rate for different time.From (4), it is known that the relationship between the additional pressure drop and heaving acceleration is linear. The additional pressure drop will decrease while the heaving acceleration decreases, which causes the decrease in the total pressure drop. It can be seen from Figure 3 that the heaving movement has a certain effect on the critical pressure drop (CPD) and critical mass flow rate (CMFR), which is discussed in the following section.


### 3.2. Effect of Heaving Period

Based on (16) and (17), the effects of heaving period on the critical pressure drop (CPD) and critical mass flow rate (CMFR) are analyzed, which are shown in Figures 4 and 5. Where the U-tube length is 0.7l0, the amplitude of heaving acceleration is 0.5 m/s2, and the heaving periods are 5, 10, 15, and 20 s, respectively.Figure 4
Variation of critical mass flow rate with time for different heaving period.Figure 5
Variation of critical pressure drop with time for different heaving period.From Figures4 and 5 it can be seen that the CPD and CMFR change periodically under the heaving movement and the CMFR increases at first and last quarter time, but the variation of CPD is opposite to it. The fluctuation period is the same as the heaving period, but the heaving periods have little influence on the fluctuation magnitude of the CMFR and CPD. The data in Figures 4 and 5 are fitted; then the variation of the CMFR and CPD with time can be expressed as follows:
(18)m˙c=m˙c,0[1+0.03294sin(2πTt)],(19)Δpc=Δpc,0[1+0.06sin(2πTt)],
where m˙c,0 and Δpc,0 are the CMFR and CPD in the steady state, respectively, m˙c,0=0.00979 kg/s, and Δpc,0=-314.58932 Pa.


### 3.3. Effect of Heaving Acceleration Amplitude

Here the U-tube length is chosen as0.7l0, the heaving period is 10 s, and the amplitudes of heaving acceleration are 0.5, 1.0, 1.5, and 2.0 m/s2, respectively. The effects of heaving acceleration amplitude on the flow instability are shown in Figures 6 and 7.Figure 6
Variation of critical mass flow rate with time for different amplitude of heaving acceleration.Figure 7
Variation of critical pressure drop with time for different amplitude of heaving acceleration.It can be seen from Figures6 and 7 that the fluctuation magnitude of CMFR and CPD increases with the increase in the amplitude of heaving acceleration. The data in Figures 6 and 7 are fitted, and the relative change values of CMFR and CPD can be written as follows:
(20)|m˙c,M-m˙c,0m˙c,0|=0.06679a0,(21)|Δpc,M-Δpc,0Δpc,0|=0.11998a0,
where m˙c,M and Δpc,M are the maximum values of CMFR and CPD in the heaving period, respectively.Combining (18) and (20), the variation of CMFR with time under the heaving movement can be obtained as follows:
(22)m˙c=m˙c,0[1+0.06679a0sin(2πTt)].Combining (19) and (21), the variation of CPD with time under the heaving movement can be obtained as follows:
(23)Δpc=Δpc,0[1+0.11998a0sin(2πTt)].


### 3.4. Effect of Heaving Movement on Flow Instability in U-Tubes with Different Length

Based on (16) and (17), the flow instability in U-tubes with different length under the heaving movement is investigated, and the results are shown in Figures 8–11.Figure 8
Variation of critical mass flow rate with time for different length.Figure 9
Variation of critical pressure drop with time for different length.Figure 10
Relative change of critical mass flow rate with U-tube length.Figure 11
Relative change of critical pressure drop with U-tube length.In Figures8 and 9, the heaving period is 10 s and the amplitude of heaving acceleration is 0.5 m/s2. l=0.7l0 and l=l0 are the shortest and longest length of U-tubes in the marine steam generator, respectively. It can be seen from Figure 8 that the CMFR of both the long tube (l=l0) and short tube (l=0.7l0) is lower than that of the U-tube with the middle length (l=0.8l0).It can be seen from Figure9 that the CPD of short tube (l=0.7l0) is higher than that of others, which means that the reverse flow will occur in the short tube easily. From Figures 8 and 9, it can be also found that the effects of heaving movement on the flow instability in U-tubes with different length are similar, which means that the heaving movement has little influence on the space distribution of reverse flow.Figures10 and 11 give the relative changes of the CMFR and CPD with U-tube length. Where the heaving period is 10 s and the amplitudes of heaving acceleration are 0.5, 1.0, 1.5, and 2.0 m/s2, respectively. It can be seen from Figures 10 and 11 that the amplitude of heaving acceleration has large influence on the relative changes of CMFR and CPD; namely, the CMFR and CPD increase with the increase in the heaving acceleration amplitude. However the effect of U-tube length on the CMFR is more obvious than that on the CPD.


## 3.1. Relationships between Total Pressure Drop and Mass Flow Rate under Heaving Movement

The heaving movement period is set to be 10 s, and the acceleration amplitude is 0.5 m/s2  [3]. The variation of the heaving movement acceleration with time is shown in Figure 2. Based on (15), the variations of total pressure drop with mass flow rate for the different heaving movement time are shown in Figure 3.Figure 2
Variation of heaving movement acceleration versus time.Figure 3
Variation of total pressure drop with mass flow rate for different time.From (4), it is known that the relationship between the additional pressure drop and heaving acceleration is linear. The additional pressure drop will decrease while the heaving acceleration decreases, which causes the decrease in the total pressure drop. It can be seen from Figure 3 that the heaving movement has a certain effect on the critical pressure drop (CPD) and critical mass flow rate (CMFR), which is discussed in the following section.


## 3.2. Effect of Heaving Period

Based on (16) and (17), the effects of heaving period on the critical pressure drop (CPD) and critical mass flow rate (CMFR) are analyzed, which are shown in Figures 4 and 5. Where the U-tube length is 0.7l0, the amplitude of heaving acceleration is 0.5 m/s2, and the heaving periods are 5, 10, 15, and 20 s, respectively.Figure 4
Variation of critical mass flow rate with time for different heaving period.Figure 5
Variation of critical pressure drop with time for different heaving period.From Figures4 and 5 it can be seen that the CPD and CMFR change periodically under the heaving movement and the CMFR increases at first and last quarter time, but the variation of CPD is opposite to it. The fluctuation period is the same as the heaving period, but the heaving periods have little influence on the fluctuation magnitude of the CMFR and CPD. The data in Figures 4 and 5 are fitted; then the variation of the CMFR and CPD with time can be expressed as follows:
(18)m˙c=m˙c,0[1+0.03294sin(2πTt)],(19)Δpc=Δpc,0[1+0.06sin(2πTt)],
where m˙c,0 and Δpc,0 are the CMFR and CPD in the steady state, respectively, m˙c,0=0.00979 kg/s, and Δpc,0=-314.58932 Pa.


## 3.3. Effect of Heaving Acceleration Amplitude

Here the U-tube length is chosen as0.7l0, the heaving period is 10 s, and the amplitudes of heaving acceleration are 0.5, 1.0, 1.5, and 2.0 m/s2, respectively. The effects of heaving acceleration amplitude on the flow instability are shown in Figures 6 and 7.Figure 6
Variation of critical mass flow rate with time for different amplitude of heaving acceleration.Figure 7
Variation of critical pressure drop with time for different amplitude of heaving acceleration.It can be seen from Figures6 and 7 that the fluctuation magnitude of CMFR and CPD increases with the increase in the amplitude of heaving acceleration. The data in Figures 6 and 7 are fitted, and the relative change values of CMFR and CPD can be written as follows:
(20)|m˙c,M-m˙c,0m˙c,0|=0.06679a0,(21)|Δpc,M-Δpc,0Δpc,0|=0.11998a0,
where m˙c,M and Δpc,M are the maximum values of CMFR and CPD in the heaving period, respectively.Combining (18) and (20), the variation of CMFR with time under the heaving movement can be obtained as follows:
(22)m˙c=m˙c,0[1+0.06679a0sin(2πTt)].Combining (19) and (21), the variation of CPD with time under the heaving movement can be obtained as follows:
(23)Δpc=Δpc,0[1+0.11998a0sin(2πTt)].


## 3.4. Effect of Heaving Movement on Flow Instability in U-Tubes with Different Length

Based on (16) and (17), the flow instability in U-tubes with different length under the heaving movement is investigated, and the results are shown in Figures 8–11.Figure 8
Variation of critical mass flow rate with time for different length.Figure 9
Variation of critical pressure drop with time for different length.Figure 10
Relative change of critical mass flow rate with U-tube length.Figure 11
Relative change of critical pressure drop with U-tube length.In Figures8 and 9, the heaving period is 10 s and the amplitude of heaving acceleration is 0.5 m/s2. l=0.7l0 and l=l0 are the shortest and longest length of U-tubes in the marine steam generator, respectively. It can be seen from Figure 8 that the CMFR of both the long tube (l=l0) and short tube (l=0.7l0) is lower than that of the U-tube with the middle length (l=0.8l0).It can be seen from Figure9 that the CPD of short tube (l=0.7l0) is higher than that of others, which means that the reverse flow will occur in the short tube easily. From Figures 8 and 9, it can be also found that the effects of heaving movement on the flow instability in U-tubes with different length are similar, which means that the heaving movement has little influence on the space distribution of reverse flow.Figures10 and 11 give the relative changes of the CMFR and CPD with U-tube length. Where the heaving period is 10 s and the amplitudes of heaving acceleration are 0.5, 1.0, 1.5, and 2.0 m/s2, respectively. It can be seen from Figures 10 and 11 that the amplitude of heaving acceleration has large influence on the relative changes of CMFR and CPD; namely, the CMFR and CPD increase with the increase in the heaving acceleration amplitude. However the effect of U-tube length on the CMFR is more obvious than that on the CPD.


## 4. Conclusions

Under heaving movement conditions, the critical pressure drop (CPD) and critical mass flow rate (CMFR) are given and analyzed, which relate to the occurrence of reverse flow in U-tubes. One-dimensional thermal-hydraulic model is established, and the effects of heaving movement on the flow instability in U-tubes are discussed. The detailed conclusions are drawn as follows.(1)
While the heaving acceleration decreases, the additional and total pressure drop will decrease, and the variation of total pressure drop and mass flow rate with time agrees with that of heaving acceleration. The heaving movement will affect the occurrence of reverse flow in the marine steam generator U-tubes.(2)
The fluctuation periods of CPD and CMFR are the same as the heaving period, but the heaving period has little influence on the CPD and CMFR fluctuation magnitude.(3)
The variations of CPD and CMFR with time under the heaving movement are derived, which are the linear function of the heaving acceleration amplitude.(4)
The U-tube length has little influence on the relative change of CPD and CMFR compared with the heaving acceleration amplitude, which means that the heaving movement has little influence on the space distribution of reverse flow in marine steam generator U-tubes.



---

*Source: 432817-2013-12-16.xml*