Scribd
Upload
83 views6 pages

Continuation Power Flow (CPF) : Join Github Today

This repository contains MATLAB code to calculate continuation power flow (CPF) for the IEEE 14-bus test system. The code performs CPF by using either load or voltage as the continuation parameter. It calculates the power flow using Newton-Raphson at each step, changing the parameter value and converging to find the next point on the loadability curve. The maximum loadability for each bus is around 4 times the initial load, with bus 11 reaching a slightly higher value of 4.06 times.

Uploaded by

Aditya Tiwari
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
83 views6 pages

Continuation Power Flow (CPF) : Join Github Today

This repository contains MATLAB code to calculate continuation power flow (CPF) for the IEEE 14-bus test system. The code performs CPF by using either load or voltage as the continuation parameter. It calculates the power flow using Newton-Raphson at each step, changing the parameter value and converging to find the next point on the loadability curve. The maximum loadability for each bus is around 4 times the initial load, with bus 11 reaching a slightly higher value of 4.06 times.

Uploaded by

Aditya Tiwari
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

sayonsom / continuation-power-flow

Dismiss
Join GitHub today
GitHub is home to over 40 million developers working together to host
and review code, manage projects, and build software together.

Sign up

This repository contains MATLAB Code for the calculating Continuation Power Flow (CPF) for IEEE-14 bus system.

14 commits 1 branch 0 packages 1 release 1 contributor GPL-2.0

Branch: master New pull request Find file Clone or download

sayonsom ... Latest commit 521961f on May 4, 2018

Results adding cpf repository 2 years ago

Calculate_PcalcQcalc.m adding cpf repository 2 years ago

Calculate_Ybus.m adding cpf repository 2 years ago

Jacobian_NRLF.m adding cpf repository 2 years ago

LICENSE Initial commit 2 years ago

NRLF.m adding cpf repository 2 years ago

README.md ... 2 years ago

ReadData.m adding cpf repository 2 years ago

main.m adding cpf repository 2 years ago

README.md

downloads 891/month Eclipse Marketplace DOI 10.5281/zenodo.1209882

Continuation Power Flow (CPF)


This repository contains MATLAB Code for the calculating Continuation Power Flow (CPF) for IEEE-14 bus system. The reference
book for this code was Mariesa Crow's Computational Methods for Electric Power Systems.

How to run the code?


Go to main.m , and change the bus number in line 43, for the bus for which CPF needs to be performed.

Steps to perform CPF using the developed program:


Get bus and branch data
Form Y bus matrix
Calculate scheduled values of real and reactive power

UPPER PART: LOAD AS CONTINUATION PARAMETER

/
Find out the Jacobian Matrix

Predictor Step:

Step size (sigma) = 0.1


Lambda = 0.1
Find out the voltage and angles

Corrector Step:

Perform Newton Raphson load flow by increasing scheduled real and reactive power by a factor of lambda from
predictor step.
Make sure the power flow converges within a certain number of iterations. If the power flow converges, the final value of
voltage and lambda is saved. It is used to plot the loadability curve.
If power flow converged, perform the next iteration of the predictor step (Step 5), and get the next point to plot, and
carry on until the power flow stops to converge in the predictor step. OR, If power flow did not converge, then move to
VOLTAGE AS CONTINUATION PARAMETER

VOLTAGE AS CONTINUATION PARAMETER:

Step size (sigma) := 0.005


Find the Jacobian Matrix
Predictor Step:
Lambda := last known value from Step 6.
Find out the voltage and angles
Corrector Step:
Perform Newton Raphson load flow by increasing scheduled real and reactive power by a factor of lambda from
predictor step.
Solution to V – Vpredicted = 0 is an additional criteria for convergence (Eq. 3.167 of book).
Make sure the power flow converges within a certain number of iterations. If the power flow converges, the final value of
voltage and lambda is saved. It is used to plot the loadability curve.
If power flow converged, perform the next iteration of the predictor step, and get the next point to plot, and carry on
until the power flow stops to converge in the predictor step (Step 9).
If the calculated lambda value is 75% of lambda last calculated in Step 6, CHANGE CONTINUATION PARAMETER AGAIN
to Power.

CHANGE CONTINUATION PARAMETER AGAIN

Repeat steps 5 and 6.


Plot the final loadability curve.

Table 1: Points where continuation factors need to be changed to avoid matrix singularity

Columns 2 and 3: Load as CF to Voltage as CF


Columns 4 and 5: Voltage as CF to Load as CF

Bus Number Lambda at change V(p.u.) Lambda at change V(p.u.)

4 4.0 0.7641 2.9 0.4491

5 4.0 0.7496 3.2 0.3946

7 4.0 0.8463 2.8 0.6013

9 4.0 0.7652 2.9 0.5152

10 4.0 0.7792 2.9 0.5692

11 4.0 0.9055 2.8 0.7955

14 4.0 0.7330 2.8 0.5780


/
If we increase the IEEE-14 bus load for all buses upto 4 times, power flow fails to converge. However, if we significantly decrease
the step size and change the continuation parameter the maximum loadability point is reached. For all buses, it is around 4.05,
however for bus 11, it is 4.06. The maximum loadability of all the buses are shown in Table 2.

Table 2: Maximum loadability points

Bus Number Lambda Max V(p.u.) at Lambda Max

4 4.0593 0.7041

5 4.0594 0.6796

7 4.0602 0.8013

9 4.0592 0.7102

10 4.0591 0.7292

11 4.0638 0.8805

14 4.0585 0.6880

Results

/
/
/
/

You might also like