Regular Hex Mesh fuel IFC example

From Serpent Wiki
Revision as of 15:10, 3 October 2016 by Ville Valtavirta (Talk | contribs) (Output plots)

Jump to: navigation, search


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). Fuel density is not modified, it is kept at a constant 10.1 g/cm3.

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

Output plots

input_geom1.png
input_geom2.png
input_mesh1png
input_mesh2.png
input_mesh3.png
input_mesh4.png