crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 21:55:41 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Merge pull request #1679 from paulromano/decay-energy-fix

Browse source 
Fix calculation of decay energy for depletion chains
This commit is contained in:
Paul Romano 2020-10-07 13:36:07 -05:00 • committed by GitHub
commit 0d20134017
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
3 changed files with 14 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

14
openmc/data/decay.py
View file

@ -298,6 +298,8 @@ class Decay(EqualityMixin):
applications.
decay_constant : uncertainties.UFloat
Decay constant in inverse seconds.
decay_energy : uncertainties.UFloat
Average energy in [eV] per decay for decay heat applications
half_life : uncertainties.UFloat
Half-life of the decay in seconds.
modes : list
@ -347,7 +349,7 @@ class Decay(EqualityMixin):
items, values = get_list_record(file_obj)
self.half_life = ufloat(items[0], items[1])
NC = items[4]//2
pairs = [x for x in zip(values[::2], values[1::2])]
pairs = list(zip(values[::2], values[1::2]))
ex = self.average_energies
ex['light'] = ufloat(*pairs[0])
ex['electromagnetic'] = ufloat(*pairs[1])
@ -359,7 +361,7 @@ class Decay(EqualityMixin):
ex['conversion'] = ufloat(*pairs[6])
ex['gamma'] = ufloat(*pairs[7])
ex['xray'] = ufloat(*pairs[8])
ex['Bremsstrahlung'] = ufloat(*pairs[9])
ex['bremsstrahlung'] = ufloat(*pairs[9])
ex['annihilation'] = ufloat(*pairs[10])
ex['alpha'] = ufloat(*pairs[11])
ex['recoil'] = ufloat(*pairs[12])
@ -464,6 +466,14 @@ class Decay(EqualityMixin):
mu, sigma = self.half_life
return ufloat(log(2.)/mu, log(2.)/mu**2*sigma)
@property
def decay_energy(self):
energy = self.average_energies
if energy:
return energy['light'] + energy['electromagnetic'] + energy['heavy']
else:
return ufloat(0, 0)
@classmethod
def from_endf(cls, ev_or_filename):
"""Generate radioactive decay data from an ENDF evaluation

3
openmc/deplete/chain.py
View file

@ -395,8 +395,7 @@ class Chain:
if not data.nuclide['stable'] and data.half_life.nominal_value != 0.0:
nuclide.half_life = data.half_life.nominal_value
nuclide.decay_energy = sum(E.nominal_value for E in
data.average_energies.values())
nuclide.decay_energy = data.decay_energy.nominal_value
sum_br = 0.0
for i, mode in enumerate(data.modes):
type_ = ','.join(mode.modes)

2
tests/unit_tests/test_data_decay.py
View file

@ -34,7 +34,7 @@ def u235_yields():
def test_nb90_halflife(nb90):
ufloat_close(nb90.half_life, ufloat(52560.0, 180.0))
ufloat_close(nb90.decay_constant, log(2.)/nb90.half_life)
ufloat_close(nb90.decay_energy, ufloat(2265527.5, 25159.400474401213))
def test_nb90_nuclide(nb90):
assert nb90.nuclide['atomic_number'] == 41