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 everyone, I am currently working on a burnup calculation for a fuel assembly with repeated geometric structures using MCNP6. I have defined two materials (Material 1 and Material 2) which are actually the same material but located in different positions. However, after running the calculation with the BURN card, I am encountering an issue where all burnup information(power fraction(Initial input is 1,but output file is 0), burnup, mass, etc.) for Material 2 is zero, while Material 1...
Hi everyone, I'm a complete beginner with MCNP and trying to learn how to perform burnup calculations. Right now, I'm feeling a bit lost and not sure where to start. I found the OECD-NEA Burnup Credit Calculational Criticality Benchmark (Phase I-B) and was wondering if anyone has worked through this specific benchmark using MCNP6? If so, would you be willing to share your MCNP input file for it? Seeing an actual working example would be incredibly helpful for my learning. I'd be really...

Similar threads

Back
Top