How to Perform Load Flow

Analysis of IEEE-9 Bus Test

System using MATPOWER

By: Dr. J. A. Laghari

 Load Flow Analysis of IEEE-9 Bus System Using

Newton
Gauss-Seidel
Raphson
Method
Method

Fast
Decoupled
Method
How to Download &
Install MATPOWER
Toolbox Completely
Free
By: Dr. J. A. Laghari
Built-in 59 Cases of Various Test
Systems in MATPOWER Toolbox
 Built-in Case List in MATPOWER Toolbox
 All of the functions (cases) listed in MATPOWER Toolbox can be found in
<MATPOWER>/data.
Table D-18: Small Transmission System Test Cases
Case Name Description
case4gs 4-bus example case from Grainger & Stevenson

case5 modified 5-bus PJM example case from Rui Bo

case6ww 6-bus example case from Wood & Wollenberg

case9 9-bus example case from Chow

case9Q case9 with reactive power costs

Source: MATPOWER User Manual 7.0, Page 188

 Built-in Case List in MATPOWER Toolbox
Table D-18: Small Transmission System Test Cases
Case Name Description
modified case9, target for example continuation power flow
case9target
case14 IEEE 14-bus case

case24_ieee_rts IEEE RTS 24-bus case

case30 30-bus case, based on IEEE 30-bus case

case_ieee30 IEEE 30-bus case

case30pwl case30 with piecewise linear costs

case30Q case30 with reactive power costs

Source: MATPOWER User Manual 7.0, Page 188

 Built-in Case List in MATPOWER Toolbox
Table D-18: Small Transmission System Test Cases
Case Name Description
case39 39-bus New England case

case57 IEEE 57-bus case

case_RTS_GMLC 96-machine, 73-bus Reliability Test System

case118 IEEE 118-bus case

case145 IEEE 145-bus case, 50 generator dynamic test

case
case300 IEEE 300-bus case

Source: MATPOWER User Manual 7.0, Page 188

 Built-in Case List in MATPOWER Toolbox
Table D-19: Small Radial Distribution System Test Cases
Case Name Description
case4_dist 4-bus example radial distribution system case

case18 18-bus radial distribution system from Grady,

Samotyj and Noyola
case22 22-bus radial distribution system from Raju,
Murthy and Ravindra
case33bw 33-bus radial distribution system from Baran
and Wu
case69 69-bus radial distribution system from Das

case85 85-bus radial distribution system from Das,

Kothari and Kalam
case141 141-bus radial distribution system from Khodr,
Olsina, De Jesus and Yusta
Source: MATPOWER User Manual 7.0, Page 188
 Built-in Case List in MATPOWER Toolbox
Table D-20: ACTIV Synthetic Grid Test Cases
Case Name Description
case_ACTIVSg200 200-bus Illinois synthetic model

case_ACTIVSg500 500-bus South Carolina synthetic model

case_ACTIVSg2000 2000-bus Texas synthetic model

case_ACTIVSg10k 10,000-bus US WECC synthetic model

case_ACTIVSg25k 25,000-bus US Northeast/Mid-Atlantic synthetic

model
case_ACTIVSg70k 70,000-bus Eastern US synthetic model

Case_SyntheticUSA 82,000-bus continental USA synthetic model

Source: MATPOWER User Manual 7.0, Page 189

 Built-in Case List in MATPOWER Toolbox
Table D-20: ACTIV Synthetic Grid Test Cases
Case Name Description
contab_ACTIVSg200 contingency table for case ACTIVSg200

contab_ACTIVSg500 contingency table for case ACTIVSg500

contab_ACTIVSg2000 contingency table for case ACTIVSg2000

contab_ACTIVSg10k contingency table for case ACTIVSg10k

scenarios_ACTIVSg200 one year of hourly zonal load scenarios for case

ACTIVSg200
scenarios_ACTIVSg2000 one year of hourly area load scenarios for case
ACTIVSg2000

Source: MATPOWER User Manual 7.0, Page 189

  Table D-21: Polish System Test Cases
Case Name Description
case2383wp Polish system - winter 1999-2000 peak

case2736sp Polish system - summer 2004 peak

case2737sop Polish system-summer 2004 off-peak

case2746wop Polish system - winter 2003-04 off-peak

case2746wp Polish system - winter 2003-04 evening peak

