Regular Hex Mesh IFC Example 1

From Serpent Wiki
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 s02 -s03
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: Geometry plot in XY.
input_geom2.png: Geometry plot in YZ.
input_mesh1.png: Fission rate/thermal flux mesh plot in XY.
input_mesh2.png: Fission rate/thermal flux mesh plot in YZ.
input_mesh3.png: Interface material temperature mesh plot in XY.
input_mesh4.png: Interface material temperature mesh plot in YZ.