Skip to content

Commit

Permalink
Merge branch 'master' into compathelper/new_version/2020-04-14-08-05-…
Browse files Browse the repository at this point in the history 
…29-584-1772495339
  • Loading branch information
@hakanergun
hakanergun authored Apr 14, 2020
2 parents b5f9d70 + 4b14844 commit c7ea451
Show file tree
Hide file tree
Showing 10 changed files with 247 additions and 233 deletions.
2 changes: 2 additions & 0 deletions .travis.yml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -26,6 +26,8 @@ before_install:
- julia -e 'using Pkg; Pkg.add("Ipopt")'
- julia -e 'using Pkg; Pkg.add("SCS")'
- julia -e 'using Pkg; Pkg.add("Documenter")'
- julia -e 'using Pkg; Pkg.add("Juniper")'
- julia -e 'using Pkg; Pkg.add("Cbc")'

script:
- if [[ -a .git/shallow ]]; then git fetch --unshallow; fi
Expand Down
8 changes: 3 additions & 5 deletions Project.toml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,15 +16,13 @@ InfrastructureModels = "~ 0.4"
Ipopt = "≥ 0.4"
JuMP = "~0.19.1, ~0.20, ~0.21"
Memento = "~0.10, ~0.11, ~0.12, ~0.13, ~1.0"
PowerModels = "~ 0.15.3, 0.16"
PowerModels = "~ 0.15.3"
SCS = "≥ 0.6.3"
julia = "^1"

[extras]
Cbc = "9961bab8-2fa3-5c5a-9d89-47fab24efd76"
Ipopt = "b6b21f68-93f8-5de0-b562-5493be1d77c9"
SCS = "c946c3f1-0d1f-5ce8-9dea-7daa1f7e2d13"
PowerModels = "c36e90e8-916a-50a6-bd94-075b64ef4655"
Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

[targets]
test = ["SCS", "Test", "Cbc", "Ipopt"]
test = ["Test", "PowerModels"]
28 changes: 25 additions & 3 deletions README.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -43,8 +43,8 @@ The developers thank Carleton Coffrin (LANL) for his support.

## Citing PowerModelsACDC

