NUREG/IA-0100

* International

'7 Agreement Report

Assessment of CCFL Model of

RELAP5/MOD3 Against

Simple Vertical Tubes

and Rod Bundle Tests

Prepared by

S. Cho, N. Arne, Korea Electric Power Corporation, Research Center

B. D. Chung, H. J. Kim, Korea Institute of Nuclear Safety

Korea Electric Power Corporation, Research Center

117-15 Yongjeon-Dong

Dong-Gu

Taejon, Korea

Korea Institute of Nuclear Safety

P.O. Box 16

Daeduk-Danji

Taejon, Korea 305-353

Office of Nuclear Regulatory Research

U.S. Nuclear Regulatory Commission

Washington, DC 20555

June 1993

Prepared as part of

The Agreement on Research Participation and Technical Exchange

under the International Thermal-Hydraulic Code Assessment

and Application Program (ICAP)

Published by

U.S. Nuclear Regulatory Commission
NOTICE

This report was prepared under an international cooperative

agreement for the exchange of technical information. Neither

the United States Government nor any agency thereof, or any of

their employees, makes any warranty, expressed or implied, or

assumes any legal liability or responsibility for any third party's

use, or the results of such use, of any information, apparatus product or process disclosed in this report, or represents that its use

by such third party would not infringe privately owned rights.

Available from

Superintendent of Documents

U.S. Government Printing Office

P.O. Box 37082

Washington, D.C. 20013-7082

and

National Technical Information Service

Springfield, VA 22161
NUREG/IA-0100

International

Agreement Report

Assessment of CCFL Model of

RELAP5/MOD3 Against

Simple Vertical Tubes

and Rod Bundle Tests

Prepared by

S. Cho, N. Ame, Korea Electric Power Corporation, Research Center

B. D. Chung, H. J. Kim, Korea Institute of Nuclear Safety

Korea Electric Power Corporation, Research Center

117-15 Yongjeon-Dong

Dong-Gu

Taejon, Korea

Korea Institute of Nuclear Safety

P.O. Box 16

Daeduk-Danji

Taejon, Korea 305-353

Office of Nuclear Regulatory Research

U.S. Nuclear Regulatory Commission

Washington, DC 20555

June 1993

Prepared as part of

The Agreement on Research Participation and Technical Exchange under the International Thermal-Hydraulic Code Assessment

and Application Program (ICAP)

Published by

U.S. Nuclear Regulatory Commission
NOTICE

This report is based on work performed under the sponsorship of The Korea

Advanced Energy Institute of Korea. The information in this report has been

provided to the USNRC under the terms of an information exchange agreement

between the United States and Korea (Agreement on Thermal-Hydraulic Research

between the United States Nuclear Regulatory Commission and The Korea Advanced

Energy Research Institute, May 1, 1986). Korea has consented to the publication

of this report as a USNRC document in order that it may receive the widest

possible circulation among the reactor safety community. Neither the United

States Government nor Korea or any agency thereof, or any of their employees,

makes any warranty, expressed or implied, or assumes any legal liability of

responsibility for any third party's use, or the results of such use, or any

information, apparatus, product or process disclosed in this report or

represents that its use by such third party would not infringe privately

owned rights.
Assessment of CCFL Model of RELAP5/MOD3

Against Simple Vertical Tubes and Rod Bundle Tests

Sungjae Cho and Namsung Arne

Korea Electric Power Corporation, Research Center

117-15 Yongjeon-Dong, Dong-Gu, Taejon, Korea,

Bub Dong Chung and Hho Jung Kim

Korea Institute of Nuclear Safety

P.O.Box 16, Daeduk-Danji,Taejon, Korea, 305-353

Abstract

The CCFL model used in RELAP5/MOD3 version 5m5 has been

assessed against simple vertical tubes and rod bundle tests

performed at a facility of Korea Atomic Energy Research Institute.

The effect of changes in tube diameter and nodalization of tube

section were investigated. The roles of interfacial drags on the

flooding characteristics are discussed. Differences between the

calculation and the experiment are also discussed.

A comparison between model assessment results and the test data

showed that the calculated value lay well on the experimental

flooding curve specified by user, but the pressure jump before

onset of flooding was not calculated.

iii

Acknowledgement

We are grateful to Dr. Sang-Hoon Lee, the president of the

Korea Institute of Nuclear Safety, for organizing the project.

We would iike to express our special appreciation to Dr. Jong-Hee

Cha and Soon-yun Won for their help in constructing the test

facility and for kind advice during the experiment.

V

Nomenclatures

pf= Density of Liquid

Pg= Density of Gas

/Jf= Dynamic Viscosity of Liquid

P9g= Dynamic Viscosity of Gas

a = Surface Tension of Liquid

g = Universal Acceleration Constant

D = Equivalent Hydraulic Diameter = 4A/P

Dp =Nondimensional Diameter = DH/[a/{g(pf-Pg )1 1 1

/ 2

K* = Kutateladze Number of Gas Flux = p. 5 a V /[9g(pf-p )1 0 .

25

= Kutateladze Number of Liquid = pf af V/[ga(•-p)J 0 2 s

= Nondimensional Gas Flux = p a V /[gDH(pf-p )10 5

* = Nondimensional Liquid Flux = pfafVf/[gDH )0. 5

P* = Nondimensional Pressure Gradient = IdP/dzl/[g(pf-p9 )]

Re£ = Film Reynolds Number = 4F//pf where r denotes film mass flowrate

per unit width

Vii

Contents

Abstract -------------------------------------------------- iii

Acknowledgement ------------------------------------------- v

Nomenclature ---------------------------------------------- Vii

List of Tables --------------------------------------------- X

List of Figures -------------------------------------------- i

1. Introduction -------------------------------------------- 1

2. Sample Problem and Description of Test Facility ------------ 2

2.1 Sample. Problems ---------------------------------------- 2

2.2 Rod Bundle Tests --------------------------------------- 3

3. Code and Input Model Description -------------------------- 8

3.1 The Code Features -------------------------------------- 8

3.2 The Calculation Input Model ------------------------------ 9

4. The Base Case Calculation -------------------------------- 17

4.1 Simple Vertical Tubes ----------------------------------- 17

4.2 Rod Bundle Tests -------------------------------------- 19

5. Sensitivity Results and Discussion ------------------------- 41

5.1 Tube Diameter ----------------------------------------- 41

5.2 Nodalization ------------------------------------------- 44

6. Run Statistics ------------------------------------------ 63

7. Conclusion --------------------------------------------- 64

References ----------------------------------------------- 66

Appendices

A. Experimental Data for Rod Bundle CCFL Test -------------- A.1

B. Input Decks for Calculations ---------------------------- B.1

ix
List of Tables

Table 5.1 Test Cases for the simple'vertical tubes ----------- 41

Table 6.1 CPU time required in the simulation --------------- 63

X
List of Figures

Fig. 2.1 Schematic Diagram of Air-Water Counter-current Flooding

Test Loop with 3x3 Tube Bundle ---------------------- 6

Fig. 2.2 Three Type of 3x3 Tube Bundle Test Section ----------- 7

Fig. 3.1 RELAP5 Nodalization for Air/Water Test Problem -------- 13

Fig. 3.2 RELAP5 Nodalization for Type-1 Test Bundle ------------ 14

Fig. 3.3 RELAP5 Nodalization for Type-2 Test Bundle ------------ 15

Fig. 3.4 RELAP5 Nodalization for Type-3 Test Bundle ------------ 16

Fig. 4.1 Air Inlet Flowrate at the Junction 4000000 -------------- 24

Fig. 4.2 Downward Liquid Flowrate at Junction 4000000 ---------- 24

Fig. 4.3 Flow Regime Change at the Test Section Junction ------- 25

Fig. 4.4 Gas Void Fraction at the Test Section Junction --------- 25

Fig. 4.5 Interfacial Drag Coefficient at the Test Section Junction- 26

Fig. 4.6 Pressure Drop in the Test Section -------------------- 26

Fig. 4.7 Mass Flux through the Junction 4000000 ---------------- 27

Fig. 4.8 Calculated Flooding Data at the Junction 4000000 -------- 27

Fig. 4.9 Air Inlet Flowrate at the Junction 4000000 -------------- 28

Fig. 4.10 Downward Liquid Flowrate at Junction 4000000 ---------- 28

Fig. 4.11 Flow Regime Change at the Test Section Junction ------- 29

Fig. 4.12 Gas Void Fraction at the Test Section Junction --------- 29

Fig. 4.13 Interfacial Drag Coefficient at the Test Section Junction- 30

Fig. 4.14 Pressure Drop in the Test Section -------------------- 30

Fig. 4.15 Mass Flux through the Junction 4000000 ---------------- 31

Fig. 4.16 Calculated Flooding Data at the Junction 4000000 -------- 31

Fig. 4.17 Air Inlet Flowrate at the Junction 4000000 -------------- 32

Fig. 4.18 Downward Liquid Flowrate at Junction 4000000 ---------- 32

Fig. 4.19 Calculated Flooding Data at the Junction 4000000 ------- 33

Fig. 4.20 Pressure Drop in the Test Section -------------------- 33

Fig. 4.21 Flow Regime Change at the Test Section Junction ------- 34

Fig. 4.22 Gas Void Fraction at the Test Section Junction --------- 34

Fig. 4.23 Interfacial Drag Coefficient at the Test Section Junction- 35

xi
Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

Fig.

4.24

4.25

4.26

4.27

4.28

4.29

4.30

4.31

4.32

4.33

4.34

5.1

5.2

5.3

5.4

5.5

5.6

5.7

5.8

5.9

5.10

5.11

5.12

5.13

5.14

5.15

5.16

5.17

5.18

5.19

5.20

5.21

Air Inlet Flowrate at the Junction 4000000 ------------- 35

Downward Liquid Flowrate at Junction 4000000 --------- 36

Calculated Flooding Data at the Junction 4000000 ------- 36

Pressure Drop in the Test Section -------------------- 37

Flow Regime Change at the Test Section Junction ------ 37

Gas Void Fraction at the Test Section Junction --------- 38

Interfacial Drag Coefficient at the Test Section Junction-38

Calculated Flooding Data at the Junction 3030000 ------- 39

Pressure Drop in the Test Section 39

Calculated Flooding Data at the Junction 3010000 ------- 40

Pressure Drop in the Test Section ------------------- 40

Calculated Flooding Data without CCFL option ( D = 0.5 cm )

Pressure Drop in the Test Section ( D = 0.5 cm ) ------ 45

Calculated Flooding Data with CCFL option (D = 0.5cm) - 46

Pressure Drop in the Test Section ( D = 0.5 cm ) ------ 46

Calculated Flooding Data with CCFL option (D = 1.0cm) - 47

Pressure Drop in the Test Section ( D = 1.0 cm ) ------ 47

Calculated Flooding Data with CCFL option (D = 1.0cm) - 48

Pressure Drop in the Test Section C D = 1.0 cm ) ------ 48

Calculated Flooding Data without CCFL option (D=3.0cm) 49

Pressure Drop in the Test Section C D = 3.0 cm ) ------ 49

Calculated Flooding Data with CCFL option (1 = 3.0cm) - 50

Pressure Drop in the Test Section ( D = 3.0 cm ) ------ 50

Calculated Flooding Data without CCFL option (D=5.Ocm) 51

Pressure Drop in the Test Section ( D = 5.0 cm ) ------ 51

Calculated Flooding Data with CCFL option (D = 5.0cm) - 52

Pressure Drop in the Test Section ( D = 5.0 cm ) ------ 52

Calculated Flooding Data without CCFL option (D=10.cm) 53

Pressure Drop in the Test Section ( D = 10. cm ) ------ 53

Calculated Flooding Data with CCFL option (D = 10.cm) - 54

Pressure Drop in the Test Section ( D = 10. cm ) ------ 54

RELAP5 Fine Nodaiization for Air/Water Test Problem --- 55

xii
Fig. 5.22 RELAP5 Fine Nodalization for Type-1 Test Bundle ----- 56

Fig. 5.23 Calculated Flooding Data with Fine Nodalization ------- 57

Fig. 5.24 Pressure Drop with Fine Nodalization ----------------- 57

Fig. 5.25 Flow Regime Change at the Test Section Junction ----- 58

Fig. 5.26 Interfacial Drag Coefficient at the Test Section Junction 58

Fig. 5.27 Gas Void Fraction at the Test Section Junctton ------- 59

Fig. 5.28 Mass Flux through the junction 4000000 -------------- 59

Fig. 5.29 Calculated Flooding Data with Fine Nodalization ------- 60

Fig. 5.30 Pressure Drop with Fine Nodalization ----------------- 60

Fig. 5.31 Flow Regime Change at the Test Section Junction ----- 61

Fig. 5.32 Interfacial Drag Coefficient at the Test Section Junction6l

Fig. 5.33 Gas Void Fraction at the Test Section Junction ------- 62

Xiii

1. Introduction

RELAP5/MOD2 [1] is an advanced computer code developed at the

Idaho National Engineering Laboratory (INEL) for thermal-hydraulic

analyses of Pressurized Water Reactor (PWR) system transients. It

was produced by improving and extending the modeling base that had

been established with the release of RELAP5/MOD1 [2] in December

1980, and was released as a 'frozen' code in April 1985. Numerous

assessments have been carried out worldwide under the

International Code Assessment and Application Program (ICAP), and

a number of deficiencies in the code [3] have been identified.

Evaluation of these deficiencies and their regulatory significance

was providing the basis for improvements for a new version of

RELAP5, namely RELAP5/MOD3. The RELAP5/MOD3 code was frozen at

version 5m5 and released in May 1990. The code improvements were

made in many areas including a CCFL model [4] and a new interphase

drag model [5].

In the present study the CCFL model of RELAP5/MOD3/5m5 is

assessed against *a simple tube as a sample problem. It is also

assessed against the rod bundle tests which were carried out in a

test facility at the Korea Atomic Energy Research Institute (KAERI)

to check interactions between the CCFL model and interphase drag

model under countercurrent condition in rod bundle.

Descriptions of sample problems and test facility are provided

in Chapter 2. A description of code and input model is given in

Chapter 3. Various assessments against the sample problems and

rod bundle tests are provided in Chapter 4. Sensitivity studies

including fine nodalization model and RELAP5/MOD2 calculation are

discussed in Chapter 5. Run statistics are given in Chapter 6.

General conclusions are drawn in Chapter 7.

Appendix A contains the time averaged data for the rod bundle

CCFL experiments at KAERI. Appendix B contains the complete set

of input lists.

1
2. Sample Problems and Description of Test Facility

The phenomenon of Counter Current Flow Limitation (CCFL) has

been extensively studied experimentally. Considering the

geometries and layouts of the nuclear reactor, there could be

some different CCFL conditions which need careful examinations to

identify the limiting condition. They are at steam generator

tubes, at the portion of the hot legs, at annulus of the vessel

downcomer, at the upper tie plate of core and the fuel bundle.

Most of experimental works have been performed in simple tubes

using air/water flow, thus simple tube tests were selected as

sample problems. The rod bundle CCFL tests have been completed by

Korea Institute of Nuclear Safety (KINS) under contract to Korea

Electric Power Co. (KEPCO) at the end of 1989. The tests were

also selected as assessment tests to check the CCFL model and

interphase drag model under countercurrent flow condition in rod

bundles.

2.1 Sample Problems

An extensive experiment was performed by Wallis and Makkenchery

[6] with various diameter pipes. In the experiments the liquid

was drained into the the test tube with various liquid entry

design and the air was introduced by a surge tank into the test

tube. The experimental results can be correlated by the well

known Wallis correlation.

ig •i* 1/2+ +mijf .* 1/2 = C (2.1)

The experimental constant varied due to the channel inlet

shape, length or viscosity effect. For most cases, m has a value

between 0.8 to 1.0 and C ranges from 0.7 to 1.0. Simple vertical

tubes were chosen as sample problems to assess the RELAP5 model.

2
The inside diameter of tube is varied from 0.5 cm to 10 cm to

check applicability of RELAP5 CCFL model.

2.2 Rod Bundle Tests

The rod bundle CCFL tests have been completed by Korea

Institute of Nuclear Safety (KINS) under contract to Korea

Electric Power Co. (KEPCO) at the end of 1989. The tests were

also selected as assessment tests to check the CCFL model and

interphase drag model under countercurrent flow condition in rod

bundles.

The test were performed using air and water flow under the

atmospheric pressure condition. The test facility consists of

test section, water and air supply system, associated sensors and

devices to measure flowrates of water and air, differential

pressure, water levels, and data recording systems as shown in

Fig. 2.1.

The major components of the test section, as shown in Fig.2.2,

are (1) a lower plenum, (2) 4x4 cm square channel (145 cm in

height) made of transparent lucite to allow visual observations,

and (3) an upper plenum. The lower plenum has two inlets, one for

water and the other for air. The two inlets are designed in such a

way that direct interaction between the incoming air and water and

the effects of entrance and entrainment at the inlet can be

avoided. The 3x3 stainless steel tube bundle is placed vertically

at the center of the lucite square channel. The tube array has the

same geometrical dimensions of the typical 17x17 PWR fuel bundle.

The outside diameter of each tube and pitch are 0.95 cm and 1.26

cm, respectively.

In order to obtain experimental data on the effects of spacer

grids and multi-channel flow interactions on the CCFL, three

tube-bundle types were used as shown in Fig.2.2 : The first type

3
is a 3x3 tube bundle with no spacer grid, the second type is a 3x3

tube bundle with one spacer grid attached at 68.3 cm axially from

the top (at "B" in Fig.2.2), and the third type is a 3x3 tube

bundle with two spacer grids attached at 17.5 cm and 68.3 cm from

the top respectively. For convenience in discussion, these

bundles are designated here as type-1 (no spacer-grid), type-2

(one-spacer grid), and type-3 (two spacer grids), respectively.

The grid assembly consists of individual slotted straps

interlocked and brazed in an egg-crate arrangement to join the

straps permanently. The test section was equipped with two

pressure taps. An air-water separator for measuring the

entrainment rate and two water level gauges for measuring the film

downflow and entrainment rates were also provided.

The water is first introduced into the inlet located at the

bottom of the lower plenum, then it rises up to the top of the

tube bundle flowing inside of each tube. The water flowing out of

the nine tubes at the top (point "A" in Fig.2.2) is then directed,

via a nozzle, to flow in reverse direction downward by gravity

over the outer surface of the 9-tubes forming a continuous annular

water film. The nozzle is placed 30 cm above the top inlet (point

"B" in Fig.2.2) of the test channel to form a fully developed

flow. The more detailed description for the test facility is

given in Reference [7] and [8].

The flooding data were obtained by varying the air flow rate

while the water .flow rate was fixed. Major quantities measured

for each test were (1) liquid penetration rate, (2) liquid

carry-over rate, (3) air inlet flow rate, and (4) test channel

pressure drop. Actually, to obtain a single CCFL data point with a

given test bundle, from 9 to 19 runs were made increasing the air

velocity in step until the downward water film flow rate becomes

zero (i.e., the zero penetration rate of water) while maintaining

a constant water input rate. Although most of data have

oscillatory behavior, the time averaged value during 100 seconds

4
are given in Appendix A. The measurement uncertainties consists

of the error of pressure sensor, A/D ( Analog to Digital )

conversion error, and calibration error. The calculated error of

channel pressure drop is about 5 mm H20. The estimated error for

air flow measurement is about 650 cc/sec and for water flow it is

about 1.7 % of liquid injection flow.

5
1) Water Storage Tank

2) Water Pump

3) Regulating Valve

4) Water -Flow Meter

5) Drain Line

6) Test Section

7),8) Pressure Tap

13) Air Flow Meter

14) Air-Water Separatoi

15) Air Exhaust

16),17) D/P Transducer

18) A/D Converter

19) Personal Computer

1 2 •/-) 

... _l1 13, . ,

Fig. 2.1 Schematic Diagram of Air-Water Counter-current Flooding

Test Loop with 3x3 Tube Bundle
9

T 000 000 %.d H

0 l.e14

3x3 Tube Bundle

Cross Section

-Water Nozzle

at "A"

PJ -Water Inlet

"A"

4- T

-J

Unit mm

Fi. rst,

Grid

Second

Grid

L C.

"Tube Bundle

with Twb Spacer

Grids (Type-3)"

"Tube Bundle

with One Spacer

Grid (Type-2)"

Air Flow I I i Air Flow

Water Flow

"Tube Bundle

with No Spacer

Grid (Type-i)"

Fig. 2.2 Three Type of 3x3 Tube Bundle Test Section
3. Code and Input Model Description

The assessment calculation for sample problems and rod bundle

CCFL experiment was done using RELAP5/MOD3/5m5, code received at

the beginning of May 1990. The descriptive document available for

this code at the time of preparing this report are scheduled to be

released later. However document of a CCFL model and a new

interphase drag model tested in this assessment are available in

reference [4] and (5].

3.1 Code Model Description

A general CCFL model, which was proposed by Bankoff [9], was

implemented in the 5m5 version of RELAP5/MOD3 that allows the user

to select the Wallis form, the Kutateladze form, or a form in

between Wallis and Kutateladze forms. The general form using

dimensionless flux is given by

H 1/2+ m H 1/2 = C (3.1)

Hg9 gV g[Pg /gw(Pf-Pg)]W/2

Hf C- fVf[pf/gw(pf-p g)]

where

w w=jD 1Dpi

I = [./9(Pf-Pg)] 11 2

In the above equation, p can be a number from 0 to 1. For P=0,

the Wallis form of CCFL equation is obtained, and for P=1, the

Kutateladze form of CCFL equation is obtained. The user should

make a decision for the input values of m,C,and p using the

appropriate experimental results.

With regard to the solution method, if the CCFL model is

8
requested by the user, the coding checks if counter current flow

exits and if the liquid downflow exceeds the limit imposed by

eq.(2.2). If this is true, the sum momentum equation and flooding

equation (2.2) are applied as the phasic momentum equations. Thus

the difference equation, which contains the Interfacial friction,

is replaced with flooding limit equation (2.2).

Although the RELAP5/MOD3 has a user controlled junction

flooding model for flooding at a singularity in the channel

geometry, the counter current flow is much related to interfacial

drag model. As a part of UK's contribution to the code

improvement, the Interfacial friction model suggested by CEGB [5]

was implemented in a released version of RELAP5/MOD3. The new

interfacial friction model was developed for the vertical bubbly

and slug flow regimes, and a CCFL correlation appropriate for a

straight uniform flow channel was also adopted for the calculation

of drift flux used for interfacial friction models.

3.2 The Calculation Input Model

3.2.1 The Sample Problems

The sample problems used in the present calculation is closely

related to the test facilities used in many previous experiment

for a vertical single tube. The specific nodalization is given by

Fig.3.1.

Most experiments were performed using air and water. The water

was drained into the test section from the top plenum, forming

downward film flow. Air was introduced by a surge tank, flowing

up into the test tube. The liquid inlet was modeled as a time

dependent junction to specify the constant liquid injection. The

air inlet junction from a surge tank was modeled as single

junction with a smooth option to avoid the inlet contraction

9
effect of phasic momentum equation. However, Most experimental

data show that the flooding occurs at the air inlet junction due

to its geometric singularity, and thus the CCFL option was used

at the air inlet junction (JUN-4000000). The air inlet flow was

modeled as a time dependent junction. For most experiments the

experimental data have oscillatory behavior during- the flooding

process, and thus the flooding data are presented as the time

averaged value during a test period for a given air and liquid

injection flow. A steady state calculation for a given air and

liquid injection flow will be very time consuming work because

there are many different flooding points. Therefore as test

sample problems, the liquid flow was fixed and the air flow was

increased slowly until the flooding is completed. After some

trial and error calculations, we find a 200 second duration is

enough to neglect the transient effect of increasing air flow.

For sample problems, the code input of the flooding curve was

specified as follows ; the simple Wallis form.

. +1I = 1 (3.2)

In view of the geometries and layout of nuclear reactor, the

specific flooding location for a single tube may be steam

generator U-tubes. A typical inside diameter of the PWR tubes is

about 1.7 cm. Thus a flooding test in 2 cm diameter tubes was

selected as a base test case, and calculations for 5 different

tube diameters ( 0.5 cm, 1 cm, 3 cm, 5 cm, and 10 cm) have been

performed as sensitivity studies. The results were discussed in

Section 5.1.

The test section pipe was modeled as 6 uniform volumes and

liquid drain tank was modeled as a relatively large single volume.

In order to compare calculated values with flooding curves

specified by user, the liquid velocity and air velocity at the

specific junction were converted to nondimensional gas flux using

10
junction void fraction and velocity as "control variables" as

shown below ;

J* = p1/2 (xV/[gD(pf-p)]W/ 2 

-_ 0.0112 D-1/2a V (3.3) J 9 V ' 9 gi gi

Jf = Pf 2afVf/[gD(pf-p )Wl/2 L 0.3195 D- /2

f f 9 afj YfJ

The physical properties of each phase were calculated at the

experimental condition 1 atm and 280 K.

3.2.2 Rod Bundle Tests

The specific nodalization and modeling were much similar to the

sample problems as shown in Fig. 3.2. But the rod bundle form of

the new interfacial friction model was activated by specifying the

third flag (b) of an extended volume control flag, bfe, to be

unity in the pipe volume of test section. The model geometry and

elevation were related to the test facilities described in section

2.2. In the rod bundle CCFL experiment, three test bundles were

used, and it was found experimentally that the flooding location

is "air inlet" for type-1 test, "spacer grid" for type-2 test, and

"first spacer grid" for type-3 test. Therefore three different

locations for CCFL junction were examined as shown in Fig. 3.2,

3.3, and 3.4 corresponding to each type of test section.

For moderate high liquid flux, the experimental flooding curves

of Kutateladze form are as follows;

K*( 1/2 + *K 1/2 D* 1/4(

* Kf = D= 1.44) for type-1 (3.4)

1/2 + K* 1/2 = 1.38 for type-2 9 f

1/2 + K* 1/2 = 1.23 for type-3 9 f

Therefore the code input value for flooding curve was specified

11
as above equations (Eq.3.4) for each type of tests.

In order to compare the calculated values with flooding curves

specified by user, the liquid velocity and air velocity at the

specific junction were converted to nondimensional gas flux using

junction void fraction and velocity as "control variables" as

shown below;

= P9 g gg / Pf-pg L,- 0.214 a V (3.5)

Kf = pf afVf/[g (pj )]/

Kf p /: rf ffcrp-p) 1 1 6.1 ctfjVfj

12
A

a

g

C

C

C

C,

4

I

I3v-' I TDv-1 ELEVATION

WVater Inlet TDJ-2 0.00 M

-0.30

. --- 0.60 M

I

U

-- -0.90 m

U

-1.50 r

CCFL Test JUN-4 --- 1.80 n

Low Plenum v

SNV-5 -/

Air Intet TDJ-6 i I - -1.93m

ILiquid TankN- lAir TV- Source 

-- 2.43 m

NodalUzo'ton for Air/Vacer Test Problem

Fig. 3.1 RELAP5 Nodalization for Air/Water Test Problem

13
TDJ-2

J,

:3

U

a)

1AW,

ELEVATION

- 0.00 n

- -02.5 ft

-- 0.50 ft

-- 0.75 r

-- 1.00

---1.25 a

-- 1.50 n

-1.-m3 r

-2.13 m

Nodoazation for Type-1 Test Bundte

Fig. 3.2 RELAP5 Nodalization for Type-1 Test Bundle

14
Fig. 3.3 RELAP5 Nodalization for Type-2 Test Bundle

15
DV-11 j I TDV-1 I ELEVATION

V:ater Inlet TDJ-2 - 0.00 m

-025 m

:• 3-02

C 0 -- 0.50 m

• 3-03 -p

D/- --- 0.75 m

• P 3-04 -J

M LA-

- - .? --- 1.00 M

3-05

--- 1.25 m

3-06

-- 1.50 M

SLow Plenum SNV-5 •lil %T]_

Air Inlet TI3J-6 I --- 1.63 m

Liquid Tank Air Source

SNV-9 ITDV-7

-- 2.13 rm

Nbdo•lizotion for Type-3 Test BundLe

Fig. 3.4 RELAP5 Nodalization for Type-3 Test Bundle

16
4. The Base Case Calculation

4.1 Simple Vertical Tubes

As discussed in section 3.2.1, the air/water countercurrent

flow in a tube of 2 cm inside diameter was simulated with

RELAP5/MOD3/5m5. Two calculations were performed with and without

a CCFL junction option. Two results were compared in terms of air

inlet flow and liquid downward flow at a CCFL Junction, flooding

cu rve, total mass flux, gas voi d fraction, flow regime,

interfaciai friction,and pressure drop in the test section. The

time step size was fixed as 0.01 second during the whole

transient. If the time step size became larger than 0.01 seconds,

it was experienced that the calculation was unstable or failed in

some cases.

Fig. 4.1 shows the gas flowrate at junction 4000000 with no

CCFL option. The negative value means the upward direction as

shown in Fig. 3.1. The gas flow was increased slowly during 200

seconds until the flooding was completed. During the first 10

seconds, the steady state calculation was performed to get the

fully developed free falling annular film flow at zero air inlet

flow. Although the air inlet flow through junction 6000000 was set

to increase monotonously, the air flow into the test section

through junction 4000000 had a little fluctuation at the start of

restriction of liquid down flow (' 120 sec.) and at the end of the

flooding (~ 170 sec.). Fig 4.2 shows the restriction of the

downward liquid flow. This restriction comes from the increase of

the interfacial drag as air flow increases. As shown in Fig.4.3,

the junction flow regime is annular-mist( regime 6) over the whole

flooding range and becomes dispersed flow regime (regime 7) as the

gas void fraction becomes nearly to unity ( > 0.9999). The gas

void fraction is shown in Fig.4.4. As shown in Fig.4.4 the liquid

film thickness increase with increasing air flow, and then

17
decrease after the flooding occurs. The liquid sweep from the

bottom to top section was calculated. The trends of the

interfacial friction, shown in Fig. 4.5, has the same trends of

liquid void fraction, because the interfacial friction factor,

which is described by a Bharathan correlation [10), increases with

increasing film thickness. The channel pressure drop trend is

similar to that of liquid void fraction, as shown in Fig. 4.6.

The new interfacial friction which is based on drift flux ,V gmay

not be applied because the mass flux is lower than 50 kg/m.sec

over most of the flooding range as shown in Fig. 4.7. If we draw

a flooding curve using the upward gas flux and the downward liquid

flux values calculated during the test transient, it is possible

to compare the calculation with the experimental results obtained

in single tube tests. Most experimental results can be correlated

by the well known Wallis correlation as discussed in section 2.2.

The experimental constant in Eq.(2.1) varied due to the channel

inlet shape, length or viscosity effect. For most cases, the

slope (m) has a value between 0.8 to 1.0 and the constant (C)

ranges from 0.7 to 1.0. However the slope of calculated flooding

curve (m) is about 0.5, although the CCFL point (C) is about 1.0,

as shown in Fig.4.8. Thus the flooding occurs at much higher air

velocity ( j = 0.7 ) than the most experimental data. It may be 9 due to lack of CCFL model applicability at the geometric

singularity in test section.

Therefore the CCFL option was applied at the junction 4000000

(air inlet). A simple Wallis form (Eq.3.2) was selected as a

input value for the flooding curve.

The inlet air flow and downward liquid flow are fluctuated

continuously after the flooding was reached (' 70 seconds) as

shown in Fig.4.9 and 4.10. The flow regime in the test section

changes between annular-mist (regime 6) and slug (regime 5) during

the flooding process, and becomes the dispersed flow regime

(regime 7) as the gas void fraction becomes near unity ( >

18
0.9999). The gas void -fraction In the test channel is shown in

Fig.4.12. The bottom junction restricts the downward flow after

70 seconds by the CCFL curve provided by user, and thus the liquid

fraction in upper junctions increases to nearly 0.2. The

increased liquid fraction results the change of flow regime; slug

regime from annular mist flow. The interfacial friction, As shown

in Fig. 4.13, increases sharply in slug flow regime. Thus the

liquid carry-over occurs through the whole test channel. It

reduces the liquid void fraction again, and the above cycling

process is continued until the flooding is completed. The channel

pressure drop, shown in Fig. 4.14, increases sharply after the

onset of flooding (70 sec.) and decreases slowly during the

flooding process. These oscillatory behavior and pressure trends

are much similar to the observed phenomena in most CCFL

experiments. The total inlet mass flux also oscillates during the

flooding. Although the air inlet flow and the downward liquid

flow oscillate, their oscillations lie near the CCFL curve

specified by user to satisfy the difference momentum equation of

RELAP5, as shown in Fig. 4.16.

4.2 Rod Bundle Tests

As described :n section 3.2.2, the calculations for each test

type were performed and compared with the experimental data. The

input data for calculations were prepared based on the specific

experimental condition. The liquid injection flow rate was fixed

to 2600 cc/min. The corresponding mass flux is about 35 kg/m 2 sec

in the case of no air flow. The gas flow was increased slowly

enough to neglect the transient effect of increasing the air flow.

4.2.1 Type-1 Test Bundle

19
For test type-i, two calculations were performed with or

without CCFL option. Two calculation results were compared with

experimental data obtained from the test facility described in

section 2.2.

The air inlet flow was increased during 200 seconds until the

flooding is completed, as shown in Fig. 4.17. The downward liquid

film flow was restricted at 110 seconds by the increased

interphase drag force (Fig. 4.18). A little fluctuation occurs at

the early stage of the flooding. The calculated flooding curve

using the above air flow and liquid flow at junction 4000000 was

presented in Fig. 4.19 as a Kutateladze form. The calculated

flooding occurs at much higher air velocity than the experimental

data ( more than 20 % ). The complete flooding occurs at K*= 3.6, g

although it was not presented in figure. The channel pressure

drop was presented in Fig. 4.20, with respect to the

nondimensional air velocity (K*) in order to compare the 9

experimental data. As shown in the figure, the pressure

continuously increases with increasing air flow, while the

experimental data shows the sudden jump after the onset of

flooding. The flow regime changes from annular mist flow (mode

6) to slug flow (mode 5) as shown in Fig. 4.21. The gas void

fraction decreases with increasing air flow rate and maintained

nearly constant during the flooding process (Fig. 4.22). The

interfacial friction factor increases somewhat slowly, not

dramatically (Fig. 4.23). However the interphase drag force was

sufficient to carry over the water flow due to the high air

velocity. All discrepancies may be resulted from the lack of

applicability of CCFL model at the geometric singularity in test

section.

In the experiment of type-1 test, the flooding occurs at the

air inlet junction. Therefore the experimental flooding curve

,Eq. (3.4), was specified as code input value at the air inlet

junction. The calculated air inlet flow and the downward liquid

20
flow at junction 4000000 are shown in Fig. 4.24 and 4.25. These

calculated values were used to generate the flooding curve, i.e.

the Fig. 4.26. The calculated flooding values lie well on the

specified CCFL curve which was obtained from the experimental

data. The pressure drop in the test channel is shown in Fig.

4.27. The experimental data of the liquid injection 2600 cc/min

are also presented in the figure. In the experiment, it was shown

that the static pressure drop in the channel is nearly

proportional to the liquid void fraction. The pressure drop

increased suddenly before the onset of flooding and oscillated

continuously. After the liquid carry over, the pressure decreased

again with decreasing liquid fraction. There are two

discrepancies between calculation results and experimental

results. First, a sudden jump occurred at a higher gas velocity

compared with experiment. In the experiment, the liquid downward

flow was restricted by the singularity at the air inlet with K*=

g 0.9, and the liquid fraction increased. Although the liquid flow

oscillated due to flow restriction, the total time averaged liquid

flow remained constant. Sometimes it was called "the onset of

liquid mixing". It resulted a sudden increase of pressure drop

before onset of flooding. In the simulation a sudden pressure

drop occured after "the onset of flooding" ( K =1.1), because the

liquid holdup in the channel is possible after "the onset of

flooding". Second, the peak pressure drop was lower than the

experimental results and increased with the air flow. It may be

resulted from the interphase drag for rod bundle. The transition

to slug flow (mode 5) from annular mist flow (mode 6) occurred

after the onset of flooding, as shown in Fig. 4.28, and thus a

sharp increase of interfacial friction coefficient was also

calculated (Fig. 4.29). The liquid film thickness increased with

increasing air flow and jumped suddenly at the point of flooding

due to restriction of the downward liquid flow at the air inlet

junction (Fig. 4.30). Although the liquid fraction decreased

21
slowly during the flooding process, the calculated pressure drop

increased due to increase of interphase drag force with increase

of air flow as shown in Fig. 4.27.

4.2.2. Type-2 and Type-3 Test Bundle

In the experiment for the other type of test bundle, the main

characteristics of flooding phenomena are almost identical to the

characteristics of type-1 test bundle. However the flooding

locations due to the singularity of geometry were quite different.

Thus a different location of CCFL junction was examined as

presented by the nodalizations shown in Fig. 3.3 and 3.4,

following the experimental observation. The flooding curve was

specified according to the equation (3.4) which was determined by

the experiments.- The calculated flooding curve at the location

of the spacer grid ( junction 3-030000 ) was compared for type-2

test bundle in Fig. 4.31. The calculated values lie well on the

curve specified by user and the channel pressure drop during

flooding process agrees well with the experimental data, as shown

in Fig. 4.32. But a sudden jump of pressure occurs just after the

onset of flooding due to the same reason described in the previous

section.

The calculated flooding curve at the location of the first

spacer grid ( junction 3-010000 ) was compared with experimental

data for type-3 test bundle also. The calculation results lie

well on the curve specified by user, but the results of CCFL at a

higher air velocity than the complete flooding point are somewhat

curious (Fig. 4.33). It may be resulted from a wrong

implementation or coding error of flooding model in 5m5 version of

RELAP5/MOD3 code. Thus no further discussions for the curious

results are given here. The pressure drop agreed well with

experimental data, except for the point of sudden jump, as shown

in Fig. 4.34. The curious pressure increase after the complete

22
flooding was considered to be a result of coding error of 5m5

version also.

23
RELAP5/MOD3/5M5 without CCFL Model

Test Cme (Jfeod - 0.16. 0 - 2 Cm)

C

2

0

I:

0

-0.1

-0.2

--0.3

-0.4

-0.5

-0.6

-0.7

-0.8

-O.9

-1

--1.1

--1.2

a 20 40 60 80 100 120 140 160 180 200

Fig. 4.1 Air Inlet Flowrate at the Junction 4000000

relap5/mod3/5m5 without CCFL Model

Test ~s (Limed - 0.16. 0 - 2 Cmn)

0

3

a.

0.17

0.16

0.15

0.14

0.13

0.12

0.11

0.1

0.09

0.07

O.O5

0.05

0.04

0.03

0.02

0.01

0

0 20 40 60 8o 100 120

Tkne (Sem )

140 180 180 200

Fic. 4.2 Downward Liquid FRowrate at Junction 4000000

24
E

z

E

8

7.8

7.6

7.4

7.2

7

6.8

6.6

6.4

6.2

6

5.8

5.6

5.4

5.2

5

RELAP5/MOD3/5M5 without CCFL Model

Test Cac (J.eed -- 0.1. 0 - 2 cm)

6 ~543 fiL

5 5 5 I p I I i I 5 I I i I I I

Fig. 4.3

0 20 40 60 80 100 120 140 160

•,,M, ( sec=on )

Flow Regime Change at the Test Section

(1 - jun.3010000 2 - jun.3020000 : 3 -

4 - jun.3040000 5 - jun.3050000 : 6 -

180 2W0

Junction

jun.3030000

jun.4000000)

RELAP5/MOD3/5M5 without CCFL Model

Toot c"o (J4feed - 0.16. D - 2 cn)

0.9

C

0

0.7

0. -

;::s. 4.4

0 20 40 60 80 100 120 140 160 180 200

Time (Seconds)

Gas Void Fraction at the Test Section Junction

(1 - jun.3010000 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 5 - jun.3050000 : 6 - jun.4000000)

25
RELAP5/MOD3/5M5 without CCFL Model

Test Com U amd - 0.16. D - 2 amn)

10

9

a

E

Cl

S

z

S

U

a

.3

C

7

.6

4

3

2

0 20 so a0 100 120 140 180 180 200

Fig. 4.5 Inter-facial Drag Coefficient at the Test

(1 - jun.3010000 : 2 - jun.3020000 : 3

4 - jun.3040000 : 5 - jun.3050000 : 6

Section Junction - jun.3030000

- jun.4000000)

relap5/mod3/5m5 without CCFL Model

Test Come (ljeed - 0.16. D - 2 cm)

INg

a.

0 -

CC

Ca

0

9-

7

7

6-

5-

4-

3-

2-

~~I ~ 0 .r-- I~~ ~I I I I I I I 20i a i 80 100 120 1.0 6 0 C) 20 40 60 150 100 120 140 160 180 200

Time (aO ,a )

Fig. 4.6 Pressure Drop in the Test Section

26
relop5/mod3/5m5 without CCFL Model

Týt Cae (L few! - 0.16. D - 2 cm)

0

100

90

80

70

60

50

40

30

20

10

0

-10

-20

-30

--40

-- 50

0 20 40 60 80 100 120 140 160 180

Ilfm6

200

Fig. 4.7 Mass Flux through the Junction 4000000

RELAP5/MOD3/5M5 without CCFL Model

Te Cmý (Jfd - 0.16. 0 - 2 cm)

0r

0

*1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 CA4 06

Uqtdd Flo-u ta [ (.ii)]

0.8

Fig. 4.8 Calculated Flooding Data at the Junction 4000000

27
RELAP5/MOD3/5M5 using CCFL Model

Test co (J_feed - 0.16. D - 2 cm)

4

C

0

0

0

-0.1

-0.2

-0-3

-0.4

-0.5

-0.6

-0.7

-0.8

-0.9

-1

--1.1

--1.2

0 20 40 60 80 100 120 140 160 180 200

Frtne ( Secotds )

Fig. 4.9 Air Inlet Flowrate at the Junction 4000000

relap5/mod3/5m5 using CCFL Model

Test Cose (jfeed - 0.16. D - 2 cm)

e

0

0

*0

3

0*

0.3

0.28

0.26

0.24

0.22

0.2

0.18

0.16

0.14

0.12

0.1

0.06

0.06

0.04

0.02

0

13 1

131

13 13

13

C

0 C

13 m

I-,

0 20 40 60 80 100 120 140 160 180 200

Time (seonds)

Fic. 4.10 Downward Liquid Flowrate at Junction 4000000

28
RELAP5/MOD3/5M5 using 6CFL Model

Test Come (J feed - 0.18. D - 2 cm)

E

E

X

a

7.5

7

6.5

6

5.5

5

4.5

4

b

0 20 40 60 80 100 120 140 160 180 200

Tie (Seond)

Fig. 4.11 Flow Regime Change at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 using CCFL Model

Test Ccme (J teed - 0.16. D - 2 cm)

oC

0

0,

a

L.

0.9

0.8

0.7

0.6

0.5

0 20 40 60 80 100 120 140 160 180 200

rime ( Seconds )

Fig. 4.12 Gas Void Fraction at

(1 - jun.3010000 : 2

4 - jun.3040000 : 5

the Test Section

- jun.3020000 3

- jun.3050000 :

Junction

- jun.SO300000

5 - jun.4000000)

29
RELAP5/MOD3/5M5 using CCFL Model

Test Cam (JLfd - 0.16, D - 2 Cm)

-S

E"S

6

0

0

h.

C

a

C

C

50

40

30

20

10

0

Fig. 4.13 Inte

(1

4

o 20 40 60 so 100 120 140 160 18o 200

rfacial Drag Coefficient ait ttie Test Section Junc

- jun.3010000 2 - jun.3020000 3 - jun.3030000

- jun.3040000 5 - jun.3050000 6 - jun.400000C

tion

relap5/mod3/5m5 using CCFL Model

Test CAse (Lf•,d - 0.16. 0 - 2 cm)

CL

02

C.

*0

C

10

9

8

7

6

5

3

2

1

0

0 20 40 60 80 100 120 140 160 180 200

Tnm-m (a.,orts)

Fig. 4.14 Pressure Drop in the Test Section

30
relap5/mod3/5rm5 using CCFL Model

Toet Cas (L.•e• - 0.16. D - 2 am)

CE

N,

,C

130

120

110

100

90

a0

70

60

50

40

30

20

10

0

-10

-20

i

g kh lIi

r R I I I I

4

0 20 40 60 80 100 120 140 160 180 200

hml (sooonds)

Fig. 4.15 Mass Flux through -the" Jun6ction 40000000

RELAP5/MOD3/5M5 using CCFL Model

Tea Case (4 _feed - 0.16. D - 2 an)

a--

0

2

1.1

1

0.9

0.O

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.2 0.4 0.6 0.8

LWqu Flowra[ owQtCLjJ)

Fig. 4.16 Calculated Flooding Data at the Junction 4000000

31
RELAP5/MOD3/5M5 without CCFL Model

Rod Bundle Type 1 ( 01=2600 cc/min )

V

a--

0

0

-0.1

-0.2

-0.3

-0.4

-0.5

-0.6

-0.7

-0.8

-0.9

-i

-. 2

-2.'

-2.4

- 1.,5

-1,.6

-1.9

-2..

-2.2

0 20 40 60 80 100 120 140 160 180

Time ( Seconds )

.17 Air Inlet Flowrate at the Junction 4000000

200

Fig. 4.

relop5/mod3/5m5 without CCFL Model

Rod Bundle Type1 ( 0-2600 cc/min)

0.24

0.22 -

0.2 -

7 --

:2

0.18 -

C

0

0

0

0.16

0.14

0.12

0.1

0.08

0.06 -

0.04 -

0.02

0

0 20 40 60 80 100 120 '40 160 1.80 200

Tnme (seconds)

Fig. 4.18 Downward Liquid Flowrate at Junction 4000000

32
RELAP5/MOD3/5M5 without CCFL Model

Rod Bundle Typal ( CM-2600 cc/ml)

t

a.

a

a

V

3

0

S

0

a

1.6

1.5

1.4

1.3

1.2

1.1

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

SI

fb

jqa

*1

0- I I I

0 0.2 0.4 0.6 0.8

Uquud nowru•t [ Pqrt(KI ]

Fig. 4.19 Calculated Flooding Data at-th.e Junction

( o - Experiment ; + - Calculation )

4000000

relap5/mod3/5m5 without CCFL Model

Rod Bundle Typel ( 01-2600cc/mmn )

10 .7

9.

8-

a

a.

0~.~

'-a

*0

s-Ia

C

C

a

-J

7

6

S

4

.3

2

00

T'7

0

0

, T , ; I A I . i 4

0.2 0.4 0.6 0.8 1 i.2

A go ~ K_9

* i i i

1.'. 1.6 1.8 2 2.2

Fig. 4.20 Pressure Drop in the Test Section

( n - Experiment ; + - Calculation )

33
RELAP5/MOD3/5M5 without CCFL Model

Rod Burxde Typal ( 0-2600cc/mn )

E

n

z

Ea

a

7.5

7

8.5

6

5.5

5

4.5

4

0 20 40 60 80 100 120

Th (Scond)

140 160 180 200

Fig. 4.21 Flow Regime Change at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 without CCFL Model

Rod Bundl Typal ( CA-2600 cc/mn )

0.9-

- 1=

C.

0

0

0.7 -

0.6

v.5 I i I . I i 0 20 40 150 i C I I I I

80 100 120 140 160 180

Trme ( Sed3 )

21

2.00

Fic. 4.22 Gas Void Fraction

(1 - jun.3010000

4 - jun.3040000

at the Test Section

2 - jun.3020000 : 3

5 - jun.3050000 : 6

Junction

- jun.3030000

- jun.4000000)

34
RELAP5/MOD3/5M5 without CCFL Model

Rod BundM Typal ( QI-2600 cc/mm )

t'4I

2!

C

50

40

30

20

10

0

0 20 40 60 80 100 120 140 160 180 200

Thn ( Sacnde )

Fig. 4.23 Interfacial Drag Coefficient at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 using CCFL Model

Rod Bun Type 1 ( 04-2600 cc/mr )

U,

2

0

0

-0.1

-0.2

-0.3

-0.4

-0.5

-0.6

-0.7

-0.8

-0.9

-1

-1.1

-1.2

-1.3

-1.4

-1.5

-1.6

-1.7

-1.8

-1.9

-2

-2.1

0 20 40 60 80 100 120

Tifm ( seOnda)

140 160 18o 200

Fig. 4.24 Air Inlet Flowrate at the Junction 4000000

35
relap5/mod3/5m5 usitng CCFL Model

Rod Be Tpa1 ( @-200 -C/m"*)

9-'

E

U.

0.18

0.16

0.14

0.12

041

0.06

0.06

0.04

0.02

0

0 20 40 60 DO 100 120 140 160

Time ,==c )

180 200

Fig. 4.25 Downward Liquid Flowrate at Junction 4000000

RELAP5/MOD3/5M5 using CCFL Model

Rod Bundie Typel ( 04-2600 cc/rrm )

t

a.

S

9

a

0

a

a

C

1.6

1.5

1.4

1.3

1.2

1.1

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.2 0.4 0.6 0.8

UquW Fk-rk t [ ,qrt•KJ_]

1

Fig. 4.26 Calculated Flooding Data at the Junction 4000000

( E - Experiment ; + - Calculation )

36
relop5/mod3/5m5 using CCFL Model

Rod B TYP,1 ( C-260-0c/nn )

9

4 -

(L

C

2-

++

I1-

0 I -

v I I A 4. I I .. . . . . . . . .

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.5 2

Fig. 4.27

KF[low [ KgJ

Pressure Drop in the Test Section

( o - Experiment ; + - Calculation

RELAP5/MOD3/5M5 using CCFL Model

Rod Bundle Typel ( QI-2600or/min )

)

z

0

i:-

8

7.5

7

6.5

6

5.5

5

4.5

4

Fig. 4.28

0 20 40 60 80 100 120 140 160 180 200

Time ( Second. )

Flow Regime Change at the Test Section Junction

(1 - jun.3010000 2 - jun.3020000 3 - jun.3030000

4 - jun.3040000 5 - jun.3050000 6 - jun.4000000)

37
RELAP5/MOD3/5M5 using CCFL Model

Rod Bundri Typel ( QC-2600 cc/min )

E

C

X

40

30

20

10

0

0 20 40 60 80 100 120 '140 160 180 200

Fig. 4.29 Interfacial Drag Coefficient at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 using CCFL Model

Rod Bundle Typel ( 01-2600 cc/mnm )

1

0.9

c

02

0.7

0.6

0.5

0 20 40 60 80 100 120 140 160 180 200

Twne C Seco•s )

Fig. 4.30 Gas Void Fraction

(1 - jun.3010000

4 - jun.3040000

at the Test Section Junction

2 - jun.3020000 3 - jun.3030000

5 - jun.3050000 6 - jun.4000000)

38
RELAP5/MOD3/5M5 using CCFL Model

Rod Bundl Type2 ( 01-2600 cc/ml)

V

V.

a

0

0

U

0

0

1.5

1.4

1.3

1.2

1.1

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.2 0.4 0.6 0.8 1

Uquid FoAt, E oqrt(K.. )]

Fig. 4.31 Calculated Flooding Data at the Junction 3030000

o - Experiment; + - Calculation)

relap5/mod3/5m5 using CCFL Model

Rod Bundle TypA2 ( I-2600=cc/an )

5

4

0-

'0

*0

C

C

3 -

2 -

00113 130 

0 - - . -

0

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2

A~r Flownt~e C Kg_ ]

Fig. 4.32 Pressure Drop in the Test Section

( C - Experiment ; + - Calculation )

39
RELAP5/MOD3/5M5 using CCFL Model

Pod BundMe Type3 ( CA-2600 cc/min)

0

1.5

1.4

1.3

1.2

1.1

I

0.9

0.8

0.7

0.8

0.5

0.4

0.3

0.2

0.1

0

Fig. 4.33

0 0.2 0.4 0.6 C.A

Lq -d , [: ,E ot(K.f) I

Calculated Floodirig Data at the Junction 3010000

( o - Experiment ; + - Calculation )

relop5/mod3/Sm5 using CCFL Model

-d fBA. j;.3 (.:-2oo•c.An ) 5

4

IL

2a

C 3

2

0

Fi

.0 0.2 0.4 0.W 0.8 1 1.2 1.4 1.6 1.8 2

9. 4.34 Pressure Drop in the Test Section

( o - Experiment ; + - Calculation )

40
5. Sensitivity Results and Discussion

5.1 Tube Diameter

The base calculations for the sample problem of a simple

vertical tube the inside diameter of which is 2 cm has been

performed. Although the inside diameter was selected for the

reason of the similarity of the steam generator U-tubes, there are

many other possibilities for the flooding in the tubes which have

a different size of the diameter. The different hydraulic

diameter causes the different interfacial drag and may result in

different flooding characteristics. The calculations have been

performed for a wide range of diameter size. A pair of

calculations were done for each case of a tube diameter, with and

without a CCFL junction option. The test cases and specified CCFL

curves are summarized in Table 5.1

Table 5.1 Test cases for the simple vertical tubes.

Case Tube Diameter CCFL curves specified

at junction 4000000 by user's input

1 0.5 cm J / . f = 1.0

2 1.0 cm j* jf = 1.0

**: 1/2+ .* 1/2

3 2.0 cm J* jf = 1.0 ( Base Case )

*:, 1/2+ .* 112

4 3.0 cm *g , f+ = 1.0

5 5.0 cm * 1/2j = 1.0

K•12+K /2 1/2

6 10.0 cm K 1/2+ K* 1 = 3.2 ( Large Tube ) g f

For the 0.5cm tube diameter ( case 1 ), the flooding occurs

much earlier than the flooding curve specified by user. Therefore

there are no differences between the results obtained with a CCFL

41
option and those obtained without a CCFL option, as shown in Fig.

5.1, 5.2, 5.3 and 5.4. The calculated pressure drop increased

monotonously with increasing air velocity, while most experiments

showed a sudden jump of pressure drop. The dispersed values were

due to oscillation before the reaching of the steady state for the

free falling film flow.

For the 1 cm tube diameter ( case 2 ), the flooding curve and

the pressure drop in the test channel calculated without a CCFL

option are presented in Fig. 5.5 and 5.6. With a specified CCFL

option chosen for test junction, the downward liquid flow was

restricted by the CCFL curve or the interfacial drag force in case

the interfacial drag force resulted in a lower liquid down flow.

In this case the flooding data would lie not on the CCFL curve but

on the restriction line determined by interfacial drag, in the

regions between j *./2_ 0.05 and 1/2= 0.2. The calculated f f

flooding curve and pressure drop are presented in Fig. 5.7 and

5.8. The pressure drop increase suddenly just after the flooding

condition is reached due to the increase of liquid void fraction

resulted from the restriction of downward liquid flow.

For the 2 cm tube diameter ( case 3), the results were fully

discussed as a base case in the previous section.

For the 3 cm tube diameter ( case 4), a flooding occured at a

much higher gas velocity when a CCFL junction option was not

chosen and the pressure drop increased monotonously with

increasing air velocity, as shown in Fig. 5.9 and 5.10. With a

CCFL junction option, the flooding data lay well on the CCFL curve

specified by user, and a sudden jump of pressure drop was

predicted as shown in Fig. 5.11 and 5.12. The Fig. 5.12 shows a

continuous oscillation of pressure drop during the flooding

process. The above trend of pressure drop is a phenomenon which

was observed in most experiments.

For the 5 cm tube diameter ( case 5 ), the calculated flooding

and pressure characteristics without a CCFL junction option have

42
almost similar trends of the case 4, except for the high

oscillation due to flow instability ( Fig.5.13, Fig. 5.14 ). This

flow instability might be resulted from the reduced interfacial

drag force and wall friction with increasing hydraulic diameter.

When a CCFL option was applied at the junction, the flooding data

lay on the specified flooding curve. However the high

oscillations were calculated in the region between jif = 0.1 and

f -2= 0.3. The liquid flow oscillated from the countercurrent

flow to cocurrent flow, while it was restricted by the CCFL curve,

as shown in Fig. 5.15. The calculated pressure drop also

oscillated continuously during the flooding process (Fig. 5.16).

Since the oscillatory behavior are found in most experiments, the

CCFL curve to be specified for the calculation was determined by

the time averaging of the oscillatory values. Thus the calculated

values oscillation should fall on the specified CCFL curve. The

oscillations were resulted from the fact that the 5m5 code select

the lower value of the two flows determined by a normal junction

equation and the other by a CCFL junction equation.

For the 10 cm tube diameter ( case 6 ) , the Wallis type of

flooding curve is not valid any more due to the characteristics of

a large diameter. Thus all calculated flooding data were

presented by the form of Kutateladze form in Fig. 5.17 for the

comparison of the experimental results. The calculation results

show that high oscillations occurred and a complete flooding

occurred at the K* = 4.8 whereas the experimental results show that

g

a complete flooding occurred at K = 3.2. The pressure drop, shown g

in Fig. 5.18, oscillated highly also. When the CCFL option was

specified at the junction, the calculated flooding data lay well

on the specified CCFL curve, although the downward flow oscillated

between countercurrent flow and cocurrent flow as shown in

Fig.5.19 and 20. The curious results were produced above the

complete flooding, K*=3.2. This might have been caused by a

coding error or incorrect implementation of a CCFL model in

43
version 5m5. Therefore no further discussions are provided for

the curious calculated results.

In summary of sensitivity studies using various tube diameters,

the high interfacial drags of the code for the small diameter

tubes can restrict the countercurrent flow before the CCFL

condition is actually reached. As the tube diameter increases,

the calculated downward liquid flow becomes unstable because the

calculated flow regime oscillated between the annular mist flow

and slug flow. The flooding occurs at a much higher air velocity

than predicted by the experiments, and thus the application of the

CCFL option at the specific junction is needed.

5.2 Nodalization

The flooding characteristics may be affected by the

nodalization. The major concern for a nodalization effect was in

the test section area in which the actual countercurrent flow

existed. The figures of the fine nodalization are presented in

Fig. 5.21 and 5.22. The results are shown in Fig. 5.23 to Fig.

33. The major calculated values of fine nodalization are not

quite different from the results of the coarse nodalization as

shown in the section 4, and thus it is concluded that flooding

characteristics are not sensitive to the number of nodalization.

44
RELAP5/MOD3/5M5 without CCFL Model

Test Cone (J feed - 0.1C. D - 0.5 Cm)

0.9 L

0.8

0.7

v 0.6 Cr 0

0.5

o.0.4

S o 0.3 -

0.2

0.1

0

0

Fig. 5.1 Calc'ulated

0.2 0.4 0.6 0.8 1

U FL~..'•,'R- [ ]qrtLf)

Flooding Data without CCFL option ( D = 0.5 cm )

relop5/mod3/5m5 without CCFL Model

Test C-,. (Lhed - 0.16. D - 0.5 cm)

E-.

C

.C

OC

10

8-

7-

6-

5

4

3

2

103

0 -r * * -,-" r• T , 6 I •I I I I

0 20 40 60 80 100 120 140 160 l1o 200

PieDme (e---o)

Fig. Pressure Drop in the Test Section (D 0.5 cm) 5.2

45
RELAP5/MOD3/5M5 using CCFL. Model

... - Cm (-•,."d - O.16..0 - 4.i c •)

I

9~~

t

CU

0

I:

U

0.6

0.4-

Oz

Q3

tic rr

01

0.2

O.*t

U

0 WEMN 8 1

0 02 (L4 0.6 Lsm nlfflt E etu

. . .1

.03

I* . - I

1

Fig. 5.3 Calculated Flooding Data with CCFL option ( D = 0.5 cm )

relop5/rnod3/5m5 using CCFL Model

Twt Com (Lfmd - 0.16, D - 0-5 mn)

10 I -

9-

8-

n.

C

CL

C

'C

7-

6-

5

4

3

2

1- 103

0 . . I

0 20 40 60 80 100 120 140 160 180 200

YWO (8.&)d*)

Fig. 5.4 Pressure Drop in the Test Section ( D = 0.5 cm )

46
•RELAP5/MOD3/5M5 without CCFL Model

Taut C=a (J•,md - 0.16..D - 1 cm)

UD

aU

s-U

U

2

C

Er

S

a

a

0.9

0.8

0.7

0.6

0.5

0.4

0.1

0

Fig. 5.5 Calculated

0.2 0.4 0.6 0.8

Foqdi 1 d lwauta -itO Cf o

Flooding Data w'ito6ut CCFL op~

4

1

tion ( D = 1.0 cm)

relap5/mod3/5m5 without CCFL Model

TGA Come (Lf-• - 0.16. v - 1 cm)

10

9

10

CL

U

C

0

8-

7

6

5-

4-

3-

2-

I

0 * l mlmm

0 20 40 60 80 100 120 140 160 180 200

Tkme (00.W4)

Fig. 5.6 Pressure Drop in the Test Section ( D = 1.0 cm )

47
RELAP5/MOD3/5M5 using CCFL Model

Test Case (Jfeed 0.16. D - I cm)

0~

t

C

0

0

a

0

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.-

0

1

ab

a

0 0.2 0.4 0.6

L.ud F•,o•-te s[rtjf) 1

Fig. 5.7 CalculatediFlooding Data with CCFL

0.8

option ( D = 1.0 cm )

IL

0~

6-0

o -C

Cc

C

0

10

9

8-

7

6

5

4

2

17

relop5/mod3/5m5 using CCFL Model

Test (3m C- e - 0.16, 0 1 -Cm)

I 0 I . " i I i I I I . I o I I I I I i I

0 20 40 60 80 100 120 '40 160 180 200

ý-. (,conds)

Fig. 5.8 Pressure Drop in the Test Section ( D = 1.0 cm )

48
RELAP5/MOD3/5M5 without CCFL Model

Test Cse (J...ed - 0.16. D - 3 am)

0

t

5

a

C

U

C

0

1

0.9 -

0.8 -

0.7 -

0.6 -

0.5 -

0.4 -

0.3 -

0.2 -

0.1 -

0r

10

En I

0

Fig. 5.9 Calculated

6 1 & i I 1 1

0.2 0.4 0.6 0.8

Liqud no-St- E qrt(LO ]

Flooding Data without CCFL option

I

( D = 3.0 cm )

relap5/mod3/5m5 without CCFL Model

Test Come (Lfeed - 0.16. D - 3 cm)

10

9-

8-

0.

C

7

5

4

2

1-

illillillillilllllllllm llllllliiiiiiiiillllIIIIIIIIIIIIIIIIIIIIII1111111111111p ý a %

a .1 --. . . . . i. .. . . .S S S S S S S I • i S

20 40 60 ,80 100 120 140 160 110 200

Time (seconds)

Fig. 5.10 Pressure Drop in the Test Section ( D = 3.0 cm )

49
RELAP5/MOD3/5M5 using CCFL Model

Toat Cc" (Jeed - 0.16. D - 3 am)

0

0

t

S.

.3

C

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.2 0.4 O.6 0.C

LIqIM Fkxuý E *at(L)]I

1

Fig. 5.11 Calculated Flooding

relop5/n

Tee 10

91

Data with CCFL option ( D = 3.0 cm )

nod3/5m5 using CCFL Model

Came (Lfeed - 0.16. V - 3 cm)

0

0~

0.

.0

*0

C

C

a

£

C,

8

7

6-

5

4.

2-

q3

b

A I ..... ý - . I a 4 & 6 1 1 1 1 1 1 . I 1 1 - 1 1 0 20 40 60 o0 100 120 140 160 180 200

Tons (second*)

Fig. 5.12 Pressure Drop in the Test Section ( D = 3.0 cm )

50
1

0.9

0.8

RELAP5/MOD3/5M5 without CCFL Model

Test Ccs (J4jeed - 0.16, D - 5 cm)

03 1

03

al

C.

S

C

a

0

L.

C

0

0.7 "4

0.6 -

0.5

0.4

0.3-

0.2 -

0 -

0 0.2 0.4 0.6

Uquid Fwith o sqtCFL)

9 Diata Without CCFL

I

0.8

Fig. 5.13 Calculated Floodin option ( D = 5.0 cm )

relop5/mod3/5m5 'without CCFL Model

Test Cas Qjeed - 0.16. D - 5 cm)

10 -

9- a

8-

0

0.

0.

0~~

Sc

50

C. g..

0. -~

S

C

C

V

U

7

6-

5-

4

3

2

1-

0

a

a a 0

0 0

] n I_ I

0 -A

0 20 40 60 80 100 120 140 160 180 200

Tme (ewconds)

the Test Section (D =5.0 cm) Fig. 5.14 Pressure Drop in

51
RELAP5/MOD3/5M5 using CCFL Model

Test Cme (Jjfeed - 0.16. D - 5 ci)

0

0*

a

a

2

0

ii:

a

0

0

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.4 0.6

LFqlid lwit h cA C

Fig. 5.15 Calculated Flooding Data with CCFL

0.8

option ( D = 5.0 cm )

