*                                                                             
***********************************************************************         
*BWR/5 core having 764 9x9 bundles, 3323 MWt, UO2 fueled (1/8th is modeled)
*
*No PLFRs are modeled, i.e. all rods are FLFRs (Hench-Gillis corr development did not account for this design feature)
*Lumping approach used, i.e. each bundle is a single macro-rod
*Cruciform control rod channels are lumped in macro-channels
*
*The two water rods actually present in the real bundles are modeled as a single larger hollow tube, through which no
*coolant is assumed to flow. The total core flow entered here is equal to the actual core flow minus 4% (4% is the flow
*percentage that in a real plant flows through water rods)
*
*Core radial bundle PF=1.51
*Axial PF= 1.47
*Local PF expressed by means of max J1 factor (J1=1.198, corresponding to a max local PF of 1.26)
*
*(Hench-Gillis correlation used to compute CPR)
*
*Because of the use of the Hench-Gillis correlation, core total length and heated length are set equal. When
*this corr is used and the two length differ, VIPRE does not predict correctly the bundle power. By setting
*them equal, MCPR is predicted correctly, while core pres drop is underestimated by a factor of about 1.046
*
1,0,0 
UO2 Core  
*
*
***********************
*GROUP GEOM
*
geom,126,126,30,0,7,0                         *geom.1: 126 channels, 30 ax.nodes, 7 channel types
145.98, 0.0,0.00                              *geom.2: total length
*
90, 14.4357,129.3078,102.290                  *geom.5: number of channels, flow area, wetted and heated perimetr (full bundles)
3,5,6,10,11,12,14,15,16,17,22,23,24,25,26,28, *geom.6: channel # 
29,30,31,32,33,39,40,41,42,43,44,45,47,48,49,50, 
51,52,53,54,61,62,63,64,65,66,67,68,69,71,72,73, 
74,75,76,77,78,79,80,87,88,89,90,91,92,93,94,95, 
96,97,98,99,100,101,102,103,104,105,106,111,112,113,114,115, 
116,117,118,119,120,121,122,123,124,125,
*
11,  7.2178, 64.6539, 51.145                 *geom.5: number of channels, flow area, wetted and heated perimetr (half bundles along the core diagonal)
1,4,7,13,18,27,34,46,55,70,81,               *geom.6: channel # 
*
6, 11.6006, 9.4762, 9.4762           *geom.5: num. of chan, flow area, wet. and heat. per. (1st type of cruc.control rod insertion chan) 
8,19,35,56,82,107,                   *geom.6: channel #
*
1,  4.9442, 4.0612, 4.0612           *geom.5: num. of chan, flow area, wet. and heat. per. (2nd type of cruc.control rod insertion chan)
2,                                   *geom.6: channel #
*
4, 6.6778, 5.415, 5.415              *geom.5: num. of chan, flow area, wet. and heat. per. (2nd type of cruc.control rod insertion chan)
9,21,38,60,                          *geom.6: channel # 
*
13, 13.3556, 10.83, 10.83            *geom.5: num. of chan, flow area, wet. and heat. per. (3rd type of cruc.control rod insertion chan)
20,36,37,57,58,59,83,84,85,86,108,109,110,    *geom.6: channel # 
*
1, 32.50,35.1975,35.1975             *geom.5: num. of chan, flow area, wet. and heat. per. (Channel between outermost budle ring and core baffle)
126,                                 *geom.6: channel #
*
41,0,1,                              *geom.7: number of gaps
1,2,8,0.8775,1.00                    *geom.11: gap #, from chan, to chan, gap width, centroid distance
2,8,9,1.1700,1.00 
3,8,19,0.8775,1.00 
4,9,20,1.1700,1.00 
5,19,20,1.1700,1.00 
6,19,35,0.8775,1.00 
7,20,21,1.1700,1.00 
8,20,36,1.1700,1.00 
9,21,37,1.1700,1.00 
10,35,36,1.1700,1.00 
11,35,56,0.8775,1.00 
12,36,37,1.1700,1.00 
13,36,57,1.1700,1.00 
14,37,38,1.1700,1.00 
15,37,58,1.1700,1.00 
16,38,59,1.1700,1.00 
17,56,57,1.1700,1.00 
18,56,82,0.8775,1.00 
19,57,58,1.1700,1.00 
20,57,83,1.1700,1.00 
21,58,59,1.1700,1.00 
22,58,84,1.1700,1.00 
23,59,60,1.1700,1.00 
24,59,85,1.1700,1.00 
25,60,86,1.1700,1.00 
26,82,83,1.1700,1.00 
27,82,107,0.8775,1.00 
28,83,84,1.1700,1.00 
29,83,108,1.1700,1.00 
30,84,85,1.1700,1.00 
31,84,109,1.1700,1.00 
32,85,86,1.1700,1.00 
33,85,110,1.1700,1.00 
34,86,126,2.3400,1.00 
35,107,108,1.1700,1.00 
36,107,126,0.8775,1.00 
37,108,109,1.1700,1.00 
38,108,126,1.1700,1.00 
39,109,110,1.1700,1.00 
40,109,126,1.1700,1.00 
41,110,126,2.3400,1.00 
*
**************************
*GROUP PROP
*
prop,0,1,2,1                    *prop.1: EPRI water properties
*
**************************
*GROUP RODS
*
rods,1,303,1,3,1,0,0,4,0,-1,0   *rods.1: num.of ax profiles, numb.of rods, conduct.model, 3 rod geom types, 1-extra mater prop, CPR corr
145.98,0.0,0,0                  *rods2 (ax. power distrib. starts at z=0.0 in)
27                              *rods.3: number of ax locations
    0.0,0.000,    2.9196,0.307,   8.7588,1.1011,   14.598,1.3797      *rods.4: z, peaking factor
  20.4372,1.4669, 26.2764,1.4684, 32.1156,1.4391,  37.9548,1.4008
  43.794,1.3614,  49.6332,1.3233, 55.4724,1.2868,  61.3116,1.2513
  67.1508,1.2155, 72.99,1.1779,  78.8292,1.1363,  84.6684,1.0876
  90.5076,1.0358, 96.3468,0.9328,  102.186,0.8836,  108.0252,0.8324
 113.8644,0.7808,  119.7036,0.7239, 125.5428,0.6388, 131.382,0.5119
 137.2212,0.2121,   143.0604, 0.0445,  145.98, 0.00
*
101,101,101                     *rods.29: number of canisters, water rods and pin rods
1,2,                            *rods.30: inside and outside channel # for each canister
3,2,8, 
4,2,8, 
5,8, 
6,8,9, 
7,9, 
10,8,19, 
11,8,9,19 
12,9,20, 
13,9,20, 
14,19, 
15,19,20, 
16,20, 
17,20,21, 
18,21, 
22,19,35, 
23,19,20,35 
24,20,36, 
25,20,21,36 
26,21, 
27,21,37, 
28,35, 
29,35,36, 
30,36, 
31,36,37, 
32,37, 
33,37,38, 
34,38, 
39,35,56, 
40,35,36,56 
41,36,57, 
42,36,37,57 
43,37,58, 
44,37,38,58 
45,38,59, 
46,38,59, 
47,56, 
48,56,57, 
49,57, 
50,57,58, 
51,58, 
52,58,59, 
53,59, 
54,59,60, 
55,60, 
61,56,82, 
62,56,57,82 
63,57,83, 
64,57,58,83 
65,58,84, 
66,58,59,84 
67,59,85, 
68,59,60,85 
69,60,86, 
70,60,86, 
71,82, 
72,82,83, 
73,83, 
74,83,84, 
75,84, 
76,84,85, 
77,85, 
78,85,86, 
79,86, 
80,86,126, 
81,126, 
87,82,107, 
88,82,83,107 
89,83,108, 
90,83,84,108 
91,84,109, 
92,84,85,109 
93,85,110, 
94,85,86,110 
95,86,126, 
96,86,126, 
97,107, 
98,107,108, 
99,108, 
100,108,109, 
101,109, 
102,109,110, 
103,110, 
104,110,126, 
105,126, 
106,126, 
111,107,126, 
112,107,108,126 
113,108,126, 
114,108,109,126 
115,109,126, 
116,109,110,126 
117,110,126, 
118,110,126, 
119,126, 
120,126, 
121,126, 
122,126, 
123,126, 
124,126, 
125,126,
*
11,0.5,0.5                     *rods.31: 11 canisters are modeled as half-canisters
1,3,6,10,15,21,28,36,45,55,66, *rods.32: canister # 
0
*
2,-1                           *rods.33 (canisters are geom 2 type)
0
*
1,-1                           *rods.35
0 
0,0.0                          *rods.37
*
0.5,11                         *rods.39 
1,4,7,13,18,27,34,46,55,70,81, *rods.40
*
1.0,90 
3,5,6,10,11,12,14,15,16,17,22,23 
24,25,26,28,29,30,31,32,33,39,40,41 
42,43,44,45,47,48,49,50,51,52,53,54 
61,62,63,64,65,66,67,68,69,71,72,73 
74,75,76,77,78,79,80,87,88,89,90,91 
92,93,94,95,96,97,98,99,100,101,102,103 
104,105,106,111,112,113,114,115,116,117,118,119 
120,121,122,123,124,125 
0
*
3,-1                           *rods41 (water tubes are geom 3 type)
0 
1,-1                           *rods43
0 
0,0.0                          *rods45
37.0,11                        *rods.47: each half-bundle contains 37 fuel pins
1,4,7,13,18,27,34,46,55,70,81, 
74.0,90                        *rods.47: each full-bundle contains 74 fuel pins 
3,5,6,10,11,12,14,15,16,17,22,23 
24,25,26,28,29,30,31,32,33,39,40,41 
42,43,44,45,47,48,49,50,51,52,53,54 
61,62,63,64,65,66,67,68,69,71,72,73 
74,75,76,77,78,79,80,87,88,89,90,91 
92,93,94,95,96,97,98,99,100,101,102,103 
104,105,106,111,112,113,114,115,116,117,118,119 
120,121,122,123,124,125 
0 
1,101                          *rods49 (lumped bundles are geom 1 type)
203,204,205,206,207,208,209,210,211,212,213,214    *rods50
215,216,217,218,219,220,221,222,223,224,225,226 
227,228,229,230,231,232,233,234,235,236,237,238 
239,240,241,242,243,244,245,246,247,248,249,250 
251,252,253,254,255,256,257,258,259,260,261,262 
263,264,265,266,267,268,269,270,271,272,273,274 
275,276,277,278,279,280,281,282,283,284,285,286 
287,288,289,290,291,292,293,294,295,296,297,298 
299,300,301,302,303 
0 
1,-1                           *rods51
0                           
1                              *rods53
0.9,1.15,1.17,1.18,1.02,1.02,1.02,1.02,0.96,1.01,0.45,1.2        *rods.56: bundle core radial peaking factors
1.05,1.45,0.95,1.28,1.51,0.99,1.46,1.33,1.02,1.42,1.14,1.5 
1.18,1.45,1.06,0.98,1.3,1.46,1.28,1.24,1.02,1.22,1.44,1.0 
1.02,1.26,1.41,0.92,1.37,1.15,0.99,1.09,1.43,1.09,0.95,1.1 
0.97,1.39,1.23,1.4,1.15,1.38,1.02,1.26,0.94,1.28,1.09,1.33 
0.99,1.34,1.0,1.27,1.02,0.76,0.8,1.04,1.22,0.85,0.85,1.0 
1.22,1.03,0.84,0.43,0.74,0.79,0.92,0.79,0.83,1.08,0.93,0.72 
0.43,0.26,0.73,0.73,0.76,0.67,0.73,0.75,0.58,0.39,0.35,0.35 
0.34,0.32,0.35,0.33,0.27 
25                                  *rods59
2,8,9,19,20,21,35,36,37,38,56,57    *rods60 
58,59,60,82,83,84,85,86,107,108,109,110 
126 
0,1.198,0                      *rods.61: max J1 factor for each bundle
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0 
0,1.198,0
*
1,nucl, 0.4400, 0.376, 12, 0.0000, 0.028     *rods62
0, 0, 0, 0, 0,1000.0, 0.960, 0.0000          *rods63
*
2,wall,21.66,0.0000,1                        *rods68
3,1,0.100,0.0                                *rods69 (zircaloy is material 1)
*
3,tube, 1.3859,1.3011,1 
3,1,0.0424,0.0
*
1,1,409.7                                   *rods.70              
662.0,0.076,10.05                           *rods.71: zircalloy properties
*
*******************************
*GROUP OPER
*
oper,1,1,0,1,0,0,0,0,0                      *oper.1: inlet conditions as inlet T and flow rate
-1.0,0.0,2.0,0.0050,0                       *oper.2: all channels have same pressure drop, 2% heat generated in the coolant
1,                                          *oper.3
126,                                        *oper.4: channel 126 is excluded from pressure drop uniformity among channels
*
1035.0000,533.000,3070.217,50, 0.00,   *oper.5:P(psi),inlet T(oF),1/8th core flow(-4% to account for wtr rods),power per fuel rod (kW/rod)
0,0,0,0,0,0                                 *oper.12
*
*******************************
*GROUP CORR
*
corr,-2,2,0                                 *corr.1
epri,epri,epri,none,                        *corr.2 
0.2,                                        *corr.3
ditb,chen,chen,epri,cond,g5.7,              *corr.6 
0,0, 0.0000                                 *corr.16
hnch                                        *corr.18
*
*******************************
*GROUP DRAG
*
drag,2,0,1,                                 *drag.1
0.15139,-0.18,0,94.0485,-1,0,               *drag.2 (Cheng&Todreas coeff for bundles)
0.184,-0.2,0,64,-1,0,                       *drag.2 (McAdams corr for turb flow, T&K page 379)
25,2,                                       *drag.3
2,8,9,19,20,21,35,36,37,38,56,57,58,59,60,82  *drag.4
83,84,85,86,107,108,109,110,126 
0 
0.5,0.562,                                  *drag.5
*
*******************************
*GROUP GRID
*
grid,0,6                            *grid.1: 6 types of form losses (lower tie plate, grids, upper tie plate, central bundle orif.coeff,
*                                            peropheral bundle orif. coeff, crucifotm control rod insertion channel orif. coeff)
9.4609, 1.104, 0.3751,21.089,182.049,270.000   *grid.2: values of the form losses
*
89,10                                       *grid.4: 89 central bundles having form losses at 10 ax. locations
1,3,4,5,6,7,10,11,12,13,14,15,16,17,18,22   *grid.5: channel #
23,24,25,26,27,28,29,30,31,32,33,34,39,40,41,42 
43,44,45,46,47,48,49,50,51,52,53,54,55,61,62,63 
64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79 
80,87,88,89,90,91,92,93,94,95,97,98,99,100,101,102 
103,104,111,112,113,114,115,116,117
*
0.000,4,  7.3,1, 19.5,2, 39.0,2, 58.5,2, 78.0,2, 97.5,2,117.0,2   *grid6 (LTP raised to 7.3", UTP lowered to 143", to account for tot.length=heated length)                                                                  
136.5,2,145.0,3
*
12,10                                                               *grid4 (for peripheral bundles)
81,96,105,106,118,119,120,121,122,123,124,125                       *grid.5: channel #
0.000,5,  7.3,1, 19.5,2, 39.0,2, 58.5,2, 78.0,2, 97.5,2,117.0,2     *grid.6: ax.location, form loss type # 
139.5,2,145.0,3
*
25,1                                             *grid4 (for cruc. control rod channels)
2,8,9,19,20,21,35,36,37,38,56,57,58,59,60,82     *grid.5: channel #
83,84,85,86,107,108,109,110,126 
0.000,6                                          *grid.6: ax.location, form loss type # 
0
*
*******************************
*GROUP CONT
cont                                                                *cont.1
0.0,0,1000,1000,2,0,   *iterative solution                          *cont.2
0.10,0.00001,0.001,0.05,0.01,0.9,1.5,1.0,                           *cont.3
5,0,0,0,0,0,1,1,0,0,0,1,0,0,                                        *cont.6
500.0,0.0,0.0,0.0,0.0,0.0,                                          *cont.7
endd 
0
