crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Ensure that unstable nuclides with no listed decay modes still decay

Browse source
This commit is contained in:
Paul Romano 2020-06-02 07:05:37 -05:00
parent ffe0f0283a
commit cff1262624
1 changed files with 4 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

10
openmc/deplete/chain.py
View file

@ -490,16 +490,14 @@ class Chain:
fission_yields = self.get_default_fission_yields()
for i, nuc in enumerate(self.nuclides):
if nuc.n_decay_modes != 0:
# Decay paths
# Loss
# Loss from radioactive decay
if nuc.half_life is not None:
decay_constant = math.log(2) / nuc.half_life
if decay_constant != 0.0:
matrix[i, i] -= decay_constant
# Gain
# Gain from radioactive decay
if nuc.n_decay_modes != 0:
for _, target, branching_ratio in nuc.decay_modes:
# Allow for total annihilation for debug purposes
if target is not None: