2D PWR assembly burnup example (division with div-card)
Jump to navigation
Jump to 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