├── CHANGELOG.md ├── LICENSE ├── MODEL.pdf ├── MODEL.tex ├── README.md ├── case24_7_jb.m ├── case3120_5_he.m ├── case39_10_he.m ├── case5_3_he.m ├── case67.m └── nem_2000bus_hvdc.m /CHANGELOG.md: -------------------------------------------------------------------------------- 1 | PGLib OPF HVDC Benchmarks Change Log 2 | ==================================== 3 | ### v23.09 4 | - Two new test cases added: case67.m, nem_2000bus_hvdc.m 5 | - Version number update for all other test cases, further no change 6 | ### v19.08 7 | - Initial release 8 | -------------------------------------------------------------------------------- /LICENSE: -------------------------------------------------------------------------------- 1 | Data License: 2 | 3 | Creative Commons Attribution 4.0 International license 4 | http://creativecommons.org/licenses/by/4.0/ 5 | 6 | See data files for specific copyright holders and attribution details. 7 | 8 | 9 | Software License: 10 | 11 | MIT License 12 | Copyright (c) 2019 A Library of IEEE PES Power Grid Benchmarks 13 | 14 | Permission is hereby granted, free of charge, to any person obtaining a copy 15 | of this software and associated documentation files (the "Software"), to deal 16 | in the Software without restriction, including without limitation the rights 17 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 18 | copies of the Software, and to permit persons to whom the Software is 19 | furnished to do so, subject to the following conditions: 20 | 21 | The above copyright notice and this permission notice shall be included in all 22 | copies or substantial portions of the Software. 23 | 24 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 25 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 26 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 27 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 28 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 29 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 30 | SOFTWARE. 31 | 32 | -------------------------------------------------------------------------------- /MODEL.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/power-grid-lib/pglib-opf-hvdc/f928f95cde06d035554d92854f03950e074850e0/MODEL.pdf -------------------------------------------------------------------------------- /MODEL.tex: -------------------------------------------------------------------------------- 1 | % render with latexit and then export as png, default settings 2 | % 18.00 pt (std display) 3 | % 9.00 pt (retina display) 4 | \documentclass{article} 5 | \usepackage{amssymb, amsmath, bm} 6 | \begin{document} 7 | 8 | 9 | 10 | % render with latexit and then export as png, default settings 11 | % 18.00 pt (std display) 12 | % 9.00 pt (retina display) 13 | % \usepackage{amssymb, amsmath, bm} 14 | \begin{align} 15 | % 16 | \mbox{\bf data:} & \nonumber \\ 17 | & \bm {S^{gl}}_k, \bm {S^{gu}}_k \;\; \forall k \in G \nonumber \\ 18 | & \bm {S^{cl}}_l, \bm {S^{cu}}_l \;\; \forall l \in C \nonumber \\ 19 | & \bm c_{2k}, \bm c_{1k}, \bm c_{0k} \;\; \forall k \in G \nonumber \\ 20 | & \bm a_l, \bm b_l, \bm c_l \;\; \forall l \in C \nonumber \\ 21 | & \bm {v^l}_i, \bm {v^u}_i \;\; \forall i \in N \nonumber \\ 22 | & \bm {v^l}_i, \bm {v^u}_i \;\; \forall i \in N_{dc} \nonumber \\ 23 | & \bm S^d_i, \bm Y^s_{i} \;\; \forall i \in N \nonumber \\ 24 | & \bm P^d_i, \bm Y^s_{i} \;\; \forall i \in N_{dc} \nonumber \\ 25 | & \bm Y_{ij}, \bm {b^c}_{ij}, \bm{T}_{ij} \;\; \forall (i,j) \in E \nonumber \\ 26 | & \bm Y_{ij} \;\; \forall (i,j) \in E_{dc} \nonumber \\ 27 | & \bm {s^u}_{ij}, \bm {\theta^{\Delta l}}_{ij}, \bm {\theta^{\Delta u}}_{ij} \;\; \forall (i,j) \in E \nonumber \\ 28 | & \bm {p^u}_{ij} \;\; \forall (i,j) \in E_{dc} \nonumber \\ 29 | & \bm r \nonumber \\ 30 | & p_{dc} \nonumber \\ 31 | % 32 | \mbox{\bf variables: } & \nonumber \\ 33 | & S^g_k \;\; \forall k\in G \nonumber \\ 34 | & S^c_l \;\; \forall l\in C \nonumber \\ 35 | & P^{c, dc}_l \;\; \forall l\in C \nonumber \\ 36 | & V_i \;\; \forall i\in N \nonumber \\ 37 | & V_i \;\; \forall i\in N_{dc} \nonumber \\ 38 | & S_{ij} \;\; \forall (i,j) \in E \cup E^R \nonumber \\ 39 | % 40 | \end{align} \\ 41 | \begin{align} 42 | \mbox{\bf minimize: } & \sum_{k \in G} \bm c_{2k} (\Re(S^g_k))^2 + \bm c_{1k}\Re(S^g_k) + \bm c_{0k} \\ 43 | % 44 | \mbox{\bf subject to: } & \nonumber \\ 45 | & \angle V_{\bm r} = 0 \\ 46 | & \bm {S^{gl}}_k \leq S^g_k \leq \bm {S^{gu}}_k \;\; \forall k \in G \\ 47 | & \bm {v^l}_i \leq |V_i| \leq \bm {v^u}_i \;\; \forall i \in N \\ 48 | & \sum_{\substack{k \in G_i}} S^g_k + \sum_{\substack{l \in C_i}} S^c_l - {\bm S^d_i} - \bm Y^s_{i} |V_i|^2 = \sum_{\substack{(i,j)\in E_i \cup E_i^R}} S_{ij} \;\; \forall i\in N \\ 49 | & S_{ij} = \left( \bm Y^*_{ij} - \bm i\frac{\bm {b^c}_{ij}}{2} \right) \frac{|V_i|^2}{|\bm{T}_{ij}|^2} - \bm Y^*_{ij} \frac{V_i V^*_j}{\bm{T}_{ij}} \;\; \forall (i,j)\in E \\ 50 | & S_{ji} = \left( \bm Y^*_{ij} - \bm i\frac{\bm {b^c}_{ij}}{2} \right) |V_j|^2 - \bm Y^*_{ij} \frac{V^*_i V_j}{\bm{T}^*_{ij}} \;\; \forall (i,j)\in E \\ 51 | & |S_{ij}| \leq \bm {s^u}_{ij} \;\; \forall (i,j) \in E \cup E^R \\ 52 | & \bm {\theta^{\Delta l}}_{ij} \leq \angle (V_i V^*_j) \leq \bm {\theta^{\Delta u}}_{ij} \;\; \forall (i,j) \in E \\ 53 | & \bm {S^{cl}}_l \leq S^c_l \leq \bm {S^{cu}}_l \;\; \forall l \in C \\ 54 | & \sum_{\substack{k \in G_i}} P^g_k + \sum_{\substack{l \in C_i}} P^{c, dc}_l - {\bm P^d_i} - \bm Y^s_{i} |V_i|^2 = \sum_{\substack{(i,j)\in E_{i, dc} \cup E_{i, dc}^R}} P_{ij} \;\; \forall i\in N_{dc} \\ 55 | & P_{ij} = p_{dc} \bm Y_{ij} \cdot( V_i^2 - V_i V_j) \;\; \forall (i,j)\in E_{dc} \\ 56 | & |P_{ij}| \leq \bm {p^u}_{ij} \;\; \forall (i,j) \in E_{dc} \cup E_{dc}^R \\ 57 | & P^c_l + P^{c, dc}_l = a + b |I^c_l| + c |I^c_l|^2 \;\; \forall l \in C \\ 58 | & |V_i|^2 |I^c_l|^2 = (S^c_l)^2 \;\; \forall l \in C_i \;\; \forall i \in N \\ 59 | % 60 | \end{align} 61 | 62 | \end{document} -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # Power Grid Lib - Optimal Power Flow with HVDC Lines 2 | 3 | This benchmark library is curated and maintained by the [IEEE PES Task Force on Benchmarks for Validation of Emerging Power System Algorithms](https://power-grid-lib.github.io/) and is designed to evaluate a variation of the Optimal Power Flow problem with HVDC lines. Specifically, these cases are designed for benchmarking algorithms that solve the following the Non-Convex Nonlinear Program described in the formulation [PDF document](MODEL.pdf). An extensive description of the case data and optimization model is provided in the [transactions paper](https://ieeexplore.ieee.org/document/8636236). 4 | 5 | All of the cases files are curated in an extended version of the [MatACDC](https://www.esat.kuleuven.be/electa/teaching/matacdc/MatACDCManual) data format. An open-source reference implementations are available in [PowerModelsACDC.jl](https://github.com/hakanergun/PowerModelsACDC.jl). 6 | 7 | ## Problem Overview 8 | 9 | These cases are useful for benchmarking solution methods for a variant of the optimal power flow problem common in the academic literature. The features of this model are: 10 | * Respresentation of the HVDC converter station including converter transformer, harmonic filter and quadratic losses 11 | * Point to point and fully meshed HVDC configuations 12 | * Optimization of the ac side active and reactive power set points and dc side active power set points for all converters 13 | * Parametrized model to omit converter transformers or filters 14 | 15 | 16 | ## Case File Overview 17 | 18 | * CASE_5_3 Power flow data for modified 5 bus, 5 gen, 3 bus dc case based on PJM 5-bus system 19 | * CASE_24_7 based on the IEEE reliability test system with HVDC grid connecting three zones 20 | * CASE_39_10 Power flow data for 39 bus New England system with additional 10 converter stations and 12 dc branches 21 | * CASE_3120_5 Power flow data for Polish system - summer 2008 morning peak, with 5 additional converter stations and 5 dc branches 22 | 23 | ## Example script to run test case 24 | ``` 25 | using PowerModelsACDC, PowerModels, Ipopt 26 | 27 | 28 | ipopt = with_optimizer(Ipopt.Optimizer, tol=1e-6, print_level=0) 29 | s = Dict("output" => Dict("branch_flows" => true), "conv_losses_mp" => true) 30 | 31 | file5_3 = "case5_3_he.m" 32 | data = PowerModels.parse_file(file5_3) 33 | data["convdc"] = data["dcconv"] 34 | data["busdc"] = data["dcbus"] 35 | data["branchdc"] = data["dcbranch"] 36 | PowerModelsACDC.process_additional_data!(data) 37 | result5_3 = run_acdcopf(data, ACPPowerModel, ipopt; setting = s) 38 | ``` 39 | ## Contributions 40 | 41 | All case files are provided under a [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which allows anyone to share or adapt these cases as long as they give appropriate credit to the original author, provide a link to the license, and indicate if changes were made. 42 | 43 | Community-based recommendations and contributions are welcome and encouraged in all PGLib repositories. Please feel free to submit comments and questions in the [issue tracker](https://github.com/power-grid-lib/pglib-uc/issues). Corrections and new network contributions are welcome via pull requests. All data contributions are subject to a quality assurance review by the repository curator(s). 44 | 45 | 46 | ## Citation Guidelines 47 | 48 | This repository is not static. Consequently, it is critically important to indicate the version number when referencing this repository in scholarly work. 49 | 50 | Users of this these cases are encouraged to cite the original source documents mentioned in this overview document. 51 | 52 | 53 | 54 | -------------------------------------------------------------------------------- /case24_7_jb.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | %%%% %%%%% 3 | %%%% IEEE PES Power Grid Library - Optimal Power Flow - v23.09 %%%%% 4 | %%%% (https://github.com/power-grid-lib/pglib-opf) %%%%% 5 | %%%% Benchmark Group - Typical Operations %%%%% 6 | %%%% 10 - May - 2019 %%%%% 7 | %%%% %%%%% 8 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 | 10 | % CASE24_IEEE_RTS1996_3zones Power flow data for system based on the 11 | % IEEE RELIABILITY TEST SYSTEM. 12 | % Please see CASEFORMAT for details on the case file format. 13 | % 14 | % This system data is based on the MATPOWER case file CASE24_IEEE_RTS 15 | % which is based on the IEEE RELIABILITY TEST SYSTEM 16 | % 17 | % The data has been adopted to corresponding with the 18 | % IEEE Two Area RTS-96 data from... 19 | % IEEE Reliability Test System Task Force of Applications of 20 | % Probability Methods Subcommittee, "IEEE reliability test system-96," 21 | % IEEE Transactions on Power Systems, Vol. 14, No. 3, Aug. 1999, 22 | % pp. 1010-1020. 23 | % 24 | % The IEEE Two Area RTS-96 network has been extended and now includes 3 25 | % asynchronous zones (node numbers 1xx, 2xx and 3xx). 26 | % Data on zone 1 and 2 taken from the IEEE Two Area RTS-96 with following 27 | % adaptations: 28 | % - nodes renumbered according to IEEE Two Area RTS-96 data 29 | % - gen U100 at node 107 disabled (commented) 30 | % - gen U76 at node 201 disabled (commented) 31 | % - slack node zone 2: node 213 32 | % - lines 107-203, 113-215, 123-217 removed 33 | % Data on zone 3 added: 34 | % - nodes 301 and 302 35 | % - gen at node 302 36 | % - line 301-302 37 | % 38 | % MATPOWER case file data provided by Bruce Wollenberg 39 | % (MATPOWER file case24_ieee_rts.m) and adapted for use with MatACDC 40 | % by Jef Beerten. 41 | 42 | 43 | function mpc = case24_3zones_acdc() 44 | %% MATPOWER Case Format : Version 1 45 | mpc.version = '1'; 46 | 47 | %%----- Power Flow Data -----%% 48 | %% system MVA base 49 | mpc.baseMVA = 100; 50 | 51 | %% bus data 52 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 53 | mpc.bus = [ 54 | 101 2 108 22 0 0 11 1.035 0 138 1 1.05 0.95; 55 | 102 2 97 20 0 0 11 1.035 0 138 1 1.05 0.95; 56 | 103 1 180 37 0 0 11 1 0 138 1 1.05 0.95; 57 | 104 1 74 15 0 0 11 1 0 138 1 1.05 0.95; 58 | 105 1 71 14 0 0 11 1 0 138 1 1.05 0.95; 59 | 106 1 136 28 0 1.00 12 1 0 138 1 1.05 0.95; 60 | 107 2 125 25 0 0 12 1.025 0 138 1 1.05 0.95; 61 | 108 1 171 35 0 0 12 1 0 138 1 1.05 0.95; 62 | 109 1 175 36 0 0 11 1 0 138 1 1.05 0.95; 63 | 110 1 195 40 0 0 12 1 0 138 1 1.05 0.95; 64 | 111 1 0 0 0 0 13 1 0 230 1 1.05 0.95; 65 | 112 1 0 0 0 0 13 1 0 230 1 1.05 0.95; 66 | 113 3 265 54 0 0 13 1.02 0 230 1 1.05 0.95; 67 | 114 2 194 39 0 0 13 0.98 0 230 1 1.05 0.95; 68 | 115 2 317 64 0 0 14 1.014 0 230 1 1.05 0.95; 69 | 116 2 100 20 0 0 14 1.017 0 230 1 1.05 0.95; 70 | 117 1 0 0 0 0 14 1 0 230 1 1.05 0.95; 71 | 118 2 333 68 0 0 14 1.05 0 230 1 1.05 0.95; 72 | 119 1 181 37 0 0 13 1 0 230 1 1.05 0.95; 73 | 120 1 128 26 0 0 13 1 0 230 1 1.05 0.95; 74 | 121 2 0 0 0 0 14 1.05 0 230 1 1.05 0.95; 75 | 122 2 0 0 0 0 14 1.05 0 230 1 1.05 0.95; 76 | 123 1 0 0 0 0 13 1.05 0 230 1 1.05 0.95; 77 | 124 1 0 0 0 0 14 1 0 230 1 1.05 0.95; 78 | 201 2 108 22 0 0 11 1.035 0 138 2 1.05 0.95; 79 | 202 2 97 20 0 0 11 1.035 0 138 2 1.05 0.95; 80 | 203 1 180 37 0 0 11 1 0 138 2 1.05 0.95; 81 | 204 1 74 15 0 0 11 1 0 138 2 1.05 0.95; 82 | 205 1 71 14 0 0 11 1 0 138 2 1.05 0.95; 83 | 206 1 136 28 0 1.00 12 1 0 138 2 1.05 0.95; 84 | 207 2 125 25 0 0 12 1.025 0 138 2 1.05 0.95; 85 | 208 1 171 35 0 0 12 1 0 138 2 1.05 0.95; 86 | 209 1 175 36 0 0 11 1 0 138 2 1.05 0.95; 87 | 210 1 195 40 0 0 12 1 0 138 2 1.05 0.95; 88 | 211 1 0 0 0 0 13 1 0 230 2 1.05 0.95; 89 | 212 1 0 0 0 0 13 1 0 230 2 1.05 0.95; 90 | 213 3 265 54 0 0 13 1.02 0 230 2 1.05 0.95; 91 | 214 2 194 39 0 0 13 0.98 0 230 2 1.05 0.95; 92 | 215 2 317 64 0 0 14 1.014 0 230 2 1.05 0.95; 93 | 216 2 100 20 0 0 14 1.017 0 230 2 1.05 0.95; 94 | 217 1 0 0 0 0 14 1 0 230 2 1.05 0.95; 95 | 218 2 333 68 0 0 14 1.05 0 230 2 1.05 0.95; 96 | 219 1 181 37 0 0 13 1 0 230 2 1.05 0.95; 97 | 220 1 128 26 0 0 13 1 0 230 2 1.05 0.95; 98 | 221 2 0 0 0 0 14 1.05 0 230 2 1.05 0.95; 99 | 222 2 0 0 0 0 14 1.05 0 230 2 1.05 0.95; 100 | 223 2 0 0 0 0 13 1.05 0 230 2 1.05 0.95; 101 | 224 1 0 0 0 0 14 1 0 230 2 1.05 0.95; 102 | 301 1 0 0 0 0 14 1 0 230 3 1.05 0.95; 103 | 302 3 0 0 0 0 14 1.025 0 230 3 1.05 0.95; 104 | ]; 105 | 106 | %% generator data 107 | % 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 % Unit Code 108 | mpc.gen = [ 109 | 101 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 110 | 101 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 111 | 101 76 14.1 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 112 | 101 76 14.1 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 113 | 102 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 114 | 102 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 115 | 102 76 7.0 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 116 | 102 76 7.0 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 117 | 107 80 17.2 60 0 1.025 100 1 100 25 0 0 0 0 0 0 0 0 0 0 0; % U100 118 | 107 80 17.2 60 0 1.025 100 1 100 25 0 0 0 0 0 0 0 0 0 0 0; % U100 119 | 113 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 120 | 113 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 121 | 113 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 122 | 114 0 13.7 200 -50 0.980 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; % SynCond 123 | 115 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 124 | 115 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 125 | 115 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 126 | 115 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 127 | 115 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 128 | 115 155 0.05 80 -50 1.014 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 129 | 116 155 25.22 80 -50 1.017 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 130 | 118 400 137.4 200 -50 1.05 100 1 400 100 0 0 0 0 0 0 0 0 0 0 0; % U400 131 | 121 400 108.2 200 -50 1.05 100 1 400 100 0 0 0 0 0 0 0 0 0 0 0; % U400 132 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 133 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 134 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 135 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 136 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 137 | 122 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 138 | 123 155 31.79 80 -50 1.05 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 139 | 123 155 31.79 80 -50 1.05 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 140 | 123 350 71.78 150 -25 1.05 100 1 350 140 0 0 0 0 0 0 0 0 0 0 0; % U350 141 | 201 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 142 | 201 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 143 | 201 76 14.1 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 144 | 202 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 145 | 202 10 0 10 0 1.035 100 1 20 16 0 0 0 0 0 0 0 0 0 0 0; % U20 146 | 202 76 7.0 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 147 | 202 76 7.0 30 -25 1.035 100 1 76 15.2 0 0 0 0 0 0 0 0 0 0 0; % U76 148 | 207 80 17.2 60 0 1.025 100 1 100 25 0 0 0 0 0 0 0 0 0 0 0; % U100 149 | 207 80 17.2 60 0 1.025 100 1 100 25 0 0 0 0 0 0 0 0 0 0 0; % U100 150 | 207 80 17.2 60 0 1.025 100 1 100 25 0 0 0 0 0 0 0 0 0 0 0; % U100 151 | 213 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 152 | 213 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 153 | 213 95.1 40.7 80 0 1.020 100 1 197 69 0 0 0 0 0 0 0 0 0 0 0; % U197 154 | 214 0 13.68 200 -50 0.980 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; % SynCond 155 | 215 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 156 | 215 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 157 | 215 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 158 | 215 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 159 | 215 12 0 6 0 1.014 100 1 12 2.4 0 0 0 0 0 0 0 0 0 0 0; % U12 160 | 215 155 0.048 80 -50 1.014 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 161 | 216 155 25.22 80 -50 1.017 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 162 | 218 400 137.4 200 -50 1.05 100 1 400 100 0 0 0 0 0 0 0 0 0 0 0; % U400 163 | 221 400 108.2 200 -50 1.05 100 1 400 100 0 0 0 0 0 0 0 0 0 0 0; % U400 164 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 165 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 166 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 167 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 168 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 169 | 222 50 -4.96 16 -10 1.05 100 1 50 10 0 0 0 0 0 0 0 0 0 0 0; % U50 170 | 223 155 31.79 80 -50 1.05 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 171 | 223 155 31.79 80 -50 1.05 100 1 155 54.3 0 0 0 0 0 0 0 0 0 0 0; % U155 172 | 223 350 71.78 150 -25 1.05 100 1 350 140 0 0 0 0 0 0 0 0 0 0 0; % U350 173 | 302 150 10 150 -25 1.05 100 1 350 140 0 0 0 0 0 0 0 0 0 0 0; % U350 174 | ]; 175 | 176 | %% branch data 177 | % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 178 | mpc.branch = [ 179 | 101 102 0.003 0.014 0.461 175 250 200 0 0 1 -360 360; 180 | 101 103 0.055 0.211 0.057 175 208 220 0 0 1 -360 360; 181 | 101 105 0.022 0.085 0.023 175 208 220 0 0 1 -360 360; 182 | 102 104 0.033 0.127 0.034 175 208 220 0 0 1 -360 360; 183 | 102 106 0.050 0.192 0.052 175 208 220 0 0 1 -360 360; 184 | 103 109 0.031 0.119 0.032 175 208 220 0 0 1 -360 360; 185 | 103 124 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 186 | 104 109 0.027 0.104 0.028 175 208 220 0 0 1 -360 360; 187 | 105 110 0.022 0.088 0.024 175 208 220 0 0 1 -360 360; 188 | 106 110 0.014 0.061 2.459 400 193 200 0 0 1 -360 360; % limit changed for feasibility of OPF 189 | 107 108 0.016 0.061 0.017 175 208 220 0 0 1 -360 360; 190 | 108 109 0.043 0.165 0.045 175 208 220 0 0 1 -360 360; 191 | 108 110 0.043 0.165 0.045 175 208 220 0 0 1 -360 360; 192 | 109 111 0.002 0.084 0 400 510 600 1.03 0 1 -360 360; 193 | 109 112 0.002 0.084 0 400 510 600 1.03 0 1 -360 360; 194 | 110 111 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 195 | 110 112 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 196 | 111 113 0.006 0.048 0.100 500 600 625 0 0 1 -360 360; 197 | 111 114 0.005 0.042 0.088 500 625 625 0 0 1 -360 360; 198 | 112 113 0.006 0.048 0.100 500 625 625 0 0 1 -360 360; 199 | 112 123 0.012 0.097 0.203 500 625 625 0 0 1 -360 360; 200 | 113 123 0.011 0.087 0.182 500 625 625 0 0 1 -360 360; 201 | 114 116 0.005 0.059 0.082 500 625 625 0 0 1 -360 360; % x value changed 202 | 115 116 0.002 0.017 0.036 500 600 625 0 0 1 -360 360; 203 | 115 121 0.006 0.049 0.103 500 600 625 0 0 1 -360 360; 204 | 115 121 0.006 0.049 0.103 500 600 625 0 0 1 -360 360; 205 | 115 124 0.007 0.052 0.109 500 600 625 0 0 1 -360 360; 206 | 116 117 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 207 | 116 119 0.003 0.023 0.049 500 600 625 0 0 1 -360 360; 208 | 117 118 0.002 0.014 0.030 500 600 625 0 0 1 -360 360; 209 | 117 122 0.014 0.105 0.221 500 600 625 0 0 1 -360 360; 210 | 118 121 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 211 | 118 121 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 212 | 119 120 0.005 0.040 0.083 500 600 625 0 0 1 -360 360; 213 | 119 120 0.005 0.040 0.083 500 600 625 0 0 1 -360 360; 214 | 120 123 0.003 0.022 0.046 500 600 625 0 0 1 -360 360; 215 | 120 123 0.003 0.022 0.046 500 600 625 0 0 1 -360 360; 216 | 121 122 0.009 0.068 0.142 500 600 625 0 0 1 -360 360; 217 | 201 202 0.003 0.014 0.461 175 250 200 0 0 1 -360 360; 218 | 201 203 0.055 0.211 0.057 175 208 220 0 0 1 -360 360; 219 | 201 205 0.022 0.085 0.023 175 208 220 0 0 1 -360 360; 220 | 202 204 0.033 0.127 0.034 175 208 220 0 0 1 -360 360; 221 | 202 206 0.050 0.192 0.052 175 208 220 0 0 1 -360 360; 222 | 203 209 0.031 0.119 0.032 175 208 220 0 0 1 -360 360; 223 | 203 224 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 224 | 204 209 0.027 0.104 0.028 175 208 220 0 0 1 -360 360; 225 | 205 210 0.022 0.088 0.024 175 208 220 0 0 1 -360 360; 226 | 206 210 0.014 0.061 2.459 400 193 200 0 0 1 -360 360; % line limit changed for feasibility 227 | 207 208 0.016 0.061 0.017 175 208 220 0 0 1 -360 360; 228 | 208 209 0.043 0.165 0.045 175 208 220 0 0 1 -360 360; 229 | 208 210 0.043 0.165 0.045 175 208 220 0 0 1 -360 360; 230 | 209 211 0.002 0.084 0 400 510 600 1.03 0 1 -360 360; 231 | 209 212 0.002 0.084 0 400 510 600 1.03 0 1 -360 360; 232 | 210 211 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 233 | 210 212 0.002 0.084 0 400 510 600 1.015 0 1 -360 360; 234 | 211 213 0.006 0.048 0.100 500 600 625 0 0 1 -360 360; 235 | 211 214 0.005 0.042 0.088 500 625 625 0 0 1 -360 360; 236 | 212 213 0.006 0.048 0.100 500 625 625 0 0 1 -360 360; 237 | 212 223 0.012 0.097 0.203 500 625 625 0 0 1 -360 360; 238 | 213 223 0.011 0.087 0.182 500 625 625 0 0 1 -360 360; 239 | 214 216 0.005 0.059 0.082 500 625 625 0 0 1 -360 360; % x value changed 240 | 215 216 0.002 0.017 0.036 500 600 625 0 0 1 -360 360; 241 | 215 221 0.006 0.049 0.103 500 600 625 0 0 1 -360 360; 242 | 215 221 0.006 0.049 0.103 500 600 625 0 0 1 -360 360; 243 | 215 224 0.007 0.052 0.109 500 600 625 0 0 1 -360 360; 244 | 216 217 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 245 | 216 219 0.003 0.023 0.049 500 600 625 0 0 1 -360 360; 246 | 217 218 0.002 0.014 0.030 500 600 625 0 0 1 -360 360; 247 | 217 222 0.014 0.105 0.221 500 600 625 0 0 1 -360 360; 248 | 218 221 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 249 | 218 221 0.003 0.026 0.055 500 600 625 0 0 1 -360 360; 250 | 219 220 0.005 0.040 0.083 500 600 625 0 0 1 -360 360; 251 | 219 220 0.005 0.040 0.083 500 600 625 0 0 1 -360 360; 252 | 220 223 0.003 0.022 0.046 500 600 625 0 0 1 -360 360; 253 | 220 223 0.003 0.022 0.046 500 600 625 0 0 1 -360 360; 254 | 221 222 0.009 0.068 0.142 500 600 625 0 0 1 -360 360; 255 | 301 302 0.000 0.001 0.000 500 600 625 0 0 1 -360 360; 256 | ]; 257 | %% dc grid topology 258 | %colunm_names% dcpoles 259 | mpc.dcpol=2; % numbers of poles (1=monopolar grid, 2=bipolar grid) 260 | 261 | %% bus data 262 | %column_names% busdc_i grid Pdc Vdc basekVdc Vdcmax Vdcmin Cdc 263 | mpc.dcbus = [ 264 | 1 1 0 1 150 1.1 0.9 0; 265 | 2 1 0 1 150 1.1 0.9 0; 266 | 3 1 0 1 150 1.1 0.9 0; 267 | 4 2 0 1 300 1.1 0.9 0; 268 | 5 2 0 1 300 1.1 0.9 0; 269 | 6 2 0 1 300 1.1 0.9 0; 270 | 7 2 0 1 300 1.1 0.9 0; 271 | ]; 272 | 273 | %% converters 274 | %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 275 | mpc.dcconv = [ 276 | 1 107 2 1 0 50 0 1 0.001 0.10 1 1 0.09 0 0.0001 0.16 0 138 1.2 0.9 1.1 1 1.103 0.887 2.885 4.371 0 0 1 0 200 -200 200 -200; 277 | 2 204 1 2 75.3 -50 0 1 0.001 0.10 1 1 0.09 0 0.0001 0.16 0 138 1.2 0.9 1.1 1 1.103 0.887 2.885 4.371 0 0 1 0 200 -200 200 -200; 278 | 3 301 1 1 -141.9 130 0 1 0.001 0.05 1 1 0.045 0 0.0001 0.08 0 138 1.2 0.9 2.2 1 2.206 0.887 1.442 2.1850 0 0 1 0 200 -200 200 -200; 279 | 4 113 2 1 131.5 75.9 0 1 0.0005 0.05 1 1 0 0 0.0001 0.08 0 345 1.2 0.5 2.2 1 2.206 1.8 5.94 9 0 0 1 0 200 -200 200 -200; 280 | 5 123 1 1 -61.7 0 0 1 0.001 0.10 1 1 0 0 0.0001 0.16 0 345 1.2 0.5 1.1 1 1.103 1.8 11.88 18 0 0 1 0 200 -200 200 -200; 281 | 6 215 1 2 -123.4 -10 0 1 0.0005 0.05 1 1 0 0 0.0001 0.08 0 345 1.2 0.5 2.2 1 2.206 1.8 5.94 9 0 0 1 0 200 -200 200 -200; 282 | 7 217 1 1 50 20 0 1 0.001 0.10 1 1 0 0 0.0001 0.16 0 345 1.2 0.5 1.1 1 1.103 1.8 11.88 18 0 0 1 0 200 -200 200 -200; 283 | ]; 284 | 285 | %% branches 286 | %column_names% fbusdc tbusdc r l c rateA rateB rateC status 287 | mpc.dcbranch = [ 288 | 1 3 0.0352 0 0 100 100 100 1; 289 | 2 3 0.0352 0 0 100 100 100 1; 290 | 4 5 0.0828 0 0 100 100 100 1; 291 | 4 7 0.0704 0 0 100 100 100 1; 292 | 4 6 0.0718 0 0 100 100 100 1; 293 | 5 7 0.0760 0 0 100 100 100 1; 294 | 6 7 0.0248 0 0 100 100 100 1; 295 | ]; 296 | 297 | %% generator cost data 298 | % 2 startup shutdown n c(n-1) ... c0 299 | mpc.gencost = [ 300 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 301 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 302 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 303 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 304 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 305 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 306 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 307 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 308 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 309 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 310 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 311 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 312 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 313 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 314 | 2 1500.0 0.0 3 0.000000 0.000000 0.000000; 315 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 316 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 317 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 318 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 319 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 320 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 321 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 322 | 2 1500.0 0.0 3 0.000213 4.423100 395.374900; 323 | 2 1500.0 0.0 3 0.000213 4.423100 395.374900; 324 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 325 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 326 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 327 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 328 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 329 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 330 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 331 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 332 | 2 1500.0 0.0 3 0.004895 11.849500 665.109400; 333 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 334 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 335 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 336 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 337 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 338 | 2 1500.0 0.0 3 0.000000 130.000000 400.684900; 339 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 340 | 2 1500.0 0.0 3 0.014142 16.081100 212.307600; 341 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 342 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 343 | 2 1500.0 0.0 3 0.052672 43.661500 781.521000; 344 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 345 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 346 | 2 1500.0 0.0 3 0.007170 48.580400 832.757500; 347 | 2 1500.0 0.0 3 0.000000 0.000000 0.000000; 348 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 349 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 350 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 351 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 352 | 2 1500.0 0.0 3 0.328412 56.564000 86.385200; 353 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 354 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 355 | 2 1500.0 0.0 3 0.000213 4.423100 395.374900; 356 | 2 1500.0 0.0 3 0.000213 4.423100 395.374900; 357 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 358 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 359 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 360 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 361 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 362 | 2 1500.0 0.0 3 0.000000 0.001000 0.001000; 363 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 364 | 2 1500.0 0.0 3 0.008342 12.388300 382.239100; 365 | ]; 366 | -------------------------------------------------------------------------------- /case39_10_he.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | %%%% %%%%% 3 | %%%% IEEE PES Power Grid Library - Optimal Power Flow - v23.09 %%%%% 4 | %%%% (https://github.com/power-grid-lib/pglib-opf) %%%%% 5 | %%%% Benchmark Group - Typical Operations %%%%% 6 | %%%% 10 - May - 2019 %%%%% 7 | %%%% %%%%% 8 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 | % 10 | % CASE39 Power flow data for 39 bus New England system. 11 | % 12 | % Data taken from [1] with the following modifications/additions: 13 | % 14 | % - renumbered gen buses consecutively (as in [2] and [4]) 15 | % - added Pmin = 0 for all gens 16 | % - added Qmin, Qmax for gens at 31 & 39 (copied from gen at 35) 17 | % - added Vg based on V in bus data (missing for bus 39) 18 | % - added Vg, Pg, Pd, Qd at bus 39 from [2] (same in [4]) 19 | % - added Pmax at bus 39: Pmax = Pg + 100 20 | % - added line flow limits and area data from [4] 21 | % - added voltage limits, Vmax = 1.06, Vmin = 0.94 22 | % - added identical quadratic generator costs 23 | % - increased Pmax for gen at bus 34 from 308 to 508 24 | % (assumed typo in [1], makes initial solved case feasible) 25 | % - re-solved power flow 26 | % 27 | % Notes: 28 | % - Bus 39, its generator and 2 connecting lines were added 29 | % (by authors of [1]) to represent the interconnection with 30 | % the rest of the eastern interconnect, and did not include 31 | % Vg, Pg, Qg, Pd, Qd, Pmin, Pmax, Qmin or Qmax. 32 | % - As the swing bus, bus 31 did not include and Q limits. 33 | % - The voltages, etc in [1] appear to be quite close to the 34 | % power flow solution of the case before adding bus 39 with 35 | % it's generator and connecting branches, though the solution 36 | % is not exact. 37 | % - Explicit voltage setpoints for gen buses are not given, so 38 | % they are taken from the bus data, however this results in two 39 | % binding Q limits at buses 34 & 37, so the corresponding 40 | % voltages have probably deviated from their original setpoints. 41 | % - The generator locations and types are as follows: 42 | % 1 30 hydro 43 | % 2 31 nuke01 44 | % 3 32 nuke02 45 | % 4 33 fossil02 46 | % 5 34 fossil01 47 | % 6 35 nuke03 48 | % 7 36 fossil04 49 | % 8 37 nuke04 50 | % 9 38 nuke05 51 | % 10 39 interconnection to rest of US/Canada 52 | % 53 | % This is a solved power flow case, but it includes the following 54 | % violations: 55 | % - Pmax violated at bus 31: Pg = 677.87, Pmax = 646 56 | % - Qmin violated at bus 37: Qg = -1.37, Qmin = 0 57 | % 58 | % References: 59 | % [1] G. W. Bills, et.al., "On-Line Stability Analysis Study" 60 | % RP90-1 Report for the Edison Electric Institute, October 12, 1970, 61 | % pp. 1-20 - 1-35. 62 | % prepared by E. M. Gulachenski - New England Electric System 63 | % J. M. Undrill - General Electric Co. 64 | % "generally representative of the New England 345 KV system, but is 65 | % not an exact or complete model of any past, present or projected 66 | % configuration of the actual New England 345 KV system. 67 | % [2] M. A. Pai, Energy Function Analysis for Power System Stability, 68 | % Kluwer Academic Publishers, Boston, 1989. 69 | % (references [3] as source of data) 70 | % [3] Athay, T.; Podmore, R.; Virmani, S., "A Practical Method for the 71 | % Direct Analysis of Transient Stability," IEEE Transactions on Power 72 | % Apparatus and Systems , vol.PAS-98, no.2, pp.573-584, March 1979. 73 | % URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4113518&isnumber=4113486 74 | % (references [1] as source of data) 75 | % [4] Data included with TC Calculator at http://www.pserc.cornell.edu/tcc/ 76 | % for 39-bus system. 77 | % [5] DC grid data added from: Rimez, J., 2014. Optimal Operation of 78 | % Hybrid AC/DC Meshed Grids (Optimale uitbating van hybriede vermaasde 79 | % AC/DC netwerken), PhD Thesis, KU Leuven. 80 | % 81 | % 82 | % 83 | % Created by Hakan Ergun in 2019. 84 | % 85 | % Copyright (c) 1989 by The Institute of Electrical and Electronics Engineers (IEEE) 86 | % Licensed under the Creative Commons Attribution 4.0 87 | % International license, http://creativecommons.org/licenses/by/4.0/ 88 | % 89 | % Contact M.E. Brennan (me.brennan@ieee.org) for inquries on further reuse of 90 | % this dataset. 91 | % 92 | function mpc = case39_10_he 93 | mpc.version = '2'; 94 | mpc.baseMVA = 100.0; 95 | 96 | %% bus data 97 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 98 | mpc.bus = [ 99 | 1 1 97.6 44.2 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 100 | 2 1 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 101 | 3 1 322.0 2.4 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 102 | 4 1 500.0 184.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 103 | 5 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 104 | 6 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 105 | 7 1 233.8 84.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 106 | 8 1 522.0 176.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 107 | 9 1 6.5 -66.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 108 | 10 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 109 | 11 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 110 | 12 1 8.53 88.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 111 | 13 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 112 | 14 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 113 | 15 1 320.0 153.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 114 | 16 1 329.0 32.3 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 115 | 17 1 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 116 | 18 1 158.0 30.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 117 | 19 1 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 118 | 20 1 680.0 103.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 119 | 21 1 274.0 115.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 120 | 22 1 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 121 | 23 1 247.5 84.6 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 122 | 24 1 308.6 -92.2 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 123 | 25 1 224.0 47.2 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 124 | 26 1 139.0 17.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 125 | 27 1 281.0 75.5 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 126 | 28 1 206.0 27.6 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 127 | 29 1 283.5 26.9 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 128 | 30 2 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 129 | 31 3 9.2 4.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 130 | 32 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 131 | 33 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 132 | 34 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 133 | 35 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 134 | 36 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 135 | 37 2 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000; 136 | 38 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000; 137 | 39 2 1104.0 250.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000; 138 | ]; 139 | 140 | %% generator data 141 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 142 | mpc.gen = [ 143 | 30 520.0 270.0 400.0 140.0 1.0 100.0 1 1040.0 0.0; % NUC 144 | 31 323.0 100.0 300.0 -100.0 1.0 100.0 1 646.0 0.0; % COW 145 | 32 362.5 225.0 300.0 150.0 1.0 100.0 1 725.0 0.0; % COW 146 | 33 326.0 125.0 250.0 0.0 1.0 100.0 1 652.0 0.0; % COW 147 | 34 254.0 83.5 167.0 0.0 1.0 100.0 1 508.0 0.0; % COW 148 | 35 343.5 100.0 300.0 -100.0 1.0 100.0 1 687.0 0.0; % COW 149 | 36 290.0 120.0 240.0 0.0 1.0 100.0 1 580.0 0.0; % COW 150 | 37 282.0 125.0 250.0 0.0 1.0 100.0 1 564.0 0.0; % COW 151 | 38 432.5 75.0 300.0 -150.0 1.0 100.0 1 865.0 0.0; % COW 152 | 39 550.0 100.0 300.0 -100.0 1.0 100.0 1 1100.0 0.0; % COW 153 | ]; 154 | 155 | %% generator cost data 156 | % 2 startup shutdown n c(n-1) ... c0 157 | mpc.gencost = [ 158 | 2 0.0 0.0 3 0.000000 6.724778 0.000000; % NUC 159 | 2 0.0 0.0 3 0.000000 14.707625 0.000000; % COW 160 | 2 0.0 0.0 3 0.000000 24.804734 0.000000; % COW 161 | 2 0.0 0.0 3 0.000000 34.844643 0.000000; % COW 162 | 2 0.0 0.0 3 0.000000 24.652994 0.000000; % COW 163 | 2 0.0 0.0 3 0.000000 32.306483 0.000000; % COW 164 | 2 0.0 0.0 3 0.000000 18.157477 0.000000; % COW 165 | 2 0.0 0.0 3 0.000000 31.550181 0.000000; % COW 166 | 2 0.0 0.0 3 0.000000 22.503168 0.000000; % COW 167 | 2 0.0 0.0 3 0.000000 27.434444 0.000000; % COW 168 | ]; 169 | 170 | %% branch data 171 | % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 172 | mpc.branch = [ 173 | 1 2 0.0035 0.0411 0.6987 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 174 | 1 39 0.001 0.025 0.75 1000.0 1000.0 1000.0 0.0 0.0 1 -30.0 30.0; 175 | 2 3 0.0013 0.0151 0.2572 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0; 176 | 2 25 0.007 0.0086 0.146 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0; 177 | 2 30 0.0 0.0181 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0; 178 | 3 4 0.0013 0.0213 0.2214 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0; 179 | 3 18 0.0011 0.0133 0.2138 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0; 180 | 4 5 0.0008 0.0128 0.1342 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 181 | 4 14 0.0008 0.0129 0.1382 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0; 182 | 5 6 0.0002 0.0026 0.0434 1200.0 1200.0 1200.0 0.0 0.0 1 -30.0 30.0; 183 | 5 8 0.0008 0.0112 0.1476 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 184 | 6 7 0.0006 0.0092 0.113 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 185 | 6 11 0.0007 0.0082 0.1389 480.0 480.0 480.0 0.0 0.0 1 -30.0 30.0; 186 | 6 31 0.0 0.025 0.0 1800.0 1800.0 1800.0 1.07 0.0 1 -30.0 30.0; 187 | 7 8 0.0004 0.0046 0.078 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 188 | 8 9 0.0023 0.0363 0.3804 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 189 | 9 39 0.001 0.025 1.2 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 190 | 10 11 0.0004 0.0043 0.0729 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 191 | 10 13 0.0004 0.0043 0.0729 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 192 | 10 32 0.0 0.02 0.0 900.0 900.0 2500.0 1.07 0.0 1 -30.0 30.0; 193 | 12 11 0.0016 0.0435 0.0 500.0 500.0 500.0 1.006 0.0 1 -30.0 30.0; 194 | 12 13 0.0016 0.0435 0.0 500.0 500.0 500.0 1.006 0.0 1 -30.0 30.0; 195 | 13 14 0.0009 0.0101 0.1723 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 196 | 14 15 0.0018 0.0217 0.366 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 197 | 15 16 0.0009 0.0094 0.171 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 198 | 16 17 0.0007 0.0089 0.1342 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 199 | 16 19 0.0016 0.0195 0.304 600.0 600.0 2500.0 0.0 0.0 1 -30.0 30.0; 200 | 16 21 0.0008 0.0135 0.2548 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 201 | 16 24 0.0003 0.0059 0.068 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 202 | 17 18 0.0007 0.0082 0.1319 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 203 | 17 27 0.0013 0.0173 0.3216 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 204 | 19 20 0.0007 0.0138 0.0 900.0 900.0 2500.0 1.06 0.0 1 -30.0 30.0; 205 | 19 33 0.0007 0.0142 0.0 900.0 900.0 2500.0 1.07 0.0 1 -30.0 30.0; 206 | 20 34 0.0009 0.018 0.0 900.0 900.0 2500.0 1.009 0.0 1 -30.0 30.0; 207 | 21 22 0.0008 0.014 0.2565 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0; 208 | 22 23 0.0006 0.0096 0.1846 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 209 | 22 35 0.0 0.0143 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0; 210 | 23 24 0.0022 0.035 0.361 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 211 | 23 36 0.0005 0.0272 0.0 900.0 900.0 2500.0 0.0 0.0 1 -30.0 30.0; 212 | 25 26 0.0032 0.0323 0.531 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 213 | 25 37 0.0006 0.0232 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0; 214 | 26 27 0.0014 0.0147 0.2396 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 215 | 26 28 0.0043 0.0474 0.7802 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 216 | 26 29 0.0057 0.0625 1.029 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 217 | 28 29 0.0014 0.0151 0.249 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0; 218 | 29 38 0.0008 0.0156 0.0 1200.0 1200.0 2500.0 1.025 0.0 1 -30.0 30.0; 219 | ]; 220 | 221 | %% dc grid topology 222 | %colunm_names% dcpoles 223 | mpc.dcpol=2; 224 | % numbers of poles (1=monopolar grid, 2=bipolar grid) 225 | %% bus data 226 | %column_names% busdc_i grid Pdc Vdc basekVdc Vdcmax Vdcmin Cdc 227 | mpc.dcbus = [ 228 | 1 1 0 1 345 1.1 0.9 0; 229 | 2 1 0 1 345 1.1 0.9 0; 230 | 3 1 0 1 345 1.1 0.9 0; 231 | 4 1 0 1 345 1.1 0.9 0; 232 | 5 1 0 1 345 1.1 0.9 0; 233 | 6 1 0 1 345 1.1 0.9 0; 234 | 7 1 0 1 345 1.1 0.9 0; 235 | 8 1 0 1 345 1.1 0.9 0; 236 | 9 1 0 1 345 1.1 0.9 0; 237 | 10 1 0 1 345 1.1 0.9 0; 238 | ]; 239 | 240 | % %% converters 241 | %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 242 | mpc.dcconv = [ 243 | 1 2 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 244 | 2 9 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 245 | 3 10 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 246 | 4 18 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 247 | 5 26 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 248 | 6 29 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 249 | 7 24 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 250 | 8 14 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 251 | 9 23 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 252 | 10 13 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.1033 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 100 -100; 253 | ]; 254 | 255 | %% branches 256 | %column_names% fbusdc tbusdc r l c rateA rateB rateC status 257 | mpc.dcbranch = [ 258 | 1 2 0.01 0 0 100 100 100 1; 259 | 2 3 0.01 0 0 100 100 100 1; 260 | 1 4 0.01 0 0 100 100 100 1; 261 | 2 4 0.01 0 0 100 100 100 1; 262 | 2 4 0.01 0 0 100 100 100 1; 263 | 1 5 0.01 0 0 100 100 100 1; 264 | 5 6 0.01 0 0 100 100 100 1; 265 | 5 7 0.01 0 0 100 100 100 1; 266 | 7 4 0.01 0 0 100 100 100 1; 267 | 4 8 0.01 0 0 100 100 100 1; 268 | 8 9 0.01 0 0 100 100 100 1; 269 | 8 10 0.01 0 0 100 100 100 1; 270 | ]; 271 | 272 | % INFO : === Translation Options === 273 | % INFO : Phase Angle Bound: 30.0 (deg.) 274 | % INFO : Gen Active Cost Model: stat 275 | % INFO : Setting Flat Start 276 | % INFO : 277 | % INFO : === Generator Classification Notes === 278 | % INFO : NUC 1 - 3.97 279 | % INFO : COW 9 - 96.03 280 | % INFO : 281 | % INFO : === Generator Active Cost Stat Model Notes === 282 | % INFO : Updated Generator Cost: NUC - 0.2 0.3 0.01 -> 0 6.72477811338 0 283 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 14.7076252515 0 284 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 24.8047338223 0 285 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 34.84464286 0 286 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 24.6529938713 0 287 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 32.306483114 0 288 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 18.1574766212 0 289 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 31.5501806971 0 290 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 22.5031675377 0 291 | % INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 27.4344440931 0 292 | % INFO : 293 | % INFO : === Generator Bounds Update Notes === 294 | % INFO : 295 | % INFO : === Base KV Replacement Notes === 296 | % INFO : 297 | % INFO : === Transformer Setting Replacement Notes === 298 | % WARNING : Transformer 23-36 connects the same voltage levels (345.0, 345.0) and has no phase shift, changing tap ratio 1.0 => 0.0 299 | % INFO : 300 | % INFO : === Line Capacity Monotonicity Notes === 301 | % INFO : 302 | % INFO : === Voltage Setpoint Replacement Notes === 303 | % INFO : Bus 1 : V=1.0393836, theta=-13.536602 -> V=1.0, theta=0.0 304 | % INFO : Bus 2 : V=1.0484941, theta=-9.7852666 -> V=1.0, theta=0.0 305 | % INFO : Bus 3 : V=1.0307077, theta=-12.276384 -> V=1.0, theta=0.0 306 | % INFO : Bus 4 : V=1.00446, theta=-12.626734 -> V=1.0, theta=0.0 307 | % INFO : Bus 5 : V=1.0060063, theta=-11.192339 -> V=1.0, theta=0.0 308 | % INFO : Bus 6 : V=1.0082256, theta=-10.40833 -> V=1.0, theta=0.0 309 | % INFO : Bus 7 : V=0.99839728, theta=-12.755626 -> V=1.0, theta=0.0 310 | % INFO : Bus 8 : V=0.99787232, theta=-13.335844 -> V=1.0, theta=0.0 311 | % INFO : Bus 9 : V=1.038332, theta=-14.178442 -> V=1.0, theta=0.0 312 | % INFO : Bus 10 : V=1.0178431, theta=-8.170875 -> V=1.0, theta=0.0 313 | % INFO : Bus 11 : V=1.0133858, theta=-8.9369663 -> V=1.0, theta=0.0 314 | % INFO : Bus 12 : V=1.000815, theta=-8.9988236 -> V=1.0, theta=0.0 315 | % INFO : Bus 13 : V=1.014923, theta=-8.9299272 -> V=1.0, theta=0.0 316 | % INFO : Bus 14 : V=1.012319, theta=-10.715295 -> V=1.0, theta=0.0 317 | % INFO : Bus 15 : V=1.0161854, theta=-11.345399 -> V=1.0, theta=0.0 318 | % INFO : Bus 16 : V=1.0325203, theta=-10.033348 -> V=1.0, theta=0.0 319 | % INFO : Bus 17 : V=1.0342365, theta=-11.116436 -> V=1.0, theta=0.0 320 | % INFO : Bus 18 : V=1.0315726, theta=-11.986168 -> V=1.0, theta=0.0 321 | % INFO : Bus 19 : V=1.0501068, theta=-5.4100729 -> V=1.0, theta=0.0 322 | % INFO : Bus 20 : V=0.99101054, theta=-6.8211783 -> V=1.0, theta=0.0 323 | % INFO : Bus 21 : V=1.0323192, theta=-7.6287461 -> V=1.0, theta=0.0 324 | % INFO : Bus 22 : V=1.0501427, theta=-3.1831199 -> V=1.0, theta=0.0 325 | % INFO : Bus 23 : V=1.0451451, theta=-3.3812763 -> V=1.0, theta=0.0 326 | % INFO : Bus 24 : V=1.038001, theta=-9.9137585 -> V=1.0, theta=0.0 327 | % INFO : Bus 25 : V=1.0576827, theta=-8.3692354 -> V=1.0, theta=0.0 328 | % INFO : Bus 26 : V=1.0525613, theta=-9.4387696 -> V=1.0, theta=0.0 329 | % INFO : Bus 27 : V=1.0383449, theta=-11.362152 -> V=1.0, theta=0.0 330 | % INFO : Bus 28 : V=1.0503737, theta=-5.9283592 -> V=1.0, theta=0.0 331 | % INFO : Bus 29 : V=1.0501149, theta=-3.1698741 -> V=1.0, theta=0.0 332 | % INFO : Bus 30 : V=1.0499, theta=-7.3704746 -> V=1.0, theta=0.0 333 | % INFO : Bus 31 : V=0.982, theta=0.0 -> V=1.0, theta=0.0 334 | % INFO : Bus 32 : V=0.9841, theta=-0.1884374 -> V=1.0, theta=0.0 335 | % INFO : Bus 33 : V=0.9972, theta=-0.19317445 -> V=1.0, theta=0.0 336 | % INFO : Bus 34 : V=1.0123, theta=-1.631119 -> V=1.0, theta=0.0 337 | % INFO : Bus 35 : V=1.0494, theta=1.7765069 -> V=1.0, theta=0.0 338 | % INFO : Bus 36 : V=1.0636, theta=4.4684374 -> V=1.0, theta=0.0 339 | % INFO : Bus 37 : V=1.0275, theta=-1.5828988 -> V=1.0, theta=0.0 340 | % INFO : Bus 38 : V=1.0265, theta=3.8928177 -> V=1.0, theta=0.0 341 | % INFO : Bus 39 : V=1.03, theta=-14.535256 -> V=1.0, theta=0.0 342 | % INFO : 343 | % INFO : === Generator Setpoint Replacement Notes === 344 | % INFO : Gen at bus 30 : Pg=250.0, Qg=161.762 -> Pg=520.0, Qg=270.0 345 | % INFO : Gen at bus 30 : Vg=1.0499 -> Vg=1.0 346 | % INFO : Gen at bus 31 : Pg=677.871, Qg=221.574 -> Pg=323.0, Qg=100.0 347 | % INFO : Gen at bus 31 : Vg=0.982 -> Vg=1.0 348 | % INFO : Gen at bus 32 : Pg=650.0, Qg=206.965 -> Pg=362.5, Qg=225.0 349 | % INFO : Gen at bus 32 : Vg=0.9841 -> Vg=1.0 350 | % INFO : Gen at bus 33 : Pg=632.0, Qg=108.293 -> Pg=326.0, Qg=125.0 351 | % INFO : Gen at bus 33 : Vg=0.9972 -> Vg=1.0 352 | % INFO : Gen at bus 34 : Pg=508.0, Qg=166.688 -> Pg=254.0, Qg=83.5 353 | % INFO : Gen at bus 34 : Vg=1.0123 -> Vg=1.0 354 | % INFO : Gen at bus 35 : Pg=650.0, Qg=210.661 -> Pg=343.5, Qg=100.0 355 | % INFO : Gen at bus 35 : Vg=1.0494 -> Vg=1.0 356 | % INFO : Gen at bus 36 : Pg=560.0, Qg=100.165 -> Pg=290.0, Qg=120.0 357 | % INFO : Gen at bus 36 : Vg=1.0636 -> Vg=1.0 358 | % INFO : Gen at bus 37 : Pg=540.0, Qg=-1.36945 -> Pg=282.0, Qg=125.0 359 | % INFO : Gen at bus 37 : Vg=1.0275 -> Vg=1.0 360 | % INFO : Gen at bus 38 : Pg=830.0, Qg=21.7327 -> Pg=432.5, Qg=75.0 361 | % INFO : Gen at bus 38 : Vg=1.0265 -> Vg=1.0 362 | % INFO : Gen at bus 39 : Pg=1000.0, Qg=78.4674 -> Pg=550.0, Qg=100.0 363 | % INFO : Gen at bus 39 : Vg=1.03 -> Vg=1.0 364 | % INFO : 365 | % INFO : === Writing Matpower Case File Notes === -------------------------------------------------------------------------------- /case5_3_he.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | %%%% %%%%% 3 | %%%% IEEE PES Power Grid Library - Optimal Power Flow - v23.09 %%%%% 4 | %%%% (https://github.com/power-grid-lib/pglib-opf) %%%%% 5 | %%%% Benchmark Group - Typical Operations %%%%% 6 | %%%% 10 - May - 2019 %%%%% 7 | %%%% %%%%% 8 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 | % 10 | % CASE5_3 Power flow data for modified 5 bus, 5 gen, 3 bus dc case based on PJM 5-bus system 11 | % Please see CASEFORMAT for details on the case file format. 12 | % 13 | % Based on data from ... 14 | % F.Li and R.Bo, "Small Test Systems for Power System Economic Studies", 15 | % Proceedings of the 2010 IEEE Power & Energy Society General Meeting 16 | % and dc system form ... 17 | % J. Beerten, D. Van Hertem and R. Belmans, "VSC MTDC systems with a distributed DC voltage control - A power flow approach," 18 | % 2011 IEEE Trondheim PowerTech, Trondheim, 2011, pp. 1-6. doi: 10.1109/PTC.2011.6019434 19 | % 20 | % Created by Hakan Ergun in 2019. 21 | % 22 | % Copyright (c) 2010 by The Institute of Electrical and Electronics Engineers (IEEE) 23 | % Licensed under the Creative Commons Attribution 4.0 24 | % International license, http://creativecommons.org/licenses/by/4.0/ 25 | % 26 | % Contact M.E. Brennan (me.brennan@ieee.org) for inquries on further reuse of 27 | % this dataset. 28 | 29 | function mpc = case5_3_he() 30 | mpc.version = '2'; 31 | mpc.baseMVA = 100.0; 32 | 33 | %% area data 34 | % area refbus 35 | mpc.areas = [ 36 | 1 4; 37 | ]; 38 | 39 | %% bus data 40 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 41 | mpc.bus = [ 42 | 1 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000; 43 | 2 1 300.0 98.61 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000; 44 | 3 2 300.0 98.61 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000; 45 | 4 3 400.0 131.47 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000; 46 | 5 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000; 47 | ]; 48 | 49 | %% generator data 50 | % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 51 | mpc.gen = [ 52 | 1 20.0 0.0 30.0 -30.0 1.0 100.0 1 40.0 0.0; 53 | 1 85.0 0.0 127.5 -127.5 1.0 100.0 1 170.0 0.0; 54 | 3 260.0 0.0 390.0 -390.0 1.0 100.0 1 520.0 0.0; 55 | 4 100.0 0.0 150.0 -150.0 1.0 100.0 1 200.0 0.0; 56 | 5 300.0 0.0 450.0 -450.0 1.0 100.0 1 600.0 0.0; 57 | ]; 58 | 59 | %% generator cost data 60 | % 2 startup shutdown n c(n-1) ... c0 61 | mpc.gencost = [ 62 | 2 0.0 0.0 3 0.000000 14.000000 0.000000; 63 | 2 0.0 0.0 3 0.000000 15.000000 0.000000; 64 | 2 0.0 0.0 3 0.000000 30.000000 0.000000; 65 | 2 0.0 0.0 3 0.000000 40.000000 0.000000; 66 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 67 | ]; 68 | 69 | %% branch data 70 | % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 71 | mpc.branch = [ 72 | 1 2 0.00281 0.0281 0.00712 400.0 400.0 400.0 0.0 0.0 1 -30.0 30.0; 73 | 1 4 0.00304 0.0304 0.00658 426 426 426 0.0 0.0 1 -30.0 30.0; 74 | 1 5 0.00064 0.0064 0.03126 426 426 426 0.0 0.0 1 -30.0 30.0; 75 | 2 3 0.00108 0.0108 0.01852 426 426 426 0.0 0.0 1 -30.0 30.0; 76 | 3 4 0.00297 0.0297 0.00674 426 426 426 0.0 0.0 1 -30.0 30.0; 77 | 4 5 0.00297 0.0297 0.00674 240.0 240.0 240.0 0.0 0.0 1 -30.0 30.0; 78 | ]; 79 | 80 | 81 | %% dc grid topology 82 | %colunm_names% dcpoles 83 | mpc.dcpol=2; 84 | % numbers of poles (1=monopolar grid, 2=bipolar grid) 85 | %% bus data 86 | %column_names% busdc_i grid Pdc Vdc basekVdc Vdcmax Vdcmin Cdc 87 | mpc.dcbus = [ 88 | 1 1 0 1 345 1.1 0.9 0; 89 | 2 1 0 1 345 1.1 0.9 0; 90 | 3 1 0 1 345 1.1 0.9 0; 91 | ]; 92 | 93 | %% converters 94 | %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 95 | mpc.dcconv = [ 96 | 1 2 1 1 -60 -40 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -58.6274 1.0079 0 100 -100 50 -50; 97 | 2 3 2 1 0 0 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0070 21.9013 1.0000 0 100 -100 50 -50; 98 | 3 5 1 1 35 5 0 1 0.0015 0.1121 1 1 0.0887 1 0.0001 0.16428 1 345 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 36.1856 0.9978 0 100 -100 50 -50; 99 | ]; 100 | 101 | %% branches 102 | %column_names% fbusdc tbusdc r l c rateA rateB rateC status 103 | mpc.dcbranch = [ 104 | 1 2 0.052 0 0 100 100 100 1; 105 | 2 3 0.052 0 0 100 100 100 1; 106 | 1 3 0.073 0 0 100 100 100 1; 107 | ]; 108 | 109 | % INFO : === Translation Options === 110 | % INFO : Phase Angle Bound: 30.0 (deg.) 111 | % INFO : Line Capacity Model: stat 112 | % INFO : Setting Flat Start 113 | % INFO : Line Capacity PAB: 15.0 (deg.) 114 | % INFO : 115 | % INFO : === Generator Bounds Update Notes === 116 | % INFO : 117 | % INFO : === Base KV Replacement Notes === 118 | % INFO : 119 | % INFO : === Transformer Setting Replacement Notes === 120 | % INFO : 121 | % INFO : === Line Capacity Stat Model Notes === 122 | % INFO : Updated Thermal Rating: on line 1-4 : Rate A, Rate B, Rate C , 9900.0, 0.0, 0.0 -> 426 123 | % INFO : Updated Thermal Rating: on line 1-5 : Rate A, Rate B, Rate C , 9900.0, 0.0, 0.0 -> 426 124 | % INFO : Updated Thermal Rating: on line 2-3 : Rate A, Rate B, Rate C , 9900.0, 0.0, 0.0 -> 426 125 | % INFO : Updated Thermal Rating: on line 3-4 : Rate A, Rate B, Rate C , 9900.0, 0.0, 0.0 -> 426 126 | % INFO : 127 | % INFO : === Line Capacity Monotonicity Notes === 128 | % INFO : 129 | % INFO : === Voltage Setpoint Replacement Notes === 130 | % INFO : Bus 1 : V=1.0, theta=0.0 -> V=1.0, theta=0.0 131 | % INFO : Bus 2 : V=1.0, theta=0.0 -> V=1.0, theta=0.0 132 | % INFO : Bus 3 : V=1.0, theta=0.0 -> V=1.0, theta=0.0 133 | % INFO : Bus 4 : V=1.0, theta=0.0 -> V=1.0, theta=0.0 134 | % INFO : Bus 5 : V=1.0, theta=0.0 -> V=1.0, theta=0.0 135 | % INFO : 136 | % INFO : === Generator Setpoint Replacement Notes === 137 | % INFO : Gen at bus 1 : Pg=40.0, Qg=0.0 -> Pg=20.0, Qg=0.0 138 | % INFO : Gen at bus 1 : Vg=1.0 -> Vg=1.0 139 | % INFO : Gen at bus 1 : Pg=170.0, Qg=0.0 -> Pg=85.0, Qg=0.0 140 | % INFO : Gen at bus 1 : Vg=1.0 -> Vg=1.0 141 | % INFO : Gen at bus 3 : Pg=323.49, Qg=0.0 -> Pg=260.0, Qg=0.0 142 | % INFO : Gen at bus 3 : Vg=1.0 -> Vg=1.0 143 | % INFO : Gen at bus 4 : Pg=0.0, Qg=0.0 -> Pg=100.0, Qg=0.0 144 | % INFO : Gen at bus 4 : Vg=1.0 -> Vg=1.0 145 | % INFO : Gen at bus 5 : Pg=466.51, Qg=0.0 -> Pg=300.0, Qg=0.0 146 | % INFO : Gen at bus 5 : Vg=1.0 -> Vg=1.0 147 | % INFO : 148 | % INFO : === Writing Matpower Case File Notes === 149 | -------------------------------------------------------------------------------- /case67.m: -------------------------------------------------------------------------------- 1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 | %%%% %%%%% 3 | %%%% IEEE PES Power Grid Library - Optimal Power Flow with HVDC Lines - v23.09 %%%%% 4 | %%%% (https://github.com/power-grid-lib/pglib-opf-hvdc) %%%%% 5 | %%%% Benchmark Group - Typical Operations %%%%% 6 | %%%% 23 - August - 2023 %%%%% 7 | %%%% %%%%% 8 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9 | % 10 | % AC/DC grid OPF test case based on: 11 | % Sass, F., Sennewald, T., Marten, A.-K. and Westermann, D. (2017), 12 | % Mixed AC high-voltage direct current benchmark test system for security constrained optimal power flow calculation. 13 | % IET Gener. Transm. Distrib., 11: 447-455. https://doi.org/10.1049/iet-gtd.2016.0993 14 | % 15 | % Copyright (c) 2023 by Hakan Ergun, Vaishally Bhardwaj ({hakan.ergun, vaishally.bhardwaj}@kuleuven.be) 16 | % Licensed under the Creative Commons Attribution 4.0 17 | % International license, http://creativecommons.org/licenses/by/4.0/ 18 | % 19 | 20 | function mpc = case67 21 | mpc.version = '2'; 22 | mpc.baseMVA = 100.0; 23 | 24 | %% bus data 25 | % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 26 | mpc.bus = [ 27 | 1 3 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 28 | 2 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 29 | 3 1 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 30 | 4 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 31 | 5 2 0.0 0.0 0.0 0.0 1 1.0526 0 380.0 1 1.10000 0.90000; 32 | 6 1 191.0 76.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 33 | 7 1 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 34 | 8 1 287.0 73.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 35 | 9 1 186.0 74.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 36 | 10 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 37 | 11 1 271.0 55.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 38 | 12 1 171.0 87.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 39 | 13 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 40 | 14 1 199.0 60.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 41 | 15 1 113.0 52.5 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 42 | 16 1 38 7.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 43 | 17 1 275.0 106.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 44 | 18 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 45 | 19 1 165.0 46.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 46 | 20 1 178.0 82.5 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 47 | 21 1 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 48 | 22 1 30.0 7.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 49 | 23 1 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 50 | 24 1 32.0 7.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 51 | 25 2 0.0 0.0 0.0 0.0 1 1.0 0 380.0 1 1.10000 0.90000; 52 | 26 1 395.0 89.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 53 | 27 1 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 54 | 28 1 665.0 99.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 55 | 29 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 56 | 30 1 266.0 100.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 57 | 31 1 845.0 119.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 58 | 32 1 332.0 137.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 59 | 33 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 60 | 34 1 540.0 158.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 61 | 35 1 460.0 97.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 62 | 36 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 63 | 37 1 451.0 190.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 64 | 38 1 150.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 65 | 39 1 629.0 87.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 66 | 40 1 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 67 | 41 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 68 | 42 1 859.0 180.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 69 | 43 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 70 | 44 1 474.0 92.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 71 | 45 1 668.0 109.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 72 | 46 1 614.0 95.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 73 | 47 1 81.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 74 | 48 1 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 75 | 49 1 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 76 | 50 2 0.0 0.0 0.0 0.0 2 1.0 0 380.0 1 1.10000 0.90000; 77 | 51 1 430.0 123.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 78 | 52 1 309.0 102.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 79 | 53 1 100.0 30.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 80 | 54 1 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 81 | 55 1 303.0 110.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 82 | 56 2 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 83 | 57 1 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 84 | 58 1 324.0 157.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 85 | 59 2 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 86 | 60 1 115.0 42.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 87 | 61 1 187.0 75.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 88 | 62 1 319.0 95.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 89 | 63 2 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 90 | 64 2 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 91 | 65 1 315.0 97.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 92 | 66 2 0.0 0.0 0.0 0.0 3 1.0 0 380.0 1 1.10000 0.90000; 93 | 67 2 0.0 0.0 0.0 0.0 4 1.0 0 380.0 1 1.10000 0.90000; 94 | 95 | ]; 96 | 97 | %% generator data 98 | % 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 alpha 99 | mpc.gen = [ 100 | 1 700.0 23.0 1000.0 -500.0 1.0526 100.0 1 1000.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 101 | 2 1500.0 100.0 0.0 0.0 1.0526 100.0 1 1500.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 102 | 4 523.0 140.0 350.0 -350.0 1.0526 100.0 1 560.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 103 | 5 1200.0 100.0 000.0 000.0 1.0526 100.0 1 1200.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 104 | 10 436.0 105.0 350.0 -350.0 1.0526 100.0 1 560.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 105 | 13 541.0 117.0 300.0 -300.0 1.0526 100.0 1 630.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 106 | 18 681.0 -35.0 400.0 -400.0 1.0263 100.0 1 720.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 107 | 25 469.0 59.0 250.0 -250.0 1.0395 100.0 1 560.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 108 | 29 500.0 101.0 350.0 -350.0 1.0263 100.0 1 630.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 109 | 33 496.0 306.0 500.0 -500.0 1.0263 100.0 1 850.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 110 | 36 512.0 249.0 400.0 -400.0 1.0263 100.0 1 720.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 111 | 41 350.0 238.0 450.0 -450.0 1.0395 100.0 1 850.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 112 | 43 574.0 223.0 500.0 -250.0 1.0263 100.0 1 720.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 113 | 50 581.0 150.0 400.0 -400.0 1.0395 100.0 1 720.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 114 | 56 496.0 56.0 250.0 -250.0 1.0395 100.0 1 560.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 115 | 59 431.0 206.0 350.0 -350.0 1.0447 100.0 1 720.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 116 | 63 488.0 152.0 250.0 -300.0 1.0447 100.0 1 520.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 117 | 64 300.0 61.0 250.0 -400.0 1.0474 100.0 1 560.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 118 | 66 537.0 143.0 300.0 -400.0 1.0447 100.0 1 630.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 119 | 67 800.0 0.0 0.0 0.0 1.0526 100.0 1 800.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 120 | 121 | ]; 122 | 123 | %% generator data 124 | % 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 alpha 125 | %mpc.gen = [ 126 | % 1 700.0 23.0 1000.0 -500.0 1.0526 100.0 1 1000.0 400.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 127 | % 2 1500.0 100.0 0.0 0.0 1.0526 100.0 1 1500.0 1500.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 128 | % 4 523.0 140.0 350.0 -350.0 1.0526 100.0 1 560.0 220.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 129 | % 5 1200.0 100.0 000.0 000.0 1.0526 100.0 1 1200.0 1200.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 130 | % 10 436.0 105.0 350.0 -350.0 1.0526 100.0 1 560.0 220.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 131 | % 13 541.0 117.0 300.0 -300.0 1.0526 100.0 1 630.0 250.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 132 | % 18 681.0 -35.0 400.0 -400.0 1.0263 100.0 1 720.0 300.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 133 | % 25 469.0 59.0 250.0 -250.0 1.0395 100.0 1 560.0 220.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 134 | % 29 500.0 101.0 350.0 -350.0 1.0263 100.0 1 630.0 250.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 135 | % 33 496.0 306.0 500.0 -500.0 1.0263 100.0 1 850.0 350.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 136 | % 36 512.0 249.0 400.0 -400.0 1.0263 100.0 1 720.0 300.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 137 | % 41 350.0 238.0 450.0 -450.0 1.0395 100.0 1 850.0 350.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 138 | % 43 574.0 223.0 500.0 -250.0 1.0263 100.0 1 720.0 220.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 139 | % 50 581.0 150.0 400.0 -400.0 1.0395 100.0 1 720.0 300.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 140 | % 56 496.0 56.0 250.0 -250.0 1.0395 100.0 1 560.0 220.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 141 | % 59 431.0 206.0 350.0 -350.0 1.0447 100.0 1 720.0 300.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 142 | % 63 488.0 152.0 250.0 -300.0 1.0447 100.0 1 520.0 250.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 143 | % 64 300.0 61.0 250.0 -400.0 1.0474 100.0 1 560.0 250.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 144 | % 66 537.0 143.0 300.0 -400.0 1.0447 100.0 1 630.0 250.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 145 | % 67 800.0 0.0 0.0 0.0 1.0526 100.0 1 800.0 800.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1; 146 | % 147 | %]; 148 | 149 | %% branch data 150 | % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 151 | mpc.branch = [ 152 | 1 5 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 153 | 1 7 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 154 | 1 8 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 155 | 1 14 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 156 | 2 3 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 157 | 2 9 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 158 | 2 12 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 159 | 3 4 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 160 | 3 10 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 161 | 3 12 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 162 | 3 9 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 163 | 4 14 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 164 | 4 19 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 165 | 5 6 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 166 | 5 7 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 167 | 5 8 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 168 | 6 7 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 169 | 7 15 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 170 | 7 16 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 171 | 8 9 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 172 | 10 11 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 173 | 10 22 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 174 | 11 12 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 175 | 11 13 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 176 | 12 13 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 177 | 13 53 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 178 | 14 15 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 179 | 14 18 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 180 | 16 17 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 181 | 16 18 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 182 | 17 24 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 183 | 18 24 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 184 | 19 20 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 185 | 19 23 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 186 | 20 21 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 187 | 21 25 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 188 | 21 22 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 189 | 21 23 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 190 | 22 25 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 191 | 18 20 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 192 | 24 49 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 193 | 25 43 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 194 | 26 27 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 195 | 26 31 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 196 | 26 40 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 197 | 27 28 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 198 | 28 35 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 199 | 28 37 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 200 | 29 39 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 201 | 29 44 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 202 | 30 31 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 203 | 31 27 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 204 | 30 26 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 205 | 32 40 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 206 | 41 40 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 207 | 43 44 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 208 | 33 51 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 209 | 33 34 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 210 | 34 51 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 211 | 35 33 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 212 | 35 36 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 213 | 35 47 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 214 | 29 35 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 215 | 36 37 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 216 | 36 38 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 217 | 37 38 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 218 | 39 40 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 219 | 39 43 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 220 | 41 42 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 221 | 42 43 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 222 | 42 49 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 223 | 43 49 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 224 | 44 45 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 225 | 44 48 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 226 | 45 46 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 227 | 45 50 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 228 | 47 48 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 229 | 46 48 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 230 | 47 50 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 231 | 47 51 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 232 | 47 59 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 233 | 52 53 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 234 | 52 54 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 235 | 63 55 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 236 | 22 56 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 237 | 54 65 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 238 | 55 57 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 239 | 58 61 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 240 | 56 59 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 241 | 57 58 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 242 | 56 58 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 243 | 58 60 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 244 | 62 66 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 245 | 61 62 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 246 | 52 64 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 247 | 62 63 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 248 | 59 60 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 249 | 63 57 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 250 | 65 66 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 251 | 66 54 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 252 | 66 64 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 253 | 30 32 0.00207756 0.01800554 0.61242207 900.0 900.0 900.0 0.0 0.0 1 -60.0 60.0; 254 | ]; 255 | 256 | %% dc grid topology 257 | %colunm_names% dcpoles 258 | mpc.dcpol=2; 259 | % numbers of poles (1=monopolar grid, 2=bipolar grid) 260 | %% bus data 261 | %column_names% busdc_i grid Pdc Vdc basekVdc Vdcmax Vdcmin Cdc area 262 | mpc.busdc = [ 263 | 1 1 0 1 500 1.1 0.9 0 1; 264 | 2 1 0 1 500 1.1 0.9 0 2; 265 | 3 1 0 1 500 1.1 0.9 0 1; 266 | 4 1 0 1 500 1.1 0.9 0 1; 267 | 5 1 0 1 500 1.1 0.9 0 2; 268 | 6 1 0 1 500 1.1 0.9 0 3; 269 | 7 1 0 1 500 1.1 0.9 0 3; 270 | 8 1 0 1 500 1.1 0.9 0 2; 271 | 9 1 0 1 500 1.1 0.9 0 4; 272 | ]; 273 | 274 | %% converters 275 | %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 276 | mpc.convdc = [ 277 | 1 7 2 1 -577.5 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -465.9871 0.9999 0 2000 -2000 1000 -1000; 278 | 2 40 3 1 1000 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 500.0000 0.9947 0 2000 -2000 1000 -1000; 279 | 3 3 3 1 -550 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -517.1051 1.0022 0 2000 -2000 1000 -1000; 280 | 4 23 3 1 -600 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -560.8855 0.9972 0 2000 -2000 1000 -1000; 281 | 5 48 3 1 1000 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 1103.897 0.9921 0 2000 -2000 1000 -1000; 282 | 6 54 3 1 50 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 51.31034 0.9963 0 2000 -2000 1000 -1000; 283 | 7 57 3 1 -550 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -516.4622 0.9957 0 2000 -2000 1000 -1000; 284 | 8 27 3 1 1000 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 1111.644 0.9900 0 2000 -2000 1000 -1000; 285 | 9 67 1 2 -800 0 0 1 0.01 0.01 1 1 0.01 0 0.01 0.01 0 500 1.1 0.9 1.1 1 1.103 0.887 2.885 2.885 0.0050 -731.3881 1.0066 0 1000 -1000 1000 -1000; 286 | ]; 287 | 288 | %% branches 289 | %column_names% fbusdc tbusdc r l c rateA rateB rateC status 290 | mpc.branchdc = [ 291 | 1 2 0.0012 0 0 1575 1575 1575 1; 292 | 3 4 0.0012 0 0 1575 1575 1575 1; 293 | 4 5 0.0012 0 0 1575 1575 1575 1; 294 | 6 7 0.0012 0 0 1575 1575 1575 1; 295 | 1 3 0.0012 0 0 1575 1575 1575 1; 296 | 2 8 0.0012 0 0 1575 1575 1575 1; 297 | 8 5 0.0012 0 0 1575 1575 1575 1; 298 | 4 6 0.0012 0 0 1575 1575 1575 1; 299 | 2 4 0.0012 0 0 1575 1575 1575 1; 300 | 5 7 0.0012 0 0 1575 1575 1575 1; 301 | 3 9 0.0012 0 0 1575 1575 1575 1; 302 | ]; 303 | 304 | %% generator cost data 305 | % 2 startup shutdown n c(n-1) ... c0 306 | mpc.gencost = [ 307 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 308 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 309 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 310 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 311 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 312 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 313 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 314 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 315 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 316 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 317 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 318 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 319 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 320 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 321 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 322 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 323 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 324 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 325 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 326 | 2 0.0 0.0 3 0.000000 10.000000 0.000000; 327 | 328 | ]; 329 | 330 | --------------------------------------------------------------------------------