Regular Hex Mesh coolant IFC example

From Serpent Wiki
Revision as of 11:23, 8 March 2018 by Ville Valtavirta (Talk | contribs) (coolant.ifc)

Jump to: navigation, search

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
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  hexxc 0 0 30
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

% --- 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