Flash Unit and HDA Flowsheet Tutorial Revisions to reflect changes to initialization and scaling updates as well as some overall improvements

.github/workflows/core.yml

@@ -1,5 +1,4 @@
 name: Tests
-
 on:
   push:
     branches:
@@ -88,7 +87,7 @@ jobs:
           - win64
         include:
           - os: linux
-            runner-image: ubuntu-20.04
+            runner-image: ubuntu-latest
           - os: win64
             runner-image: windows-2022
     steps:
@@ -109,10 +108,10 @@ jobs:
         run: |
           pwd
           ls idaes_examples
-          pytest -v idaes_examples --ignore=idaes_examples/notebooks/docs/surrogates/sco2/alamo/ 
+          pytest -v idaes_examples --ignore=idaes_examples/notebooks/docs/surrogates/sco2/alamo/
       - name: Upload pytest-xdist worker logs
         if: success() || failure()
-        uses: actions/upload-artifact@v3
+        uses: actions/upload-artifact@v4
         with:
-          name: pytest_worker_logs
+          name: pytest-xdist-logs-${{ matrix.python-version }}-${{ runner.os }}
           path: "tests_*.log"

idaes_examples/mod/hda/hda_ideal_VLE.py

@@ -23,6 +23,7 @@
     log,
     NonNegativeReals,
     Var,
+    value,
     Set,
     Param,
     sqrt,
@@ -45,6 +46,7 @@
     LiquidPhase,
     VaporPhase,
 )
+from idaes.core.util.model_statistics import degrees_of_freedom
 from idaes.core.util.constants import Constants as const
 from idaes.core.util.initialization import fix_state_vars, solve_indexed_blocks
 from idaes.core.initialization import InitializerBase
@@ -548,10 +550,10 @@ def fix_initialization_states(blk):
         # Fix state variables
         fix_state_vars(blk)
 
-        # Also need to deactivate sum of mole fraction constraint
-        for k in blk.values():
-            if not k.config.defined_state:
-                k.equilibrium_constraint.deactivate()
+        # # Also need to deactivate sum of mole fraction constraint
+        # for k in blk.values():
+        #     if not k.config.defined_state:
+        #         k.equilibrium_constraint.deactivate()
 
 
 @declare_process_block_class("IdealStateBlock", block_class=_IdealStateBlock)
@@ -1011,7 +1013,7 @@ def calculate_dew_point_pressure(self, clear_components=True):
     # as needed
     def _temperature_bubble(self):
         self.temperature_bubble = Param(
-            initialize=33.0, units=pyunits.K, doc="Bubble point temperature"
+            initialize=298.15, units=pyunits.K, doc="Bubble point temperature"
         )
 
     def _temperature_dew(self):
@@ -1039,7 +1041,7 @@ def rule_psat_dew(b, j):
             )
 
             def rule_temp_dew(b):
-                return (
+                ans = (
                     b.pressure
                     * sum(
                         b.mole_frac_comp[i] / b._p_sat_dewT[i]
@@ -1048,6 +1050,7 @@ def rule_temp_dew(b):
                     - 1
                     == 0
                 )
+                return ans
 
             self.eq_temperature_dew = Constraint(rule=rule_temp_dew)
         except AttributeError: