# Effects of Diameters on Countercurrent Flow Limitation at a Square Top End in Vertical Pipes

Authors: Michio Murase; Koji Nishida; Toshihide Torige; Toshiya Takaki; Raito Goda; Akio Tomiyama
Journal: Science and Technology of Nuclear Installations
 (2018)

---

## Abstract

The falling liquid flow rate under flooding conditions is limited at a square top end of a vertical pipe in the pressurizer surge line with the diameter of about 300 mm that consists of a vertical pipe, a vertical elbow, and a slightly inclined pipe with elbows. In this study, therefore, we evaluated effects of diameters on countercurrent flow limitation (CCFL) at the square top end in vertical pipes by using existing air-water data in the diameter range ofD = 19-250 mm. As a result, we found that there was a strong relationship between the constant CK and the slopem in the Wallis-type correlation where the Kutateladze parameters were used for the dimensionless gas and liquid velocities. The constant CK and the slopem increased when the water level is increased in the upper tankh. CCFL at the square top end of the vertical pipes could be expressed by the Kutateladze parameters with CK = 1.53±0.11 andm = 0.97 forD ≥ 30 mm. The CK values were smaller forD = 19-25 mm than those forD ≥ 30 mm.


---

## Body

## 1. Introduction

Under postulated accident conditions such as loss-of-coolant accidents (LOCAs) in pressurized water reactors (PWRs), steam and condensate water form countercurrent flows in a hot leg (consisting of a horizontal pipe, a 50-deg vertical elbow, and a short inclined pipe) and a pressurizer surge line (consisting of a slightly inclined pipe with elbows, a vertical elbow, and a vertical pipe), and flooding may occur. For transient and accident analyses of PWRs, characteristics of countercurrent flow limitation (CCFL) should be considered strictly, where CCFL is defined by the relationship between the time-averaged gas superficial velocity and falling liquid superficial velocity,JG and JL, under flooding conditions, because the falling liquid flow rate affects the water mass in the reactor core and cooling of the fuel rods.Previously, we (Murase et al. [1]) developed a one-dimensional computation method with parameters adjusted from CCFL data in hot leg and pressurizer surge line models (Mayinger et al. [2]; Minami et al. [3]; Futatsugi et al. [4]), and we could predict CCFL in nearly horizontal pipes for the hot leg and pressurizer surge line within a practical uncertainty. There are only a few CCFL experiments simulating the pressurizer surge line with a diameter of about 300 mm. In available data (Takeuchi et al. [5]; Futatsugi et al. [4]; Yu et al. [6, 7]), the pipe diameters wereD = 30-90 mm, so that the CCFL data were not sufficient to evaluate effects of the diameters. In the pressurizer surge line, the falling liquid flow rate is limited at the square top end of the vertical pipe, and so we (Yamamoto et al. [8]) derived a correlation for CCFL-U at the square top end for the pressurizer surge line by using existing data forD = 19-140 mm (Richter [9]; Doi et al. [10]). However, CCFL-U behavior is very complex and there are some disagreements among experimental data. Therefore, it is important to evaluate effects of the diameters on CCFL-U at the square top end of vertical pipes and to apply these results to the pressurizer surge line.In this study, we evaluated effects of diameters on CCFL-U at the square top end in vertical pipes by using existing air-water data reported by Richter [9] forD = 19-140 mm, Doi et al. [10] forD = 30-60 mm, Wallis and Kuo [11] forD = 19-145 mm, Bharathan et al. [12] forD = 19-250 mm, and Matsumura and Kaminaga [13] forD = 20 mm. From data by Richter [9], Wallis and Kuo [11], and Bharathan et al. [12], CCFL-U for large diameters was evaluated. The effects of the water level in the upper tank were evaluated from data mainly by Doi et al. [10], and the difference between data by Richter [9] and Doi et al. [10] was discussed. From the results, CCFL-U characteristics were classified depending on the diameter, the shape of the top end, and the water level in the upper tank, and features of CCFL-U were summarized. All data used in this study were obtained under air-water conditions at atmospheric pressure.


## 2. Previous Studies on CCFL at the Square Top End

### 2.1. General Form of CCFL Correlation

For safety analyses during transients and accidents in nuclear power plants, the CCFL correlation by Wallis [15] has often been applied to evaluate the falling liquid flow rate. The general form of the Wallis correlation is given by(1)HG∗1/2+mHL∗1/2=Ci,i=KorW(2)Hk∗=JkρkgwρL-ρG1/2,k=GorLwhere g [m/s2] is the gravitational acceleration,J [m/s] is the superficial velocity, H∗ [-] is the dimensionless velocity, w [m] is the characteristic length, andρ [kg/m3] is the density. The slopem and the constantCi are determined from experiments, where the subscripts K and W show the Kutateladze parameter and Wallis parameter for Hk∗, respectively. Bankoff et al. [16] defined the characteristic length w by(3)w=D1-βLβ,0≤β≤1,(4)L=σgρL-ρG1/2whereD [m] is the diameter,L [m] is the Laplace capillary length, andσ [N/m] is the surface tension. In (1), H∗ indicates the Wallis parameter J∗ at β = 0, while H∗ is the Kutateladze parameter K∗ at β = 1. J∗ and K∗ can be converted to each other by using D∗, as follows:(5)J∗=K∗D∗1/2,D∗=DL.Equations (2) and (3) show primary parameters for the length scale (D orL) and fluid properties (ρ andσ). However, it is well known that w in (3) changes depending onD in vertical pipes (Wallis and Makkenchery [17]) and that the liquid viscosity (which is not included in (2) and (4)) affects CCFL in vertical pipes (Wallis [15]).


### 2.2. CCFL Characteristics and Technical Issues

