2D CANDU fuel bundle geometry

<nowiki>

% --- CANDU cluster ------------------------------------------

set title "CANDU"

% --- Fuel pin:

pin 1 fuel 0.6122 clad 0.6540 cool

% --- Lattice (type = 4, 4 rings, 3rd ring rotated 15 deg.):

lat 10 4 0.0 0.0 4

1  0.0000   0.0  1
6  1.4885   0.0  1 1 1 1 1 1

12 2.8755 15.0 1 1 1 1 1 1 1 1 1 1 1 1 18 4.3305 0.0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

% --- Surfaces (core pitch = 18.191 cm):

surf 1 cyl 0.0 0.0 5.16890  % Pressure tube inner wall surf 2 cyl 0.0 0.0 5.60320  % Pressure tube outer wall surf 3 cyl 0.0 0.0 6.44780  % Calandria tube inner wall surf 4 cyl 0.0 0.0 6.58750  % Calandria tube outer wall surf 5 sqc 0.0 0.0 9.09570  % Outer boundary

% --- Cells:

cell 1 0 fill 10 -1  % Pin lattice cell 2 0 tube 1 -2  % Pressure tube cell 3 0 void 2 -3  % Void between tubes cell 4 0 caltube 3 -4  % Calandria tube cell 5 0 moder 4 -5  % Moderator channel cell 6 0 outside 5  % Outside world

% --- Fuel (UO2, natural uranium, 0.7% U-235):

mat fuel -10.4375010

8016.09c    -1.18473E+1 

92235.09c -6.27118E-1 92238.09c -8.75256E+1

% --- Cladding:

mat clad -6.44 25055.06c -1.60000E-1 28000.06c -6.00000E-2 24000.06c -1.10000E-1 40000.06c -9.97100E+1

5010.06c    -5.7409e-05
5011.06c    -2.5259E-04

% --- Pressure tube:

mat tube -6.57 40000.06c -9.75000E+1

5010.06c    -3.8889E-05
5011.06c    -1.7111E-04

% --- Calandria tube:

mat caltube -6.44 25055.06c -1.60000E-1 28000.06c -6.00000E-2 24000.06c -1.10000E-1 40000.06c -9.97100E+1

5010.06c    -5.7409e-05
5011.06c    -2.5259E-04

% --- Coolant water:

mat cool -0.812120 moder lwtr 1001 moder hwtr 1002

8016.06c    -7.99449E-1 
1002.06c    -1.99768E-1
1001.06c    -7.83774E-4

% --- Moderator water:

mat moder -1.082885 moder lwtr 1001 moder hwtr 1002

8016.06c    -7.98895E-1
1002.06c    -2.01016E-1
1001.06c    -8.96000E-5

% --- Thermal scattering data for light and heavy water:

therm lwtr lwj3.11t therm hwtr hwj3.11t

% --- Cross section data library file path:

set acelib "sss_jeff311u.xsdata"

% --- Periodic boundary condition:

set bc 3

% --- Neutron population and criticality cycles:

set pop 2000 500 20

% --- Geometry and mesh plots:

plot 3 500 500 mesh 3 500 500

% ------------------------------------------------------------