ANL-CEEESA/UnitCommitment.jl
1
0
Fork 0
mirror of https://github.com/ANL-CEEESA/UnitCommitment.jl.git synced 2025-12-06 00:08:52 -06:00
Code Issues Packages Projects Releases Wiki Activity

enforce_transmission: Add constraint refs to model

Browse Source
This commit is contained in:
Alinson S. Xavier
2024-08-05 20:03:38 -05:00
parent bb880575ef
commit 7897adbee1
1 changed files with 34 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

60
src/solution/methods/XavQiuWanThi2019/enforce.jl
View File

@@ -26,59 +26,67 @@ function _enforce_transmission(;
isf::Matrix{Float64},
lodf::Matrix{Float64},
)::Nothing
instance = model[:instance]
limit::Float64 = 0.0
overflow = model[:overflow]
net_injection = model[:net_injection]
lm = violation.monitored_line
lc = violation.outage_line
t = violation.time
eq_flow_ub = _init(model, :eq_flow_ub)
eq_flow_lb = _init(model, :eq_flow_lb)
eq_flow_def = _init(model, :eq_flow_def)
eq_idx = (
sc.name,
lm.name,
lc === nothing ? "Base" : lc.name,
t,
)
if violation.outage_line === nothing
limit = violation.monitored_line.normal_flow_limit[violation.time]
if lc === nothing
limit = lm.normal_flow_limit[t]
@info @sprintf(
" %8.3f MW overflow in %-5s time %3d (pre-contingency, scenario %s)",
violation.amount,
violation.monitored_line.name,
violation.time,
lm.name,
t,
sc.name,
)
else
limit = violation.monitored_line.emergency_flow_limit[violation.time]
limit = lm.emergency_flow_limit[t]
@info @sprintf(
" %8.3f MW overflow in %-5s time %3d (outage: line %s, scenario %s)",
violation.amount,
violation.monitored_line.name,
violation.time,
violation.outage_line.name,
lm.name,
t,
lc.name,
sc.name,
)
end
fm = violation.monitored_line.name
t = violation.time
flow = @variable(model, base_name = "flow[$fm,$t]")
v = overflow[sc.name, lm.name, t]
flow = @variable(model, base_name = "flow[$eq_idx]")
eq_flow_ub[eq_idx] = @constraint(model, flow <= limit + v)
eq_flow_lb[eq_idx] = @constraint(model, -flow <= limit + v)
v = overflow[sc.name, violation.monitored_line.name, violation.time]
@constraint(model, flow <= limit + v)
@constraint(model, -flow <= limit + v)
if violation.outage_line === nothing
@constraint(
if lc === nothing
eq_flow_def[eq_idx] = @constraint(
model,
flow == sum(
net_injection[sc.name, b.name, violation.time] *
isf[violation.monitored_line.offset, b.offset] for
net_injection[sc.name, b.name, t] *
isf[lm.offset, b.offset] for
b in sc.buses if b.offset > 0
)
)
else
@constraint(
eq_flow_def[eq_idx] = @constraint(
model,
flow == sum(
net_injection[sc.name, b.name, violation.time] * (
isf[violation.monitored_line.offset, b.offset] + (
net_injection[sc.name, b.name, t] * (
isf[lm.offset, b.offset] + (
lodf[
violation.monitored_line.offset,
violation.outage_line.offset,
] * isf[violation.outage_line.offset, b.offset]
lm.offset,
lc.offset,
] * isf[lc.offset, b.offset]
)
) for b in sc.buses if b.offset > 0
)