Syncing GreenSteel and GreenHEART features

examples/reference_plants/01-onshore-steel-mn/onshore-steel-mn.py

@@ -57,4 +57,4 @@
 
     print("LCOE: ", lcoe, "[$/MWh]")
     print("LCOH: ", lcoh, "[$/kg]")
-    print("LCOS: ", lcos, "[$/metric-tonne]")
+    print("LCOS: ", lcos, "[$/metric-tonne]")

examples/reference_plants/02-onshore-ammonia-tx/onshore-ammonia-tx.py

@@ -57,4 +57,4 @@
 
     print("LCOE: ", lcoe, "[$/MWh]")
     print("LCOH: ", lcoh, "[$/kg]")
-    print("LCOA: ", lcoa, "[$/kg]")
+    print("LCOA: ", lcoa, "[$/kg]")

examples/reference_plants/03-offshore-hydrogen-gom/offshore-hydrogen-gom.py

@@ -23,6 +23,7 @@
 
 # ORBIT imports
 from ORBIT.core.library import initialize_library
+
 initialize_library(os.path.join(os.getcwd(), "./input/"))
 
 # run the stuff
@@ -32,7 +33,7 @@
     filename_turbine_config = "./input/turbines/" + turbine_model + ".yaml"
     filename_floris_config = "./input/floris/floris_input_osw_17MW.yaml"
     filename_hopp_config = "./input/plant/hopp_config_gom.yaml"
