Skip to content

Update MT-HVDC quadratic converter model #103

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
acostarelli wants to merge 9 commits into
base: main
Choose a base branch
from
Open

Update MT-HVDC quadratic converter model #103

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
35145a8
claude hvdc
Apr 29, 2026
b39c63d
formatting
Apr 29, 2026
b745791
fix linear loss adding;fix hdf output
Apr 29, 2026
4c8c8fe
formatting
Apr 30, 2026
abbc792
move _voltage helpers to constructor
Apr 30, 2026
baa0288
Add HVDCTwoTerminalVSC two-terminal HVDC formulation (#119)
acostarelli May 18, 2026
f0a590a
Simplify HVDC converter loss modeling; consolidate HVDC tests
May 18, 2026
79200d9
Tidy quadratic converter loss helpers
May 19, 2026
046c508
use jump utils; fix docstring
May 19, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 6 additions & 6 deletions Project.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,18 +23,18 @@ TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
[weakdeps]
PowerFlows = "94fada2c-0ca5-4b90-a1fb-4bc5b59ccfc7"

[sources]
InfrastructureSystems = {rev = "IS4", url = "https://github.com/Sienna-Platform/InfrastructureSystems.jl"}
PowerSystems = {rev = "psy6", url = "https://github.com/Sienna-Platform/PowerSystems.jl"}
InfrastructureOptimizationModels = {rev = "main", url = "https://github.com/Sienna-Platform/InfrastructureOptimizationModels.jl"}

[extensions]
PowerFlowsExt = "PowerFlows"

[sources]
InfrastructureSystems = {url = "https://github.com/NREL-Sienna/InfrastructureSystems.jl", rev = "IS4"}
PowerSystems = {url = "https://github.com/NREL-Sienna/PowerSystems.jl", rev = "psy6"}
InfrastructureOptimizationModels = {url = "https://github.com/NREL-Sienna/InfrastructureOptimizationModels.jl", rev = "lk/pom-test-fixes"}

[compat]
Dates = "1"
DocStringExtensions = "~0.8, ~0.9"
InfrastructureOptimizationModels = "0.1"
InfrastructureOptimizationModels = "0.1.0"
InfrastructureSystems = "3"
InteractiveUtils = "1.11.0"
JuMP = "^1.28"
Expand Down
16 changes: 3 additions & 13 deletions src/PowerOperationsModels.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -505,24 +505,12 @@ export PostContingencyActivePowerReserveDeploymentVariable
# HVDC Variables
export DCVoltage
export DCLineCurrent
export ConverterPowerDirection
export ConverterCurrent
export SquaredConverterCurrent
export InterpolationSquaredCurrentVariable
export InterpolationBinarySquaredCurrentVariable
export ConverterPositiveCurrent
export ConverterNegativeCurrent
export SquaredDCVoltage
export InterpolationSquaredVoltageVariable
export InterpolationBinarySquaredVoltageVariable
export AuxBilinearConverterVariable
export AuxBilinearSquaredConverterVariable
export InterpolationSquaredBilinearVariable
export InterpolationBinarySquaredBilinearVariable
export ConverterCurrentDirection
export HVDCFlowDirectionVariable
export HVDCLosses
export ConverterDCPower
export ConverterCurrentDirection

Comment thread
acostarelli marked this conversation as resolved.
# Load Variables
export ShiftUpActivePowerVariable
Expand Down Expand Up @@ -774,6 +762,8 @@ export HVDCTwoTerminalLCC
# Converter Formulations
export LosslessConverter
export LinearLossConverter
export AbstractQuadraticLossConverter
export Bin2QuadraticLossConverter
export QuadraticLossConverter

# DC Line Formulations
Expand Down
73 changes: 0 additions & 73 deletions src/core/constraints.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -577,19 +577,6 @@ struct DCLineCurrentConstraint <: ConstraintType end

struct NodalBalanceCurrentConstraint <: ConstraintType end

"""
Struct to create the constraints that compute the converter DC power based on current and voltage.

The specified constraints are formulated as:
```math
\\begin{align*}
& p_c = 0.5 * (γ^sq - v^sq - i^sq), \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& γ_c = v_c + i_c, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
\\end{align*}
```
"""
struct ConverterPowerCalculationConstraint <: ConstraintType end

"""
Struct to create the constraints that decide the balance of AC and DC power of the converter.

Expand All @@ -603,66 +590,6 @@ The specified constraints are formulated as:
"""
struct ConverterLossConstraint <: ConstraintType end

"""
Struct to create the McCormick envelopes constraints that decide the bounds on the DC active power.

The specified constraints are formulated as:
```math
\\begin{align*}
& p_c >= V^{min} i_c + v_c I^{min} - I^{min}V^{min}, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& p_c >= V^{max} i_c + v_c I^{max} - I^{max}V^{max}, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& p_c <= V^{max} i_c + v_c I^{min} - I^{min}V^{max}, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& p_c <= V^{min} i_c + v_c I^{max} - I^{max}V^{min}, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
\\end{align*}
```
"""
struct ConverterMcCormickEnvelopes <: ConstraintType end

"""
Struct to create the Quadratic PWL interpolation constraints that decide square value of the voltage.
In this case x = voltage and y = squared_voltage.
The specified constraints are formulated as:
```math
\\begin{align*}
& x = x_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& y = y_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& z_k \\le \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
& z_k \\ge \\delta_{k+1}, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
\\end{align*}
```
"""
struct InterpolationVoltageConstraints <: ConstraintType end

"""
Struct to create the Quadratic PWL interpolation constraints that decide square value of the current.
In this case x = current and y = squared_current.
The specified constraints are formulated as:
```math
\\begin{align*}
& x = x_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& y = y_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& z_k \\le \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
& z_k \\ge \\delta_{k+1}, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
\\end{align*}
```
"""
struct InterpolationCurrentConstraints <: ConstraintType end

"""
Struct to create the Quadratic PWL interpolation constraints that decide square value of the bilinear variable γ.
In this case x = γ and y = squared_γ.
The specified constraints are formulated as:
```math
\\begin{align*}
& x = x_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& y = y_0 + \\sum_{k=1}^K (x_{k} - x_{k-1}) \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\} \\\\
& z_k \\le \\delta_k, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
& z_k \\ge \\delta_{k+1}, \\quad \\forall t \\in \\{1,\\dots, T\\}, \\forall k \\in \\{1,\\dots, K-1\\} \\\\
\\end{align*}
```
"""
struct InterpolationBilinearConstraints <: ConstraintType end

"""
Struct to create the constraints that set the absolute value for the current to use in losses through a lossy Interconnecting Power Converter.
The specified constraint is formulated as:
Expand Down
22 changes: 20 additions & 2 deletions src/core/formulations.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -201,9 +201,27 @@ Linear Loss InterconnectingConverter Model
struct LinearLossConverter <: AbstractConverterFormulation end

"""
Quadratic Loss InterconnectingConverter Model
Abstract supertype for InterconnectingConverter formulations with quadratic losses.
"""
struct QuadraticLossConverter <: AbstractConverterFormulation end
abstract type AbstractQuadraticLossConverter <: AbstractConverterFormulation end

"""
Quadratic Loss InterconnectingConverter using the Bin2 separable bilinear approximation
(`v·i = ½((v+i)² − v² − i²)`) with a SOS2-based PWL approximation for x².
"""
struct Bin2QuadraticLossConverter <: AbstractQuadraticLossConverter end

"""
Quadratic Loss InterconnectingConverter using exact bilinear (v·i) and quadratic (i²)
products. Requires an NLP-capable solver (e.g., Ipopt).

Attributes (set on `DeviceModel`):
- `use_linear_loss`: Include the linear loss term `b·(i⁺ + i⁻)` (default false). The
linear term requires an absolute-value reformulation that introduces a binary
variable, so this is only solvable by MINLP solvers (e.g., Gurobi); pure NLP
solvers like Ipopt cannot solve the model with this option enabled.
"""
struct QuadraticLossConverter <: AbstractQuadraticLossConverter end

############################## HVDC Lines Formulations ##################################
abstract type AbstractDCLineFormulation <: AbstractBranchFormulation end
Expand Down
79 changes: 0 additions & 79 deletions src/core/variables.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -150,30 +150,12 @@ Docs abbreviation: ``i_l^{dc}``
"""
struct DCLineCurrent <: VariableType end

"""
Struct to dispatch the creation of Squared Voltage Variables for DC formulations
Docs abbreviation: ``v^{sq,dc}``
"""
struct SquaredDCVoltage <: VariableType end

"""
Struct to dispatch the creation of DC Converter Current Variables for DC formulations
Docs abbreviation: ``i_c^{dc}``
"""
struct ConverterCurrent <: VariableType end

"""
Struct to dispatch the creation of DC Converter Power Variables for DC formulations
Docs abbreviation: ``p_c^{dc}``
"""
struct ConverterDCPower <: VariableType end

"""
Struct to dispatch the creation of Squared DC Converter Current Variables for DC formulations
Docs abbreviation: ``i_c^{sq,dc}``
"""
struct SquaredConverterCurrent <: VariableType end

"""
Struct to dispatch the creation of DC Converter Positive Term Current Variables for DC formulations
Docs abbreviation: ``i_c^{+,dc}``
Expand All @@ -192,67 +174,6 @@ Docs abbreviation: `\\nu_c``
"""
struct ConverterCurrentDirection <: VariableType end

"""
Struct to dispatch the creation of Binary Variable for Converter Power Direction
Docs abbreviation: ``\\kappa_c^{dc}``
"""
struct ConverterPowerDirection <: VariableType end

"""
Struct to dispatch the creation of Auxiliary Variable for Converter Bilinear term: v * i
Docs abbreviation: ``\\gamma_c^{dc}``
"""
struct AuxBilinearConverterVariable <: VariableType end

"""
Struct to dispatch the creation of Auxiliary Variable for Squared Converter Bilinear term: v * i

Docs abbreviation: ``\\gamma_c^{sq,dc}``
"""
struct AuxBilinearSquaredConverterVariable <: VariableType end

"""
Struct to dispatch the creation of Continuous Interpolation Variable for Squared Converter Voltage

Docs abbreviation: ``\\delta_c^{v}``
"""
struct InterpolationSquaredVoltageVariable <: InterpolationVariableType end

"""
Struct to dispatch the creation of Binary Interpolation Variable for Squared Converter Voltage

Docs abbreviation: ``z_c^{v}``
"""
struct InterpolationBinarySquaredVoltageVariable <: BinaryInterpolationVariableType end

"""
Struct to dispatch the creation of Continuous Interpolation Variable for Squared Converter Current

Docs abbreviation: ``\\delta_c^{i}``
"""
struct InterpolationSquaredCurrentVariable <: InterpolationVariableType end

"""
Struct to dispatch the creation of Binary Interpolation Variable for Squared Converter Current

Docs abbreviation: ``z_c^{i}``
"""
struct InterpolationBinarySquaredCurrentVariable <: BinaryInterpolationVariableType end

"""
Struct to dispatch the creation of Continuous Interpolation Variable for Squared Converter AuxVar

Docs abbreviation: ``\\delta_c^{\\gamma}``
"""
struct InterpolationSquaredBilinearVariable <: InterpolationVariableType end

"""
Struct to dispatch the creation of Binary Interpolation Variable for Squared Converter AuxVar

Docs abbreviation: ``z_c^{\\gamma}``
"""
struct InterpolationBinarySquaredBilinearVariable <: BinaryInterpolationVariableType end

#########################################################
#########################################################

Expand Down
Loading
Loading