def test_enth_mol_phase(m): for j in m.params.component_list: m.params.get_component(j).config.enth_mol_liq_comp = dummy_call m.params.get_component(j).config.enth_mol_ig_comp = dummy_call m.params.get_component(j).config.dens_mol_liq_comp = dummy_call2 assert str(Ideal.enth_mol_phase( m.props[1], "Vap")) == ('props[1].mole_frac_phase_comp[Vap,a]*42.0 + ' 'props[1].mole_frac_phase_comp[Vap,b]*42.0 + ' 'props[1].mole_frac_phase_comp[Vap,c]*42.0') assert str(Ideal.enth_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) + (m.props[1].pressure - m.params.pressure_ref) / m.props[1].dens_mol_phase["Liq"])
def test_enth_mol_phase(m): m.props[1].enth_mol_phase_comp = Var(m.params.phase_list, m.params.component_list) for p in m.params.phase_list: assert str(Ideal.enth_mol_phase(m.props[1], p)) == str( sum(m.props[1].mole_frac_phase_comp[p, j] * m.props[1].enth_mol_phase_comp[p, j] for j in m.params.component_list))