Figure1 shows CCFL characteristics at the square top end in vertical pipes. CCFL data reported by Richter [9] and Doi et al. [10] are well expressed by the Kutateladze parameters. From the data ofD =19-140 mm in Figure 1(a), the Wallis-type correlation (6) was derived using the least-square method (Yamamoto et al. [8]).(6)KG∗1/2+0.9KL∗1/2=1.5±0.1Figure 1(b) shows the CCFL constant CK, which was obtained for each experimental case using the least-square method. Data by Wallis and Kuo [11] are not the CCFL constant but the zero water penetration (ZWP), which generally agrees with CK. Technical issues are disagreements between data by Richter [9] and Doi et al. [10] forD = about 30 mm, and data by Richter [9] and Wallis and Kuo [11] forD = about 140 mm. The CK values by Richter [9] and the ZWP values by Wallis and Kuo [11] suggest that CCFL characteristics might be expressed by the Wallis parameter (CW = 0.7) for the region of small diameters ofD < 40 mm, but this should be confirmed from CCFL data.Figure 1
CCFL characteristics at square top ends: (a) CCFL expressed by Kutateladze parameters and (b) CCFL constant (broken-line box: disagreement between data).
(a) (b)


### 2.3. Experimental Conditions in Previous Studies

Figure2 shows the shapes near the top ends of the vertical pipes used in previous experiments, and Table 1 lists major experimental conditions. Wallis and Kuo [11] used a protruding top end (Figure 2(a)), and the water level in the upper tank was low (h < 7 mm), where the gas core continued from the lower tank to the atmosphere. Bharathan et al. [12] used a square top end (Figure 2(b)) and also a protruding top end, and the water level in the upper tank was midlevel, where the upper part of the top end was covered by the air-water mixture. Richter [9] referred to CCFL data by Wallis and Makkenchery [17], but we could not find the CCFL data in Wallis and Makkenchery [17]. Therefore, we do not know the shape near the top end of the vertical pipe reported by Richter [9]. From the comparison of data by Richter [9] with data by Wallis and Kuo [11] and Bharathan et al. [12], we assumed that a protruding top end was used and the water level was low. Matsumura and Kaminaga [13] used a square top end with a horizontal plate (Figure 2(c)) to ensure clear observations. Doi et al. [10] used a square top end with three upper tanks (Figure 2(d)): a rectangular tank (UT1), a circular tank (UT2), and a cylindrical tank (UT3). All experiments [9–13] listed in Table 1 were done under air-water conditions at the room temperature and atmospheric pressure.Table 1
Experimental conditions (D, diameter; h, water level; ∗our judgement of the condition).
Reference Top end D [mm] h [mm] Wallis and Kuo [11] Protruding 6.4, 19, 31.2, 38.2, < 7 63.5, 113, 145 Bharathan et al. [12] Protruding 19, 148, 250 Mid-level Square 19, 25, 51, 148 Richter [9] (Protruding)∗ 19, 25, 32, 38, 44, 51, (< 7)∗ 64, 76, 89, 100, 140 Matsumura and Square 12, 20 50, 100, 200 Kaminaga [13] Doi et al. [10] Square 30, 45, 60 (UT1) 200 30, 45, 60 (UT2) 200 30 (UT3) 100, 200, 300, 450, 600Figure 2
Shapes near the top end of the vertical pipe used by (a) Wallis and Kuo [11], (b) Bharathan et al. [12], (c) Matsumura and Kaminaga [13], and (d) Doi et al. [10].
(a) (b) (c) (d)


## 2.1. General Form of CCFL Correlation

For safety analyses during transients and accidents in nuclear power plants, the CCFL correlation by Wallis [15] has often been applied to evaluate the falling liquid flow rate. The general form of the Wallis correlation is given by(1)HG∗1/2+mHL∗1/2=Ci,i=KorW(2)Hk∗=JkρkgwρL-ρG1/2,k=GorLwhere g [m/s2] is the gravitational acceleration,J [m/s] is the superficial velocity, H∗ [-] is the dimensionless velocity, w [m] is the characteristic length, andρ [kg/m3] is the density. The slopem and the constantCi are determined from experiments, where the subscripts K and W show the Kutateladze parameter and Wallis parameter for Hk∗, respectively. Bankoff et al. [16] defined the characteristic length w by(3)w=D1-βLβ,0≤β≤1,(4)L=σgρL-ρG1/2whereD [m] is the diameter,L [m] is the Laplace capillary length, andσ [N/m] is the surface tension. In (1), H∗ indicates the Wallis parameter J∗ at β = 0, while H∗ is the Kutateladze parameter K∗ at β = 1. J∗ and K∗ can be converted to each other by using D∗, as follows:(5)J∗=K∗D∗1/2,D∗=DL.Equations (2) and (3) show primary parameters for the length scale (D orL) and fluid properties (ρ andσ). However, it is well known that w in (3) changes depending onD in vertical pipes (Wallis and Makkenchery [17]) and that the liquid viscosity (which is not included in (2) and (4)) affects CCFL in vertical pipes (Wallis [15]).


## 2.2. CCFL Characteristics and Technical Issues

Figure1 shows CCFL characteristics at the square top end in vertical pipes. CCFL data reported by Richter [9] and Doi et al. [10] are well expressed by the Kutateladze parameters. From the data ofD =19-140 mm in Figure 1(a), the Wallis-type correlation (6) was derived using the least-square method (Yamamoto et al. [8]).(6)KG∗1/2+0.9KL∗1/2=1.5±0.1Figure 1(b) shows the CCFL constant CK, which was obtained for each experimental case using the least-square method. Data by Wallis and Kuo [11] are not the CCFL constant but the zero water penetration (ZWP), which generally agrees with CK. Technical issues are disagreements between data by Richter [9] and Doi et al. [10] forD = about 30 mm, and data by Richter [9] and Wallis and Kuo [11] forD = about 140 mm. The CK values by Richter [9] and the ZWP values by Wallis and Kuo [11] suggest that CCFL characteristics might be expressed by the Wallis parameter (CW = 0.7) for the region of small diameters ofD < 40 mm, but this should be confirmed from CCFL data.Figure 1
CCFL characteristics at square top ends: (a) CCFL expressed by Kutateladze parameters and (b) CCFL constant (broken-line box: disagreement between data).
(a) (b)


