How does MNCP calculate an F6 Tally?

  • Thread starter Thread starter MadGander
  • Start date Start date
  • Tags Tags
    Mcnp
AI Thread Summary
MCNP calculates an f6 tally by summing the energy deposited in a cell, which is then expressed in units of energy per cell volume. To compare an attenuated energy spectrum with an f6 tally, one approach is to calculate the average energy of the spectrum and divide it by the mass of the f6 tally cell, but discrepancies may arise due to scattered X-rays and electrons not reaching the detector. The average energy in the attenuated spectrum should not exceed that of the initial spectrum, indicating a potential error in calculations. It is essential to ensure that the total energy conservation is maintained across all energy bins. Considering an f8 tally may provide additional insights, and sharing MCNP input/output could help identify any mistakes in the analysis.
MadGander
Messages
28
Reaction score
1
TL;DR Summary
How does MNCP calculate an F6 Tally?
I'm wondering how exactly MCNP calculates an f6 tally? I'm trying to compare a theoretical result with an MCNP f6 tally (MeV/g). I have an initial energy spectrum and a thin layer of lead that attenuates the x-rays. Using the attenuation coefficient at each energy (bin width of 0.5 kev from 0-100 kev) and intensity at each energy, I produced an attenuated energy spectrum.

How do I compare this attenuated spectrum to an f6 tally? My thought process is to find the average energy within the spectrum and then divide by the mass of the f6 tally cell. However, this is a couple orders of magnitude off the MCNP ouput so I'm guessing it's wrong.

Another thing to note is that my average energy in the attenuated spectrum is higher than the average energy in an initial spectrum. The attenuator is a 0.08 cm thick lead sample, so there's no way the final average energy should be higher than the initial, right?
 
Last edited by a moderator:
Engineering news on Phys.org
Your simulation might need to account for scattered X-rays and electrons that don't make it to the detector. Orders of magnitude error makes me suspect a math mistake though.

Total Initial energy = Total Final energy + Total F6 + Total scattered (that does not hit the detector). Multiplying each energy bin by the midpoint of it's energy and summing the total should produce simple values that add up correctly. Individual bins don't have to obey conservation of energy but the total should.

I assume the F6 tally sums the energy loss in the cell and that value already has units of energy per cell volume (since twice the cell would tally twice the amount), that multiplied by the density would result in an answer that is in MeV/g.

All tallies are per source particle of course. Either multiply up the F6/scattered values or divide down the other two. If you have a largely monoenergetic source and you want it to look like a book value dividing down might make the most sense. Do each step one at a time and keep track of your units.

I'm wondering if an F8 tally would be better and I'm finding it hard to think clearly about this. If you are able to share your results and/or your MCNP input/output and/or your working out, that might give you the best chance of someone spotting the problem.
 
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...
Back
Top