case3012wp Polish system - winter 2007-08 evening peak

case3120sp Polish system - summer 2008 morning peak

case3375wp Polish system plus - winter 2007-08 evening

peak
Source: MATPOWER User Manual 7.0, Page 189
 Built-in Case List in MATPOWER Toolbox
Table D-22: PEGASE European System Test Cases
Case Name Description
case89pegase 89-bus portion of European transmission
system from PEGASE project
case1354pegase 1354-bus portion of European transmission
system from PEGASE project
case2869pegase 2869-bus portion of European transmission
system from PEGASE project
case9241pegase 9241-bus portion of European transmission
system from PEGASE project*
case13659pegase 13659-bus portion of European transmission
system from PEGASE project

Source: MATPOWER User Manual 7.0, Page 190

  Table D-23: RTE French System Test Cases
Case Name Description
case1888rte 1888-bus snapshot of VHV French transmission
system from RTE
case1951rte 1951-bus snapshot of VHV French transmission
system from RTE
case2848rte 2848-bus snapshot of VHV French transmission
system from RTE
case2868rte 2868-bus snapshot of VHV French transmission
system from RTE
case6468rte 6468-bus snapshot of VHV and HV French
transmission system from RTE
case6470rte 6470-bus snapshot of VHV and HV French
transmission system from RTE
case6495rte 6495-bus snapshot of VHV and HV French
transmission system from RTE
case6515rte 6515-bus snapshot of VHV and HV French
transmission system from RTE
Source: MATPOWER User Manual 7.0, Page 190
  How to Load data of any case
 Syntax to Load a Case:
 mpc = loadcase(casefilename);
 Where mpc stands for matpower case.

 Syntax to Load IEEE-9 Bus Test System Case:

 The M-file name for IEEE-9 Bus test name is case9 as mentioned in the
previous Table. Therefore, its syntax will become:
 mpc = loadcase(case9);
 This will open the structure file of IEEE 9 bus test system with its all data.

Source: MATPOWER User Manual 7.0, Page 20

  How to Open M-file of a case:
 All of the functions (cases) listed in MATPOWER Toolbox can be found in
<MATPOWER>/data.
 Go to <MATPOWER>/data folder in MATLAB. Double click on any M-file of a
case. It will be opened.

Source: MATPOWER User Manual 7.0, Page 20

 Single Line Diagram of IEEE 9 Bus Test System:
Bus2 T2 Bus7 Bus8 Bus9 T3 Bus3
250 MVA 150 MVA
Gen 2 19 kV / 230 kV 13.8 kV / 230 kV

Gen 3

Bus5 Bus6

Bus4

T1
100 MVA
16.5 kV / 230 kV

Bus1

Gen1
  Line Data of
IEEE 9 Bus
Test System:
 It is also known as Western
System Coordinating
Council (WECC)
3-Machine, 9-bus system

Source: Power System Control and Stability

by P. M. Anderson, Chapter 2. Page No. 38
  Line Data
of IEEE 9
Bus Test
System:

Source: Power System Dynamics and

Stability by Peter W. Sauer and M. A. Pai,
Chapter 7. Example 7.1, Page No. 164-165.
  Bus Data of IEEE-9 Bus System
Bus No. Bus Code Voltage Generation Load
(P.U.) MW Mvar MW Mvar
1 3 (Slack) 1.04 0 0 0 0
2 2 (PV) 1.025 163 6.7 0 0
3 2 (PV) 1.025 85 -10.9 0 0
4 1 (PQ) 1 0 0 0 0
5 1 (PQ) 1 0 0 125 50
6 1 (PQ) 1 0 0 90 30
7 1 (PQ) 1 0 0 0 0
8 1 (PQ) 1 0 0 100 35
9 1 (PQ) 1 0 0 0 0
 In MATPOWER Toolbox, Bus Code 1 = PQ bus; Bus Code 2 = PV Bus
 Bus Code 3 = Slack Bus
  Line or Branch Data of IEEE-9 Bus System
Line Line
R X ½B B
From To
1 4 0 0.0576 0 0
4 5 0.01 0.085 0.088 0.176
4 6 0.017 0.092 0.079 0.158
6 9 0.039 0.17 0.179 0.358
5 7 0.032 0.161 0.153 0.306
9 3 0 0.0586 0 0
7 2 0 0.0625 0 0
9 8 0.0119 0.1008 0.1045 0.209
7 8 0.0085 0.072 0.0745 0.149
  How to Perform Load Flow Analysis in MATPOWER
 Syntax to Run the load flow analysis of a Case:
 runpf('casefilename');
 Syntax to Run Load Flow Analysis of IEEE9 Bus Test System Case:
 The M-file name for IEEE-9 Bus test name is case9 as mentioned in the
previous Table. Therefore, its syntax will become:
 runpf('case9');
 By default it will perform the Load Flow Analysis using Newton Raphson
Method.

Source: MATPOWER User Manual 7.0, Page 20

  Voltage, Power
Flows of IEEE
9 Bus Test
System:

Source: Power System Control and Stability

by P. M. Anderson, Chapter 2. Page No. 38
 Voltage, Power Flows of IEEE 9 Bus Test System:

Source: Power System Dynamics and

Stability by Peter W. Sauer and M. A. Pai,
Chapter 7. Example 7.1, Page No. 164-165.
  Load Flow Solution of IEEE-9 Bus System
Bus No. Generation Load Bus Voltage
MW Mvar MW Mvar Voltage(p.u.) Angle
1 71.6 27.0 0 0 1.04 0
2 163 6.7 0 0 1.025 9.3°
3 85 -10.9 0 0 1.025 4.7°
4 0 0 0 0 1.026 -2.2°
5 0 0 125 50 0.996 -4.0°
6 0 0 90 30 1.013 -3.7°
7 0 0 0 0 1.026 3.7°
8 0 0 100 35 1.016 0.7°
9 0 0 0 0 1.032 2.0°
 How to Perform Load Flow
Analysis using Gauss-Seidel Method
  How to Perform Load Flow Analysis in MATPOWER
 Syntax to Change the Type of Load Flow Method of a Case:
 mpopt = mpoption('pf.alg', 'LFName');
 results = runpf('casefilename', mpopt);
 Where mpopt = matpower option.
 Pf.alg = power flow analysis algorithm.
 Lfname = name of the load flow analysis method.

Source: MATPOWER User Manual 7.0, Page 20

  Table 4-2: Power Flow Options
Name Default Description

model 'AC' AC vs. DC modeling for power flow and OPF formulation
'AC' --- use AC formulation and corresponding alg options
'DC' ---- use DC formulation and corresponding alg options
pf.alg 'NR' AC Power Flow algorithm:
'NR' --- Newtons's method (formulation depends on values of
pf.current_balance and pf.v_cartesian options)
'NR-SC' --- Newtons's method (power mismatch, cartesian)
'NR-IP' --- Newtons's method (current mismatch, polar)
'NR-IC' --- Newtons's method (current mismatch, cartesian)
'FDXB' --- Fast-Decoupled (XB version)
'FDBX' --- Fast-Decouple (BX version)
'GS' --- Gauss-Seidel

Source: MATPOWER User Manual 7.0, Page 45

 How to Perform Load Flow Analysis in MATPOWER

 Syntax to Change the Load Flow Method to Gauss

Seidel Method for IEEE-9 Bus Test System:

 mpopt = mpoption('pf.alg', 'GS');

 results = runpf('case9', mpopt);
How to Perform Load Flow Analysis
using Fast Decoupled Method
 How to Perform Load Flow Analysis in MATPOWER

 Syntax to Change the Load Flow Method to Fast Decoupled

Method (BX Version) for IEEE-9 Bus Test System:

 mpopt = mpoption('pf.alg', 'FDBX');

 results = runpf('case9', mpopt);
 How to Perform Load Flow Analysis in MATPOWER

 Syntax to Change the Load Flow Method to Fast Decoupled

Method (XB Version) for IEEE-9 Bus Test System:

 mpopt = mpoption('pf.alg', 'FDXB');

 results = runpf('case9', mpopt);
How to Perform Load Flow Analysis
using Newton Raphson Method
 How to Perform Load Flow Analysis in MATPOWER

 Syntax to Change the Load Flow Method to Newton

Raphson Method for IEEE-9 Bus Test System:

 mpopt = mpoption('pf.alg', 'NR');

 results = runpf('case9', mpopt);
 Performance Comparison of Different Load Flow
Methods
THANK YOU VERY MUCH