## 2.3. Experimental Conditions in Previous Studies

Figure2 shows the shapes near the top ends of the vertical pipes used in previous experiments, and Table 1 lists major experimental conditions. Wallis and Kuo [11] used a protruding top end (Figure 2(a)), and the water level in the upper tank was low (h < 7 mm), where the gas core continued from the lower tank to the atmosphere. Bharathan et al. [12] used a square top end (Figure 2(b)) and also a protruding top end, and the water level in the upper tank was midlevel, where the upper part of the top end was covered by the air-water mixture. Richter [9] referred to CCFL data by Wallis and Makkenchery [17], but we could not find the CCFL data in Wallis and Makkenchery [17]. Therefore, we do not know the shape near the top end of the vertical pipe reported by Richter [9]. From the comparison of data by Richter [9] with data by Wallis and Kuo [11] and Bharathan et al. [12], we assumed that a protruding top end was used and the water level was low. Matsumura and Kaminaga [13] used a square top end with a horizontal plate (Figure 2(c)) to ensure clear observations. Doi et al. [10] used a square top end with three upper tanks (Figure 2(d)): a rectangular tank (UT1), a circular tank (UT2), and a cylindrical tank (UT3). All experiments [9–13] listed in Table 1 were done under air-water conditions at the room temperature and atmospheric pressure.Table 1
Experimental conditions (D, diameter; h, water level; ∗our judgement of the condition).
Reference Top end D [mm] h [mm] Wallis and Kuo [11] Protruding 6.4, 19, 31.2, 38.2, < 7 63.5, 113, 145 Bharathan et al. [12] Protruding 19, 148, 250 Mid-level Square 19, 25, 51, 148 Richter [9] (Protruding)∗ 19, 25, 32, 38, 44, 51, (< 7)∗ 64, 76, 89, 100, 140 Matsumura and Square 12, 20 50, 100, 200 Kaminaga [13] Doi et al. [10] Square 30, 45, 60 (UT1) 200 30, 45, 60 (UT2) 200 30 (UT3) 100, 200, 300, 450, 600Figure 2
Shapes near the top end of the vertical pipe used by (a) Wallis and Kuo [11], (b) Bharathan et al. [12], (c) Matsumura and Kaminaga [13], and (d) Doi et al. [10].
(a) (b) (c) (d)


## 3. CCFL Characteristics at the Square Top End

### 3.1. CCFL in Vertical Pipes of Large Diameters

Our main interest is CCFL in large diameters for application to pressurizer surge lines withD = about 300 mm. Bharathan et al. [12] reported some air-water data with such large diameters. Figure 3 compares CCFL data in large diameters with (6). The single data point by Wallis and Kuo [11] is for the zero water penetration with the protruding top end. The data by Richter [9] and Wallis and Kuo [11] agreed well with data by Bharathan et al. [12] with the protruding top end. Hence the data by Richter [9] might be measured with the protruding top end. This was the reason why we judged the top end in the experiment by Richter [9] to be “protruding” in Table 1. The difference between data points by Richter [9] and Wallis and Kuo [11] forD = 140-145 mm in Figure 1(b) was due to the small slope in the middle range of KG∗1/2 and the large slope in the large range of KG∗1/2.Figure 3
CCFL characteristics for large diameters and effects of protruding end (A: smallCK value by Richter [9] in Figure 1(b)).Our interest is not CCFL with the protruding top end but CCFL with the square top end.KL∗1/2 values with the square top end are larger in the wide range of KG∗1/2 than those with the protruding top end. Figure 3 shows that (6) can be applied to large diameters up toD = 250 mm (Richter et al. [14]).


### 3.2. Effects of Water Level in Upper Tank on CCFL

We (Doi et al. [10]) measured CCFL characteristics at the square top end using the three upper tanks shown in Figure 2(d). Figure 4 shows the effects of the upper tanks and water levels there on the CCFL constant CK and the slopem. As shown in Figure 4(a), the shape of the upper tanks affected CK andm due to different sloshing behavior there and pressure fluctuation in the lower tank, which affected fluctuation of the air flow rate into the vertical pipe from the lower tank. Clear effects of the diameters on CK andm were not observed, and effects of the upper tanks on CK andm were not significant.Figure 4
Effects of (a) upper tank by Doi et al. [10] (UT1-UT3: cf. Figure 2(d)) and (b) water levelh on CK andm.
(a) (b)As shown in Figure4(b), CK andm withD = 30 mm increased with increasing the water levelh in the upper tank (Doi et al. [10]). As decreasingh, CK andm by Doi et al. [10] approached those withD = 32 mm by Wallis and Kuo [11] with a low water level and Richter [9]. This was the reason we judged that the water level in the experiment by Richter [9] was low in Table 1. Matsumura and Kaminaga [13] reported that clear differences were not measured forD = 20 mm andh = 0.05-0.2 m. Therefore, the effects of the water levels may depend on the shape of the upper tank and pipe diameter. Figure 4 suggests a clear interrelation between CK andm, which will be discussed later.


### 3.3. CCFL in Vertical Pipes of Small Diameters

As shown in Figure1(b), there was clear disagreement between the data by Richter [9] and Doi et al. [10] forD = about 30 mm, and CK values by Richter [9] and ZWP by Wallis and Kuo [11] suggested that CCFL characteristics might be expressed by the Wallis parameter (CW = 0.7) in the region ofD < 40 mm. The large difference of CK values by Doi et al. [10] forD = about 30 mm in Figure 1(b) was due to effects of the water level in the upper tank as shown in Figure 4(b).Figure5 shows CCFL characteristics for (a) medium and (b) small diameters. For the medium diameters ofD = 44-51 mm, all data points were within the uncertainty of (6) as shown in Figure 5(a). ForD = 51 mm, the protruding or square end and the water level at low or midlevel (cf. Table 1) did not affect CCFL. For the small diameters ofD = 19-25 mm, CCFL characteristics were clearly different from (6) as shown in Figure 5(b). The protruding or square end did not affect CCFL from the 19 mm diameter data by Bharathan et al. [12]. The slope of data with the low water level by Richter [9] (D = 19 mm) was clearly smaller than that of other data with the water level at midlevel. This suggests that the small CK values for the small diameters ofD = 19-38 mm by Richter [9] and Wallis and Kuo [11] in Figure 1(b) might be due to the low water level. For data with the square top end, the slope was similar to that of (6) and there were no clear trends for diameters.Figure 5
CCFL characteristics in (a) medium diameters and (b) small diameters with water levels at midlevel (exceptD = 19 mm by Richter [9]).
(a) (b)Figure6 shows CCFL characteristics in small diameters with low water levels by Richter [9]. As shown in Figure 6(a), the slopem could be expressed by an exponent function of the dimensionless diameter D∗ asm = 0.16D∗0.58. Figure 6(b) shows CCFL data plotted for JG∗1/2 and D∗0.58JL∗1/2. From the data points in Figure 6(b), the CCFL correlation was derived using the least-square method as(7)KG∗1/2+0.16D∗0.58KL∗1/2=0.71±0.05All data points in the diameter range ofD = 19-44 mm were within the small uncertainty of ±0.05. The data points for the diameter ofD = 44 mm were also within the uncertainty of (6) as shown in Figure 5(a).Figure 6
CCFL characteristics for small diameters with low water levels by Richter [9]: (a) slopem, and (b) CCFL.
(a) (b)


### 3.4. Features of CCFL at the Square Top End

Figure7 shows the constant CK and the slopem in (1) for CCFL at the square top end. The CK andm values differed largely depending on the shape of the top end (square or protruding), the diameter, and the water level in the upper tank. The difference was mainly due to the shape of the top end (square or protruding) in the large diameters of D∗>37 (D > 100 mm) and the water level in the upper tank in the small diameters of D∗<16 (D < 44 mm). This means that the shape of the top end and water level did not significantly affect CCFL in the diameter range ofD = 44-100 mm (D∗ = 16-37).Figure 7
CCFL characteristics at the square top end: (a) CCFL constantCK and (b) slopem.
(a) (b)CCFL characteristics depended on the shape of the top end (square or protruding), diameters, and the water level in the upper tank, and Table2 lists classifications of CCFL. Figure 8 shows features of CCFL characteristics at the square top end (including the protruding top end). As shown in Figure 8(a), a clear relationship between CK andm was found except for the region “R-IV” (small diameters ofD = 19-25 mm with the water level at midlevel). The relationship between CK andm was expressed by(8)CK=0.52ln⁡m+1.54±0.051.19≤CK≤1.7,0.5≤m≤1.4.In region “R-I”, the CCFL characteristics (CK = 1.5±0.1 andm = 0.88±0.13) were measured withD = 30-100 mm andh ≤ 200 mm or midlevel (cf. “Basic” in Table 2). On increasing the diameter or the water level with the square top end, CK andm increased (region “R-II”). In region “R-III”, all data were obtained with the protruding top end. CK andm in region “R-IV” with the small diameters ofD = 19-25 mm and the water level at midlevel were outside the range of values calculated by using (8). Our major interest is CCFL in large diameters with the square top end, and therefore all data in Regions “R-I” and “R-II” and data with large diameters in region “R-III” are plotted in Figure 8(b). From data points in Figure 8(b), we derived the CCFL correlation using the least-square method as follows: (9)KG∗1/2+0.97KL∗1/2=1.53±0.1130mm≤D≤250mm,11≤D∗≤94.95 % of the 290 data points were within the uncertainty of ±0.11. Some data points at large JG∗1/2 withD = 148 mm and the square top end and some data points at medium JG∗1/2 with the protruding top end, which were larger and smaller than (9), respectively, were outside the uncertainty of ±0.11. Though the difference between (6) and (9) was small, we recommend using (9) because of the wide range of CCFL data used for its derivation.Table 2
Classifications of CCFL at the square top end (Region “R-”: cf. Figure8(a)).
Classification Top end D [mm] h [mm] Data base Basic, pe (R-I) Protruding 44, 51, 64, 76, 89, 100 Low Richter [9]  Basic, se (R-I) Square 30, 45, 60 100, 200 Doi et al. [10] 51 Mid-level Bharathan et al. [12]  Highh, se (R-II) Square 30 300, 450, 600 Doi et al. [10]  LargeD, se (R-II) Square 148 Mid-level Bharathan et al. [12]  LargeD, pe (R-III) Protruding 140 Low Richter [9] 148 Mid-level Bharathan et al. [12] 250 Mid-level Richter et al. [14]  Lowh, pe (R-III) Protruding 19, 25, 32, 38 Low Richter [9]  SmallD, pe (R-IV) Protruding 19 Mid-level Bharathan et al. [12]  SmallD, se (R-IV) Square 19, 25 Mid-level Bharathan et al. [12] 20 50, 100, 200 Matumura and Kaminaga [13]Figure 8
Features of CCFL at the square top end: (a) relationship betweenCK andm (classifications: cf. Table 2) and (b) CCFL for medium and large diameters.
(a) (b)Table3 summarizes effects of parameters on CCFL at the square top end. Effects of “Large”, “Medium”, and “Small” were determined by comparison with conditions of the square top end and water level at midlevel forD = 30-100 mm where the uncertainty was ±0.1 (6.7%) for CK = 1.5. Our target to improve a CCFL correlation was uncertainty smaller than ±10% for the CCFL constant CK, and “Medium” was within the target of ±10%. “Large” was uncertainty larger than ±10% and we recommend using a different CCFL correlation for “Large”. However, the data with the protruding end forD = 140-250 mm were included to derive (9) due to limited data for large diameters.Table 3
Effects of parameters on CCFL at square top end comparing with conditions of the square top end and water level at midlevel forD = 30-100 mm (“Medium” is uncertainty range of about 5-10 % for CK, for “Large” different correlation is recommended, bold area; database used for (9)).
D [mm] 19-25  30-38 44-100 140-250 Diameter(Square end) Large  Small (base case) Medium  Protruding end  Small  Medium Large  Low water level  Large  SmallFor low water levels, we recommend using the minimum of (7) for small diameters and (9) for medium to large diameters. For the small diameters of 19-25 mm with a square end, we recommend using a correlation specific for a diameter, because CCFL characteristics were different among the databases as shown in Figure 5(b).


