Cobra 20XTR Bedienungsanleitung

-WOSUB-A SUBCHANNEL CODE FOR STEADY-STATE ANDTRANSIENT THERMAL-HYDRAULIC ANALYSIS OFBWR FUEL PIN BUNDLESVOLUME IIUSER'S MANUALbyL. Guillebaud, A

viii2.6.11 Listing of Subroutine CONTI 2282.6.12 Listing of Subroutine SWEEP 2322.6.13 Listing of Subroutine SOLMAT 2412.5.14 Listing of Su

-66-1.4.7 Listing of Subroutine SWEEP-- ---- -- -- -- ------ -----

-87-SUPOUIJT I.JE qwEcPC :- SCHEME FOR SUPCHiNNELcCLCJLATION IN STEA)y STATF A TNISIFEJTL0OGICAL IImF.V,IF.--VIOCIP(20) .T-E(20),FQE(20),iI(2n},iP~-N(

-88-CaLL \'"-X 0000100001S A r,=. 00001I~~~~~~~~~~~ ~~~~~00001i 1O=1',v T CJ. 0 0010r103 I=l iC'tLC~ 00001t~aAI, (t ) 00001aJIN=i

-89-SAV V (I) =S Av; 'XT ( I ) =x TarJ ( I ) =I-I EJ ( I ) = I WI Cf!T TIN JF.4 \/ =o4 V/S AVC TEPAT TN C-E,'EC)n 12I=I - I&quo

-' I00001,a jAP T-JI A f ) oorJLT ( I ) =_ .r) ( °T D ( T ) = -L ( I n o oEG S ( i ) = I ,LT: .Io TEPE. iV(T)=I ii) O OlJ TFL~~~~~~~~~~~ii ~~~F

-92.-(' 16 () J= ,J Ar X 000 .'Zin -! ' r, T = -rP,T+. f=)P( (v) , ~ T I ~ 000'P I k 2=2 000'I T=' oin000'I TI&ap

-92-1.4 . Listing of Subroutines VMIX and IXIN--- - -- -- ---- ---

-93-SJU2ROJTT,'IE VWOIX-.:--:I G'NE.AL SChEvmE FOP vAO(n IXIN.V, IN SJbCHAN\JFLSL)OG C L ,!' ~1D TTCOmON/rEJ/AOD(2h) ,Po(4.20) ~F(5 I02

3a~~~~ ~~~94-, Ia ) L T;' t. Is I rl(i . 1Al '; MO)0F-LIF, ] C- (i icCq "N)q;? i-' TC ,-,! .i' .^ .,, / ·i T /)

-95-Fe E:' T UOLGOj TO 101" *** fIVEPTF)O L',i SLIP F.CTORFITOY SLIP- ( I J)C -T, ISVEPS"AL SLTP MODElt!.C7=.?11 SLIP(I)=1-./A&a

ixLIST OF FIGURESFigure Page1 Input Deck Structure After First Control Card 62 Input Deck Structure After Second Control Card 83 Input Deck Struct

-96-1.4.9 Listing of Subroutines SPLITD and SOLMAT.. - - I- .-. . ._.. - - - --.-- --·- -I -----I I -

-97-SIJ3PUITTNE SDLITO(J)i *'' S:'LITS TrHF IIVFRTEF) 'N,! INrT VADOR 4N;) LIutlID FnLOW. FVAL!IATES° C '~:' >T &apos

-98-C :"::' t-VFi FL,, FrNTr4ALv))'A 1 =1 ,*cJ %= .,c i ' ( )1-' S^VE=SVE+i,,,t(JJ)*O K(IK)rr) TO 1!1n ~,cAVE=S V-+LI(I) &quo

-99-SUOl.PUTriNF SOL '4T(TW. .4)r *" ZSOLVES THE MATOIX FOtJATTN 0r CSS FLOWSr)T'.4FNSION rw3(0,4) .W(PO) S() .VEC(40)Co-kMON!/GE0Mv/JO

-100-1 .4 .10 Listing of Subroutine CONTI

-101-SUPOtJTr riF CNTIr 0 SOLVES COJNTTt)ITY FOUATTDNS IN EACHt F TrE SURCHAJ\INELSLOGICL 8,l FPV, IFREV, I r4, ILnCCO N l /UNT/I'DZ .DT V7EO, N

--102 -t JV=U .O=j. v i.+: / /AI F- I Ii A n ) "= q + V Z r)OVD M'"F '0C, I =/r7A J=AJI , Z' 3P S I +OI/ CLT ( IT;)

-10 3-1.4.11 Listing of Subroutine CHF--Y-IIP IIIXII IXI- __

,T L ,)~ITI\" CALCTJ ,Ar, F r-, CTITICAL riAT FltX hNr CTTICFL; Fr j)4C TJF l",J.)LF / ,C'" u"I '4 TO) F

-105-GAP=XGTOT*I 00Q0000./( Y A-A/144. )(3';;.3Ak=GA / O(n r).00CIE CoPPELATIOJPC4 I T=3204.PE F)=XP/PCQ I TaC= (t-o E) / ( ( I ·356c'G&apo

xLIST OF TABLESTab le Page1 Card Image of WOSUB-I Input for Sample Problem 392 Map of Input Data 453 Input Display L6Printout of R

-106-·;,) TC i -n7' I (,V X.GT. P) ,n TO AnC-'ITJL=l0 0n0.o(1.63-(0-e70>r,-qt) -(4.71 AX) Ni)!,C TO 9n:: C-iTJL=1 0rn00nl.*

Go TO ?2;" -107-21' CDITi=CHFR(J)2? CN')iT Ut UEIF(OATA(6) .LE..) CITR=C4TF3/31AQ9.731vt;ITC (6. 23) CPIT923 FOR'AT(//Inx.q4eNFTT

-108-1.4 .12 Listing of Subroutine CHEN-'"'

-109-ClJOjtu TI E rFjCO).,4O/ /tLJh!T/DZ,T,v7FO O.AN,AN1 .7E,: -L,FVCPa.7uP R,P-CGrC',UN /PR P/P9ROO\tAPH\AP,HATAL4MTLAAMgDHSATGAMMA,GAr4! A VRnV

