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