Need help with this MCNPX Error

  • Thread starter Thread starter unique empire
  • Start date Start date
  • Tags Tags
    Mcnp5 Mcnpx
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 3K views
unique empire
Messages
1
Reaction score
0
TL;DR
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
1000011063.webp
 

Attachments

Last edited by a moderator:
Physics news on Phys.org
Hi @Unique empire, welcome to PhysicsForums.

Just a quick attempt. Cell 12 is inside cz 0.4369, and this is infinite. It's later limited to 2 meters long I think. Work out the volume and put it in the cell definition.

Code:
12 7 -10.17 -51 vol= u=4  imp:n=1

Don't just read the command line messages, read the output file.
 
Reply
  • Like
Likes   Reactions: unique empire and berkeman