UGM 2016 PB2 steadystate
Revision as of 14:51, 29 November 2017 by Antti Rintala (talk | contribs)
Contents
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
