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

updated depletion test to perform loading of new nuclides from the chain

Browse source
This commit is contained in:
Jingang Liang 2019-02-08 12:03:50 -08:00
parent 0adad96bf0
commit f8c8e4e972
3 changed files with 5 additions and 13 deletions
Download patch file Download diff file Expand all files Collapse all files

4
tests/regression_tests/deplete/example_geometry.py
View file

@ -47,10 +47,8 @@ def generate_initial_number_density():
fuel_dict['U238'] = 2.21371e22
fuel_dict['O16'] = 4.62954e22
fuel_dict['O17'] = 1.127684e20
fuel_dict['I135'] = 1.0e10
fuel_dict['Xe135'] = 1.0e10
fuel_dict['Xe136'] = 1.0e10
fuel_dict['Cs135'] = 1.0e10
fuel_dict['Gd156'] = 1.0e10
fuel_dict['Gd157'] = 1.0e10
# fuel_dict['O18'] = 9.51352e19 # Does not exist in ENDF71, merged into 17
@ -62,10 +60,8 @@ def generate_initial_number_density():
fuel_gd_dict['Gd156'] = 3.95517E+10
fuel_gd_dict['Gd157'] = 1.08156e20
fuel_gd_dict['O16'] = 4.64035e22
fuel_dict['I135'] = 1.0e10
fuel_dict['Xe136'] = 1.0e10
fuel_dict['Xe135'] = 1.0e10
fuel_dict['Cs135'] = 1.0e10
# There are a whole bunch of 1e-10 stuff here.
# Concentration to be used for cladding

BIN
tests/regression_tests/deplete/test_reference.h5
View file

Binary file not shown.

14
tests/unit_tests/test_deplete_resultslist.py
View file

@ -20,8 +20,7 @@ def test_get_atoms(res):
t, n = res.get_atoms("1", "Xe135")
t_ref = [0.0, 1296000.0, 2592000.0, 3888000.0]
n_ref = [6.674732823364922e+08, 4.022377919003244e+14,
3.469325382312551e+14, 3.639463751297369e+14]
n_ref = [6.67473282e+08, 3.76986925e+14, 3.68587383e+14, 3.91338675e+14]
np.testing.assert_allclose(t, t_ref)
np.testing.assert_allclose(n, n_ref)
@ -32,13 +31,11 @@ def test_get_reaction_rate(res):
t, r = res.get_reaction_rate("1", "Xe135", "(n,gamma)")
t_ref = [0.0, 1296000.0, 2592000.0, 3888000.0]
n_ref = [6.674732823364922e+08, 4.022377919003244e+14,
3.469325382312551e+14, 3.639463751297369e+14]
xs_ref = np.array([3.272630849911638e-05, 2.663333206501429e-05,
3.378999816491878e-05, 3.277013860196171e-05])
n_ref = [6.67473282e+08, 3.76986925e+14, 3.68587383e+14, 3.91338675e+14]
xs_ref = [3.32282266e-05, 2.76207120e-05, 4.10986677e-05, 3.72453665e-05]
np.testing.assert_allclose(t, t_ref)
np.testing.assert_allclose(r, n_ref * xs_ref)
np.testing.assert_allclose(r, np.array(n_ref) * xs_ref)
def test_get_eigenvalue(res):
@ -46,8 +43,7 @@ def test_get_eigenvalue(res):
t, k = res.get_eigenvalue()
t_ref = [0.0, 1296000.0, 2592000.0, 3888000.0]
k_ref = [1.2162584685809341, 1.1021840248590875, 1.223040353281542,
1.2360979537312164]
k_ref = [1.16984322, 1.19097427, 1.03012572, 1.20045627]
np.testing.assert_allclose(t, t_ref)
np.testing.assert_allclose(k, k_ref)