Difference between revisions of "Regular Hex Mesh coolant IFC example"
From Serpent Wiki
(→coolant.ifc) |
(→Main input) |
||
Line 46: | Line 46: | ||
Sn-124.03c 2.7897e-05 | Sn-124.03c 2.7897e-05 | ||
− | mat cool sum rgb 200 200 255 | + | mat cool sum rgb 200 200 255 moder HinWater 1001 |
H-1.03c 4.9457e-02 | H-1.03c 4.9457e-02 | ||
H-2.03c 7.4196e-06 | H-2.03c 7.4196e-06 | ||
Line 52: | Line 52: | ||
O-17.03c 9.3982e-06 | O-17.03c 9.3982e-06 | ||
%O-18.03c 5.0701e-05 | %O-18.03c 5.0701e-05 | ||
+ | |||
+ | % --- Use on-the-fly interpolation for thermal scattering data | ||
+ | % lwj3.07t = 474 K / JEFF3.1.1 | ||
+ | % lwj3.09t = 524 K / JEFF3.1.1 | ||
+ | % lwj3.11t = 574 K / JEFF3.1.1 | ||
+ | |||
+ | therm HinWater 0 lwj3.07t lwj3.09t lwj3.11t | ||
%%%%%%%%%%%%%%%%%%%%%%%%%%% | %%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
Line 57: | Line 64: | ||
%%%%%%%%%%%%%%%%%%%%%%%%%%% | %%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
− | surf s01 | + | surf s01 hexyc 0 0 7 |
surf s02 pz -50 | surf s02 pz -50 | ||
surf s03 pz 50 | surf s03 pz 50 | ||
Line 126: | Line 133: | ||
plot 3 500 500 | plot 3 500 500 | ||
plot 1 500 500 | plot 1 500 500 | ||
+ | |||
+ | % --- Sample temperature and density data in 20 points on z-axis | ||
+ | |||
+ | sample 1 0 0 1 0 0 10 -49.9 49.9 | ||
% --- Fission rate / thermal flux plot | % --- Fission rate / thermal flux plot | ||
Line 139: | Line 150: | ||
% --- Include interface for assembly-wise fuel temperature | % --- Include interface for assembly-wise fuel temperature | ||
− | ifc "./coolant.ifc"</nowiki> | + | ifc "./coolant.ifc" |
+ | </nowiki> | ||
=== coolant.ifc === | === coolant.ifc === |
Revision as of 12:24, 8 March 2018
Input example for single assembly regular mesh based interface for coolant 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 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 sum rgb 200 200 255 moder HinWater 1001 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 % --- Use on-the-fly interpolation for thermal scattering data % lwj3.07t = 474 K / JEFF3.1.1 % lwj3.09t = 524 K / JEFF3.1.1 % lwj3.11t = 574 K / JEFF3.1.1 therm HinWater 0 lwj3.07t lwj3.09t lwj3.11t %%%%%%%%%%%%%%%%%%%%%%%%%%% % --- Geometry definition % %%%%%%%%%%%%%%%%%%%%%%%%%%% surf s01 hexyc 0 0 7 surf s02 pz -50 surf s03 pz 50 surf sINF inf % --- Core cell c01 0 fill lPIN -s01 s02 -s03 cell c02 0 outside s01 s02 -s03 cell c03 0 outside -s02 cell c04 0 outside s03 % --- 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 % --- Sample temperature and density data in 20 points on z-axis sample 1 0 0 1 0 0 10 -49.9 49.9 % --- 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 "./coolant.ifc"
coolant.ifc
The interface consists of five axial layers from bottom to top. Here the mesh type 4 is similar to the lattice type 2 in the input (used for the lPIN lattice), although a cartesian mesh could have also been used as there is no variation in the horizontal direction. We let the coolant density decrease from 1 g/cm3 to 0.6 g/cm3 while increasing the coolant temperature from 520 K to 540 K.
2 cool 0 4 0 0 14 -50.0 50.0 1 1 5 -1.0 520 -0.9 525 -0.8 530 -0.7 535 -0.6 540