Jd/add thermal time series

Project.toml

@@ -1,7 +1,7 @@
 name = "PowerSimulations"
 uuid = "e690365d-45e2-57bb-ac84-44ba829e73c4"
 authors = ["Jose Daniel Lara", "Clayton Barrows", "Daniel Thom", "Dheepak Krishnamurthy", "Sourabh Dalvi"]
-version = "0.28.3"
+version = "0.29.0"
 
 [deps]
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"

README.md

@@ -43,6 +43,8 @@ using PowerSimulations
 using PowerSystems
 ```
 
+For information on using the package, see the [stable documentation](https://nrel-sienna.github.io/PowerSimulations.jl/stable/). Use the [in-development documentation](https://nrel-sienna.github.io/PowerSimulations.jl/dev/) for the version of the documentation which contains the unreleased features.
+
 ## Development
 
 Contributions to the development and enhancement of PowerSimulations is welcome. Please see [CONTRIBUTING.md](https://github.com/NREL-Sienna/PowerSimulations.jl/blob/main/CONTRIBUTING.md) for code contribution guidelines.

src/PowerSimulations.jl

@@ -556,7 +556,7 @@ include("devices_models/devices/common/duration_constraints.jl")
 include("devices_models/devices/common/get_time_series.jl")
 
 # Device Modeling components
-include("devices_models/devices/interfaces.jl")
+include("devices_models/devices/default_interface_methods.jl")
 include("devices_models/devices/common/add_to_expression.jl")
 include("devices_models/devices/common/set_expression.jl")
 include("devices_models/devices/renewable_generation.jl")

src/devices_models/device_constructors/thermalgeneration_constructor.jl

@@ -56,6 +56,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, D())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -135,6 +139,17 @@ function construct_device!(
     add_constraints!(container, RampConstraint, devices, model, network_model)
     add_constraints!(container, DurationConstraint, devices, model, network_model)
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
+
     add_feedforward_constraints!(container, model, devices)
 
     objective_function!(container, devices, model, get_network_formulation(network_model))
@@ -164,6 +179,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, D())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -227,6 +246,16 @@ function construct_device!(
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
     add_constraints!(container, RampConstraint, devices, model, network_model)
     add_constraints!(container, DurationConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -256,6 +285,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalBasicUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -335,6 +368,17 @@ function construct_device!(
     )
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
+
     add_feedforward_constraints!(container, model, devices)
 
     objective_function!(container, devices, model, get_network_formulation(network_model))
@@ -361,6 +405,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalBasicUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -423,6 +471,16 @@ function construct_device!(
     )
 
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -829,6 +887,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalMultiStartUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -931,6 +993,16 @@ function construct_device!(
         network_model,
     )
     add_constraints!(container, ActiveRangeICConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -965,6 +1037,10 @@ function construct_device!(
     add_variables!(container, TimeDurationOff, devices, ThermalMultiStartUnitCommitment())
     add_variables!(container, PowerOutput, devices, ThermalMultiStartUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -1053,6 +1129,16 @@ function construct_device!(
         network_model,
     )
     add_constraints!(container, ActiveRangeICConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -1092,6 +1178,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalCompactUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -1169,6 +1259,16 @@ function construct_device!(
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
     add_constraints!(container, RampConstraint, devices, model, network_model)
     add_constraints!(container, DurationConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -1202,6 +1302,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalCompactUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -1270,6 +1374,16 @@ function construct_device!(
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
     add_constraints!(container, RampConstraint, devices, model, network_model)
     add_constraints!(container, DurationConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -1307,6 +1421,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalBasicCompactUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -1382,6 +1500,16 @@ function construct_device!(
         network_model,
     )
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 
@@ -1413,6 +1541,10 @@ function construct_device!(
 
     initial_conditions!(container, devices, ThermalBasicCompactUnitCommitment())
 
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_parameters!(container, ActivePowerTimeSeriesParameter, devices, model)
+    end
+
     add_to_expression!(
         container,
         ActivePowerBalance,
@@ -1479,6 +1611,16 @@ function construct_device!(
     )
 
     add_constraints!(container, CommitmentConstraint, devices, model, network_model)
+    if haskey(get_time_series_names(model), ActivePowerTimeSeriesParameter)
+        add_constraints!(
+            container,
+            ActivePowerVariableTimeSeriesLimitsConstraint,
+            ActivePowerRangeExpressionUB,
+            devices,
+            model,
+            network_model,
+        )
+    end
 
     add_feedforward_constraints!(container, model, devices)
 

src/devices_models/devices/common/add_to_expression.jl

@@ -45,7 +45,7 @@ function add_expressions!(
         D,
         PSY.get_name.(devices),
         time_steps;
-        meta = PSY.get_name(first(devices)),
+        meta = get_service_name(model),
     )
     return
 end

src/devices_models/devices/common/objective_function/common.jl

@@ -88,8 +88,14 @@ function _add_vom_cost_to_objective!(
     )
     for t in get_time_steps(container)
         exp =
-            _add_proportional_term!(container, T(), d, cost_term_normalized * multiplier, t)
-        add_to_expression!(container, ProductionCostExpression, exp, d, t)
+            _add_proportional_term!(
+                container,
+                T(),
+                component,
+                cost_term_normalized * multiplier,
+                t,
+            )
+        add_to_expression!(container, ProductionCostExpression, exp, component, t)
     end
     return
 end

src/devices_models/devices/common/range_constraint.jl

@@ -705,7 +705,12 @@ function lower_bound_range_with_parameter!(
     mult = get_multiplier_array(param_container)
     jump_model = get_jump_model(container)
     time_steps = axes(constraint_container)[2]
+    ts_name = get_time_series_names(model)[P]
+    ts_type = get_default_time_series_type(container)
     for device in devices
+        if !(PSY.has_time_series(device, ts_type, ts_name))
+            continue
+        end
         name = PSY.get_name(device)
         param = get_parameter_column_refs(param_container, name)
         for t in time_steps
@@ -730,8 +735,13 @@ function _add_parameterized_lower_bound_range_constraints_impl!(
     W <: AbstractDeviceFormulation,
 }
     time_steps = get_time_steps(container)
-    names = [PSY.get_name(d) for d in devices]
-
+    ts_name = get_time_series_names(model)[U]
+    ts_type = get_default_time_series_type(container)
+    names = [PSY.get_name(d) for d in devices if PSY.has_time_series(d, ts_type, ts_name)]
+    if isempty(names)
+        @debug "There are no $V devices with time series data"
+        return
+    end
     constraint =
         add_constraints_container!(container, T(), V, names, time_steps; meta = "lb")
 
@@ -834,14 +844,19 @@ function upper_bound_range_with_parameter!(
     lhs_array,
     param::P,
     devices::IS.FlattenIteratorWrapper{V},
-    ::DeviceModel{V, W},
+    model::DeviceModel{V, W},
 ) where {P <: TimeSeriesParameter, V <: PSY.Component, W <: AbstractDeviceFormulation}
     param_container = get_parameter(container, param, V)
     mult = get_multiplier_array(param_container)
     jump_model = get_jump_model(container)
     time_steps = axes(constraint_container)[2]
+    ts_name = get_time_series_names(model)[P]
+    ts_type = get_default_time_series_type(container)
     for device in devices
         name = PSY.get_name(device)
+        if !(PSY.has_time_series(device, ts_type, ts_name))
+            continue
+        end
         param = get_parameter_column_refs(param_container, name)
         for t in time_steps
             constraint_container[name, t] =
@@ -865,7 +880,13 @@ function _add_parameterized_upper_bound_range_constraints_impl!(
     W <: AbstractDeviceFormulation,
 }
     time_steps = get_time_steps(container)
-    names = [PSY.get_name(d) for d in devices]
+    ts_name = get_time_series_names(model)[P]
+    ts_type = get_default_time_series_type(container)
+    names = [PSY.get_name(d) for d in devices if PSY.has_time_series(d, ts_type, ts_name)]
+    if isempty(names)
+        @debug "There are no $V devices with time series data $ts_type, $ts_name"
+        return
+    end
 
     constraint =
         add_constraints_container!(container, T(), V, names, time_steps; meta = "ub")