3D BWR assembly transient example
(Redirected from 3D PWR assembly transient example)
Note: For transient simulations, the initial system needs to be as close to criticality as possible. The exact multiplication factor of the system may vary if different cross section libraries is used. In that case, adjust the boron content in the cool and moder materials to achieve criticality.
Contents
Generating the steady state source distributions
% --- Asymmetric BWR assembly with Gd-pins ------------------- set title "Peach Bottom 2 from UAM, source generation " include pins include geometry include materials % --- Reflective boundary condition in XY: set bc 2 2 1 % --- Neutron population and criticality cycles: set pop 80000 1000000 200 1.0 100 % --- Geometry plots: plot 3 1000 1000 0.0 plot 3 1000 1000 0.0 -7.62 7.62 -7.62 7.62 plot 2 500 1500 % --- Total power for normalization: set power 4.0000e+06 % --- Save source to "./source"-files set savesrc "./source" 1.00 0.25 1 1 1 % --- Do not generate group constants set gcu -1 % --- Do not use implicit capture, nxn, or fission set impl 0 0 0
Time dependent simulation (no reactivity insertion)
% --- Asymmetric BWR assembly with Gd-pins ------------------- set title "Peach Bottom 2 from UAM, time dependent simulation" include pins include geometry include materials % --- Reflective boundary condition in XY: set bc 2 2 1 % --- Increase size of precursor buffer set pbuf 10 % --- Total power for normalization: set power 4.0000e+06 % --- Set up time structures for simulation and detectors: % 1 time interval for simulation (no population control) tme simutime 2 1 0 1 % % 20 time bins for detector tme dettime 2 20 0 1 % --- Neutron population: set nps 200000000 40000 simutime % --- Output after each 10 cycles set outp 10 % --- Link source (use point-wise precursor tracking) set dynsrc "./source" 1 % --- Set up detectors % --- Neutron population as a function of time det 1 dr -15 void di dettime % --- Fission energy deposition as a function of time det 2 dr -8 void di dettime % --- Do not generate group constants set gcu -1 % --- Do not use implicit capture, nxn, or fission set impl 0 0 0
Additional files
pins
% --- 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
geometry
% --- Assembly lattice definition 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 4 pz -182.88 % Active stack length 365.76 cm surf 5 pz 182.88 % Active stack length 365.76 cm surf 9 cuboid -7.62 7.62 -7.62 7.62 -250.0 250.0 % --- Cell definitions: cell 3 101 fill 100 -1 % Pin lattice cell 4 101 box 1 % Channel box wall cell 5 0 fill 101 -2 4 -5 % assembly inside channel box outer wall cell 6 0 moder 2 4 -5 -9 % Water outside channel box cell 9 0 moder -4 -9 % bottom reflector cell 10 0 moder 5 -9 % top reflector cell 99 0 outside 9 % Outside world
materials
% --- Fuel materials: % Fuel1 2.93 wt % enrichment mat fuel1 -10.42 rgb 100 0 0 92235.09c -0.0258279 92238.09c -0.85567041 8016.09c -0.11850169 % Fuel2 1.94 wt % enrichment mat fuel2 -10.42 rgb 160 0 0 92235.09c -0.01710132 92238.09c -0.86441006 8016.09c -0.11848862 % Fuel3 1.69 wt % enrichment mat fuel3 -10.42 rgb 200 75 75 92235.09c -0.0148976 92238.09c -0.86661708 8016.09c -0.11848532 % Fuel4 1.33 wt % enrichment mat fuel4 -10.42 rgb 255 150 150 92235.09c -0.01172421 92238.09c -0.86979522 8016.09c -0.11848057 % --- Fuel with Gd: % Fuel5 2.93 wt % enrichment, 3.0 wt% Gadcontent mat fuel5 -10.25 rgb 0 100 0 92235.09c -0.02505306 92238.09c -0.8300003 8016.09c -0.11891791 64152.09c -0.00005029 64154.09c -0.00055541 64155.09c -0.00379519 64156.09c -0.00528302 64157.09c -0.004065 64158.09c -0.00649317 64160.09c -0.00578667 % --- Zirc2 and channel box wall: % --- "Zircaloy-2" [PNNL-15870, Rev. 1] mat Zirc2 -6.56000E+00 rgb 200 200 200 %*OK 8016.06c -1.19376E-03 8017.06c -4.83282E-07 24050.06c -4.16117E-05 24052.06c -8.34483E-04 24053.06c -9.64457E-05 24054.06c -2.44600E-05 26054.06c -5.62862E-05 26056.06c -9.16258E-04 26057.06c -2.15389E-05 26058.06c -2.91667E-06 28058.06c -3.35317E-04 28060.06c -1.33612E-04 28061.06c -5.90496E-06 28062.06c -1.91358E-05 28064.06c -5.03067E-06 40090.06c -4.98111E-01 40091.06c -1.09835E-01 40092.06c -1.69731E-01 40094.06c -1.75753E-01 40096.06c -2.89183E-02 50112.06c -1.27668E-04 50114.06c -8.84175E-05 50115.06c -4.59485E-05 50116.06c -1.98205E-03 50117.06c -1.05596E-03 50118.06c -3.35857E-03 50119.06c -1.20129E-03 50120.06c -4.59450E-03 50122.06c -6.63830E-04 50124.06c -8.43778E-04 % Assembly box treated as zircaloy mat box -6.56000E+00 tmp 557.0 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 557.0 rgb 100 100 100 6012.03c -3.95366E-04 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 557 rgb 150 150 255 1001.03c 0.66667 8016.03c 0.33333 5010.03c 0.0000575 % --- Moderator: mat moder -0.46072 moder lwtr 1001 tmp 557 rgb 50 50 255 1001.03c 0.666667 8016.03c 0.333333 5010.03c 0.0000575 % --- Thermal scattering data for light water: %(HinH20 at 573.60K) therm lwtr lwj3.11t
