def check_NU(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.05
feed = Material(
{
922340000: 0.000055,
922350000: 0.00720,
922380000: 0.992745,
}
)
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1e-11)
# print "casc.mat_prod = " + repr(casc.mat_prod)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.05, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.0989473684211, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.901052631579, 1.0)
assert_almost_equal(casc.N / 27.183583424704818, 1.0, 4)
assert_almost_equal(casc.M / 13.387464890476533, 1.0, 4)
assert_almost_equal(casc.Mstar / 236.5621860655, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 288.62731727645644, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.761263453429, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 7.69362000806, 1.0, 5)
def check_vision(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.055
feed = Material(
{
922340000: 0.000183963025893197,
922350000: 0.00818576605617839,
922360000: 0.00610641667100979,
922380000: 0.985523854246919,
}
)
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1e-11)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.055, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.10830030583196934, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.89169969416803063, 1.0)
assert_almost_equal(casc.N / 27.38162850698868, 1.0, 2)
assert_almost_equal(casc.M / 15.09646512546496, 1.0, 2)
assert_almost_equal(casc.Mstar / 236.58177606549995, 1.0, 4)
assert_almost_equal(casc.l_t_per_feed / 326.8956175003255, 1.0, 4)
assert_almost_equal(casc.swu_per_feed / 0.85102089049, 1.0, 4)
assert_almost_equal(casc.swu_per_prod / 7.85797310499, 1.0, 4)
def check_vision(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.055
feed = Material({
922340000: 0.000183963025893197,
922350000: 0.00818576605617839,
922360000: 0.00610641667100979,
922380000: 0.985523854246919,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-11)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.055, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.10830030583196934, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.89169969416803063, 1.0)
assert_almost_equal(casc.N / 27.38162850698868, 1.0, 2)
assert_almost_equal(casc.M / 15.09646512546496, 1.0, 2)
assert_almost_equal(casc.Mstar / 236.58177606549995, 1.0, 4)
assert_almost_equal(casc.l_t_per_feed / 326.8956175003255, 1.0, 4)
assert_almost_equal(casc.swu_per_feed / 0.85102089049, 1.0, 4)
assert_almost_equal(casc.swu_per_prod / 7.85797310499, 1.0, 4)
def check_sample_feed(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.06
feed = Material(
{
922320000: 1.1 * (10.0**-9),
922340000: 0.00021,
922350000: 0.0092,
922360000: 0.0042,
922380000: 0.9863899989,
}
)
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1e-11, max_iter=100)
# print "casc.mat_prod = " + repr(casc.mat_prod)
# print "casc.mat_prod = " + repr(casc.mat_tail)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.06, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.11652173913043479, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.88347826086956527, 1.0)
assert_almost_equal(casc.N / 26.864660071132583, 1.0, 5)
assert_almost_equal(casc.M / 16.637884564470365, 1.0, 5)
assert_almost_equal(casc.Mstar / 236.57708506549994, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 357.3888391866117, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.9322804173594426, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 8.000914029577306, 1.0, 5)
def check_NU(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.05
feed = Material({
922340000: 0.000055,
922350000: 0.00720,
922380000: 0.992745,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-11)
#print "casc.mat_prod = " + repr(casc.mat_prod)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.05, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.0989473684211, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.901052631579, 1.0)
assert_almost_equal(casc.N / 27.183583424704818, 1.0, 4)
assert_almost_equal(casc.M / 13.387464890476533, 1.0, 4)
assert_almost_equal(casc.Mstar / 236.5621860655, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 288.62731727645644, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.761263453429, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 7.69362000806, 1.0, 5)
def check_sample_feed(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.06
feed = Material({
922320000: 1.1 * (10.0**-9),
922340000: 0.00021,
922350000: 0.0092,
922360000: 0.0042,
922380000: 0.9863899989,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-11, max_iter=100)
#print "casc.mat_prod = " + repr(casc.mat_prod)
#print "casc.mat_prod = " + repr(casc.mat_tail)
assert_almost_equal(casc.mat_prod.comp[922350000], 0.06, 5)
assert_almost_equal(casc.mat_tail.comp[922350000], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.11652173913043479, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.88347826086956527, 1.0)
assert_almost_equal(casc.N / 26.864660071132583, 1.0, 5)
assert_almost_equal(casc.M / 16.637884564470365, 1.0, 5)
assert_almost_equal(casc.Mstar / 236.57708506549994, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 357.3888391866117, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.9322804173594426, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 8.000914029577306, 1.0, 5)
def check_NU(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.05
feed = Material({
922340: 0.000055,
922350: 0.00720,
922380: 0.992745,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-11)
assert_almost_equal(casc.mat_prod.comp[922350], 0.05, 5)
assert_almost_equal(casc.mat_tail.comp[922350], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.0989473684211, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.901052631579, 1.0)
assert_almost_equal(casc.N / 27.1835088212, 1.0, 4)
assert_almost_equal(casc.M / 13.3875092512, 1.0, 4)
assert_almost_equal(casc.Mstar / 236.562179, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 288.627270162, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.761263453429, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 7.69362000806, 1.0, 5)
def check_sample_feed(solver):
orig_casc = enr.default_uranium_cascade()
orig_casc.x_prod_j = 0.06
feed = Material({
922320: 1.1 * (10.0**-9),
922340: 0.00021,
922350: 0.0092,
922360: 0.0042,
922380: 0.9863899989,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-11, max_iter=100)
assert_almost_equal(casc.mat_prod.comp[922350], 0.06, 5)
assert_almost_equal(casc.mat_tail.comp[922350], 0.0025, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.11652173913043479, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.88347826086956527, 1.0)
assert_almost_equal(casc.N / 26.8646352802, 1.0, 5)
assert_almost_equal(casc.M / 16.6379009423, 1.0, 5)
assert_almost_equal(casc.Mstar / 236.577085, 1.0, 5)
assert_almost_equal(casc.l_t_per_feed / 357.388791749, 1.0, 5)
assert_almost_equal(casc.swu_per_feed / 0.932280175218, 1.0, 5)
assert_almost_equal(casc.swu_per_prod / 8.0009119515, 1.0, 5)
def check_tungsten(solver):
# This test comes from 'Multicomponent Isotope Separation in Matched
# Abundance Ratio Cascades Composed of Stages with Large Separation Factors'
# by casc. von Halle, 1987.
orig_casc = enr.Cascade()
orig_casc.alpha = 1.16306
orig_casc.Mstar = 181.3
orig_casc.j = 741800
orig_casc.k = 741860
orig_casc.N = 30.0
orig_casc.M = 10.0
orig_casc.x_prod_j = 0.5109
orig_casc.x_tail_j = 0.00014
feed = Material(
{
741800000: 0.0014,
741820000: 0.26416,
741830000: 0.14409,
741840000: 0.30618,
741860000: 0.28417,
}
)
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1e-7)
assert_almost_equal(casc.mat_prod.comp[741800000], 0.5109, 5)
assert_almost_equal(casc.mat_tail.comp[741800000], 0.00014, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.0024669120526274574, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.99753308794737272, 1.0)
assert_almost_equal(casc.N / 43.557515688533513, 1.0, 2)
assert_almost_equal(casc.M / 11.49556481009056, 1.0, 2)
assert_almost_equal(casc.Mstar / 181.16425540249995, 1.0, 4)
assert_almost_equal(casc.l_t_per_feed / 96.81774564292206, 1.0, 3)
assert_almost_equal(casc.swu_per_feed / 2.22221945305, 1.0, 3)
assert_almost_equal(casc.swu_per_prod / 900.810164953, 1.0, 3)
def check_tungsten(solver):
# This test comes from 'Multicomponent Isotope Separation in Matched
# Abundance Ratio Cascades Composed of Stages with Large Separation Factors'
# by casc. von Halle, 1987.
orig_casc = enr.Cascade()
orig_casc.alpha = 1.16306
orig_casc.Mstar = 181.3
orig_casc.j = 741800
orig_casc.k = 741860
orig_casc.N = 30.0
orig_casc.M = 10.0
orig_casc.x_prod_j = 0.5109
orig_casc.x_tail_j = 0.00014
feed = Material({
741800000: 0.0014,
741820000: 0.26416,
741830000: 0.14409,
741840000: 0.30618,
741860000: 0.28417,
})
orig_casc.mat_feed = feed
casc = enr.multicomponent(orig_casc, solver=solver, tolerance=1E-7)
assert_almost_equal(casc.mat_prod.comp[741800000], 0.5109, 5)
assert_almost_equal(casc.mat_tail.comp[741800000], 0.00014, 5)
assert_almost_equal(casc.mat_feed.mass / 1.0, 1.0)
assert_almost_equal(casc.mat_prod.mass / 0.0024669120526274574, 1.0)
assert_almost_equal(casc.mat_tail.mass / 0.99753308794737272, 1.0)
assert_almost_equal(casc.N / 43.557515688533513, 1.0, 2)
assert_almost_equal(casc.M / 11.49556481009056, 1.0, 2)
assert_almost_equal(casc.Mstar / 181.16425540249995, 1.0, 4)
assert_almost_equal(casc.l_t_per_feed / 96.81774564292206, 1.0, 3)
assert_almost_equal(casc.swu_per_feed / 2.22221945305, 1.0, 3)
assert_almost_equal(casc.swu_per_prod / 900.810164953, 1.0, 3)