Finance models to account for lifetime performance of electrolyzer (follow-on to 463)

CHANGELOG.md

@@ -29,6 +29,7 @@
   - Remove unused dependencies.
   - Fixes typos for skipped folders.
   - Fixes missing dependencies for `gis` modifier used in new iron mapping tests.
+- Added `PerformanceModelBaseClass` and standardized outputs of converter performance models
 
 ## 0.5.1 [December 18, 2025]
 

docs/user_guide/how_to_set_up_an_analysis.md

@@ -142,9 +142,9 @@ h2i_model.run()
 
 # h2i_model.post_process()
 
-# Print the total hydrogen produced by the electrolyzer in kg/year
-total_hydrogen = h2i_model.model.get_val("electrolyzer.total_hydrogen_produced", units="kg/year")[0]
-print(f"Total hydrogen produced by the electrolyzer: {total_hydrogen:.2f} kg/year")
+# Print the average annual hydrogen produced by the electrolyzer in kg/year
+annual_hydrogen = h2i_model.model.get_val("electrolyzer.annual_hydrogen_produced", units="kg/year").mean()
+print(f"Total hydrogen produced by the electrolyzer: {annual_hydrogen:.2f} kg/year")
 ```
 
 This will run the analysis defined in the config files and generate the output files in the through the `post_process` method.
@@ -169,9 +169,9 @@ h2i_model.run()
 # Post-process the results
 # h2i_model.post_process()
 
-# Print the total hydrogen produced by the electrolyzer in kg/year
-total_hydrogen = h2i_model.model.get_val("electrolyzer.total_hydrogen_produced", units="kg/year")[0]
-print(f"Total hydrogen produced by the electrolyzer: {total_hydrogen:.2f} kg/year")
+# Print the average annual hydrogen produced by the electrolyzer in kg/year
+annual_hydrogen = h2i_model.model.get_val("electrolyzer.annual_hydrogen_produced", units="kg/year").mean()
+print(f"Total hydrogen produced by the electrolyzer: {annual_hydrogen:.2f} kg/year")
 ```
 
 This is especially useful when you want to run an H2I model as a script and modify parameters dynamically without changing the original YAML configuration file.

docs/user_guide/run_size_modes.md

@@ -188,10 +188,10 @@ model.run()
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -225,10 +225,10 @@ model.run()
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -270,10 +270,10 @@ model.run()
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -323,10 +323,10 @@ model.run()
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:

examples/08_wind_electrolyzer/plant_config.yaml

@@ -61,14 +61,17 @@ finance_parameters:
   finance_subgroups:
     electricity_profast:
       commodity: "electricity"
+      # commodity_stream: "wind"
       finance_groups: ["profast_model"]
       technologies: ["wind"]
     electricity_custom:
       commodity: "electricity"
+      # commodity_stream: "wind"
       finance_groups: ["custom_model"]
       technologies: ["wind"]
     hydrogen:
       commodity: "hydrogen"
+      # commodity_stream: "electrolyzer"
       commodity_desc: "produced"
       finance_groups: ["custom_model","profast_model"]
       technologies: ["wind", "electrolyzer"]

examples/08_wind_electrolyzer/user_finance_model/simple_lco.py

@@ -19,6 +19,7 @@ def initialize(self):
         self.options.declare("tech_config", types=dict)
         self.options.declare("commodity_type", types=str)
         self.options.declare("description", types=str, default="")
+        self.options.declare("commodity_stream", types=str, default="")  # unused in this model
 
     def setup(self):
         if self.options["commodity_type"] == "electricity":
@@ -41,10 +42,13 @@ def setup(self):
             else f"{LCO_base_str}_{self.options['description']}"
         )
 
+        plant_life = int(self.options["plant_config"]["plant"]["plant_life"])
+
         # add inputs for commodity production and costs
         self.add_input(
             f"total_{self.options['commodity_type']}_produced",
             val=0.0,
+            shape=plant_life,
             units=commodity_units,
         )
         tech_config = self.tech_config = self.options["tech_config"]

examples/25_sizing_modes/run_size_modes.py

@@ -82,10 +82,10 @@ def compute(self, inputs, outputs, discrete_inputs, discrete_outputs):
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -112,10 +112,10 @@ def compute(self, inputs, outputs, discrete_inputs, discrete_outputs):
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -153,10 +153,10 @@ def compute(self, inputs, outputs, discrete_inputs, discrete_outputs):
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:
@@ -198,10 +198,10 @@ def compute(self, inputs, outputs, discrete_inputs, discrete_outputs):
 for value in [
     "electrolyzer.electricity_in",
     "electrolyzer.electrolyzer_size_mw",
-    "electrolyzer.hydrogen_capacity_factor",
+    "electrolyzer.capacity_factor",
     "ammonia.hydrogen_in",
     "ammonia.max_hydrogen_capacity",
-    "ammonia.ammonia_capacity_factor",
+    "ammonia.capacity_factor",
     "finance_subgroup_h2.LCOH",
     "finance_subgroup_nh3.LCOA",
 ]:

examples/28_iron_map/driver_config.yaml

@@ -10,7 +10,7 @@ general:
 driver:
   design_of_experiments:
     flag: true
-    debug_print: true
+    debug_print: false
     generator: "csvgen"
     filename: "ned_reduced_sitelist.csv"
     run_parallel: False

examples/28_iron_map/plant_config.yaml

@@ -1,21 +1,10 @@
 name: "plant_config"
 description: "This plant is located in MN, USA..."
 
-site:
-  latitude: 41.717
-  longitude: -88.398
-
-  # array of polygons defining boundaries with x/y coords
-  boundaries: [
-    {
-      x: [0.0, 1000.0, 1000.0, 0.0],
-      y: [0.0, 0.0, 100.0, 1000.0],
-    },
-    {
-      x: [2000.0, 2500.0, 2000.0],
-      y: [2000.0, 2000.0, 2500.0],
-    }
-  ]
+sites:
+  site:
+    latitude: 41.717
+    longitude: -88.398
 
 # array of arrays containing left-to-right technology
 # interconnections; can support bidirectional connections

examples/test/test_all_examples.py

@@ -72,11 +72,11 @@ def test_steel_example(subtests):
             pytest.approx(
                 model.prob.get_val("finance_subgroup_hydrogen.LCOH_delivered")[0], rel=1e-3
             )
-            == 8.270362492342693
+            == 8.235313509720276
         )
 
     with subtests.test("Check LCOS"):
-        assert pytest.approx(model.prob.get_val("steel.LCOS")[0], rel=1e-3) == 1266.6193378846617
+        assert pytest.approx(model.prob.get_val("steel.LCOS")[0], rel=1e-3) == 1264.2821232584045
 
     with subtests.test("Check total adjusted CapEx"):
         assert (
@@ -175,15 +175,15 @@ def test_simple_ammonia_example(subtests):
     with subtests.test("Check LCOH"):
         assert (
             pytest.approx(model.prob.get_val("finance_subgroup_hydrogen.LCOH")[0], rel=1e-3)
-            == 4.025446
+            == 4.01554337
         )
 
     with subtests.test("Check price of hydrogen"):
         assert (
             pytest.approx(
                 model.prob.get_val("finance_subgroup_hydrogen.price_hydrogen")[0], rel=1e-3
             )
-            == 4.025446
+            == 4.01554337
         )
 
     # Currently underestimated compared to the Reference Design Doc
@@ -279,7 +279,7 @@ def test_ammonia_synloop_example(subtests):
     with subtests.test("Check LCOH"):
         assert (
             pytest.approx(model.prob.get_val("finance_subgroup_h2.LCOH")[0], rel=1e-6)
-            == 4.025385101169759
+            == 4.015543377027795
         )
 
     with subtests.test("Check LCOA"):
@@ -320,7 +320,7 @@ def test_co2h_methanol_example(subtests):
 
     # Check levelized cost of methanol (LCOM)
     with subtests.test("Check CO2 Hydrogenation LCOM"):
-        assert pytest.approx(model.prob.get_val("methanol.LCOM")[0], rel=1e-6) == 1.7555607442
+        assert pytest.approx(model.prob.get_val("methanol.LCOM")[0], rel=1e-6) == 1.751617212520347
 
 
 @pytest.mark.skipif(importlib.util.find_spec("mcm") is None, reason="mcm is not installed")
@@ -356,7 +356,9 @@ def test_wind_h2_opt_example(subtests):
 
     model_init.post_process()
 
-    annual_h20 = model_init.prob.get_val("electrolyzer.total_hydrogen_produced", units="kg/year")[0]
+    annual_h20 = model_init.prob.get_val("electrolyzer.annual_hydrogen_produced", units="kg/year")[
+        0
+    ]
 
     # Create a H2Integrate model
     model = H2IntegrateModel(Path.cwd() / "wind_plant_electrolyzer.yaml")
@@ -422,7 +424,7 @@ def test_wind_h2_opt_example(subtests):
     with subtests.test("Check minimum total hydrogen produced"):
         assert (
             pytest.approx(
-                model.prob.get_val("electrolyzer.total_hydrogen_produced", units="kg/year")[0],
+                model.prob.get_val("electrolyzer.annual_hydrogen_produced", units="kg/year")[0],
                 abs=15000,
             )
             == 29028700
@@ -504,7 +506,7 @@ def test_splitter_wind_doc_h2_example(subtests):
     with subtests.test("Check LCOH"):
         assert (
             pytest.approx(model.prob.get_val("finance_subgroup_hydrogen.LCOH")[0], rel=1e-3)
-            == 9.82319908
+            == 9.80590834551064
         )
 
     with subtests.test("Check LCOC"):
@@ -587,7 +589,7 @@ def test_hydrogen_dispatch_example(subtests):
                 model.prob.get_val("finance_subgroup_all_hydrogen.LCOH", units="USD/kg")[0],
                 rel=1e-5,
             )
-            == 5.380013537850591
+            == 5.6545296034469
         )
 
     with subtests.test("Check dispatched h2 LCOH"):
