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

adjust heating number when reading from ACE

Browse source
This commit is contained in:
liangjg 2019-03-21 11:38:49 -04:00
parent 90c70f910a
commit 21888e6cb5
1 changed files with 11 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

16
openmc/data/photon.py
View file

@ -513,10 +513,6 @@ class IncidentPhoton(EqualityMixin):
for mt in (502, 504, 515, 522, 525):
data.reactions[mt] = PhotonReaction.from_ace(ace, mt)
# Get heating cross sections [eV*barn] from factors [MeV per collision]
data.reactions[525].xs.y *= sum([data.reactions[mt].xs.y for mt in
(502, 504, 515, 522)]) * EV_PER_MEV
# Compton profiles
n_shell = ace.nxs[5]
if n_shell != 0:
@ -997,7 +993,17 @@ class PhotonReaction(EqualityMixin):
# Store cross section
xs = ace.xss[idx : idx+n].copy()
if mt in (502, 504, 515, 522):
if mt == 525:
# Get heating cross sections [eV barn] from factors [MeV/collision]
# by multiplying the total xs (sum of (502, 504, 515, 522))
idx = ace.jxs[1] + n
total_xs = ace.xss[idx : idx + 4*n].copy()
total_xs = total_xs.reshape((4, n))
nonzero = (total_xs != 0.0)
total_xs[nonzero] = np.exp(total_xs[nonzero])
total_xs = total_xs.sum(axis=0)
xs *= total_xs * EV_PER_MEV
else:
nonzero = (xs != 0.0)
xs[nonzero] = np.exp(xs[nonzero])
rx.xs = Tabulated1D(energy, xs, [n], [5])