## 3.1. CCFL in Vertical Pipes of Large Diameters

Our main interest is CCFL in large diameters for application to pressurizer surge lines withD = about 300 mm. Bharathan et al. [12] reported some air-water data with such large diameters. Figure 3 compares CCFL data in large diameters with (6). The single data point by Wallis and Kuo [11] is for the zero water penetration with the protruding top end. The data by Richter [9] and Wallis and Kuo [11] agreed well with data by Bharathan et al. [12] with the protruding top end. Hence the data by Richter [9] might be measured with the protruding top end. This was the reason why we judged the top end in the experiment by Richter [9] to be “protruding” in Table 1. The difference between data points by Richter [9] and Wallis and Kuo [11] forD = 140-145 mm in Figure 1(b) was due to the small slope in the middle range of KG∗1/2 and the large slope in the large range of KG∗1/2.Figure 3
CCFL characteristics for large diameters and effects of protruding end (A: smallCK value by Richter [9] in Figure 1(b)).Our interest is not CCFL with the protruding top end but CCFL with the square top end.KL∗1/2 values with the square top end are larger in the wide range of KG∗1/2 than those with the protruding top end. Figure 3 shows that (6) can be applied to large diameters up toD = 250 mm (Richter et al. [14]).


## 3.2. Effects of Water Level in Upper Tank on CCFL

We (Doi et al. [10]) measured CCFL characteristics at the square top end using the three upper tanks shown in Figure 2(d). Figure 4 shows the effects of the upper tanks and water levels there on the CCFL constant CK and the slopem. As shown in Figure 4(a), the shape of the upper tanks affected CK andm due to different sloshing behavior there and pressure fluctuation in the lower tank, which affected fluctuation of the air flow rate into the vertical pipe from the lower tank. Clear effects of the diameters on CK andm were not observed, and effects of the upper tanks on CK andm were not significant.Figure 4
Effects of (a) upper tank by Doi et al. [10] (UT1-UT3: cf. Figure 2(d)) and (b) water levelh on CK andm.
(a) (b)As shown in Figure4(b), CK andm withD = 30 mm increased with increasing the water levelh in the upper tank (Doi et al. [10]). As decreasingh, CK andm by Doi et al. [10] approached those withD = 32 mm by Wallis and Kuo [11] with a low water level and Richter [9]. This was the reason we judged that the water level in the experiment by Richter [9] was low in Table 1. Matsumura and Kaminaga [13] reported that clear differences were not measured forD = 20 mm andh = 0.05-0.2 m. Therefore, the effects of the water levels may depend on the shape of the upper tank and pipe diameter. Figure 4 suggests a clear interrelation between CK andm, which will be discussed later.


## 3.3. CCFL in Vertical Pipes of Small Diameters

As shown in Figure1(b), there was clear disagreement between the data by Richter [9] and Doi et al. [10] forD = about 30 mm, and CK values by Richter [9] and ZWP by Wallis and Kuo [11] suggested that CCFL characteristics might be expressed by the Wallis parameter (CW = 0.7) in the region ofD < 40 mm. The large difference of CK values by Doi et al. [10] forD = about 30 mm in Figure 1(b) was due to effects of the water level in the upper tank as shown in Figure 4(b).Figure5 shows CCFL characteristics for (a) medium and (b) small diameters. For the medium diameters ofD = 44-51 mm, all data points were within the uncertainty of (6) as shown in Figure 5(a). ForD = 51 mm, the protruding or square end and the water level at low or midlevel (cf. Table 1) did not affect CCFL. For the small diameters ofD = 19-25 mm, CCFL characteristics were clearly different from (6) as shown in Figure 5(b). The protruding or square end did not affect CCFL from the 19 mm diameter data by Bharathan et al. [12]. The slope of data with the low water level by Richter [9] (D = 19 mm) was clearly smaller than that of other data with the water level at midlevel. This suggests that the small CK values for the small diameters ofD = 19-38 mm by Richter [9] and Wallis and Kuo [11] in Figure 1(b) might be due to the low water level. For data with the square top end, the slope was similar to that of (6) and there were no clear trends for diameters.Figure 5
CCFL characteristics in (a) medium diameters and (b) small diameters with water levels at midlevel (exceptD = 19 mm by Richter [9]).
(a) (b)Figure6 shows CCFL characteristics in small diameters with low water levels by Richter [9]. As shown in Figure 6(a), the slopem could be expressed by an exponent function of the dimensionless diameter D∗ asm = 0.16D∗0.58. Figure 6(b) shows CCFL data plotted for JG∗1/2 and D∗0.58JL∗1/2. From the data points in Figure 6(b), the CCFL correlation was derived using the least-square method as(7)KG∗1/2+0.16D∗0.58KL∗1/2=0.71±0.05All data points in the diameter range ofD = 19-44 mm were within the small uncertainty of ±0.05. The data points for the diameter ofD = 44 mm were also within the uncertainty of (6) as shown in Figure 5(a).Figure 6
CCFL characteristics for small diameters with low water levels by Richter [9]: (a) slopem, and (b) CCFL.
(a) (b)


