I have tried to benchmark the photon transport implemented in OpenMC against MCNP and Tripoli4, with JEFF33 data. So I modeled two spheres, the inner most is in U235 (density of 19.5g/cm3) and the outer most is made of air (N14 at 80% and O16 at 20%, with a density of 0.001g/cm3). Their respective radius are 8,5cm and 15cm. Both are centered on the same point. I also placed a cylinder in the sphere of air, at a distance of the center of 10.5cm. The cylinder has a radius of 2cm and a height of 5cm. The material it is made is Si28 at 50% and C0 (C12 in MCNP) at 50% and has a density of 3.15g/cm3. All materials are at 293.15K.
I made a criticality calculation and measured the heat deposited by neutrons and photons on the cylinder. I also measured the flux discretized in 18 energy groups of neutrons and photons in the different parts/cells. I found significant differences for the heat deposition caused by photons. The figures below illustrates that issue. MCNP and Tripoli4 seems to agree on the same value while OpenMC does not. However, the mean squared error for the comparison of Tripoli4 and OpenMC against MCNP shows that the photon transport is correct since they all agree.
Energy levels in MeV are : [
The flux spectrum shown above is for the flux in the relevant cell. So yes, it does agree. In the github issue I added the output file of MCNP (“crit-ref” I could not attach it here because of the file extension) where you will find all the input data for the simulation.
Thank you for pointing out the mistake in MCNP. This will be corrected. I will repost the corrected result. I am not familiar with the delayed_photon_scaling, I will look into that.
Here are the xml’s. settings.xml (435 Bytes) materials.xml (587 Bytes) geometry.xml (543 Bytes)
Just ran this;
I’m 99.9% sure that it’s the delayed gammas.
I get a result 40% lower with delayed gammas off; in-line with your MCNP/Tripoli result.
By default MCNP does not simulate delayed gammas, only prompt gammas.
Depending on the nuclear system in question delayed gammas can be ~50% of the total gamma energy produced (or compelled) from fission.
There are delayed gamma modules in MCNP (act gam, something like that). And other means to simulate delayed gammas (R2S, etc.)
These delayed gammas can also sometimes be called capture gammas or other.