def test_sigma_a():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
sig_a = sigma_a("U238", group_struct=E_g)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a("U238", group_struct=E_g)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a("U235")
observed = (0.0 <= sig_a).all()
assert_true(observed)
def test_sigma_a():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
sig_a = sigma_a('U238', group_struct=E_g)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a('U238', group_struct=E_g)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a('U235')
observed = (0.0 <= sig_a).all()
assert_true(observed)
def test_sigma_a():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
E_n = xs_cache['E_n']
phi_n = np.ones(len(E_n) - 1)
sig_a = sigma_a('U238', E_g, E_n, phi_n)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a('U238', E_g, E_n, phi_n)
observed = (0.0 <= sig_a).all()
assert_true(observed)
sig_a = sigma_a('U235')
observed = (0.0 <= sig_a).all()
assert_true(observed)
def test_sigma_t():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
sig_t = sigma_t("U238", 600.0, E_g)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a("U238", 600.0, E_g) + sigma_s("U238", 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t("U238", 600.0, E_g)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a("U238", 600.0, E_g) + sigma_s("U238", 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t("U235")
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a("U235", 600.0) + sigma_s("U235", 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
def test_sigma_t():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
sig_t = sigma_t('U238', 600.0, E_g)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U238', 600.0, E_g) + sigma_s('U238', 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t('U238', 600.0, E_g)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U238', 600.0, E_g) + sigma_s('U238', 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t('U235')
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U235', 600.0) + sigma_s('U235', 600.0, E_g)
assert_array_almost_equal(sig_t, expected)
def test_sigma_t():
E_g = np.array([10.0, 7.5, 5.0, 2.5, 0.1])
E_n = xs_cache['E_n']
phi_n = np.ones(len(E_n) - 1)
sig_t = sigma_t('U238', 600.0, E_g, E_n, phi_n)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U238', E_g, E_n, phi_n) + sigma_s('U238', 600.0, E_g, E_n, phi_n)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t('U238', 600.0, E_g, E_n, phi_n)
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U238', E_g, E_n, phi_n) + sigma_s('U238', 600.0, E_g, E_n, phi_n)
assert_array_almost_equal(sig_t, expected)
sig_t = sigma_t('U235')
observed = (0.0 <= sig_t).all()
assert_true(observed)
expected = sigma_a('U235') + sigma_s('U235', 600.0, E_g, E_n, phi_n)
assert_array_almost_equal(sig_t, expected)
def _get_destruction(self, nuc, decay=True):
"""Computes the destruction rate of the nuclide.
Parameters
----------
nuc : int
Name of the nuclide in question
decay : bool
True if the decay constant should be added to the returned value.
False if only destruction from neutron reactions should be considered.
Returns
-------
d : float
Destruction rate of the nuclide.
"""
xscache = self.xscache
sig_a = sigma_a(nuc, xs_cache=xscache)
d = utils.from_barns(sig_a[0], 'cm2') * xscache['phi_g'][0]
if decay and not np.isnan(data.decay_const(nuc)):
d += data.decay_const(nuc)
return d
def _get_destruction(self, nuc, decay=True):
"""Computes the destruction rate of the nuclide.
Parameters
----------
nuc : int
Name of the nuclide in question
decay : bool
True if the decay constant should be added to the returned value.
False if only destruction from neutron reactions should be considered.
Returns
-------
d : float
Destruction rate of the nuclide.
"""
xscache = self.xscache
sig_a = sigma_a(nuc, xs_cache=xscache)
d = utils.from_barns(sig_a[0], 'cm2') * xscache['phi_g'][0]
if decay and not np.isnan(data.decay_const(nuc)):
d += data.decay_const(nuc)
return d
