Updating NGCC Flowsheet

idaes_examples/mod/power_gen/gas_turbine.py

@@ -234,7 +234,7 @@ def _add_models(self):
             support_isentropic_performance_curves=True,
         )
 
-    def _add_performance_curves_gts1(self, flow_scale=0.896):
+    def _add_performance_curves_gts1(self, flow_scale=0.896*pyo.units.s/pyo.units.m**3):
         """Add isentropic head and efficiency curves for gas turbine stage 1"""
 
         @self.gts1.performance_curve.Constraint(
@@ -262,9 +262,9 @@ def head_isen_eqn(b, t):
             )
             return b.head_isentropic[t] == -(
                 -2085.1 * f**3 + 38433 * f**2 - 150764 * f + 422313
-            )
+            )*pyo.units.m**2/pyo.units.s**2
 
-    def _add_performance_curves_gts2(self, flow_scale=0.896):
+    def _add_performance_curves_gts2(self, flow_scale=0.896*pyo.units.s/pyo.units.m**3):
         """Add isentropic head and efficiency curves for gas turbine stage 2"""
 
         @self.gts2.performance_curve.Constraint(
@@ -290,9 +290,9 @@ def head_isen_eqn(b, t):
             )
             return b.head_isentropic[t] == -(
                 -1676.3 * f**3 + 34916 * f**2 - 173801 * f + 456957
-            )
+            )*pyo.units.m**2/pyo.units.s**2
 
-    def _add_performance_curves_gts3(self, flow_scale=0.896):
+    def _add_performance_curves_gts3(self, flow_scale=0.896*pyo.units.s/pyo.units.m**3):
         """Add isentropic head and efficiency curves for gas turbine stage 3"""
 
         @self.gts3.performance_curve.Constraint(
@@ -320,7 +320,7 @@ def head_isen_eqn(b, t):
             )
             return b.head_isentropic[t] == -(
                 -1373.6 * f**3 + 31759 * f**2 - 188528 * f + 500520
-            )
+            )*pyo.units.m**2/pyo.units.s**2
 
     def _add_constraints(self):
         """Add additional flowsheet constraints and expressions"""
@@ -759,10 +759,10 @@ def _set_scaling(self):
         for i in ["air05", "air07", "air09"]:
             iscale.set_scaling_factor(self.splt1.split_fraction[0.0, i], 100)
 
-        iscale.set_scaling_factor(self.valve01.control_volume.work, 1e-8)
-        iscale.set_scaling_factor(self.valve02.control_volume.work, 1e-8)
-        iscale.set_scaling_factor(self.valve03.control_volume.work, 1e-8)
-        iscale.set_scaling_factor(self.vsv.control_volume.work, 1e-8)
+        iscale.set_scaling_factor(self.valve01.control_volume.work, 1)
+        iscale.set_scaling_factor(self.valve02.control_volume.work, 1)
+        iscale.set_scaling_factor(self.valve03.control_volume.work, 1)
+        iscale.set_scaling_factor(self.vsv.control_volume.work, 1)
         iscale.set_scaling_factor(self.cmp1.control_volume.work, 1e-8)
         iscale.set_scaling_factor(self.gts1.control_volume.work, 1e-8)
         iscale.set_scaling_factor(self.gts2.control_volume.work, 1e-8)
