Difference between revisions of "Regular Hex Mesh coolant IFC example"

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/cm³ to 0.6 g/cm³ 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