Updates and Improvements to REFLO System Costing

src/watertap_contrib/reflo/analysis/case_studies/KBHDP/__init__.py

@@ -1,2 +1 @@
 from .components import *
-from .KBHDP_SOA import *

src/watertap_contrib/reflo/analysis/multiperiod/ltmed_vagmd_semibatch/MED_VAGMD_semibatch_class.py

@@ -500,6 +500,9 @@ def add_costing_packages(self):
         """
         self.costing = TreatmentCosting()
         self.costing.base_currency = pyunits.USD_2020
+        # set heat and electricity costs to be non-zero
+        self.costing.heat_cost.set_value(0.01)
+        self.costing.electricity_cost.fix(0.07)
 
         # The costing model is built upon the last time step
         blk = self.mp.get_active_process_blocks()[-1].fs

src/watertap_contrib/reflo/analysis/multiperiod/ltmed_vagmd_semibatch/test/test_MED_VAGMD_semibatch_flowsheet_multiperiod.py

@@ -260,11 +260,11 @@ def test_costing(self, MED_VAGMD_semibatch_frame):
             68657.789, rel=1e-3
         )
         assert cost_performance["Annual heat cost ($)"] == pytest.approx(
-            150540.917, rel=1e-3
+            128236.19, rel=1e-3
         )
         assert cost_performance["Annual electricity cost ($)"] == pytest.approx(
             4064.964, rel=1e-3
         )
         assert cost_performance["Overall LCOW ($/m3)"] == pytest.approx(
-            1.76369, rel=1e-3
+            1.70261, rel=1e-3
         )

src/watertap_contrib/reflo/analysis/multiperiod/vagmd_batch/test/test_VAGMD_batch_flowsheet_multiperiod.py

@@ -28,7 +28,7 @@
 
 from watertap.core.solvers import get_solver
 
-from watertap_contrib.reflo.costing import REFLOCosting
+from watertap_contrib.reflo.costing import TreatmentCosting
 from watertap_contrib.reflo.analysis.multiperiod.vagmd_batch.VAGMD_batch_multiperiod_unit_model import (
     VAGMDbatchSurrogate,
 )
@@ -248,7 +248,10 @@ def test_costing(self, VAGMD_batch_frame_AS7C15L_Closed):
         # The costing model is built upon the last time step
         vagmd = m.fs.VAGMD.mp.get_active_process_blocks()[-1].fs.vagmd
 
-        m.fs.costing = REFLOCosting()
+        m.fs.costing = TreatmentCosting()
+        # set heat and electricity costs to be non-zero
+        m.fs.costing.heat_cost.set_value(0.01)
+        m.fs.costing.electricity_cost.fix(0.07)
         m.fs.costing.base_currency = pyunits.USD_2020
 
         m.fs.VAGMD.add_costing_module(m.fs.costing)
@@ -277,7 +280,7 @@ def test_costing(self, VAGMD_batch_frame_AS7C15L_Closed):
         assert pytest.approx(151892.658, rel=1e-3) == value(
             vagmd.costing.fixed_operating_cost
         )
-        assert pytest.approx(2.777, rel=1e-3) == value(m.fs.costing.LCOW)
+        assert pytest.approx(2.500, rel=1e-3) == value(m.fs.costing.LCOW)
 
 
 class TestVAGMDbatchAS7C15L_HighSalinityClosed:

src/watertap_contrib/reflo/costing/solar/photovoltaic.py

@@ -92,6 +92,7 @@ def build_photovoltaic_cost_param_block(blk):
 )
 def cost_pv(blk):
 
+    blk.costing_package.has_electricity_generation = True
     global_params = blk.costing_package
     pv_params = blk.costing_package.photovoltaic
     pv = blk.unit_model