relop5/mod3/5m5 using CCFL Model

Twt Cam (Lfeed - 0.16. D - 5 cmi)

C

.0

10

9-

8-

7

6

5

4-

3

2

I

0 0 ~0 1

13

0 .1r 1 4 4 4 4 4 4 4 4 4 4 4 T I I 1 4 4 4 4 r - i i i i i j i I i i I I i I I

0 20 40 60 80 100 120 140 160

Time (,-ond)

i a 0I

ISO 200

Fig. 5.16 Pressure Drop in the Test Section ( D = 5.0 cm )

52
RELAP5/MOD3/5M5 without CCFL Model

Test Case (K toed - 0-966. D - 10 cmi)

01

Cr

0

2.6

2.4

2.2

2

1.8

1.6

1.4

1.2

1

0.8

0.6

0.4

0.2

0

0aEP

a

,P0

1a

n A

a VaRr 0

-FrV a

0 Za

13

0

0

' I

0.2 0.4 0.6

LuquAd no-to [ sqrt(K_f) ]

d Flooding Data without CCFL

0

O.8 i

I

Fig. 5.17 Calculate option ( D = 10. cm )

relop5/mod3/5m5 without CCFL Model

Test Cose (Ij eed - O.96. D - 10 cm)

10 --

9-

0

a0~

o -' '-U0.0

ec

.~0

I.. ~- a.

U

C

C

a

0

8

7

6

5-

4

2

1-

aa

Do

0 13

0 -mmmmCW--

0 40 80 120 160 200 240 280

Time (••ond.)

Fig. 5.18 Pressure Drop in the Test Section ( D 10. cm )

53
RELAP5/MOD3/5M5 using CCFL Model

Test Cc" (K_feed - 0-966. D - 10 cm)

C

a.

2.1

2

1.9

1.8

1.7

1.6

1.5

1.4

1.3

1.2

1.1

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

0 0.2 0.4 0.6 0.8

Uq.uid flowrote [ zqrt(Kj)]

Fig. 5.19 Caalcuiated Floodinig Data with CCFL opti<)n ( D = 10. cm )

:0

i.-

C

AC

0

10

9-

8-

7

6

5

4

3

2-

relap5/mod3/5m5 using CCFL Model

Test Case (jfeed - 0.966. 0 - 10 cm)

0

aa oo

13 0

3 0

0 9

a aaa0

0

0 20 40 60 80 100 120 140 160 180 200

ine (teest)

Fig. 5.20 Pressure Drop in the Test Section ( D =10. cm )

54
Fig. 5.21 RELAP5 Fine Nodalization for Air/Water Test Problem

55
TDV-11 TDV-1 ELEVATION

Water Inlet T-lJ-2 0.00 m

S 02

-- --- 0.25 tn

C: 03

CP 04

C o2 05 ---0.50 m

-P

U S 06

UP--0.75 m - 07

-" • 08

09

10 -- 1.25 n

12

. CCFL Test JUN-4 -- 1.50 m

SLow Plenum I r

SNV-5

Air Inlet TDJ-6

I -- L63 im

Liquid Tank Air Source

SNV-9 TDV-7 --- 133 

Nodalizaoton f or Type-1 Test Bundle

Fig. 5.22 RELAP5 Fine Nodalization for Type-1 Test Bundle

56
RELAP5/MOD3/5M5 using CCFL Model

Fkn Nodrg (JJmed - 0.16. D - 2 cm)

1.1

t

V.

6

e

0

C

0

a

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0

Fig. 5.23

0 0.2 0.4 0.6

uqtdd nFlow [ .q-ry_) 3

Calculated Flooding Data with Fine

Ox.

Nodalization

relop5/mod3/5m5 using CCFL Model

Fln NOdhn (Lfe* - 0.16. D - 2 cm)

10 ~1

9-

V

0.

0.

go

C

C

0

0

8-

7 -

6

5

4-

3

2-

1-

13

2

0 -IlP

0

. . . . . 4 .. a a a a a I r r 20 40 60 80 100 120 140 160 180 200

Fim. wi (5D nea)

Fig. 5.24 Pressure Drop with Fine Nodalization

57
RELAP5/MOD3/5M5 using CCFL Model

Fin. Noding (J4_fe - 0.16. D - 2 CM)

E

E

SY

0

a

7.5

7

6.5

6

5.5

5

4.5

4

0 20 40 60 80 100 120 140 160 180 200

rime ( Sands )

Fig. 5.25 Flow Regime Change at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 using CCFL Model

Fin. Nodn• (J food - 0.16. D - 2 ern)

E

94

2

1.5

50

40

30

20

10

0

0 20 40 60 80 100 120 140

Tvr,o C sec"Wds )

160 180 200

Fig. 5.26 Interfacial Drag Coefficient at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

58
RELAP5/MOD3/5M5 using CCFL Model

Fin oding (Jed - 0.16. D - 2 cm)

0.9

0

0

40

0.7

0.6

0.5

0 20 40 60 80 100 120 140 160 180 200

Fig. 5.27 Gas Void Fraction at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 6 - jun.4000000)

relap5/mod3/5m5 using CCFL Model

FIe NodkV (L.fed - 0.18. D - 2 cm)

CIO

E

"I

CD

X

100

90

80

70

60

50

40

30

20

10

0

-10

-20

-30

-40

-50

I.

0 20 40 60 80 100 120

irme (*econd,)

140 10 180 200

Fig. 5.28 Mass Flux through the Junction 4000000

59
RELAP5/MOD3/5M5 using CCFL Model

Fin Hoing ( QI-2600 cc/min )

t

a.

U

6

2

0

U

a

0

1.5

1.4

1.3

1.2

1.1

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0-

00.2

Fig. 5.29 Calculated

04. 0.6 0.8

Uquod inoDtuh [ e oqrt(Kjf

Flooding Data with Fine Noidalization

relOp5/mod3/5m5 using CCFL Model

F6ne Nodng ( 4-2600 cc/mn )

5

4 -

0.

C

60

z

U

3 -

2

1 -

-I-I

a -I 0 I I . I I , , J

0 0.2 . 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

Air Flow [ KSg

Fig. 5.30 Pressure Drop with Fine Nodalization

60
RELAP5/MOD3/5M5 using CCFL Model

Rod Bundle Typel ( CH-2600cc/in )

E

0

AX.

a

7.5

7

6.5

6

5.5

5

4.5

4 SI

0 20 40 60 80 100 120 140 160 180 200

Time ( Sow%)

Fig. 5.31 Flow Regime Change at the Test Section Junction

(1 - jun.3010000 : 2 - jun.3020000 : 3 - jun.3030000

4 - jun.3040000 : 5 - jun.3050000 : 6 - jun.4000000)

RELAP5/MOD3/5M5 using CCFL Model

Rod Bundle Typel ( 4-2650 cc/mln )

0 Z 0

LC

.S

0.05

0.04

0.03

0.02

0.01

0

F ig. 5.32

0 20 40 60 80 100 120 140 160 ISO 200

Time ( Secon )

Interfacial Drag Coefficient at the Test Section Junctlcn•

(1 - jun.3010000 2 - jun.3020000 3 - jun.3030000

4 - jun.3040000 5 - jun.3050000 6 - jun.4000000i

61
RELAP5/MOD3/5M5 using CCFL Model

Rod Bundie Typal ( QI-2600 cc/mli )

C. 0

0o

0 0.7

0.6

0.5

0 20 40 60 80 100 120 140 160 180 200

77me ( Seconds )

Fig. 5.33 Gas Void Fraction at the Test Section Junction

(1 - jun.3010000 2 - jun.3020000 3 - jun.3030000

4 - jun.3040000 5 - jun.3050000 6 - jun.4000000)

62
6. Run Statistics

The simulation for the CCFL model assessment is not a transient

case basically, thus it is meaningless to compare the simulation

time. However the range of air flow used to test the CCFL model

was obtained by increasing air flow during a given period.

Therefore the run time required by the simulation will be helpful

for the other assessment or determination of the transient period.

The computer used for the simulation is a CRAY-2S which has 4

central process units, and the operating system is UNICOS. The

time step size was fixed as 0.01 second to avoid the numerical

instability of the liquid film flow, discussed in the section 4.1.

The whole transient time was also fixed as 200 seconds. The

computation efficiency is summarized in Table 6.1. A plot of CPU

time was not needed because it increases linearly with the

transient time.

Table 6.1 CPU time required in the simulation

Case Transient Computing Number of Number of Total Grind

Real Time Time Total Volume Time Step Time

(sec) (CPU) (C) (DT) (CPU/C/DT)

1 200.0 151.52 11 20,000 6.89 x 10-'

2 200.0 152.34 11 20,000 6.92 x 10 4

3 200.0 156.41 11 20,000 7.11 x 10-4

4 200.0 159.83 11 20,000 7.27 x 10 4

5 200.0 263.14 17 20,000 7.74 x 10-4

" case 1 : coarse noding without CCFL option in simple tube

" case 2 : coarse noding with CCFL option in simple tube

" case 3 : coarse noding without CCFL option in rod bundle

* case 4 : coarse noding with CCFL option in rod bundle

" case 5 : fine noding with CCFL option in simple tube

63
7. Conclusion

Assessment calculations for sample problems in a simple tube

and a rod bundle were carried out using RELAP5/MOD3/5m5 code. For

the purpose of comparison, a pair of calculations were performed

with and without CCFL option. Sensitivity studies with respect to

the diameter of tube was also performed. In the small diameter

case, the flooding was calculated to occur before the CCF.L

occurred in the experiment, due to the increased interfacial drag

force. In the large diameter case, the flooding occurs at a

higher gas velocity than the experimental results, when the CCFL

model option was not selected. With applying the CCFL option to

the specific junction, the calculated values lay well on the

curve specified by user, and the pressure oscillations and jump

after onset of flooding agreed with most experimental

observations. Curious calculation results obtained by specifying

the Kutateladze CCFL curve for the large diameter tube may be

resulted from a coding error of 5m5 version.

Other assessment calculations were carried out for a rod bundle

test, which was performed by KINS( Korea Institute of Nuclear

Safety ). Calculated results were compared with experimental data

for each type of test section ; type-1 which has no spacer grid,

type-2 which has one spacer grid, and type-3 which has two spacer

grids. A flooding curve which was determined by experiments was

applied to the specific junction. Calculated values lay well on

the curve specified by user, and the pressure drop increased

suddenly just after onset of flooding due to restriction of

downward liquid flow at the specific junction. However, the

pressure jump before onset of flooding was not calculated because .

RELAP5 did not have models of "liquid mixing" which was observed

in the experiment at the transition area before "onset of

flooding".

The fine nodalization for the test section did not affect the

64
whole flooding phenomena, and it was concluded that the

interfacial drag model and CCFL curve given by user mainly

characterized the flooding process.

65
References

1. V.H. Ransom, et al.,"RELAP5/MOD2 Code Manual Volume 1: Code

Structure, System Models, and Solution Methods",NUREG/CR-4312 (1985)

2. V.H. Ransom, et al.,"RELAP5/MOD1 Code Manual Volume 1: Code

Structure, System Models, and Solution Methods",NUREG/CR-1826 (1982)

3. ICAP Relap5 Newsletter , NUREG/BR-0103,Vol.3, No.4, October (1987)

4. R.A. Riemke, "Report on CCFL and Interphase Drag Model for

RELAP5/MOD3", EGG-TFM-8012, February (1988)

5. J.M. Putney, "Development of a New Bubbly-Slug Interfacial Friction

Model for RELAP5 - Final Report", CEGB Draft Report,

PWR/HTWG/P(89)722, June (1989)

6. G.B. Wallis and S. Makkenchery, "The Hanging Film Phenomena in

Vertical Annular Two Phase Flows", Journal of Fluid Engg.,Sep., pp

297-298 (1974)

7. M.H.Chun and B.D.Chung, "Countercurrent Flow Limitation Experiments

with Vertical Multi-tube Bundle with and without Spacer Grids", Int.

Comm. Heat Mass Transfer, Vol.17, pp.237-246 (1990)

8. M.H.Chun et. al., "Experimental Investigation of Flooding in

Air-Water Countercurrent Flow with a vertical adiabatic Multi-rod

Bundle", FED-Vol.99, HTD-Vol.155, pp.237-242,Winter Annual Meeting,

ASME, Nov. (1990)

9. S.G. Bankoff, R.S. Tankin, M.C. Yuen and C.L. Hsieh,"Countercurrent

Flow of Air/Water and Steam/Water through a Horizontal Perforated

Plate", Int. J. Heat & Mass Transfer, Vol.24, No.8, pp 1381-1395

(1981)

10. D. Bharathan, G.B. Wallis and H.J. Richter,"Air-water

Countercurrent Annular Flow", EPRI NP-1165, Sep. (1979)

66
APPENDIX A. EXPERIMENTAL DATA FOR ROD BUNDLE CCFL TEST

Appendix A. Experimental Data for Rod Bundle CCFL Test

Numerical data for experiments of the three type of test

section are tabulated as following. In the following tables the

oscillation data are presented as averaged value ± 1 standard

deviation. The nondimensional conversions are performed by using

the physical and geometrical conditions of the present experiment.

The conditions used for nondimensional conversions are listed as

following.

TI= Temperature of Liquid = 282 °K

T = Temperature of Gas = 280 *K

pf= Density of Liquid = 1000 Kg/m 3

Pg = Density of Gas = 1.2713 Kg/mi3

pJf= Dynamic Viscosity of Liquid = 1303 X 10- Kg/m.sec

Pg Dynamic Viscosity of Gas = 173.9 X 10- 6Kg/m.sec

a = Surface Tension of Liquid = 73.656 X 10- 3 N/m

g = Universal acceleration Constant = 9.81 m/sec2

DH= Equivalent Hydraulic Diameter = 4A/P = 1.03 X 10-2 m

The following lists of the nomenclatures are used in tables.

0 = Water Injection Rate

K= Kutateladze Number of Gas Flux = p 05xg V /[ga(pf-pg )]0.25

K = Kutateladze Number of Liquid = Pf V f/ga(pfp 0 25

K= Kutateladze Number of Water Injection

j = Nondimensional Gas Flux = pgcagVg/[gDH(Pf-PgJ)]0.

jf = Nondimensional Liquid Flux = pfaf /[gDH(pf-p•]° 5

JI = Nondimensional Liquid Injection

P Nondimensional Pressure Gradient = IdP/dzl/[g(pf-pg)]

A.1
Re£ = Film Reynolds Number = 4r/p/ ., where r' denotes film mass flowrate

per unit width

V= Critical Nondimensional Gas Velocity

!AfV g,c 0.s 0.8 =-- (P /pf) (11N/•

where N P= Viscosity Number = pf./(PfJ'7C/(Pf-Pg)] /2

V g,c = Gas Velocity which results flooding

A.2
Test Data for

QI = 0 cc/min

Single Phase

T I= 282 *K

Air Flow in Bare Rod Bundle

T= 280 °K

9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

0.252a

5.150±

5.799±

6.825±

8.207±

9.134±

9.718±

11.240±

11.732±

0.257

0.061

0.053

0.023

0.066

0.053

0.046

0.053

0.063

0.224

4.573

5.149

6.061

7.288

8.111

8.629

9.981

10.418

1.3

16.4

20.8

25.2

34.0

38.4

43.8

77.7

79.2

±

±

±

±

±

±

±

±

±

3.7

0.0

2.8

1.2

0.3

0.0

0.5

0.5

1.3

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

Nondimensional Form of Test Data

K*= 0.0 j 0.0 Reý - 0.0 V* = 0.0 C

K P ~ K * K *i*i

9f c9 f4

1

2

3

4

5

6

7

8

9

0.049±

0.995±

1.120±

1.318±

1.586±

1.765±

1.877±

2.171:±

2.266±

0.050

0.012

0.010

0.005

0.013

0.010

0.009

0.010

0.012

0.001±

0.01 6±

0.021±

0.025±

0.034±

0.038±

0.044±

0.078±

0.079±

0.004

0.000

0.003

0.001

0.000

0.000

0.000

0.001

0.001

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.025±

0.513±

0.578±

0.680±

0.818±

0.911±

0.969±

1.120±

1.169±

0.026

0.006

0.005

0.002

0.007

0.005

0.005

0.005

0.006

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

A.3
Test Data for Two Phase Air Flow in Bare Rod Bundle

Q1= 230 cc/min T I= 282 °K T= 280 *K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

0.000±

6.764±

7.853±

9.055±

9.456±

9.844±

10.736±

11.248±

8.823±

8.448±

8.083±

0.000

0.035

0.016

0.069

0.044

0.054

0.026

0.073

0.012

0.066

0.039

0.000

6.007

6.974

8.041

8.397

8.741

9.533

9.988

7.834

7.502

7.178

26.7

28.4

33.9

82.7

94.7

110.3

80.5

82.3

91.9

99.1

111.9

±

±

±

±

±

±

±

±

±

±

±

2.3

6.0

2.7

4.1

6.8

6.9

3.5

1.9

7.7

7.3

9.2

227.5

227.5

227.5

85.8

62.2

38.6

23.6

0.0

79.4

126.6

156.7

0.0

0.0

0.0

141.6

165.3

188.9

203.9

227.5

148.1

100.9

70.8

Nondimensional Form of Test Data

K= 0.021 0.011 Re. = 43.7 V = 1.323

C

A A

K K9 9 P K*• f K Kc JgJf . .9Jf J.c

1

2

3

4

5

6

7

8

9

10

11

0.000±

1.307±

1.517±

1.749±

1.827±

1.902±

2.074±

2.173±

1.704±

1.632±

1.562±

0.000

0.007

0.003

0.013

0.009

0.010

0.005

0.014

0.002

0.013

0.008

0.027±

0.028±

0.034±

0.083±

0.095±

0.110±

0.081±

0.082±

0.092±

0.099±

0.112±

0.002

0.006

0.003

0.004

0.007

0.007

0.004

0.002

0.008

0.007

0.009

0.021

0.021

0.021

0.008

0.006

0.003

0.002

0.000

0.007

0.011

0.014

0.000

0.000

0.000

0.013

0.015

0.017

0.018

0.021

0.013

0.009

0.006

0.000±

0.674±

0.783±

0.903±

0.943±

0.981±

1.070±

1.121±

0.879±

0.842±

0.806±

0.000

0.004

0.002

0.007

0.004

0.005

0.003

0.007

0.001

0.007

0.004

0.011

0.011

0.011

0.004

0.003

0.002

0.001

0.000

0.004

0.006

0.007

0.000

0.000

0.000

0.007

0.008

0.009

0.009

0.011

0.007

0.005

0.003

A.4
Test Data for Two Phase Air Flow in Bare Rod Bundle

Q 530 cc/min T =282 °K Tg= 280 °K

g

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1 0.000± 0.000 0.000 25.6 ± 2.2 530.1 0.0

2 6.195± 0.023 5.501 55.5 ± 0.0 530.1 0.0

3 7.220± 0.072 6.411 123.3 ± 12.2 530.1 0.0

4 7.453± 0.045 6.618 130.0 ± 11.7 354.1 149.3

5 7.818± 0.045 6.942 127.7 ± 10.7 255.4 247.7

6 8.058± 0.051 7.155 122.7 ± 7.9 176.0 327.5

7 8.332± 0.020 7.399 119.2 ± 7.4 137.4 365.5

8 10.041± 0.034 8.916 87.3 ± 3.8 57.9 450.7

9 11.710± 0.072 10.399 85.9 ± 4.9 0.0 530.5

Nondimensional Form of Test Data

K= 0.048 jI= 0.025 Re£ = 100.7 V = 1.235

C

K* K K* .* *. . 9 Kf c g Jf Jc

1

2

3

4

5

6

7

8

9

0.000±

1.197±

1.395±

1.440±

1.510±

1.557±

1.610±

1.940±

2.262±

0.000

0.004

0.014

0.009

0.009

0.010

0.004

0.007

0.014

0.026±

0.056±

0.123±

0.130±

0.128±

0. 123±

0.119±

0.087±

0.086±

0.002

0.000

0.012

0.012

0.011

0.008

0.007

0.004

0.005

0.048

0.048

0.048

0.032

0.023

0.016

0.012

0.005

0.000

0.027

0.027

0.027

0.040

0.042

0.044

0.046

0.041

0.036

0.000±

0.618±

0.720±

0.743±

0.779±

0.803±

0.831 ±

1.001±

1.167±

0.000

0.002

0.007

0.005

0.004

0.005

0.002

0.003

0.007

0.025

0.025

0.025

0.016

0.012

0.008

0.006

0.003

0.000

0.000

0.000

0.000

0.007

0.012

0.015

0.017

0.021

0.025

A.5
Test Data for Two Phase Air Flow in Bare Rod Bundle

Q= 730 cc/min T =282 °K T-= 280 °K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

0.535±

6.228±

6.829±

7.287±

7.677±

8.045±

8.299±

8.472±

8.946±

9.683±

10.375±

10.807±

11.421±

8.396±

7.859±

7.176±

7.336±

7.466±

7.704±

0.134

0.024

0.023

0.060

0.064

0.047

0.027

0.068

0.046

0.038

0.030

0.042

0.043

0.047

0.020

0.080

0.020

0.048

0.074

0.475

5.530

6.065

6.471

6.817

7.144

7.370

7.523

7.944

8.599

9.213

9.597

10.142

7.456

6.979

6.372

6.515

6.630

6.841

38.4

57.3

139.3

135.5

134.7

130.7

127.5

124.7

116.4

101.3

85.4

80.5

84.7

123.2

133.5

139.0

136.4

134.1

133.7

±

±

±

±

±

±

±

±

±

±

±

±

±

±

±

±

+

±

±

0.6

1.6

11.5

12.6

11.6

11.2

8.3

6.8

4.0

2.4

4.0

2.2

2.7

6.3

10.2

11.4

11.4

11.4

10.5

727.6

727.6

706.1

530.1

321.9

216.8

176.0

133.1

118.0

57.9

0.0

0.0

0.0

137.4

236.1

510.8

491.5

431.4

313.3

-0.0

-0.0

21.5

197.5

405.6

510.8

551.6

594.5

609.5

669.6

727.6

727.6

727.6

590.2

491.5

216.8

236.1

296.2

414.2

A.6
Nondimensional Form of Test Data

K*0.066 j 0.034 Re = 138.7 V= 1.147

K * p • K* K: *: j*.•* .

9 f J9 if Jc

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

0.103± 0.026

1.203± 0.005

1.319± 0.005

1.408± 0.012

1.483± 0.012

1.554± 0.009

1.603± 0.005

1.637± 0.013

1.728± 0.009

1.871± 0.007

2.004± 0.006

2.088± 0.008

2.206± 0.008

1.622± 0.009

1.518± 0.004

1.386± 0.016

1.417± 0.004

1.442± 0.009

1.488± 0.014

0.038± 0.001

0.057± 0.002

0.139± 0.012

0.136± 0.013

0.135± 0.012

0.131± 0.011

0.128± 0.008

0.125± 0.007

0.117± 0.004

0.101± 0.002

0.086± 0.004

0.081± 0.002

0.085± 0.003

0.123± 0.006

0.134± 0.010

0.139± 0.011

0.137± 0.011

0.134± 0.011

0.134± 0.011

0.066 -0.000

0.066 -0.000

0.064

0.048

0.029

0.020

0.016

0.012

0.011

0.005

0.000

0.000

0.000

0.012

0.021

0.046

0.044

0.039

0.028

0.002

0.018

0.037

0.046

0.050

0.054

0.055

0.060

0.066

0.066

0.066

0.053

0.044

0.020

0.021

0.027

0.037

0.053± 0.013-

0.621± 0.002

0.681± 0.002

0.726± 0.006

0.765± 0.006

0.802± 0.005

0.827± 0.003

0.844± 0.007

0.892± 0.005

0.965± 0.004

1.034± 0.003

1.077± 0.004

1.138± 0.004

0.837± 0.005

0.783± 0.002

0.715± 0.008

0.731± 0.002

0.744± 0.005

0.768± 0.007

0.034 -0.000

0.034 -0.000

0.033 0.001

0.025 0.009

0.015 0.019

0.010 0.024

0.008 0.026

0.006 0.028

0.005 0.028

0.003 0.031

0.000 0.034

0.000 0.034

0.000 0.034

0.006 0.027

0.011 0.023

0.024 0.010

0.023 0.011

0.020 0.014

0.015 0.019

A.7
Test Data for Two Phase Air Flow in Bare Rod Bundle

QI 1570 cc/min Ti= 282 °K T = 280 °K

9

Air

Flowrate

(1000 cc/sec)

Channel Pressure Liquid Flowrate

Velocity

(m/sec)

Drop

(mm-H20)

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

0.052±

4.711±

5.715±

6.2144±

6.454±

6.563±

6.808±

6.892±

7.206±

7.613±

8.066±

8.290±

8.535±

9.014±

10.065±

10.657±

0.150

0.055

0.058

0.027

0.063

0.043

0.023

0.045

0.076

0.029

0.045

0.032

0.042

0.049

0.044

0.058

0.046

4.183

5.075

5.518

5.631

5.827

6.046

6.120

6.399

6.760

7.163

7.361

7.579

8.005

8.937

9.463

13.3

45.6

135.7

161.1

161.6

158.3

158.3

154.9

153.8

148.9

140.3

137.6

132.8

123.1

99.5

85.1

±

±

±

±

±

±

±

±

±

±

±

±

±

±

±

±

0.7

7.5

15.7

17.6

18.4

14.6

13.5

14.8

14.1

12.0

11.8

11.2

9.5

4.6

5.0

3.1

1571.0

1571.0

1571.0

1296.3

1098.9

1021.6

903.5

824.1

652.4

354.1

156.7

118.0

79.4

62.2

38.6

0.0

-0.0

-0.0

0.0

274.7

472.2

549.4

667.5

746.9

918.6

1216.9

1414.3

1453.0

1491.6

1508.8

1532.4

1571.0

Nondimensional Form of Test Data

K= 0.142 j 0.073 Rel = 298.3 V* = 1.015 C

K K * K K* .* .

gf c Jg f Jc

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

0.010±

0.910±

1.104±

1.200±

1.247±

1.268±

1.31 5±

1.331±

1.392±

1.471±

1.558±

1.601±

1.649±

1.741±

1.944±

2.059`+

0.029

0.011

0.011

0.005

0.012

0.008

0.004

0.009

0.015

0.006

0.009

0.006

0.008

0.009

0.009

0.011

0.01 3t

0.046±t

0. 136±t

0.161:±

0. 162±

0.158±

0.158±

0.155±

0.154±

0.149±

0.140±

0.138±

0.133±

0. 123±

0.100±

0.085±

0.001

0.007

0.016

0.018

0.018

0.015

0.014

0.015

0.014

0.012

0.012

0.011

0.010

0.005

0.005

0.003

0.142

0.142

0.142

0.117

0.099

0.092

0.082

0.074

0.059

0.032

0.014

0.011

0.007

0.006

0.003

0.000

-0.000

-0.000

-0.000

0.025

0.043

0.050

0.060

0.067

0.083

0.110

0.128

0.131

0.135

0.136

0.138

0.142

0.005±

0.470±

0.570±

0.619±

0.643±

0.654±

0.679±

0.687±

0.718±

0.759±

0.804±

0.826±

0.851±

0.899±

1.003±

1.062±

0.015

0.005

0.006

0.003

0.006

0.004

0.002

0.004

0.008

0.003

0.005

0.003

0.004

0.005

0.004

0.006

0.073

0.073

0.073

0.060

0.051

0.048

0.042

0.038

0.030

0.016

0.007

0.005

0.004

0.003

0.002

0.000

-0.000

-0.000

-0.000

0.013

0.022

0.026

0.031

0.035

0.043

0.057

0.066

0.068

0.069

0.070

0.071

0.073

A.8
Test Data for Two Phase Air Flow in Bare Rod

Q 2600 cc/min TI= 282 'K T = 280 °K

Bundle

Air

Flowrate

(1000 cc/sec)

Channel Pressure Liquid Flowrate

Velocity

(m/sec)

Drop

(mm-H20)

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

4.509±

4.898±

5.194±

5.517±

5.790±

5.995±

6.377±

6.832±

7.441±

8.136±

8.387±

9.114±

10.056±

0.120

0.065

0.059

0.063

0.081

0.103

0.080

0.027

0.045

0.061

0.042

0.061

0.040

4.004

4.349

4.612

4.900

5.142

5.324

5.663

6.066

6.608

7.225

7.448

8.093

8.929

55.8

187.7

210.4

205.5

200.4

197.9

189.8

172.9

159.3

146.2

140.5

125.5

108.8

±

±

±

±

±

±

±

±

±

±

±

±

±

1.5

30.1

34.6

31.9

32.1

29.9

25.5

18.0

14.3

12.5

12.1

7.5

7.1

2590.5

2590.5

2120.4

1807.1

1571.0

1491.6

1216.9

942.2

588.1

156.7

79.4

62.2

0.0

-0.0

-0.0

470.0

783.4

1019.4

1098.9

1373.6

1648.3

2002.4

2433.8

2511.1

2528.2

2590.5

Nondimensional Form of Test Data

K> 0.234 j 0.121 ReT = 494.0 V* = 0.794

C

K * K K .* .* .I

9 f c 9 f c

1

2

3

4

5

6

7

8

9

10

11

12

13

0.871±

0.946±

1.003±

1.066±

1.119±

1.158±

1.232±

1.320±

1.438±

1.572±

1.620±

1.761±

1.943±

0.023

0.013

0.011

0.012

0.016

0.020

0.015

0.005

0.009

0.012

0.008

0.012

0.008

0.056±

0.188±

0.211±

0.206±

0.201±

0.198±

0.190±

0.173±

0.159±

0.146±

0.141±

0.126±

0.109±

0.001

0.030

0.035

0.032

0.032

0.030

0.026

0.018

0.014

0.012

0.012

0.008

0.007

0.234

0.234

0.191

0.163

0.142

0.135

0.110

0.085

0.053

0.014

0.007

0.006

0.000

-0.000

-0.000

0.042

0.071

0.092

0.099

0.124

0.149

0.181

0.220

0.227

0.228

0.234

0.449±

0.488±

0.518±

0.550±

0.577±

0.598±

0.636±

0.681±

0.742±

0.811±

0.836±

0.908±

1.002±

0.012

0.006

0.006

0.006

0.008

0.010

0.008

0.003

0.005

0.006

0.004

0.006

0.004

0.121

0.121

0.099

0.084

0.073

0.069

0.057

0.044

0.027

0.007

0.004

0.003

0.000

-0.000

-0.000

0.022

0.036

0.047

0.051

0.064

0.077

0.093

0.113

0.117

0.118

0.121

A.9
Test Data for Two Phase Air Flow in Bare Rod Bundle

Q =3600 cc/min T 282 °K T-= 280 °K

g

Air

Flow rate

(1000 cc/sec)

Channel Pressure Liquid Flowrate

Velocity

(m/sec)

Drop

(mm-H20)

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

4.137±

4.580±

5.110±

5.355±

5.900±

6.184±

6.510±

6.972±

7.722±

8.172±

8.584±

9.521±

0.118

0.155

0.092

0.124

0.076

0.036

0.054

0.062

0.083

0.072

0.040

0.072

3.673

4.067

4.538

4.756

5.239

5.491

5.781

6.191

6.857

7.257

7.623

8.455

39.6

194.3

213.5

218.1

208.1

201.2

188.9

176.2

156.5

151.4

143.5

118.6

±

±

±

±

±

±

±

±

±

±

±

±

8.1

37.2

32.7

29.2

25.1

24.2

21.5

15.3

9.8

10.3

10.5

6.7

3612.1

3219.3

2199.9

1901.5

1455.1

1373.6

1139.6

942.2

510.8

236.1

79.4

0.0

-0.0

392.8

1412.2

1710.5

2156.9

2238.5

2472.4

2669.9

3101.3

3376.0

3532.6

3612.1

Nondimensional Form of Test Data

K 0.326 j 0.128 Re. = 684.0 V = 0.706 C

K* P= K* K * g .f *c

9 f c J f c

1

2

3

4

5

6

7

8

9

10

11

12

0.799±

0.885±

0.987±

1.035±

1.140±

1.195±

1.258±

1.347±

1.492±

1.579±

1.658±

1.839±

0.023

0.030

0.018

0.024

0.015

0.007

0.011

0.012

0.016

0.014

0.008

0.014

0.040±

0.195±

0.214±

0.218±

0.208±

0.201±

0.189±

0.176±

0.157±

0.152±

0.144±

0.119±

0.008

0.037

0.033

0.029

0.025

0.024

0.022

0.015

0.010

0.010

0.010

0.007

0.326

0.291

0.199

0.172

0.131

0.124

0.103

0.085

0.046

0.021

0.007

0.000

-0.000

0.035

0.128

0.154

0.195

0.202

0.223

0.241

0.280

0.305

0.319

0.326

0.412±

0.457±

0.509±

0.534±

0.588±

0.616±

0.649±

0.695±

0.770±

0.815±

0.856±

0.949±

0.012

0.015

0.009

0.012

0.008

0.004

0.005

0.006

0.008

0.007

0.004

0.007

0.168

0.150

0.102

0.089

0.068

0.064

0.053

0.044

0.024

0.011

0.004

0.000

-0.000

0.018

0.066

0.080

0.100

0.104

0.115

0.124

0.144

0.157

0.165

0.168

A.10
Test Data for Single Phase Air Flow

Qi 0 cc/min T= 282 0K Tg=

9

in Single Grid Rod Bundle

280 * K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

3.243±

4.597±

5.524±

6.316±

7.295±

7.842±

8.569±

8.841±

10.134±

10.757±

0.102

0.109

0.024

0.058

0.034

0.020

0.023

0.028

0.055

0.035

2.880

4.082

4.905

5.608

6.478

6.963

7.609

7.851

8.999

9.552

15.4

20.8

24.7

29.2

34.2

37.2

42.3

52.6

74.5

88.8

±

±

±

±

±

±

±

±

±

±

1.3

3.2

2.7

3.2

1.2

1.9

3.3

5.3

3.9

5.4

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

Nondimensional Form of Test Data

K = 0.0 0.0 0.0 V * = 0.0

C

K * pK * K .* .* * K Kf Jg Jf Jc

1

2

3

4

5

6

7

8

9

10

0.626±

0.888±

1.067±

1.220±

1.409±

1.515±

1.655±

1.708±

1.958±

2.078±

0.020

0.021

0.005

0.011

0.007

0.004

0.004

0.005

0.011

0.007

0.015±

0.021±

0.025±

0.029±

0.034±

0.037±

0.042±

0.053±

0.075±

0.089±

0.001

0.003

0.003

0.003

0.001

0.002

0.003

0.005

0.004

0.005

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.323±

0.458±

0.551±

0.630±

0.727±

0.782±

0.854±

0.881±

1.01 0±

1.072±

0.010

0.011

0.002

0.006

0.003

0.002

0.002

0.003

0.005

0.003

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

A. 11
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

Q0= 180 cc/min T =282 °K T_ z280 °K

1 g

Air

Flowrate

(1000 cc/sec)

Channel Pressure

Velocity Drop

(m/sec) (mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

3.434±

5.655±

7.843±

7.980±

7.978±

8.027±

8.277±

8.8199±

8.977±

9.374±

10.064±

0.033

0.045

0.020

0.051

0.059

0.050

0.036

0.015

0.041

0.042

0.042

3.050

5.022

6.965

7.086

7.084

7.128

7.350

7.831

7.972

8.325

8.937

42.7

56.8

124.4

128.5

130.2

127.9

129.4

120.9

107.2

76.1

54.7

±

±

±

±

±

±

±

±

±

±

±

3.0

0.8

3.6

5.4

6.6

4.9

6.6

5.1

5.6

3.0

3.0

180.6

180.6

157.1

117.8

94.2

110.0

62.8

15.7

4.7

0.0

0.0

0.0

0.0

23.6

62.8

86.4

70.7

117.8

164.9

175.9

180.6

180.6

Nondimensional Form of Test Data

Ki:*= 0.016 4=0.008 ReT, = 34.2 V* = 1.323 C

K P Kf K c J Jf Jc

1.

2

3

4

5

6

7

8

9

10

11

0.663±

1.092±

1.515±

1.542±

1.541±

1.551±

1.599±

1.704±

1.734±

1.811±

1.944±

0.006

0.009

0.004

0.010

0.011

0.010

0.007

0.003

0.008

0.008

0.008

0.043±

0.057±

0.125±

0.129±

0.130±

0.128±

0.130±

0.121±

0.107±

0.076±

0.055±

0.003

0.001

0.004

0.005

0.007

0.005

0.007

0.005

0.006

0.003

0.003

0.016

0.016

0.014

0.011

0.009

0.010

0.006

0.001

0.000

0.000

0.000

0.000

0.000

0.002

0.006

0.008

0.006

0.011

0.015

0.016

0.016

0.016

0.342±

0.564±

0.782±

0.795±

0.795±

0.800±

0.825±

0.879±

0.895±

0.934±

1.003±

0.003

0.005

0.002

0.005

0.006

0.005

0.004

0.001

0.004

0.004

0.004

0.008

0.008

0.007

0.005

0.004

0.005

0.003

0.001

0.000

0.000

0.000

0.000

0.000

0.001

0.003

0.004

0.003

0.005

0.008

0.008

0.008

0.008

A.12
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

I = 530 cc/min T = 282 °K T = 280 °K

I 9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

2.554±

5.132±

6.977±

7.318±

7.445±

7.865±

7.954±

8.282±

8.361±

8.819±

9.228±

0.078

0.067

0.051

0.027

0.047

0.022

0.047

0.030

0.016

0.007

0.023

2.268

4.557

6.195

6.498

6.611

6.984

7.063

7.354

7.424

7.831

8.195

53.3 ±

71.4 ±

109.7 ±

104.0±

104.8 ±

104.6 ±

88.5 ±

83.0 ±

90.1±

88.1 ±

76.1 ±

4.6

4.4

8.2

11.1

8.2

8.0

7.4

6.3

12.8

8.6

3.4

530.1

530.1

432.0

353.4

255.2

98.2

58.9

39.3

19.6

3.9

0.0

0.0

0.0

98.2

149.2

255.2

412.3

471.2

490.9

510.5

530.0

530.1

Nondimensional Form of Test Data

K/ 0.048 j 0.025 Re£ = 100.7 V = 1.191 C

K K 9* Kf ,K f K ., .* ,: KC Jg Jf Jc

1

2

3

4

5

6

7

8

9

10

11

0.493±

0.991±

1.348±

1.41 4±

1.438±

1.519±

1.537±

1.600±

1.615±

1.704±

1.783±

0.015

0.013

0.010

0.005

0.009

0.004

0.009

0.006

0.003

0.001

0.005

0.053±

0.072±

0.110±

0.104±

0.105±

0.105±

0.089±

0.083±

0.090±

0.088±

0.076±

0.005

0.004

0.008

0.011

0.008

0.008

0.007

0.006

0.013

0.009

0.003

0.048

0.048

0.039

0.032

0.023

0.009

0.005

0.004

0.002

0.000

0.000

0.000

0.000

0.009

0.013

0.023

0.037

0.043

0.044

0.046

0.048

0.048

0.255±

0.512±

0.695±

0.729±

0.742±

0.784±

0.793±

0.826±

0.833±

0.879±

0.920±

0.008

0.007

0.005

0.003

0.005

0.002

0.005

0.003

0.002

0.001

0.002

0.025

0.025

0.020

0.016

0.012

0.005

0.003

0.002

0.001

0.000

0.000

0.000

0.000

0.005

0.007

0.012

0.019

0.022

0.023

0.024

0.025

0.025

A.13
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

0I= 750 cc/min Ti= 282 °K T_= 280 OK

9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

2.849±

5.495±

6.249±

6.895±

7.052±

7.404±

7.792±

8.013±

8.312±

8.630±

9.154±

0.008

0.058

0.011

0.045

0.047

0.036

0.066

0.052

0.022

0.069

0.048

2.530

4.879

5.549

6.122

6.262

6.574

6.919

7.116

7.381

7.664

8.129

46.1 ±

53.7±

74.1 ±

107.0 ±

108.3 ±

107.1 ±

107.6 ±

105.2 ±

97.1 ±

80.8 ±

71.8 ±

7.6

14.8

5.1

9.8

8.4

9.0

7.2

6.7

8.2

2.2

5.9

746.1

746.1

742.2

490.9

451.6

294.5

106.0

58.9

19.6

3.9

0.0

0.0

0.0

3.9

255.2

294.5

451.6

628.3

687.2

726.5

742.*2

746.1

Nondimensional Form of Test Data

K= 0.067 j 0.035 Rer = 142.5 V = 1.059 C

K p .K* .* .•

f C J f c

1

2

3

4

5

6

7

8

9

10

11

0.550±

1.061±

1.207±

1.332±

1.362±

1.430±

1.505±

1.548±

1.606±

1.667±

1.768±

0.002

0.011

0.002

0.009

0.009

0.007

0.013

0.010

0.004

0.013

0.009

0.046±

0.054±

0.074±

0.107±

0.108±

0.107±

0.108±

0.105±

0.097±

0.081±

0.072±

0.008

0.015

0.005

0.010

0.008

0.009

0.007

0.007

0.008

0.002

0.006

0.067

0.067

0.067

0.044

0.041

0.027

0.010

0.005

0.002

0.000

0.000

0.000

0.000

0.000

0.023

0.027

0.041

0.057

0.062

0.066

0.067

0.067

0.284±

0.548t

0.623±

0.687±

0.703±

0. 738±

0.777±

0.799±

0.829±

0.860±+

0.913±

0.001

0.006

0.001

0.005

0.005

0.004

0.007

0.005

0.002

0.007

0.005

0.035

0.035

0.035

0.023

0.021

0.014

0.005

0.003

0.001

0.000

0.000

0.000

0.000

0.000

0.012

0.014

0.021

0.029

0.032

0.034

0.035

0.035

A.14
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

QI = 1650 cc/min T I282 °K T = 280 °K

Air

Flowrate

(1000 cc/sec)

Channel Pressure Liquid Flowrate

Velocity Drop

(m/sec) (mm-H20)

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

4.365+

4.714±

5.412±

5.599±

6.167±

6.371±

6.666±

7.070±

7.424±

7.922±

8.293±

8.587±

0.054

0.050

0.114

0.023

0.036

0.078

0.095

0.043

0.048

0.044

0.029

0.025

3.876

4.186

4.806

4.972

5.476

5.658

5.920

6.278

6.592

7.035

7.364

7.625

72.5 ±

79.8 ±

92.0±

94.0±

91.7±

105.2±

110.9 ±

112.0 ±

108.7 ±

104.9 ±

105.4 ±

79.1 ±

6.3

6.1

11.8

11.6

10.3

10.9

9.8

8.4

7.7

7.4

7.5

2.5

1649.3

1649.3

1335.2

1217.4

981.7

824.7

628.3

353.4

235.6

117.8

39.3

0.0

0.0

7.9

314.2

432.0

667.6

824.7

1021.0

1295.9

1413.7

1531.5

1610.1

1649.3

Nondimensional Form of Test Data

K= 0.149 j/ 0.077 Reý = 313.5 V = 0.750

C

9 c J9 Jf Jc

1

2

3

4

5

6

7

8

9

10

11

12

0.843±

0.911±

1.046±

1.082±

1.191±

1.231±

1.288±

1.366±

1.434±

1.530±

1.602±

1.659±

0.010

0.010

0.022

0.004

0.007

0.015

0.018

0.008

0.009

0.009

0.006

0.005

0.073±

0.080±

0.092±

0.094±

0.092±

0.105±

0.111±

0.112±

0.109±

0.105±

0.105±

0.079±

0.006

0.006

0.012

0.012

0.010

0.011

0.010

0.008

0.008

0.007

0.008

0.003

0.149

0.149

0.121

0.110

0.089

0.074

0.057

0.032

0.021

0.011

0.004

0.000

0.000

0.001

0.028

0.039

0.060

0.074

0.092

0.117

0.128

0.138

0.145

0.149

0.435±

0.470±

0.539±

0.558±

0.615±

0.635±

0.665±

0.705±

0.740±

0.790±

0.827±

0.856±

0.005

0.005

0.011

0.002

0.004

0.008

0.010

0.004

0.005

0.004

0.003

0.002

0.077

0.077

0.062

0.057

0.046

0.038

0.029

0.016

0.011

0.005

0.002

0.000

0.000

0.000

0.015

0.020

0.031

0.038

0.048

0.060

0.066

0.071

0.075

0.077

A.15
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

QI = 2800 cc/min TI= 282 °K T = 280 °K

9

Air

Flowrate

(1000 cc/sec)

Channel Pressure

Velocity

(m/sec)

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

2.545±

3.614±

3.826±

4.340±

4.640±

5.086±

5.574±

6.230±

6.604±

6.892±

7.335±

7.846±

8.316±

8.725±

9.311±

0.135

0.173

0.113

0.070

0.087

0.077

0.028

0.023

0.077

0.045

0.022

0.045

0.025

0.085

0.048

2.260

3.209

3.397

3.854

4.120

4.516

4.950

5.532

5.865

6.120

6.514

6.967

7.385

7.748

8.268

54.2 ±

84.7±

88.7±

101.2±

101.8±

97.2±

95.3±

106.8±

113.9 ±

115.4 ±

108.7 ±

92.6 ±

94.6 ±

91.9 ±

77.6 +

6.3

11.5

11.8

13.3

12.4

11.3

10.8

11.1

8.6

7.6

8.1

9.2

9.7

6.7

2.8

2827.4

2811.7

2772.5

2200.0

2042.0

1492.3

1256.6

942.5

628.3

455.5

274.9

117.8

39.3

0.0

0.0

0.0

15.7

55.0

628.3

785.4

1335.2

1570.8

1885.0

2199.1

2371.9

2552.5

2709.6

2788.1

2827.4

2827.4

Nondimensional Form of Test Data

Kl= 0.255 jI= 0.132 Reý = 532.0 V = 0.662

C

K * * K * . .* .* 9Kf c g f c

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

0.492± 0.026

0.698± 0.033

0.739± 0.022

0.838± 0.014

0.896± 0.017

0.982± 0.015

1.077± 0.005

1.204± 0.004

1.276± 0.015

1.331± 0.009

1.417± 0.004

1.516± 0.009

1.607± 0.005

1.686± 0.016

1.799± 0.009

0.054±

0.085±

0.089±

0.101±

0.102±

0.097±

0.095±

0.107±

0.114±

0.116±

0.109±

0.093±

0.095±

0.092±

0.078±

0.006

0.012

0.012

0.013

0.012

0.011

0.011

0.011

0.009

0.008

0.008

0.009

0.010

0.007

0.003

0.255

0.254

0.250

0.199

0.184

0.135

0.113

0.085

0.057

0.041

0.025

0.011

0.004

0.000

0.000

0.000

0.001

0.005

0.057

0.071

0.121

0.142

0.170

0.199

0.214

0.230

0.245

0.252

0.255

0.255

0.254±

0.360±

0.381±

0.433±

0.463±

0.507±

0.556±

0.621±

0.658±

0.687±

0.731±

0.782±

0.829±

0.870±

0.928±

0.013

0.017

0.011

0.007

0.009

0.008

0.003

0.002

0.008

0.004

0.002

0.004

0.002

0.008

0.005

0.132

0.131

0.129

0.103

0.095

0.070

0.059

0.044

0.029

0.021

0.013

0.005

0.002

0.000

0.000

0.000

0.001

0.003

0.029

0.037

0.062

0.073

0.088

0.102

0.111

0.119

0.126

0.130

0.132

0.132

A.16
Test Data for Two Phase Air Flow in Single Grid Rod Bundle

= 3900 cc/min T = 282 °K T = 280 °K

9

Air Channel Pressure Liquid Flowrate

Flowrate Velocity Drop Downward Upward

(1000 cc/sec) (m/sec) (mm-H20) (cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

3.326±

3.778±

4.199±

4.836±

5.490±

6.013±

6.403±

6.758±

7.177±

7.542±

7.913±

8.497±

9.105±

0.137

0.153

0.078

0.032

0.067

0.093

0.073

0.063

0.079

0.050

0.035

0.061

0.059

2.953

3.355

3.729

4.295

4.875

5.339

5.686

6.001

6.373

6.697

7.027

7.545

8.086

97.0±

109.3±

106.8±

103.9±

101.0±

102.6±

117.0 ±

115.4 ±

112.6 ±

106.5 ±

92.9 ±

77.5 ±

85.6 ±

15.6

14.2

13.0

11.7

10.8

12.4

9.6

9.2

8.1

8.3

8.1

3.5

5.4

3887.7

3479.3

3102.3

1885.0

1335.2

1021.0

667.6

526.2

353.4

223.8

110.0

7.9

0.0

39.3

447.7

1138.8

2042.0

2591.8

2906.0

3259.4

3400.8

3573.6

3703.2

3817.0

3919.1

3927.0

Nondimensional Form of Test Data

K/ 0.351 j 0.181 Ref = 741.0 V = 0.556

C

K* pK * K. .* .*

P Kf c Jg Jf Jc

1

2

3

4

5

6

7

8

9

10

11

12

13

0.642±

0.730±

0.811±

0.934±

1.061±

1.162±

1.237±

1.305±

1.387±

1.457±

1.529±

1.641±

1.759±

0.027

0.030

0.015

0.006

0.013

0.018

0.014

0.012

0.015

0.010

0.007

0.012

0.011

0.097±

0.109±

0.107±

0.104±

0.101±

0.103±

0.117±

0.116±

0.113±

0.107±

0.093±

0.078±

0.086±

0.016

0.014

0.013

0.012

0.011

0.012

0.010

0.009

0.008

0.008

0.008

0.003

0.005

0.351

0.314

0.280

0.170

0.121

0.092

0.060

0.048

0.032

0.020

0.010

0.001

0.000

0.004

0.040

0.103

0.184

0.234

0.262

0.294

0.307

0.323

0.334

0.345

0.354

0.355

0.332±

0.377±

0.419±

0.482±

0.547±

0.599±

0.638±

0.674±

0.715±

0.752±

0.789±

0.847±

0.908±

0.014

0.015

0.008

0.003

0.007

0.009

0.007

0.006

0.008

0.005

0.003

0.006

0.006

0.181

0.162

0.145

0.088

0.062

0.048

0.031

0.025

0.016

0.010

0.005

0.000

0.000

0.002

0.021

0.053

0.095

0.121

0.135

0.152

0.158

0.167

0.173

0.178

0.183

0.183

A.17
Test Data for Single Phase Air

Q I = 0 cc/min TI= 282 *K

Flow

T

in Double Grid Rod Bundle

280 * K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

41

5

6

7

8

9

10

11

2.706±

3.993±

5.119±

6.248±

7.205±

7.645±

8.323±

8.579±

9.093±

9.496±

9.722±

0.024

0.027

0.065

0.014

0.078

0.044

0.017

0.023

0.063

0.068

0.052

2.403

3.546

4.546

5.548

6.398

6.788

7.391

7.619

8.075

8.433

8.633

8.991

13.8 ±

22.5 ±

29.7 ±

35.8 ±

46.1 ±

55.7 ±

61.7 ±

72.4 ±

78.0 ±

87.3 ±

90.2 ±

0.4

1.7

1.7

2.2

5.5

2.7

5.4

4.6

1.4

2.5

4.9

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

0.0

12 10.125± 0.066 104.9 ± 4.6 0.0

Nondimensional Form of Test Data

K= 0.0 0.0 Re£= 0.0 V* = 0.0

C

K 9 K f Kc Jg Jf Jc

1

2

3

4

5

6

7

8

9

10

11

12

0.523±

0.771±

0.989±

1.207±

1.392±

1.477±

1.608±

1.657±

1.757±

1.834±

1.878±

1.956±

0.005

0.005

0.012

0.003

0.015

0.008

0.003

0.005

0.012

0.013

0.010

0.013

0.014±

0.023±

0.030±

0.036±

0.046±

0.056±

0.062±

0.072±

0.078±

0.087±

0.090±_

0.105±

0.000

0.002

0.002

0.002

0.005

0.003

0.005

0.005

0.001

0.002

0.005

0.005

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.270±

0.398±

0.510±

0.623±

0.718±

0.762±

0.830±

0.855±

0.906±

0.947±

0.969±

1.009±

0.002

0.003

0.006

0.001

0.008

0.004

0.002

0.002

0.006

0.007

0.005

0.007

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

0.000

A.18
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

QI =170 cc/min T 282 °K T = 280 °K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

2.847±

4.666±

5.574±

6.421±

6.756±

7.074±

7.366±

7.624±

8.110±

0.011

0.052

0.028

0.070

0.053

0.037

0.023

0.020

0.043

2.529

4.144

4.949

5.702

6.000

6.282

6.541

6.771

7.202

7.540

0.2

20.2

32.3

40.0

55.5

58.7

81.5

87.4

87.1

±

±

±

±

±

±

±

±

±

2.1

3.0

0.7

3.2

1.0

2.1

6.6

1.1

1.6

172.8

172.8

172.8

128.8

94.2

62.8

23.6

3.1

0.0

0.0

0.0

0.0

0.0

44.0

78.6

110.0

149.2

169.7

172.8

10 8.491± 0.064 66.2 ± 3.2 172.8

Nondimensional Form of Test Data

K 0.016 j 0.008 Re. -

32.3 V = 1.059

C

K* P* * K* " * "* SKf c 9 f Ic

1

2

3

4

5

6

7

8

9

10

0.550±

0.901±

1.077±

1.240±

1.305±

1.367±

1.423±

1.473±

1.567±

1.640±-

0.002

0.010

0.005

0.014

0.010

0.007

0.005

0.004

0.008

0.012

0.000±

0.020±

0.032±

0.040±

0.056±

0.059±

0.082±

0.088±

0.087±

0.066±_

0.002

0.003

0.001

0.003

0.001

0.002

0.007

0.001

0.002

0.003

0.016

0.016

0.016

0.012

0.009

0.006

0.002

0.000

0.000

0.000

0.000

0.000

0.000

0.004

0.007

0.010

0.013

0.015

0.016

0.016

0.284±

0.465±

0.556±

0.640±

0.673±

0.705±

0.734±

0.760±

0.808±

0.846+

0.001

0.005

0.003

0.007

0.005

0.004

0.002

0.002

0.004

0.006

0.008

0.008

0.008

0.006

0.004

0.003

0.001

0.000

0.000

0.000

0.000

0.000

0.000

0.002

0.004

0.005

0.007

0.008

0.008

0.008

A. 19
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

0I= 510 cc/min T I= 282 °K T = 280 *K 9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

1.922± 0.065 1.707 34.2 ± 0.8

3.906± 0.030 3.468 38.6 ± 1.3

5.134± 0.063 4.559 55.8 ± 0.6

5.848± 0.056 5.193 36.8 ± 1.4

6.191± 0.018 5.498 43.6 ± 6.0

6.515± 0.032 5.786 47.7 ± 7.0

6.792± 0.018 6.031 53.9 ± 4.7

7.102± 0.044 6.307 87.0 ± 1.8

7.314± 0.029 6.495 82.5 ± 1.3

7.594± 0.039 6.744 89.8 ± 5.3

8.042± 0.048 7.141 87.6 ± 0.8

8.554± 0.035 7.596 63.3 ± 3.7

Nondimensional Form of Test Data

510.5 0.0

510.5 0.0

502.7 7.9

361.3 149.2

255.3 255.3

176.7 333.8

137.4 373.1

58.9 451.6

39.3 471.2

0.0 510.5

0.0 510.5

0.0 510.5

K/ = 0.046 j, = 0.024 Re = 96.9 V* = 0.882 C

K • P K * K *i*i "* " "* K f c J f c

1

2

3

4

5

6

7

8

9

10

11

12

0.371±

0.754±

0.992±

1.130±

1.196±

1.259±

1.312±

1.372±

1.413±

1.467±

1.554±

1.653±

0.013

0.006

0.012

0.011

0.003

0.006

0.004

0.009

0.006

0.007

0.009

0.007

0.034±

0.039±

0.056±

0.037±

0.044±

0.048±

0.054±

0.087±

0.083±

0.090±

0.088±

0.063±

0.001

0.001

0.001

0.001

0.006

0.007

0.005

0.002

0.001

0.005

0.001

0.004

0.046

0.046

0.045

0.033

0.023

0.016

0.012

0.005

0.004

0.000

0.000

0.000

0.000

0.000

0.001

0.013

0.023

0.030

0.034

0.041

0.043

0.046

0.046

0.046

0.192±

0.389±

0.512±

0.583±

0.617±

0.649±

0.677±

0.708±

0.729±

0.757±

0.802±

0.853±

0.006

0.003

0.006

0.006

0.002

0.003

0.002

0.004

0.003

0.004

0.005

0.003

0.024

0.024

0.023

0.017

0.012

0.008

0.006

0.003

0.002

0.000

0.000

0.000

0.000

0.000

0.000

0.007

0.012

0.016

0.017

0.021

0.022

0.024

0.024

0.024

A.20
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

0 = 750 cc/min TI= 282 °K T9= 280 *K 9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

2.329± 0.232

4.261± 0.074

4.956± 0.062

5.614± 0.027

6.116± 0.044

6.312± 0.063

6.633± 0.086

7.071± 0.038

7.407± 0.030

7.862± 0.020

8.386± 0.042

5.228± 0.019

5.589± 0.022

2.068

3.784

4.401

4.985

5.431

5.605

5.890

6.279

6.578

6.981

7.447

4.643

4.963

36.5 ± 1.4 746.1 0.0

22.9 ± 1.5 746.1 0.0

33.2 ± 0.7 687.2 58.9

38.9 ± 2.2 451.6 294.5

44.7 ± 7.3 314.2 432.0

53.8 ± 5.1 216.0 530.1

57.8 ± 3.1 157.1 589.0

65.9 ± 5.3 58.9 687.2

58.6 ± 1.7 19.6 726.5

57.8 ± 1.3 0.0 746.1

68.0± 10.4 0.0 746.1

59.0 ± .2.6 628.3 117.8

58.8 ± 3.1 510.5 235.6

Nondimensional Form of Test Data

K 0.067 j= 0.035 Re. = 142.5 V = 0.838

C

K Kf Kc J9 Jf Jc

1

2

3

4

5

6

7

8

9

10

11

12

13

0.450±

0.823±

0.957±

1.084±

1.181±

1.219±

1.281±

1.366±

1.431±

1.5199±

1.620±

1.0100±

1.080±

0.045

0.014

0.012

0.005

0.009

0.012

0.017

0.007

0.006

0.004

0.008

0.004

0.004

0.037±

0.023-L

0.033±

0.039±

0.045±

0.054±

0.058±

0.066±

0.059±

0.058±

0.068±

0.059±

0.059±

0.001

0.001

0.001

0.002

0.007

0.005

0.003

0.005

0.002

0.001

0.010

0.003

0.003

0.067

0.067

0.062

0.041

0.028

0.020

0.014

0.005

0.002

0.000

0.000

0.057

0.046

0.000

0.000

0.005

0.027

0.039

0.048

0.053

0.062

0.066

0.067

0.067

0.011

0.021

0.2322±

0.425±

0.494±

0.560±

0.6100±

0.629±

0.661±

0.705±

0.738±

0.784±

0.836±

0.521±

0.557±

0.023

0.007

0.006

0.003

0.004

0.006

0.009

0.004

0.003

0.002

0.004

0.002

0.002

0.035

0.035

0.032

0.021

0.015

0.010

0.007

0.003

0.001

0.000

0.000

0.029

0.024

0.000

0.000

0.003

0.014

0.020

0.025

0.027

0.032

0.034

0.035

0.035

0.005

0.011

A.21
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

Q 1500 cc/min T =282 'K T9= 280 °K

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

3.246± 0.100

3.989± 0.023

4.391± 0.031

4.916± 0.045

5.198± 0.032

5.573± 0.027

6.201± 0.020

6.583± 0.061

7.010± 0.053

7.313± 0.030

7.723± 0.081

8.319+ 0.028

2.882

3.542

3.899

4.365

4.616

4.949

5.506

5.846

6.225

6.494

6.858

7.387

45.6

55.8

57.1

64.7

66.5

71.1

84.5

96.1

79.9

64.5

63.0

±

±

±

±

±

±

±

±

±

±

±

6.7

1.0

1.3

4.4

4.0

5.2

5.3

7.1

1.7

3.3

8.8

1492.3 0.0

1413.7 78.5

1295.9 196.3

942.5 549.8

863.9 628.3

628.3 863.9

353.4 1138.8

219.9 1272.3

78.5 1413.7

47.1 1445.1

0.0 1492.3

80.7 ± 3.9 0.0 1492.3

Nondimensional Form of Test Data

KI*= 0.135 j/= 0.070 Reý = 285.0 V* = 0.662

C

gA

K g J f Jc

1

2

3

4

5

6

7

8

9

10

11

12

0.627±

0.771±

0.848±

0.950±

1.004±

1.077±

1.198±

1.272±

1.354±

1.413±

1.492±

1.607±

0.019

0.004

0.006

0.009

0.006

0.005

0.004

0.012

0.010

0.006

0.016

0.005

0.046±

0.056±

0.057±

0.065±

0.067±

0.071±

0.085±

0.096±

0.080±

0.065±

0.063±

0.081±

0.007

0.001

0.001

0.004

0.004

0.005

0.005

0.007

0.002

0.003

0.009

0.004

0.135

0.128

0.117

0.085

0.078

0.057

0.032

0.020

0.007

0.004

0.000

0.000

0.000

0.007

0.018

0.050

0.057

0.078

0.103

0.115

0.128

0.130

0.135

0.135

0.324±

0.398±

0.438±

0.490±

0.518±

0.555±

0.618±

0.656±

0.699±

0.729±

0.770+±

0.829±

0.010

0.002

0.003

0.004

0.003

0.003

0.002

0.006

0.005

0.003

0.008

0.003

0.070

0.066

0.060

0.044

0.040

0.029

0.016

0.010

0.004

0.002

0.000

0.000

0.000

0.004

0.009

0.026

0.029

0.040

0.053

0.059

0.066

0.067

0.070

0.070

A. 22
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

Q = 2500 cc/min T,= 282 °K T = 280 °K I I 9

Air

Flowrate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

0.723± 0.152

1.722± 0.188

2.763± 0.049

3.919± 0.033

4.330± 0.070

4.731± 0.051

5.258± 0.050

5.546± 0.022

5.689± 0.053

6.170± 0.034

6.283± 0.051

6.646± 0.092

6.855± 0.023

7.386± 0.020

8.133± 0.056

0.642

1.530

2.454

3.480

3.845

4.201

4.669

4.925

5.052

5.479

5.580

5.901

6.087

6.558

7.222

33.2

51.4

44.2

56.4

58.0

60.6

73.4

77.2

78.7

90.6

103.8

103.6

96.3

87.1

±

±

±

±

±

±

±

±

±

±

±

±

±

0.6

6.1

6.4

1.4

2.3

3.6

5.1

3.6

3.6

8.7

7.9

7.6

8.0

4.3

2513.3 0.0

2513.3 0.0

2497.6 15.7

2081.3 432.0

1649.3 863.9

1453.0 1060.3

1060.3 1453.0

903.2 1610.1

785.4 1727.9

573.3 1939.9

392.7 2120.6

235.6 2277.7

196.3 2316.9

62.8 2450.4

67.1 ± 3.0 0.0 2513.3

Nondimensional Form of Test Data

K 0.227 j 0.117 Rer, = 475.0 V = 0.574 C

K * K * .K .* .*

9 Kf Kc f c

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

0.140± 0.029

0.333± 0.036

0.534± 0.009

0.757± 0.006

0.836± 0.013

0.914± 0.010

1.016± 0.010

1.071± 0.004

1.099± 0.010

1.192± 0.007

1.214± 0.010

1.284± 0.018

1.324± 0.004

1.427± 0.004

1.571± 0.011

0.033±

0.051±

0.044±

0.057±

0.058±

0.061±

0.073±

0.077±

0.079±

0.091±

0.104±

0.104±

0.096±

0.087±

0.067±

0.001

0.006

0.006

0.001

0.002

0.004

0.005

0.004

0.004

0.009

0.008

0.008

0.008

0.004

0.003

0.227

0.227

0.226

0.188

0.149

0.131

0.096

0.082

0.071

0.052

0.035

0.021

0.018

0.006

0.000

0.000

0.000

0.001

0.039

0.078

0.096

0.131

0.145

0.156

0.175

0.191

0.206

0.209

0.221

0.227

0.072±

0.172±

0.275±

0.391±

0.432±

0.472±

0.524±

0.553±

0.567±

0.615±

0.626±

0.662±

0.683±

0.736±

0.811±

0.015

0.019

0.005

0.003

0.007

0.005

0.005

0.002

0.005

0.003

0.005

0.009

0.002

0.002

0.006

0.117

0.117

0.116

0.097

0.077

0.068

0.049

0.042

0.037

0.027

0.018

0.011

0.009

0.003

0.000

0.000

0.000

0.001

0.020

0.040

0.049

0.068

0.075

0.081

0.090

0.099

0.106

0.108

0.114

0.117

A.23
Test Data for Two Phase Air Flow in Double Grid Rod Bundle

I =3700 cc/min Ti= 282 °K T = 280 °K

g

Air

Flow rate

(1000 cc/sec)

Channel

Velocity

(m/sec)

Pressure

Drop

(mm-H20)

Liquid Flowrate

Downward Upward

(cc/min) (cc/min)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

1.079±

1.661±

2.708±

3.199±

3.405±

3.888±

4.299±

4.837±

5.493±

5.809±

6.264±

6.798±

7.072±

7.438±

8.266±

0.317

0.311

0.039

0.108

0.038

0.021

0.076

0.028

0.074

0.057

0.030

0.020

0.037

0.043

0.045

0.959

1.475

2.405

2.841

3.023

3.452

3.818

4.295

4.878

5.158

5.563

6.037

6.280

6.605

7.340

59.4 ±

68.8±

55.5 ±

56.4 ±

56.6 ±

58.7 ±

60.5 ±

.68.7 ±

76.7 ±

84.5 ±

108.6 ±

101.6±

98.1 ±

82.2 ±

98.8 ±

4.5

10.7

2.9

2.1

2.5

3.4

4.0

5.4

5.5

4.3

8.6

10.1

8.4

4.8

4.8

3691.4

3675.7

3495.0

2748.9

2434.7

1885.0

1570.8

1335.2

863.9

706.9

432.0

235.6

133.5

55.0

0.0

0.0

15.7

196.3

942.5

1256.6

1806.4

2120.6

2356.2

2827.4

2984.5

3259.4

3455.8

3557.9

3636.4

3691.4

Nondimensional Form of Test Data

I•= 0.333 j = 0.172 Re£ = 703.0 V* = 0.424

C

g Kf c Jg Jf Jc

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

0.209±

0.321±

0.523±

0.618±

0.658±

0.751±

0.831±

0.934±

1.061±

1.122±

1.210±

1.313±

1.366±

1.437±

1.597±

0.061

0.060

0.008

0.021

0.007

0.004

0.015

0.005

0.014

0.011

0.006

0.004

0.007

0.008

0.009

0.059± 0.004

0.069± 0.011

0.056± 0.003

0.056± 0.002

0.057± 0.003

0.059± 0.003

0.061± 0.004

0.069± 0.005

0.077± 0.005

0.085± 0.004

0.109± 0.009

0.102± 0.010

0.098± 0.008

0.082± 0.005

0.099± 0.005

0.333

0.332

0.316

0.248

0.220

0.170

0.142

0.121

0.078

0.064

0.039

0.021

0.012

0.005

0.000

0.000

0.001

0.018

0.085

0.113

0.163

0.191

0.213

0.255

0.270

0.294

0.312

0.321

0.328

0.333

0.108±

0.166±

0.270±

0.319±

0.339±

0.388±

0.429±

0.482±

0.548±

0.579±

0.624±

0.678±

0.705±

0.741±

0.824±

0.032

0.031

0.004

0.011

0.004

0.002

0.008

0.003

0.007

0.006

0.003

0.002

0.004

0.004

0.005

0.172

0.171

0.163

0.128

0.113

0.088

0.073

0.062

0.040

0.033

0.020

0.011

0.006

0.003

0.000

0.000

0.001

0.009

0.044

0.059

0.084

0.099

0.110

0.132

0.139

0.152

0.161

0.166

0.169

0.172

A.24
APPENDIX B. INPUT DECKS FOR CALCULATIONS

Aopendix B. Input Decks for Calculations

15$$11"14 $ $5 $"$$S$S $$$$$$

ccfi model test dla=0.5 cm, jg=O - 1.0

ccfl junction = 4000000

Sccfi option = wallis type ( m=1.0 c-1.0]

000100 new transnt

0000102 sl si

0000110 air

0000115 1.0 * --------------------------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 1000 10000 * -.................................................

0010000 waterin tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0 *z --------------------------------------------------

0020000 waterIn tmdpjun

0020101 001000000 003000000 5.0-6

0020200 1

0020201 0. 6.954-4 0. 0. * jf = 0.16 *----------------------------------

0030000 section pipe

0030001 6

0030101 1.9635-5 8

0030301 0.30 6

0030601 -90.0 6

0030801 0.0 0.005 6

0031001 00000 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 8

0031301 0.0 0.0 0.0 5 .....................................................

0040000 downflow snglJun

0040101 003010000 005000000 1.9535-5 0.0 0.0 100000

0040110 0.005 0.0 1.0 1.0 *** wallis flooding

0040201 1 0.0 0.0 0.0 ----------------------------------

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0

*.................................................

0060000 alrin tmdojun

0060101 007000000 005010003 1.9635-5

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 7.000 0.0

* ----------------

0070000 airin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0 * -----------------------------------------------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0 ..............................---------------------

0090000 tanki snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0 *..............................-------------------

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0

* --------------------------------------------

0110000 upplenum tmdpvoi

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0 ----------------------------------------------

20500100 I-In mult 5.0-6 0.0 0

20500101 voidfj,002000000 velfJ,002000000 * -----------------------------------

20500200 I-downl mult 1.9635-5 0.0 0

20500201 voldfJ,004000000 velfj,004000000 *a------------------------------------------------

20500300 I-up mult 1.9635-5 0.0 0

20500301 voldfj,010000000 velfj,010000000

20500400 J-f mult 4.5180 0.0 0

20500401 voldfj,004000000 velfj,004000000

20500500 J-g mult 0.1584 0.0 0

20500501 void gJ,004000000 velgj,004000000

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3080000 -1.0 p,3010000 *------------------------------------

20500700 J-fs mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4 *-------

20500800 J-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5 * --------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 sort cntrlvar,7 *-------------------------

20501000 J-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8

*a---------------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

*a- -----------------------

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

00Q0304 cntrivar,4

0000305 cntrlvar,5

0000306 cntrlvar,6

0000307 cntrlvar,9

0000308 cntrlvar,10

0000309 cntrivar,1l1

t-----------------------------------------------

* end of Input data
sSS$SS$$$$S$$$$$$$$$$$$$$$$$$$$t$$$$$$$~SSSSS$$$S~SS~$sS$$$$$$$$$

* ccfi model test dia= I cm, jg'O - 1.0

* ccfi junction = 4000000

* ccfl option = waills type ( m=1.0 c=1.0)

000100 new transnt

0000102 sl si

0000110 air

0000115 1.0

* -------- 0------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 1000 10000

* ---------------------------------

0010000 waterin tmdpvoI

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

* -------------- ------------------

0020000 waterin tmdpjun

0020101 001000000 003000000 0.00002

0020200 t

0020201 0. 3.934-3 0. 0. * Jf = 0.18

* ---------------------------------

0030000 section pipe

0030001 6

0030101 7.854-5 6

0030301 0.30 6

0030801 -90.0 6

ED 0030801 0.0 0.01 8

tO 0031001 00000 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0,0

*----------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0W10200 4

0110201 0.0 1.0+5 282.0 1.0

* ------------------------------------

20500100 I-In mutt 0.00002 0.0 0

20500101 voldfJ,002000000 velfJ,002000000

*------------------------------------

20500200 1-down1 mutt 7.8540-5 0.0 0

20500201 voldfJ,004000000 velfj,004000000

* ------------------------------------------------

20500300 I-up muit 7.8540-5 0.0 0

20500301 voidfJ,010000000 velfj,010000000

*------------------------------------

20500400 J-f mutt 3.1947 0.0 0

20500401 voidfj,004000000 velfJ,O04000000 *- ----------- --------------

20500500 J-g mutt 0.1120 0.0 0

20500501 void gj,004000000 ve0gj,004000000

*------------------------------------

20500800 diff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3060000 -1.0 P.3010000

*------------------------------------

20500700 J-fs mutt 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 J-gs muit -1.0 0.0 0 1 0.0

20500801 cntrivar,5

* --------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

* -----------------

20501000 J-9-sart stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8

* ---------------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrivar,lO

*-------------------

0000301 cntrivar,I

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrtvar,4

0000305 cntrlvar,S

0000308 cntrlvar,5

0000307 cntrlvar,9

000•308 cntrivar,10

0000309 cntrlvar,l1

* end of Input data

0040000

0040101

0040110

0040201

0050000

0050101

0050200

downflow snglJun

003010000 005000000 7.8540-5 0.0 0.0 100000

0.01 0.0 1.0 1.0 *** waills flooding

1 0.0 0.0 0.0

lowerpi snglvol

0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

4 1.0+5 282.0 1.O

* ---------------------------------

0060000 alrin tmdpjun

0060101 007000000 005010003 7.8540-5

0060200 0

0080201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 8.780 0.0

*-r -------------------------------

0070000 alrin tmdovol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.45 282.0 1.0

* -----------------------------------------------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

* --------------------------------------------------

0090000 tanki snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0

* --------------------------------------------------
- ccfl model test dla= 2 cm , jg=O - 1.0

* ccfl junction = 4000000

* ccfl option = wallis type ( m=1.o c=1.0

000100 new transnt

0000102 sl sl

0000110 air

0000115 1.0

* -------- 0------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 11000 10000 *----------------------------------

0010000 waterln tmdDvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0

* ---------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

* ------------------------------------

20500100 I-In mult 0.00008 0.0 0

20500101 voldfj,002000000 velfj,002000000

*------------------------------------

20500200 I-downl mult 3.1416-4 0.0 0

20500201 voldfj,004000000 velfJ,004000000

* -----------------------------------

20500300 I-up mult 3.1416-4 0.0 0

20500301 voidfJ,010000006 velfj,010000000

20800007 flJ,3010000

20800008 flJ,3020000

20800009 fl),3030000

20800010 flj,3040000

20800011 flJ,3050000

20800012 flj,4000000

* ------------------

20800013 lregj,3010000

20800014 Iregj,3020000

20800015 lregJ,3030000

20800016 IregJ,3040000

20800017 Iregj,3050000

20800018 Iregj,4000000

* ------------------

20800019 florgJ,3010000

20800020 florgj,3020000

20800021 florgj,3030000

20800022 florgj,3040000

20800023 florgJ,3050000

20800024 florgJ,4000000

* ------------------

20800025 vgJj,3010000

20800026 vgjj,3020000

20800027 vgjj,3030000

20800028 vgJj,3040000

20800029 vgjj,3050000

20800030 vgjj,4000000

0020000 waterln tmdpjun

0020101 001000000 003000000 0.00008

0020200 I

0020201 0. 0.0223 0. 0. * Jf = 0.16

0030000 section pipe

0030001 6

0030101 3.1416-4 6

0030301 0.30 6

0030601 -90.0 6

0030801 0.0 0.03 8

(V 0031001 00000 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

* -----------------------------------

0040000 downflow sngljun

0040101 003010000 005000000 3.1416-4 0.0 0.0 10000

0040110 0.02 0.0 1.0 1.0 *** wallIs flooding

0040201 1 0.0 0.0 0.0

*----------------------------------

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0

0050200 4 1.0+5 282.0 1.0

*----------------------------------

0060000 airin tmdpjun

0060101 007000000 005010003 3.1416-4

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 15.00 0.0 ----------------------------------

0070000 alrin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0

0070200 4

0070201 0.0 1.+5 282.0 1.0

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 0000(

0080201 1 0.0 0.0 0.0

*- - - --- ---------

20500400 J-f mult 2.2590 0.0 0

20500401 voldfj,004000000 velfJ,004000000

-----------

0

C

20500500 J-g mult 0.07920 0.0 0

20500501 voldgJ,004000000 velgJ,004000000

*------------------------------------

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3060000 -1.0 P,3010000 *------------------------------------

20500700 J-fs mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 J-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5

* --------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 sort cntrlvar,7 * -------------------

20501000 J-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8

* ------------------

20501100 curve sum 1.0 0.0 0

0 20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

*-------------------

20501200 massf mult 3183.00 0.0 0

205012011 mflowJ,004000000

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0090305 cntrlvar,5

0000306 cntrlvar,6

O 0000307 cntrlvar,9

0000308 cntrlvar.10

0000309 cntrlvar,11

0000310 cntrlvar,12

* --------------------

20800001 voldJ,3010000

20800002 void J,3020000

20800003 voldJ,3030000

20800004 voldJ,3040000

20800005 voidJ,3050000

20800006 vold 1,4000000 *--------------------------------------

0090000 tank1 snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0 k -------------------- -----------------------------.
$s$$s$s$$sss$s$ss$sssssss$ss$ss$$$$$s$sSsss$s$$$Ssss$s$ssSsssSSsss

* ccfl model test diae 3 cm, jg=O - 1.0

* ccfl junction = 4000000

* ccfl option = wallis type ( m=1.0 c=1.0

000100 new transnt

0000102 si sl

0000110 air

0000115 1.0

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 1000 10000

* ---------------------------------

0010000 waterln tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

*----------------------------------

0020000 waterln tmdpjun

0020101 001000000 003000000 0.00018

0020200 I

0020201 0. 0.0813 0. 0. * Jf = 0.16

Ur

0030000 section pipe

0030001 6

0030101 7.0690-4 6

0030301 0.30 6

0030601 -90.0 8

0030801 0.0 0.03 6

0031001 00000 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

* -----------------------------------------------------

0040000 downflow snglJun

0040101 003010000 005000000 7.0690-4 0.0 0.0 100000

0040110 0.03 0.0 1.0 1.0 *** wallis flooding

0040201 I1 0.0 0.0 0.0

*----------------------------------

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0

*----------------------------------

0060000 alrln tmdpjun

0060101 007000000 005010003 7.0690-4

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 15.20 0.0

*-L -------------------------------

0070000 alrin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0

*----------------------------------

0100000 Upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.? 0000

0100201 1 0.0 0.0 0.0

* ---------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

* ------------------------------------

20500100 I-In mult 0.00018 0.0 0

20500101 voldfj,002000000 velfj,002000000

*------------------------------------

20500200 I-downi mult 7.0690-4 0.0 0

20500201 voidfj,004000000 velf),004000000

* -----------------------------------

20500300 I-up mult 7.0690-4 0.0 0

20500301 voldfJ,010000000 velfj,01000000o

* -----------------------------------

20500400 J-f mult 1.8444 0.0 0

20500401 voldfj,004000000 velfi,004000000

* -----------------------------------

20500500 J-9 mult 0.08488 0.0 0

20500501 voldgJ,004000000 velgJ,004000000

* -----------------------------------

20500600 dlff-p sum 1.0 0.0 0

20500801 0.0 1.0 P,3060000 -1.0 p.3010000

* -----------------------------------

20500700 J-fs mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 J-gs mult -1.0 0.0 0 1 0.0

20500801 , cntrlvar,5

* --------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 aqrt cntrlvar,7

20501000 J-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8 * ---------------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

*-------------------

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,5

0000308 cntrlvar,S

0000307 cntrlvar,9

0060308 cntrlvar,10

0000309 cntrlvar,111

-------------------

* end of input data

* -----------------------------------------------------

0080000 tankflow sngljun

0000101 005010000 009000000 0.0346

0080201 1 0.0 0.0 0.0

0.0 0.0 00000

* ---------------------------------

0090000 tanki sngivol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0
ss$$$$$$$$$$$$$s$$$$$$$$$$$$$$$$$$$s$$s$s$$$$$$$$$$s$sss$$$$$$$$$$

* ccfl model test dlaz 5 cm, jg=O - 1.0

ccfl junction = 4000000

* ccf! option = wallis type ( m:l.0 c:l.0 )

000100 new transnt

0000102 sl sl

0000110 air

0000115 1.0

* 10.----- 0------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.011 3 100 1000 10000 *----------------------------------

0100000 upflow sn9lJun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0

0110000 upnlenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

20500100 I-In mult 0.0005 0.0 0

20500101 voldfJ,002000000 velfJ,002000000

0010000 waterin tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

*.------------------------------------------------

0020000 waterln tmdpjun

0020101 001000000 003000000 0.0005

0020200 1

0020201 0. 0.219600 0. 0. * Jf 0.16 *----------------------------------

0030000 section pipe

0030001 6

0030101 1.9635-3 6

(E 0030301 0.30 6

0030601 -90.0 6

0030801 0.0 0.05 6

0031001 00000 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5 *------------------------------------

0040000 downflow sngljun

0040101 003010000 005000000 1.9635-3 0.0 0.0 100000

0040110 0.05 0.0 1.0 1.0 *** wallis flooding

0040201 1 0.0 0.0 0.0 * ------------------------------------------------

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0

* ---------------------------------------------

0060000 alrin tmdpjun

0060101 007000000 005010003 1.9635-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 19.60 0.0

0060204 300.0 0.00 29.91 0.0 *--------------------------------------------------

0070000 airin tmdovol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0

20500200 I-down1 .mult 1.9635-3 0.0 0

20500201 voldfj,004000000 velfj,004000000 *------------------------------------

20500300 i-up mult 1.9635-3 0.0 0

20500301 voldfJ,010000000 velfJ,010000000 *------------------------------------

20500400 J-f mult 1.4287 0.0 0

20500401 voldfJ,004000000 velfJ,004000000

20500500 J-g muit 0.05 0.0 0

20500501 void gj,004000000 velgj,004000000

20500600 dIlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3060000 -1.0 p.3010000 *------------------------------------

20500700 J-Is mult 1.0 0.0 0 1 0.0

20500701 cntrivar,4 * --------------

20500800 J-gs mult -1.0 0.0 0 1 0.0

20500801 cntrivar,5 * --------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrivar,7

*a--------------------------

20501000 J-g-slrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8 * ------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

* c---------------

0000301 cntrivar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,5

0000308 cntrlvar ,9

0000307 cntrlvar,9

00d0308 cntrivar,10

0000309 cntrivar.t I

* end of Input data

*.---------------------------------------------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0 *----------------------------------

0090000 tanki snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0
$$$$$$$$$$$$$$$$$$$$s~ssss$$$$$$$$$s$$$$s$s~$$$$$$sss$$s$$$$$$$$$$$

* ccfl model test tube diameter 10 cm

* ccfl junction = 4000000

* ccfl option = kutatelagze type ( m=l.0 c- sclrt[3.2])

000100 new transnt

0000102 si si

0000110 air

0000115 1.0

* 10------ 0------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 300.0 1.0-7 0.01 3 100 1000 10000

0010000 waterln tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

0020000 waterln tmdpjun

0020101 001000000 003000000 0.002

0020200 1

0020201 0. 1.244100 0. 0. * 2600 cc/min *..............................-------------------

00300C0 section pipe

0030001 6

0030101 7.8540-3 6

0030301 0.25 6

0030601 -90.0 8

0030801 0.0 0.1000 6

0031001 00000 6

0031101 00000 5

0031201 4 l.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

* -----------------------------------

Ow

0100000 upflow $ngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0 * ---------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0 * ------------------------------------------------------

20500100 I-In mult 0.0020 0.0 0

20500101 voldfj,002000000 velfj,002000000 * -----------------------------------------------------

20500200 I-downl muit 7.8540-3 0.0 0

20500201 voldfJ,004000000 vOlf J,004000000 *------------------------------------

20500300 I-up mult 7.8540-3 0.0 0

20500301 voldfj,010000000 velfJ,010000000 *------------------------------------

20500400 k-f mult 6.1000 0.0 0

20500401 voldfj,004000000 velfJ,004000OO0 * ------------------------------------

20500500 k-a mutt 0.214 0.0 0

20500501 voldgj,004000000 velgJ,004000000

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3050000 -1.0 p,3010000 * ----- ------ ---------

20500700 k-ts muit 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

20500800 k-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5

* --------------

20500900 k-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

20501000 k-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,S

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

* c---------------

0000301 cntrlvar,2

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,5

0000308 cntrcvar,8

0000307 cntrlvar,O

0000308 cntrlvar,10

0000309 cntrdvart 1 --------------------------------------

* end of Input data

0040000 downfflow sngijun

0040101 003010000 005000000 7.8540-3 0.0 0.0 100000

0040110 0.1000 1.0 1.79 1.0 *** kutateladze form

C040201 1 0.0 0.0 0.0 S....................................---............

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0 * --------------------------------------------------

0060000 airln tmdpjun

0060101 007000000 005010003 7.8540-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 20.0 0.0

0060204 300.0 0.00 30.0 0.0 * ---------------------------------

0070000 alrin tmdavol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0 * -----------------------------------------------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

0090000 tanki snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0
$$s$$s$•$ss$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$sss$s$s$$s$s$s$$1$$$$

- ccfl model test bundle type I ( no spacer grid

* ccfI option junction = 4000000

* ccfl type = kutateladze form : constant from experiment

000100 new transnt

0000102 sI sl

0000110 air

0000115 1.0

*--------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 300.0 1.0-7 0.01 3 100 1000 10000

* ---------------------------------

0010000 waterln tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

* ---------------------------------

0020000 waterln tmdpJun

0020=01 001000000 003000000 3.80-4

0020200 1

0020201 0. 0.040000 0. 0. * 2600 cc/min

* ------------------------------

0030000 section pipe

0030001 6

0030101 1.1261-3 6

0030301 0.25 6

0030601 -90.0 6

0030801 0.0 0.0103 6

0031001 00100 6

0031101 00000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

*------------------------------------

0040000 downflow sngljun

0040101 003010000 005000000 1.1261-3 0.0 0.0 100000

0040110 0.0103 1.0 1.44 1.0 *** kutateladze form

0040201 1 0.0 0.0 0.0

* ------------------------------

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0

*----------------------------------

0060000 alrln tmdpjun

0060101 007000000 005010003 1.1281-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 10.0 0.0

0069204 300.0 0.00 20.0 0.0

*----------------------------------

0070000 alrin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 I.+5 282.0 1.0

* ---------------------- -------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

* ---------------------------------

0090000 tanki snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 I 0.0 0.0 0.0

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.045 282.0 1.0

20800011 fiJ,3050000 20800012 flj,4000000

20800013 florgj,3010000

20800014 florgJ,3020000

20800015 florgJ,3030000

20800016 florgJ,3040000

20800017 florgj,3050000

20800018 florgj,4000000

* end of Input data

20500100 I-in mult 3.80-4 0.0 0 20500101 voldfJ,002000000 velf ,002000000

20500200 I-downl mult 1.1281-3 0.0 0 20500201 voldfj,004000000 velfj,004000000

X -----------------------------------------------

20500300 I-up muit 1.1261-3 0.0 0

20500301 voldfJ,010000000 velfJ,010000000

20500400 k-f mutt 6.1000 0.0 0 20500401 voldfj,004000000 velfj,004000000

20500500 k-g mult 0.214 0.0 0 20500501 voldgj,004000000 velgj,004000000

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3050000 -1.0 p,30t0000

* -----------------------------------

20500700 k-fs muit 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 k-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5

*---------------

20500900 k-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

* -----------------

20501000 k-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrIvar,8

* ------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

* ----------------

0000301 cntrlvar,l

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,5

0000308 nntrrvarS

000030? cntrlvar,9

0000308 cntrlvar,10 0000309 cntrivar,11i

20800001 20800002

20800003

20800004

20800005

20800006

20800O07

20800008

20800009

20800010

voldJ,3010000 void J,3020000

void J,3030000

voldJ,3040000

vold j,3050000

void ],4000000

fl J,3010000

f 11,3020000

f1j,3030000

f1 1,3040000
5$ $$$$$ $$$$$$$$S$$S$$$$$$$SS$$$$$SS$$$S$$$$$S$$$$$$$$$$S$$$$$$$$$

- ccfl model test bundle type 2 ( i spacer grid

, ccfl junction = 3030000 , 4000000

a ccfl kutatefadze form : constant from experiment

000100 new transnt

0000102 sI si

0000110 air

0000115 1.0

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0,01 3 100 1000 10000

* ----------------------- --------------------------

0010000 waterln tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0 *-------- --------------------------

0020000 waterln tmdpjun

0020101 001000000 003000000 3.80-4

0020200 1

0020201 0. 0.040000 0. 0. * 2600 cc/min

* ----------------------- --------------------------

0030000 section pipe

0030001 6

0030101 1.1261-3 6

0030301 0.25 6

0030601 -90.0 6

CO 0030801 0.0 0.0103 6

ý0 0031001 00100 6

0031101 000000 1

0031102 000000 2

0031103 100000 3

0031104 000000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 6

0031301 0.0 0.0 0.0 5

0031401 0.0103 1.0 1.23 1.0 1 **a Grid 01

0031402 0.0103 1.0 1.36 1.0 5 *aa Grid 02

* -----------------------------------------------------

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

* --------------------------------------------------

0090000 tank1 snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0

010000C upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0

----------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

* -.------------------- ..........-----------------------

20500100 I-In mult 3.80-4 0.0 0

20500101 voidfj,002000000 velfj,002000000

* -----------------------------------

20500200 1-downi mult 1.1261-3 0.0 0

20500201 voidfj,004000000 velfj,004000000

* -----------------------------------

20500300 I-up mult 1.1261-3 0.0 0

20500301 voldfj,010000000 velfj,0l0000000

20500400 k-f mult 6.1000 0.0 0

20500401 voldfJ, OW volfJ, o

t ------------------- Ix ---------------- . ...--- ---------------

20500500 k-g mult 0.214 0.0 0

20500501 voldgj,0cOQooJ ,•4j, velgJ 4O 100 *------.------------ 3-P -----------------------

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3050000 -1.0 p,3010000

* -----------------------------------

20500700 k-fs mult 1.0 0.0 0. 1 0.0

20500701 cntrlvar,4

* -------------------

20500800 k-ge mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5 * ---------------------

20500900 k-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

20501000 k-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8 * ------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

20501400 k-f mult 6.1000 0.0 0

20501401 voidfj,003010000 velfj,003010000

*a------------------------ 0----------------..

20501500 k-g mult 0.214 0.0 0

20501501 voidgj,003010000 velgj,003010000 *------------------------------------

20501700 k-fs mult 1.0 0.0 0 1 0.0

20501701 cntrlvar,14

* -....................

20501800 k-gs mult -1.0 0.0 0 1 0.0

20501801 cntrlvar,15

* -............. ......

20501900 k-f-eqrt stdfnctn 1.0 0.0 0

20501901 sqrt cntrlvar,17

--...........................

20502000 k-g-sqrt stdfnctn 1.0 0.0 0

20502001 sort cntrlvar,18 a ---------------------------

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,5

0000306 cntrlvar,6

0000307 cntrlvar,9

0000308 cntrlvar,10

0000309 cntrlvar,.11

0000310 cntrlvar,19

0000311 cntrlvar,20 a----------------------------------------------

a end of Input data

0040000

0040101

0040110

0040201

downflow snglJun

003010000 005000000 1.1261-3 0.0 0.0 100000

0.0103 1.0 1.44 1.0 *** kutateladze form

1 0.0 0.0 0.0

a ----------------------------------------------

0050000 lowerl snglvel

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0

0060000 airin tmdpjun

0060101 007000000 005010003 1.1261-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

*0060203 150.0 0.00 7.00 0.0

*0060204 250.0 0.00 7.00 0.0

0060203 200.0 0.00 10.0 0.0

0070000 airin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0

a ---------------------- ------------------------------

0080000 tankflOw sngljun
a0

* ccfl model test bundle type 3 ( 2 spacer grid

* ccfl Junction = 3010000, 3030000 , 4000000

* ccfI kutateladze form : constant from experiment

000100 new transnt

0000102 si si

0000110 air

0000115 1.0 .-------------------------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 1000 10000 * -------------------------------------------------

0010000 waterin tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0 .-------------------------------------------------

0020000 waterln tmdpjun

0020101 001000000 003000000 3.80-4

0020200 I

0020201 0. 0.040000 0. 0. * 2600 cc/min * --------------------------------------------------

0030000 section pipe

0030001 6

0030101 1.1261-3 6

0030301 0.25 6

0030601 -90.0 6

0030801 0.0 0.0103 6

0031001 00100 6

0031101 100000 1

0031102 000000 2

0031103 100000 3

0031104 000000 5

0031201 4 1.+5 282.0 1.000 0.0 0.0 8

0031301 0.0 0.0 0.0 5

0031401 0.0103 1.0 1.23 1.0 I ** Grid 01

0031402 0.0103 1.0 1.36 1.0 5 **G Grid 02

*------------------------------------

0040000 downflow sngijun

0040101 003010000 005000000 1.1261-3 0.0 0.0 100000

0040110 0.0103 1.0 1.44 1.0 *** kutateladze form

0040201 1 0.0 0.0 0.0 * --------------------------------------------------

0050000 lowerpl sngivol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0 .------------------------------------------------

0060000 airin tmdpJun

0060101 007000000 005010003 1.1261-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

*0060203 150.0 0.00 7.00 0.0

*0060204 250.0 0.00 7.00 0.0

0060203 200.0 0.00 10.0 0.0

* -------------------------------------------------

0070000 airin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 l.+5 282.0 1.0

*---- n-------------------

0080000 tankfiow sngijun

0090000 tankl snglvol

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0

0090200 4 1.0+5 282.0 1.0

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0 *----------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

00

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

*a ---------- ------------ ------------ ------------

* ------------------

20500100 I-In mult 3.80-4 0.0 0

20500101 voldfj,002000000 velfj,002000000

20500200 1-down1 mult 1.1261-3 0.0 0

20500201 voidfj,00OAG0000 velfj,004,G0000 * -------------------- lei -------------- A.43 --------------

20500300 i-up mult 1.1261-3 0.0 0

20500301 voldfj,010000000 velfJ,010000000

20500400 k-f mult 6.1000 0.0 0

20500401 voldfj,0oo oo0000 veifJ,00o O3000

* -----------------

20502000 k-g-sqrt stdfnctn 1.0 0.0 0

20502001 sqrt cntrlvar,18

* ----------------

0000301 cntrlvar,2

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrivar,5

0000306 cntrnvar, a

0000307 cntrlvar,0

0000308 cntrivar,10

0000309 cntrvar,111

0000310 cntrlvar,19

0000311 cntrlvar,20 --------------------------------

* end of Input data

20500500 k-g mult 0.214 0.0 0

20500501 voldg J,00O OO00 velgj,00#O0000 *------------------304. .------------ .-t.........

20500600 diff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3050000 -1.0 P,3010000

*a----------------------------------------------------

20500700 k-fS mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4 * --------------

20500800 k-ga mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5 * --------------

20500900 k-f-sqrt stdfnctn 1.0 0.0 0

20500901 sQrt cntrlvar,7 a -----------------------

20501000 k-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8 ----------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrIvar,9 1.0 cntrlvar,10 * ---................----------------------

20501400 k-f muit 6.1000 0.0 0

20e01401 voidfJ,003010000 velfj,003010000 a ------------------------------------------------

20501500 k-g muit 0.214 0.0 0

20501501 voidgJ,003010000 velgj,003010000 *----------------------------------------------------

20501700 k-fs mult 1.0 0.0 0 1 0.0

20501701 cntrlvar,14

*a---------------------

20501800 k-gs mult -1.0 0.0 0 1 0.0

20501801 cntrlvar,15 *--------

20501900 k-f-sqrt stdfnctn 1.0 0.0 0

20501901 sqrt cntrlvar,17
0l

* ccfl model test dia= 2 cm , jg=O - 1.0

* ccfl junction = 4000000

* ccfl option = wallis type ( m=l.0 c=1.0

* fine nodallzation In the test section

000100 new transnt

0000102 sl sl

0000110 air

0000115 1.0 * ---------------------------------

0000201 . 0.00 1.0-7 0.01 3 100 1000 10000

0000202 200.0 1.0-7 0.01 3 100 1000 10000

* ---------------------------------

0010000 waterIn tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0

* ---------------------------------

0020000 waterin tmdpjun

0020101 001000000 003000000 0.00008

0020200 1

0020201 0. 0.0223 0. 0. * if = 0.16

0030000 section pipe

0030001 12

0030101 3.1416-4 12

0030301 0.15 12

0030601 -90.0 12

0030801 0.0 0.03 12

0031001 00000 12

0031101 00000 iI

0031201 4 1.+5 282.0 1.000 0.0 0.0 12

0031301 0.0 0.0 0.0 11

0040000 downfiow sngljun

0040101 003010000 005000000 3.1416-4 0.0 0.0 100000

0040110 0.02 0.0 1.0 1.0 *** waills flooding

0040201 1 0.0 0.0 0.0 * --------------------------------------------------

0050000 lowerpl snglvol

0050101 0.0348 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 00

0050200 4 1.0+5 282.0 1.0 * --------------------------------------------------

0060000 alrin tmdpjun

0060101 007000000 005010003 3.1416-4

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 15.00 0.0

a -Z ------------------------------------------------

0070000 alrin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 00

0070200 4

0070201 0.0 1.+5 282.0 1.0

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 1 0.0 0.0 0.0

..................................................

0090000 tanki snglvol

0090101 0,0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0 00

0090200 4 1.0+5 282.0 1.0

* --------------------------------------------------

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0

0110000 upPlenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

20500100 I-In mult 0.00008 0.0 0

20500101 voidfJ,002000000 velfJ,002000000

*------------------------------------

20500200 I-downl mult 3.1416-4 0.0 0

20500201 voidfJ,004000000 ve|fj,004000000

* -..................-.................................

20500300 I-up muit 3.1418-4 0.0 0

20500301 voidfj,010000000 velfJ,010000000

*------------------------------------

20500400 J-f mult 2.2590 0.0 0

20500401 voidfj,004000000 velfj,004000000

a -----------------------------------------------------

20500500 J-9 mult 0.07920 0.0 0

20500501 voidgj,004000000 velgj,004000000

a ------------------------------------------------

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 P,3120000 -1.0 p,3010000

*a------------------------------------------------

20500700 j-fs mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 J-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,5

a -------------------

20500900 J-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

* -----------------

20501000 J-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrivar,8

* ------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrivar,9 1.0 cntrlvarl0

*-------------------

20501200 massf mult 3183.00 0.0 0

20501201 mflowj,004000000

* -------------------

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

0000305 cntrlvar,S

0000308 cntrlvar,6

0000307 cntrivar,9

0000308 cntrivar,10

0000309 cntrlvar,1 1

0000310 cntrlvar,12

* ...........------------------------------------

20800001 voidj,3020000

20800002 voidj,3040000

20800003 voldJ,3060000

20800004 void J,3080000

20800005 voldJ,3f00000

20800006 void j,4000000

S-------------------------

20800007 flj,3020000

20800008 fij.3040000

20800009 flJ,3080000

20800010 fij,3080000

20800011 flj,3100000

20800012 fl),4000000

S-------------------------

20800013 Iregj,3020000

20800014 IregJ,3040000

20800015 IregJ,3080000

20800018 Iregj,3080000 20800017 iregJ,3OOOO 1 20800018 Iregj,4000000

S-------------------------

20800019 florgj,3020000

20800020 florgj,3040000

20800021 florgj,3080000

20800022 florgJ,300000

20800023 florgJ,3100000

20800024 florgj,4000000 a ---------------------------

20800025 vgjj,3020000

20800028 vgjj,3040000

20800027 vgjj,3060000

20800028 vgJj,4080000 20800029 vgjj,3100000

20800030 vgjj,4000000
S$$s$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$:

* ccfl model test bundle type I ( no spacer grid

* ccfI option junction = 4000000 * ccfl type = kutateladze form : constant from experiment

* file nodalization In the test section

000100 new transnt

0000102 si SI

0000110 air

0000115 1.0

* -------- 0------------------------------

0000201 10.00 1.0-7 0.01 3 100 1000 10000

0000202 300.0 1.0-7 0.01 3 100 1000 10000

0010000 waterln tmdpvol

0010101 1.0 0.1 0. 0. 0. 0. 0. 0. 00

0010200 3

0010201 0. 1.0+5 282.0 *----------------------------------

0020000 waterln tmdpJun

0020101 001000000 003000000 3.80-4

0020200 1

0020201 0. 0.040000 0. 0. * 2600 cc/min

0030000 section pipe

0030001 12

0030101 1.1261-3 12

0030301 0.125 12

0030601 -90.0 12

0030801 0.0 0.0103 12

0031001 00100 12

0031101 00000 11

0031201 4 1.+5 282.0 1.000 0.0 0.0 12

0031301 0.0 0.0 0.0 I1 * -----------------------------------

0090101 0.0346 0.5 0.0 0.0 -90.0 -0.5 0.0 0.0

0090200 4 1.0+5 282.0 1.0 * ---------------------------------

0100000 upflow sngljun

0100101 003000000 011000000 0.0 0.1 0.1 0000

0100201 1 0.0 0.0 0.0 * ---------------------------------

0110000 upplenum tmdpvol

0110101 0.125 0.44 0. 0. 0. 0. 0. 0. 00

0110200 4

0110201 0.0 1.0+5 282.0 1.0

G0 20800009 flJ,3060000

20800010 flJ,3080000

20800011 flj,3100000

20800012 fl1,4000000

20800013 florgj,3020000

20800014 florgj,3040000

20800015 florgJ,3060000

20800018 florgJ,3080o00

20800017 florgJ,3100000

20800018 florgJ,4000000

* end of Input data

U-,

20500100 I-In mult 3.80-4 0.0 0

20500101 voldfJ,002000000 velfj,002000000 * -----------------------------------

20500200 1-downl mult 1.1261-3 0.0 0

20500201 voldfJ,004000000 velfj,004000000 -----------------------------------

20500300 I-up mult 1.1261-3 0.0 0

20500301 voldfJ,010000000 velfJ,010000000

* -----------------------------------

20500400 k-f mult 6.1000 0.0 0

20500401 voldfj,004000000 velfj,004000000

20500500 k-g mult 0.214 0.0 0

20500501 voldgj,004000000 velgJ,004000000

20500600 dlff-p sum 1.0 0.0 0

20500601 0.0 1.0 p,3100000 -1.0 p,3020000

* -----------------------------------

20500700 k-fs mult 1.0 0.0 0 1 0.0

20500701 cntrlvar,4

* --------------

20500800 k-gs mult -1.0 0.0 0 1 0.0

20500801 cntrlvar,S *..........----------

20500900 k-f-sqrt stdfnctn 1.0 0.0 0

20500901 sqrt cntrlvar,7

20501000 k-g-sqrt stdfnctn 1.0 0.0 0

20501001 sqrt cntrlvar,8 *•---------------------------

20501100 curve sum 1.0 0.0 0

20501101 0.0 1.0 cntrlvar,9 1.0 cntrlvar,10

* ----------------

0000301 cntrlvar,1

0000302 cntrlvar,2

0000303 cntrlvar,3

0000304 cntrlvar,4

00(0305 cntrlvar,5

0000306 cntrlvar,6

0000307 cntrlvar,9

0000308 cntrlvar,10

0000309 cntrlvar,11I *--------------------------------------

0040000

0040101

0040110

0040201

downflow sngljun

003010000 005000000 1.1261-3 0.0 0.0 100000

0.0103 1.0 1.44 1.0 *** kutateladze form

1 0.0 0.0 0.0

0050000 lowerpl snglvol

0050101 0.0346 0.13 0.0 0.0 -90.0 -0.13 0.0 0.0 0

0050200 4 1.0+5 282.0 1.0 *.............................--------------------

0060000 airln tmdpjun

0060101 007000000 005010003 1.1261-3

0060200 0

0060201 0.0 0.0 0.0 0.0

0060202 10.0 0.0 0.0 0.0

0060203 200.0 0.00 10.0 0.0

0060204 300.0 0.00 20.0 0.0

* --------------------------------------------------

0070000 alrin tmdpvol

0070101 1.0 0.1 0.0 0.0 0.0 0.0 0.0 0.0 0

0070200 4

0070201 0.0 1.+5 282.0 1.0

* -----------------------------------

0080000 tankflow sngljun

0080101 005010000 009000000 0.0346 0.0 0.0 00000

0080201 I 0.0 0.0 0.0 * --------------------------------------------------

0090000 tankl snglvol

)0

O0

20800001

20800002

20800003

20800004

20800005

20800006

20800007

20800008

voldJ,3020000

voldJ,3040000

vold j,3080000

void j,3080000

voldJ,3100000

voldJ,4000000

flj,3020000

fi j,3040000

NRC FORM 335 U.S. NUCLEAR REGULATORY COMMISSION 1. REPORT NUMBER

.12-89) (A.i|ne~d by NRC. Add Vol.. Sum., Rev.,

NRCM. 1102. SHE AddeIndu Numbert, If sny,) 3201.3202 BIBLIOGRAPHIC DATA SHEET

ISO instruc-ons on the ,evesej NUREG/IA-0100

2. TITLE AND SUBTITLE

Assessment of CCFL Model of RELAP5/MOD3 3. DATE REPORT PUBLISHED

Against Simple Vertical Tubes and Rod Bundle Tests MONT --- YEAR

June 1993

4. FIN OR GRANT NUM8ER

L2245

5. AUTHORIS) 6. TYPE OF REPORT

S. Cho, N. Arne, Korea Electric Power Corp., Research Center Technical

B.D. Chung, H.J. Kim, Korea Instittue of Nuclear Safety 7. PERIODCOVERED ,,ncuusoa Cteu)

a. rn r tI I ~ ~U flU4P4u~s .'~up~iur,~ I i IVL. rvvoe vwo. u,,c ar egin. .~.,vuwarntggary ~mnuw1. M m,,,g aore: a cotrator OrvaI

a,.,,. ndnwdMkvadd~reXJ

Korea Electric Power Corporation, Research Center

Korea Institute of Nuclear Safety

117-15 Yongjeon-Dong, Dong-Gu, Taejon, Korea

P.O. Box 16, Daeduk-Danji, Taejon, Korea, 305-353

9. SPONSORING ORGANIZATION - NAME AND ADDRESS (It NRC. type "Sa,. jabovee'; if contractrt, provide NRC Division, Office or Region. UoS. NucIlea Regulatory Commi/on.

and maiMA addras& I

Office of Nuclear Regulatory Research

U.S. Nuclear Regulatory Commission

Washington, D.C. 20555

10. SUPPLEMENTARY NOTES

11. ABSTRACT 1200 wo,,' or iesi

The CCFL model used in RELAP5/MOD3 version 5m5 has been assessed against simple vertical tubes and rod

bundle tests performed at a facility of Korea Atomic Energy Research Institute.

The effect of changes in tube diameter and nodalization of tube section were investigated. The roles of interfacial

drags on the flooding characteristics are discussed. Differences between the calculation and the experiment are also

discussed.

A comparison between model assessment results and the test data showed that the calculated value lay well on the

experimental flooding curve specified by user, but the pressure jump before onset of flooding was not calculated.

12. KEY WORDS/DESCR!PTORS (List woros oap/ram; that wittist researAe" In Ifcadng the t.1 13. AVAI LABILITY STATEMENT

ICAP Program Unlimited

RELAP5/MOD3 14. SECURITY CLASSIFICATION

CCFL Model ,rIT Page

Unclassified

(Tihk Revorit

Unclassified

15. NUMBER OF PAGES

16. PRICE

NRC FORM 335 12491


Federal Recycling Program

UNITED STATES

NUCLEAR REGULATORY COMMISSION

WASHINGTON, D.C. 20555-0001

SPECIAL FOURTH-CLASS RATE

POSTAGE AND FEES PAID

USNRC

PERMIT NO. G-67

OFFICIAL BUSINESS

PENALTY FOR PRIVATE USE, $300