Why is My MCNP Burnup Calculation Failing?

  • Thread starter Thread starter Rafimah
  • Start date Start date
  • Tags Tags
    Mcnp Simulation
AI Thread Summary
The user is encountering a geometry error in their MCNP burnup calculation, specifically receiving a "no cell found" message, which results in lost particles and no final keff estimate. The issue likely stems from the void specification, where the user needs to ensure that the geometry correctly encompasses all necessary space around the pellet. It is suggested to use the union operator for the void cell specification to address this problem. Additionally, for repeated geometry, placing the pellet in a reflecting box is recommended. Clarification on these geometry specifications is crucial for successful calculations in MCNP.
Rafimah
Messages
13
Reaction score
1
TL;DR Summary
I'm trying to learn to calculate burnup in MCNP but I'm having some difficulties
Hi everyone, I'm a newbie to MCNP, I'm trying to calculate burnup for this pellet I include here in a PWR in an infinitely repeated geometry, but it seems to be failing for some reason. I get the error message:

ctm = 0.00 nrn = 0
dump 1 on file runtpp.h5 nps = 0 coll = 0
source distribution written to file srcte cycle= 0
xact is done
cp0 = 0.18
Geometry error: no cell found.
run terminated because 10 particles got lost.
=====> 4.05 M histories/hr (based on wall-clock time in mcrun)
warning. 10 particles got lost
warning. no final keff estimate made because no active cycles run.
ctm = 0.00 nrn = 189
dump 2 on file runtpp.h5 nps = 10 coll = 0
mcrun is done

Can anyone help me out with this? I imagine it has to do with the "geometry error: no cell found" line but I don't understand the issue. The cell seems alright to me and the ksrc point is inside the cylinder.

Any help would be very appreciated!
 

Attachments

Engineering news on Phys.org
"20 0 1 2 -3"
This is the void, so it needs to be all space you don't care about. Everything above AS WELL AS everything below AS WELL AS everything outside the cylinder. So you need the union operator :
Code:
1:2:-3
As a rule of thumb, expect to use the union operator a lot in the void cell specification. I've not checked anything else and this is not yet a repeated geometry, the simplest way to do that would be to put the pellet in a reflecting box.

I was going to welcome you to PF, but unusually for an MCNP question you have other posts. Welcome to Nuclear Engineering anyway @Rafimah !
 
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

Replies
8
Views
1K
Replies
13
Views
3K
Replies
1
Views
4K
Replies
20
Views
18K
Replies
4
Views
3K
Back
Top