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

Expand activation test to cover all cases for FluxCollapseHelper

Browse source
This commit is contained in:
Paul Romano 2020-09-03 14:44:00 -05:00
parent 35f8f03f10
commit 7e48637c1d
1 changed files with 11 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

15
tests/unit_tests/test_deplete_activation.py
View file

@ -40,8 +40,16 @@ def model():
return model
@pytest.mark.parametrize("reaction_rate_mode", ["direct", "flux"])
def test_activation(run_in_tmpdir, model, reaction_rate_mode):
ENERGIES = np.logspace(log10(1e-5), log10(2e7), 100)
@pytest.mark.parametrize("reaction_rate_mode,reaction_rate_opts", [
("direct", {}),
("flux", {'energies': ENERGIES}),
("flux", {'energies': ENERGIES, 'reactions': ['(n,gamma)']}),
("flux", {'energies': ENERGIES, 'reactions': ['(n,gamma)'], 'nuclides': ['W186']}),
])
def test_activation(run_in_tmpdir, model, reaction_rate_mode, reaction_rate_opts):
# Determine (n.gamma) reaction rate using initial run
sp = model.run()
with openmc.StatePoint(sp) as sp:
@ -56,12 +64,11 @@ def test_activation(run_in_tmpdir, model, reaction_rate_mode):
chain.export_to_xml('test_chain.xml')
# Create transport operator
energies = np.logspace(log10(1e-5), log10(2e7), 100)
op = openmc.deplete.Operator(
model.geometry, model.settings, 'test_chain.xml',
normalization_mode="source-rate",
reaction_rate_mode=reaction_rate_mode,
reaction_rate_opts={'energies': energies},
reaction_rate_opts=reaction_rate_opts,
)
# To determine the source rate necessary to reduce W186 density in half, we