Peak load management heuristic control

.pre-commit-config.yaml

@@ -27,7 +27,7 @@ repos:
       - id: yamlfix
         # Exclude YAML files that are used as inputs for testing or examples, or ones for desired schedule in HOPP as they
         # expect misformatted YAML files.
-        exclude: ^(h2integrate/core/test/inputs/.*|examples/11_hybrid_energy_plant/tech_inputs/desired_schedules/.*)$
+        exclude: ^(h2integrate/core/test/inputs/.*|examples/11_hybrid_energy_plant/tech_inputs/desired_schedules/.*|examples/33_hybrid_energy_plant/demand_profiles/.*)$
   - repo: https://github.com/astral-sh/ruff-pre-commit
   # Ruff version.
     rev: v0.8.1

examples/33_peak_load_management/33_peak_load_management.yaml

@@ -0,0 +1,5 @@
+name: H2Integrate_config
+system_summary: Peak load management dispatch
+driver_config: driver_config.yaml
+technology_config: tech_config.yaml
+plant_config: plant_config.yaml

examples/33_peak_load_management/driver_config.yaml

@@ -0,0 +1,5 @@
+name: driver_config
+description: Driver configuration for multivariable streams example
+general:
+  folder_output: outputs
+  create_om_reports: false

examples/33_peak_load_management/plant_config.yaml

@@ -0,0 +1,11 @@
+name: plant_config
+description: Demonstrates multivariable streams with a gas combiner
+plant:
+  plant_life: 30
+  simulation:
+    n_timesteps: 168
+    dt: 3600
+    timezone: -6
+    start_time: 2025/07/01 00:00:00
+technology_interconnections:
+  - [battery, grid_buy, [electricity_out, electricity_in]]

examples/33_peak_load_management/run_peak_load_management.py

@@ -0,0 +1,72 @@
+"""
+Example 33: Peak load management dispatch
+
+This example demonstrates:
+1. Peak load management dispatch open loop control
+2. Battery charging without an input stream, assuming purchase from the grid
+
+"""
+
+import numpy as np
+import matplotlib.pyplot as plt
+
+from h2integrate.core.utilities import build_time_series_from_plant_config
+from h2integrate.core.file_utils import load_yaml
+from h2integrate.core.h2integrate_model import H2IntegrateModel
+
+
+# Create and setup the H2Integrate model
+model = H2IntegrateModel("33_peak_load_management.yaml")
+
+model.setup()
+
+model.run()
+
+plant_config = load_yaml("plant_config.yaml")
+supervisor_demand = np.asarray(
+    load_yaml("demand_profiles/demand_profile_supervisor.yaml"), dtype=float
+)
+secondary_demand = np.asarray(
+    load_yaml("demand_profiles/demand_profile_secondary.yaml"), dtype=float
+)
+
+time_series = build_time_series_from_plant_config(plant_config)
+
+# Example profiles may be shorter than the simulation horizon; plot over shared length.
+n_plot = min(len(time_series), len(supervisor_demand), len(secondary_demand))
+time_plot = time_series[:n_plot]
+
+fig, ax = plt.subplots(4, 1, sharex=True)
+
+ax[0].plot(time_plot, supervisor_demand[:n_plot] * 1e-3, label="Supervisory demand (MW)")
+ax[0].plot(time_plot, secondary_demand[:n_plot] * 1e-3, label="Secondary demand (MW)")
+ax[0].set_ylabel("Power (MW)")
+ax[0].legend(loc="upper right")
+
+ax[1].plot(time_plot, model.prob.get_val("battery.SOC", units="percent"))
+ax[1].set(ylabel="SOC")
+
+ax[2].plot(time_plot, secondary_demand[:n_plot] * 1e-3, label="Original demand (MW)")
+ax[2].plot(
+    time_plot,
+    model.prob.get_val("battery.electricity_out", units="MW"),
+    label="Battery charge/discharge",
+)
+ax[2].set(ylabel="Power (MW)")
+ax[2].legend()
+
+ax[3].plot(time_plot, secondary_demand[:n_plot] * 1e-3, label="Original demand (MW)")
+ax[3].plot(
+    time_plot,
+    model.prob.get_val("battery.unmet_electricity_demand_out", units="MW"),
+    label="New demand profile",
+)
+ax[3].set(ylabel="Power (MW)")
+ax[3].legend()
+
+
+ax[3].tick_params(axis="x", labelrotation=90)
+
+# import pdb; pdb.set_trace()
+plt.tight_layout()
+plt.show()

