3D BWR assembly transient example

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 2

% --- 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:

% 10000 neutrons in 10 batches (1000 neutrons per batch)

set nps 4000000000 800000 simutime

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