-1i0 -T C=r.] , TrFC=Ac(n r.,_5-(r. n1-.' T~:,i) )P .'F s = . ! . a -:, r P ( ( 1 r . / ( T ';,,, + 4 7 3 . ) ) :1 . O -51-741 .9

-111-2? IF (HTCO(J r ) .T.).)KF tr,(J, Pi) =KFL4G(JqMrD) +1HTCA (.J ,( -IPOD) =HTC VG .J, 4OD) +HTCn (J I )TraVG (J.NIRO) =T=WAVG (J,.'P ) +T

-112-1.4 .1 Listing of Subroutines FUEL and FVS12

-113-SlIUPROUT T 'tF FV 1 2 (FK. X V)IFr() ln02!),20O 1 ' c=-1.GO TO 30?n C= .3' I'C( K- l )40, 50 V=A RS(F-. 1 ( X .- 3

-1 4 -O 2 =Oad~ 3 .-I ~)' ( 3r. a4 ,.''' )2 i~ ~: ~- i i Q T (,t q) D aT. T 1 . C,r: F TU, ti T ! ( i ( n E 1 ?. ) )F

-115-1.4.14 Listing of Subroutine EDIT

CHAPTER 1USER'S MANUAL OF THE COMPUTER CODEWOS UB -I1_____11__1_111111_1_11111.____·1__1_· _

SI; -DIT FD T Tr..>c pi'T ')T t'T T IFLnIi AL I CmlC V '. A r) r'! j g ' ) / \ 7 . '? T. , 7 , r i&apo

-117-STPI;,~& () =nOnI ro? TI',(])P! STQ IG(p) =' C:,FST 3 \Ir-(3)STq I fN) ( ;Q) =FCIONF*SfT I t ( )ST T ',qP ( ' ) =c)O

-118-T -( .1' i T o'i]aT_, m = , ( .T .I,]'!-(-. -2', :-,T:T, ,'b)

-119-1.4.15 Listings of Subroutines WATER and HYDP

' -120-S.JuPU TIE W'TE (P): i-- . i.TR. Dh..-.'TIFq FITTTNr, FUNCTTOJS_~ C~O:'AO. /UT/f37,r.Dv O ,," N1. ,.7E,?'L

