Example #1
0
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)
Example #2
0
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)
Example #3
0
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)
Example #4
0
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)
Example #5
0
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)
Example #6
0
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)
Example #7
0
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)
Example #8
0
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)
Example #9
0
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)
Example #10
0
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)