% --- NEA BUC Criticality Safety Benchmark Phase III-C
set title "NEA/NSC/R/(2015)6"
% -----------------------------------------------------------------------------
% --- Geometry ---------------------------------------------------------------
% --- Fuel Pin definitions:
pin 1
fuel1 0.48
clad 0.55
cool40
pin 2
fuel2 0.48
clad 0.55
cool40
pin 3
fuel3 0.48
clad 0.55
cool40
pin 4
fuel4 0.48
clad 0.55
cool40
pin 5
fuel5 0.48
clad 0.55
cool40
pin G
fuelG 0.48
clad 0.55
cool40
% --- Empty lattice position:
pin W
cool40
% --- Lattice (type = 1, pin pitch = 1.45):
lat 10 1 0.0 0.0 11 11 1.45
W W W W W W W W W W W
W 5 4 3 2 2 W W W W W
W W 3 G 2 G W W W W W
W W W 1 2 2 W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
W W W W W W W W W W W
% --- 45 degree symmetry
set usym 10 3 2 0.0 0.0 225 45
% --- Outer channel & box wall (assembly pitch = 15.24):
surf 1 sqc 0.0 0.0 6.700
surf 2 sqc 0.0 0.0 6.950
surf 3 sqc 0.0 0.0 7.620
% --- Water channel:
surf 4 sqc 0.0 0.0 1.850
surf 5 sqc 0.0 0.0 1.925
% --- Cell definitions:
cell 1 0 moder -4 % Inner water channel
cell 2 0 box 4 -5 % Inner water channel wall
cell 3 0 fill 10 -1 5 % Pin lattice
cell 4 0 box 1 -2 % Channel box wall
cell 5 0 moder 2 -3 % Outer water cnallel
cell 99 0 outside 3 % Outside world
% --- Reflective boundary condition:
set bc 2
% -----------------------------------------------------------------------------
% --- Material compositions ---------------------------------------------------
% --- Fuel type 1: 4.9wt% U235
mat fuel1 sum burn 1 rgb 227 7 19
92234.09c 8.5649E-6
92235.09c 1.1221E-3
92238.09c 2.1495E-2
8016.09c 4.5252E-2
% --- Fuel type 2: 4.4wt% U235
mat fuel2 sum burn 1 rgb 238 130 98
92234.09c 7.6910E-6
92235.09c 1.0076E-3
92238.09c 2.1609E-2
8016.09c 4.5249E-2
% --- Fuel type 3: 3.9wt% U235
mat fuel3 sum burn 1 rgb 220 227 131
92234.09c 6.8170E-6
92235.09c 8.9315E-4
92238.09c 2.1723E-2
8016.09c 4.5247E-2
% --- Fuel type 4: 3.4wt% U235
mat fuel4 sum burn 1 rgb 212 190 202
92234.09c 5.9431E-6
92235.09c 7.7865E-4
92238.09c 2.1838E-2
8016.09c 4.5244E-2
% --- Fuel type 5: 2.1wt% U235
mat fuel5 sum burn 1 rgb 182 178 157
92234.09c 3.6708E-6
92235.09c 4.8094E-4
92238.09c 2.2134E-2
8016.09c 4.5238E-2
% --- Fuel type G: 3.4wt% U235, 5.0wt% Gd
mat fuelG sum burn 1 rgb 7 227 19
92234.09c 5.5429E-6
92235.09c 7.2622E-4
92238.09c 2.0367E-2
8016.09c 4.4678E-2
64154.09c 3.6051E-5
64155.09c 2.4475E-4
64156.09c 3.3852E-4
64157.09c 2.5881E-4
64158.09c 4.1079E-4
64160.09c 3.6481E-4
% --- Division into burnup zones
div fuel1 sep 1
div fuel2 sep 1
div fuel3 sep 1
div fuel4 sep 1
div fuel5 sep 1
div fuelG sep 1 subr 10 0.0 0.48
% --- Fuel volumes:
set mvol
fuel1 1 2.8953
fuel2 5 2.8953
fuel2 4 5.7906
fuel2 3 5.7906
fuel2 2 2.8953
fuel2 1 5.7906
fuel3 2 2.8953
fuel3 1 5.7906
fuel4 1 5.7906
fuel5 1 2.8953
fuelG 20 0.2895
fuelG 19 0.2895
fuelG 18 0.2895
fuelG 17 0.2895
fuelG 16 0.2895
fuelG 15 0.2895
fuelG 14 0.2895
fuelG 13 0.2895
fuelG 12 0.2895
fuelG 11 0.2895
fuelG 10 0.5791
fuelG 9 0.5791
fuelG 8 0.5791
fuelG 7 0.5791
fuelG 6 0.5791
fuelG 5 0.5791
fuelG 4 0.5791
fuelG 3 0.5791
fuelG 2 0.5791
fuelG 1 0.5791
% --- Cladding and channel box wall:
mat clad sum tmp 559 rgb 255 251 220
50000.03c 4.9797E-4
26000.03c 9.1782E-5
24000.03c 7.5861E-5
28000.03c 4.0314E-5
40000.03c 4.2465E-2
mat box sum tmp 559 rgb 255 251 220
50000.03c 4.9797E-4
26000.03c 9.1782E-5
24000.03c 7.5861E-5
28000.03c 4.0314E-5
40000.03c 4.2465E-2
% --- Coolant with different void fractions
mat cool0 sum tmp 559 moder lwtr 1001 rgb 255 255 255
1001.03c 4.9316E-2
8016.03c 2.4658E-2
mat cool40 sum tmp 559 moder lwtr 1001 rgb 255 255 255
1001.03c 3.0588E-2
8016.03c 1.5294E-2
mat cool70 sum tmp 559 moder lwtr 1001 rgb 255 255 255
1001.03c 1.6542E-2
8016.03c 8.2712E-2
% --- Moderator:
mat moder sum tmp 559 moder lwtr 1001 rgb 161 218 247
1001.03c 4.9316E-2
8016.03c 2.4658E-2
% --- Thermal scattering data for light water:
therm lwtr 559 lwj3.09t lwj3.11t
% -----------------------------------------------------------------------------
% --- Options & data libraries ------------------------------------------------
% --- Libraries:
set acelib "sss_jeff311u.xsdata"
set declib "sss_jeff31.dec"
set nfylib "sss_jeff31.nfy"
% --- Neutron population and criticality cycles:
set pop 2000 500 20
% --- Geometry and mesh plots:
plot 3 500 500
mesh 3 500 500
% -----------------------------------------------------------------------------
% --- Burnup calculation ------------------------------------------------------
% --- Power density for normalization
set powdens 25.3E-3
% --- Predictor-corrector algorithm & sub-steps
set pcc LELI 10 10
% --- Depletion history
dep butot
0.2 0.5
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
20 25 30 35 40 45 50
set inventory
922340
922350
922360
922380
932370
942380
942390
942400
942410
942420
952410
952430
962440
380900
420950
430990
441010
451030
471090
531290
541310
551330
551340
551370
621470
621490
621500
621510
621520
631530
631540
631550
641550
641560
641570
641580
% -----------------------------------------------------------------------------
