def test_dens_mol_phase_liq(m):
for j in m.params.component_list:
m.params.get_component(j).config.dens_mol_liq_comp = dummy_call
assert str(Ideal.dens_mol_phase(m.props[1], "Liq")) == str(
sum(m.props[1].mole_frac_phase_comp["Liq", j] * 42
for j in m.params.component_list))
def test_dens_mol_phase_sol(m_sol):
for j in m_sol.params.component_list:
m_sol.params.get_component(j).config.dens_mol_sol_comp = dummy_call
m_sol.params.get_component(j).config.dens_mol_liq_comp = dummy_call
assert str(Ideal.dens_mol_phase(m_sol.props[1], "Sol")) == str(
1 / sum(m_sol.props[1].mole_frac_phase_comp["Sol", j] / 42
for j in m_sol.params.component_list))
def test_dens_mol_phase_liq(m):
for j in m.params.component_list:
m.params.get_component(j).config.dens_mol_liq_comp = dummy_call
assert str(Ideal.dens_mol_phase(m.props[1], "Liq")) == str(
1 / (1 / 42 * m.props[1].mole_frac_phase_comp["Liq", "a"] +
1 / 42 * m.props[1].mole_frac_phase_comp["Liq", "b"] +
1 / 42 * m.props[1].mole_frac_phase_comp["Liq", "c"]))
def test_dens_mol_phase_invalid_phase(m_sol):
with pytest.raises(PropertyNotSupportedError):
Ideal.dens_mol_phase(m_sol.props[1], "Sol")
def test_dens_mol_phase_vap(m):
assert str(Ideal.dens_mol_phase(
m.props[1],
"Vap")) == ('props[1].pressure/(kg*m**2/J/s**2*(8.314462618*(J)/mol/K)'
'*props[1].temperature)')