Cobra-Iv-I: Thermal-Hydraulic Nuclear: Cobra of Bundle Fuel and Cores
Cobra-Iv-I: Thermal-Hydraulic Nuclear: Cobra of Bundle Fuel and Cores
Cobra-Iv-I: Thermal-Hydraulic Nuclear: Cobra of Bundle Fuel and Cores
March 1976
.-.. ~~
C. W. Stewart
R. J. Cena
D. S. Rowe
A. M. Sutey
March 1976
Approved by
D. S, T r e n t ,'Managcr '
F l u i d and Thermal Engineering S e c t i o n
Engineering Technology Department
BATTELLE
P a c i f i c Northwest L a b o r a t o r i e s
Richland, Washington 99352
{h
The COBRA-IV-I computer code uses t h e subchannel a n a l y s i s approach t o
'determine t h e e n t h a l p y and f l o w d i s t r i b u t i o n i n r o d bundles' f o r both steady-
s t a t e and t r a n s i e n t c o n d i t i o n s . The s t e a d y - s t a t e and t r a n s i e n t s o l u t i o n
schemes used i n COBRA-IIIC a r e s t i l l a v a i l a b l e i n COBRA-IV-I as t h e i m p l i c i t
s o l u t i o n scheme o p t i o n . I n a d d i t i o n t o these techniques, a new e x p l i c i t
s o l u t i o n scheme i s now a v a i l a b l e which a l l o w s t h e c a l c u l a t i o n o f severe
t r a n s i e n t s i n v o l v i n g f l o w r e v e r s a l s , r e c i r c u l a t i o n s , e x p u l s i o n and r e e n t r y
flows, w i t h a pressure o r f l o w boundary c o n d i t i o n s p e c i f i e d . Significant
s t o r a g e compaction and reduced r u n n i n g times have been achieved t o a l l o w
t h e c a l c u l a t i o n o f problems i n v o l v i n g hundreds o f subchannels.
THIS PAGE
WAS INTENTIONALLY
LEFTBLANK
CONTENTS
ABSTRACT ............................... i
LIST OF FIGURES . . . . . . . . . . . . . . . . . . . . . . . . . . . . v
LIST OF TABLES . . . . . . . . . . . . . . . . . . . . . . . . . . . . v.
1.0 INTRODUCTION. . . . . . . . . . . . . . . . . . . . . . . . . . . 1
2.0 SUMMARY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
6.1 ORGANIZATION . . . . . . . . . . . . . . . . . . . . . . . . 29
iii
CONTENTS ( C o n t i n u e d )
8.0 .
NOMENCLATURE . .... . . ................... 77
APPENDIX D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . D-1
APPENDIX E . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . E-1
APPENDIX F . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . F-1
APPENDIX G . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . G-1
LIST OF FIGURES
Figure Page
Number Number
COBRA- I V - I F l ow Diagram
Regional A p p l i c a t i o n o f Heat T r a n s f e r C o r r e l a t i o n s
LIST OF TABLES
3 Storage L o e a t i o n o r EQUIVALENCE V a r i a b l e s
1.0 INTRODUCTION
COBRA-IV-I i n c l u d e s a l l t h e p r e v i o u s c a p a b i l i t y o f COBRA-IIIC p l u s
s i g n i f i c a n t l y expanded computational c a p a b i l i t y d e s c r i b e d below.
Continuity
ap
A - am + [DC] Tw = O
+ -
at ax
Energy
A x i a l Momentum
Transverse Momentum
* av w
aw
- au w [DC] P = c
dt ax ay
I n a d d i t i o n t o t h e c o n s e r v a t i o n equations, i t i s necessary t o s p e c i f y
f l u i d p r o p e r t i e s and c o n s t i t u t i v e equations t o form a c l o s e d s e t o f equa-
t i o n s f o r so!ution. I n COBRA, s p e c i f i c a t i o n s a r e r e q u i r e d f o r a x i a l and
l a t e r a l f r i c t i o n o r form losses, heat t r a n s f e r c o e f f i c i e n t s , c o e f f i c i e n t s
f o r t u r b u l e n t exchange, s l i p r a t i o , and an e q u a t i o n o f s t a t e f o r an
incompressible, b u t t h e r m a l l y expanding f l u i d . T h i s i.s w r i t t e n as:
where P* = r e f e r e n c e pressure.
T
'[Dc] and [ D C ] ~ r e p l a c e t h e n o t a t i o n [S] and [S] i n Reference 1.
TABLE 1. Major Assumptions and R e s t r i c t i o n s
Common t o b o t h numerical a l g o r i t h m s :
Energy
I.I No s o n i c p r o p o g a t i o n
Only mechanisms f o r i n t e r n a l
energy change a r e heat t r a n s f e r ,
convection, -and t u r b u l e n t
Neglect viscous d i s s i p a t i o n
N e g l e c t work a g a i n s t t h e g r a v i t y f i e l d
Gravitational acceleration i s the only
'body f o r c e
Energy
Energy
.
Momentum
I mixing
Electromagnetic f o r c e s n o t
considered
Transverse body f o r c e .neglected Transverse momentum No l a t e r a l buoyancy e f f e c t s
( r e s t r i c t s o n l y when flow
channel n o t v e r t i c a l )
Neglect surface tension contributions Momentum Requires equal phase pressures
F l u i d drag f o r c e s represented by w a l l Momentum F l u i d - f l u i d shear n e g l e c t e d
f r i c t i o n and form drag ( t u r b u l e n t momentum exchange i s
modeled, however)
~ e g l e c ttwo o f t h e s i x terms i n t h e ~ r a n s v e r s emomentum A1 lows a r b i t r a r y l a t e r a l
t r a n s v e r s e momentum e q u a t i o n (see coup1 i n g o f channels ~
Appendix C f o r c l a r i f i c a t i o n )
Flow area v a r i e s l i n e a r l y between No s t e p area changes
axial locations
For t h e e x p l i c i t method:
Homogeneous f l u i d Momentum and energy Equal phase v e l o c i t i e s , i.e.,
u = U,
9
vg = v1
Courant r e s t r i c t i o n Maximum t i m e s t e p l i m i t e d t n
At (Ax/u
For t h e i m p l i c i t method:
S p e c i f i e d a x i a l s l i p i s an o p t i o n ; Momentum and energy ugf U1' Vg = V1
however, t r a n s v e r s e s l i p i s n e g l e c t e d
No r e v e r s e ' f 1ow A1 1 Stability
3.2 IMPLICIT NUMERICAL SOLUTION
Since t h e f i n i t e d i f f e r e n c e f o r m u l a t i o n s o f COBRA-IV-I d i f f e r b u t
s l i g h t l y f r o m those g i v e n i n ~ e f e r e n c e1, t h e new f o r m u l a t i o n s a r e o n l y
s u m a r i zed here.
Continuity
Energy
A'; (hj - -h . ) m
j-1
( h . - hj-,) T T
At
J + Ax + [DC] h P j - hj [DC] w j
where
J J
Using equation (6) and the following definitions and identities,
where F . consists of the terms contained inside the { I , and R j contains the
J
terms inside the [I in equation ( 9 ) .
Equation (1 0) i s rearranged and written as
where r 1
and
't
Using equation (1 1 ) , the pressure difference between subchannel s at el eva-
tion j-1 can be written as .
Combining equation (15) wit'h (12) results in the set of simultaneous equa-
tions of the form
where
[mj] = [DC] [R~] Axt [RP]
and
{bj} = [DC] {Pj - F AX} - {FP}
Encrgy ( Impl i ci t )
A X I A L MOMENTUM
L A T E R A L MOMENTUM
2,
I
In t h i s case, E . represents the volume d i l a t i o n or volumetric e r r o r which
J
i s t o be corrected by pressure adjustment. The time derivative i s absent
because an Eulerian c e l l does n o t change i t s volume.
To form the pressure correction equation, the momentum equations (21)
and ( 2 2 ) can be substituted f o r the flows in equation (24). The derivative
of the r e s u l t with respect t o c e l l pressure i s
the d i l a t i o n , E i s s u f f i c i e n t l y small.
j '
T H I S PAGE
WASINTENTIONALLY
L E F T BLANK
4.0 THERMAL TRANSPORT
\
I n ^ a d d i t i o n t o t h e thermal t r a n s p o r t mechanisms a v a i 1a b l e i n COBRA-I II C Y
t h r e e new models have been added t o C O B R A - I V - I . These c o n s i s t o f :
a = 1 Planar
a = 2 Cylindrical
F i r s t , t h e r a d i a l c o o r d i n a t e i s approximated by t h e orthogonal c o l l o c a -
t i o n method o f weighted r e s i d u a l s . I n t h i s method a polynomial expansion
o f t h e dependent v a r i a b l e i s s u b s t i t u t e d i n t o t h e d i f f e r e n t i a l e q u a t i o n and
solved a t N c o l l o c a t i o n points. I n t h e p r e s e n t a p p l i c a t i o n , a symmetric
polynomial i s used:
or, i n m a t r i x notation,
where
and 1
1
which,, w r i t t e n i n m a t r i x n o t a t i o n w i t h { d l r e p l a c e d by [Q]- { e l , gives
v2{ e l = [ ~ l [ ~ l - ' { e l= [ B l i e l
where
dr2i-2
Cji -- dr /rj and Dji = V,2 ( r 2 i - Z j
Fuel I n t e r i o r ( i = 1 t o N-1)
Fuel Surface ( i = N)
- 1 - -+ -1 t c l ad
H h I(,
gap gap clad
Clad Surface ( i = N+1)
H CR*
+ S!JL-- (T/+I -; : IT 1 - -H s u r f n+l
( T ~ + l' T f l u i d 1
t c l ad . , tclad
where A denotes t h e p r e v i o u s i t e r a t i o n v a l u e , G i s d e f i n e d by e q u a t i o n ( 2 8 ) ,
and G ' i s t h e d e r i v a t i v e . o f G w i t h r e s p e c t t o T. The term c o n t a i n i n g 8 is
j
s u b s t i t u t e d on t h e i m p l i c i t s i d e o f t h e boundary equations w i t h t h e o t h e r
two terms e v a l u a t e d e x p l i c i t l y f r o m t h e p r e v i o u s i t e r a t i o n s o l u t i o n . Once
t h e 8 s o l u t i o n i s determined, an i n v e r s e t r a n s f o r m a t i o n i s performed t o
o b t a i n t h e t r u e temperature s o l u t i o n .
ROD RADIUS, R
CLADDING
hk = t h e subchannel s u r f a c e heat t r a n s f e r c o e f f i c i e n t as c a l c u l a t e d
C
by s u b r o u t i n e s PROP o r HTCOR. /
-Ayk
- - an i n p u t heat conduction parameter
kw
I n t h e i m p l i c i t scheme, e q u a t i o n ( 4 0 ) i s r e w r i t t e n i n t h e f o l l o w i n g
form
W I DTH OF WALL
CONNECT! NG
SUBCHANNELS
i AND j
25
where
{ b ) = v e c t o r o f f o r c i n g f u n c t i o n s composed of t h e terms on t h e r i g h t
s i d e o f e q u a t i o n s (42) and ( 7 ) .
S e l e c t i o n o f t h e a p p r o p r i a t e h e a t t r a n s f e r c o r r e l a t i o n i s based on,
several c r i t e r i a . The f i r s t i s an e v a l u a t i o n o f .the p h y s i c a l s t a t e o f t h e
f l u i d t o determine which regime o f heat t r a n s f e r i s o c c u r r i n g (e.g.,
forced convection, f i l m b o i l i n g , e t c . ) . The second i s t h e e v a l u a t i o n o f
t h e c r i t i c a l h e a t f l u x (CHF) and t h e c l a d temperature t h a t would e x i s t a t
t h i s h e a t f l u x t o determine whether t h e pre-CHF o r p o s t - C H F ' c o r r e l a t i o n s
are appropriate. The t h i r d c r i t e r i o n i s t h e comparison o f two o r more
calculated heat transfer c o e f f i c i e n t s t o s e l e c t the proper c o e f f i c i e n t i n
r e g i o n s where c o r r e l a t i o n s o v e r l a p .
Another a d d i t i o n t o COBRA i s t h e i n c l u s i o n o f a r o u t i n e t o c a l c u l a t e
t h e p r o p e r t i e s of superheated steam t o a l l o w c o n ~ p u t a t i o n so v e r d wider
range o f condi t i o n s .
6.1 ORGANIZATION
. a. The r e s u l t s o f t h e s t e a d y - s t a t e c a l c u l a t i o n s a r e p r i n t e d by s u b r o u t i n e
/
schemes.
G (READ INPUT)
TRANS lENT
replacement o f - m a n y e x t e r n a l s u b r o u t i n e c a l l s w i t h more e f f i c i e n t
i n t e r n a l statement f u n c t i o n s ;
replacement o f t h e i m p l i c i t m a t r i x c r o s s f l o w s o l u t i o n w i t h an
i t e r a t i v e s o l uti'on;
1 oops.
The layout of the common storage area LARGE 1 , which contains most of
the variables discussed in Table 2, is shown in Table 3.
TABLE 2. EQUIVALENCE R e l a t i o n s h i p s
1 DR(1) )I(' !
-
Dimensioned by Rods
MR 7
HINLET(L) VPA(1)
Oimensioned by Channels
FINLET(1) VISC(1)
RC = 18
TItILET(1) VISW(1)
'U(1.1) SAVEAL(1)
UOLD(1 , I ) SAFEAL(Z5) BLOCK 2
BETA(1 .I) SAFEAL(49) ' A x i a l l y Dependent
BSAVE(1 . I ) SAVEAL(73) Variables
.................................. o.----.---
imensionea by Rods
---- --------
TROD(1 .I, I ) SAVEAL(265) and Fuel Nodes (MN = 6 )
.-.-.-------------------------------
I6~~~!s1~1~!0!---~~---
SAVEAL(373)
SAVCiIL( 331 )
SAVEAL( 409)
SAVEAL(427)
SAVEAL(445)
S A V E A L ( J ~ ~ ) Oimensioned by Channels
SAVEAL(48I ) MC = 18
SAVEAL( 499)
SAVEAL(517)
SAVEAL(535)
SAVEAL(553)
SAVEAL(571)
SAVFAI.(SRq)
BETA(l.1)
BSAVE(l.1)
PSI(1.1)
PSI(1.1)
TRF(l.1)
CHFR(1 .I)
P(1.1)
STUW(1 ,I)
BSAVE(l.2)
STUW(1.2)
BSAVE(l.3) BLOCK 3
3Tud(l,3) Tal~uur.ary Vat ~alrlc Stoi.ouc
VISC(1)
VISCH(1) (In order as variables are
vrzcwll~ 5torcd in GOMMOtl, see
Table 3)
CON(1)
CON(1)
CP( I J
CP(1)
FSP(1)
FSP(1)
FSP(1)
FSP(1)
FSP(1)
WP(1)
UP(1)
WP(1)
WP(1)
dAA(1.l)
Another s i g n i f i c a n t compaction o f s t o r a g e i s t h e r e d u c t i o n o f t h e
c r o s s f l o w c o e f f i c i e n t m a t r i x AAA(K,K) t o an e q u i v a l e n t AAA(L,K) matrix,
where K i s t h e number of gap connections and L i s e i t h e r 7 f o r square
&. a r r a y s o r 5 f o r hexagonal a r r a y s . . This was made p o s s i b l e by r e p l a c i n g t h e
m a t r i x ~ s o l u t i o na l g o r i t h m w i t h an i t e r a t i v e s o l u t i o n technique. Any row o f
AAA has no more than 7 nonzero elements. The diagonal elements AAA(K, K)
a r e t r a n s f e r r e d t o AAA(1,l) and succeeding nonzero v a l ues f i11 i n rows 2
through 7 by u s i n g t h e gap l o c a t i o n a r r a y LOCA.
6.2.3 E x t e r n a l Devices
The c e n t r a l memory r e q u i r e m e n t f o r t h e p u b l i s h e d v e r s i o n i s
ZOZ,0008 words. I n addition, peripheral devices are required f o r input,
o u t p u t and f o u r s c r a t c h u n i t s . The l a r g e s t v e r s i o n o f t h e code executed
t o d a t e was dimensioned f o r a f u l l 217-pin wire-wrapped b u n d l e r e q u i r i n g
2 7 5 ~ 0 0 0words
~ o f c e n t r a l memory. The code i s a l s o s e t up t o u t i l i z e l a r g e
- c o r e memory (LCM) on CDC-7600 machines so t h a t c e n t r a l memory r e q u i r e m e n t s
do n o t r e p r e s e n t a p r a c t i c a l l i m i t t o problem s i z e .
THIS PAGE
WAS INTENTIONALLY
LEFT BLANK
7.0 COBRA INPUT DESCRIPTION .
The i n p u t i n s t r u c t i o n s i n t h i s s e c t i o n a r e ordered as t h e i n p u t i s
encountered by t h e code. An i d e n t i f y i n g l a b e l i s assigned t o each d a t a c a r d
o r c a r d t y p e ( i f m u l t i p l e cards r e a d ) . A l i s t o f t h e v a r i a b l e s and t h e
format under which t h e y a r e read f o l l o w s . The l a b e l i s composed o f t h e
name o f t h e s u b r o u t i n e where t h e read takes p l a c e and t h e r e l a t i v e o r d e r o f
occurrence i n t h a t subroutine, (e.g., COBRA.l). 0ptional:data are included
i n t h e o r d e r t h e y a r e encountered, w i t h t h e c r i t e r i o n f o r i n c l u s i o n i n t o t h e
d a t a deck s t a t e d i n t h e d e s c r i p t i o n .
Several u s e r - s p e c i f i e d e m p i r i c a l c o r r e l a t i o n s a r e i n p u t t o COBRA. To
c l a r i f y which c o r r e l a t i o n i s desired, t h e e m p i r i c a l c o r r e l a t i o n s a r e des-
c r i b e d i n Appendix B. To f u r t h e r a s s i s t t h e user i n s e t t i n g up t h e i n p u t ,
many o f t h e c a l c u l a t i o n parameters have been s e t t o a d e f a u l t v a l u e which
w i l l be used by t h e code if t h e u s e r leaves t h e parameter i n p u t f i e l d b l a n k .
The parameters which have d e f a u l t e d values a r e s p e c i f i e d i n t h e i n p u t i n s t r u c -
t i o n s and t h e d e f a u l t v a l u e i s given.
., Since most parameters w i t h d e f a u l t
values a f f e c t t h e " s t a b i l i t y " o f t h e s o l u t i o n , i t i s recommended t h a t t h e
d e f a u l t v a l u e be used u n t i l f a m i l i a r i t y w i t h each parameter i s reached. As
an a i d t o t h e user, a d e s c r i p t i o n o f t h e i m p o r t a n t code v a r i a b l e s i s pre-
sented i n Appendix F. F i n a l l y , a number o f sample problems have been
i n c l u d e d i n Appendix G showing many o f t h e o p t i o n s d e s c r i b e d i n t h i s r e p o r t .
*
I t i s hoped t h a t t h e i n f o r m a t i o n p r o v i d e d w i l l a i d t h e u s e r i n u t i l i z i n g
COBRA-IV-I t o i t s f u l l e s t potential.
INPUT INSTRUCTIONS
Must be t h e f i r s t d a t a c a r d o f t h e i n p u t
deck. Read o n l y once.
Optional i n p u t - MAXT n e g a t i v e
R e s t a r t o p t i o n s , where:
NA = Number o f a d d i t i o n a l i t e r a t i o n s f o r
i m p l i c i t s o l u t i o n scheme. May be r e s e t
i n setup i f NJUMP = 1.
= 1, b e g i n a t r a n s i e n t s o l u t i o n a t time
zero from a p r e v i o u s s t e a d y - s t a t e s o l u -
t i o n , must r e a d i n t r a n s i e n t d a t a from
setup. For o t h e r r e s t a r t cases, I T = 0.
= P r i n t oflion f o r i n p u t data. J1 0,
p r i n t o n l y new i n p u t data. J1 = 1,
p r i n t a l l i n p u t data. J1 = 2 , p r i n t
only operating conditions. J1 = 10,
p r i n t , a l l i n p u t data, then s t o p .
DATA CARD OR GROUP CONTROL CARD ( C o n t i n u e d )
TEXT = O u t p u t t e x t f o r problem i d e n t i f i c a t i o n ,
maximum 68 c h a r a c t e r s .
= NPROP, number o f p r o p e r t y c a r d s t o be
read.
PP = system p r e s s u r e (psia)
TT = Temperature (OF)
VV F = L i q u i d s p e c i f i c volume (ft3/lb)
-
SETUP.1.2 DTMAX Format (F5.3)
Optional Input - N2 = 1 and N3 = 1
44
DATA CARD OR GROUP CONTROL CARD (Continued)
= NVISCW. Wall V i s c o s i t y C o r r e l a t i o n .
N4 = 0, no heated w a l l c o r r e l a t i o n i n
f r i c t i o n factor. N4 = 1, ' i n c l u d e
heated w a l l c o r r e l a t i o n i n f r i c ' t i o n
factor.
= LAMNF, l a m i n a r f r i c t i o n f a c t o r .
N5 = 0, no l a m i n a r c o r r e l a t i o n .
N 5 = 1, read i n up t o 4 s e t s o f
laminar f r i c t i o n f a c t o r c o r r e l a t i o n
constants.
Optional i n p u t - N5 = 1.
Laminar f r i c t i o n f a c t o r c o e f f i c i e n t s .
-
SETUP.2.4
'
NV, AV(1)
. ..
( I = 1 t o 7)
Optional I n p u t NZ > 4.
Format (15, 7E10.5)
S l ip r a t i o and v o i d f r a c t i o n c o r r e l a t i o n , where:
= polynomial ~ r d e rfor' l t ~ ev u i d f r a c t i o n
c o r r c l s t l o n 5 f N2 = 6. I f N2 = . 5 , NV
i s n o t used.
= I f N2 = 6, AV(1) . i s J. s i n g l c term
s l i p ratio. I f N2 = 6, read A V ( I ) ,
I = 1 t o NV c o e f f i c i e n t s f o r up t o
6 t h o r d e r polynomial f u n c t i o n o f v o i d
f r a c t i o n versus steam q u a l i t y .
I
-
SETUP.2.5 NF, AF(I)
Optional Input
( I = 1 t o 7)
- N3 > 4
Format ( 1 5 , . ' 7 ~ 1 0 . 5 )
AF = Constants f o r up t o 6 t h o r d e r p o l y -
nomial f u n c t i o n o f t h e two-phase .-,
Y = R e l a t i v e p o s i t i o n (X/L) a t which h e a t
f l u x i s given, where L i s t h e t o t a l
l e n g t h o f t h e bundle. Must i n c l u d e
0.0 and 1.0 as end p o i n t s .
= T o t a l number o f subchannels, r e g a r d l e s s
o f N1.
= O p t i o n t o s p e c i f y N3 r e l a t i v e d i r e c t i o n
angles f o r c r o s s f l o w between a d j a c e n t
subchannel s.
= L o g i c a l u n i t from which c a r d s e t
SETUP.4.1 i s t o be read. N5 = 0,
r e a d from i n p u t deck. N5 > 0, r e a d
cards from l o g i c a l u n i t N5. Recommend
N5 = 9 i f cards a r e n o t s p e c i f i e d .
= L o g i c a l u n i t from which c a r d s e t
SETUP.4.3 i s t o be read. N6 = 0,
r e a d from i n p u t deck. N6 > 0, r e a d
cards from l o g i c a l u n i t N6. Recommend
N6 = 9 if cards a r e n o t s p e c i f i e d ,
= Suhchannel i d e n t i f i c a t i o n number.
2
= Nominal subchannel area ( i n ) .
= Nominal subchannel w e t t e d p e r i m e t e r
(in).
DATA CARD OR GROUP CONTROL CARD (Continued)
= Adjacent subchannel i d e n t i f i c a t i o n
number, f o r up t o 4 subchannels
a d j a c e n t t o subchannel I. I f t h e
subchannels a r e i n p u t w i t h ascending
i d e n t i f i c a t i o n numbers, o n l y connec-
t i o n s LC(1,L) > I need be read i n .
DIST = C e n t r o i d - t o - c e n t r o i d d i s t a n c e between
t h e a d j a c e n t subchannels s p e c i f i e d by
LC ( i n ) . ~ e q u i r e do n l y i f N1 = 3
i n Group 10 (SETUP.lO).
Optional I n p u t - N4 > 0
- . .
?See Appendix B f o r d e f i n i t i o n s .
DATA CARD OR GROUP CONTROL CARD (Continued)
RWALL(1 ,KW),
i
= NAXI.,, n~.rrnb~r
o f axial locations for
subchannel a r e a v a r i a t i o n .
= NARAMP, t h e number o f i t e r a t i o n s f o r
gradual i n s e r t i o n o f area v a r i a t i o n s .
T f b l a n k o r x r n , NARAMP - 1. For a
" r e s t a r t " case, NARAMP must be r e r e a d
i f desired.
DATA CARD OR GROUP CONTROL CARD (Continued)
= I d e n t i f i c a t i o n Number o f a subchannel
f o r wh.ich area v a r i a t i o n s a r e being
specified. Read a c c o r d i n g t o (15)
format, then s k i p t o t h e n e x t c a r d and
read a complete s e t o f f a c t o r s (AFACT)
corresponding t o t h e AXL l o c a t i o n s .
. Repeat u n t i . 1 f a c t o r s f o r N1 subchannels
a r e read.
= NGXL, number o f a x i a l l o c a t i o n s f o r
GAP v a r i a t i o n .
DATA CARD OR GROUP CONTROL CARD ( C o n t i n u e d )
T a b l e o f a x i a l l o c a t i o n s , where:
F o r N1 GAPS r e a d GAP v a r i a t i o n s a t N2 a x i a l l o c a t i o n s
(GAPXL) , where:
,
= GAP i d e n t i f i c a t i o n number o f GAP t o be
varied. Read K [ F o r m a t ( I 5 ) ] , then s k i p t o
n e x t c a r d and r e a d N2 GAP v a r i a t i o n
factors. Repeat u n t i l f a c t o r s f o r N1
GAPS a r e read.
= 56. N1 = 1 s p e c i f i e s w i r e wrap
input only. N1 = 2 s p e c i f i e s g r i d
spacer i n p u t o n l y . N1 = 3 s p e c i f i e s
b o t h w i r e wrap ir ~ p u Ldrld g r i d spacer-
loss coefficient input.
*
g r i d spacers.
= L o g i c a l u n i t t h a t c a r d s e t s SETUP.7.2
and SETUP.7.3 a r e t o be read. N6 = 0,
r e a d from i n p u t deck. N6 > 0, r e a d
. .from l o g i c a l u n i t N6. Recommend
N6 = 9 i f cards a r e n o t s p e c i f i e d .
\
, -
SETUP. 7.1 PITCH, DIA., THICK
Optional Input -
' Format (3E10.5)
N1 = 1 o r 3
-
. ,
Optional Input - N1 = 1 o r 3 ,
K = GAP number.
DATA CARD OR GROUP CONTROL CARD (Continued)
= The e f f e c t i v e f r a c t i o n o f a p i t c h l e n g t h
f o r f o r c i n g crossflow. Recommended
value i s d e r i v e d from t h e equation '
DUR = AXIPITCH
where
AX = t h e a x i a l node l e n g t h ( i n ) .
PITCH = t h e w i r e wrap p l tcti ( i l l ) .
Optional I n p u t - N1 = 1 . o r 3
Wrap i n v e n t o r y where:
channels a r e read.
DATA CARD OR GROUP CONTROL CARD (Continued)
-
SETUP.7.4 GRIDXL(I)., IGRID(1)
Optional I n p u t -
( I = 1 t o N3)
N1 = 2 o r 3
Format [ 6 ( ~ 5 . 2 . 1 5 ) ]
A x i a l l o c a t i o n of g r i d s and t t i e t y p e o f g r i d a t each a x i a l
location
IGRID = The g r i d t y p e a t a x i a l l o c a t i o n
, . GRIDXL.
-
SETUP.^.^ . r ( J y CD(Jy1)y Ky FXFLOW(K,I)
Format (15, E5.2; 15, E5.2)
Optional I n p u t - N1 = 2 o r 3
, I 1 = 1 t o NCHNL), I = 1 t o N41
Read t h e l o s s c o e f f i c i e n t and f o r c e d c r o s s f l o w , i f
desired, i n each subchannel f o r eash g r i d type.
f o r a1 1 subchannel s f o r g r i d t y p e 1 ,' t h e n
r e p e a t f o r succeeding g r i d types.
= 1 d e n t i f i c a . t i o n number o f subcha.nne1 f o r
which d a t a a r e being s u p p l i e d .
= LOSS c o e f f i c i e n t i r ~subchanr~el J f o r
g r i d type I .
= GAP a s s o c i a t e d w i t h subchannel Jy
through which f o r c e d d i v e r s i o n cross-
flow i s specified. I f b l a n k , no
forced crossflow.
FXFLOW = F r a c t i o n o f a x i a l f l o w which i s
d i v e r t e d through GAP K f o r g r i d t y p e I .
I f b l a n k , no f o r c e d c r o s s f l o w .
DATA CARD OR GROUP CONTROL CARD ( C o n t i n u e d )
= The number o f c a r d s of r o d d a t a t o be
read. One c a r d f o r each r o d modeled.
For RESTART cases o r m u l t i p l e case
r u n s , i t i s o n l y necessary t o r e a d
new r o d i n p u t .
= NROD, t h e t o t a l number o f r o d s t o be
model ed r e g a r d l e s s o f N1 .
= NCy o r d e r of a p p r o x i m a t i o n used i n
f u e l mode.1 . N3 = 0, no f u e l model.
N3 = 2 , 2nd o r d e r c o l l o c a t i o n s o l u t i o n .
N3 = 3, 3 r d o r d e r c o l l o c a t i o n . N3 > 3
o r N3.= 1 i s i l l e g a l .
= NFUELT, t h e number o f f u e l m a t e r l a l s
f o r w h i c h thermal p r o p e r t . i e s a r e t o
he s p e c i f i e d . Nnt a p p l i c a b l e i f N3 = 0.
I f b l a n k , NFUELT = 1 I s assigned.
= NCHFy c r i t i c a l h e a t f l u x o p t i o n .
NS = O, No ChF c a l c u l a t i o n s a r e
performed. N5 = 1, BAN-2 c o r r e l a t i o n
i s used. N5 = 2 , t h e W-3 c o r r e l a t i o n
i s used.
56
DATA CARD OR GROUP CONTROL CARD-(Continued)
= NHTC , h e a t t r a n s f e r c o r r e l a t i o n o p t i o n .
N7 = 0, D i t t u s and B o e l t e r o r o t h e r
c o r r e l a t i o n s p e c i f i e d i n SETUP.2.2
i n p u t i s used throughout. N7 = 1,
use a complete b o i 1ing and nonboi 1 ing
c o r r e l a t i o n package s i m i l a r t o t h a t
o f t h e RELAP-4 computer code.
= NRODTP, o p t i o n f o r a x i a l l y v a r y i n g
fuel material. N8 = 0, each fuel r o d
i s constructed o f a single material
and no a x i a l l y v a r y i n g d a t a a r e read.
N8 > 0, must r e a d f u e l zone . i n f o r m a t i o n
.(SETUP.8.4), f o r each f u e l t y p e N.
= L o g i c a l u n i t from which c a r d s e t
SETUP.8.1 i s t o be read. N9 = 0, r e a d
from i n p u t deck. N9 > 0, r e a d from
l o g i c a l u n i t N9. Recommend N9 = 9
when n o t r e a d i n g from cards.
For N8 = 0 ( A x i a l l y uniform f u e l )
ABS (N) corresponds t o one o f N4 mate-
r i a l s (sETUP.8.2) o f which Rod I
i s made.
= Rod i d e n t i f i c a t i o n number.
.A,
Card Label variables Format and E x p l a n a t i o n
CCLAD = S p e c i f i c h e a t o f c l a d (Btu/lb-OF).
A.
Card Label . Variables Format and Explanation
IZTYP = Type of material in fuel zone ending
at ZEND. Each IZTYP corresponds to a
material specifTed in SETUP.8.2. Only
'
= 55, an o p t i o n f o r r a d i a l thermal
conduction mixing. N3 = 0, no thermal
conduction. N3 = 1, read thermal
geometry f a c t o r , GK f o r r a d i a l thermal
conduction.
Optional I n p u t - N2 > 1
SETUP. 10.3 -
GK Format (F5.3)
Optional I n p u t - N3 = 1
= IH, o p t i o n f o r s p e c i f i e d i n l e t enthalpy
o r temperature. N1 = 0, H I N (SETUP.ll. 1 )
i s t h e i n l e t enthalpy. N1 = 1, H I N i s
t h e i n l e t temperature. N1 = 2', read i n
an i n l e t enthalpy f o r each subchannel.
N1 = 3, read i n an i n l e t temperature f o r
each subchannel.
68
DATA CARD OR GROUP. CONTROL CARD (Continued)
69
DATA CARD OR GROUP CONTROL CARD (Continued)
i
= Transient f o r c i n g f u n c t i o n f o r e x i t
enthalpy. Read i n NHX p a i r s o f values
f o r a t a b l e o f e x i t enthalpy f a c t o r
versus time.
PEX I T = system
T ~ Q pr~sqllre, (psis).
' . -
SETUP.11.2 '. HINLET(1)
Optional I n p u t
( I = 1 t o NCHANL)
N1 = 2 o r 3
Format (12E5.0)
.-.
subchannel .
-
SETUP.11.3 FINLET(1)
Optional I n p u t
( I = = 1 TO NCHNL)
- N2 > 2
Format (12E5.0)
I n l e t f l o w r a t e , where:
S E T J J P . ~.4
~ YP(I), FP(1) ( I - = 1 t o N3.) Format (12E5.0)
_____t
Optional I n p u t - N3 > 1
Pressure t r a n s i e n t t a b l e , where:
,
,
-
SETUP.11.5 YH(I), FH(1)
. Optional I n p u t
( I = 1 TO N4)
N4 > 1
Format (12E5.0)
- . . .....
?
* .
or t h e f r a c t i o n o f i n l e t temperature
(N1 = 1 o r 3 ) a t t h e t r a n s i e n t time
(YH).
71
DATA CARD OR GROUP CONTROL CARD (Continued)
-
SETUP.11.6 YG(I), FG(1)
Optional I n p u t
( I = 1 t o N5)
- N5 > 1
Format (12E5.0)
- .
p s i (N7
(YG).
-
(N7 = 1) o r FG i s t h e pressure drop i n
0) a t t h e t r a n s i e n t time
-
SETUP.11.7 YQ(I), FQ(1)
Optional I n p u t
(I= 1 t o N6)
N6 > 1
Format ( 1 2 ~ 5 . 0 )
-. .. .
Heat f l u x t r a n s i e n t table.
Optional I n p u t N8 > 1
. E x i t enthalpy t r a n s i e n t t a b l e .
N1 = NOUT, p r i n t option. N1 = 0, p r i n t
subchannel data o n l y . N1 = 1, p r i n t
subchannel data and crossflows o n l y .
N1 = 2, p r i n t subchannel data and f u e l
r o d temperatures only. N1 = 3, p r i n t
subchannel data, c r o s s f l ows and f u e l
r o d temperatures.
= NPROD, an o p t i o n f o r f u e l r o d heat f l u x
and/or temperature p r i n t o u t . N3 = 0,
data f o r a l l rods a r e p r i n t e d i f c a l l e d
f o r by N1. N3 > 0, read i n N3 r o d
i d e n t i f i c a t i o n numbers o f rods t o be
printed. I f NCHF (SETUP.8) i s >O, CHF
data i s a l s o p r i n t e d along w i t h t h e r o d
data.
DATA CARD OR GROUP CONTROL CARD (Continued)
printed. N2 = 0, p r i n t . a l 1 subchannels.
HIP(1). 2
) Mass flux (million lb/hr-ft )
ZPT(1)
. ,
) Enthal py (1 000 Btu/l b)
ZPT(2)
HIP(3)
) Pressure Drop (psi)
ZPT(3)
HIP(4)
} , Crossflow (1 b/sec-ft)
ZPT(~)
Mass flow rate, enthalpy and pressure
drop are plotted for each subchannel
specified (see N6 and N7). Crossflow
1s plotted Tor each GAP specified
(see NG and NS). If HIP(1) = 0, that
quantity is not plotted.
8. o NOMENCLATURE+'
Wette,d perimeter, ( f t )
P r a n d t l number
Heat r a t e , ( ~ t u l s e c - f t ' )
Reynol ds number
Radi a1 coordinate
Temperature, (OF)
Time, (sec)
A x i a l subchannel ~ e l o c i ~ t y(,f t l s e c )
M i x t u r e s p e c i f i c volume, ( f t 3 / 1b )
Q u a l i t y , mg/(mg + mf) ,
A x i a l node length, ( f t )
AY Width of thermal ly-conducting wall, ( f t )
u Viscosity, (1bf/ft-sec)
Superscripts
*' Convected quantity
- Previous time quantity
Axial levels
REFERENCES
SUBROUTINE DESCRIPTION
Subroutine BVOID ( I )
PROGRAM COBRA
FX = Q u a n t i t y t o ' be Tound
X =. Independent v a r i a b l e
F = I n p u t a r r a y o f , values of t h e dependent v a r i a b l e
Y
= I n p u t a r r a y of values of t h e independent v a r i a b l e
=.Number o f F values i n t a b l e
N
J = Error signal
ISAVE = Table search switch. For ISAVE = 1, a complete t a b l e search
on t h e independent v a r i a b l e i s done. For ISAVE = 2, t h e
l o c a t i o n i n t h e t a b l e which brackets t h e independent v a r i a b l e
i s known from a previous c a l l t o curve and t h e t a b l e search i s
n o t performed..
Subroutine 'DIFFER (IPART, J, JX)
Subroutine DIFFER i s d i v i d e d i n t o 4 p a r t s as i n d i c a t e d by t h e v a r i a b l e
IPART.
I - channel number
N - .rod number
J - calculation level
JX - axial location
HTC - heat t r a n s f e r c o e f f i c i e n t c a l c u l a t e d i n HTCOR
1MODE.- f l a g t.0 designate regions on t h e b o i l i n g curve
NSS -
conduction model. NSS i n d i c a t e s from which r o u t i n e HEAT i s c a l l e d .
0, SCHEME, NSS = 1, XSCIIEM. NSS =. 2 , i n i t i a l i z a t i o n i n COBRA.
JB = t h e c a l c u l a t i o n l e v e l and JX3 i s t h e a x i a l l o c a t i o n .
HTCOR c a l c u l a t e s a heat t r a n s f e r c o e f f i c i e n t by t h e c o r r e l a t i o n
specified i n the c a l l l i s t .
The arguments i n t h e c a l l l i s t are:
HTC - heat t r a n s f e r c o e f f i c i e n t
TSURF - r o d surface temperature
I - channel number
N - r o d number
J - calculation level
RE - Reynolds number
DE - h y d r a u l i c diameter
IPART - f l a g - t o designate a p a r t i c u l a r ' c o r r e l a t i o n
Subroutine ISWAP ( I , J, K l
Subroutine LOAD (X, Y, Z, MIN, MAX, LIMIT, STEP, IMAGE, CARD, LU) I
X = F i r s t v a r i a b l e t o be loaded. (Step = 1 )
Y = Second v a r i a b l e t o be loaded. (Step = 2 ) .
Z = T h i r d v a r i a b l e t o be loaded. (Step = 3 )
STEP = The number o f v a r i a b l e s t o be loaded s e q u e n t i a l l y . (1 cSTEP53)
MIN = The maximum allowable s e t s o f (X,Y ,Z) which can be loaxed.
MAX = I d e n t i f i e s t h e number of s e t s of (X,Y ,Z) t h e user attempts
t o load.
I f (MAX>MIN), then a dummy read i s used t o account f o r remaining
cards. This allows subsequent i n p u t t o be edited.
LIMIT = I d e n t i f i e s t h e maximum number o f data values per card which can
be loaded v i a t h e forrnattcd IMAGE.
IMAGE = I d e n t i f i e s t h e v a r i a b l e FORMAT t o be used f o r l o a d i n g i n p u t
data.
CARD = I n p u t card counter; f o r l a b e l i n g which card may be i n e r r o r .
LU = I n p u t device used f o r reading i n p u t data.
Subroutine LOADL (X, Y, Z, MIN, MAX, LIMIT,. STEP, IMAGE, CARD, LU)
LOADL i s i d e n t i c a l i n f u n c t i o n t o subroutine LOAD, except t h a t (X,Y,Z)
a r e v a r i a b l e s which may be a l l o c a t e d t o Level 2, ( l a r g e core memory) on a
CYBER-7600 system. This r o u t i n e i s used regardless o f t h e system on which
COBRA i s being r u n and t h e necessary storage a l l o c a t i o n s a r e provided by
Program SPECSET.
Subroutine LIMITS (NUM, MIN, MAX, GROUP, CARD, ERROR)
This r o u t i n e i s designed f o r use w i t h i n subroutine SETUP t o enhance
t h e e d i t i n g of i n p u t data. I t s f u n c t i o n i s t o guarantee t h a t t h e number
o f values read i n t o an a r r a y are w i t h i n t h e dimensioned l i m i t s o f t h e
code. I f not, t h e parameter, NUM, i s changed o n l y f o r t h e convenience o f
e d i t i n g ; however, t h e case w i l l be terminated a f t e r t h e i n p u t i s edited.
I n a d d i t i o n , a - l i m i t o f t w e n t y - f i v e accumulated e r r o r s are allowed, a f t e r
which t h e e d i t i n g w i l l cease.
Subroutine ROLLIT (J, JX, SAVEA1, SAVEA2, SAVEA3, NWR, NDXP1, LUO, LUI)
The c a l c u l a t i o n a l procedure i s b e s t i l l u s t r a t e d by t h e f l o w c h a r t
shown i n Figure A-2. F i r s t t h e enthal py h ( x ) i s c a l c u l a t e d and an estimate
f o r m(x) i s made f o r t h e f i r s t i t e r a t i o n only. The crossflow w(x) and
a x i a l f l o w m(x) a r e then calculated. Any value of w(x) o r m(x) n o t con-
verged t o w i t h i n a selected t o l e r a n c e i s s u f f i c i e n t t o s e t t h e nonconverg-
ence f l a g JUMP t o 1. Following t h e flow and crossflow c a l c u l a t i o n s , t h e
pressure P(x) and pressure d i f f e r e n c e [D]P(x) i s calculated. The procedure
i s repeated stepwise up t h e bundle u n t i l t h e e x i t i s reached. I f t h e flows
and crossflows have converged (JUMP = 2) o r t h e maximum allowable number
of i t e r a t i o n s i s reached, c o n t r o l i s r e t u r n e d t o t h e main program. Other-
P
wise t h e bundle i s swept again using new estimates f o r flows i n t h e
calculations.
FROM SUBROUTINE SETUP
DO NOT READ ANY FOLLOWING RETURN FROM
ADDITIONAL INPUT
BEGIN A TRANSIENT
( NJUMP = 2 \r FROM T l M E ZERO
FROM A PREVIOUS
STEADY STATE SOLUTION
CALCULATE A NEW
AS FIRST GUESS
I CONTINUE
lTERATI ONS ON
PREV lOUS STEADY
STATE SOLUTION
/1 CONTINUE
TIME STEPS
ON PREVIOUS
TRANSIENT SOLUTl ON
11
I CONTINUE TI ME STEPS
ON A PREVIOUS TRANSIENT
SOLUTION. WITH ADDITIONAL
DATAREAD 11
\
1 CALCULATE
PARAMETERS
ATx=O 1
SET NEW VALUE OF x
CALCULATE h(x)
CALCULATE m(x)
- -
END OF CHANNEL
1 YES
A- 9
Regardless o f convergence, a minimum number o f i t e r a t i o n s i s r e q u i r e d
depending on t h e f u e l model options employed:
-- . - . . -.
M i nimum Number \
2 No f u e l model
JUMP- Convergence i n d i c a t o r
Subroutlne SETUP
Subroutine V O I D ,
Subroutine V O I D c a l c u l a t e s t h e subcooled v o i d f r a c t i o n , b u l k v o i d
f r a c t i o n , d e n s i t y , e f f e c t i v e s p e c i f i c volume f o r momentum, two-phase
f r i c t i o n g r a d i e n t m u l t i p l i e r v e l o c i t y and energy t r a n s p o r t v e l o c i t y .
Several c o r r e l a t i o n s a r e i n c l u d e d i n t h i s subroutine t h a t . t h e user can
s e 1 e c t . b ~option'. These a r e provided as an example w i t h t h e thought t h a t
t h e users w i l l s e t up c o r r e l a t i o n s t h a t a r e most a p p l i c a b l e t o t h e i r
p a r t i c u l a r problems. This subroutine i s n o t c a l l e d from t h e e x p l i c i t
s o l u t i o n scheme.
Subroutine XSCHEM c a r r i e s o u t t h e e x p l i c i t t r a n s i e n t s o l u t i o n . I t s
f u n c t i o n i s e q u i v a l e n t t o SCHEME f o r t h e i m p l i c i t s o l u t i o n . A normal c a l l
t o XSCHEM (IPART = 2 ) advances t h e s o l u t i o n one complete time step and
computes t h e d u r a t i o n o f t h e n e x t time step. For t h e running s t a r t options,
t h e steady-state i m p l i c i t s o l u t i o n i s f i r s t c a l c u l a t e d . The r e s u l t i n g
subchannel flows a r e then used t o determine t h e maximum time step allowed
f o r t h e e x p l i c i t t r a n s i e n t by a c a l l t o XSCHEM (IPART = 1 ) . Otherwise, an
i n p u t value f o r t h e maximum a l l o w a b l e time step i s used f o r t h e f i r s t time
step.
FUEL TEMPERATURES
+
I BEGIN CHANNEL SWEEP I
+
'COMPUTE TENTATIVE F L W S
EXPANSI W TERM
YES L
MAXIMUM COURAM NUMBER
CONSTITUTIVE RELATIONSHIPS
I n a d d i t i o n t o t h e s t a t e d c o r r e l a t i o n s , a minimum heat t r a n s f e r c o e f f i -
c i e n t o f 5 ~ t u / h r - f t z - " ~i s assumed.
CHF
I n ,q
A - B Forced convection
C - D T r a n s i t i o n b o i l i n g and t r a n s i t i o n
pool' b o i l i n g
!
- D - E F i l m b o i l i n g , low pressure f i l m
b o i l i n g and pool f i l m b o i l i n g
c' -
-.*, FIGURE B-1. Regional A p p l i c a t i o n o f Heat Transfer C o r r e l a t i o n s
POST CHF
transition boiling - McDonough, M i l i c h and in^(^)
stable f i l m b o i l i n g - ~roeneveld'~)
low pressure f i l m b o i ' l i n g - Dougall and oh sen ow(^)
pool f i l m b o i l i n g - ~erenson'~),
t r a n s i t i o n pool b o i l i n g ( h )
(a) M. Jacob, Heat Transfer, v o l . 1, New York, Wiley &. Sons, 1957.
(b) J. S. Thom, et.al., ' B o i l i n g i n Subcooled Water During Flow Up Heated Tubes
o r A n n u l i ' , Proc. I n s t n . Mech. Enqrs., v o l . 180, p a r t 36, 1966, pp. 226-246.
(c) V. E. Schrock and L. M. Grossman, Forced Convection B o i l i n g Studies, F i n a l
Report on Forced Convection V a p o r i z a t i o n P r o j e c t , TID-14632, 1959.
(d) J. 6. McDonough. W. M i l i c h , and E. C. King, P a r t i a l F i l m B o i l i n g w i t h Water
a t 2000 p s i q i n a Round V e r t i c a l Tube. MSA Research Corp., Technical Report
62, NP-6976, 1958.
(e) D. C. Groeneveld, An I n v e s t i g a t i o n o f Heat T r a n s f e r i n t h e L i q u i d D e f i c i e n t
Regime, AECL-3281 (Rev.) December 1968; r e v i s e d August 1969.
(f) R. q . n o n g a l l and W. M. Rohcenow, F i l m - B s i l i n q on t h e I n s i d e o f V c r t f c a l
Tubes w i t h Upward Flow o f t h e F l u i d a t Low Q u a l i t i e s , MIT-TR-9079-26, 1963. .
(9) P. J. Berenson, ' F i l m - B o i l i n g Heat T r a n s f e r from a H o r i z o n t a l Surface,'
J. o f Heat T r a n s f e r , v o l . 83, August 1961, pp. 553-558.
(h) K. V. Moore and W. H. R e t t i g . 'RELAP-4: A Computer Program f o r T r a n s i e n t
Thermal-Hydraulic A n a l y s i s ' A e r o j e t Nuclear Co., ANCR-1127, December 1973.
(i) E. D. Hughes, A C o r r e l a t i o n o f Rod Bundle C r i t i c a l Heat F l u x f o r Water i n
t h e Pressure Range 150-725 p s i a , IN-1412, J u l y 1970.
(j) P. G. B a r n e t t , A Correlarisrl o f Burnout Data f o r U n f f o n l y Heated A n n u l i and
I t s Use f o r P r e d i c t i n g Burnout i n U n i f o r m l y Heated Rod Bundles, AEEW-R 463, 1966
(k) $1. S. Gr?l.lerstedt e t . , a l . , ' C o r r e l a t i o n o f C r i t i c a l Heat F l u x i n a Bundle
Cooled by Pressurized Water; Two Phase Flow and Heat T r a n s f e r i n Rod Bundles,
Symposium proceedings o f Winter Annual Meeting o f t h e American S o c i e t y o f
Mechanical Engineers, Los Angeles, C a l i f o r n i a , November 1966, pp. 63-71.
COMPUTER PROGRAM CORRELATIONS
F r i c t i o n Factor
-
.rhi s c o r r e c t i o n i s based on t h e assumption t h a t t h e t o t a l perimeter* con-
s i s t s o f two regions--one heated and t h e o t h e r unheated--and t h a t t h e
heated p o r t i o n has uniform heat f l u x . t h e heat t r a n s f e r c o e f f i c i e n t i s
c a l c u l a t e d from
Two-Phase F r i c t i o n Mu1t i p 1 i e r
Armand ( 7 )
Polynomial Function
where t h e c o e f . f i c i e n t s a r e supplied as i n p u t .
Spacer Loss C o e f f i c i e n t
Void ~ r a c t i o n
S l ip Model
a = 0.
'
where y i s a s p e c i f i e d s l i p r a t i o .
a = 0. X -
< 0.
Subc0.01 ed Void F r a c t i o n
x = o xe < xd
X = X, - Xd exp (-
Xd
Xe
- 1) Xe/Xd < 1
Y~
AT = - ~ Q ( P , + l o g (1 + Pr(T - 1))) - 30
5 < YB.<
Phh
1
Several. forms of equations f o r specifying t h e t u r b u l e n t crossflow
a r e included. The p r e s e n t l y a v a i l a b l e forms i n COBRA f o r c a l c u l a t i n g
w ' incl.ude:
? --- .----..-- .
w l k = B sk 6 (8-20)
E a
-
where Re =
U
-
and a and b a r e i n p u t constants. Since a d e f i n i t i v e mixing c o r r e l a t i o n
does n o t e x i s t and o t h e r forms a r e a v a i l a b l e , (11 ,l2,13) t h e user should
s e t up . c o r r e l a t i o n s o f h i s choice.
Two-Phase Turbulent M i x i n q
Complete i n f o r m a t i o n concerning n i x i n g d u r i n g b o i l i n g i s n o t a v a i l -
able. It i,s known,,however, t h a t mixing i s s t r o n g l y dependent on q u a l i t y ;
t h e r e f o r e , COBRA i s s e t . u p t o accept B as a t a b u l a r f u n c t i o n o f q u a l i t y .
When t h e q u a l i t y o f two adjacent subchannels i s different, t h e c a l c u l a -
t i o n s use a q u a l i t y c a l c u l a t e d from t h e mean mixed enthalpy o f t h e two sub-
channel s.
Transient Correlations
C r i t i c a l Heat F l ux C o r r e l a t i o n s
and the entire integral taken as a summation over the increments of AX from
REFERENCES
where
Flux = I;[ ,,
f l u x o f K momentum through gap L can be computed by
W;
3
,o,e,
5
I n r e g u l a r c a r t e s i a n geome ry, t h i s o p e r a t i o n forms t h e exact centered
f i n i t e d i f f e r e n c e . a n a l o g t o apv lay. I n a r o d bundle, i t i s o n l y an approxi-
mation which i s most accurate when t h e gap space i s c l o s e t o t h e subchannel
h y d r a u l i c diameter. Test runs have shown t h a t t h e use o f t h e transverse
momentum term produces an e f f e c t in. t h e r i g h t d i r e c t i o n b u t t h e r e s u l t s have
n o t been compared t o experimental data.
C-2.
FIGURE. .-. Example of Gap D i r e c t i o n Angles
APPENDIX D.
1. t o a u t o m a t i c a l l y s e t up a c o n s i s t e n t s e t o f dimensions f o r COBRA
compatible w i t h t h e u s e r - s p e c i f i e d problem size;
2. t o c a l c u l a t e t h e r e l a t i v e storage l o c a t i o n o f t h e v a r i a b l e s
equivalenced t o t h e v e c t o r SAVEAL;
ME = MX i f no storage on p e r i p h e r a l s
ME - 3 f o r storage on p e r i p h e r a l s ( t h i s i m p l i e s t h a t a x i a l values w i l l
be r o l l e d i n COBRA)
*AF N E k P L t Y A N K $ % $
*CO SPtCS
YC"
HGr
CI X a
Ft Nt
HH 0
WE+
PPn
"LO
M 2.
MY.
Hva
nI a
no=
NIr
N(E
Cl Kn
R4E
nta
PSP
?I 10
*CD ,SPEC 1 -
DIPENSIC& A C ( ~ C ) ~ P ~ ( M C ) ~ P ~ ( ~ C ) ~ C C ( ~ C ] ~ O R ( ~ R ) ~
1 NI~LET(PC)~FINLET(~C)I~~~L€T(~CI~
2 LC(MC,U),GAPS("C,Y),DIST(MC,U), ,
3 D A l Ih(hZ1
EQ~I~ALENcE ( D A T I ~ ( ~ ) ~ ~ R ( ~ ~ ~ E ) ~ ~ ~ I ) ~ T I C ( ~ ~ ~ ~ C ~ ~ ~ ~ I
1 (P*(l)rSTU~(lrhEl),(FH(l),rP(l,rE)),
2 ~ ~ ~ ~ l l ~ ~ P ~ ~ ~ ~ l ~ ~ t ~ ~ ~ ~ l ) ~ ~ ~
3 ( F I ~ ~ E ~ ( ~ ) I V ~ S C ( ~ ) ) , ( T I N L E T ~ ~ ) ~ ~ I ~ C ~ ~ ~ ~ ~ ~
U ( L C ( l ~ ~ l ~ ~ F I L ~ ( l ) l r ( G ~ ~ ~ ~ l ~ f ~ t P E ~ I ~ ( l ~ ~ ~
5 (01ST(l,1)1~AA(1~~3))
CLCM ~ ~ V EL.~,~~~I~,AC,P~,~H,DC,O~,MI?~L~T,FI~LE~~TI~LET,LC~
CLCH 1 GAPS,UIST
*CD SPEC2
COHHCN ~U(~!G~~J~(PG),X(~%)~A~GL€~?G)PFAC~O~(~G)~LENG~~(~~G)~
1 A X I A L ( ~ P ] , Y ( ~ P ) ~ U ~ ( ~ ~ C ) ~ A L P ~ A ( P C ) , C C ' A L ( ~ C I ~ ~ ~ T C ~ T ~ ~ ( ~ C ) ~
2 ' S~HPF(hC)rdIl>T~(t~'k),IK~(fih),Jk~(plN],UkCP(~,~W),
3 t f i ~ ~ ~ i i I, . C
i (d ~ ~A I . L ( Z , ~ I ~ ) , C . . ~ ~ @( I. -~ C
u *PI )
I K I J , A F L L X ~ Z , T ~ E ~ A ~ P I ~ ~ A X , I ~ ~ I ~ I ~ C ' R ~ G ~ ~ , J ~ , J ~ ~ J ~ ~ J ~ ,
1 J~~ATC~~L,G~X,EL~V~FE~R~R,ITE~A~~YLOX~DTGCIC~GC~I~TGC~GC
2 3L~CX~SLUXI~SLGOX,SLG~X~,O%UI,SLCXOlrO~2~DkCT,PI~~PCUt~F~T,
3 F ~ ~ L ~ ~ K ~ ~ ~ ~ ~ ~ A P , ~ Z ~ ~ ~ ~ Y ~ E Q ~ ~ I U U ~ L U I ~
0 S ~ D T ~ S ~ D T Z ~ S ~ G C I A P ~ I I ~ T I D X I G T ~ ~ ~ R ~ ~ E R H X ~ ~ € ~ W Y ~ I ~ U Y ~ I T
5 A C C C L Y , A C C ~ ~ ~ , A C C E L ~ , I ~ T ~ Q E , N A A ~ , ~ A ~ A P ~ , A ~ ~ E L F ~ ~ A M P ~ G ,
b ~ A A ~ ~ ~ A A Y P ~ , H ~ Y I Y ~ U S O O ~ I I T L , P ~ O ~ ~ ~ F C ~ K B F C ~ U ~ ~
*,EhOl
COPKGN /LCCP/ ~ C M ~ ~ L , N K I N D X P ~ ~ K T R Y X ~ N T ~ ~ Y Y ~ ~ ~ I ~ . K , ~ R U D I ~ ~ O ~
,E h 0 2
COkrCtJ / L b R G € l / S A V E A L ( ~ ~ ~ , Q P R I ~ ( F C ~ ~ V ( P C ~ ~ V P ~ ~ ( ~ C ) , V I S C ~ H C ~ ~ ~
1 V I ~ C ~ ~ ~ C ~ ~ H F T L ~ ( H C ) ~ ~ ~ ~ ( ~ C ) ~ C P ( ~ C ) , F S P ( ~ C ) ~
a PEPIH(PC),PPERIM(P:C~,~P~PC),AA~LMU,~O)
*rkbO3
CO!'HCN /LAI?GEJ/ A N ( ~ C ) , O H Y D ~ ( * C ) ~ L O ~ P ( ~ ' ~ ~ ~ G ) ~ N ~ H A P ( ~ C J ~
1 N~QAPS(HC~~VP(*C),T(~~)PILOCS(~III~C)
*,Eh04
INTEGER L b E V A S ( k 2 1
E ' J L I V P L E ~ c E( S A V E A L I L A E V A S )
REAL I D T G C , L E N G f d r K I J
. .
I .
1 ' (RHO(l,Il,SAVEbL(~C))t(~~LI:(1,1),9AVEAL(~~)),
2 ( A ( I ~ ~ ) # S A V E P L ( Y S ) ) ~ ( ~ ~ Y U ( ~ , ~ ) I ~ A \ ~ ~ ~ L ( ~ L ] ) *
,
5 ( ~ ( ~ ~ ~ I , ~ A ~ E ~ L ( N ~ ) I , ( F C L U ( ~ ~ ~ I I S ~ V E ~ L ~ ~ X ~ I ~
(1 ( ~ ~ ~ U ~ L D ( ~ ~ ~ ) ~ S ~ V E ~ L ( ~ ~ ~ ~ ~ ( R ~ ( C ~ ) A ~ ( ~ . I ~ ~
S ~ L ~ ( ~ , I ) , S ~ ~ E P L ( ~ ~ ) ) I ( ~ T . U * ~ ~ , ~ ~ ~ S A V E A L ( N ~ ) ~ ~
b (AP(l,l),SAbEAL(NZ))
CLCM LEVEL 21' P
CLCM LEVEL 2r P
CLC* LEVEL 2, F
CLCH LEVEL 28. QWL'
CLCV LEVEL 2r bCLO
CLCM LEVEL 2r FCLC
CLCH LFVkl. Zr QKOCLO
CCCV Lf'JkL 2, A
,CLCM LEvtL 2r rbYC
CLCY LEIEL 2, LB
CLCH LEVEL 2r STUa
CLCM. LEVEL 2, b~
CLC* LEVEL 2,
ECcCflAE
*CD S P E C u .
' DIPEkSICb ~ ( M ~ , ~ E ) , * C L O ( ~ ~ , I Y E ) , ~ ) E T A ( V ~ ~ ~ E I ,
1 U S ~ V E ( ~ ~ t ~ E l ~ F S . I ( ~ 1 ~ ~ E l t T k F ~ ~ l ~ ~ E ) ~
2 G P P ( . r l , n E l # T*ALL(.pl,C'E) t T ~ C L O ( *, ?~ E l
EQLIVALEkCE ( * ( ~ ~ I I ~ S A V E A ~ ( ~ ~ ) , ( ~ ~ ~ L U ~ ~ , ~ ~ , ~ A V € ~ L ~ ~ ~ ) ~
1 ~~ET4(l~ll~S~V~EAL(~U)~,i(~S~VE(ltl)~S6.btAL(HS)l~
2 ( ~ S I ( ~ ~ ~ ~ ~ S ~ V E A L ( ~ ~ ) ) , ( G ~ ~ ~ ~ I ~ ~ ~ S
3 ( T ~ ~ F ( ~ ~ ~ ) ~ S A ~ E A L ( ~ B ) ) , ( T ~ A L ~ ( I , ! ) ~ ~ A ~ E A ~ ( ~ Y ) ) #
1) (T~CLG(lrl),SAVEAL(No)l
CLCV L E V E L 2r
CLCt' LFIEL 2, kCL0
c L c n LEVF.L 2, RETA
C L C ~L E V E L a, PSAVE,PSI
CLCH L E V E L 2 , GAP, T k F
CLCY LEVEL Z,T*ALLPT~OLC
rCD S P E C S
DIFEhSICN cPFH(Yl,~'E),CCh~'dL(~-l,rE),FLl~X(~lrPE)t
1 ~ ~ U ~ ( ~ N , N U ~ M E ) ~ H S U ~ F ( ~ ~ ~ ~ ~ ~ ~ ) , I ~ T T P ( * ~ I C E ) P
2 ~cG€(r\,rE)
EGLIvALENCi: ( C ~ F ~ ( ~ ~ ~ ) ~ S ~ V ~ ~ A L ( ~ ~ ) ) ~ ( C . C ~ A ~ L ( ~ , I ) ~ S A ~ E A L
1 ( F L U ~ ( ~ ~ ~ ~ ~ S ~ V E A L ( ~ D ~ ~ ~ ( ~ ! ~ ~ ~ F ~ ~ I ~ ~ I S ~ ~ E A
' 2 ( ~ ~ U E ( ~ , ~ ~ ~ S A ~ € A L ( ~ Y ~ ~ ~ ( ~ K C G ~ ~ , ~ ~ ~ ~ , ~ ~ V E A L
3 (IF~YP(l~l)~ChFH(lrll)
1PTEGE.R CCkANL,
c L c n L ~ ~ 2E, LC ~ F M
C i c M L E b k L 2, CCw.4NL
CLCH L E V E L 2 , F L U X
CLC" L F V E L 2, T G C C
CLCH L E V E L 2, HSLidF
CLCM L E V E L 2, I F T Y P
C L C n L E V E L i?r v C C E
*CO S P E r t
LOGICAL GGID
ccmrub / P G Y X O / I G R I D ( ~ Z ) ~ ~ ~ G R ~ I O , M G R I D T , J ~ , ~ G T Y F E , N ~ ~ A ~ P ~ T ~ I I C K ~ D I C ,
1 GUIOxL(~Z)rCD(~C,hK),FCX\JtEcG),FXFL.Ch~MG~MK),,
2 O U R ( * G ) ~ P I T C H , ~ B A R ( ~ ~ ) I X C H C S ~ ( ~ G ~ ~ ) ~ G H ~ O
*, E h 0 5
TABLE D-1 . . Continued
OI~E~SIGH CNRFLX(HX)
E~LIVALE~CE (~NbFLX(l)rnfjAR(l))
c c * r r - w / I c ~ h T / ICUU!VT
. . . . . .. . .. .. - . .
CLCM L k ' E L 2, ChbFLX . . .
CLCM LEVEL 2 1 I G H I O
4co s~tt7
.
COPeCN / g ~ f i E A / A F ~ c T ( ~ L , P A ) , G F ~ C T ( ~ L I ~ ~ ) , ~ X L ( ~ ' L ) , C ~ P X L ( K L ) #
1 ~AWARU,%~FACT,NAX~,NGXL,~GAPS~YC~(MA),NGAP(~S)
*,EhDb
CGr',!CN / L A E G E ~ / I D A ~ E A ( M C ) ~ I U G A ~ ( P G ) ~ G ~ P N ( M )
*,Eh07
CLCH L E V E L 2, ICARE~
*CD SPEC(!
D1"ENSICbi D P D X ( P C ~ ~ C Y D X ( ~ C ) ~ C F I ~ X ( ~ C ) ~ W ~ T ~ ( P ~ ~ E )
ECUIUALENCE ( C ~ P ~ ( ~ ) ~ C P ( ~ ~ ) , ( C F D % ( ~ ) I F ~ P ~ ~ ) ) ~ ~ O P O X ~ ~ ~ ~
1 (P(l~1),~9T~(l,l))
CLCH -LEVEL 2 , C P O X P ~ ~ ~ X , O F O X I R S T ~
* G O SPEC9
corr!CN / P Q C P l / P P ( ~ P I , ~ ~ ( ~ ~ ~ F ) I v V O ( ~ P ~ ~ VH V~F (-C PF ) ,~ ~ P .~ ~ ~ C C O ~ F ~
1 KKF(~P),SSIG~A(KP),AA(~)~~M(~),CC(U),UUF(MP),
2 ~GPT(I~P),T~'CUT(PP),~PPUCT(~P)~KFI:~T~~P)IUP~(,T(FP),
s ~ P O ~ T ~ ~ P ) , O T ~ A ~ , P ~ ~ F C L , G C W I T ~ A A L ~ ~ ~ ~ ~
U ~IPHOP,NVXSC~,NT~PE(HC),ISTEA)L,L~~NF
h,EhOd
RE.AL KPCUT,KKF
*CD S P E C 1 0
CCPYGN /PGCP2/ P W E P , ~ F ~ V F , V G ~ ~ ~ ~ H C ~ I ; P ~ K F , S I G ~ A I ~ ' F G I V F G ~ R H G F ~
1 RHOG,AF(~),~V(~),NV,~F,J~UIL(~C),AH(~)
*, EhDQ
RkAL K F
b C D SPEC1 1
COrnCN / n I X f / A ~ E T A ~ B B ~ ~ A , X ~ U A L ( M P ) ~ ~ X ( ~ P ) , C ~ , ~ T M , ~ H ~ ~ ~ ~ J ~ , ~ S
a,EhDlo
*CO SPECIZ
cOFt1CN / F b E L / KFU~L(~TJ,~CL~U(~T),~FUEL(~T),RCLA~(~T)~CF~EL(~T),
1 CCLA~(~T)~TCL~~(~T)~TFLLIU~~GAP(~T),LY(KU,~)~
2 . P H I ( ~ U ~ ~ ) , ~ ~ O ~ A L ( H ~ ~ , O ( C R ) , P @ ~ € ~ ~ G F ~ E L ~ ~ ~
J FLXX(~R),~P~F,ZEND(MT,~Y)~TCLOR(C~,~~)~IDFUEL(HR),
Q N ~ J A ~ P rNn T
C r ~P c n p ) , t i < r i N E ( w T l ~ I Z T Y P ( ~ T , H Y ) , N R C D T P
+rEIDll
REDL %FUEL,KCLAD
CLCH L E V E L 20 %FUEL
*CD SPEC13
CCIPNCN /BCHF/ ~ C H F R ( * X ) , M C ~ F R C ( P X ) , ~ C H F ~ R ( H X ) , ~ C H F , ~ C H F X
*,EhOli!
- REAL. PCIiFR .
*CD B ~ E c I ~
bInEhsICN ~ S ( H ~ , M E ) , ~ S * ( M ~ ~ M ~ ~ I ~ E ~ ~ ~ ( ~ C ~ ~ T F : P ~ ~ ( ~ C J ,
1 CPK(PC),A~P(PI,M~)~~LC~(CC~~~(P~)~~IJ(M~)D
2 ~ A ~ ~ ( ~ ~ ~ P E ) t ~ P ( ~ l ~ U C ) ~ O P ~ ~ ( ~ C ) ~ I F U I ~ ( U l i l ~
5 ISTAEi(!4G),FJCC(!$C)pbAb(hf), IcCC(~B)~~PP(MG)
EQLIV.ALEF;cE ~~J~(l~l)~US~~t~l),~~I(l~l))~(Tt~~C1(~),VISC~(~l),
1 ~ T ~ ~ ~ 2 ( l ~ ~ S T ~ ~ ( l ~ l ) ~ ~ ~ O ~ ~ ~ l ~
Z (~PP(t,lltA~~l,l)),(PLCrJ~ll,8l~~+(lt~)~1
3 (~(l)~~P(l)~C1J(l)r~P~~~l))~~MbRU(l~l),BET~(l~l))~
i) ~ S P ~ l ~ l ~ t ~ U ~ ~ l ~ l ~ ~ ~ ~ ~ ~ ~ ~
5 (IL~C(l)~ILCCS(1~1))~~OP~P(l)~FSF(I),HJ~O(~)l~
b (ISTA~(l~~IPCIV(l),F~OIV(l))
CLCM LEVEL 2 , U ~ , ~ ~ ~ ~ T E ~ ~ ~ , T E ~ ~ ~ ~ , ~ P K , A P P , F L C ~ ~ ~ ~ C ~ J , E ~ R K , ~ P P
CLCM LEVEL 2 , SF,AAC,ILCL,OPWP,~J~U
CLcn L E V E L 7, ISTA",IFDIV
TABLE D-1 . Continued
t C 0 SPEC15 .
91PEiiSIc~ SAVEA~(H~),SALEAZ(M~),SAVEA~(K~),
sAVdkS(!.iS)
EnLIbALEKCE (SAVF'A~(~),SALEJL(~)I,(SAVEA~(~),~SAVE.AL(NK)),
1 ( ~ A v E A ~ ~ ~ ) , ~ A ~ E A ' L ( ~ ~ ) ) , ( s P v H E s ( ~ ) , s P v E A ~ ( ~ ) )
cLcr L E ~ & L 2, S ~ V E A ~ , S A V E A Z , S P ~ E A S
CLCh LEVEL 2, 3AL2kS .
' *CO SPECI!,
!-jlrt'hSrCt.: S I ! ~ ' ~ ( I ~ C ) , S I J ~ L ( ~ C ) , S ~ I * ~, S- (Trrcj () n ~ ) , d T b J r i ( n ~ )
E'.:LIVALEbCE ( S ~ ~ ~ L ( ~ ) I V I S C ( I ) I ~ ( S C , ~ ~ H ( ~ ~ ~ C P ~ ~ ) ) ~
1 (S~~X(l)~vISC~(l)),(STHJ[l)rFSP(l))r
i! ( 3 l ~ J ~ l ( l ) ~ C t i h ( l ) )
CLCH LE'+EL 2 , S L P . L , S L @ ' P , S O " % , S T h J t S T ~ J H l
*CO Y P E C 1 7
D T " E h S 1 C V AAAPk:l(t;7) t L F k ~ ( i " 2 ), L @ C A l ( h b )
E ~ ~ ~ I V A L E ~( CA A
E A(l,l)nAAAP~~(t)),(fi(l,l),~Pfi).r(l))t
1 ( L C C ~ ( ~ , ~ ) , L C C A ~ ( ~ I )
cLcn L E ~ E L~ , A A A P R M , ~ P W ~ , L O C A ~
OIrth3ICN CS~~~~(~G),UZPUM(PG)ILIEIFWI*(~GI~SPP~~C~ML~~~
1 h C i L n F l q ( * G l ,Ufjf'R+{(WC), APPFfit'(pC) ,kSAVk ( k l i ) ,
2 CTII("G)rIITJJ("G)
ErJLIVALENCE. (~~SAVE(l),dSAvE(1,3)),(DTII(l)tEEA~E(l~l~~~
1 (OlJJ(ll~tiS~Vt(1~2))~~UY~FR~~l),LS1*~l~1ll~
2 ~ ~ ~ S P K t ~ ~ l ~ ~ ~ S ~ l ~ l ~ ~ , ~ ~ E T P k ~ ~ l ~ ~ B E l
s ( s ~ P ~ ~ ( ~ ) , s P ( ~ ~ ~ ) I , ( ~ c L u P. c ( ~ ~ ~ ) ) ~
1 ( u r ? ~ r ( l ) r u ~ ( l ~ l ) ~ ) , ( ~ P P ~ ~ ~ ( l ) , ~ ~ ~ ( l ; l ) )
CLCM CEbEL Z , , u s * P W M , ( ! S F 4 ~ , H E T P R r , s P C ) R M
CLCN L E V E L 2, ~ C L C P ? , U ~ P R Y D A P P P ~ ~ # ~ ~ Y A V ~
CLCM L E b E L 2 , o T I I , U T J J
. *CO S P E C 1 8
O I Y E h S I ~ N~ l ( f l 2 ) , ~ A ~ l ( ~ 7 ' ) ~ e s 4 v E l ( ~ ~ ) , ~ ~ I 1 ~ ~ ~ l , L ~ C ~ 1 ( ~ 6 ) r
I P("G)rSETAl(Ht)
EnbIvALENCE (Al(l),A(trl)),(~AAl(l),AA~(l~l))r
1 ( H S ~ V E 1 ~ 1 ) , ~ S ~ ~ F . ( l ~ l ~ ) ~ ~ ~ S I l ( l ~ ~ P S I ~ 1 ~ 1 ) ) ~
Z ( Y ( ~ ) ~ ~ P ( ~ ) ) ~ ( L C C ~ ~ ( ~ ~ ~ L C C A ~ ~ ~ ~ ~ ~ I ~ B E . T
CLCH L E V E L ~ , A ~ ~ A A A ~ ~ ~ S A V E I , P S I ~ ~ L O C J ~ ~ B , H F T A ~
*CO S D E C 1 9
c C ~ ' I * C ) ~ / I N P A I ~ /C I ~ ~ ~ P U T ( ~ ~ ) ~ ? R I ~ T . ( ! ~ ! ~ , T E X T ( ~ ~ ) ~ C A, .~ E ( ~ ) P ~ I ~ E . ( ~ ) ,
1 T D ~ ~ Y ( ~ O ) ~ P ~ T U T , ~ T T I ~ & I ~ T I " ~ , ~ I F ( ~ ~ , Z P T ( ~ ) ~
2 YP(~P~,FP(rP),YC(~P)nFG(~P),YP(YP~,Fk(MPlt
3 YO(P~I,F~(~O),YT(P?),F~(~~),F~)I(*!P~,Y~'X(PP),
11 k~bT,kFX~hCUT,hPCFAFi,tjPC4.Pt~Fh;CCE,hPRU@t~PLTCht
5 .. P ~ I N ~ ~ ~ ( ~ ~ ) , ~ ~ L T G P , ~ ~ S Y ~ ' L T , ~ S ) ~ I P T , N S K I P X , M N , ~ R ,
6 ~ C , ~ G , ~ P , ~ X I ~ ~ * ~ H Y ~ ~ E I I * Z , ~ ~ , ~ ~ , ~ P , P ~ I , ~ ! T , H L ~ ~ P ,
7 ~!G,*;orhH,IC6T,IH(VS)~J~~y.S)!~AxT,~JuYP~
8 I ~ ~ P ( ~ O ) , J ~ P ~ ~ U ~ , I P , I ~ ~ I ~ I I F ~ ~ H , I F ~ ] L T ~ ~ K ~ O ~
0 hDTr:JuTrl~rASE~~T~IES~ITSTkPtT~X~TPD~~~OFLC~J~CFLO
+,HI~,TIN,CIN,LZ,OXX
* , E ~ ~ ~ ~ , ~ I G ~ A L ( ~ ~ ) , F ~ , H Z , + ~ , H ~ , H ~ , ) . ~ I H ~ , ~ ~
C C P ~ C/ L~A E G E Z / PHINTC(~C),PRINTG(~G),PRINTR(PF),~P~TCH(MCI~
1 IPLTCP(~C)~~EXIT(nc)~b9AR~rc~
',EhUlU
If+IEGER P Q I ~ T C , ~ R I C T G , P R I ~ T . \ r , ~ P H I % T A
LCGICAL PSINT
CbCM L E V t L 2,FRIhJTC *
'CD SPEC20
DIrEhSICh S A V ~ A ~ ( ~ ~ ) ~ S P V E ~ ~ ( M ~ ) , S A V E ~ ~ ( Y I )
CLCH LEVEL Z I S A V E A ~ , S A V E A ~ , S A V E A ~
-..
TABLE D-1 .. Continued
*CD SPEC21
O I ~ E ~ S I CP ~n l t r c )
ECLI~ALE~Ct(PhI(l)rFSP(I))
C L C M LEVEL 2,PhI
*CO SPEC22
O I ~ E ~ S I CF ~( ~ P ) I Y PI
*COHDECK S P E C 2 3
CLCM L E V E L 2 ~ 3 Q C 2 3 r l ~ ~ C 2 3 1 S ~ C 2 ~ ~ I S ~ c 2 4 ~ 8 P C b l ~ I S P C b l
C L C H LEVEL ~ , S P C ~ ~ ~ I S P C ~ ~ ~ S F C ~ ~ I I ~ P C ~ ~ I S F C ~ O I I S P C ~ O
*CD 3PkC2b
CLC3 LEVEL 2 , X , l r Z
*CD SPEC25
CLCq L k v E L Z I ' J A K ~
*cn S P E ~ Z ~
CLC? LEVEL ~ ~ I s ~ C R E , S T ~ J U E I S T C R E ~ I L O G I C L
*CD SPEC17
CLCH L E V k L i ? r A , 0
LAST
APPENDIX E
PROGRAM GEOM
19 ROD
37 ROD
Card 5 (4(2F5.3))
1. i n t e r i o r - i n t e r i o r connections
2. i n t e r i o r - s i d e connections
3: side - s i d e connections
4. corner s i d e connections
Card 6 (611)
V a r i a b l e (Dimension) - -
Subroutine D e i i n it i o n
52 . input Subccoled v o i d o p t i o n . (Group 2, N1)
= 0, no S U ~ C O O ~ E : \lords
= 1, Levy s u b c ~ o l e dv o i d model
SPEC 2 input Option f o r bulk void correlation. (Group 2, N2)
= 0, homcgenouj b u l k v o i d model
= 1, m o d i f i e d Arzand, b u l k v o i d model
= 5 , .homogenous mcdel w i t h s p e c i f i e d s l i p
ratio
= 6, b u l k v o i d f r a c t i o n i n p u t as an Nth
o r d e r polynomizl f u n c t i o n of qua1 it y
SPEC 2 input I n p u t data o p t i o n f o r two-phase f r i c t i o n -
mu1 t i p l i e r , (groiip 2, Ii3)
= 0, homc)eneoirs i ~ 0 d 2 1
= 1, Ani1andtf1-iction m u l t i p l i e r
= 5 , f r i c t i o n ~ u l t l p i l e ri r ~ p u las ail Nth
o r d e r polynomial f u n c t j o n o f q u a l i t y
input R a d i a l t h e r z a l c o t i d u ~ ~ l oupL
n lull. (01i r ~ p10. 113)
= 0, conduction absent
= 1, c o n d ~ c t i o np r e s e n t
SPEC 6 input Flags type u f spacer design. (Group 7, N1)
= 1, w i r e wraps
= 2, g r i d spacers
= 3, both
SPEC 10 PROP I n d i c a t e s a x i a l l o c a t i o n o f onset o f sub-.
c o o l e d b o i l i n g i n Subchannel I.
SPEC 2 SETUP see I %
SPEC. 2 SETUP see IK'II
SPEC 19 SETUP see 11.: .
SPEC 19 RESULT see If.:?
SPEC 19 input Determines t h e t y p e o f AP t r a n s i e n t . (Group 11, N7)
= 0, prcssi!rs r l 1 - n ~i s (liven i n p s i
= 1, pressure drop t r a n s i e n t i s a f r a c t i o n o f
t h e s t e a d y - s t a t e pressure d r o p
Numeri z s l s o l u t i o n o p t i o n . (f;rolrp Y, ~ 3 )
= 0, i m p l i c i t s o l u t i o n
= 1 , expl i c i t s o l u t i o n ; pressure boundary ,
condition, running s t a r t
= 2, e x p l i c i t s o l u t i o n , pressure boundary
condition, standing s t a r t
= 3, e x p l i c i t s o l u t i c n , i n l e t f l o w speci-
f i e d running s t z r t
= 4, e x p l i c i t s o ' l u t i o n , i r l l e r f l a w s p r c l -
f i e d standing s t a r t
'
K12 SPEC 2 SETUP S e t t o one i f gap o r i e n t a t i o n angles a r e used.
KASC SPEC 19 input Case number, (case c o n t r o l c a r d )
KCLAD (MT) SPEC 12 input Thermai c o n d g c t i v i t y o f c l a d d i n g m a t e r i a l ,
Btu/sec-f t - O F . (Group 8)
SPEC 10 PROP- S a t u r a t e d l i q u i d thermal c o n d u c t i v i t y a t t h e
r e f e r e n c e pressure, Btu/sec-ft-OF.
COBRA- I V - I VARIAaLES DESCRIPT 1011 ( S o n t i nued)
Variable (Dimension.1 - -k
Colndec Subroutine Definition
NSCBC SPEC 11 input ~ u ~ c o o l emixing
d option, (Group 10, Nl )
= 0, 6 = W/SG = a
= 1, U/GS - a(Re)b
= 2, W/GD - a ( ~ e ) b
=. 3, W/GS - (S/ZIJ) a ( ~ e ) b
NSKIPT SPEC 19 input Output option. P r i n t every NSKIRT time steps.
(Group 9, N1)
NSKIPX SPEC 19 input Output option. P r i n t every NSKIPX a x i a l nodes.
(Group 9, R1)
NSKPLT SPEC 19 input Output option. P l o t every NSKPLT time steps.
(Group 12, N6)
NTRIES SPEC 19 input External i t e r a t i o n 1i m i t f o r imp1 i c i t s o l u t i o n .
(Group 9)
NTRYX SPEC 2 input It e r d t i o n l i m i t f o r e x p l i c i t flow f i e l d s o l u t i o n .
(Group 9)
NI'KYY SPEC 2 inpu't I t e r a t i o n 1i m i t f o r e x p l i c i t energy s o l u t i o n .
(Group 3 )
NTYPC(MC) SPEC 9 Input Flag f o r channel type f o r f r i c t i o n f a c t o r
NV SPEC 10 input Number o f c o e f f i c i e n t s of polynomial function
expressing two-phase f r i c t i o n m u l t i p l i e r ; see AV.
NVISCW SPEC 9 input Option f o r w a l l v i s c o s i t y c o r r e c t i o n t o f r i c t i o n
f a c t o r . (Group 2, N4)
= 0, no c o r r e c t i o n
= 1, c o r r e c t i o n used
NWK SPEC 2 input Number o f w a l l connections. (Group 4, R4)
NWR SPEC 2 COBRA Maximum number of.connections t o a channel
(thermal + - f l o w ) .
SPEC 2 input I n i t i a l wrap inventory f o r each subchannel.
NZONE(MT) SPEC 12 input Number of fuel zones.
OUTPUT ( 12 ) SPEC 19 RESULT Array use? t o p r i n t . output.
P(M1 ,ME) SPEC 3 SCHEME R e l a t i v e pressure a t a x i a l l o c a t i o n J i,n
subchannel I,l b / f t 2 .
PERIM(MC ) SPEC 2 COBRA Wetted perimeter, PW/12, ft.
PEXIT SPEC 2 input System operating pressure, psia. (Gro.up 11 )
PH(MC) SPEC 1 input Subchannel heated perimeter, i n . (Group 4)
PHI (MC) SPEC 21 VOID Two-phase f r i c t i o n pressure g r a d i e n t m u l t i p l i e r .
P H I (MK,ti) SPEC 12 . input ~ r a c t i o no f the heated perimeter. o f rod N f a c i n g
subchannel I. (Group 8)
PHTOT SPEC 19 SETUP T o t a l heated perimeter, ft.
PI SPEC 2 COBRA Constant I?, PI .= 355./113
PIN SPEC 2 COBRA Pressure drop ( p s i ) . f o r pressure boundary condition.
PITCH SPEC' 6 input P i t c h l e n g t h o f w i r e wrap, ft. '(Group 7 )
POUT SPEC 2 XSCEM E x i t pressure, l b / f t
2
.
POWER SPEC 2 COBRA Value o f FQ a t appropriate t r a n s i e n t time.
PP(MP) SPEC 9 input ~ r e i s u r ea r r a y i n the f l u i d droperty t a h l e , psia.
(Group 1 )
PREF SPEC 10 input System reference pressure, psia,. (Group 11 )
COBRA-. I V - 1' VARIABLES DESCRIPTION (Continued)
. . SAMPLE PROBLEMS
. O A T & F O H t r Y L r ~ I lS O L U T ~ O N
t X T L i i h d L I T E Y A ~ I G cL I n 1 1 24 HIfvInuM A v I k h b A L A T t W P 1 I C h S 5
I t t l k R N q L k T E w ~ 1 1 C uL I M I T 36 F R A C T l U h UUhOR CLLL U B l A k $0,0000
C i 1 b V E Q G E h C C FbCTCH3 ACC€bEnbT l U h f AClLA6
t x v t ~ ~ s( ao ~~ j h i etbuo CaC38FLCc SGLCVlOh 4,6000
11.fERb~L ( o e t k ) ~ O - O ~ O L A T E R A L C ~ ILA.P ,no00
f LO# (OFIF1 OOOUO $LO* I ,ocloo
. .
R I X I ~ ~ GC C ~ U E L A T ~ C ~ S
JbfiCclCLCn m I u I h G 1 BET4 6 ,.0200
G) H C l L I h G f'IXIh.G, bET? 1 8 B S ~ U W E D ' ~ ~A8
M LS U B C ~ G L E C
I
m
O P E A A T I ~ ' C~ O t ~ g I f I t i . h . 5
SYSTE* ?WESSUPE a 6 0 . 0 PSXA
IhLkY EhTnrLVy C 9968,,rd- B l U / L R
b v G , r A 3 3 VLLCCITY C 3,OCO M I L L l O P L M / [ H R i S C B ? )
, I s L t T T E ~ P E R A T U ~ C8 800.,O O t G H k . E S F
A v G , q t A l FLUX e ,5D001:'0 * ? L L I n N t!TU/[HR=SGFl)
U h I F f i a ~ I N L E T TEhPEEAlIJRi .
IJKJFL'kM XIVLET M A S 3 WELOCIlV
7 - P I N WIRE-WRAPPED BUNDLE SAMPLE PROBLEM RESULTS (Cont inued)
,lU905E408 bH/6EC
P A S 2 FLD* CLT .lu905L t o 8 LB/SkC
rbas FLcw EwwuW m p 2 P 7 J d E e 1 0 LR/GCC
tCUBL
olrab 1 T V
O,(rO6
GaUOO
OrOOO
(0,000
OeOOO
0 e 0 0 i)
oso00
C a OOiJ
0 ~ 0 0 0
OaOOO
U,OOO
6 ~ U O O
051.51' a L r ~ ' 9 i f4?0hls I-'~I)o~o @*Of
9SLSt8 ' 1601"!1 bb"OIS 6ZOC0 OeSF
QSIS?e 1119'90 ILOLOS iinl)v" Oenl
QSLS?' SSh7'90 OhmhOS SLlbO aeE6
'J:LS2° nb19@ut a l 0 l o ~ zbOz0 I 0 O ?l
9SLSP' k$Ris91 Ob"6bh 01S5'1 o 0tr
951St' RSSS'Q~ E9'0bfI oZ.Eho % 0 0.5
QSLSt?' SS?5'90 SF.'? b n .@Enlmt OOb?
o;15t0 Lhhhe9! 0O0wC)h 25nnot OORZ!
QGLS?' 1 7 v n m ~ 1 Obshdh hYSLm2 l"1.7
9SLS?' !rCbm9t CSe90~ 2 8 5 ~h U0S?
,PSlc;?' rF?k'Qt 9Z0w~n OeQhes 0 O 57
o;l.sts 9S~C'91 ae'b~ti t l ~ g ~ h OOR?
Q;LC,t?* ZRkC'91 IL'IL~ lZbooh 1)'fi'
QGLG?' OSl~'Q1 )rv"t~~v 3562'~ n
Oil.S?' saPtfQ1 9 1 '5011 ~ 9 b ~ 0 n U"I?
JSlcl?' , 5bSF'QO b0OfPV alob"n 0°1;7
95LST' O1fEm'JI tPOtiSh' bnli.'es tObt
Qs/.5?s bf.0t090 S t OSSn 1115°S b06t
0 5 1 SF' ~SL0'90 ~Om?s0 bOf d e s UOLL
75LSZ' bSh!OP! ~ ~ ' ~ I I ~ lVn F I D o C03t cn
3il.'i?' ~SO1'91 f ~ ' 5 b n ~ L ~ O O Y OOSl 0
0
05/.SP0 kLROOQ! veeirnr, D>SL*J (! n 1
OiLqt?' 9653O91 ab'bS4 9 C ~ i ~ ~ l cer\
9 ; l G?' 9?fdm91 lL'Stf9 Z*Cf ' 1 ? ~ ? l
Q;1.Grm hS00'91 nnm?<v .at&3'1 G'lt
Q;lc;Fr hHLbmS! olebZn LC*O.L . L'*r:1
94152' ?1Sbm6l hho<-?n EQbZee (! * b
9;Lsi" ?h?b'St ?9'?Zn Ibbi'4 u'd
~5/.5?' 0Lb~'S1 nZ"bln 9Slb0c t oL
9il'it' tola-st L O ' O I ~ b~17'b 3'u
OsLq?' tFh?'SI !tb*t?IV (;1\5@b T'S
191J'SO ZS'bdn b4b d e b (8 t
*;l=#?@
Q~Ls?' I LRL'SO S?'ctOb 8:-51'71 U'F.
9SLSZ' Z?01'51 ~h*?Ofl ~ J L S l ~ ~ W 5'2
QSLS?' tSfl'S1 Olebbf b21Lei;1 0'1
OSLS?' raoLmso ;:.*;LC ~Lb3'Il 0'5
(ur-069 (93611.4) ( ~ ~ I Q I H )( 1 ~ ('li1)
V3bV Al13313h AdJVM&Nq ,d*B173iJ 33\01615
7 - P I N WE2E-WRAPPED BUNDLESAMPLE PROBLEM RESULTS (Continued)
CHANIiEL . ~ ~ E s ~ ~ T s
CASE 'I' frP1N N I S E * H A P ~ A ~ P LP AEW L E Y OPTE OUIZII~~ o ttnk i~,uu,ii I
cnrcehEL ~ I F . , ~ ~ I L Y ~
CASE I' 7mPrb WRAP ~ A P P L LPFUILEP
EGUILe VCID
Q U ~ L ~ T BYR A C l l O N
OaOOO 0,000
O8l)0O 0,000
0,ooo 0,000
6,000 0,000
OSOOO 0,000
0 a 000 0;OOO
0,000 0,000
OaOOO oaooo
0,UOO oaooo
0,000 Oa00O
0~000 0 a 000
0,000 OaOOO
Osir0O 0,000
0;ooo osooo
OIO0U 0,000
Oefi00 OeOUO
OsOO0 0,000
0,000 Us000
oa000 00000
,
1) 1) 0 0 Oe01Q
0. UOO 00000
oa0ou 0,000
00000 OaOOO
0.000 oaooo
0,000 0,000
os000 OeOOO
o,aoo oeouo
OaOOO. 0~000
0.000 01000
0.000 0,000
OeOOO 0.000
00000 OIOOO
o,oo(! 00000
osooo OIOOO
osoou 0 I 000
oaoou 00000
0,000 OeOOO
D BUNDLE SAMPLE PROBLEM RESULTS ( C o n t i n ~ e d l
o
1959
1903
I bSS
01Tlu
, la30
I.riij9O
3mlb17
3,0819
1,092U
2,9403
,00300 00000 o 1471 2,9137
1)001)0 0.000 e lubd 2.6966
OoOOO 0,000 ,1552 2.79@8
0,OO~ 00000 ? IbD9 2.9622
00000 00000 01352 2, 6 7 8 1
Don00 00000 0 1299 209699
30093 OOOUO 1,!ic ( 4 J.(rlUJ
lo000 0,000 l1 8 IS lmO12b
,30000 - 00000 ll a b 2 3.09111
Ooi)OU Oo0i)0 .la21 3,0263
OoOOq OoOOO o 1492 3?~519
Om000 00000 ' olblb 2o91bir
OOOOO 00000 01471 il09UlO
0.00~ ooooo i0 6 8 2069b7
OmOOO Om000 s 155.4 2o8U21
0#000 00000 r1blZ 209080
0 ; 000 00000 (1136 i!m86Ud
L,~)OU ooooo ,~ecu 2.4970
C8000 OoOOO 1 . 3.UL87
O o ~ ~ O O OoOOO ole26 3o11317
OoOOO 0,000 o 1879 So0i!27
Om000 OoOOO 0 184,l 3 o 3591
00000 01000 ? l r \ l 3 3668
Oo00O Go000 o 1672 3ai)lbJ
OmOOO OoOOO ,lb'Jl 2o9600
. > a r
7 - P I N WIRE-WRAPPED BUNDLE SAMPLE PROBLEM RESULTS (Continued)
7 - P I N W IRE-WRAPPED BUNDLE SAMPLE PROBLEM RESULTS (Continued )
C( 5 1 ( 6, 1 2 ) * ( 7, 8)
*,OOObl -,0090E ,0 1277
' 0m21107 r0147b4 3 Z 4 5 &,2
ma21397 ~059Fi1 1 13 2 7 1
-.44i7l -,zos~~ ,19355
- , ~ ~ o J z -,20t102 . ,19216
r13Q2d ' rr4h75b a18199
#Oh219 wp3SHU1 ,1 5 1 9 2
rZ15S3 m12929 $52825
~ Z l t ~ l Z o 05230 ,43432
m19H'JQ m l 1.02789
r 17bOQ
D 12700
a 12ubC
.
a 2 1 3'46
lub5;3
D 17530
~7961t1
i;ssrqa
~bi)'lU?
-,i!OOhd e 12795 ~51386
r.dCfJ63 r l23i(b 037580
0m4751J -0231Jb , 19tU2.
-~3bl)& m~2Ob?O ,1027F
,13U3b ma07638 s lb9OC
.
.m05729
21552
,i?1 b 6 5
0,36277
alL61u
e05S74
s 13Hlll
05i!i)$7 .
,11271 1
c'?OObU m21557 BeOlb2l
o 17711 021762 ' e78019
r 12997 eZli29U 0qu'39! .
m 12524 ,17990 ,bdULL
-,2olSU elZd'i9 057 187
*e2G472 9 12414 i31413
mr~75S7 ae20051 6 19060
ofi36141 -eZ04149 ,19113
a13059 -rU7&21 e.14986
,0 5 5 7 8 e m55931 .;I l d 6 4
,2 1Qh9 r 13179 ,.s 18 3 b
,21&27 ,05817 i42UbU
,19358 eEl3I1 l;?l049
r 15 8 8 2 aZ1UUI 6.760'15
,133al 1726Q i40920
THERMAL MODELING
a ROD D n A I L 4082STEELCLAD
0.836 DIAMETER U02 FUEL
GAP COND = 1000 Btulhr-ft2
STEEL
3PACEH O A T A
$PbC€:h TYPE .KC, ! 2 3
LCCAfIbfi (x:/Ll ? 0- 1 3 ,729 ?013
SPACEY T Y P E 1
C W A ~ I ? E LD Q G CHANNEL DRAG CHANNEL O R b G C H A L b C & b8AG
Ni!, CCEFF, hoe , C C E F F I hO, CCEFF, NCg COEFP,
1 ,500 i! ,750' 3' ,750. '4 mT?o
5 2,000 b ZmCi3 7 2,009
SPACFH TYPE a
C h A h l h t L 0466 CHANNEL DYAG CHANNEL DRAG CHbhkEL DRAG
NO, CCEFFB NOd CcEFF8 kO! COLFF, COkFF,
t)
1 ,500 2. ,750 3 1750 . 4 8 150
I 5 Oe3O0 Q 0,006 7 OeOOO
N
d
S P A C F R TVPk 5
C*4AlhtL C d L C Ch4lVhEl. O k A G CHAhNtL 0HbG CRAhhtL DRAG
Vll, CCPFFc NO, CteFFa hue C(ltFFP NCe . COEPf
I o 600 a o,oao 3 O ~ O O O Y 0~000"
S lo800 b 1 ,000 7 leOOO
ROD f k P U 1 DATA
AOn TYPE O l A ~ A D I A L P O ~ E R S H A C T I C ~ UCP POWER 10 4 c J A C E h I C B A N ~ E C S (ADJ, CHAhhEL h0,B
NO, NCI, ( 1 FACTCR
a I f,OCOO 1,5000 ,bbb76 1 ) ,4166( ,41b7(
8 3 ) 4 ) a o e 0 0 0 0 ( - 0 ) m0,0000( -0) e080000( -0)
z i loorno isoooo ,!bb7( 1) ,41b1( a) 6 6 4 1 -0,ooooc -01 ~ o , o o o o ~m o ) Q O , O O O Q (- 0 1
3 1 1,0an0 ,500~ ,9bb7$ 0) ,4gbb9 ,2) ,Ql67( 3 ) mO,OOOO( r O ) -0,00000 -0) ~ O d O O O O- ~0 )
L ~ P L I C I T S C L L T I C h * I T H PRE8SbRL DRGP S P t C I F I E C
RL h \ I h G 3 T A a l =QC.M 1 " v L l C l E S t L U T f O h
OLTA F C Q I ~ P L l C I f:3OLUTICk
E X l E G v b L I T ~ Q P T ' I C * IL I P I T 20
I ~ I E ~ ~ ~ A I I. T F ~ A I ' I C ~ t . 1 ~ 1 ~ 20
C " Y ~ ~ C E . \ I . E FnC'rKQ?
~ X T ~ R ~(nwlr) A L , 1000
1hTtkhbL ( G I / * ) ,0010
G) FLC* (OF/?) ,FlOO
I
N
N 0 4 1 4 + U R EXPLICIT SULUIIC*
' 1 4 r r lNl,r 1 1 , 9~1 1 ~ ~ .UlUO 6ECONDa ACCELEQAllCh P A C l L 6 1,0000
COURAkT h U b 3 E Q ,4306 PW'ACTIOh OF LAST C ~ r k t i t -0,OOUO
X T E U A T X V ~c t v ~ c m r c ep u T o n a IIEW4TICS L I k l T S
VCLUPE t c a l n I O V / V ) ,0100 PRESSLkt l l E . k I T l O h 30
Ehlt!Al.Pv EiWUY (13Mm) .0001 Eh'EHGV t U u A ' 1 L h 20
rc(rI\l; C C a d t ; ~ l i c k 3
.pb.iCr.l.lED ~ l X I f . C , tit'TA 8 ,0200
b C I L 1 4 G P I U f h C , UCTA I d A S S ~ W L O$AWE AS S'JRCOCLEC
C C ~ O ~ C l l l irhl l l h t , G E C * C l R i FACICR 8 ,5020
C P E R A l l \ G CC41:1TIChY
3VSTE- P E t 3 9 i r f i E o 15.0 P S ~ A
INLLT Ehl*ALPV 8 1230.0 H l l l / L 4
AVG. * ~ S S r E ~ r c l T Yo ,010 MILLIO~ L B / ( n Q - S C ? l )
I~~.F.T T Z ~ P E & A T & ~I! € 0 . 0 DEGWtES C
AVG, t-EbT FLUX 8 ,120030 h I L L l O b 01U/(Hfi-30FT)
EXIT € h l r b L p v 91b.0 8 1 U I L 4
1 4 D I b I D L A L SLlCbbhhCL 1 E Y P E ~ A I U R ESPECIFIED (ChANhELmlEPPERlT~~E)
f I -lloo,a) ( a -iioo,o) ( 3 -1100.0) ( u -tioo,o) ( 9 -1o00,o) ( b -1000,o) ( 7 mtooo,o?
FLfl-S S @ L I l F C Q E i U A L FdES3kRt GRADIthT fCubhhEC-FLCb)
. THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
b;hL 1 i.b
FACTCH
. -
t ,0000
1 gOOOU '
0.0000
F u ~ c T I t l iF:R HEAT FLUX
H E 4 1 FLUK
PACTCH
1.0000
6,0090
0,0000
H A X I P L M EHROR
INTEfihAL EXTERNAL
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
i
UAlf OY/2I/Tb 0 lIHt 0 1
E N E l i l i Y EdLAkCk .m
P L C N ENERGY I N ,l O S 4 5 E t O 5 BTU;~~C
E N E R G Y ADDED , I J B O b E + O 3 BTU/SEC
FLOW EbiERGY CUT ,10597€+0'5 BTbf 8EC
f NEfiGY ERRPA m.85668E602 BTU/StC
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
EGUIL. VGIO
QUALXTY F c n C t t O ~
0.000 0,000
0,000 0.000
Or000 0,000
OIOOO OqOOO
0,600 0.000
OSGOU 0.000
O'eOOO Om000
01000 08000
0,000 0.000
OrGO0 OaOOO
01000 0.000
08300 00000
O,OO(r OIOOO
o0ooo 0.000
08000 OaOOO
Oe0,00 Om000
OmOOO OeOOO
0,000 0.000
0.000 0,000
06QO(' 0mOOC
0.000 0.000
GaOOO 01000
OeOOO 0,000
0.000 0.000
0.000 00000
, .
TIME .,
0,OOono SECOkD9 PAL:SSlrGE R
TE-PERATURE
15,'0 PSIA
C d P t i h E C fiESUL.19
C A.SE 1 FUEL THE MA^ hOOEL TksT
CHANNEL 3
.45
I *I ** I* *I *I * I
*+t+++t+++++t+t+++t
~t+++*tt+t+++++++++
*+tt++t+t+++t+t+t*t
*+ttt+t++t+tt+++t++
I+ttt++++tt++ttt+tt
*+t+tt+tt*+tt+++t+t
*+tttttttt+t++++tt+
*tttttttt*tttt+tt+t
*+ttt++t++ttt+ttt++
*+tttttt+t++ttttt++
*tttttttt++tt+tt++t
Ittt+tttttt+++ttt+tt
*t++t++~t++tt+ttt+t+
*+tttt+ttttttttttttt
*+ttttt+t++ttttt+ttt
*+t++ttt+++tt+ttt+tt
~+tt+tt+tt+t+t+++tt+
It*ttttt+++*ttt+t+t+
*+tttttttttt++tttt++
*t+tt++++t+tt+ttt+tt
a++t+++~+t+rttt**~t*
~t+t+tt+ttttttttttt+
*tt+ttttt+++t+t+tt++
Itt+tttt+tttt+ttt+t+
*+4tt+ttt+++++t+++t*
I * I I L L b I* I* **
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
. .
Cb!Ah!+.EL hE8lrLTg
C a SF, I FUEL TMEMAL wUGEL T E S T
. .
TICF a ,ZiI00(1 SECChD5 PPkS8LRE 8 15.0 PSI? O?Th FOR CHAhNEL 1
ECUILo
QlrbLITY
s.co0
0,000
PbOO!,
0,000
00000
0,000
O*OOO
00000
0,000
0,ooo
~ ~ 0 0 0
oo500
0.603
0,000
0*000
00000
0,000
0;ooo
30000
rJ.001)
3.0PO
OI000
poooo
Oo0OY
0,000
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued )
CHA'4kEL HE3l!L13
CIA? 1 FIJEL T H E V P L MODEL TEST
-
cast: 1 FLIEL Q H L M ~ LP - ~ O ~1 L~ 3 1 DAlt ,04/21/7b 0 TIM5 lZa,ZUrlU 1
CLAD
1099,8
IO9Pab
1099,U
1099,3
1107al
1 1 1 0 a'&
1113e?
1 1 178 1
1 1 19.5
1120.7
1 1 2 1 .e
112284
1122e1
1 l Z l a0
1119,i
.1117,C
1109.0
1108aL
1108a3
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
EXPLICIT SCLUTIONaaea
THERMAL MODEL SAMPLE PROBLEM RESULTS (Co7tinued)
C ? A L ~ V t Lk t 311~13
cA ~ E 1 F U E L T ~ S H A LNOEEL T I ~ B I D A T ~ O U / Z I / ~ ~ o ~ i r t LZ,ZM,LO 1
CHANNEL 7
050
qt +I *z 10 r t *r t ~ * I
*
I
I
t
t
(I
*
C
I
I
I
I
I.
I
I
t
THERMAL MODEL SAMPLE PROBLEM RESULTS (Conti nued)
D A T A COG CHANNEL 1
VCID FLOW
FRACIION (L8/8EC 1
0 a 000 l.dU34
0 ,0 0 0 1,4965;
0,000 t ,5082
d,OOQ 1,5201
OIOOO 1;53fL
0,oCIo !, 5 4 7 6
0,000 la5b2M~
0,OCO l ,578J
0 a 000 1,593b
0,000 1 6084
Oa01r0 1,6225
OeOOO 1,6558.
0.000 9,047?
0 a 000 l eb58?
OaDOU ,
I b682!
OaOOO 1p 6 7 6 i
0,OOO l ab(r29'
OeOOO 1,689L
oaooo !,705b,
0 0 OOU 1,7965
0,000 t ,b e 5 5 1
UaOOO I aA97L
00000 O a8942
0,UOO 1 s8ti67'
09000 8 a8757'
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
C:lA'u'itr'L .wt3ULT3
C.A St: I FUEL fhEMAL MODEL T E S T
FLOW ~ S F SL C $
(Le/sEc) (rbti/niw1a)
3,4504 ,b90P
3,4857. q b9bb
3m516d ,7052 .
3,5520 ,7i~?ti
3,5348 ;71h4
34blb8 ,7226
Sabu7T ,729~
3,6772 ,7309
3.7052 ,7405
307315 ,7457
3;7559 rn 7SOb
3e7763 el551
3e7Ybb ,7591
3,8160 ,762d
3,6331 ,.7660
S98U74 o 7b89
3,8595 a7713
f ,8b99 ,7734
3mM737 a7741
3iP5(1% . c7'4A7
388272 ab!lY
3 0828r? ~6120
38L1321 ~6128
308373 h 135'
3,8417 ,6143
THERMAL MODEL SAMPLE PROBLEM RESULTS ( c o n t i rued)
C*A:!.XEL iiEYttLT3
CASE. 1 F U E L THEMAL ~ O C ~T ~L S T
. . DATE 04/21/10 0 TlNt 12aZbadb 1
5, 1 ) nt 60 71
eCL413 ,012lU
eU3loO a08537
)O!lOU 8!lIb49
8OJOI4 eOA5Oh
802Yb4 gi1147M
a62910 #UlU5',
eOZhbU ,0145U
s.li2793 a U 1395
06270s - a U 1j 5 U
,025YJ #cjl2'lr(
0ozi(5e so~2c?7
,OiZY9 .~~104@
e02014 8 ill ULjb
;u!ors ,OO'~OII
e0878b eUtJt5B
eO1SUU eU375U
eOldOt #iIObSlj
eb1AU7 8U05S4
800957 ,UUUE(I
eOU857 . o(r3Ui"i
pOGSb7 ;0ttdd~
,OU400 suUdil~
@UC$bb au313r
e OJlZO e~OI)bO
THERMAL MOt3EL SAMPLE PROBLEM RESULTS (Cor;.tinued)
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
0 4 T A F W O r EXpL:CIT SLLUIION~~~,
. 3 Z a1689
? 1.9597
31ahbbb'
31r343.l
.30au553
3Ga53'31
.3Oa04l0
29r5391
'ita9723
ia. 3 9 a 2
27a76O3
.i7, 1.197
ib.02lb
25.7253
2Q.9oS1
Zba263U
i3aS22S
22.7762
2Za~tbU
i 1 a25U8
THERMAL MODEL SAMFLE PROBLEM RESULTS (Conti n u ~ d )
TEMPERATUPt4 ( F 1
[RELATIVE HAClUY(W/H))
T( 1) T( i l
.
1( 3 )
(1a000) ( ,bt7) ( ,333)
lOQ9.8 109qa8 1099,8
1099#b 109F.b kU90Db
1099a5 iOYYcb 10'49,b
1104a9 11OTe0 llOBab
1299,7 1900',3 2460a0
136609 2321ab flJ91,2
iUJZ.9 Lbd7a5 3h58,9
iU90e1 30Szab 4151a4
1532a6 ' 32J4.3. 4370a5
1534,l 3237.4 U373,b
155581 3239.1 &47ueh
1535,9 3219,2 U373a9
.:.1504aZ 30bOmk U157.3
luJ9.8 269CB1 3b6dab
1375ag 253&,1 Jl0Ua&
1309,0.,i,, 1 9 7 3 I S ZY95*3
iiiuae iiibr9 1114.5
l l ~ 8 a O ii-GF,,a i l C C e i
I i o e a s Ilo8.e 110d.b
v so il*
+++++++++++++++++++*
+*++++*++*+++++++++*
++++++++++'+++++++++V
+++++++++++++++b+++?
+++++++++++++++++++*
....................
+++++++++++++++++++*
....................
++++++++++++++*++++r
. ....................
+++++++++++++++++++*
++++++++++++++++++*
. ++++++++++++++++++*
++++++++++++++++++I
++*+++++*++++++++**
++++++++++++++++4+* +++*+b+*++**++'+++*v c,
++++++++++++++++++* +++++,+*++++++*++*. t,
+4++++4+++4+++++++r . +++*+b*+++++++++*t.
++++++++++++++~+++v ++++++++++++++++ t
+*++**+*++++*+++++* +++++b+*++*+++++*** 1
** I* *P 24
** v * LC C* ** * * *Iv * 1, +*
Sk' Sh'
k 13~slVY3 1 73YNV!43 .
TEERMAL MODEL SAMPLE PROBLEM RESULTS (Cortinued)
rib
a 16
a 15
a15
0 19
s 14
a19
a 14
a 13
o 13
812
a12
a11
,it
a 12
t i3
013
a l 4
a10
a15
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued )
C M A t I b i i r L iiESULT$
C &St: 1 F O E L T F E ~ A L~ O C € L I L S I OATE U4/Zl/70 0 TlMt 12a2Ua30 1
C I V E H Y I O h CnC'3SFLb* d E T * E E l u b O J A C L N I C H b W E L 3 e k ( 1 e J l o ( L 8 / S E C r F T ) ,
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
'$:. a !+
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
VOID
f f i A C 1 ION
oeooo
OeOOO
oeooo
0.000
0.0~0
OeOUO
OeOOO
On000
ceouo
(i.000
OeGclO
0,000
0.000
6e000
OeU00
CeOOO
0 e 000
o e o o o
ceooo
0.000
UaCiriI
a,. 0 0.0
0 I 000
u.000
0.OtlO
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued )
nxraL
0.0 -
2.u
U#O
-,
--
---
ha(!
A,.,
10,o
12.0
t4.Q --
Ih,O
1 --
I
23.9
22*ll -
5: %fl,OD
--
---
,???.I)
2!.,(l
. d : ,
32,\1
---
$.~,n
35,o
JII,O
il'!.O -
$do()
4h.O -
r
THERMAL MODEL SAMPLE PROBLEM RESULTS (Continued)
CLAD
DAMPED U-TUBE MANOMETER OSCILLATION
One of the most important new c a p a b i l i t i e s of COBRA-IV i s t h a t of con-
sidering density gradients.approaching steam-water interface conditions.
The. modeling of a U-tube o s c i l l a t i o n i s an ideal way t o demonstrate t h i s
capability. The model a l s o employs the code's new a b i l i t y to a r b i t r a r i l y
block l a t e r a l flows in t r a n s i e n t calculations.
The geometry chosen f o r the U-tube i s i l l u s t r a t e d in Figure G-3. The
flow area of 36 square inches, the 6-inch gap and the axial node length of
3 inches were chosen to provide the same area f o r l a t e r a l flow as f o r axial
flow and thus a uniform velocity a1 1 the way around the "U". A1 1 gaps
except the f i r s t two a r e closed via input. The only damping present consists
o f a drag coefficient of 2.0 placed a t the 14.4-inch level in channel 1 .
One channel i s i n i t i a l l y f i l l e d with saturated liquid and the other with
saturated vapor a t 60 psia. The calculation i s carried ntrt t o 4 . 0 seconds
which allows two f u l l o s c i l l a t i o n s (natural period i s 1.92 seconds.)
The Sntermediate output f o r t h i s problem shows "surface level" in
addition t o the normal t r a n s i e n t output. This i s a post-calculated quantity
which i s determined by integrating the previously calculated density and i s
not used elsewhere i n the code. Refer t o the channel 'results page f o r the .
actual density gradient calculated.
LIQUID
E X P L I C I T Stil.UTI(:h N I T W I k L E l F L t j h S 8 ~ t C I . F I t D
S T d t . ~ l l r . r ; 9 1 6 G r Flc1.lr.l Z E H C FI!LL#IWY
CHANNEL b E h 6 T H
N u r b E R UF n x 1 A L h t i L t 3
AXIAL hGDE L t h C l h
TOTAL T U A h 9 1 k ~ l V l H k
NUMEER OF T l n c S T E P S
~ O M I N A L l l k t BTLP
ACCELEHATIGh FAClGW
P R A C T I O h OF L A S T C b A t v L t
ITERATION L I H I l S
PRESSbfit I l t M P l l b h
EhkNGY L L U A T I L h
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT ( ~ o n t i n u e c )
~ ' i A t l \ t . L ktSt1,L T $
C4YE 1 IYI.~-HE'HM~L lJ-f JHE l l r h f J e c ~ ~Dk C ~ C ~ ~ S T R A T ~ O N
5 104 ,001
SiikF.5CE L C V t L 72.09l)
3 ,aotrln , ~ \ !
SJdFsCE LEVEL 72,UflU
Z ; , sOdll2 ,1123
SUVPhCE 1.Evt.~ 72.01~~
2 .o:;llc ,n53
S ~ J ~ F L CLEi v E l . 72,udi)
2 ,OIJ1?U ,0110'
Stj9F:CE l.6 VkL 72,OClU
2 8Uc?tZ5 e@~S.
E 1 ~ ~ F a CLEVEL t 72,000
2 . eI?f!13a roll8
JdriF:,Ci: LEVEI 71.@999
2 ,03135 ,05d
Ssi'F i C E I.k.'ubL 71,099
2 ,nil iu o ,051
ZJugcCt LEVEL 71,998
2 0111 U 5 ,052
Z U ~ F A C EL E V E L 71.996
2 ,31)15l ,052.
S0J+FbCE CEVKL 71,993.
2 eOOlSh ,051
Slik'caCt C ~ V F L 71,9P9
2 ,aOlbi! .05 1
5~'irFbl;k. L ~ V E L . lle9ilu
2 a0Qlh9 e OSq
S U ~ F A C ELEVEL 71,970
2 ,00lj5 ,ous
S C E F - A C E LEVEL 71,9tb
3 ,o:llYZ roue
- f U h p A C E LEVEL 71 r95U
3 ,00189 ,04A
SURPJLF LP.V.EL 71,938
3 ,On196 ,OU7
S ~ ~ F A CL E V E L 71.C19
3 ,602n3 ,oh6
SU9F6CE L E V E L 71,897
.3 *(lll?l 1 ,946
S k k F L C E CEVt'L 71,872
3 ,0lid19 ( 0 45
SllgFACF. L E V E L 7l.,OU3
3 ~03227 4 005
$UQFICE LELKL 71,el0
n ,00236 ,005
SukFbCE LEVtL 71,775
If ,002U5 , 006
SU2FALE LEVEL 71.737
u , f?I!ZSU ,047
SvkF&CE L E V E L 7 1 ,C9h
U ,1!?2hU ,lltlg
SukFltF LEVtL 71.t52
5 ,u327~ ,65 1
S U q F A C E C?.VCL 71 i6nS
b ,(:02AU ,055
2.92 F\i'Cf 1. fib? L , '7 1'; 5 5.5
? ,ii0295
- . . -
-. . ,059
DAMPED U-TUBE MAYOMETER PROBLEM OUTPUT ( C o n ~ i n u e d )
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Con.t.inuec)
D A T A F O f i CbANNtL 2
O'tl
bWb4
O*JV
(lei'(
0*0f
U*LG
O'bS
F * 1s
CWHo
V*S')
(#*ell
OObF
(1
* 0F
Il*%F
(I * Il r
o0L?
U'n?
(I* ti!
o'w1
C.51
@'PI
6
(1 *9
0's
3'4:lZ
DAMPED U-TUBE MANOMETER PROBLEM.OUSPUT (Continued)
. . c a r 4 Fllllr EXPLICIT YCLU~I~I~,,.,
voLurt
ERROR
j\*d
3 j> *
34, o
3J2 ,0
' I S ,9
.J ,11
fl
SI*$
s 11
>r.:1
.
.:i
b t~ ,;r
s*0
hh i)
bq,~)
7249
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT ( C o n t i n u e d l
~ N ? H A L ~ YTkn*Pt:AAIUHE CL'~t811Y
(i3TU?L9) . 10tG*F) (LB/CU-FT)
1 174adO
250.20
ZSt,C0
2 6 1 #r;O
a
9.7,H'i
.
8 -a
HC I, 21
217,ilum7~
I I-!i,JSh4L
0,00300
o,oa~oo
0 i j 1) '10 9
,
$1 i j 0 '1 u I)
O.dli31)O
O,?)IIJOO
,
0 0 1) -1 0 0
d,uOJuO
,
0 3 1) ) 0 1)
ii, 3 3 ii)u
, 5 3 .j I, 9
,
0 1) 0 .) 1) 0
0,Du3JO
i1.001)Or)
l J , i)1)1)03
i1,ir1)$3u
0,301)01)
0 q 1) 0 :) 0 0
I ) * 1)I)l)l~O
ci,oirn~i)
9,001)OO
rr,1)01ion
DAMPED U-TLIBE MANOMETER PROBLEM OUTPUT (Continued)
.
1 : I. .; , TI*$ CUUHA~;T VEL~~CITCYOMPL~LA~S TCTAL
~31r-3 STEP ht~I:~~aii) L'/DX V/DY (b/cr*v~oY) TIME
-
3'JiiF .ALE L E V E L
5uGFdf:E LEVEL
3 ,00~22
,
05.253
u u ,gu,?
,202
26,730
27,1u5
Ll?,b3
!J 1 ,tiGu27 ,1 9 9 ub,~Z
3 i J ? c 4 ( . E LEVEL 3'7,544 32,'03h
a1 ,C1)02i! ,191 Un.hZ
SJCFJCS LEVEL 3q.13a 32.~05
- i ~ . ,09a2? *'?I) 1 47,ilJ
3 Q G F b C E LE>EL 3*,731 33,250
Uo , ,6Cd2!3 ,200 Ohaq5
S c . ~ n ~ ~ C E C t ~ E3 3L , 3 ? 5 33 0650
20 *09?r29 e 1 917 L O , ii8
S'JUFAC' LEVkL 37,923 54,057
3 '1 ,(!GQ~O c 1qR rib 2 5
SiJxFnCk LEVEL 37,519 JUa'JbO
13 *UilbjS a 1J Y U5aQ6
3 G G ; ~ i t LLVEL 37,115 3 4 , daU
3J 9 J O U 3 5 a 199 OSe72
S<JKFACE LkrEL 33,711) 35.2i8
3Q n v t ~ j 7 144 US.SZ
S . l + ' J i F : LEVEL 3ha305 55a673
39 .Oil037 *to1 U5,ii.
3,JrFJCE LEVEL 3SaF99 3t,U79
i .? .CljU36 *21)t 'r 5 ,9'7
52-FnC E L E q E L 3 5 , ui*5 3 6 ,U82
!1 a Oou3@ o.198 . 05.41
B J ~ F A C EL E V E L 35, 0'40 3b,8F3
S7 0!!1ruZ ,1 9 7 UU,QU
S ~ ~ ~ F J LEVEL C E 34,693 37,iku
51 *00UUll sl9Q u~;dl
S a J ? F i r . 5 LEVEL 3d.291 37 6 8 0
57 ,f l i l u u 7 0 197 q u au 2
~ ; J U P ~L ~ E VE E L ss,riae 35 948
Ib ' *OOu50 ,197 4~,11
i u i F . t i E CtvF.1. 33,485 3rruYO
5 ?I a Ouu5 1 ~ Z d u 00 32
S U ~ F A C EL E V t L
; ;j . a:'~so
S:JhF.+CE L E ~ ~ L 3 2 . h 7 9
33rOH1
,a91
3AaAQ!1
391295
44.0~
2s a00893
"
S ~ ? F ~ CL EE V ~ I 2 7 ,.sou
DAMPED'U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
D A T A FOR CHANNEL 1
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
*;. I
TI*f: COu44~1 V E L L C I T Y COpPChEhTS TCTAL
rslts Y TE? YUH~EQ li/!?x v/ 0 (b/cx+r/oY) TlrE
59,531
3.04
ENTHALPY
(HllJ/L9)
I 174,ni)
23r1,2n
(?50,20
<50,20
253820
251J,20
250, lo.
250.1Q
259.19
250,ZO
250.20
2so.20 .
250,2U
25UoZO
250,zO
250 20
250 022
253.52
C51.1~0
25.9a 5 2
360sSR
977.07
Ill'Ja00
I 17i1rSr)
117Uebl'
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (continued)
oarA F R C * 'XDLIC[T SLLUTIO~~,,.,
.
;
~ ! J ? F : \ C E L k V=L 22.211. 40ob74
49 .uo714 .a04 28 09
S 1 J i F 4 C F LEVEL 22.638 ~ 9 02116
27 001)795 ,208 Zh.8d
3iJ*F A C E L E V E L 23.0b6 ~8.a17
; 5 ,:!rlh90 ,209 2q.66
~ u ~ G . I C F C5UE!. 23.091 uP,.3q1
; 'j a113hC12 0205 29.67
~ ~ . J * F A C LEE v < L 23.911 U?.971
LI cj ,1:fi670 ,205 33;irl
SU.+FACF L E v c L ZU.330 07.551
. j .fi)Oh70 ,203 39,Ob
S ! J ~ F ' A C EL E V E L 2U07U9 07.132
44 00Obhb 0202 '39,ZL
j1.J$iciE LEVEL .25.1b8 b6a712
'i4 00flhhZ 0202 3:).Ul
3UGFACE LEVEL 25.598 lJb02Q2
a5 ' .iJI!b52 ,rob 31.10
SIJGFACE LEVEL ' 26,?0B 115.871
33 r UOhuO .a@A 31 .au
SIJHFACE L t v k L 26,025 45. U53
ti 3 o'JOh5h oL"J2 31.42
SUGk A C E 'LEbEL %o.e37 u'5.0Crl .
.'r3 . .01!b35 0 21) 1 31.50
S1JkFpT.E LEVEL 27.2d8, ::&.hill
LJ m .I)Gh3S .20 1 31.h2
SUSFdtE L E V E L 27.bh0 UUa217
2 11 000632 0201 31.70
SUFFJCk L E V E L ZPoC7U &30i?03
rt ir *ll(!b!l .291 3.! # 77
S1.JMFACE LEVkl. 2*.ui(8 U3a.3R!4
2u 000h27 0263 32.11
JiJsFICE C E V ' L 2!!,93,? 02,975
" J .l1iae!p ,i l i h 32.*'5
SbUF A C E I t v ? ~ 29,315 ~~,SOII
,!
S ~ S F A C EL E ~ E L
2 1
.OOhlU
.\J?b 1 d
2s. 72.3
0203
.2cIO
.
42 l 5 c
32.22
32.54
S!JkFACE. LEVEL
.-..' 3J
., - C L . -.
;qnhi3
30.129
a'.*=.
. 198
4
Ula7U5
-. *
32.43
..
u3 ~Jfiblb 199 3i?ouO
S~JQFACF. L E V ~ L 3 a . u ~ ~ uoa'92q
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
VOID
FliACTION
0,000
0;ooo
0.000
OaOOO
oaouo
0.000
OaOOO
OaOOO
eooo
a009
199
,660
,679,
,957
,986
1996
.9F9
1 a000
1 a000
1 aouo
i, 0 0 0
I,OOO
1 a000
1.000
l a000
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
C*4*rP.EL YtSlILT3
CASE 1 IS;IlnEfi!*&L Cl-Tt;t)E ?dAl ~JI;ETER DEVCh,$TRbT I O N OAIk 04/22/76 0 Tlfik 12,1~,31 1
SIJIFICE
;.
9U4FbC.C L E V E L
1
LElvEL
.
' ? 3 ht L
32.963
33.385
0 198
38.5~~
$9e~*h
32.53
i~1 eD1'414 ,199 3's 3 9
~ U U F A T EL E V E L 33,739 J2aOA1
-11 ,05617 a 139 322“
S ! J ~ ~ ~L EC VEE L 3u.194 37,676
31 eOJh?O a19h ' 3 2 0 10
a 036*3
39.805
e 195
e lf45
32 059 .
30.17
20.83
' S I J ~ F L C E LEVEL UO.ZOS 3l.hhl
35 mil0670 e 175 29.39
3 0 4 F i C F LEV.EL OiJ,hOZ J l rZh2
S LA .00673 a 199 21.03
SllsF&CF:LEvEL ol,(ii~l 3 n ,Hh+?
33 .OOh70 .202 2 0 ,?b
3USF i f E LEVEL ul13Q4 3 0 a bhS
32 aOo67b ,197 29.36
s u ~ r ~ rL EE V E L
3 ;,
SURFACE L t V E L
.onbe5
s i @7hu
bL,l7R
,
195
30,074
Z91bfiU
26e91
.-
'
F Q f i CHANNEL 1
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
3 ~ b7r r r ~ E X P L T C I T s r L b T I a N . , , ,
rle,il4u
.
191
25.e40
25,415
23.?7 ,000u
VOID PLOW
FkACTION (L8/SEC)
3.0 ?,0(10 ?, 0600
3.0 0,000 ,511b
h , fl 0,000 -14e21b5
9.0 0,000 alUe23'Jl
12.0 O*fIOO alua2519
l5.n 0,050 -18,2696
IA.11 0,000 .IV*Z~bd
Z I. n 0.000 al4e303J
2 4 ,0 OaOOO ~lur3188
c 7.0 0,030 -14r3.530
~ I ~ * O OaOOO mI4aJd57
55.4 0,OOO -14,3567
.5 h ,il 0,OoO -14,3650
$1'>
.,
s q , 11
J C il
0
0,006
eOu0
,OUO
-14r3728
-1uo37bO
-1uou53
I I ~ I . ~ - ,003 -1j.eoul
51 .I1 ,040 -7,4118
'>;I .n 0U85 -,9073
5 1.1, ,939 0,2796
bib.
.--.
,983 -,bCCL
h{.d ,996 ,028 7
h h . ij 1 0000 -a0291
hQ, n 1 to00 a.0292
12.9 1 ,000 -,OZVL
DAMPED U-TUBE MANOMETER PROBLEM OUTPUT (Continued)
BNWL-1962
DISTRIBUTION UC-32
No. of No. of
Copies Copies
OFFS ITE
/
A. A. Churm 5 U.S. Nuclear R e g u l a t o r y Commission
ERDA Chicago P a t e n t Group ATTN: C h i e f , Systems
Chicago Operations O f f i c e E n g i n e e r i n g Branch
9800 South Cass Avenue D i v i s i o n o f Reactor S a f e t y Research
v Argonne, I L 60439 Was h i ngton, DC 20555
27 ERDA Technical I n f o r m a t i on E. Davidson
Center U.S. Nuclear Regulatory Commission
2 R. J. B a r i b o l d i D i v i s i o n o f Reactor S a f e t y Research
ERDA Chicago Operations Office Washington, DC 20555
9800 South Cass Avenue
Argonne, I L 60439 S. Fabic
U.S. N u c l e a r Regulatory Commission
U. S. Energy Research & Develop- D i v i s i o n o f Reactor S a f e t y Research
ment. A d m i n i s t r a t i o n Washington, DC 20555
ATTN : A c t i n g Di r e c t o r , Program
Analysis & Evaluation D i v i s i o n S. I s r a e l
Idaho Operations O f f i c e U. S. Nuclear Regulatory Commission
550 2nd S t r e e t D i v i s i o n o f Technical Review
Idaho F a l l s , I D 83401 Reactor Systems Branch
H. C u l l i n g f o r d Washington, DC 20555
U.S. Energy Research and
Development A d m i n i s t r a t i o n P. N o r i a n
Engi n e e r i ng Design and Component .-
U. S Nuclear ~ e ~story u l Commi s s i o n
Development Branch D i v i s i o n o f Technical Review -
D i v i s i on o f C o n t r o l 1 ed Core Performance Branch
Thermonuclear Research Was h i ngton, DC 20555
Washington, DC 20545
J . Hopenfeld 201 t a n Rostoczy
U.S. Nuclear Regulatory Commission
U.S. Energy Research and D i v i s i o n o f Technical Review
Development A d m i n i s t r a t i o n Reactor Systems Branch
Componen t s E n g i n e e r i n g and Washington, DC 20555
Development Branch
D i v i s i o n o f Reactor Research D. F. Ross
and Development U. S. Nuclear Regulatory Commission
Washington, DC 20545 D i v i s i o n o f Technical Review
B. Twining Core Performance Branch
U.S. Energy Research and Washington, DC 20555
Devel opment Admi n i s t r a t i on
Systems and A p p l i c a t i o n 2 Union Carbide Corp.
Studies Branch Oak Ridge N a t i o n a l L a b o r a t o r y
Division o f Controlled P,O. Box Y
Oak Ridge, Tennessee 37830
I*, Thermonuclear Research
Washington, DC 20545 M. Fontana
J. Wan t l and
.(
No. of No. o f
Copies Copies
L. S h o t k i n 2 A e r o j e t Nuclear Company
U.S. Nuclear Regulatory Com- ATTN: Mgr, Systems A n a l y s i s
mission Branch
D i v i s i o n o f Reactor S a f e t y 550 2nd S t r e e t
Research Idaho F a l l s , I D 83401
Washington, DC 20555
2 A e r o j e t - N u c l e a r Company
L. S. Tong ATTN: Mgr, LOFT F a c i l i t y
U.S. Nuclear Regulatory Com- Division
mission 550 2nd S t r e e t
D i v i s i o n o f Reactor S a f e t y Idaho F a l l s , I D 83401
Research
Washington, DC 20555 . 2 A e r o j e t Nuclear Company
ATTN: Mgr, PBF F a c i l i t y D i v i s i o n
A e r o j e t Nucl e a r Company 550 2nd S t r e e t
ATTN: Mgr, Systems S a f e t y Idaho F a l l s , I D 83401
Research D i v i s i o n
550 2nd S t r e e t 2 A e r o j e t Nuclear Company
Idaho F a l l s , I D 83401 ATTN: Mgr, Reactor S a f e t y
Technical Support D i v i s i o n
2 A e r o j e t Nuclear Company . 550 2nd S t r e e t
ATTN: Mgr, Model V e r i f i c a t i o n Idaho F a l l s , I D 83401
Branch
550 2nd S t r e e t 2 Aero j e t Nucl ear Company
Idaho F a l l s , I D 83401 ATTN: Mgr, Thermal Fuels
BehavioF Program
2 A e r o j e t Nuclear Company 550 2nd S t r e e t
ATTN: Mgr, A n a l y t i c a l Model ldaho Fa'lls, ID 83401
Development Branch
550 2nd S t r e e t 2 A e r o j e t Nuclear Company
Idaho F a l l s , I D 83401 ATTN r Mgr , Semi scal e Program
D i v i s.ion
2 A e r o j e t Nuclear Company 550 2nd S t r e e t
ATTN: Mgr, A n a l y s i s A p p l i c a - Idaho F a l l s , I D 83401
t i o n s Branch
550 2nd S t r e e t 2 A e r o j e t Nuclear Company
Idaho F a l l s , I D 83401 ATTN: Mgr, LOFT Program
550 2nd S t r e e t
2 A e r o j e t Nuclear Company Idaho F a l l s , I D 83401
ATTN: Mgr, Experiment
S p e c i f i c a t i o n and A n a l y s i s P. L o t t e s
Br'drlch Argonnc Nat iunal LdburVatory
550 2nd S t r e e t 9700 South Cass Avenue
Idaho F a l l s , I D 83401 Argonne, I L 60439
No. of No. o f
Copies Copi es
R. Henry A. Gopalakrishnan
Argonne N a t i o n a l L a b o r a t o r y E l e c t r i c Power Research
9700 South Cass Avenue Institute
Argonne, I L 60439 3412 H i l l v i e w Avenue
P. 0. Box 10412
W. W u l f f Pal o A1 t o , CA 94304
Brookhaven N a t i o n a l L a b o r a t o r y
Upton, Long I s l a n d , NY 11973 J. T. A. ~ o b e r t s
E l e c t r i c Power Research
J. McCluskey Institute
Brown & Root, I n c . 341 2 H i 11view Avenue
P. 0. Box 3 P. 0. Box 10412 . '
M. P l e s s e t C. S u l l i v a n
Cal i f o r n i a I n s t i t u t e o f E l e c t r i c Power. Research
Techno1ogy Institute
Department of Appl i e d Mechanics 341 2 H i 11view Avenue
Pasadena; CA. 91109 P. 0. Box 10412
Palo A l t o , CA 94304
M r . J. Block
CREARE, I n c . General E l e c t r i c Company
Hanover, New Hampshire 03755 ATTN : Mgr, Reactor Systems
Safety Development
Consolidated Edison Co. o f San Jose, CA. 95114
New York, I n c .
ATTN : S e n i o r Vice P r e s i d e n t General E l e c t r i c . C o m p a n y
4 I r v i n g Place ATTN: P. McGee, Mgr o f Design
New York, NY 10003 Engineering
Breeder Reactor Development
G. B. W a l l i s Operation .
Dartmouth C o l l ege 310 Deguigne D r i v e
Hanover, NH 03755 Sunnyvale, CA 94086
A. E. B e r g l e s N. E. Todreas
Department o f Mechanical Department o f Nucl e a r
Engineering Engineering
Iowa .State Uni vers it y Massachusetts I n s t i t u t e o f
Ames, I A 50011 Technology
Cambridge, MA 02139
C. T r u e s d e l l
Department o f R a t i o n a l D. M. Chapin
Mechanics MPR Associates
Johns Hopkins U n i v e r s i t y 1140 Connecticut Avenue
B a l t i m o r e , MD 21218 Washington, D. C. 20036
J. Chen H. Estrada ,\
Department o f Mechanical . MPR Associates
Eligineet.irlg 'I 140 C o n n e c t i c u t Avenue
Lehigh U n i v e r s i t y Washington, D. C. 20036
.Bet h l ehem, PA 1801 5
P. Lax
C. W. H i r t Courant I n s t i t u t e o f Math Sciences
Group Leader T-3 New York U n i v e r s i t y
Los Alamos S c i e n t i f i c 251 Mercer S t r e e t
Laboratory New York, N. Y . 10012
Los Alamos, NM 87545
S. G. B a n k o f f
K. D. L a t h r o p Department o f Chemical
A s s i s t a n t D i v i s i o n Leader Engineering
T-Division Northwestern U n i v e r s i t y
Los Alamos S c i e n t i f i c Evanston, I L 60201
Laboratory
Los Alamos, NM 87545 B. I. S p i n r a d
Department of Nuclear
P. G r i f f i t h Engi neeri.ng
Department of Nucl e a r Oregon S t a t e U n i v e r s i t y
Engineering Corva 11 i s , OR 97331
Massachusetts I n s t i t u t e o f
Technol ogy S. L. Lee
Cambridge, MA 021 39 Department o f Mechanical
Engineering
K. F. Hansen S t a t e U n i v e r s i t y o f New York
Department o f Nucl e a r S toriy Brook, NY '1.1 194 -
Engineering
Massachusetts I n s t i t u t e o f K. M i l l e r
Technol ogy Department o f Mathematics
Cambri dge , MA 021 39 U n i v e r s i t y o f Cal i f o r n i a
a t Berkeley -
Berkeley, CA 94720
No. o f No. o f
Copi,es Copies
Corpora t i o n
J. E. C a s t e r l i n ' e Waltz M i l l S i t e
Department o f Chemical P. 0. Box 158
Engineering Madison, PA 15663
Col umbia U n i v e r s i t y
.New York, NY 10027 C. L. Caso
Westinghouse E l e c t r i c
' C. Morgan Corporation
Babcock & Wilcox P. 0. Box 355
. a P. 0. Box 1260 ' P i t t s b u r g h , PA 15230
Lynchburg, VA 24505
,. ONS ITE
... P. L o r e n z f u i
B e t t i s Atomic Power ERDA R i c h l a n d o p e r a t i o n s O f f i c e
Laboratory
P. 0. Box 79 B. J. Melton
West M i f f l i n , PA' 15122
4 Westinghouse - Hanford
G. diLauro J. Hanson
Combustion Engineering, I n c ,
A. E. W a l t e r
Windsor, CT 06095
3. F. Wett
J. M. Yatabe
L. P. Bupp
E x x ~ nNuclear C o r p o r a t i o n
75 B a t t e l l e-Northwest
21 10 Horn Rapid.s Road
Richland, WA 99352
J. M. Bates
S. H. Bush
G. A. Sofer N. E. Garter
Exxon Nuclear C o r p o r a t i o n
R. J. Cena
21 10 Horn Rapids Road
Richland, WA 99352
J . M. Creer
No. o f
Copies
ONS ITE
J. M. Cuta
D. E. Fitzsimmons
C. R. Hann
B. M. Johnson
D. K, Kreid
W. W. Laity
C. A. McMonagle
D. E. Olesen
M. G. Patrick
L. T. Pedersen
n. S. Rnwe
C. W. Stewart
A. M. Sutey
D. S. T r e n t
C. L. Wheeler (50)
R. D. W i d r i g
F. R. Zaloudek
Technical I n f o r m a t i o n ( 3 )
Technical Pub1 i c a t i ons