## 3.4. Features of CCFL at the Square Top End

Figure7 shows the constant CK and the slopem in (1) for CCFL at the square top end. The CK andm values differed largely depending on the shape of the top end (square or protruding), the diameter, and the water level in the upper tank. The difference was mainly due to the shape of the top end (square or protruding) in the large diameters of D∗>37 (D > 100 mm) and the water level in the upper tank in the small diameters of D∗<16 (D < 44 mm). This means that the shape of the top end and water level did not significantly affect CCFL in the diameter range ofD = 44-100 mm (D∗ = 16-37).Figure 7
CCFL characteristics at the square top end: (a) CCFL constantCK and (b) slopem.
(a) (b)CCFL characteristics depended on the shape of the top end (square or protruding), diameters, and the water level in the upper tank, and Table2 lists classifications of CCFL. Figure 8 shows features of CCFL characteristics at the square top end (including the protruding top end). As shown in Figure 8(a), a clear relationship between CK andm was found except for the region “R-IV” (small diameters ofD = 19-25 mm with the water level at midlevel). The relationship between CK andm was expressed by(8)CK=0.52ln⁡m+1.54±0.051.19≤CK≤1.7,0.5≤m≤1.4.In region “R-I”, the CCFL characteristics (CK = 1.5±0.1 andm = 0.88±0.13) were measured withD = 30-100 mm andh ≤ 200 mm or midlevel (cf. “Basic” in Table 2). On increasing the diameter or the water level with the square top end, CK andm increased (region “R-II”). In region “R-III”, all data were obtained with the protruding top end. CK andm in region “R-IV” with the small diameters ofD = 19-25 mm and the water level at midlevel were outside the range of values calculated by using (8). Our major interest is CCFL in large diameters with the square top end, and therefore all data in Regions “R-I” and “R-II” and data with large diameters in region “R-III” are plotted in Figure 8(b). From data points in Figure 8(b), we derived the CCFL correlation using the least-square method as follows: (9)KG∗1/2+0.97KL∗1/2=1.53±0.1130mm≤D≤250mm,11≤D∗≤94.95 % of the 290 data points were within the uncertainty of ±0.11. Some data points at large JG∗1/2 withD = 148 mm and the square top end and some data points at medium JG∗1/2 with the protruding top end, which were larger and smaller than (9), respectively, were outside the uncertainty of ±0.11. Though the difference between (6) and (9) was small, we recommend using (9) because of the wide range of CCFL data used for its derivation.Table 2
Classifications of CCFL at the square top end (Region “R-”: cf. Figure8(a)).
Classification Top end D [mm] h [mm] Data base Basic, pe (R-I) Protruding 44, 51, 64, 76, 89, 100 Low Richter [9]  Basic, se (R-I) Square 30, 45, 60 100, 200 Doi et al. [10] 51 Mid-level Bharathan et al. [12]  Highh, se (R-II) Square 30 300, 450, 600 Doi et al. [10]  LargeD, se (R-II) Square 148 Mid-level Bharathan et al. [12]  LargeD, pe (R-III) Protruding 140 Low Richter [9] 148 Mid-level Bharathan et al. [12] 250 Mid-level Richter et al. [14]  Lowh, pe (R-III) Protruding 19, 25, 32, 38 Low Richter [9]  SmallD, pe (R-IV) Protruding 19 Mid-level Bharathan et al. [12]  SmallD, se (R-IV) Square 19, 25 Mid-level Bharathan et al. [12] 20 50, 100, 200 Matumura and Kaminaga [13]Figure 8
Features of CCFL at the square top end: (a) relationship betweenCK andm (classifications: cf. Table 2) and (b) CCFL for medium and large diameters.
(a) (b)Table3 summarizes effects of parameters on CCFL at the square top end. Effects of “Large”, “Medium”, and “Small” were determined by comparison with conditions of the square top end and water level at midlevel forD = 30-100 mm where the uncertainty was ±0.1 (6.7%) for CK = 1.5. Our target to improve a CCFL correlation was uncertainty smaller than ±10% for the CCFL constant CK, and “Medium” was within the target of ±10%. “Large” was uncertainty larger than ±10% and we recommend using a different CCFL correlation for “Large”. However, the data with the protruding end forD = 140-250 mm were included to derive (9) due to limited data for large diameters.Table 3
Effects of parameters on CCFL at square top end comparing with conditions of the square top end and water level at midlevel forD = 30-100 mm (“Medium” is uncertainty range of about 5-10 % for CK, for “Large” different correlation is recommended, bold area; database used for (9)).
D [mm] 19-25  30-38 44-100 140-250 Diameter(Square end) Large  Small (base case) Medium  Protruding end  Small  Medium Large  Low water level  Large  SmallFor low water levels, we recommend using the minimum of (7) for small diameters and (9) for medium to large diameters. For the small diameters of 19-25 mm with a square end, we recommend using a correlation specific for a diameter, because CCFL characteristics were different among the databases as shown in Figure 5(b).


## 4. Discussion

Our final target is to evaluate CCFL characteristics in an actual pressurizer surge line. Hence we have to evaluate effects of pressure and temperature on CCFL and also system characteristics of the pressurizer surge line. However, there are only a few CCFL experiments simulating the pressurizer surge line with the diameter of about 300 mm. In available data (Takeuchi et al. [5]; Futatsugi et al. [4]; Yu et al. [6, 7]), the pipe diameters wereD = 30-90 mm, and CCFL data were not sufficient to evaluate effects of the diameters and fluid properties. Therefore, the features of CCFL at the square top end discussed in Section 3.4 are important when we evaluate CCFL characteristics in the pressurizer surge line.


