MCNP: Critical Mass of UO2 (10% Enrichment)

AI Thread Summary
A program was developed to calculate the critical mass of uranium oxide (UO2) with 10% enrichment, yielding a keff value near 1 for specific volume and density. The user inquired about automating the program to determine density and volume without manual input, noting that UO2 density varies between 10.4 and 10.9 g/cm^3. Suggestions included using a trick where the product of density and radius remains constant, allowing for adjustments without rerunning the program. The discussion emphasized the complexity of running kcode in MCNP, which primarily searches for k rather than solving for it directly. The conversation concluded with a focus on the importance of understanding the relationship between density, radius, and mass in critical mass calculations.
Tema3212
Messages
1
Reaction score
0
TL;DR Summary
Critical mass UO2 with an enrichment of 10%.
I wrote a program to determine the critical mass of uranium oxide with an enrichment of 10%. I got a keff approximately equal to 1 with the selected volume and density (attached a file). Is it possible to somehow run the program without writing the initial density and volume into the conditions? So that the program itself finds these values in the process of operation. After all, the practical density of UO2 varies from 10.4 to 10.9 gr/cm^3 and the volume and mass change accordingly
 

Attachments

  • 1637093757577.png
    1637093757577.png
    41 KB · Views: 212
  • 1637093806949.png
    1637093806949.png
    23.9 KB · Views: 183
  • inp111.txt
    inp111.txt
    929 bytes · Views: 172
Last edited:
Engineering news on Phys.org
:welcome:

On PF, we use only English. Please use only English in future posts.

translate.google.com translated your post as:

I wrote a program to determine the critical mass of uranium oxide with an enrichment of 10%. It turned out keff approximately equal to 1 for the selected volume and density (attached file). Is it possible to somehow run the program without writing the initial density and volume into the conditions? So that the program itself finds these values in the process. After all, the practical density of UO2 ranges from 10.4 to 10.9 g / cm ^ 3, and the volume and mass change accordingly.
 
  • Like
Likes Alex A, berkeman and Tema3212
I couldn't get the file to run without adding imp:n=0 to the void cell. I've simplified the geometry, removing an undeclared region and a surface.

When doing a kcode, the program isn't solving for k so much as it is searching for k. I don't know any way of getting mcnp to do what you want directly.

For such a simple problem there is a trick. Density * Radius appears to be an invariant. So if you run for Density=10.4 and Radius = 66.5 you have the same answer for all problems where Density * Radius = 691.6

So if the fuel has a density of 10.8 then 691.6/10.8 to gives you the radius without having to run again. Volume of a sphere * density to get the mass - this will change.

Code:
Critical mass uranium dioxide
c Cell Card                                                    
1 1  -10.4  -1 imp:n=1 
999 0        1 imp:n=0

c Surface Card
1 so  66.5

c Data Card
m1   92235 0.08815 92238 0.79333
     8016 0.11852                                                        
c Criticaly Control Card                                      
kcode 1000 1.0  50  600                                                
ksrc    0  0  0                                                        

print Table 50
 
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

Replies
10
Views
16K
Replies
2
Views
3K
Back
Top