Need help with this MCNPX Error

  • Thread starter Thread starter unique empire
  • Start date Start date
  • Tags Tags
    Mcnp5 Mcnpx
Click For Summary
SUMMARY

The forum discussion focuses on troubleshooting an MCNPX error related to cell definitions and geometry specifications in a provided MCNP code snippet. The user is advised to ensure that the volume for cell 12 is correctly defined, as it is currently set to infinite, which may lead to errors during simulation. The suggestion includes calculating the volume and updating the cell definition accordingly. Additionally, users are reminded to review the output file for detailed error messages that can provide further insights into the issue.

PREREQUISITES
  • Understanding of MCNPX (Monte Carlo N-Particle Transport Code) version 2.7 or later.
  • Familiarity with geometry definitions in MCNPX, including cell and surface specifications.
  • Knowledge of how to interpret command line messages and output files in MCNPX.
  • Basic skills in calculating volumes for geometric shapes used in simulations.
NEXT STEPS
  • Review the MCNPX documentation on cell definitions and geometry specifications.
  • Learn how to calculate volumes for different geometric configurations in MCNPX.
  • Explore common MCNPX error messages and their resolutions in the output files.
  • Practice writing and debugging MCNPX code snippets to enhance understanding of geometry and cell definitions.
USEFUL FOR

This discussion is beneficial for nuclear engineers, radiation physicists, and researchers using MCNPX for particle transport simulations who need to troubleshoot geometry-related errors in their models.

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:
Engineering 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.
 
  • Like
Likes   Reactions: unique empire and berkeman

Similar threads

  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 8 ·
Replies
8
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 7 ·
Replies
7
Views
3K
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K