Regular Hex Mesh fuel IFC example
Input example for regular mesh based interface using hexagonal mesh.
Remember to add cross section libraries to the main input using set acelib
Files
Main input
set title "Regular hex mesh interface input"
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Material definitions %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mat fuel -10.31341 rgb 255 150 150
U-234.03c 5.0131e-06
U-235.03c 5.7503e-04
U-238.03c 2.2625e-02
O-16.03c 4.5895e-02
O-17.03c 1.7482e-05
%O-18.03c 9.4313e-05
mat zirc -6.55 rgb 100 100 100
O-16.03c 3.0743e-04
O-17.03c 1.1711e-07
%O-18.03c 6.3176e-07
Cr-50.03c 3.2962e-06
Cr-52.03c 6.3564e-05
Cr-53.03c 7.2076e-06
Cr-54.03c 1.7941e-06
Fe-54.03c 8.6699e-06
Fe-56.03c 1.3610e-04
Fe-57.03c 3.1431e-06
Fe-58.03c 4.1829e-07
Zr-90.03c 2.1827e-02
Zr-91.03c 4.7600e-03
Zr-92.03c 7.2758e-03
Zr-94.03c 7.3734e-03
Zr-96.03c 1.1879e-03
Sn-112.03c 4.6735e-06
Sn-114.03c 3.1799e-06
Sn-115.03c 1.6381e-06
Sn-116.03c 7.0055e-05
Sn-117.03c 3.7003e-05
Sn-118.03c 1.1669e-04
Sn-119.03c 4.1387e-05
Sn-120.03c 1.5697e-04
Sn-122.03c 2.2308e-05
Sn-124.03c 2.7897e-05
mat cool -0.740582 rgb 200 200 255
H-1.03c 4.9457e-02
H-2.03c 7.4196e-06
O-16.03c 2.4672e-02
O-17.03c 9.3982e-06
%O-18.03c 5.0701e-05
%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Geometry definition %
%%%%%%%%%%%%%%%%%%%%%%%%%%%
surf s01 cylz 0 0 30
surf s02 pz -50
surf s03 pz 50
surf sINF inf
% --- Core
cell c01 0 fill lASS -s01
cell c02 0 outside s01 s02 -s03
cell c03 0 outside -s02
cell c04 0 outside s03
% --- Assembly lattice
lat lASS 3 0.0 0.0 5 5 15.0
uWAT uWAT uWAT uWAT uWAT
uWAT uWAT lPIN lPIN uWAT
uWAT lPIN lPIN lPIN uWAT
uWAT lPIN lPIN uWAT uWAT
uWAT uWAT uWAT uWAT uWAT
% --- Water assembly
cell cWAT uWAT cool -sINF
% --- Fuel assembly pin lattice
lat lPIN 2 0.0 0.0 15 15 1.22
ww ww ww ww ww ww ww ww ww ww ww ww ww ww ww
ww ww ww ww ww ww ww PF PF PF PF PF PF PF ww
ww ww ww ww ww ww PF PF PF PF PF PF PF PF ww
ww ww ww ww ww PF PF PF PF PF PF PF PF PF ww
ww ww ww ww PF PF PF PF PF PF PF PF PF PF ww
ww ww ww PF PF PF PF PF PF PF PF PF PF PF ww
ww ww PF PF PF PF PF PF PF PF PF PF PF PF ww
ww PF PF PF PF PF PF IT PF PF PF PF PF PF ww
ww PF PF PF PF PF PF PF PF PF PF PF PF ww ww
ww PF PF PF PF PF PF PF PF PF PF PF ww ww ww
ww PF PF PF PF PF PF PF PF PF PF ww ww ww ww
ww PF PF PF PF PF PF PF PF PF ww ww ww ww ww
ww PF PF PF PF PF PF PF PF ww ww ww ww ww ww
ww PF PF PF PF PF PF PF ww ww ww ww ww ww ww
ww ww ww ww ww ww ww ww ww ww ww ww ww ww ww
% --- Pin definitions
% Empty lattice position
pin ww
cool
% Fuel pin
pin PF
void 0.07
fuel 0.3765
void 0.3865
zirc 0.4575
cool
% Instrumentation tube
pin IT
cool 0.3850
zirc 0.4400
cool
%%%%%%%%%%%%%%%%%%%%%%%%
% --- Some run options %
%%%%%%%%%%%%%%%%%%%%%%%%
set pop 5000 1000 50
% --- Cross section libraries
%set acelib ""
% --- Geometry plot
plot 3 500 500
plot 1 500 500
% --- Fission rate / thermal flux plot
mesh 3 500 500
mesh 1 500 500
% --- Interface temperature plot
mesh 10 3 500 500
mesh 10 1 500 500
% --- Include interface for assembly-wise fuel temperature
ifc "./assemblywise.ifc"
assemblywise.ifc
The interface consists of two axial layers, first bottom then top. Here the mesh type 5 is similar to the lattice type 3 in the input (used for the lASS lattice).
2 fuel 1
assemblywise.out 2 -50.0 50.0 1
5
0 0 15 -50.0 50.0 5 5 2
-10.1 510 -10.1 530 -10.1 550 -10.1 570 -10.1 590
-10.1 610 -10.1 630 -10.1 650 -10.1 670 -10.1 690
-10.1 710 -10.1 730 -10.1 750 -10.1 770 -10.1 790
-10.1 810 -10.1 830 -10.1 850 -10.1 870 -10.1 890
-10.1 910 -10.1 930 -10.1 950 -10.1 970 -10.1 990
-10.1 1010 -10.1 1030 -10.1 1050 -10.1 1070 -10.1 1090
-10.1 1110 -10.1 1130 -10.1 1150 -10.1 1170 -10.1 1190
-10.1 1210 -10.1 1230 -10.1 1250 -10.1 1270 -10.1 1290
-10.1 1310 -10.1 1330 -10.1 1350 -10.1 1370 -10.1 1390
-10.1 1410 -10.1 1430 -10.1 1450 -10.1 1470 -10.1 1490
