unique empire
- 1
- 0
- TL;DR Summary
- I need help,please!!!
I tried to write the MCNP code for a South Korean SMART reactor and wanted to do a burnup analysis of it via MCNPX, but something went wrong
This is my mcnp code:
Code:
1 4 -0.7119 1 -2 3 -4 5 -6 fill=1 imp:n=1
26 4 -0.7119 11 -12 13 -14 5 -6 u=12 lat=1 fill=13 imp:n=1
20 4 -0.7119 46 -47 -45 #1 fill=12 imp:n=1
19 0 #20 imp:n=0
3 4 -0.7119 11 -12 13 -14 u=1 lat=1 imp:n=1 fill=-4:4 -4:4 0:0
12 2r 8 9 8 12 2r
12 12 9 10 11 10 9 12 12
12 9 10 7 6 7 10 9 12
8 10 7 6 7 6 7 10 8
9 11 6 7 6 7 6 11 9
8 10 7 6 7 6 7 10 8
12 9 10 7 6 7 10 9 12
12 12 9 10 11 10 9 12 12
12 2r 8 9 8 12 2r
18 4 -0.7119 41 -42 43 -44 u=6 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 33r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 2r 4 2 8r 4 2 2r
2 7r 3 2 7r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 33r
2 2 4 2 3 4 2 2 5 2 2 4 3 2 4 2 2
2 33r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 7r 3 2 7r
2 2r 4 2 8r 4 2 2r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 33r
21 4 -0.7119 41 -42 43 -44 u=7 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 33r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 2r 4 2 8r 4 2 2r
2 5r 3 2 2 2 3 2 5r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 3r 3 2 6r 3 2 3r
2 16r
2 2 4 2 2 4 2 2 5 2 2 4 2 2 4 2 2
2 16r
2 3r 3 2 6r 3 2 3r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 5r 3 2 2 2 3 2 5r
2 2r 4 2 8r 4 2 2r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 33r
22 4 -0.7119 41 -42 43 -44 u=8 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 33r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 2r 4 2 8r 4 2 2r
2 16r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 33r
2 2 4 2 2 4 2 2 5 2 2 4 2 2 4 2 2
2 33r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 16r
2 2r 4 2 8r 4 2 2r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 33r
23 4 -0.7119 41 -42 43 -44 u=9 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 33r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 2r 4 2 8r 4 2 2r
2 7r 3 2 7r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 33r
2 2 4 2 3 4 2 2 5 2 2 4 3 2 4 2 2
2 33r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 7r 3 2 7r
2 2r 4 2 8r 4 2 2r
2 2 3 2 2 4 2 2 4 2 2 4 2 2 3 2 2
2 33r
24 4 -0.7119 41 -42 43 -44 u=10 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 16r
2 3r 3 2 6r 3 2 3r
2 4r 4 2 2 4 2 2 4 2 4r
2 2r 4 2 2 3 2 2r 3 2 2 4 2 2r
2 3 2 12r 3 2
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 2r 3 2 2 3 2 2 2 3 2 2 3 2 2r
2 16r
2 2 4 2 2 4 2 2 5 2 2 4 2 2 4 2 2
2 16r
2 2r 3 2 2 3 2 2 2 3 2 2 3 2 2r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 3 2 12r 3 2
2 2r 4 2 2 3 2 2r 3 2 2 4 2 2r
2 4r 4 2 2 4 2 2 4 2 4r
2 3r 3 2 6r 3 2 3r
2 16r
25 4 -0.7119 41 -42 43 -44 u=11 lat=1 imp:n=1 fill=-8:8 -8:8 0:0
2 16r
2 2r 3 2 8r 3 2 2r
2 4r 4 2 2 4 2 2 4 2 4r
2 3 2 4 2 2r 3 2 3 2 2r 4 2 3 2
2 3r 3 2 6r 3 2 3r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 5r 3 2 2 2 3 2 5r
2 2r 3 2 8r 3 2 2r
2 2 4 2 2 4 2 2 5 2 2 4 2 2 4 2 2
2 2r 3 2 8r 3 2 2r
2 5r 3 2 2 2 3 2 5r
2 2 4 2 2 4 2 2 4 2 2 4 2 2 4 2 2
2 3r 3 2 6r 3 2 3r
2 3 2 4 2 2r 3 2 3 2 2r 4 2 3 2
2 4r 4 2 2 4 2 2 4 2 4r
2 2r 3 2 8r 3 2 2r
2 16r
4 1 -10.286 -21 u=2 imp:n=1
5 5 -0.001625 21 -22 u=2 imp:n=1
6 3 -6.56 22 -23 u=2 imp:n=1
7 4 -0.7119 23 u=2 imp:n=1
8 2 -10.017 -21 u=3 imp:n=1
9 5 -0.001625 21 -22 u=3 imp:n=1
10 3 -6.56 22 -23 u=3 imp:n=1
11 4 -0.7119 23 u=3 imp:n=1
12 7 -10.17 -51 u=4 imp:n=1
13 6 -7.9 51 -52 u=4 imp:n=1
14 4 -0.7119 52 -31 u=4 imp:n=1
30 3 -6.56 31 -32 u=4 imp:n=1
31 4 -0.7119 32 u=4 imp:n=1
15 4 -0.7119 -31 u=5 imp:n=1
16 3 -6.56 31 -32 u=5 imp:n=1
17 4 -0.7119 32 u=5 imp:n=1
27 4 -0.7119 -48 u=13 imp:n=1
28 3 -6.56 48 -49 u=13 imp:n=1
29 4 -0.7119 49 u=13 imp:n=1
1 px -96.3747
2 px 96.3747
3 py -96.3747
4 py 96.3747
5 pz -100
6 pz 100
11 px -10.7083
12 px 10.7083
13 py -10.7083
14 py 10.7083
21 cz 0.4096
22 cz 0.41875
23 cz 0.475
31 cz 0.5615
32 cz 0.612
41 px -0.6299
42 px 0.6299
43 py -0.6299
44 py 0.6299
45 cz 136.2944
46 pz -120
47 pz 120
48 cz 0.0001
49 cz 0.0002
51 cz 0.4369
52 cz 0.48385
ksrc 1 1 0
BURN TIME=3
365 365 365
MAT=1
POWER=330
PFRAC=1 1 1
AFMIN=1e-12
OMIT=54135
MATVOL=1514593.481
BOPT=1 -11
m1 92235 -0.03 92238 -0.85 8016 -0.12
m2 92235 -0.04 92238 -0.77 8016 -0.12 64160 -0.07 $IFBA rod
m3 40091 -0.982 50119 -0.015 26056 -0.002 24052 -0.001
m4 1001 2 8016 1
m5 2004 1
m6 26056 -0.7 24052 -0.18 28059 -0.08 6012 -0.02
14028 -0.01 15031 -0.02 16032 -0.01
m7 47107 -0.8 49115 -0.15 48112 -0.05 $Ag-In-Cd
print
Attachments
Last edited by a moderator: