The following issue shows up in burnup calculations with edepmode 3 activated, at some point in the depletion calculation (2.1.31 version).
It happens only when the burnable regions contain isotopes that are also present in non-burnable regions (here, as an example, Zr, present in both the fuel and the cladding).
It doesn't show up for the first transport calculation, but only later.
Attached is an input based on the 2D PWR assembly case, taken from one the examples, and the result:
Code: Select all
***** Fri Nov 1 15:31:34 2019: - MPI task = 0 - OpenMP thread = 3 - RNG parent seed = 1572643595011 - RNG history seed = 2506919101001949955 - RNG history idx = 5280000 Fatal error in function TrackMode: Total exceeds majorant (E = 5.109989E-01, f = 1.714415E-04, mat = fuelz36) Simulation aborted.
The case runs fine when edepmode is set to 0,1 or 2. So I am guessing it might be due to the photon transport part.
I tried disabling the delta tracking (set dt 0.) but the same issue showed up.
Using optimization mode 1 made it to work, but the calculation becomes really slow. Optimization modes 2,3,4 lead to the same failure.
Do you have any suggestions on how to get the case to run with a better optimization mode than 1?