@@ -852,8 +852,10 @@ def initialize(
         outlvl=idaeslog.NOTSET,
         solver=None,
         optarg=None,
-        load_from="gas_turbine_init.json.gz",
-        save_to="gas_turbine_init.json.gz",
+        # load_from="gas_turbine_init.json.gz",
+        # save_to="gas_turbine_init.json.gz",
+        load_from=None,
+        save_to=None,
     ):
         """Initialize the gas turbine flowsheet
 
@@ -889,6 +891,7 @@ def initialize(
         self.feed_fuel1.initialize(outlvl=outlvl, solver=solver, optarg=optarg)
         # compressor
         propagate_state(self.air01)
+        self.vsv.valve_opening.fix(1-1e-6)
         self.vsv.initialize(outlvl=outlvl, solver=solver, optarg=optarg)
         propagate_state(self.air02)
         self.cmp1.initialize(outlvl=outlvl, solver=solver, optarg=optarg)

idaes_examples/mod/power_gen/hrsg.py

@@ -93,16 +93,19 @@ def _add_unit_models(self):
         prop_water = self.prop_water
         prop_gas = self.prop_gas
 
+
         ######### LP Section ###########
         self.econ_lp = BoilerHeatExchanger(
             doc="LP Economizer",
             delta_temperature_callback=delta_temp_cb,
             cold_side_name="tube",
             hot_side_name="shell",
             tube={"property_package": prop_water,
-                  "has_pressure_change": False,},
+                  "has_pressure_change": False,
+                  },
             shell={"property_package": prop_gas,
-                   "has_pressure_change": False,},
+                   "has_pressure_change": False,
+                   },
             has_holdup=False,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
             tube_arrangement=TubeArrangement.inLine,
@@ -115,10 +118,12 @@ def _add_unit_models(self):
             property_package=prop_water,
             momentum_mixing_type=MomentumMixingType.none,
             inlet_list=["econ_lp", "Preheater"],
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.drum_lp = HelmPhaseSeparator(
             doc="Phase separator for LP evaporator (partial evaporator)",
             property_package=prop_water,
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.evap_lp = HeatExchanger(
             doc="LP evaporator heat exchanger section",
@@ -135,6 +140,7 @@ def _add_unit_models(self):
             property_package=prop_water,
             momentum_mixing_type=MomentumMixingType.none,
             inlet_list=["main", "soec_makeup"],
+            property_package_args={'has_phase_equilibrium':False}            
         )
         self.split_fg_lp = Splitter(
             doc="LP superheater flue bypass gas splitter",
@@ -166,14 +172,19 @@ def _add_unit_models(self):
             doc="LP liquid split to IP and HP pumps",
             property_package=prop_water,
             outlet_list=["toIP", "toHP"],
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.pump_ip = WaterPump(
             doc="Intermediate pressure pump",
             property_package=prop_water,
+            has_phase_equilibrium=False,
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.pump_hp = WaterPump(
             doc="High pressure pump",
             property_package=prop_water,
+            has_phase_equilibrium=False,
+            property_package_args={'has_phase_equilibrium':False}
         )
         ######### IP Section ###########
         self.econ_ip1 = BoilerHeatExchanger(
@@ -182,9 +193,9 @@ def _add_unit_models(self):
             hot_side_name="shell",
             cold_side_name="tube",
             tube={"property_package": prop_water,
-                  "has_pressure_change": True,},
+                  "has_pressure_change": True},
             shell={"property_package": prop_gas,
-                   "has_pressure_change": True,},
+                   "has_pressure_change": True},
             has_holdup=False,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
             tube_arrangement=TubeArrangement.inLine,
@@ -195,21 +206,23 @@ def _add_unit_models(self):
             doc="IP economizer hot water split for natural gas preheater",
             property_package=prop_water,
             outlet_list=["toIP_ECON2", "toNGPH"],
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.econ_ip2 = BoilerHeatExchanger(
             doc="IP ecomonmizer part 2",
             delta_temperature_callback=delta_temp_cb,
             hot_side_name="shell",
             cold_side_name="tube",
             tube={"property_package": prop_water,
-                  "has_pressure_change": True,},
+                  "has_pressure_change": True},
             shell={"property_package": prop_gas,
-                   "has_pressure_change": True,},
+                   "has_pressure_change": True},
             has_holdup=False,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
             tube_arrangement=TubeArrangement.inLine,
             cold_side_water_phase="Liq",
             has_radiation=False,
+            # has_phase_equilibrium=False,
         )
         self.evap_ip = HeatExchanger(
             doc="IP evaporator (total evaporator)",
@@ -219,16 +232,17 @@ def _add_unit_models(self):
             tube={"property_package": prop_water},
             delta_temperature_callback=delta_temperature_lmtd_callback,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
+
         )
         self.sh_ip1 = BoilerHeatExchanger(
             doc="IP superheater 1",
             delta_temperature_callback=delta_temp_cb,
             hot_side_name="shell",
             cold_side_name="tube",
             tube={"property_package": prop_water,
-                  "has_pressure_change": True,},
+                  "has_pressure_change": True},
             shell={"property_package": prop_gas,
-                   "has_pressure_change": True,},
+                   "has_pressure_change": True},
             has_holdup=False,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
             tube_arrangement=TubeArrangement.inLine,
@@ -241,21 +255,24 @@ def _add_unit_models(self):
             property_package=prop_water,
             momentum_mixing_type=MomentumMixingType.none,
             inlet_list=["sh_ip1", "Cold_reheat"],
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.splitter_ip2 = HelmSplitter(
             doc="IP Splitter 2, for ejector, reclaimer and dryer",
             property_package=prop_water,
             outlet_list=["Cold_reheat", "toEjector", "toReclaimer", "toDryer"],
+            # has_phase_equilibrium=False,
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.sh_ip2 = BoilerHeatExchanger(
             doc="IP superheater 2",
             delta_temperature_callback=delta_temp_cb,
             hot_side_name="shell",
             cold_side_name="tube",
             tube={"property_package": prop_water,
-                  "has_pressure_change": True,},
+                  "has_pressure_change": True},
             shell={"property_package": prop_gas,
-                   "has_pressure_change": True,},
+                   "has_pressure_change": True},
             has_holdup=False,
             flow_pattern=HeatExchangerFlowPattern.countercurrent,
             tube_arrangement=TubeArrangement.inLine,
@@ -355,6 +372,7 @@ def _add_unit_models(self):
         self.evap_hp_valve = HelmValve(
             doc="HP evaporator valve",
             property_package=prop_water,
+            property_package_args={'has_phase_equilibrium':False}
         )
         self.evap_hp_valve.pressure_flow_equation.deactivate()
         self.evap_hp = HeatExchanger(
@@ -436,7 +454,7 @@ def _add_flowsheet_constraints(self):
         def ip_sat_vap_eqn(b, t):
             return (
                 b.tube.properties_out[t].enth_mol / 1e4
-                == (b.tube.properties_out[t].enth_mol_sat_phase["Vap"] + 30) / 1e4
+                == (b.tube.properties_out[t].enth_mol_sat_phase["Vap"] + 30*pyo.units.J/pyo.units.mol) / 1e4
             )
 
         @self.evap_hp.Constraint(
@@ -445,7 +463,7 @@ def ip_sat_vap_eqn(b, t):
         def hp_sat_vap_eqn(b, t):
             return (
                 b.tube.properties_out[t].enth_mol / 1e4
-                == (b.tube.properties_out[t].enth_mol_sat_phase["Vap"] + 30) / 1e4
+                == (b.tube.properties_out[t].enth_mol_sat_phase["Vap"] + 30*pyo.units.J/pyo.units.mol) / 1e4
             )
 
         @self.mixer1.Constraint(self.config.time, doc="Mixed state pressure eqn.")
@@ -1020,6 +1038,8 @@ def _set_scaling_factors(self):
             if hasattr(unit, "deltaP_tube_uturn"):
                 iscale.set_scaling_factor(unit.deltaP_tube_uturn, 1e-2)
 
+        # iscale.set_scaling_factor(self.sh_ip2.N_Pr_tube_eqn, 1e9)
+        # iscale.set_scaling_factor(self.sh_hp4.N_Pr_tube_eqn, 1e9)
         iscale.set_scaling_factor(self.evap_lp.shell.heat, 1e-8)
         iscale.set_scaling_factor(self.evap_lp.tube.heat, 1e-8)
         iscale.set_scaling_factor(self.evap_lp.tube.heat, 1e-8)
@@ -1058,8 +1078,10 @@ def initialize(
         outlvl=idaeslog.NOTSET,
         solver=None,
         optarg=None,
-        load_from="hrsg_init.json.gz",
-        save_to="hrsg_init.json.gz",
+        # load_from="hrsg_init.json.gz",
+        # save_to="hrsg_init.json.gz",
+        load_from=None,
+        save_to=None,
     ):
         """Initialize the HRSG flowsheet
 

idaes_examples/mod/power_gen/ngcc.py

@@ -360,7 +360,7 @@ def reboiler_duty_eqn(b, t):
 
         @self.Constraint(self.config.time)
         def net_power_constraint(b, t):
-            return b.net_power_mw[t] / 100.0 == -b.net_power[t] / 1e6 / 100.0
+            return b.net_power_mw[t] / 100.0 == pyo.units.convert(-b.net_power[t] / 100.0,to_units=pyo.units.MW)
 
         @self.Constraint(self.config.time)
         def lp_steam_temperature_eqn(b, t):
@@ -418,8 +418,10 @@ def initialize(
         outlvl=idaeslog.NOTSET,
         solver=None,
         optarg=None,
-        load_from="ngcc_init.json.gz",
-        save_to="ngcc_init.json.gz",
+        # load_from="ngcc_init.json.gz",
+        # save_to="ngcc_init.json.gz",
+        load_from=None,
+        save_to=None,
     ):
 
         init_log = idaeslog.getInitLogger(self.name, outlvl, tag="flowsheet")
@@ -438,23 +440,24 @@ def initialize(
             self.cap_additional_reboiler_duty.fix()
             self.fuel_lhv.fix()
             self.fuel_hhv.fix()
-
             self.gt.initialize(
-                load_from="gas_turbine_init.json.gz",
-                save_to="gas_turbine_init.json.gz",
+                # load_from="gas_turbine_init.json.gz",
+                # save_to="gas_turbine_init.json.gz",
+                load_from=None,
+                save_to=None
             )
             propagate_state(self.g08a)
             self.fg_translate.initialize()
             propagate_state(self.g08b, overwrite_fixed=True)
             self.hrsg.initialize(
-                load_from="hrsg_init.json.gz",
-                save_to="hrsg_init.json.gz",
+                load_from=None,
+                save_to=None
             )
             self.hrsg.sh_hp4.shell_inlet.unfix()
             propagate_state(self.t05a, overwrite_fixed=True)
             self.st.initialize(
-                load_from="steam_turbine_init.json.gz",
-                save_to="steam_turbine_init.json.gz",
+                load_from=None,
+                save_to=None
             )
 
             init_log.info(f"Open tears")