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