PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
MUTAH UNIVERSITY
College of Engineering
Department of Electrical Engineering
Power System Simulation Lab
Dr. Khaled Alawasa
LAB No. 1 Load Flow Using MATLAB
1. Objective
In this experiment, you will have the opportunity to become familiar with MATLAB software
and to know how to prepare and conduct power flow analysis using this software.
2. Pre-Lab Preparation
1. Revise the power flow (Load flow) concept.
2. Refer to the Textbook Power System Analysis by Hadi Saadat. Go through Section
6.9 and Example 6.9 of the textbook to learn how a power flow solution is done using
MATLAB (and the provided m-files by Saadat’s textbook). Learn the format of input
data and MATLAB commands used for a power flow solution.
3. Lab Procedure and Practice
Task 1 (Will be explained during the lab):
1. Down load the MATLAB file that contain the developed MATLAB commands by Prof.
Hadi Saadat.
2. Run a sample case Example 6.9.
3. Explore how the power flow can be conducted.
Page 1 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
Loads: Lines:
On bus 2, 21.7 MW and 12.7 Mvar; Bus1 to Bus2, Series impedance = 0.0192+j0.0575 pu, shunt cap admittance = 0.0528 pu
On bus 3, 2.4 MW and 1.2 Mvar Bus1 to Bus3, Series impedance = 0.0452+j0.1852 pu, shunt cap admittance = 0.0408 pu
On bus 4, 7.6 MW and 1.6 Mvar …
…
Summary of steps to prepare the MATLAB code for a power flow solution:
Values for basemva, accuracy, accel, and maxiter must be specified.
The input matrices busdata and linedata must be prepared according to the following MATLAB format.
Then, the following MATLAB code will do a power flow solution by Gauss-Seidel.
Figure below shows the general structure and steps for MATLAB command :
Page 2 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
% lf_gauss1
clear
… Note Slack Bus code 1; Load Bus code 0; Voltage-Controlled Bus code 2.
…
basemva = 100;
accuracy = 0.001;
accel = 1.8; maxiter = 100;
lfybus
lfgauss
busout
lineflow
Page 3 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
Task #2
Write and run a MATLAB code to solve the power flow problem of Example 6.8 shown below and compare
your results with the textbook solution (Use three methods), and use these setup data.
basemva = 100; accuracy = 0.001; accel = 1.8; maxiter = 100;
1) Run the MATLAB code in Gauss-Seidel and compare your results with the textbook answers.
Describe the information given by the output of the powerflow program.
2) Change the power flow method to Newton-Raphson and note the differences observed between
the Gauss-Seidel and Newton-Raphson methods.
3) Change the power flow method to Fast-Decoupled and note the differences observed.
Page 4 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
Task #2 Results: Tabulate your results in the following Tables
A- Input Data – Bus Data
Bus Bus V- V- Load Generator Static Mvar
No. code Mag. Angle +Qc/-Ql
MW MVAR MW MVAR Qmin Qmax
B- Input Data – Line Data
Bus Bus Resistance Reactance Susceptance (B) tr. tap at bus nl
From To R nr p.u. X p.u. 1/2 B p.u. (> 1 or <)
Page 5 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
C- Power Flow results: Run the code for each power flow method and tabulate the
results:
1- Power Flow Solution by Gauss-Seidel Method
Maximum Power Mismatch = …….
No. of Iterations = …….
Bus V-Mag. V- Load Generator Injected
No. Angle
MW MVAR MW MVAR MVAR
Line Flow and Losses
Line Power at bus & line flow Line losses Transformer
From to MW Mvar MVA MW Mvar tap
MW Mvar
Total losses
Page 6 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
2- Power Flow Solution by Newton-Raphson Method
Maximum Power Mismatch = …….
No. of Iterations = …….
Bus V-Mag. V- Load Generator Injected
No. Angle
MW MVAR MW MVAR MVAR
Line Flow and Losses
Line Power at bus & line flow Line losses Transformer
From to MW Mvar MVA MW Mvar tap
MW Mvar
Total losses
Page 7 of 8
PSSL-LAB#1: Load Flow Using MATLAB Summer Semester 2020/2021
3- Power Flow Solution by Fast Decoupled Method
Maximum Power Mismatch = …….
No. of Iterations = …….
Bus V-Mag. V- Load Generator Injected
No. Angle
MW MVAR MW MVAR MVAR
Line Flow and Losses
Line Power at bus & line flow Line losses Transformer
From to MW Mvar MVA MW Mvar tap
MW Mvar
Total losses
4. Lab Report
The lab report should include as follows:
Part 1: Answer the following questions
1. What does mean by power flow analysis and what is it needed?
2. What essential data are needed to describe a bus, a generator a load
3. Why does a power system need a slack bus?
4. What is the effect of changing the Reactive power limit setting (Qmin and Qmax) on the results?
Repeat Task 2 Qmin and Qmax ( -150, 400 ) with Newton-Raphson method and observe the
results. Discus how this limit is determined in real power system.
Part 2: fully completed tables listed in Task 2 , discussion and conclusions.
Page 8 of 8
