Gauss-Seidel Method of Load Flow Analysis: Algorithm Flowchart Problems Advantages & Disadvantages
Gauss-Seidel Method of Load Flow Analysis: Algorithm Flowchart Problems Advantages & Disadvantages
Gauss-Seidel Method of Load Flow Analysis: Algorithm Flowchart Problems Advantages & Disadvantages
Analysis
Contents:
Algorithm
Flowchart
Problems
Advantages & Disadvantages
Algorithm
Step 1:
Assume a flat voltage profile 1+j0 for all buses except slack bus
Step 2:
Assume a suitable value of convergence criterion ε
Step 3:
Set iteration count k=0 and assume V10 V20 V30 ………. Vn0 except
slack bus.
Step 4:
Set bus count p=1
Step 5:
Check for slack bus. If it is slack bus then go to step-12,
otherwise go to next step.
Contd…
Step 6: Check for generator bus. If it is generator bus go to next step,
otherwise go to step 9
Step 7:
Set |Vpk|=|Vp|spec calculate the reactive power by,
Contd…
If calculated Q violates the specified limit then treat this bus as
load bus
if Qp,cal k+1 < Qp,min then Qp = Qp,min
Qp,calk+1 > Qp,max then Qp = Qp,max
go to step-9
Step 8:
For generator bus the voltage magnitude is constant. The
phase of bus voltage calculated by,
Contd…
Step 9: For the load bus the value of voltage can be calculated by,
𝑃𝑝−𝑄𝑝
𝑘+1 1 𝑝−1 𝑘+1 𝑛 𝑘
𝑉𝑝 =𝑌 ∗ − 𝑌 𝑉
𝑞=1 𝑝𝑞 𝑞 − 𝑌 𝑣
𝑞=𝑝+1 𝑝𝑞 𝑞
𝑝𝑝 𝑉𝑝𝑘
Step 10: An acceleration factor α can be used for faster convergence.
Vp,acck+1 = Vpk+α(Vpk+1- Vpk)
Then set,
Vpk+1 = Vp,acck+1 α=1.6
Step 11:
Calculate, ∆Vpk+1=Vpk+1- Vpk
Step 12:
Repeat steps 5 to 11 until all the bus voltages have been calculated .
Continue until bus count is n.
Contd…
Step 13:
Find the largest of the absolute value of
change in voltage. |∆Vmax|<ε then move to next step.
Otherwise increment the iteration count and go to
step-4.
Step-14
Calculate the line flows and slack bus power
using bus voltages.
The above algorithm has been explained in the next slide with an
example.
3 − 𝑗9 −2 + 𝑗5 −1 + 𝑗4
𝑌𝑏𝑢𝑠 = −2 + 𝑗5 5 − 𝑗14 −3 + 𝑗9
−1 + 𝑗4 −3 + 𝑗9 4 − 𝑗13
Step 2: Initial voltages are considered as 1 p.u for all buses.
𝑉20 = 1.2+j0 𝑉30 = 1+j0
Bus 1 is slack bus so its voltage remains constant for all iterations.
𝑉10 =𝑉11 = 𝑉12 = 𝑉13 = 1.05+j0
Step 3: Start iteration count, for first iteration set k=0
Step 4: Bus 1 is slack bus p=1 ∴ 𝑽𝟏
𝟏 = 1.05+j0
1
𝑄2 =(-1) im 𝑉20 ∗
∗ 𝑌21 𝑉11 + 𝑌22 𝑉20 + 𝑌23 𝑉30
= (-1) im 1.2 − 𝑗0 ∗ (−2 + 𝑗5) 1.05 + 𝑗0 + 5 − 𝑗14 1.2 + 𝑗0 + −3 +
𝑗9 (1 + 𝑗0)
=(-1)im 𝑗1.08 − 𝑗3.06 ∴ 𝑸𝟏
𝟐 = 3.06 p.u 0 < Q < 3.
The Q-limit is violated which is greater than the upper limit. Set 𝑄2 = 3 𝑝. 𝑢 Therefore the
bus acts as load bus. ∴ 𝑽𝟎
𝟐 = 1+j0
For load bus calculate the value of voltage magnitude and phase angle,
1 𝑃 𝑝 −𝑗𝑄 𝑝 𝑝 −1 𝑛
𝑉𝑝𝑘, +1 =𝑌 ∗ − 𝑞 =1 𝑌𝑝𝑞 𝑉𝑞𝑘 +1 − 𝑞 =𝑝 +1 𝑌𝑝𝑞 𝑣𝑞𝑘
𝑝𝑝 𝑉𝑝𝑘
1 𝑃 2−𝑗𝑄 2
𝑉21 = ∗ − 𝑌21 𝑉11 − 𝑌23 𝑣3
0
𝑌22 𝑉20
1 3+𝑗 3
=5−𝑗 14 − −2 + 𝑗5 (1 + 𝑗0) – −3 + 𝑗9 (1 − 𝑗0)
1−𝑗 0
1 −4+𝑗 2
= − −1 + 𝑗4 (1.05 + 𝑗0) – −3 + 𝑗9 (1.023 + 𝑗0.063)
4−𝑗 13 1−𝑗 0
Result:
The voltages at the end of first iteration is,
𝑉11 = 1.05+j0 p.u
𝑉21 =1.023 + j0.063p.u
𝑉31 = 0.803 – j0.194 p.u
References:
1. Hadi Saadat, ‘Power System Analysis’, Tata McGraw Hill Education Pvt.
Ltd., New Delhi, 21st reprint, 2010.
2. Kundur P., ‘Power System Stability and Control, Tata McGraw Hill Education
Pvt. Ltd., New Delhi, 10th reprint, 2010.
3. Pai M A, ‘Computer Techniques in Power System Analysis’, Tata Mc Graw-Hill
Publishing
Company Ltd., New Delhi, Second Edition, 2007.
4. J. Duncan Glover, Mulukutla S. Sarma, Thomas J. Overbye, ‘ Power System
Analysis & Design’, Cengage Learning, Fifth Edition, 2012.
5. Olle. I. Elgerd, ‘Electric Energy Systems Theory – An Introduction’, Tata
McGraw Hill Publishing Company Limited, New Delhi, Second Edition, 2012.
6. C.A.Gross, “Power System Analysis,” Wiley India, 2011.
7. M.Jeraldin Ahila “Power System Analysis”, Lakshmi Publications, Chennai,
Eleventh Edition 2017.
8. Other Web Resources