UGM 2016 PB2 steadystate

From Serpent Wiki
Revision as of 14:51, 29 November 2017 by Antti Rintala (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search


Non-coupled

Main input

% --- Set simulation title

set title "Peach Bottom 2, steady state"

% --- Link XS-libraries etc.



% --- Include other input files

include "./include/pins"
include "./include/geometry"
include "./include/materials"

% --- Reflective boundary condition in XY

set bc 2 2 1

% --- Switch group constant generation off

set gcu -1

% --- Neutron population and criticality cycles:
%     (too low for physical results)

set pop 5000 100 50

% --- Total power for normalization:

set power 3.9429e+06

% --- Some plots

plot 2 500 1500
plot 3 500 500

mesh 3 500 500
mesh 2 500 1000

mesh 10 3 500 500
mesh 10 2 500 1000

Included files

To be stored in a subfolder called ./include/

include/pins

% --- Fuel Pin definitions:

% --- Empty lattice position:

pin 99
cool

% --- Fuel pins

pin 1
fuel4  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 2
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 3
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 4
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 5
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 6
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 7
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 8
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 9
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 10
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 11
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 12
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 13
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 14
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 15
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 16
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 17
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 18
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 19
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 20
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 21
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 22
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 23
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 24
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 25
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 26
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 27
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 28
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

include/geometry

% --- Pin lattice

lat 100 1  0.23876  0.23876 9 9 1.875
99 99 99 99 99 99 99 99 99
99  1  2  3  4  5  6  7 99
99  2  8  9 10 11 12 13 99
99  3  9 14 15 16 17 18 99
99  4 10 15 19 20 21 22 99
99  5 11 16 20 23 24 25 99
99  6 12 17 21 24 26 27 99
99  7 13 18 22 25 27 28 99
99 99 99 99 99 99 99 99 99

% --- Outer channel (assembly pitch = 15.375):

surf  1  sqc   0.23876  0.23876  6.70306
surf  2  sqc   0.23876  0.23876  6.90626
surf  8  sqc   0.23876  0.23876  6.90626
surf  9  cuboid  -7.62 7.62 -7.62 7.62 -182.88 182.88

% --- Cell definitions:

cell  3  101  fill 100 -1         % Pin lattice
cell  4  101  box       1         % Channel box wall
cell  5  0  fill 101   -8 -9      % assembly inside channel box outer wall
cell  6  0  moder       8 -9      % Water outside channel box

cell 99  0  outside   9           % Outside world

include/materials

% --- Fuel materials:

% Fuel1 2.93 wt % enrichment
mat fuel1   -10.42 tmp 900 rgb 100 0 0
92235.03c   -0.0258279
92238.03c   -0.85567041
8016.03c    -0.11850169

% Fuel2 1.94 wt % enrichment
mat fuel2   -10.42 tmp 900 rgb 160 0 0
92235.03c   -0.01710132
92238.03c   -0.86441006
8016.03c    -0.11848862

% Fuel3 1.69 wt % enrichment
mat fuel3   -10.42 tmp 900 rgb 200 75 75
92235.03c   -0.0148976
92238.03c   -0.86661708
8016.03c    -0.11848532

% Fuel4 1.33 wt % enrichment
mat fuel4   -10.42 tmp 900 rgb 255 150 150
92235.03c   -0.01172421
92238.03c   -0.86979522
8016.03c    -0.11848057

% --- Fuel with Gd:

% Fuel5 2.93 wt % enrichment, 3 wt% Gadcontent
mat fuel5   -10.25 tmp 900 rgb 0 100 0
92235.03c   -0.02505306
92238.03c   -0.8300003
8016.03c    -0.11891791
64152.03c   -0.00005029
64154.03c   -0.00055541
64155.03c   -0.00379519
64156.03c   -0.00528302
64157.03c   -0.004065
64158.03c   -0.00649317
64160.03c   -0.00578667

% --- "Zircaloy-2" [PNNL-15870, Rev. 1]
% Xe18 commented because of lack of XS data

mat Zirc2 -6.56 tmp 600 rgb 200 200 200
 8016.03c  -1.19376E-03
 8017.03c  -4.83282E-07
24050.03c  -4.16117E-05
24052.03c  -8.34483E-04
24053.03c  -9.64457E-05
24054.03c  -2.44600E-05
26054.03c  -5.62862E-05
26056.03c  -9.16258E-04
26057.03c  -2.15389E-05
26058.03c  -2.91667E-06
28058.03c  -3.35317E-04
28060.03c  -1.33612E-04
28061.03c  -5.90496E-06
28062.03c  -1.91358E-05
28064.03c  -5.03067E-06
40090.03c  -4.98111E-01
40091.03c  -1.09835E-01
40092.03c  -1.69731E-01
40094.03c  -1.75753E-01
40096.03c  -2.89183E-02
50112.03c  -1.27668E-04
50114.03c  -8.84175E-05
50115.03c  -4.59485E-05
50116.03c  -1.98205E-03
50117.03c  -1.05596E-03
50118.03c  -3.35857E-03
50119.03c  -1.20129E-03
50120.03c  -4.59450E-03
50122.03c  -6.63830E-04
50124.03c  -8.43778E-04

% --- Assembly box treated as zircaloy
mat box      -6.56000E+00 tmp 556 rgb 200 200 200
 8016.03c  -1.19376E-03
 8017.03c  -4.83282E-07
24050.03c  -4.16117E-05
24052.03c  -8.34483E-04
24053.03c  -9.64457E-05
24054.03c  -2.44600E-05
26054.03c  -5.62862E-05
26056.03c  -9.16258E-04
26057.03c  -2.15389E-05
26058.03c  -2.91667E-06
28058.03c  -3.35317E-04
28060.03c  -1.33612E-04
28061.03c  -5.90496E-06
28062.03c  -1.91358E-05
28064.03c  -5.03067E-06
40090.03c  -4.98111E-01
40091.03c  -1.09835E-01
40092.03c  -1.69731E-01
40094.03c  -1.75753E-01
40096.03c  -2.89183E-02
50112.03c  -1.27668E-04
50114.03c  -8.84175E-05
50115.03c  -4.59485E-05
50116.03c  -1.98205E-03
50117.03c  -1.05596E-03
50118.03c  -3.35857E-03
50119.03c  -1.20129E-03
50120.03c  -4.59450E-03
50122.03c  -6.63830E-04
50124.03c  -8.43778E-04

% --- "Steel, Stainless 304" [PNNL-15870, Rev. 1]

mat ssteel -8.00000E+00 tmp 556 rgb 100 100 100
 6012.03c  -3.95366E-04
% 6013.03c  -4.63372E-06
14028.03c  -4.59332E-03
14029.03c  -2.41681E-04
14030.03c  -1.64994E-04
15031.03c  -2.30000E-04
16032.03c  -1.42073E-04
16033.03c  -1.15681E-06
16034.03c  -6.75336E-06
16036.03c  -1.68255E-08
24050.03c  -7.93000E-03
24052.03c  -1.59029E-01
24053.03c  -1.83798E-02
24054.03c  -4.66139E-03
25055.03c  -1.00000E-02
26054.03c  -3.96166E-02
26056.03c  -6.44901E-01
26057.03c  -1.51600E-02
26058.03c  -2.05287E-03
28058.03c  -6.21579E-02
28060.03c  -2.47678E-02
28061.03c  -1.09461E-03
28062.03c  -3.54721E-03
28064.03c  -9.32539E-04

% --- Coolant (40% void fraction):

mat cool     -0.46072  moder lwtr 1001 tmp 556 rgb 150 150 255
 1001.03c     0.66667
 8016.03c     0.33333

% --- Moderator:

mat moder    -0.46072  moder lwtr 1001 tmp 556 rgb 50 50 255
 1001.03c     0.666667
 8016.03c     0.333333

% --- Thermal scattering data for light water:

%(HinH20 at 573.60K)
therm lwtr lwj3.11t

Coupled

Main input

% --- Set simulation title

set title "Peach Bottom 2 with FINIX, steady state"

% --- Link XS-libraries etc.



% --- Include other input files

include "./include/pins"
include "./include/geometry"
include "./include/materials"

% --- Reflective boundary condition in XY

set bc 2 2 1

% --- Switch group constant generation off

set gcu -1

% --- Neutron population and criticality cycles:
%     (too low for physical results)

set pop 5000 100 50

% --- Number of coupled calculation iterations

set ccmaxiter   3

% --- Use fission source passing with 10 inactive cycles between iterations

set fsp 1 10

% --- Total power for normalization:

set power 3.9429e+06

% --- FINIX file names

set finrodfile "./include/rod_library.txt"
set finoptionsfile "./include/options_library.txt"
set finscenariofile "./include/scenario_library.txt"

% --- FINIX initialization interface file has to be created

%set fininitfile "./Finix0.ifc"

% --- FINIX power binning definitions

finbin powerbins1 0.0 0.60579 5 -182.88 182.88 24

% --- FINIX solver definitions

finrod 1  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 2  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 3  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 4  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 5  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 6  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 7  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 8  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 9  rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 10 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 11 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 12 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 13 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 14 rod_pb2_yesgd options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 15 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 16 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 17 rod_pb2_yesgd options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 18 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 19 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 20 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 21 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 22 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 23 rod_pb2_yesgd options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 24 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 25 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 26 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 1
finrod 27 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 2
finrod 28 rod_pb2_nogd  options_pb2_all scenario_pb2_all powerbins1 nrods 1

% --- Some plots

plot 2 500 1500
plot 3 500 500

mesh 3 500 500
mesh 2 500 1000

mesh 10 3 500 500
mesh 10 2 500 1000

Included files

To be stored in a subfolder called ./include/

include/pins

% --- Fuel Pin definitions:

% --- Empty lattice position:

pin 99
cool

% --- Fuel pins

pin 1
fuel4  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 2
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 3
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 4
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 5
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 6
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 7
fuel3  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 8
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 9
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 10
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 11
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 12
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 13
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 14
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 15
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 16
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 17
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 18
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 19
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 20
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 21
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 22
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 23
fuel5  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 24
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 25
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 26
fuel1  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 27
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

pin 28
fuel2  0.60579
void   0.62103
Zirc2  0.71501
cool

include/geometry

% --- Pin lattice

lat 100 1  0.23876  0.23876 9 9 1.875
99 99 99 99 99 99 99 99 99
99  1  2  3  4  5  6  7 99
99  2  8  9 10 11 12 13 99
99  3  9 14 15 16 17 18 99
99  4 10 15 19 20 21 22 99
99  5 11 16 20 23 24 25 99
99  6 12 17 21 24 26 27 99
99  7 13 18 22 25 27 28 99
99 99 99 99 99 99 99 99 99

% --- Outer channel (assembly pitch = 15.375):

surf  1  sqc   0.23876  0.23876  6.70306
surf  2  sqc   0.23876  0.23876  6.90626
surf  8  sqc   0.23876  0.23876  6.90626
surf  9  cuboid  -7.62 7.62 -7.62 7.62 -182.88 182.88

% --- Cell definitions:

cell  3  101  fill 100 -1         % Pin lattice
cell  4  101  box       1         % Channel box wall
cell  5  0  fill 101   -8 -9      % assembly inside channel box outer wall
cell  6  0  moder       8 -9      % Water outside channel box

cell 99  0  outside   9           % Outside world

include/materials

% --- Fuel materials:

% Fuel1 2.93 wt % enrichment
mat fuel1   -10.42 tft 556 2500 rgb 100 0 0
92235.03c   -0.0258279
92238.03c   -0.85567041
8016.03c    -0.11850169

% Fuel2 1.94 wt % enrichment
mat fuel2   -10.42 tft 556 2500 rgb 160 0 0
92235.03c   -0.01710132
92238.03c   -0.86441006
8016.03c    -0.11848862

% Fuel3 1.69 wt % enrichment
mat fuel3   -10.42 tft 556 2500 rgb 200 75 75
92235.03c   -0.0148976
92238.03c   -0.86661708
8016.03c    -0.11848532

% Fuel4 1.33 wt % enrichment
mat fuel4   -10.42 tft 556 2500 rgb 255 150 150
92235.03c   -0.01172421
92238.03c   -0.86979522
8016.03c    -0.11848057

% --- Fuel with Gd:

% Fuel5 2.93 wt % enrichment, 3 wt% Gadcontent
mat fuel5   -10.25 tft 556 2500 rgb 0 100 0
92235.03c   -0.02505306
92238.03c   -0.8300003
8016.03c    -0.11891791
64152.03c   -0.00005029
64154.03c   -0.00055541
64155.03c   -0.00379519
64156.03c   -0.00528302
64157.03c   -0.004065
64158.03c   -0.00649317
64160.03c   -0.00578667

% --- "Zircaloy-2" [PNNL-15870, Rev. 1]
% Xe18 commented because of lack of XS data

mat Zirc2 -6.56000E+00 tft 556 1000 rgb 200 200 200
 8016.03c  -1.19376E-03
 8017.03c  -4.83282E-07
24050.03c  -4.16117E-05
24052.03c  -8.34483E-04
24053.03c  -9.64457E-05
24054.03c  -2.44600E-05
26054.03c  -5.62862E-05
26056.03c  -9.16258E-04
26057.03c  -2.15389E-05
26058.03c  -2.91667E-06
28058.03c  -3.35317E-04
28060.03c  -1.33612E-04
28061.03c  -5.90496E-06
28062.03c  -1.91358E-05
28064.03c  -5.03067E-06
40090.03c  -4.98111E-01
40091.03c  -1.09835E-01
40092.03c  -1.69731E-01
40094.03c  -1.75753E-01
40096.03c  -2.89183E-02
50112.03c  -1.27668E-04
50114.03c  -8.84175E-05
50115.03c  -4.59485E-05
50116.03c  -1.98205E-03
50117.03c  -1.05596E-03
50118.03c  -3.35857E-03
50119.03c  -1.20129E-03
50120.03c  -4.59450E-03
50122.03c  -6.63830E-04
50124.03c  -8.43778E-04

% --- Assembly box treated as zircaloy
mat box      -6.56000E+00 tmp 556 rgb 200 200 200
 8016.03c  -1.19376E-03
 8017.03c  -4.83282E-07
24050.03c  -4.16117E-05
24052.03c  -8.34483E-04
24053.03c  -9.64457E-05
24054.03c  -2.44600E-05
26054.03c  -5.62862E-05
26056.03c  -9.16258E-04
26057.03c  -2.15389E-05
26058.03c  -2.91667E-06
28058.03c  -3.35317E-04
28060.03c  -1.33612E-04
28061.03c  -5.90496E-06
28062.03c  -1.91358E-05
28064.03c  -5.03067E-06
40090.03c  -4.98111E-01
40091.03c  -1.09835E-01
40092.03c  -1.69731E-01
40094.03c  -1.75753E-01
40096.03c  -2.89183E-02
50112.03c  -1.27668E-04
50114.03c  -8.84175E-05
50115.03c  -4.59485E-05
50116.03c  -1.98205E-03
50117.03c  -1.05596E-03
50118.03c  -3.35857E-03
50119.03c  -1.20129E-03
50120.03c  -4.59450E-03
50122.03c  -6.63830E-04
50124.03c  -8.43778E-04

% --- "Steel, Stainless 304" [PNNL-15870, Rev. 1]

mat ssteel -8.00000E+00 tmp 556 rgb 100 100 100
 6012.03c  -3.95366E-04
% 6013.03c  -4.63372E-06
14028.03c  -4.59332E-03
14029.03c  -2.41681E-04
14030.03c  -1.64994E-04
15031.03c  -2.30000E-04
16032.03c  -1.42073E-04
16033.03c  -1.15681E-06
16034.03c  -6.75336E-06
16036.03c  -1.68255E-08
24050.03c  -7.93000E-03
24052.03c  -1.59029E-01
24053.03c  -1.83798E-02
24054.03c  -4.66139E-03
25055.03c  -1.00000E-02
26054.03c  -3.96166E-02
26056.03c  -6.44901E-01
26057.03c  -1.51600E-02
26058.03c  -2.05287E-03
28058.03c  -6.21579E-02
28060.03c  -2.47678E-02
28061.03c  -1.09461E-03
28062.03c  -3.54721E-03
28064.03c  -9.32539E-04

% --- Coolant (40% void fraction):

mat cool     -0.46072  moder lwtr 1001 tmp 556 rgb 150 150 255
 1001.03c     0.66667
 8016.03c     0.33333

% --- Moderator:

mat moder    -0.46072  moder lwtr 1001 tmp 556 rgb 50 50 255
 1001.03c     0.666667
 8016.03c     0.333333

% --- Thermal scattering data for light water:

%(HinH20 at 573.60K)
therm lwtr lwj3.11t

include/rod_library.txt

!!!!!!                                             !!!!!!
!!!!!! Peach Bottom 2 BWR rod from UAM description !!!!!!
!!!!!! No gadolinia                                !!!!!!
!!!!!!                                             !!!!!!

begin rod_pb2_nogd

! --- Basic rod geometry

pellet_inner_radius = 0.0
pellet_outer_radius = 6.0579e-03
clad_inner_radius = 6.2103e-03
clad_outer_radius = 7.1501e-03
fuel_axial_length = 3.6576
clad_axial_length = 3.6576
plenum_length = 0.4013

! --- Fuel information

pellet_roughness = 2.0e-6
fractional_density = 0.951
gadolinia_weight_fraction = 0.0

! --- Clad information

clad_roughness = 0.5e-6
coldwork = 0.5
clad_oxygen_concentration = 0.0012
fast_neutron_fluence = 0.0

! --- Fill gas pressure (Pa)
fill_gas_pressure = 0.69e6

! --- Initial (fill) plenum temperature (K)
fill_gas_temperature = 300.0

! --- Fill gas fractions
gas_fraction_He = 1.0
gas_fraction_Ar = 0.0
gas_fraction_Kr = 0.0
gas_fraction_Xe = 0.0
gas_fraction_H2 = 0.0
gas_fraction_N2 = 0.0
gas_fraction_H2O = 0.0

! --- Rod pitch
pitch = 18.75e-3

! --- Number of rods in unit cell (1.0 in square, 0.5 in triangular lattice)
unitcell = 1.0

end rod_pb2_nogd


!!!!!!                                             !!!!!!!
!!!!!! Peach Bottom 2 BWR rod from UAM description !!!!!!!
!!!!!! 3.0 wt-% gadolinia                          !!!!!!!
!!!!!!                                             !!!!!!!


begin rod_pb2_yesgd

! --- Basic rod geometry

pellet_inner_radius = 0.0
pellet_outer_radius = 6.0579e-03
clad_inner_radius = 6.2103e-03
clad_outer_radius = 7.1501e-03
fuel_axial_length     = 3.6576
clad_axial_length = 3.6576
plenum_length = 0.4013

! --- Fuel information

pellet_roughness = 2.0e-6
fractional_density = 0.951
gadolinia_weight_fraction = 0.03

! --- Clad information

clad_roughness = 0.5e-6
coldwork = 0.5
clad_oxygen_concentration = 0.0012
fast_neutron_fluence = 0.0

! --- Fill gas pressure (Pa)
fill_gas_pressure = 0.69e6

! --- Initial (fill) plenum temperature (K)
fill_gas_temperature = 300.0

! --- Fill gas fractions
gas_fraction_He = 1.0
gas_fraction_Ar = 0.0
gas_fraction_Kr = 0.0
gas_fraction_Xe = 0.0
gas_fraction_H2 = 0.0
gas_fraction_N2 = 0.0
gas_fraction_H2O = 0.0

! --- Rod pitch
pitch = 18.75e-3

! --- Number of rods in unit cell (1.0 in square, 0.5 in triangular lattice)
unitcell = 1.0

end rod_pb2_yesgd

include/options_library.txt

!!!!!!                                                 !!!!!!
!!!!!! Basic options for the Peach Bottom 2 simulation !!!!!!
!!!!!!                                                 !!!!!!

begin options_pb2_all

! --- Number of nodes in the thermal/mechanical solution

axial_nodes = 24
pellet_radial_nodes = 51
clad_radial_nodes = 21

! --- Boundary condition is fixed cladding outer temperature
!     NB: This must be given in the scenario

boundary_option = 0

end options_pb2_all

!!!!!!                                                                    !!!!!!
!!!!!! Basic options for the Peach Bottom 2 simulation with 5 axial zones !!!!!!
!!!!!!                                                                    !!!!!!

begin options_pb2_5axial

! --- Number of nodes in the thermal/mechanical solution

axial_nodes = 5
pellet_radial_nodes = 51
clad_radial_nodes = 21

! --- Boundary condition is fixed cladding outer temperature
!     NB: This must be given in the scenario

boundary_option = 0

end options_pb2_5axial

include/scenario_library.txt

!!!!!!
! --- Basic scenario for Peach Bottom 2 simulations using clad surface temperature
!     as a boundary condition
!!!!!!

begin scenario_pb2_all

! --- Clad temperature history 

! --- Single axial zone
 
clad_temperature_history_zones = 3.6576

! --- Temperature, time pairs for first axial zone
!     Only one time point now, only one axial zone

clad_temperature_history(1) = 557.0, 0.0

end scenario_pb2_all