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