-    filename_orbit_config= "./input/plant/orbit-config-"+turbine_model+"-gom.yaml"
+    filename_orbit_config = "./input/plant/orbit-config-" + turbine_model + "-gom.yaml"
     filename_greenheart_config = "./input/plant/greenheart_config_offshore_gom.yaml"
 
     config = GreenHeartSimulationConfig(
@@ -61,4 +62,4 @@
     lcoh = prob.get_val("lcoh", units="USD/kg")
 
     print("LCOE: ", lcoe, "[$/MWh]")
-    print("LCOH: ", lcoh, "[$/kg]")
+    print("LCOH: ", lcoh, "[$/kg]")

examples/reference_plants/04-offshore-hydrogen-ny/offshore-hydrogen-ny.py

@@ -23,6 +23,7 @@
 
 # ORBIT imports
 from ORBIT.core.library import initialize_library
+
 initialize_library(os.path.join(os.getcwd(), "./input/"))
 
 # run the stuff
@@ -32,7 +33,7 @@
     filename_turbine_config = "./input/turbines/" + turbine_model + ".yaml"
     filename_floris_config = "./input/floris/floris_input_osw_15MW.yaml"
     filename_hopp_config = "./input/plant/hopp_config_ny.yaml"
-    filename_orbit_config= "./input/plant/orbit-config-"+turbine_model+"-ny.yaml"
+    filename_orbit_config = "./input/plant/orbit-config-" + turbine_model + "-ny.yaml"
     filename_greenheart_config = "./input/plant/greenheart_config_offshore_ny.yaml"
 
     config = GreenHeartSimulationConfig(
@@ -61,4 +62,4 @@
     lcoh = prob.get_val("lcoh", units="USD/kg")
 
     print("LCOE: ", lcoe, "[$/MWh]")
-    print("LCOH: ", lcoh, "[$/kg]")
+    print("LCOH: ", lcoh, "[$/kg]")

examples/reference_plants/05-offshore-hydrogen-ca/offshore-hydrogen-ca.py

@@ -23,6 +23,7 @@
 
 # ORBIT imports
 from ORBIT.core.library import initialize_library
+
 initialize_library(os.path.join(os.getcwd(), "./input/"))
 
 # run the stuff
@@ -32,7 +33,7 @@
     filename_turbine_config = "./input/turbines/" + turbine_model + ".yaml"
     filename_floris_config = "./input/floris/floris_input_osw_15MW.yaml"
     filename_hopp_config = "./input/plant/hopp_config_ca.yaml"
-    filename_orbit_config= "./input/plant/orbit-config-"+turbine_model+"-ca.yaml"
+    filename_orbit_config = "./input/plant/orbit-config-" + turbine_model + "-ca.yaml"
     filename_greenheart_config = "./input/plant/greenheart_config_offshore_ca.yaml"
 
     config = GreenHeartSimulationConfig(
@@ -61,4 +62,4 @@
     lcoh = prob.get_val("lcoh", units="USD/kg")
 
     print("LCOE: ", lcoe, "[$/MWh]")
-    print("LCOH: ", lcoh, "[$/kg]")
+    print("LCOH: ", lcoh, "[$/kg]")

greenheart/simulation/greenheart_simulation.py

@@ -34,6 +34,7 @@
 import greenheart.tools.eco.utilities as he_util
 import greenheart.tools.eco.hydrogen_mgmt as he_h2
 
+
 @define
 class GreenHeartSimulationConfig:
     """
@@ -157,10 +158,10 @@ def __attrs_post_init__(self):
             self.orbit_config["plant"]["num_turbines"] = int(
                 self.wind_rating * 1e-3 / self.turbine_config["turbine_rating"]
             )
-            self.hopp_config["technologies"]["wind"]["num_turbines"] = (
-                self.orbit_config["plant"]["num_turbines"]
-            )
-        
+            self.hopp_config["technologies"]["wind"][
+                "num_turbines"
+            ] = self.orbit_config["plant"]["num_turbines"]
+
         if self.grid_connection != None:
             self.greenheart_config["project_parameters"][
                 "grid_connection"
@@ -238,38 +239,74 @@ def setup_greenheart_simulation(config: GreenHeartSimulationConfig):
     # run orbit for wind plant construction and other costs
     ## TODO get correct weather (wind, wave) inputs for ORBIT input (possibly via ERA5)
     if config.design_scenario["wind_location"] == "offshore":
-
-        if config.orbit_config["plant"]["num_turbines"] != config.hopp_config["technologies"]["wind"]["num_turbines"]:
+
+        if (
+            config.orbit_config["plant"]["num_turbines"]
+            != config.hopp_config["technologies"]["wind"]["num_turbines"]
+        ):
             config.orbit_config["plant"].update(
-                {"num_turbines": config.hopp_config["technologies"]["wind"]["num_turbines"]}
+                {
+                    "num_turbines": config.hopp_config["technologies"]["wind"][
+                        "num_turbines"
+                    ]
+                }
             )
-            warnings.warn(f"'num_turbines' in the orbit_config was {config.orbit_config['plant']['num_turbines']}, but 'num_turbines' in" 
-                    f"hopp_config was {config.hopp_config['technologies']['wind']['num_turbines']}. The 'num_turbines' value in the orbit_config"
-                    "is being overwritten with the value from the hopp_config", UserWarning)
-
-        if config.orbit_config["site"]["depth"] != config.greenheart_config["site"]["depth"]:
+            warnings.warn(
+                f"'num_turbines' in the orbit_config was {config.orbit_config['plant']['num_turbines']}, but 'num_turbines' in"
+                f"hopp_config was {config.hopp_config['technologies']['wind']['num_turbines']}. The 'num_turbines' value in the orbit_config"
+                "is being overwritten with the value from the hopp_config",
+                UserWarning,
+            )
+
+        if (
+            config.orbit_config["site"]["depth"]
+            != config.greenheart_config["site"]["depth"]
+        ):
             config.orbit_config["site"].update(
                 {"depth": config.greenheart_config["site"]["depth"]}
             )
-            warnings.warn(f"site depth in the orbit_config was {config.orbit_config['site']['depth']}, but site depth in" 
-                    f"greenheart_config was {config.greenheart_config['site']['depth']}. The site depth value in the orbit_config"
-                    "is being overwritten with the value from the greenheart_config", UserWarning)
-
-        if config.orbit_config["plant"]["turbine_spacing"] != config.greenheart_config["site"]["wind_layout"]["turbine_spacing"]:
-            config.orbit_config["plant"].update(
-                {"turbine_spacing": config.greenheart_config["site"]["wind_layout"]["turbine_spacing"]}
+            warnings.warn(
+                f"site depth in the orbit_config was {config.orbit_config['site']['depth']}, but site depth in"
+                f"greenheart_config was {config.greenheart_config['site']['depth']}. The site depth value in the orbit_config"
+                "is being overwritten with the value from the greenheart_config",
+                UserWarning,
             )
             warnings.warn(f"'turbine_spacing' in the orbit_config was {config.orbit_config['plant']['turbine_spacing']}, but 'turbine_spacing' in" 
                     f"greenheart_config was {config.greenheart_config['site']['wind_layout']['turbine_spacing']}. The 'turbine_spacing' value in the orbit_config"
                     "is being overwritten with the value from the greenheart_config", UserWarning)
 
         if config.orbit_config["plant"]["row_spacing"] != config.greenheart_config["site"]["wind_layout"]["row_spacing"]:
             config.orbit_config["plant"].update(
-                {"row_spacing": config.greenheart_config["site"]["wind_layout"]["row_spacing"]}
+                {
+                    "turbine_spacing": config.greenheart_config["site"]["wind_layout"][
+                        "turbine_spacing"
+                    ]
+                }
+            )
+            warnings.warn(
+                f"site depth in the orbit_config was {config.orbit_config['plant']['turbine_spacing']}, but 'turbine_spacing' in"
+                f"greenheart_config was {config.greenheart_config['site']['wind_layout']['turbine_spacing']}. The 'turbine_spacing' value in the orbit_config"
+                "is being overwritten with the value from the greenheart_config",
+                UserWarning,
+            )
+
+        if (
+            config.orbit_config["plant"]["row_spacing"]
+            != config.greenheart_config["site"]["wind_layout"]["row_spacing"]
+        ):
+            config.orbit_config["plant"].update(
+                {
+                    "row_spacing": config.greenheart_config["site"]["wind_layout"][
+                        "row_spacing"
+                    ]
+                }
+            )
+            warnings.warn(
+                f"site depth in the orbit_config was {config.orbit_config['plant']['row_spacing']}, but 'row_spacing' in"
+                f"greenheart_config was {config.greenheart_config['site']['wind_layout']['row_spacing']}. The 'row_spacing' value in the orbit_config"
+                "is being overwritten with the value from the greenheart_config",
+                UserWarning,
             )
-            warnings.warn(f"site depth in the orbit_config was {config.orbit_config['plant']['row_spacing']}, but 'row_spacing' in" 
-                    f"greenheart_config was {config.greenheart_config['site']['wind_layout']['row_spacing']}. The 'row_spacing' value in the orbit_config"
-                    "is being overwritten with the value from the greenheart_config", UserWarning)
 
         wind_config = he_fin.WindCostConfig(
             design_scenario=config.design_scenario,
@@ -368,7 +405,7 @@ def run_simulation(config: GreenHeartSimulationConfig):
             hopp_config=config.hopp_config,
             greenheart_config=config.greenheart_config,
             turbine_config=config.turbine_config,
-            hopp_interface=hopp_results["hopp_interface"]
+            hopp_interface=hopp_results["hopp_interface"],
         )
 
         wind_cost_results = he_fin.run_wind_cost_model(
@@ -457,13 +494,14 @@ def energy_internals(
         )
 
         # compressor #TODO size correctly
-        h2_transport_compressor, h2_transport_compressor_results = (
-            he_h2.run_h2_transport_compressor(
-                config.greenheart_config,
-                electrolyzer_physics_results,
-                design_scenario,
-                verbose=verbose,
-            )
+        (
+            h2_transport_compressor,
+            h2_transport_compressor_results,
+        ) = he_h2.run_h2_transport_compressor(
+            config.greenheart_config,
+            electrolyzer_physics_results,
+            design_scenario,
+            verbose=verbose,
         )
 
         # transport pipeline
@@ -814,6 +852,9 @@ def simple_solver(initial_guess=0.0):
             if config.verbose:
                 print("Running ammonia\n")
 
+            if "hydrogen_cost" not in config.greenheart_config['ammonia']['costs']['feedstocks']:
+                config.greenheart_config['ammonia']['costs']['feedstocks']['hydrogen_cost'] = lcoh
+
             # use the hydrogen amount from the electrolyzer physics model if it is not already in the config
             if (
                 "hydrogen_amount_kgpy"
@@ -916,7 +957,13 @@ def simple_solver(initial_guess=0.0):
             platform_results = platform_results
         )
 
-def run_sweeps(simulate=False, verbose=True, show_plots=True, use_profast=True, output_dir="output/"):
+def run_sweeps(
+    simulate=False,
+    verbose=True,
+    show_plots=True,
+    use_profast=True,
+    output_dir="output/",
+):
 
     if simulate:
         verbose = False
@@ -942,7 +989,8 @@ def run_sweeps(simulate=False, verbose=True, show_plots=True, use_profast=True,
                     )
                     print(lcoh_array)
                 np.savetxt(
-                    output_dir + "data/lcoh_vs_rating_%s_storage_%sMWwindplant.txt"
+                    output_dir
+                    + "data/lcoh_vs_rating_%s_storage_%sMWwindplant.txt"
                     % (storage_type, wind_rating),
                     np.c_[ratings, lcoh_array],
                 )
@@ -1033,7 +1081,11 @@ def run_sweeps(simulate=False, verbose=True, show_plots=True, use_profast=True,
 
 
 def run_policy_options_storage_types(
-    verbose=True, show_plots=False, save_plots=False, use_profast=True, output_dir="output/"
+    verbose=True,
+    show_plots=False,
+    save_plots=False,
+    use_profast=True,
+    output_dir="output/",
 ):
 
     storage_types = ["pressure_vessel", "pipe", "salt_cavern", "none"]
@@ -1066,7 +1118,11 @@ def run_policy_options_storage_types(
 
 
 def run_policy_storage_design_options(
-    verbose=False, show_plots=False, save_plots=False, use_profast=True, output_dir="output/"
+    verbose=False,
+    show_plots=False,
+    save_plots=False,
+    use_profast=True,
+    output_dir="output/",
 ):
 
     design_scenarios = [1, 2, 3, 4, 5, 6, 7]
@@ -1156,7 +1212,11 @@ def run_policy_storage_design_options(
 
 
 def run_design_options(
-    verbose=False, show_plots=False, save_plots=False, incentive_option=1, output_dir="output/"
+    verbose=False,
+    show_plots=False,
+    save_plots=False,
+    incentive_option=1,
+    output_dir="output/",
 ):
 
     design_options = range(1, 8)  # 8
@@ -1232,7 +1292,7 @@ def run_storage_options(output_dir="output/"):
             storage_type=storage_type,
             output_level=4,
             grid_connection=False,
-            output_dir=output_dir
+            output_dir=output_dir,
         )
         lcoe_list.append(lcoe)
         lcoh_list.append(lcoh)
@@ -1248,7 +1308,7 @@ def run_storage_options(output_dir="output/"):
             storage_type=storage_type,
             output_level=4,
             grid_connection=True,
-            output_dir=output_dir
+            output_dir=output_dir,
         )
         lcoh_with_grid_list.append(lcoh_with_grid)
         lcoh_grid_only_list.append(lcoh_grid_only)

greenheart/tools/eco/electrolysis.py

@@ -3,6 +3,7 @@
 import matplotlib.pyplot as plt
 from matplotlib import ticker
 import os
+import warnings
 from greenheart.tools.eco.utilities import ceildiv
 
 # import hopp.tools.hopp_tools as hopp_tools
@@ -87,6 +88,9 @@ def run_electrolyzer_physics(
         "turndown_ratio": greenheart_config["electrolyzer"]["turndown_ratio"],
     }
 
+    if "time_between_replacement" in greenheart_config['electrolyzer']:
+        warnings.warn("`time_between_replacement` as an input is depricated. It is now calculated internally and is output in electrolyzer_physics_results['H2_Results']['Time Until Replacement [hrs]'].")
+
     H2_Results, h2_ts, h2_tot, power_to_electrolyzer_kw = run_h2_PEM(
         electrical_generation_timeseries=energy_to_electrolyzer_kw,
         electrolyzer_size=electrolyzer_size_mw,
@@ -392,9 +396,10 @@ def run_electrolyzer_cost(
 
             pem_offshore = PEMCostsSingliticoModel(elec_location=offshore)
 
-            electrolyzer_capital_cost_musd, electrolyzer_om_cost_musd = (
-                pem_offshore.run(P_elec, RC_elec)
-            )
+            (
+                electrolyzer_capital_cost_musd,
+                electrolyzer_om_cost_musd,
+            ) = pem_offshore.run(P_elec, RC_elec)
 
             electrolyzer_total_capital_cost = (
                 electrolyzer_capital_cost_musd * 1e6