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))
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))
# 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])
# 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])