Fatal error in function ProcessPhotonProd

Compiler issues (Intel compiler, mpicc, etc.), data libraries and everything needed to get Serpent running is discussed here
Post Reply
gtdrew
Posts: 2
Joined: Tue Jun 05, 2018 9:46 pm
Security question 1: No
Security question 2: 93

Fatal error in function ProcessPhotonProd

Post by gtdrew » Tue Jul 09, 2019 12:10 am

Dear all,

I am attempting to get coupled neutron-photon transport simulations working in Serpent, and I am encountering the following error: "Fatal error in function ProcessPhotonProd, No fission channels found, Simulation aborted". In addition, I get the warning: "Warning message from function ProcessPhotonProd: CDF(NP - 1) = 0.000000E+00". I have the MCNP photon data libraries in my acelib path, and I have the auxiliary photon data in my pdatadir path. I have attempted runs with both the 2.1.30 version (with source code changes from viewtopic.php?f=25&t=2839&p=8402#p8402) and 2.1.31 version, both yielding the same error. I feel like the issue is stemming from an issue with my pdatadir files, but if anyone could provide some wisdom regarding this issue and possible solutions, it would be much appreciated.

User avatar
Jaakko Leppänen
Site Admin
Posts: 2167
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: Fatal error in function ProcessPhotonProd

Post by Jaakko Leppänen » Tue Jul 09, 2019 10:39 am

The subroutine actually handles the photon production reactions in neutron data. Which neutron data libraries are you using?
- Jaakko

MohammedQah
Posts: 10
Joined: Wed Mar 29, 2017 2:06 am
Security question 1: No
Security question 2: 64

Re: Fatal error in function ProcessPhotonProd

Post by MohammedQah » Wed Aug 28, 2019 7:50 pm

I got this error when I use KERMA data:

Fatal error in function ReadPendfData:

Mismatch in lowest energy of ACE and KERMA datas for nuclide 90232.03c (mt=3318)
ACE: 9.99999999999999939497E-12 KERMA: 3.99441299999999960973E-03


and also many warning messages of negative xs points or negative non-fission KERMAs in data.
any suggestion on that error.

Riku Tuominen
Posts: 36
Joined: Tue Feb 27, 2018 10:52 am
Security question 1: No
Security question 2: 93

Re: Fatal error in function ProcessPhotonProd

Post by Riku Tuominen » Thu Aug 29, 2019 2:00 pm

This error has nothing to do with processphotonprod, so next time please use the correct section of the forum and make a new thread if necessary. In order to fix the error remove the following lines of code from readpendfdata.c:

lines 419-422:
else
Die(FUNCTION_NAME, "Mismatch in lowest energy of ACE and KERMA datas "
"for nuclide %s (mt=%ld) \nACE: %.20E KERMA: %.20E",
GetText(nuc + NUCLIDE_PTR_NAME), mt, E[0], E0[0]);

lines 426-429:
else
Die(FUNCTION_NAME, "Mismatch in highest energy of ACE and KERMA datas "
"for nuclide %s (mt=%ld) \nACE: %.20E KERMA: %.20E",
GetText(nuc + NUCLIDE_PTR_NAME), mt, E[NES - 1], E0[np - 1]);

The issue of negative kerma coefficients is briefly discussed in section 2.2 of the energy deposition paper.

Post Reply