def test_build(): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.properties = SaponificationParameterBlock() m.fs.reactions = SaponificationReactionParameterBlock( default={"property_package": m.fs.properties}) m.fs.pfr = PFR( default={ "property_package": m.fs.properties, "reaction_package": m.fs.reactions, "has_equilibrium_reactions": False, "has_heat_transfer": True, "has_pressure_change": True }) assert hasattr(m.fs.pfr, "inlet") assert len(m.fs.pfr.inlet.vars) == 4 assert hasattr(m.fs.pfr.inlet, "flow_vol") assert hasattr(m.fs.pfr.inlet, "conc_mol_comp") assert hasattr(m.fs.pfr.inlet, "temperature") assert hasattr(m.fs.pfr.inlet, "pressure") assert hasattr(m.fs.pfr, "outlet") assert len(m.fs.pfr.outlet.vars) == 4 assert hasattr(m.fs.pfr.outlet, "flow_vol") assert hasattr(m.fs.pfr.outlet, "conc_mol_comp") assert hasattr(m.fs.pfr.outlet, "temperature") assert hasattr(m.fs.pfr.outlet, "pressure") assert hasattr(m.fs.pfr, "performance_eqn") assert hasattr(m.fs.pfr.control_volume, "heat") assert hasattr(m.fs.pfr, "heat_duty") assert hasattr(m.fs.pfr, "deltaP")
def test_config(): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.properties = PhysicalParameterTestBlock() m.fs.reactions = ReactionParameterTestBlock( default={"property_package": m.fs.properties}) m.fs.unit = PFR(default={ "property_package": m.fs.properties, "reaction_package": m.fs.reactions }) # Check unit config arguments assert len(m.fs.unit.config) == 19 assert m.fs.unit.config.material_balance_type == \ MaterialBalanceType.componentPhase assert m.fs.unit.config.energy_balance_type == \ EnergyBalanceType.enthalpyTotal assert m.fs.unit.config.momentum_balance_type == \ MomentumBalanceType.pressureTotal assert not m.fs.unit.config.has_heat_transfer assert not m.fs.unit.config.has_pressure_change assert not m.fs.unit.config.has_equilibrium_reactions assert not m.fs.unit.config.has_phase_equilibrium assert not m.fs.unit.config.has_heat_of_reaction assert m.fs.unit.config.property_package is m.fs.properties assert m.fs.unit.config.reaction_package is m.fs.reactions assert m.fs.unit.config.length_domain_set == [0.0, 1.0] assert m.fs.unit.config.transformation_method == "dae.finite_difference" assert m.fs.unit.config.transformation_scheme == "BACKWARD" assert m.fs.unit.config.finite_elements == 20 assert m.fs.unit.config.collocation_points == 3
def sapon(self): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.properties = SaponificationParameterBlock() m.fs.reactions = SaponificationReactionParameterBlock(default={ "property_package": m.fs.properties}) m.fs.unit = PFR(default={"property_package": m.fs.properties, "reaction_package": m.fs.reactions, "has_equilibrium_reactions": False, "has_heat_transfer": True, "has_heat_of_reaction": True, "has_pressure_change": True}) return m
def test_report(): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.properties = SaponificationParameterBlock() m.fs.reactions = SaponificationReactionParameterBlock( default={"property_package": m.fs.properties}) m.fs.pfr = PFR( default={ "property_package": m.fs.properties, "reaction_package": m.fs.reactions, "has_equilibrium_reactions": False, "has_heat_transfer": False, "has_pressure_change": False }) m.fs.pfr.report()
def test_initialize(): m = ConcreteModel() m.fs = FlowsheetBlock(default={"dynamic": False}) m.fs.properties = SaponificationParameterBlock() m.fs.reactions = SaponificationReactionParameterBlock( default={"property_package": m.fs.properties}) m.fs.pfr = PFR( default={ "property_package": m.fs.properties, "reaction_package": m.fs.reactions, "has_equilibrium_reactions": False, "has_heat_transfer": False, "has_pressure_change": False }) m.fs.pfr.inlet.flow_vol.fix(1.0) m.fs.pfr.inlet.conc_mol_comp[0, "H2O"].fix(55388.0) m.fs.pfr.inlet.conc_mol_comp[0, "NaOH"].fix(100.0) m.fs.pfr.inlet.conc_mol_comp[0, "EthylAcetate"].fix(100.0) m.fs.pfr.inlet.conc_mol_comp[0, "SodiumAcetate"].fix(0.0) m.fs.pfr.inlet.conc_mol_comp[0, "Ethanol"].fix(0.0) m.fs.pfr.inlet.temperature.fix(303.15) m.fs.pfr.inlet.pressure.fix(101325.0) m.fs.pfr.control_volume.length.fix(0.5) m.fs.pfr.control_volume.area.fix(0.1) assert degrees_of_freedom(m) == 0 m.fs.pfr.initialize(outlvl=5, optarg={'tol': 1e-6}) assert (pytest.approx(101325.0, abs=1e-2) == m.fs.pfr.outlet.pressure[0].value) assert (pytest.approx(303.15, abs=1e-2) == m.fs.pfr.outlet.temperature[0].value)