2D PWR assembly burnup example (division with div-card)

From Serpent Wiki
Revision as of 13:27, 29 November 2017 by Antti Rintala (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search
set title "PWR Burnup Calculation Based on NEA Benchmark"

% --- Fuel pins:

pin 1
UO2     0.4025
clad    0.4750
water

% --- Gd-pins:

pin 2
UO2Gd   0.4025
clad    0.4750
water

% --- Guide tube:

pin 3
water   0.5730
tube    0.6130
water

% --- Empty lattice position:

pin 4
water

% --- Pin lattice:

lat 110 1 5.06 5.06 9 9 1.265
                   
3 4 4 4 4 4 4 4 4  
1 1 4 4 4 4 4 4 4  
1 1 2 4 4 4 4 4 4  
3 1 1 3 4 4 4 4 4  
1 1 1 1 1 4 4 4 4  
2 1 1 1 1 3 4 4 4  
3 1 1 3 2 1 1 4 4  
1 1 1 1 1 1 1 1 4  
1 1 1 1 1 1 1 1 1  

% --- Symmetry:

set usym 110 3 2 0.0 0.0 45 45

% --- assembly data:

surf  1000  sqc  0.0 0.0 10.752
surf  1001  sqc  0.0 0.0 10.806

cell 110  0  fill 110   -1000
cell 111  0  water       1000 -1001
cell 112  0  outside     1001

% --- Materials:

div UO2 sep 1

mat UO2    6.7402E-02  burn 1
92234.09c  9.1361E-06
92235.09c  9.3472E-04
92238.09c  2.1523E-02
 8016.09c  4.4935E-02

div UO2Gd sep 1 subr 10 0.0 0.4025

mat UO2Gd  6.8366E-02  burn 1
92234.09c  4.2940E-06
92235.09c  5.6226E-04
92238.09c  2.0549E-02
64154.09c  4.6173E-05
64155.09c  2.9711E-04
64156.09c  4.1355E-04
64157.09c  3.1518E-04
64158.09c  4.9786E-04
64160.09c  4.3764E-04
 8016.09c  4.5243E-02

mat clad   3.8510E-02
26000.06c  1.3225E-04
24000.06c  6.7643E-05
40000.06c  3.8310E-02

mat tube   4.3206E-02
26000.06c  1.4838E-04
24000.06c  7.5891E-05
40000.06c  4.2982E-02

mat water  7.2216E-02  moder lwtr 1001
 1001.06c  4.8132E-02
 8016.06c  2.4066E-02
 5010.06c  3.6487E-06
 5011.06c  1.4686E-05

therm lwtr lwj3.11t

% --- Cross section library file path:

set acelib "/xs/sss_endfb7u.xsdata"

% --- Periodic boundary condition:

set bc 3

% --- Neutron population and criticality cycles:

set pop 6000 500 20

% --- Geometry and mesh plots:

plot 3 500 500
mesh 3 500 500

% --- Decay and fission yield libraries:

set declib "/xs/sss_endfb7.dec"
set nfylib "/xs/sss_endfb7.nfy"

% --- Irradiation cycle:

set powdens 38.6E-3

dep butot

    0.10000
    0.50000
    1.00000
    1.50000
    2.00000
    2.50000
    3.00000
    3.50000
    4.00000
    4.50000
    5.00000
    5.50000
    6.00000
    6.50000
    7.00000
    7.50000
    8.00000
    8.50000
    9.00000
    9.50000
   10.00000
   10.50000
   11.00000
   11.50000
   12.00000
   12.50000
   13.00000
   13.50000
   14.00000
   14.50000
   15.00000
   17.50000
   20.00000
   22.50000
   25.00000
   27.50000
   30.00000
   32.50000
   35.00000
   37.50000
   40.00000

% --- Nuclide inventory:

set inventory

 922340
 922350
 922360
 922370
 922380
 922390
 932360
 932370
 932380
 932390
 942360
 942380
 942390
 942400
 942410
 942420
 942430
 952410                                       
 952420
 952430
 952440
 952421
 962420
 962430
 962440
 962450
 962460
 962470
 962480
 962490
 972490
 972500
 982490
 982500
 982510
 982520
 360830
 451030
 451050
 471090
 531350
 541310
 541350
 551330
 551340
 551350
 551370
 561400
 571400
 601430
 601450
 611470
 611480
 611490
 611481
 621470
 621490
 621500
 621510
 621520
 631530
 631540
 631550
 631560
 641520
 641540
 641550
 641560
 641570
 641600

% --- Material volumes:

set mvol

UO2    36   2.03583
UO2    35   4.07166
UO2    34   4.07166
UO2    33   4.07166
UO2    32   4.07166
UO2    31   4.07166
UO2    30   4.07166
UO2    29   4.07166
UO2    28   2.03583
UO2    27   2.03583
UO2    26   4.07166
UO2    25   4.07166
UO2    24   4.07166
UO2    23   4.07166
UO2    22   4.07166
UO2    21   4.07166
UO2    20   2.03583
UO2    19   2.03583
UO2    18   4.07166
UO2    17   4.07166
UO2    16   4.07166
UO2    15   4.07166
UO2    14   4.07166
UO2    13   4.07166
UO2    12   4.07166
UO2    11   2.03583
UO2    10   4.07166
UO2     9   4.07166
UO2     8   4.07166
UO2     7   2.03583
UO2     6   4.07166
UO2     5   4.07166
UO2     4   4.07166
UO2     3   2.03583
UO2     2   2.03583
UO2     1   2.03583
UO2Gd  30   0.40717
UO2Gd  29   0.40717
UO2Gd  28   0.40717
UO2Gd  27   0.40717
UO2Gd  26   0.40717
UO2Gd  25   0.40717
UO2Gd  24   0.40717
UO2Gd  23   0.40717
UO2Gd  22   0.40717
UO2Gd  21   0.40717
UO2Gd  20   0.20358
UO2Gd  19   0.20358
UO2Gd  18   0.20358
UO2Gd  17   0.20358
UO2Gd  16   0.20358
UO2Gd  15   0.20358
UO2Gd  14   0.20358
UO2Gd  13   0.20358
UO2Gd  12   0.20358
UO2Gd  11   0.20358
UO2Gd  10   0.20358
UO2Gd   9   0.20358
UO2Gd   8   0.20358
UO2Gd   7   0.20358
UO2Gd   6   0.20358
UO2Gd   5   0.20358
UO2Gd   4   0.20358
UO2Gd   3   0.20358
UO2Gd   2   0.20358
UO2Gd   1   0.20358