def model(self): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.params = CubicParameterBlock(default={"valid_phase": "Vap"}) m.fs.params.a = Component() m.fs.params.b = Component() m.fs.params.cubic_type = CubicEoS.SRK m.fs.params.gas_const = Param(default=8.314462618) m.fs.params.pressure_crit = Param( m.fs.params.component_list, initialize={'a': 5e6, 'b': 4e6}) m.fs.params.temperature_crit = Param( m.fs.params.component_list, initialize={"a": 500, "b": 600}) m.fs.params.omega = Param( m.fs.params.component_list, initialize={"a": 0.2, "b": 0.2}) m.fs.params.kappa = Param( m.fs.params.component_list, m.fs.params.component_list, initialize={('a', 'a'): 0.0, ('a', 'b'): 0.0, ('b', 'a'): 0.0, ('b', 'b'): 0.0}) return m
def test_build_default(self): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.params = CubicParameterBlock() assert m.fs.params.state_block_class is CubicStateBlock assert m.fs.params.config.valid_phase == ('Vap', 'Liq') assert isinstance(m.fs.params.phase_list, Set) assert len(m.fs.params.phase_list) == 2 for p in m.fs.params.phase_list: assert p in ["Vap", "Liq"]