examples/33_peak_load_management/tech_config.yaml

@@ -0,0 +1,61 @@
+name: technology_config
+description: This plant charges a battery from the grid to reduce peak demand
+technologies:
+  battery:
+    performance_model:
+      model: StoragePerformanceModel
+    cost_model:
+      model: ATBBatteryCostModel
+    control_strategy:
+      model: PeakLoadManagementOpenLoopStorageController
+    model_inputs:
+      shared_parameters:
+        commodity: electricity
+        commodity_rate_units: kW
+        max_charge_rate: 2500.0  # kW/time step, 1, 2.5, or 5 MW
+        max_capacity: 10000.0  # kWh, 80 MWh
+        max_soc_fraction: 0.9  # percent as decimal
+        min_soc_fraction: 0.1  # percent as decimal
+        init_soc_fraction: 0.9  # percent as decimal
+        demand_profile: !include demand_profiles/demand_profile_secondary.yaml
+        charge_efficiency: 1.0  # percent as decimal
+        discharge_efficiency: 1.0  # percent as decimal
+      control_parameters:
+        max_discharge_rate: 2500.0  # kW/time step
+        charge_equals_discharge: true
+        demand_profile_supervisor: !include demand_profiles/demand_profile_supervisor.yaml
+        dispatch_priority_demand_profile: demand_profile_supervisor
+        max_supervisor_events: 2
+        max_supervisor_event_period: W
+        peak_range:
+          start: 12:00:00
+          end: 17:00:00
+        advance_discharge_period:
+          units: h
+          val: 2
+        delay_charge_period:
+          units: h
+          val: 4
+        allow_charge_in_peak_range: false
+        min_peak_proximity:
+          units: h
+          val: 4
+      cost_parameters:
+        cost_year: 2024
+        energy_capex: 408 # $/kWh  # conservative case for 2024 ATB utility scale batteries
+        power_capex: 379 # $/kW  # conservative case for 2024 ATB utility scale batteries
+        opex_fraction: 0.025  # 2.5% percent of capex as per 2024 ATB for utility scale batteries
+  grid_buy:
+    performance_model:
+      model: GridPerformanceModel
+    cost_model:
+      model: GridCostModel
+    model_inputs:
+      shared_parameters:
+        interconnection_size: 100000
+      cost_parameters:
+        cost_year: 2024
+        fixed_interconnection_cost: 0.0
+        interconnection_capex_per_kw: 0.0
+        interconnection_opex_per_kw: 0.0
+        electricity_buy_price: 0.09

h2integrate/control/control_strategies/storage/openloop_storage_control_base.py

@@ -12,17 +12,18 @@ class StorageOpenLoopControlBaseConfig(BaseConfig):
      Attributes:
         commodity (str): Name of the commodity being stored (e.g., "hydrogen").
         commodity_rate_units (str): Rate units of the commodity (e.g., "kg/h" or "kW").
-        demand_profile (int | float | list): Demand values for each timestep, in
+        demand_profile (int | float | list | dict): Demand values for each timestep, in
             the same units as `commodity_rate_units`. May be a scalar for constant
-            demand or a list/array for time-varying demand.
+            demand or a list/array/dict for time-varying demand. If a dict is provided, it
+            it should have two keys: "time_date" and "demand".
         commodity_amount_units (str | None, optional): Units of the commodity as an amount
             (i.e., kW*h or kg). If not provided, defaults to `commodity_rate_units*h`.
 
     """
 
     commodity: str = field()
     commodity_rate_units: str = field()
-    demand_profile: int | float | list = field()
+    demand_profile: int | float | list | dict = field()
     commodity_amount_units: str = field(default=None)
 
     def __attrs_post_init__(self):
@@ -48,9 +49,11 @@ def setup(self):
 
         commodity = self.config.commodity
 
+        demand_data = self.config.demand_profile
+
         self.add_input(
             f"{commodity}_demand",
-            val=self.config.demand_profile,
+            val=demand_data if not isinstance(demand_data, dict) else demand_data["demand"],
             shape=self.n_timesteps,
             units=self.config.commodity_rate_units,
             desc=f"Demand profile of {commodity}",