@@ -851,7 +853,7 @@ def test_wind_solar_electrolyzer_example(subtests):
                 model.prob.get_val("finance_subgroup_hydrogen.LCOH", units="USD/kg")[0],
                 rel=1e-5,
             )
-            == 5.3277923
+            == 5.3063358
         )
 
     wind_generation = model.prob.get_val("wind.electricity_out", units="kW")
@@ -886,9 +888,9 @@ def test_electrolyzer_om_example(subtests):
     with subtests.test("Check LCOE"):
         assert pytest.approx(lcoe, rel=1e-4) == 39.98869
     with subtests.test("Check LCOH with lcoh_financials"):
-        assert pytest.approx(lcoh_with_lcoh_finance, rel=1e-4) == 13.0858012
+        assert pytest.approx(lcoh_with_lcoh_finance, rel=1e-4) == 16.9204156301
     with subtests.test("Check LCOH with lcoe_financials"):
-        assert pytest.approx(lcoh_with_lcoe_finance, rel=1e-4) == 7.9935907
+        assert pytest.approx(lcoh_with_lcoe_finance, rel=1e-4) == 10.3360027653
 
 
 def test_wombat_electrolyzer_example(subtests):
@@ -917,7 +919,7 @@ def test_wombat_electrolyzer_example(subtests):
     with subtests.test("Check LCOH from custom  model"):
         assert pytest.approx(lcoh_with_custom_model, rel=1e-5) == 4.19232346
     with subtests.test("Check LCOH from ProFAST model"):
-        assert pytest.approx(lcoh_with_profast_model, rel=1e-5) == 5.32632237
+        assert pytest.approx(lcoh_with_profast_model, rel=1e-5) == 5.308630730506
     with subtests.test("Check LCOE from custom model"):
         assert pytest.approx(lcoe_with_custom_model, rel=1e-5) == 51.17615298
     with subtests.test("Check LCOE from ProFAST model"):
@@ -1243,7 +1245,7 @@ def test_csvgen_design_of_experiments(subtests):
             min_lcoh_case_num = i
 
     with subtests.test("Min LCOH value"):
-        assert pytest.approx(min_lcoh_val, rel=1e-6) == 4.468258
+        assert pytest.approx(min_lcoh_val, rel=1e-6) == 4.6630144
 
     with subtests.test("Min LCOH case number"):
         assert min_lcoh_case_num == 6
@@ -1310,9 +1312,9 @@ def test_sweeping_solar_sites_doe(subtests):
     for ci, case in enumerate(cases):
         solar_resource_data = case.get_val("site.solar_resource.solar_resource_data")
         lat_lon = f"{case.get_val('site.latitude')[0]} {case.get_val('site.longitude')[0]}"
-        solar_capacity = case.get_design_vars()["solar.system_capacity_DC"]
-        aep = case.get_val("solar.annual_energy", units="MW*h/yr")
-        lcoe = case.get_val("finance_subgroup_electricity.LCOE_optimistic", units="USD/(MW*h)")
+        solar_capacity = case.get_design_vars()["solar.system_capacity_DC"][0]
+        aep = case.get_val("solar.annual_electricity_produced", units="MW*h/yr")[0]
+        lcoe = case.get_val("finance_subgroup_electricity.LCOE_optimistic", units="USD/(MW*h)")[0]
 
         site_res = pd.DataFrame(
             [aep, lcoe, solar_capacity], index=["AEP", "LCOE", "solar_capacity"], columns=[lat_lon]
@@ -1395,7 +1397,8 @@ def test_floris_example(subtests):
     with subtests.test("Distributed wind plant capacity"):
         assert (
             pytest.approx(
-                h2i.prob.get_val("distributed_wind_plant.total_capacity", units="MW"), rel=1e-6
+                h2i.prob.get_val("distributed_wind_plant.rated_electricity_production", units="MW"),
+                rel=1e-6,
             )
             == 66.0
         )
@@ -1405,7 +1408,7 @@ def test_floris_example(subtests):
             pytest.approx(
                 np.sum(
                     h2i.prob.get_val(
-                        "distributed_wind_plant.total_electricity_produced", units="MW*h/yr"
+                        "distributed_wind_plant.total_electricity_produced", units="MW*h"
                     )
                 ),
                 rel=1e-6,
@@ -1433,7 +1436,8 @@ def test_floris_example(subtests):
     with subtests.test("Utility wind plant capacity"):
         assert (
             pytest.approx(
-                h2i.prob.get_val("utility_wind_plant.total_capacity", units="MW"), rel=1e-6
+                h2i.prob.get_val("utility_wind_plant.rated_electricity_production", units="MW"),
+                rel=1e-6,
             )
             == 120.0
         )