Problem with photonuclear physics phys:p in MCNP

  • Thread starter Thread starter Nour_wahban
  • Start date Start date
  • Tags Tags
    Mcnp
AI Thread Summary
The discussion revolves around issues encountered while using MCNP to activate cobalt-59 with neutrons from californium-252. Users are receiving warnings indicating that photonuclear physics may be necessary and that material 3 is not utilized in the problem. It is suggested that photoneutron reactions from gamma emissions could be relevant, particularly since thresholds for such reactions typically exceed 7 MeV. The complexities of the PHYS:p card in MCNP are highlighted, with recommendations for settings that enable photonuclear physics without significantly impacting results. Additionally, it is clarified that the user should not be experiencing the material 3 error, as it is indeed being used in the simulation.
Nour_wahban
Messages
5
Reaction score
0
TL;DR Summary
I'm trying to activate cobalt 59 to cobalt 60 and obtain its peaks, but the phys:p is not working.
I'm working on a very basic setup to activate cobalt 59 with neutrons emitted from californium 252. MCNP is giving me the following warnings:
warning. photonuclear physics may be needed (phys:p).
warning. material 3 is not used in the problem.
could anyone help me understand where the problem is. I've attached the input file.
 

Attachments

Last edited:
Engineering news on Phys.org
Nour_wahban said:
TL;DR Summary: I'm trying to activate cobalt 59 to cobalt 60 and obtain its peaks, but the phys:p is not working.

warning. photonuclear physics may be needed (phys:p).
I suspect that one may have to consider photoneutron reactions from the gammas emitted by 252Cf. Most photoneutron thresholds are above 7 MeV.
See - https://journals.aps.org/prc/abstract/10.1103/PhysRevC.10.795

This is a useful reference: https://www.osti.gov/servlets/purl/1323135
MatMCNP: A Code for Producing Material Cards for MCNP

About photoneutrons
https://www.sciencedirect.com/topics/physics-and-astronomy/photoneutrons
A gamma photon with energy sufficiently large to overcome the neutron binding energy (about 7 MeV in most nuclides) may cause a (γ, n) reaction.
 
You will see a lot of warnings and you should chase them. Warnings often show up a bigger problem. That said, this one is usually a comment, not a warning. It can usually be ignored, unless you know it's a common event in the problem, but again I don't think I've seen it as an actual warning, so I don't know. Might be related to the use of a distribution function for your neutron source energy.

phys:p is a more complicated card that doesn't just turn on and off photonuclear physics, but a lot of other stuff. It also tends to change between versions so you have to check the right manual to make sure the switch is the one you want.
The default for 6.3, according to the manual is PHYS:p 100 0 0 0 0 J 0
I am guessing the safest way to turn on photonuclear physics would be,
PHYS:p 100 0 0 -1 0 J 0
That according to the manual is an unbiased 'sample every collision' photonuclear treatment.
This might slow things down and probably won't affect the results much.

I do not see an issue with material 3. In the def there are only 3 spaces for the extra lines, but & is used to continue the previous line, so that is absolutely fine. Material 3 is used in cell 3. You should not be getting that error with that input file as best as I can tell.

Hope that helps!

Edit, Just to check, you are after the fast activation spectrum (n,gamma) right? You are not going to see Co-60 decay, because the program doesn't handle time like that.
 
Last edited:
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