If you find PowerModelsACDC useful in your work, we kindly request that you cite the following [publication](https://ieeexplore.ieee.org/document/8636236):

If you find PowerModelsACDC useful in your work, we kindly request that you cite the following publications:
[AC/DC OPF Core](https://ieeexplore.ieee.org/document/8636236):
```
@ARTICLE{8636236,
author={H. {Ergun} and J. {Dave} and D. {Van Hertem} and F. {Geth}},
Expand All @@ -59,7 +59,29 @@ doi={10.1109/TPWRS.2019.2897835},
ISSN={0885-8950},
month={July},}
```

[TNEP Extension 1](https://digital-library.theiet.org/content/journals/10.1049/iet-gtd.2019.0383):
```
@ARTICLE{
iet:/content/journals/10.1049/iet-gtd.2019.0383,
author = {Jay Dave},
author = {Hakan Ergun},
author = {Ting An},
author = {Jingjing Lu},
author = {Dirk Van Hertem},
keywords = {power systems;meshed HVDC grids;increased utilisation;presented formulations;convex formulations;second-order cone convex relaxation;multiple HVDC links;linear approximation;dc grids;transmission network expansion planning problem;high-voltage direct current;traditional ac grid;TNEP problem;nonlinear formulation;},
ISSN = {1751-8687},
title = {TNEP of meshed HVDC grids: ‘AC’, ‘DC’ and convex formulations},
journal = {IET Generation, Transmission & Distribution},
issue = {24},
volume = {13},
year = {2019},
month = {December},
pages = {5523-5532(9)},
publisher ={Institution of Engineering and Technology},
copyright = {© The Institution of Engineering and Technology},
url = {https://digital-library.theiet.org/content/journals/10.1049/iet-gtd.2019.0383}
}
```
## License

This code is provided under a BSD license.
8 changes: 8 additions & 0 deletions ReleaseNotes.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,8 @@
v0.3.0 Release Notes

- Added TNEP problems in single and multiperiod
- TNEP formulated in DCP, LPAC, SOCBFM, SOCBIM, SDP, QC and ACP
- TNEP test cases and unit tests integrated
- Update code base for compliance with PowerModels v0.15.3
-- Variable definitions & reporting
-- Conductors removed
48 changes: 26 additions & 22 deletions src/core/variabledcgrid.jl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -189,6 +189,8 @@ function variable_dcgrid_voltage_magnitude_sqr_ne(pm::_PM.AbstractPowerModel; nw
JuMP.set_upper_bound(wdcr_ne[br], bus_vdcmax[bi_bp[br][1]] * bus_vdcmax[bi_bp[br][2]])
end
end
report && _PM.sol_component_value(pm, nw, :busdc_ne, :wdc_ne, _PM.ids(pm, nw, :busdc_ne), wdc_ne)
report && _PM.sol_component_value(pm, nw, :busdc_ne, :wdcr_ne, _PM.ids(pm, nw, :busdc_ne), wdcr_ne)
end

function variable_dcgrid_voltage_magnitude_sqr_du(pm::_PM.AbstractPowerModel; nw::Int=pm.cnw, bounded::Bool = true, report::Bool=true) # this has to to every branch, different than its counterpart(Wdc_fr) since two candidate branches can be connected to same node and two duplicate variables will be needed
Expand Down Expand Up @@ -217,43 +219,45 @@ function variable_dcgrid_voltage_magnitude_sqr_du(pm::_PM.AbstractPowerModel; nw
JuMP.set_upper_bound(wdcr_frto_ne[i], 1.21)
end
end

report && _PM.sol_component_value(pm, nw, :busdc_ne, :wdc_du_fr, _PM.ids(pm, nw, :busdc_ne), wdc_fr_ne)
report && _PM.sol_component_value(pm, nw, :busdc_ne, :wdc_du_to, _PM.ids(pm, nw, :busdc_ne), wdc_to_ne)
report && _PM.sol_component_value(pm, nw, :busdc_ne, :wdcr_du, _PM.ids(pm, nw, :busdc_ne), wdcr_frto_ne)
end

"variable: `p_dcgrid[l,i,j]` for `(l,i,j)` in `arcs_dcgrid`"
function variable_active_dcbranch_flow_ne(pm::_PM.AbstractPowerModel; nw::Int=pm.cnw, bounded::Bool = true, report::Bool=true)
if bounded
_PM.var(pm, nw)[:p_dcgrid_ne]= JuMP.@variable(pm.model,
[(l,i,j) in _PM.ref(pm, nw, :arcs_dcgrid_ne)], base_name="$(nw)_pdcgrid_ne",
lower_bound = -_PM.ref(pm, nw, :branchdc_ne, l, "rateA"),
upper_bound = _PM.ref(pm, nw, :branchdc_ne, l, "rateA"),
start = _PM.comp_start_value(_PM.ref(pm, nw, :branchdc_ne, l), "p_start", 1.0)
)
else
_PM.var(pm, nw)[:p_dcgrid_ne] = JuMP.@variable(pm.model,
p = _PM.var(pm, nw)[:p_dcgrid_ne]= JuMP.@variable(pm.model,
[(l,i,j) in _PM.ref(pm, nw, :arcs_dcgrid_ne)], base_name="$(nw)_pdcgrid_ne",
start = _PM.comp_start_value(_PM.ref(pm, nw, :branchdc_ne, l), "p_start", 1.0)
)

if bounded
for arc in _PM.ref(pm, nw, :arcs_dcgrid_ne)
l,i,j = arc
JuMP.set_lower_bound(p[arc], -_PM.ref(pm, nw, :branchdc_ne, l)["rateA"])
JuMP.set_upper_bound(p[arc], _PM.ref(pm, nw, :branchdc_ne, l)["rateA"])
end
end

report && _PM.sol_component_value_edge(pm, nw, :branchdc_ne, :pf, :pt, _PM.ref(pm, nw, :arcs_dcgrid_from_ne), _PM.ref(pm, nw, :arcs_dcgrid_to_ne), p)
end

"variable: `ccm_dcgrid[l]` for `(l)` in `branchdc`"
function variable_dcbranch_current_sqr_ne(pm::_PM.AbstractPowerModel; nw::Int=pm.cnw, bounded::Bool = true, report::Bool=true)
vpu = 0.8
cc= _PM.var(pm, nw)[:ccm_dcgrid_ne] = JuMP.@variable(pm.model,
[l in _PM.ids(pm, nw, :branchdc_ne)], base_name="$(nw)_ccm_dcgrid_ne",
start = (_PM.comp_start_value(_PM.ref(pm, nw, :branchdc_ne, l), "p_start", 0.0) / vpu)^2
)

if bounded
_PM.var(pm, nw)[:ccm_dcgrid_ne] = JuMP.@variable(pm.model,
[l in _PM.ids(pm, nw, :branchdc_ne)], base_name="$(nw)_ccm_dcgrid_ne",
lower_bound = 0,
upper_bound = (_PM.ref(pm, nw, :branchdc_ne, l, "rateA")/vpu)^2,
start = (_PM.comp_start_value(_PM.ref(pm, nw, :branchdc_ne, l), "p_start", 0.0) / vpu)^2
)
else
_PM.var(pm, nw)[:ccm_dcgrid_ne] = JuMP.@variable(pm.model,
[l in _PM.ids(pm, nw, :branchdc_ne)], base_name="$(nw)_ccm_dcgrid_ne",
start = (_PM.comp_start_value(_PM.ref(pm, nw, :branchdc_ne, l), "p_start", 0.0)/vpu)^2,
lower_bound = 0,
)
for (l, branchdc) in _PM.ref(pm, nw, :branchdc_ne)
JuMP.set_lower_bound(cc[l], 0)
JuMP.set_upper_bound(cc[l], (branchdc["rateA"] / vpu)^2)
end
end

report && _PM.sol_component_value(pm, nw, :branchdc_ne, :ccm, _PM.ids(pm, nw, :branchdc), cc)
end

"variable: `0 <= convdc_ne[c] <= 1` for `c` in `candidate converters"
Expand Down
4 changes: 1 addition & 3 deletions test/REQUIRE
View file
Original file line number Diff line number Diff line change
@@ -1,6 +1,4 @@
julia 1.1


Ipopt 0.4
SCS 0.4
Cbc 0.4
PowerModels 0.15.3
22 changes: 3 additions & 19 deletions test/data/tnep/case6_test.m
View file
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
% Case to test two connected components in the network data
% the case is two replicates of the case3 network

% Garver system - Transmission Network Estimation Using Linear Programming, IEEE trans. on power appratus and sys
% Power system transmission network expansion planning using AC model by Rider, MJ and Garcia, AV and Romero, R
%modification: gen cost is changed. A low value is set to penalize power losses but not dominate the objective.
function mpc = case6
mpc.version = '2';
mpc.baseMVA = 100.0;
Expand Down Expand Up @@ -57,31 +57,15 @@
%% candidate branches
%column_names% fbusdc tbusdc r l c rateA rateB rateC status cost
mpc.branchdc_ne = [
1 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.1;
1 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.1;
1 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.1;
2 3 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.2;
2 3 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.2;
2 3 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.2;
2 5 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.3;
2 5 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.3;
2 5 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.3;
2 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.4;
2 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.4;
2 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.4;
3 4 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.5;
3 4 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.5;
3 4 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.5;
3 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.6;
3 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.6;
3 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.6;
4 5 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.7;
4 5 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.7;
4 5 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.7;
4 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.8;
4 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.8;
4 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.8;
5 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.9;
5 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 2.9;
5 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 3.9;
];
Expand Down
136 changes: 21 additions & 115 deletions test/data/tnep/case9_test.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -12,35 +12,35 @@
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;

%modified case to motivate expansion: L45 deleted and L56s rating reduced. gen cost = 0 to check only expansion cost
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 3 0 0 0 0 1 1 0 345 1 1.1 0.9;
2 2 0 0 0 0 1 1 0 345 1 1.1 0.9;
3 2 0 0 0 0 1 1 0 345 1 1.1 0.9;
4 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
5 1 270 90 0 0 1 1 0 345 1 1.1 0.9;
5 1 90 30 0 0 1 1 0 345 1 1.1 0.9;
6 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
7 1 300 105 0 0 1 1 0 345 1 1.1 0.9;
7 1 100 35 0 0 1 1 0 345 1 1.1 0.9;
8 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
9 1 375 150 0 0 1 1 0 345 1 1.1 0.9;
9 1 125 50 0 0 1 1 0 345 1 1.1 0.9;
];

%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf
mpc.gen = [
1 0 0 900 900 1 100 1 750 10 0 0 0 0 0 0 0 0 0 0 0;
2 163 0 900 900 1 100 1 900 10 0 0 0 0 0 0 0 0 0 0 0;
3 85 0 900 900 1 100 1 810 10 0 0 0 0 0 0 0 0 0 0 0;
1 0 0 300 -300 1 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0;
2 163 0 300 -300 1 100 1 300 10 0 0 0 0 0 0 0 0 0 0 0;
3 85 0 300 -300 1 100 1 270 10 0 0 0 0 0 0 0 0 0 0 0;
];

%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
1 4 0 0.0576 0 250 250 250 0 0 1 -360 360;
4 5 0.017 0.092 0.158 250 250 250 0 0 1 -360 360;
5 6 0.039 0.17 0.358 150 150 150 0 0 1 -360 360;
% 4 5 0.017 0.092 0.158 50 250 250 0 0 1 -360 360;
5 6 0.039 0.17 0.358 50 150 150 0 0 1 -360 360;
3 6 0 0.0586 0 300 300 300 0 0 1 -360 360;
6 7 0.0119 0.1008 0.209 150 150 150 0 0 1 -360 360;
7 8 0.0085 0.072 0.149 250 250 250 0 0 1 -360 360;
Expand All @@ -54,124 +54,30 @@
% 1 startup shutdown n x1 y1 ... xn yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 1500 0 3 0 0.01 0; % modified
2 1500 0 3 0 0.01 0; % modified
2 1500 0 3 0 0.01 0; % modified
2 1500 0 3 0 0.0 0; % modified
2 1500 0 3 0 0.0 0; % modified
2 1500 0 3 0 0.0 0; % modified
];


%% candidate dc bus data
%column_names% busdc_i grid Pdc Vdc basekVdc Vdcmax Vdcmin Cdc
mpc.busdc_ne = [
1 1 0 1 345 1.1 0.9 0;
2 1 0 1 345 1.1 0.9 0;
3 1 0 1 345 1.1 0.9 0;
4 1 0 1 345 1.1 0.9 0;
5 1 0 1 345 1.1 0.9 0;
6 1 0 1 345 1.1 0.9 0;
7 1 0 1 345 1.1 0.9 0;
8 1 0 1 345 1.1 0.9 0;
9 1 0 1 345 1.1 0.9 0;
1 1 0 1 345 1.1 0.9 0;
2 1 0 1 345 1.1 0.9 0;
];


%% candidate branches
%column_names% fbusdc tbusdc r l c rateA rateB rateC status cost
mpc.branchdc_ne = [
1 2 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.1;
1 3 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.2;
1 5 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.3;
1 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.4;
1 7 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.5;
1 8 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.6;
1 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.7;
2 3 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.8;
2 4 0.01 0.00 0.00 100.0 0.0 0.0 1.0 1.9;
2 5 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.1;
2 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.2;
2 7 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.3;
2 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.4;
3 4 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.5;
3 5 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.6;
3 7 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.7;
3 8 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.8;
3 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 2.9;
4 6 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3;
4 7 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.1;
4 8 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.2;
5 7 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.3;
5 8 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.4;
5 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.5;
6 8 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.6;
6 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.7;
7 9 0.01 0.00 0.00 100.0 0.0 0.0 1.0 3.8;

1 2 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.1;
1 3 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.2;
1 5 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.3;
1 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.4;
1 7 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.5;
1 8 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.6;
1 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.7;
2 3 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.8;
2 4 0.01 0.00 0.00 200.0 0.0 0.0 1.0 4.9;
2 5 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.0;
2 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.1;
2 7 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.2;
2 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.3;
3 4 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.4;
3 5 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.5;
3 7 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.6;
3 8 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.7;
3 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.8;
4 6 0.01 0.00 0.00 200.0 0.0 0.0 1.0 5.9;
4 7 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.0;
4 8 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.1;
5 7 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.2;
5 8 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.3;
5 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.4;
6 8 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.5;
6 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.6;
7 9 0.01 0.00 0.00 200.0 0.0 0.0 1.0 6.7;

1 2 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.1;
1 3 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.2;
1 5 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.3;
1 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.4;
1 7 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.5;
1 8 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.6;
1 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.7;
2 3 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.8;
2 4 0.01 0.00 0.00 400.0 0.0 0.0 1.0 7.9;
2 5 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8;
2 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.1;
2 7 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.2;
2 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.3;
3 4 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.4;
3 5 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.5;
3 7 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.6;
3 8 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.7;
3 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.8;
4 6 0.01 0.00 0.00 400.0 0.0 0.0 1.0 8.9;
4 7 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9;
4 8 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.1;
5 7 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.2;
5 8 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.3;
5 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.4;
6 8 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.5;
6 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.6;
7 9 0.01 0.00 0.00 400.0 0.0 0.0 1.0 9.7;
];
1 2 0.017 0.00 0.00 250 0.0 0.0 1.0 1.2;
1 2 0.017 0.00 0.00 250 0.0 0.0 1.0 1.3;
];

%% candidate converters
%column_names% busdc_i busac_i type_dc type_ac P_g Q_g islcc Vtar rtf xtf transformer tm bf filter rc xc reactor basekVac Vmmax Vmmin Imax status LossA LossB LossCrec LossCinv droop Pdcset Vdcset dVdcset Pacmax Pacmin Qacmax Qacmin cost
mpc.convdc_ne = [
1 1 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.1;
2 2 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.2;
3 3 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.3;
4 4 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.4;
5 5 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.5;
6 6 2 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.6;
7 7 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.7;
8 8 1 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.8;
9 9 2 1 -360 -1.66 0 1.0 0.001 0.001 1 1 0.001 1 0.001 0.001 1 345 1.1 0.9 10 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 700 -700 700 -700 3.9;
mpc.convdc_ne = [
1 4 1 1 -360 -1.66 0 1.0 0.001 0.1 1 1 0.08 1 0.001 0.09 1 345 1.1 0.9 100 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 1000 -1000 500 -500 4.5;
2 5 1 1 -360 -1.66 0 1.0 0.001 0.1 1 1 0.08 1 0.001 0.09 1 345 1.1 0.9 100 1 1.1033 0.887 2.885 2.885 0.0050 -52.7 1.0079 0 1000 -1000 500 -500 5;
];
Loading

0 comments on commit c7ea451

Please sign in to comment.