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