Need help with this MCNPX Error

  • Thread starter Thread starter unique empire
  • Start date Start date
  • Tags Tags
    Mcnp5 Mcnpx
AI Thread Summary
The discussion centers around a provided MCNP code, detailing various cell definitions and parameters for a simulation. Key points include the importance of accurately defining cell volumes and the implications of infinite versus limited geometries, specifically mentioning a cell that is initially infinite but later constrained to 2 meters. Participants emphasize the need to analyze output files thoroughly, not just command line messages, to ensure proper understanding and results. Additionally, there are references to material compositions and burn time settings, which are critical for the simulation's accuracy. Understanding these elements is essential for effective modeling in MCNP simulations.
unique empire
Messages
1
Reaction score
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
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 unique empire and berkeman
Hello, I'm currently trying to compare theoretical results with an MCNP simulation. I'm using two discrete sets of data, intensity (probability) and linear attenuation coefficient, both functions of energy, to produce an attenuated energy spectrum after x-rays have passed through a thin layer of lead. I've been running through the calculations and I'm getting a higher average attenuated energy (~74 keV) than initial average energy (~33 keV). My guess is I'm doing something wrong somewhere...

Similar threads

Back
Top