def test_wakao_model():
    mat = Material(k=1 * units.watt / units.meter / units.kelvin,
                   cp=1 * units.joule / units.kg / units.kelvin,
                   mu=2 * units.pascal * units.second)
    h_wakao = ConvectiveModel(mat=mat,
                              m_flow=1 * units.kg / units.g,
                              a_flow=1 * units.meter**2,
                              length_scale=1 * units.meter,
                              model='wakao')
    assert_equal(h_wakao.mu, 2 * units.pascal * units.second)
    rho = 100 * units.kg / units.meter**3
    assert_equal(h_wakao.h(rho, 0 * units.pascal * units.second),
                 h_wakao.h(rho, 2 * units.pascal * units.second))
Beispiel #2
0
def test_wakao_model():
    mat=Material(k=1*units.watt/units.meter/units.kelvin,
                 cp=1*units.joule/units.kg/units.kelvin,
                 mu=2*units.pascal*units.second
                 )
    h_wakao = ConvectiveModel(mat=mat,
                              m_flow=1*units.kg/units.g,
                              a_flow=1*units.meter**2,
                              length_scale=1*units.meter,
                              model='wakao')
    assert_equal(h_wakao.mu,
                 2*units.pascal*units.second)
    rho = 100*units.kg/units.meter**3
    assert_equal(h_wakao.h(rho, 0*units.pascal*units.second),
                 h_wakao.h(rho, 2*units.pascal*units.second))
Beispiel #3
0
# mesh size for the fuel pebble FVM calculation
l = 0.0005 * units.meter
comp_list = mod.mesh(l)
comp_list.extend(fuel.mesh(l))
comp_list.extend(shell.mesh(l))
pebble = th.THSuperComponent('pebble', t_shell, comp_list, timer=ti)
# Add convective boundary condition to the pebble
pebble.add_conv_bc('cool', h=h_cool)

cool = th.THComponent(name="cool",
                      mat=cool,
                      vol=vol_cool,
                      T0=t_cool,
                      alpha_temp=alpha_cool,
                      timer=ti)
# The coolant convects to the shell
cool.add_convection('pebble', h=h_cool, area=a_pb)
cool.add_advection('cool', m_flow / n_pebbles, t_inlet, cp=cool.cp)

components = []
for i in range(0, len(pebble.sub_comp)):
    components.append(pebble.sub_comp[i])
components.extend([pebble, cool])

uncert = [
    alpha_cool, alpha_fuel, k_mod, k_fuel, k_shell, cp_mod, cp_fuel, cp_shell,
    cp_cool,
    h_cool.h()
]
uncertainty_param = np.array([o.magnitude for o in uncert])
Beispiel #4
0
# Add convective boundary condition to the pebble
pebble.add_conv_bc('cool', h=h_cool)

cool = th.THComponent(name="cool",
                      mat=cool,
                      vol=vol_cool,
                      T0=t_cool,
                      alpha_temp=alpha_cool,
                      timer=ti)
# The coolant convects to the shell
cool.add_convection('pebble', h=h_cool, area=a_pb)
cool.add_advection('cool', m_flow/n_pebbles, t_inlet, cp=cool.cp)

components = []
for i in range(0, len(pebble.sub_comp)):
    components.append(pebble.sub_comp[i])
components.extend([pebble, cool])

uncert = [
    alpha_cool,
    alpha_fuel,
    k_mod,
    k_fuel,
    k_shell,
    cp_mod,
    cp_fuel,
    cp_shell,
    cp_cool,
    h_cool.h()]
uncertainty_param = np.array([o.magnitude for o in uncert])