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