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

Update dictionary of secondary particles for Chain

Browse source
This commit is contained in:
Andrew Johnson 2019-08-23 11:43:19 -05:00
parent 6c25dbf33b
commit c8e205aff8
No known key found for this signature in database
GPG key ID: 253418E91B7F6FEB
1 changed files with 12 additions and 13 deletions
Download patch file Download diff file Expand all files Collapse all files

25
openmc/deplete/chain.py
View file

@ -104,17 +104,17 @@ def replace_missing(product, decay_data):
return product
_secondary_particles = {
"p": ["H1"], "d": ["H2"], "t": ["H3"], "3He": ["He3"], "a": ["He4"],
"2nd": ["H2"], "na": ["He4"], "n3a": ["He4"] * 3, "2na": ["He4"],
"3na": ["He4"], "np": ["H1"], "n2a": ["He4"] * 2,
"2n2a": ["He4"] * 2, "nd": ["H2"], "nt": ["H3"],
"nHe-3": ["He3"], "nd2a": ["H2", "He4"], "nt2a": ["H3", "He4", "He4"],
"2np": ["H1"], "3np": ["H1"], "n2p": ["H1"] * 2,
"2a": ["He4"] * 2, "3a": ["He4"] * 3, "2p": ["H1"] * 2,
"pa": ["H1", "He4"], "t2a": ["H3", "He4", "He4"],
"d2a": ["H2", "He4", "He4"], "pd": ["H1", "H2"], "pt": ["H1", "H3"],
"da": ["H2", "He4"]}
_SECONDARY_PARTICLES = {
"(n,p)": ["H1"], "(n,d)": ["H2"], "(n,t)": ["H3"], "(n,3He)": ["He3"],
"(n,a)": ["He4"], "(n,2nd)": ["H2"], "(n,na)": ["He4"], "(n,3na)": ["He4"],
"(n,n3a)": ["He4"] * 3, "(n,2na)": ["He4"], "(n,np)": ["H1"],
"(n,n2a)": ["He4"] * 2, "(n,2n2a)": ["He4"] * 2, "(n,nd)": ["H2"],
"(n,nt)": ["H3"], "(n,nHe-3)": ["He3"], "(n,nd2a)": ["H2", "He4"],
"(n,nt2a)": ["H3", "He4", "He4"], "(n,2np)": ["H1"], "(n,3np)": ["H1"],
"(n,n2p)": ["H1"] * 2, "(n,2a)": ["He4"] * 2, "(n,3a)": ["He4"] * 3,
"(n,2p)": ["H1"] * 2, "(n,pa)": ["H1", "He4"],
"(n,t2a)": ["H3", "He4", "He4"], "(n,d2a)": ["H2", "He4", "He4"],
"(n,pd)": ["H1", "H2"], "(n,pt)": ["H1", "H3"], "(n,da)": ["H2", "He4"]}
class Chain(object):
@ -559,8 +559,7 @@ class Chain(object):
bad_sums = {}
# Secondary products, like alpha particles, should not be modified
secondary = _secondary_particles.get(
reaction[reaction.index(",") + 1:-1], [])
secondary = _SECONDARY_PARTICLES.get(reaction, [])
# Check for validity before manipulation