-121-SiJ3pOOiJT INF HYOP (EY,' )T H ,UG)' ::--' :' EV4tUT$ YDRULIC PPOP:-RTIECL'GICAL i,. TRV. IVIFPEv\. I TA, ILnC

1.4.16 Listings of Functions TIIriEF FFACT,T'PFi and TFLM"- I-L II-~--- I --- XCIYC~--~-I1"P~--911~_ Ip-ll___·__··-----·(-m·

-12 3-Fl .'.!JT I O, T T '. - F ( * AX T )C ':- lt.. f-IT\TNr, )- I"T; .Nl T I [EGOr TO ( I ,. 3) , T141=1K= IK = K+IT :1R L FF(T

-12 4-I CI ,.qFg~i=~,:,n~~~ -, ., T ( !-J Y -:. 'h~ ..,t. )r .-F-4Ci'F,CT~ TI, C,,-,~F r 1 T],' Ph/q ;"L_ rL.! riJLTItI.[,

-125-1.4.17 Listings of Subroutines INIT, INVERTand TRANSI

-c -1.0 IntroductionWOSUB-I is a subchannel commuter code for the steady-stateand transient thermal-hydraulic analysis of BWR fuel rodbundles. T

-126-I T I t ! I T-"- T T -t. "' : CO'\- 1.' TR C: ", .-'/ r.q :t ./ , 9 tq 2n ,); ( ( 5 V) 3; i I Ni L F T ( G i?.

-127-Rt!QR Olt T \'' - " ,iV;PT (! ' , T i V,K.' IP )C : : VV r r~ T T' * C~"!. h CT[ r !)''~&a

q :fi;4L~JITT T ' I ^>J (XCL! ,' i)~--- .r'LUATF , i r )AI _ CO CC) FFTC T TSC :t- 'I u 17J / / 7 .4 7 -7,.C , t ? FVC, 7 tR

-129-,- y( I (Ntr~ )*J2 +K) =V0+CKHO: 10J,:t (I ) =TJH(;. T-) 1 31-~ IF ( I-'.~C-~,) 17.1 i , 1Q17 I (I-1) 20,20,211': I 1- t A4 ( L

-130-,.( T ) ='?., T i Is = ( ' (L ) - I 1 ::DELF ) /FrI_= ( r '- (L) - ( T i ) -r) F) / lRF I (T.'T. ) GO Tn 2nT-.

-131-1.5 Some Additional Remarks on the Volumetric CrossflowCoefficient Matrix and the Recirculation Loop Concept1.5.1 Subroutine MATSETThe subrouti

-132-NCHA =number of subchannels in thebundleIOCIR = number of recirculation loocsin the bundleMATDIT'.is then the nmber of boundariesfor

-133-The variable NTW (I,J.) is the gap widthbetween channels I and J (1 I, J < NCHAN).Simultaneously, the DO loop 6005 sets up the

-1 34-1.,4 --?-$ byanne -Z a e wh - t- ax z' -' -- <r~~~r~~~,,,,;,3 V~~~~~I --.6Note: 3 Recrculation loops around A., , 8 Subchanne lsIAT

-135-__ _ _ I s I _ _ _ ._.. ~ -- Pein7 P:t7 IFig. 17: 10-Subchannels arraysymmetryNote: 3 RecirculationB, C, E)10 Subchannels(with axisofl

-3-specifically to the code's capability to analyze BWR fuelpin bundles where it is expected to give more realisticanswers than the conventional

Note that for the square geometry arrange-ment, a center-type subchannel has 4 neighbors. Inthe subroutine MATSET there is a point in the

-137-domain centered around the principaldiagonal.In analyzing the matrix more closely, one deducesthat the Gaussian Elimination method is the onl

-138-1.5.3 Recirculation Loop Concept1.5.3.- Se -Up of the ProblemOne of the original features in WOSUB is theuse of the recirculation loop

-139-As shown in volume 1 [ 3 ] these two sets of condi-tions can be written in matrix form as:MR = QTo illustrate the nature of the co

-140Table 7 Matrix I II1 2 3 4 1 5 .6 .7 !9 o1 1 .0 0 0 0 0 i o2 -1 O 1 i O O ; O O i O C3~~~ ~~~~ 0 0 ; j I ~~~~~ ~ ~~~~~, .0 3 0 0 o -1 0 1 1 0 0 0

Note that:a. The array qiJ is initialized to zero beforethe calculation.b. The array qi is initialized and filled with thepreviously discusse

-142-1.5.3.2 Subroutine SOLMAT: Solution for theVolumetric CrossflowsOnce the matrix TMAS (I,J) has been in-verted, subroutine SOLMAT solves

-143-1. 6 Limitations and Execution TimeThis section is devoted to an analysis of the main limitationsof WOSUB-1 with respect to the calculation

-144-This might result, for large cases of matrices, in longcomputer- . calculation time, and therefore costly runs.One should note though, this

-145-CHAPTER 2USER'S MANUAL OF THE COMPUTER CODEWOSUB-II

-4-1.1 Input Data DescripticnAll the data of a problem are read in as a single floatingpoint vector; conversion to integer variables is performe

-146-2.1 IntroductionThe code WOSUB-II extends the calculational capabilitiesof the WOSUB-I code to a maximum of 45 subchannels, 36 typesof rods,

2.2 Implementation of ExtensionsThe WOSUB-II code is defined by the following series ofmodifications which have been brought to the original version

-145-2.2.1 List cf P'odified ArraysThe modified arrays are given in detail in the fcllow-ing tables for the purpose of easy reference in c

-149-MAINIMODIFIED ARRAYSRERGH (45)NAROD (45), KROD (45)HF (36)VINLET (45), GINLET (45)HF (100,-)CHANN (45), A (45)FFLC (-,45)HYD (45), HYPER (45

-150-DISPLAYODIFID ARRAYSiMODIFIED ARRAYSJOIN (45), NJOIN(45,-),MCIRC(28SL (-, 45)RODN (36)RERGH (45)NAROD (45),KROD (45)VINLET (45),GINLET (45

-151-IMODIFIED ARRAYSJOIN (45),SL (45)NJOIN(45,-),MCIRC(28)RODN (36)RERGH (45)NAROD (45),KROD(-,45),VINLET(45)GINLET (45)AF (-,36)HF (100, -)CHA

-152-1MODIFIED ARRAYSNAROD(45),KROD(-,45),VINLET(45)GINLET (45)HF (-,36)HF (100,-)CHANN (45), A(45),FFLC(-,45)HYD(45),HPER(-,45)JOIN(45),NJOIN(45,-

-15 3-TMAT (72,-),N(72,-),VIV2 (72)TMAT (-,72),V (-, 7 2)c.INVERT0090312(72),453xxxC.I.~_ _ _ I

-15 4-IM- -T Tr-nA7R A V NAROD (45) ,KROD(- ,5),VINLT(45)GINLET (45)HF (-,36)HF (00,-)CHANN (45), A(45), FFLC (-,45)HYD (45), HPER (-,45)JOIN

-155-TRANSMODIFIED ARRAYSNAROD (45),KROD(-,45),VINLETGINLET (45)(455HF (-, 36)HF (100,-)CHANN (45),A (45), FFLC (-,45)HYD (45), HPER (-,45)PEAK (36)A

-5-Group 10 Peaking factorsGroup 11 Flux shapeGroup 12 Blockage specificationsGroup 13 Transient specificationsGroup 14 Perturbation timetablesGro

i --156-MODIFIED ARRAYSAJT(45) ,AJVAP945),FIHLIN (45),QT(45), QVT (45), QLT (45)EXCHM (45),TDH(45), GFAJN(45), AJHLIN (45),HDIV(45)L( 45),VFI(

-157-IMODIFIED ARRAYSm ., - ,,,,,SLIP(45),VR7 (45,4)JOIN(45),NJOIN(45,-),MCIRC(28)SL (-,45)RODN (36)NTW (45,-)TMATI (72,-), SRMIX (72,-)SALDD (

-i, S -MIXIN (cont'd.)MODIFTED ARRAYSNAROD (45), KROD (-,45)HF (-,36)HF (100,-)VINLET (45), GINLET (45)CHANN(45), A(45),FFLC(-,HYD(45), HPER(-,

-159-'IMODIFIED ARRAYSAJT(45) ,AJVAP(45) ,FIHLIN(45)QLT(45)QT(45), QUT(45),EXCHM(45) ,HLIN(45) ,HDIV(45)TDH(45),GFL(45),VFI(45)AJN(45),AJVN(45)AJ

-160-(cont'd.)MODIFIED1ARRAYSJOIN(45.),NJOIN(45,-),MCIRC(28)SL(-,45)RODN (36)NTW (45,-)TMATI (72,-),SRMIX(72,-)SALDD (72,-)TMATI (-,72), SRMIX (

-161-1MODIFIED ARRAYSTW3(45,-),W(45)VEC (72), RHS(72)JOIN(45),NJOIN(45,-)MCIRC (28), SL (-,45)RODN (36)NTW (45,-)TMATI (72,-),TMAT(72,-),XLU(72,-T

-162-MODIFIED ARRAYS_iiii i RERGH (45)WATERCARD #1583231x5II . . . . I I4 llxl-lc`"^lsIl"IIIVYI-Y-Y"~"" -··---

-163-IMODIFIED ARRAYSRERGH (45)HYDCARD172834i1x5CC(CCIII'.I__2 1- ----------

-i64-MODIFIED ARRAYSRERGH (45)CONTICARD #16222345xII.-. _ · _ ___ __ ____ __I__ __ ___ ___ __ __ __IIIL

-165-MODIFIED ARRAYSi .IB(45), IRE(45),IFRE(45)IBN (45), IREN(45)RERGH (45) --AJT(45),AJVAP(45),FIHLIN(45)QT (45), QUT(45),QLT(45)EXCHM( 45)TDH(

-6-12.111i jCon t rolF a SRi TFDR m PtI5vInputKMA TE( 2.6)Ida ta Ca rdFc*.PATrLZ .12M npu, daac23 24FIG. 1:Input Deck StructureAfter 1st Control Card

-166-SWEEP (cont'd.)MODIFIED ARRAYS GR (100,-)QTR (-,45),ALF(-,45),HEL(-,45)GR (- ,45)TQUAL (100,-),XQUAL (100,-)TQUAL (-,45), XQUAL (-,45)NAR

C-167-EDIT EDITCARD #176217631764176717681769177017711772177317742IMODIFIED ARRAYSRERGH (45)NAROD (45),KROD(-,45),VINLET(45)GINLET (45)HF (-,36)HF

I-16 8-EDIT EDIT 2 (cont'dlMODIFIED ARRAYSRODN (36)ENI (45)CARD #177517763 .452x1xI iI .I 4J,3(it._- _ ___ _---------

-169-2.3 Input DataThe increase in the array sizes as well as the set up ofnew arrays have resulted in changes in the input, although theinpu

- 70 -additional explanations are given vy placing numbers intcthe column "'Rearks,' wn"hi ch refers te user to Section 1.1

-171-.4i...o ioE-4 (4-14-C-w '"o Q~ra8 r4ao ~ ~ OU v'm 04-i (:0 W s .y8 0 5 cZ, o i A o .p , 4- 4-4 8 .I 0 0 4'-'-u-

-172-azNH o tCl) * o Q) boo t IICD f-.lbO-H I a 4q- a'4.4 o M-ior -- oF4*nh w144-J8( tha) - ColaCI II10 w ( V 2 a)H-CD-14- ru 40) 0

-173-laa4J4 J4-I (~~~~~~~~ ~ ~ ~ 4- 0-w .S <, .S H .S 0 0t, , 4 W 4 nv <w r .w@'x I E S 4.J g ft*H: OHl eC- o-*IIrr ·rA zwr-4 H r4f-IIi-I

-174a lo o ( a g v-I ° iJ C/) w i g 0 X E M tco U; cn UM 4 UJ U) Cn ) 08~ 4 0) c mcoa, -4..*CD f00aS) 00 0 . *CDIraI_·Ii3jIiIII.IIIiII

i f -175-It~~~~~~~~~~~-.,4',~o .orX,~~~~~~~~~~~~I- 4'1 H i M ! i~~ a a t. toc.- O( r l~ ~4~ .. S~i 4 q - -·~ ~~~~·rq1 ~ ·f~ r( ~ om ~4 kC

-7-Therefore, the two remaining zones on the second inputdata card will be ignored by the code.2. If a second set of input data contains

-176-LI'4~4-J 4- LrO "0 0 0o 0 0 _-l*Hi -*H *-1 04-III--r--q..- r-r- %-0N C C H I II C) L4H CN ((D (D CD(0 (0 (.1-I r-%_ CT) _CD D CD

-177-0 0H o HoZO Z O )U (4R 4 .4_ _-v H4-,.N .--bo bO~U) 0 tv 4 -4 1 H, ( eH 4 O 2Q ) O Z. E+ o n )C '- n~x 'N*E -H HIr-8 .tL- I

Is~~~.4-)~~~~~~4 ~ ~ '.4# rX~ g 3e T B H-4 t-H >I; ' *r-4 40~~ r-.44-4~~~P1 0 t4-) Ls A40 (4 44-- 4 024- C/o'I ~H c

C-179-0p .Oto440 00 I aC) > a-4HH *,1 > H 0}c CO 4 -H0 00C4 S C awH0=:Ezcoc4tD.C(IrltO4-H- r-,--04o~.~,-IO 0(O l OC0 -<004-4Oit >. -Fr-X

01 O %b .-I Z 04 1-4 -4.74 iO F- 4-J --%la :$ L1toxaC 30 $4-cc3k 4i -% H:3 S 4 a) 1-4 p -O ) I-TI 4 4 0 v a o sm D CD-r 1- cer3 rz -4)

-181-W tO-Z 4 OO - ;H 4- 4 OE-,4-H [1 4- IIP a +- P -c a -40c 3. CC.) 4 3II0,140vz--4.0C.CDC"V.1WCMIC(tHlH4-H .-q*r ' *H H 00 4J

-182-00 0)0H I- 0Pe EH O : L)0 b00 OUEzi 4-0H O0 U0ro0-O Z0C) 40 C.,- aUO uC.)a4-0 rP O* 04*r OH HU 0_ 1U-H CNN c E- 4 W0,-E0C)zX::DS tH- H r

pC.04oCD0I0 II- 03- 0rd 4Jar)04Imto4ioW >E t 00 0 r-4O a .a HC C.4-iE- > 3 C96 o04 Oo0 -0o- 4 .I-,4-A.01rU003 .03).H~~i(Iz0E-'-~-InI3IE

-184 -2.4 Output DescriptionThe organization of the output f iOSUB-II is essentiallythe same as that used by WOSUB-I. Therefore, the user shoul

-185-2.6 Listing of WOSUB-II

lwow W

1iCon trotCaxdIFIG. 2: Input Deck StructureAfter 2nd Control Card1d 2i-8-- T (Ew Plo. 6ONPIT (2.-ACOT-HA CElI6)4 a Zs cN- ---------- -- -- -

-186-2.6.1 L sting of M'AT

M. - " el M eO'1COOO00000000o000000I.l4 L b. -¢ Z 3.3 Z " ,r"-~- U'U -0a 3 mn a CC Z t- ::> E ~ u - -*aP-*C_ :

O0 O Q O C-1.88-- . 0 o e S j- a'- r -" Ln T,. "000000::oCooOo'000000z - 0 0 0 0 0a c 00000000o 0O0000000 O 00 Co 0 0

-189-v- NC" C0000000 '0aC,0000.0i!0cr 000000000000OCNI _ 0000 C 0000000000co0o CCCC :t : · Ir2 2 -OSn%w'w0,5-'C iF4'®IIC.)

0 00 0 0 0 0 0 0 0 0 0 0 0 0 0-190-o}~~~~~~1 tVS%O L!(IJ ENC CC, 0C .CC ;*\ 01cacV)LIIaI-'Cc VLA FC i ' 4- 0 0 '

-191-r r-O CCOo00N,I IM U -to r Wrv r r r , IrCO0COO0C 0o00C C C C) oooIIIC0 0!C N t0 0c COocCC0 C C 00 C 00000 0 00000o o o C o( r a: Ma C 0000 C

0 0 9 o 9 oI§~~~ .~-192-O 0 @ o 0 'o. 0-Un 0i 0 0-:v * * C * %IM0 0 0 0 00 0 a 0 0C 0 o o0 o C 0 0 C0 0 0 0 0I* _CI'

I~~ .r -0~*t. i~-0 0 0o 0 0o 0 0 0o 0 0 0'd-cl3 O O O660O O O9~~~~~~~~~~~~~~~~1 \. r..;C42 0r.4 ·cO r~Or0 .J ,.; O-* f ." ~~~ ~M

2.6.2 Listing of DISPL.-_111----  .,1---------11---z 1-----·--- ----- c-·Y---------- --pl--PI--------U--·-L·_l-·ll··· L

O0%000 0co00002 -CE Z42U34V2 00..LU sS| *. 3 C,· ;: E-.C M - :N rr% 3 < ua:4 r SC) ~ o= 44.4 M-t -s? l "0 ZI- M! -t) =4 *

j-9-l 2I1o 1o1 V2120 0La<, Contro( CmatLouAs IrA Pvk(CAAFIG. 3: Input Deck StructureAfter Last Control CardI1A21 22 t 2.tz olo o_ _ _1 __

196-0 C '" ..000E-.C .Omn~~~- 4IIC4_ eC- O C-tl,C· CC _CIa: F4 (*c _ SC-OCO OC CC O00O_O COCC C0 C0 CC,F--OU..-t'-4C;0C.-4-e

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0I-197-C- O P" C" o o o o0CC G OO OCCt C Ck C a: * . .400.3 i O 0_C/- es0C. O*r0 . s-0f,-

-198-2.6.3 Listing of EDIT

0 0 0 00 0 0 0 0 0.0 0 a 0 ·i!f!II-199.i.ii IIiiit!iiIsi1I.4 o-0, e, -D. _ ._ 40: * ,e ...-O '* ,2 ..I. I( AC-01 U, M *t *: *C- -IO C.C c

-200-C' C o -C o _:'LiCr 00TC , C Cc C 0 CCOOtf- - -C C0o0=COU. ,u f- ,V- 0 a-C 0 O CC4 CN -. C4. C (" O : .C C O C

... _-- -_- ...0 0 0 0 0 0 0 0 0 0 0 0 0 Q 0 0 0 0C . Co W c c . Z( iii: ii io000cco000o000000-4U-wI,40:3r b 3U-AIi npL -

-202-2.6.4 Listings of Functions ROFT UN andSubroutine GEOIRY------

-203-U roo oo oO0o co C00o o000C 0 C-.-Q.bCLIC.4.u._ .. +L~cL.C -·C;Z ,' -- C-CLL ' .. ... - CQ 9 0 @ a 0- 0 0 @ 0 G It

@ @ * @@, §@, * 0 0 0 0 0 0-20--2 0 4-4 u) Io 0 c r - MT ZN N N N N ft n FI ^ m r? c,o0 0 0 000000o oo o0000000000000COOC o C C C

0 O, 0 0 0 0 0 0 0 O -205-00 00C C000000 000C CC a a O C O0COCC coo-4UCC000000000000OOOOOOOOOOCO0000000000000OOOCCO OOOO0a 0% o- N M U On r- No %o

-10-3. Large flexibility of inputting a set of input dataof any size since the control card does not operate with afixed field value.1.

-2 06-2.6.5 Listing of Subroutine STEADY

**< ****0 3 0 * 4-207-qf3 3 Q 3 ,, , 0 o h a uCo000 C 000 C 0 000000COOCOCC000C0000 o o C o C o o o o o 0 0 0 oaC .> CC .

'~ 0 0 @ @ *-0 0 * 0 * * * 0 0 0 -20 8-r- M C' o -N r. tr WD , & -Ln r u nm m LnO O O O ; C C C C,0000000000000 00 00000OCCC00000CO

-209-2.6.6 Listing of Subroutine TRANS

L -210-8 .j-210-0000C00c NCD CC 000 C C CO o3 . I.00 C000000 % r00O 00 Ca 00O C0 0O0 O0 o oo o0 o0 000C >LAL . L. Z au (>- Co CI Ja 4

00,ID<8 tR-M- * ' *. .f * I -*0f---2l -*-211-N 4 Ln Z rl- Cx 0a C O C O a oa C o C 0 0 a a=0Q0 0 a a CN-C 0 00 00O0 00

Q 0 0 0 0* 3 0 0@@ @ * -0 0 -2 12-_ ` O 0 5 30.0000000000000CO000000000000000CO C:CI-!-.-JI;-J(._IC.CCCCJ~~~~~f'~ ~ ~ ~ ~ ~ ~ +-,m-,Y~~b_I

-213-2.6.7 Listings of Subroutines VMEX and MIXIN

0 .) * * * * * @ @ 0 0 -21 4-\C o rr IU r) X U) £C > ac 6 C _~nLu 3 eS t 3 3 3 3 3 e ;t n LE Vnm Fl rn c1 rn .n rn ?

o0 a a 0 0 0 0 0 @ 0 a a @ 0 aG-215-CCOO C0CCOCCCOCCo C C= a o= C CX d, LU_.. > rL 4* .J- CC* ,L.= e,-.Z _4 P) O-_ C- S Z -- U -<

-11-r-Ir- c',z aOI--H ( r1H,: 4P .II IIo r1-0o .zC) 'J- )c hao H rHC) C ) 4 U)C)rO 0~* ~ * ' (a C) M l *H4Q0 0 c +-4 -H 4-s J Q a

,.. --,c ,,L 7ror, *'' '-' ' 7 ; .' C .v -.-' ;- __Z 0 @ a , @0000 B **9 , o o @ @0-9-2 1 6-I C _ '- N

m IN --W~~~~~~r --- ---irris 0 0 a a a 0 0 0 0 0 0 0 0-217---- - M NNMMMMM0000C0000000000coocc o4000COQ'._e+ I- ,.,.%, ~-. l--+ - -+

-218-2.6.8 Lsting of Subroutine SPLI7TD~-cl~~~-`---slYI-'- ----- -- -·----I^C·--)-·lll-I -- ··l·----·lgllP--·-X--PII--CI-·

-rI M I f -i -[rA ,@ o,.*v* ** 0 0*@-219-P"0000a-LAZ L L .LA.- lt L 4h'.4- _J LP tL 7_ -aaU1·SUN r t3 * W-I'3¢ _ -In U&qu

I~r·;~t~, ~ i~~i-j-:10 a a - O 0 O-i 0 - -o I O·--·x-- Y"-: '-·-.Ur-- -- ··- I-1- ----*?;I r. -··^--·)r"·x"2-r

-221-cll Cc C- c I D P-O- . --. -O-" 0-Ln n nX . Ln- - -l - - -O 0 0 C 0 COO C0 00 COCCO0OCC 0 OO O C¢ air I_ _CO0 0CC0 U cJ nt

-222-2.6.9 Listings of Subroutines ATER and HYDP-I __·-_I_ C--- -

iM·* cvuc c-qc;r· ·'"`"rP'VL*·4;T.I :2 ;;*$-bP;;.ia `i·Si6'-kZ 74i'L·F 44Liu:e-"" j

- '; sC ~l3 F t;,-~- 'i: -"-"'- --s-22 4-0, o- -c LU tr, CC -o0 oo CC o0-- 0 000CLL- .LL.C >CJ9 I--.m* -- -. --

-225-2.6.10 Listings of Functions TIMEF, FFACT,TPFM and TFLM

-12 -Q)N.-O CU) . .C r''C C4,I 0)II II !I IIO: -oH ~'"'C-H-1H 4O> ) CIC) r-4-)C.CO T rC CG) - Q) III Q IC O

-22 6-r t LC . ¢% ", ° ,NMn U) L L \ tn Ln D D OC C C O 0 CU00O0000000C C C C C C C G OO:C C C C C O C O0 CICcr--fD ,-C

TZ~~~~~~~~~~~w-227-- --22 7-0% o- N .t Ln I ne e m t t It m.o Co C Coo o000 C 0000 000 CCC 0 orLr:.Cb...o '- -. e-,L L..I L+_ _c: _* I i

-22 8-2 .. 11 Listing of Subroutine CONiTI"---------P-------IIIIIIIPIII'UI-rX 111--·--···_··11-···IIII_IX

-229-4 n e in \ c: ;r o _ N 0 -e;n 4o ; Ls A Qa o _ L ^ O ur 9 > s-,.O .& 'D~ ^0 ~J t (M c' D -' X ( M . 4 . ., .

' '.~ .. Y " "''-'...' .. ' '~ :~ '~ -~ ~' " ' ., '~.. _ _ A

-231-0 00 t 0S00 . 0 0 ar,- M C4 X -Z rl w O O o o O o C o0 oIN,_!I I- :- I,_ _ P. P,_i-- -4 - 4 - -I, .- - I, fl >C 0 0000

-2 32-2.6.12 Listing of Subroutine SWEEPI --

90 -* * 0,0O0O0@0 00 0 00-2 0 0-2 33-Ln Dor r-0o cc o0po40- -4 ..-- I-,000o 0o oo aoLOcxOC-,C.-JO9in4?0WAUnU.t'-4.,--* 4 4 -* . -( c*.i Lt

-2 .3 4... ' ~ ~ ~ ,~-*.~'" ~ ;; ~' :L~ ~ . ~' ~ -,~~rt,~k'~ ~ '~*-,:~:~.~..~,.~.' .'~ ?&apo

I(, @0,0 @ 0 * 0 00 0 0 0 0-235-n 0D 1 0 OM o- I~,-N M, LtD I- uO 0 -N M.tn %C!'- C- , ?, C D o N I 4 I 4 4 0 UM Ln I Un UM Ln tr U

-13-bO bOd i LUv v84)J E4-H .r-o -HCDW Ev 4-4 2 44CN cr_S _ _e4Q)o ,- .oI -Q) 0 -x) a )O W) -0o4Z v-0bJ.i8C eLD ( (D OD ri dC- C

-~ ~ ~ ~ ~ ~ ~~~~·j -·C·);TI .r," 00@ @00 o Q @ 0 0-236Ljr C I - C . Q -U ", Lr ! %- C: , C'"ll-4 - -" """

:-a '- tS· :u·~~~;·vu,.t'.rr~Lnrrrcau_~- '--I-~-Y . _s-- > -r * - --···w, ·-L·4;rs.rirt·ft gcidrir ·dir.4-zc L:S y

- '-- ;C--'· s -4S -. rs - 4·-Y·-_ - -+sg- . -. M ... ·-- r·:-:-·' ._ ...- .-. ·-_,.f' ir ·- -. . --'-·-

0 ~ 0,, 0 0@ 0 9 0 0 0 O-239 -239-00000 rl r en e f C -N -.- n -o r- , c- C- C CN CM . _" _* _e _ "" A0 -4 ; # N N N N..

-240-2 .6 .13 Listing of Subroutine SOLNrAT- -------- DI ·-- Yl -·l s1·- ·^~)~- -~W ILI~- ~ ---- II-I - ILL -I--I(q - _ · ~ ll_ l^ P

0 0 0 0 0 0 0 0 0 0 0 0 0 -241-I C -- , t;C -F, C . , . r _- - -C . ... r *. I C *Z.? ' C -_ .,- 7C CYt n .- '- .'*

-24 2-2.6.1Ci Listing of Subroutine MATSET

~~~- -~~~~~~~~~~~~~~~~~~~- .* , , 0 ·0 0-243-a C C o CD000000000000000000000000o oc)oo000D c0O O _ _.C00 0000000000000000:> pM pr 4 100000000

:) 0 O Ql $ @ G : O .' ~n O O -Q -244-'C -a O -' 0" C' C: LC 0 L a;O 0 0 -: t- t-C U N ' .- 0 r' C 0 C 0 C 0

a ' 0@ ···0 0S 0 @ a-245-C o -U· o '- r -0 00000a O C C C0 00 00 0 00O000 0 0 00C 0, ° O-N o Lco Cm r X 0 0 _;X a D .3 X 3 9

4Z sio o obD t. tb3 Uc C rb t3- C b% rr oz Cr C a; H O c60O C C cl-i C *H 4O O O O O 'E £ wr t o 2 C r: TJ C l -C r_, 0 ,0 c0 0 ·r , 0 *C ::a Q)

-246-2.6.1: Listing of Subroutine iN\RT- ------ -- I

* @00 @0°* 000 Oo -247-L. I Cc CC0 00000 0000 I OO OOO0 D 0 0 CD O - NO OC OC O"_ , - ,_ 0..0 QOOQ -.!__. G 0 0 'G--C0T..CcP a:y -

I··P

.AwtREFERENCES[1] L. J. Guillebaud: Vapor Diffusion in BWR Fuel RodBundles, N. E. and M. S. Thesis, Department ofNuclear Engineering, M.I.T., Januar

-15-.4-'.-8 's rQ o6o::3 o .ooC) .r W3o t @ oi-4o (r C) Q..#H:s -r-4 U X; H aCNf-l0c" J a)8 C o 0- O o o---4 *,-4 t~ 0C

-16-rj- 0O ' Ur > 9E U JLO0... r-4 IZ, ~ TaL C" O 4 a) LI 4O *3 0 0 0O)- C 0 ' itV1 1 1 1I-, l- .- ,-L3 w w

ao oH-)r--0 0H ICH za CW9 HZ O-4--C0)(nJ r-4'a') bO4~-412z O0co)030)HD DN N0I!I-IO1r-CN CN C,Ir rI IH- H- H--: U3 O(O (D COo o c&ap

-WOSUB-A SUBCHANNEL CODE FOR STEADY-STATE ANDTRANSIENT THERMAL-HYDRAULIC ANALYSIS OFBWR FUEL PIN BUNDLESVOLUME IIUSER'S MANUALbyL. Guillebaud, A

-18-.H C * -4-) 44- a) ·O IC 4-i U C0 O O.a o a Q)4.) r U C) ( I C) h 4-:Ez 30n:O;3 ; H :-3(o (0 LoCN CN CZri I Ir rl CNC)

-1.9-V4-004H: 3:0 0~O o a> OOHHk e O·X -a}X C?) > @H0Cn *H > H W 4J (O =H P: C i a< > 0P P0z:(oC%-WD IO Ho I.03oo0E-

i-20O--2H 00 -PH OiO + .- *0 m o N: ~1 H 0 0X~~~~~~~~~~~~~~~~~~~b 4J O 4 U; >4 -:. 4-J O· ~ ~,--~ ~ ~{~--~ ~ ~::~::~ o ~ I { 0 : ::

-2 1-)o 0 aFe O ) mO O ., .~ I0 -r-q @- IEn , c r-i 4- >H 4-' I ra 0 Pz la-4 ; 4- Q)a 4 c) r_ : a. z a)Ca)0a) r0OaC) 4 3 H P -U. , _E-0A0CHzHC

a-22-0OZ 4- Uo a) H H I)H ~aCL. EH O toW- CC)-t; U tC) aW OUa) .laCE-H 4-4HO 00 HOor OO -C)0Ou 0 O04-i*H 0au E0 Om .HN N No t.oCZ CD (0.- q

CI1r4cnto ~ ~ ~ ~ 3o0I~( . 0Z ,. 0 /0.O( *S .r4a 4 + ,j0O-1 H0 , aa) ·>%CJ3 r CO tbO X0 IC.) C3 3o o ofWWX4O r-I H 9 e0 0to0(13·r 3r a)n C nH HO CK

-2 4-a) HZ C)C) 4D3 a)ch c W)1 04)'- 5) Ce Q):5 a) CePCd .L> S v C, U O4.. )a' aC0 -4 a), aCTC aa ,X F.444CO-d 5 0HH

-25-1.1.3.2 Recommendations for Input Data Selection1. NCHAN = number of subchannel types in thebundle.A subchannel is defined by its geometry

-26-One has:16 cards of GROUP 3 or from DATA 11 to DATA 9016 cards of GROUP 5 or from DATA 401 to DATA 4803 cards of GROUP 6 or from DATA 601 to DAT

:)-27-..) , .. .. .~~~~~~~~~~. rjI.CI(TYPICAL) (TYPICAL)Fig. 4 Nine-Rod Bundle with Uniform PowerWC.DSDistribution

iiA CKNOWLE DGEMENTSThis research was initiated about July 1976 when LouisGuillebaud pursued the assessment and extension of the MATTEO-code,

(-2 8-CCCFIG. 5: 9-Rod Bundle with Non-Uniform Power Distribution(from GE-13049 )4Al -_ -%llk l-^-----C^ ---------·I --·---p--------- ---nrul

-29-iFIG. 6: 8 Subchannel Bundle Case (Symmetric withRespect to the Second Column of Rods)IIll-RXZ6L r a .C

-30-FIG. 7: 10 Subchannel Bundle Case (Symmetric withRespect to the Principal Diagonal)3ss -----·--·.--·IIIC-···-----·111- -·-- IIIIICI---··--r·

-31-3. ICON = Key for selecting input system of units. SeeSection 12 for a complete discussion of the procedures used by thecode in both cases of uni

-32-Wv;ck BFig. 8: Si.mple subchannel layout to describeinput data #11 through #20·1(II______________ -·illC··-----··lll-P·ll I___ :.Lt A.

-33-am: Sapie ubchannei numering scheme usedfor the explanation of the concept ofrecirculation loopLoop around rod 1: 01050602Loop 2 around rod 2: 020

-34-Therefore, the input data for subchannel 2 are:JOIN (2) = 2NAROD (2) = 2NJOIN (2,1) = 1 NJOIN (2,2) = 3 UIJCIJN (2,3) = 0NJTN (2,4) = 0KROD (1,2)

-35-7, 8 NOCIR and MCIRC(K) Recirculation loop specificationsWOSUB-I used the concept of the recirculation loop infuel rod bundles. This concept

I-36-9, 10 CHAJN (I) and RODS (L)For the cse of 8 subchannel-s as shao-,,in Fig. 6, theseinput data read if' entered by the first option:

-37-1.1.3.4 Sample InputFig. 10 shows the sample case layout. Thedimensions are given in inches. Other data of importanceare:PressureTotal bundle p

iiiall the pieces of information and to implement them into thecode. His first achievement is the addition of a novel fuel pintemperat

-38-Fig. 10: Sample Case Layoutc ic.nc.e--- tc. 14

I * *** I ,' M 3Lot _ OLF1 104 .11.I .25.1.2 (2.30311.1I .40A.'f)3.604.2 .704.71514.4..aO4.9qS.10 3 (' r4 U 10I ..7;'4110 I-,i-i

1.1.1200on0.1.51 000n.9. Table 1 (Continued)1.1.1.qr .41 .-.i -; - :·~ ] 1.1.1 , 1. I 2i I ? ) * C I2I 0 ;i1.is21420 1?0 I2n;l)? 0 I?? o O "

1.2 Output DescriptionThe printout of the input data and the results areorganized by a hierarchy of three levels.1. After having read

1.2.1 List of Some Important Notations Used in theOutput of the CodeGTOT Total mass flow into the bundleFMOUT Total mass flow out of the bun

-43-PDROP Pressure drop across one nodeITER Number of iterations to achieve convergenceAVG HTC Average heat transfer coefficientAVG TW Average clad

-4 4-1.2.2 Sample Output~ ~c19"1~1c~s~-~~~"~~~1~1~~~"~~1"~"11111 l1sli

O C0 00o 00 0o 0o 00 0 0 0 0 0 0 0 0 0 0 0 ¢N 0 N0 0 0 6 01' 0 0 0; 0 0 0 * 0 0 0o 0 0 0 0 o 0 0 0 0 0 0

-46-0000ho4xCaCV0Zo000a'0000 0 0000aaa100A000oo0oao0 N E N _,r o 100-- -- e , U n o:.00000000 i,C;SU1~~~9~EOC'I4tc00 A~C 0 (92Eo2hiii

-47--0la0Isco..-I00014U.Pobem:nn o.-4o '0IiIiiiIIiC'.4 .0r · I. .r4 N-,,'"CU,I. =o o m."t.C'O o _II.* Ce "U, 2 ZC:

ABSTRA.CTThe WOSUB-codes are spin-offs and extensions of the MATTEO-code 2 ]. The series of reports describe WOSUB-Iand WOSUB-II in their

4 8-L)04I30Z (N=000o (o00(N 040 "i, > 0oir30003a N 0C Otwar.0- -ooo mOml::* 0CaC :*FCa 0ao00' S 0 0 0a -0r· ~ UC C9 _~-3

-49- .co 0E vC 9 0I 'a a JIl a' ifl 0 II. II * a' ' a 0 in in in p aM O 0 ak 'C 0 U- a' %a 3 p. a ~o in

CNOE- NE-C):. rO OCiLr~,OCSC)OP CC 1 -C)l, C1(7 Ia: - CHC)_ C)* C-'O zF- _ 'a_ C).-50-C 0 C 0 C) 00 0 0t 3 = In ' I I I L _

-51-C '4N N 'N N , N N N N C (N c N N N N N m 'x --q.. 2 , n c c'H, C ) O 0 r~ ,D O 0 I O-000 000 000000O *-4O O

-52-1.3 Organization of the Code1.3.1 Flow ChartThe flow chart of the overall program organizationis given in Fig. 11. MAIN calls eith

-53-VMIX Calculation of vapor mixing in subchannelsMIXIN Vapor mixing modelentrySLIPF Transverse slip modelentryREYSF Turbulent shear stress calcul

GENERALMAIN PROGRAMSTARTREAD AND PRINTDATA VECTORSET OPTIONS ANDI GEOMETRY DATAROFUN;WATER- DENSITYWATER:WATERPROPERT-I ESI SET TRANSIENT DATAIFFINA

ZO JJ Lz -L iU UD LU-- -F >4za0 CZ HO kF 3r- -LLr,I a: .~~~~~~~~~'-4-55-0-JI- a.4CL0Zc xo0<w0:(n0eOI0H00JLUIHLUIO-Jz0wLUJJ-jzHLUV)

TRANSI EN TT= T + TL i ~ iiP VARIESnoyesy e sNEW TIMESTEP4FUEL PIN Lno Fl NAIyesRETURNTO MAINCHiyesy EDIT:SUMMARYPRINTOUTEDIT 2:DETAILED O

CALCULAT IONSUBROUTIN ESCHEMESWEEPITERATIVE(SEE DETAILLOOPNEXT PAGE)D2-PHASE FLOWMULTIPLIERNEW MIIX rNGCOE FFI CIENTSFig. : Subroutine S&apos

VTABLE OF CONTENTSPageChapter 1. WOSUB-I1.0 Introduction 21.1 Input Data Description 41.1.1 Ordering scheme for the input card deck 41.1

iTERATIVEOF SUBROUTINEITMAX = 10REL- RELI(I.) first value of relaxation parameterI I]I SUBCHANNELS CAL-CULATIONS (DO LOOP)AVERAGE pitno /

-59-1.3.1.2 List of common blocksThe common blocks of the WOSUB-II codeare listed in Table 6 where it is shown whatsubroutine needs what

-60-EIuQ _ __X_ x~Vi~n~ 1X ~= S D9nr x In -.x Ix xx xI_ _ _ __ -- -_ -x x _ _os___ = - _ x _, --- _ _ _XIt x x > 8eXiXM X X | X 8dSII x xLNIUd 8x

-61-1.3.2 Parameters Specified in the CodeThis section deals with two different kinds ofparameters set internally by the code.1.3.2.1 Parameters Used

-62-The built-in conversion factors in MAIN andEDIT/EDIT 2 are:HCONF = .4304 (For enthalpy)GCONF = .007373 (For mass flow rates)FICONF = .0

-63-REL 1 = 1.REL 2 = .7Moreover the code uses with this option:GRAD = 980. *ABS ( CGRAV)ANI = 1. /ANANM = 1. /(AN 1.)

-6 4 -1.4 Listing of WOSUB-I1.4.1 Listing of MAIN

-65-.:.. w- THc- T'I3T ,",ll PCOC-'S T!F OPTTDIY4S "r i. 5 ''I-'"Ivv CALCJLtATTqn'AtR cAF)FL!! 7GICaL [C

T ="C ' r ' = ,rL ( LI )T L 4Ar V="', -I -,'I '=r) 1 ?) ,7 + P, ! j -- - l ( L C ))Tcr () ) TJ T

I =1." 1 = 1 .A1 -1 . INj;I'(I= rJT (arT(T (L+ ) -1001 )' =-,CO (I = l,T (OaTa L+?) t' 001 I=,_ I ', L .:1 ': _

viPage1.4 Listing of WOSUB-I 641.4.1 Listing of AIN 641.4.2 Listing of subroutine DISPL 701.4.3 Listing of subroutine GEOMRY 74i..4 Li

·. WI~~~~~w w w mw- I·~ IIII---1··~~ LII~YILIUC ~ ~ IL· C~lel)IL I C--6 8-'K U,. ., =7TT . . - "1 t h zq H/Fr (. ) =!F I i:

1~l-. -~i--;3ascP I JI -- _, _ ,,,.t --1 0-69-1 71 FLC ( IL . I) =OAT (LI)1A" 1 L=L+20I°(IPI-R)(O TO 000I;'(nATA(901) .GT..)ANI=!)ATAA(901)7

iT ('r (32' ) :1.L T I 7-jvJ =Ie' T ( ,S0T. t n:7]): -in: I ( ;. T L ( Pi3 q) -: C:F;C TA ( ^ )3 ). 4 TA ( )r )PC =- A T

1.4.2 Listing of Subroutine DISPL---------- ----- ---

---- ----- *------ -·- -- ---- ---- ------- ----·-- -------- --- --------- ·------ --·- ---------- ·------- -- ·- -.. U

~~~~w .~~~-71-,S, ' "~Ji~T T E n I D' '. ::* T :P T O I DL. Y:r -:L KF .KCij 3!,T'CL IC)ON. IF)DIFTt-)- ".J;&apo

I'( ';:);F ) TO -2-- r _ (~. "'' T) 7 )':'T iEi J, TIcTa F T PJT , A S S Tt F) O TOTTF, YF)A T T .P ;&

-73-I OAP CO)"IOUC TA"CF I: C ,c .n.4 .)i,- .r: ! ,C -- -- -I- - ---- I --

-7L-1.4 .3 Listing of Subroutine GECTIRY

-75--;:-*:-;:- t I ) r , IrTTY FITTT\IG : . ICll!T ') THj ' PRESSIi: '= l | '= 1 *7 -f !%7; POUT TI Grn-.N ¥r ..- .- T

viiPage1.5.3 Recirculation loop concept 1381.5.3.1 Set up of the problem 1381.5.3.2 Subroutine SOLMAT: Solution forvolumetric crossflow 1421.6 Lim

- ...: .,,J! ...I;~( tJ ' : *:/ J¥I I T(I;T T-=T .T'¾-r-L (K . Ij =1 ._-I , ( "= I] ~ c 'V ~ | = a -L( r + , I1 ,&apos

-77-1.4.4 Listing of Subroutine MATSET

-78-q-".) tT TT r- T I T7 Tr:-T lTipC -^ A&T?.]'A .-C.LS CFl'W'E-CTI)'~ A-ND LO',it r .'h~, () , . .rT r).9&ap

-79-, no7 ""tI (6.hrn) (TJ(J) j=..*<)' ic (=IJ( I )'~.=IJ(I)0o I = 1 ."AT IMJ) =JOITN.)-0 n V009 J=1 ,JJIrF(JOf~,J(J)O.

~~~~aan ~ ·-. --. --I- nssfl W WI··~r*u·r ~ur···· lr~-flfl --- -·IU·;~·I~.Y~r~C·nHIl~ l~-80 ---------- --I--- ----1-- _-I... _ .._I -._ ____ ___

-81-1.4.5 Listing of Subroutine STEADY111-...__ ... ____

-82-r .-:-:.-.- -T;LLJY $TATF 313TI". %CT5 I'.ITAt. C3I.lDTiOS SL*,^, " DH 4 /L J" .?/ T ,¢TV7'pi h, ,7; ? V P 7^, 7! R D.. .

-83-1.4.6 Listing of Subroutine TRANS--- 1 1-11---

-84-S;j 'i O;T'T T.E T ^J'S~-:- -::-:I: T'.SIiilTjTrIE. I -TST iL;LET CONflITiJ)OJS F FCH TImF STCP, --: COn ;rTOi q P>IT

-85-2 C .;T 1, "J 'EIf(I P;L..-'.f) f T 3V I.! =V .1v! T"!=T IF (T IVTNLCVI'i ,vIZERoQTR.),,cA==F-XP f -Sr)QT(VI,/ (f).1 )V