Parameterize 'Peaking Boiler Cost per KW' to address NREL#392

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -1111,6 +1111,18 @@ def __init__(self, model: Model):
             ErrMessage="assume default peaking boiler efficiency (85%)",
             ToolTipText="Peaking boiler efficiency"
         )
+        self._default_peaking_boiler_cost_USD_per_kw = 65
+        self.peaking_boiler_cost_per_kw = self.ParameterDict[self.peaking_boiler_cost_per_kw.Name] = floatParameter(
+            "Peaking Boiler Cost per KW",
+            DefaultValue=self._default_peaking_boiler_cost_USD_per_kw,
+            Min=0,
+            Max=1000,
+            UnitType=Units.ENERGYCOST,
+            PreferredUnits=EnergyCostUnit.DOLLARSPERKW,
+            CurrentUnits=EnergyCostUnit.DOLLARSPERKW,
+            Required=False,
+            ToolTipText="Peaking boiler cost per KW of maximum peaking boiler demand"
+        )
         self.dhpipingcostrate = self.ParameterDict[self.dhpipingcostrate.Name] = floatParameter(
             "District Heating Piping Cost Rate",
             DefaultValue=1200,
@@ -1689,13 +1701,18 @@ def __init__(self, model: Model):
             PreferredUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR,
             CurrentUnits=CurrencyFrequencyUnit.MDOLLARSPERYEAR
         )
+
         # district heating
         self.peakingboilercost = self.OutputParameterDict[self.peakingboilercost.Name] = OutputParameter(
             Name="Peaking boiler cost",
             UnitType=Units.CURRENCY,
             PreferredUnits=CurrencyUnit.MDOLLARS,
-            CurrentUnits=CurrencyUnit.MDOLLARS
+            CurrentUnits=CurrencyUnit.MDOLLARS,
+            ToolTipText=f'Default cost: ${self._default_peaking_boiler_cost_USD_per_kw}/KW '
+                        f'of maximum peaking boiler demand. '
+                        f'Provide {self.peaking_boiler_cost_per_kw.Name} override the default.'
         )
+
         self.dhdistrictcost = self.OutputParameterDict[self.dhdistrictcost.Name] = OutputParameter(
             Name="District Heating System Cost",
             UnitType=Units.CURRENCY,
@@ -2386,8 +2403,13 @@ def Calculate(self, model: Model) -> None:
                 self.Cplant.value = 1.12 * 1.15 * self.ccplantadjfactor.value * 250E-6 * np.max(
                     model.surfaceplant.HeatExtracted.value) * 1000.
 
-                self.peakingboilercost.value = 65 * model.surfaceplant.max_peaking_boiler_demand.value / 1000  # add 65$/KW for peaking boiler
-                self.Cplant.value += self.peakingboilercost.value  # add peaking boiler cost to surface plant cost
+                # add 65$/KW for peaking boiler
+                self.peakingboilercost.value = (self.peaking_boiler_cost_per_kw.quantity()
+                                                .to('USD / kilowatt').magnitude
+                                                * model.surfaceplant.max_peaking_boiler_demand.value / 1000)
+
+                # add peaking boiler cost to surface plant cost
+                self.Cplant.value += self.peakingboilercost.value
 
 
         else:  # all other options have power plant

src/geophires_x_schema_generator/geophires-request.json

@@ -1861,6 +1861,15 @@
       "minimum": 0,
       "maximum": 1
     },
+    "Peaking Boiler Cost per KW": {
+      "description": "Peaking boiler cost per KW of maximum peaking boiler demand",
+      "type": "number",
+      "units": "USD/kW",
+      "category": "Economics",
+      "default": 65,
+      "minimum": 0,
+      "maximum": 1000
+    },
     "District Heating Piping Cost Rate": {
       "description": "District heating piping cost rate ($/m)",
       "type": "number",

tests/geophires_x_tests/test_economics.py

@@ -7,6 +7,7 @@
 # ruff: noqa: I001  # Successful module initialization is dependent on this specific import order.
 from geophires_x.Model import Model
 from geophires_x.Economics import CalculateFinancialPerformance
+from geophires_x_client import GeophiresXResult, GeophiresXClient, GeophiresInputParameters
 from tests.base_test_case import BaseTestCase
 
 
@@ -115,3 +116,24 @@ def _new_model(self) -> Model:
         os.chdir(stash_cwd)
 
         return m
+
+    def test_peaking_boiler_cost(self):
+        def _get_result(peaking_boiler_cost_: int) -> GeophiresXResult:
+            return GeophiresXClient().get_geophires_result(
+                GeophiresInputParameters(
+                    from_file_path=self._get_test_file_path('../examples/example12_DH.txt'),
+                    params={
+                        'Peaking Boiler Cost per KW': peaking_boiler_cost_,
+                    },
+                )
+            )
+
+        def _lcoh_pbc(r: GeophiresXResult) -> tuple[float, float]:
+            return (
+                r.result['SUMMARY OF RESULTS']['Direct-Use heat breakeven price (LCOH)']['value'],
+                r.result['CAPITAL COSTS (M$)']['of which Peaking Boiler Cost']['value'],
+            )
+
+        lcoh, peaking_boiler_cost = _lcoh_pbc(_get_result(0))
+        self.assertLess(lcoh, 13.19)
+        self.assertEqual(0, peaking_boiler_cost)