### 4.1. Effects of Pressure and Temperature

Ilyukhin et al. [18, 19] reported CCFL data withD = 20 mm at pressures ofP = 0.6-4.1 MPa and withD = 30 and 40 mm atP = 0.3 -1.6 MPa, respectively. The data points forD = 30 and 40 mm atP = 0.3 -1.6 MPa obtained by Ilyukhin et al. [19] were plotted on the diagram of Kk∗1/2D∗1/8(ρG/ρL)0.05. Since the data points for each pressure were not separated, we cannot evaluate Kk∗1/2 values from the data. They measured CCFL data with the protruding top and bottom ends at a high water level in the upper tank. Figure 8(a) and Table 3, therefore, may be useful for evaluating effects of fluid properties from the data by Ilyukhin et al. [19] withD = 20 mm atP = 0.6-4.1 MPa.


### 4.2. Computations with the Annular Flow Model

CCFL characteristics in vertical pipes can be computed by using the annular flow model with the maximum flow condition of∂JL∗/∂α=0 or ∂JG∗/∂α=0 and correlations for the interfacial friction factor and the wall friction factor (Richter [9]; Bharathan and Wallis [20]; Sudo [21]; Yamamoto et al. [8]). In the computations, the correlation for the interfacial friction factor is important. Sudo [21] improved the correlation for the interfacial friction factor by Bharathan and Wallis [20] to obtain good agreement between computed results and data by Richter [9]. As shown in Figure 7(a), however, the CK value by Richter [9] had a maximum at an intermediate diameter and decreased with increasing the diameter toward large diameters due to the protruding end. The correlation for the interfacial friction factor improved by Sudo [21] gave a small CK value for a large diameter (Yamamoto et al. [8]). Therefore, the correlation for the interfacial friction factor should be improved to apply it to a large diameter without a protruding end.


### 4.3. System Characteristics

Figure9 shows CCFL system characteristics of the pressurizer surge line (PSL) models. The dimensionless liquid velocity KL∗ of APEX (Takeuchi et al. [5]) (withD = 89 mm and steam-water) and the 1/10-scale model by Futatsugi et al. [4] (withD = 30 mm and air-water), which simulated the pressurizer surge line models in the AP600 and a PWR, respectively, is smaller than that of (9) for the square top end in vertical pipes. From the integration test data in Figure 9, the following CCFL correlation was derived using the least-square method (Yamamoto et al. [8]):(10)KG∗1/2+0.94KL∗1/2=1.24±0.1.Though the reason is not clear as to why the falling water velocity of the integration test data is smaller than that of (9), both sets of integration test data (Takeuchi et al. [5]; Futatsugi et al. [4]) can be well expressed by (10). Data used for deriving (10) are limited. Yu et al. [6, 7] measured CCFL using the 1/4-scale model (D = 90 mm) of the pressurizer surge line model of the AP1000 with air-water and steam-water at atmospheric pressure. They changed the water level in the pressurizer simulator up toh = 900 mm, and it affected CCFL. The effects ofh on CK andm in their experiments [6, 7] showed trends similar to those by Doi et al. [10] shown in Figure 4(b). Therefore, Figure 4(b), Figure 8(a), and Table 3 may be useful, when we evaluate data by Yu et al. [6, 7].Figure 9
System characteristics of CCFL in a pressurizer surge line.


### 4.4. Technical Issues

Gas-liquid countercurrent flows in vertical pipes have been studied for over a half century, and it is well known that CCFL characteristics strongly depend on the shapes at the top and bottom ends (Wallis [15]; Bankoff and Lee [22]). We (Yamamoto et al. [23]) classified CCFL in vertical pipes depending on the limiting locations at the square bottom end (CCFL-L) and the square top end (CCFL-U) and inside the vertical pipe (CCFL-P), and we found that the characteristic lengths defined by (3) areD,L, andD0.5L0.5 for CCFL-L, CCFL-U, and CCFL-P, respectively. However, (1)-(4) only show effects of primary parameters, and actual CCFL characteristics are very complex and depend on some secondary parameters such as the diameters (cf. Table 3), water levels, and fluid properties. Moreover, experimental data of CCFL differ among experimental facilities and conditions, and we have to consider some uncertainty among experimental data. Experiments at high pressure and temperature are few in number and effects of fluid properties as secondary parameters are still not clear. Wallis [15] showed that the liquid viscosity (which is not included in (1)-(4)) affects CCFL and he used the dimensionless viscosity number. We (Murase et al. [24]) used the viscosity ratio of the gas and liquid phases for CCFL-L, CCFL-U, and CCFL-P, but Ilyukhin et al. [18, 19] used the density ratio of the gas and liquid phases for CCFL with the protruding top and bottom ends. Therefore, CCFL correlations, which include effects of fluid properties, should be carefully used considering the application range. As the uncertainty of CCFL will increase on increasing the number of experimental data due to the inherent differences among the experimental facilities, the CCFL correlation should be derived using experimental data suitable for the system of concern.Moreover, CCFL characteristics from integration tests expressed by (10) are different from those of (9) as shown in Figure 9. Therefore, engineering judgement is required to choose an appropriate CCFL correlation. At present, we recommend using the correlation equation (10), for the actual pressurizer surge line.


## 4.1. Effects of Pressure and Temperature

Ilyukhin et al. [18, 19] reported CCFL data withD = 20 mm at pressures ofP = 0.6-4.1 MPa and withD = 30 and 40 mm atP = 0.3 -1.6 MPa, respectively. The data points forD = 30 and 40 mm atP = 0.3 -1.6 MPa obtained by Ilyukhin et al. [19] were plotted on the diagram of Kk∗1/2D∗1/8(ρG/ρL)0.05. Since the data points for each pressure were not separated, we cannot evaluate Kk∗1/2 values from the data. They measured CCFL data with the protruding top and bottom ends at a high water level in the upper tank. Figure 8(a) and Table 3, therefore, may be useful for evaluating effects of fluid properties from the data by Ilyukhin et al. [19] withD = 20 mm atP = 0.6-4.1 MPa.


## 4.2. Computations with the Annular Flow Model

CCFL characteristics in vertical pipes can be computed by using the annular flow model with the maximum flow condition of∂JL∗/∂α=0 or ∂JG∗/∂α=0 and correlations for the interfacial friction factor and the wall friction factor (Richter [9]; Bharathan and Wallis [20]; Sudo [21]; Yamamoto et al. [8]). In the computations, the correlation for the interfacial friction factor is important. Sudo [21] improved the correlation for the interfacial friction factor by Bharathan and Wallis [20] to obtain good agreement between computed results and data by Richter [9]. As shown in Figure 7(a), however, the CK value by Richter [9] had a maximum at an intermediate diameter and decreased with increasing the diameter toward large diameters due to the protruding end. The correlation for the interfacial friction factor improved by Sudo [21] gave a small CK value for a large diameter (Yamamoto et al. [8]). Therefore, the correlation for the interfacial friction factor should be improved to apply it to a large diameter without a protruding end.


## 4.3. System Characteristics

Figure9 shows CCFL system characteristics of the pressurizer surge line (PSL) models. The dimensionless liquid velocity KL∗ of APEX (Takeuchi et al. [5]) (withD = 89 mm and steam-water) and the 1/10-scale model by Futatsugi et al. [4] (withD = 30 mm and air-water), which simulated the pressurizer surge line models in the AP600 and a PWR, respectively, is smaller than that of (9) for the square top end in vertical pipes. From the integration test data in Figure 9, the following CCFL correlation was derived using the least-square method (Yamamoto et al. [8]):(10)KG∗1/2+0.94KL∗1/2=1.24±0.1.Though the reason is not clear as to why the falling water velocity of the integration test data is smaller than that of (9), both sets of integration test data (Takeuchi et al. [5]; Futatsugi et al. [4]) can be well expressed by (10). Data used for deriving (10) are limited. Yu et al. [6, 7] measured CCFL using the 1/4-scale model (D = 90 mm) of the pressurizer surge line model of the AP1000 with air-water and steam-water at atmospheric pressure. They changed the water level in the pressurizer simulator up toh = 900 mm, and it affected CCFL. The effects ofh on CK andm in their experiments [6, 7] showed trends similar to those by Doi et al. [10] shown in Figure 4(b). Therefore, Figure 4(b), Figure 8(a), and Table 3 may be useful, when we evaluate data by Yu et al. [6, 7].Figure 9
System characteristics of CCFL in a pressurizer surge line.


## 4.4. Technical Issues

Gas-liquid countercurrent flows in vertical pipes have been studied for over a half century, and it is well known that CCFL characteristics strongly depend on the shapes at the top and bottom ends (Wallis [15]; Bankoff and Lee [22]). We (Yamamoto et al. [23]) classified CCFL in vertical pipes depending on the limiting locations at the square bottom end (CCFL-L) and the square top end (CCFL-U) and inside the vertical pipe (CCFL-P), and we found that the characteristic lengths defined by (3) areD,L, andD0.5L0.5 for CCFL-L, CCFL-U, and CCFL-P, respectively. However, (1)-(4) only show effects of primary parameters, and actual CCFL characteristics are very complex and depend on some secondary parameters such as the diameters (cf. Table 3), water levels, and fluid properties. Moreover, experimental data of CCFL differ among experimental facilities and conditions, and we have to consider some uncertainty among experimental data. Experiments at high pressure and temperature are few in number and effects of fluid properties as secondary parameters are still not clear. Wallis [15] showed that the liquid viscosity (which is not included in (1)-(4)) affects CCFL and he used the dimensionless viscosity number. We (Murase et al. [24]) used the viscosity ratio of the gas and liquid phases for CCFL-L, CCFL-U, and CCFL-P, but Ilyukhin et al. [18, 19] used the density ratio of the gas and liquid phases for CCFL with the protruding top and bottom ends. Therefore, CCFL correlations, which include effects of fluid properties, should be carefully used considering the application range. As the uncertainty of CCFL will increase on increasing the number of experimental data due to the inherent differences among the experimental facilities, the CCFL correlation should be derived using experimental data suitable for the system of concern.Moreover, CCFL characteristics from integration tests expressed by (10) are different from those of (9) as shown in Figure 9. Therefore, engineering judgement is required to choose an appropriate CCFL correlation. At present, we recommend using the correlation equation (10), for the actual pressurizer surge line.


## 5. Conclusions

In this study, we evaluated effects of diameters on CCFL at the square top end in vertical pipes with and without the protruding end in the upper tank by using existing air-water data in the diameter range ofD = 19-250 mm and with the water level in the upper tank up toh = 600 mm. The results we obtained are as follows.(1)
We found that there was a strong relationship between the constantCK and the slopem in the Wallis-type CCFL correlation except for small diameters ofD = 19-25 mm with the water level at midlevel, where CK increased with increasingm.(2)
On increasing the water level in the upper tank,CK andm in the Wallis-type CCFL correlation increased forD = 30 mm. However, effects of the water level were complex and depended on the range of diameters and water levels.(3)
CCFL characteristics at the square top end were expressed by the Kutateladze parameters, and the constant and slope in the Wallis-type CCFL correlation wereCK = 1.53±0.11 andm = 0.97, respectively, in the diameter range ofD = 30-250 mm. The CK values were smaller than CK = 1.53 in the diameter range ofD = 19-25 mm.The remaining technical issues for CCFL characteristics at the square top end are effects of fluid properties and application of the results to evaluate CCFL characteristics in the pressurizer surge line.



---

*Source: 1426718-2018-12-02.xml*