1import os
2import sys
3
4import numpy as np
5
6from openmc import Material
7from openmc.data import NATURAL_ABUNDANCE, atomic_mass
8
9
10def test_enrichment():
11    # This test doesn't require an OpenMC run.  We just need to make sure the
12    # element.expand() method expands Uranium to the proper enrichment.
13
14    uranium = Material()
15    uranium.add_element('U', 1.0, 'wo', 4.95)
16    densities = uranium.get_nuclide_densities()
17
18    sum_densities = 0.
19    for nuc in densities.keys():
20        assert nuc in ('U234', 'U235', 'U236', 'U238')
21        sum_densities += densities[nuc][1]
22
23    # Compute the weight percent U235
24    enrichment = densities['U235'][1] / sum_densities
25    assert np.isclose(enrichment, 0.0495, rtol=1.e-8)
26
27    # Compute the ratio of U234/U235
28    u234_to_u235 = densities['U234'][1] / densities['U235'][1]
29    assert np.isclose(u234_to_u235, 0.0089, rtol=1.e-8)
30
31    # Compute the ratio of U236/U235
32    u236_to_u235 = densities['U236'][1] / densities['U235'][1]
33    assert np.isclose(u236_to_u235, 0.0046, rtol=1.e-8)
34