Load Flow Analysis

3.
LOAD FLOW ANALYSIS

3.1. INTRODUCTION
Successful operation of electrical power systems requires that:
Generation must supply the demand (load) plus the losses,
Bus voltage magnitudes must remain close to rated values,
Generators must operate within specified real and reactive power limits,
Transmission lines and transformers should not be overloaded for long periods.
Therefore it is important that voltages and power flows in an electrical system can be
determined for a given set of loading and operating conditions. This is known as the load
flow problem. The study of various methods of solution to a complex power system network
is referred to as load flow study. The solution provides the voltages at various buses, power
flowing in various lines and line losses. The main information obtained from a load flow
study are the magnitude and phase angle of the voltage at each bus and the real and reactive
power flowing in each line. The load flow solution also gives the initial conditions of the
system when the transient behaviour of the system is to be studied. The load flow study of a
power system is essential to decide the best operation of existing system, for planning the
future expansion of the system and for designing a new power system.
A load flow study of a power system generally requires the following steps
i.
Representation of the system by single line diagram
ii.
Determine the impedance (admittance) diagram using the information in

single line diagram
iii.
Formulation of network equations
iv.
Solution of network equations
3.2. BUS CLASSIFICATION

In a power system the buses are meeting points of various components. The
generators will feed energy to buses and loads will draw energy from buses. In the network of
a power system the buses becomes nodes and so a voltage can be specified for each bus.
Therefore each bus in a power system, is associated with four quantities and they are real
power, reactive power, magnitude of voltage and phase angle of voltage. In a load flow
problem two quantities (out of four) are specified for each bus and the remaining two
quantities are obtained by solving the load flow equations. The buses of a power system can
be classified into following three types based on the quantities being specified for the buses
Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam
Page 1
i) Load bus or PQ bus

The bus is called load bus, when real and reactive components of power are specified
for the bus. The load flow equations can be solved to find the magnitude and phase of bus
voltage. In a load bus the voltage is allowed to vary within permissible limits, for example
5% .
ii) Generator bus or voltage controlled bus or PV bus
The bus is called generator bus, when real power and magnitude of bus voltage are
specified for the bus. The load flow equations can be solved to find the reactive power and
phase angle of bus voltage. Usually for generator buses, reactive power limits will be
specified.
iii) Slack bus or swing bus or reference bus
The bus is called slack bus if the magnitude and phase angle of bus voltage are
specified for the bus. The slack bus is the reference bus for load flow solution and usually one
of the generator buses is selected as the slack bus.
The following table gives the summary of the above classifications.
Bus type
Quantities specified
Load bus
P,Q
Generator bus
P,|V|
Q,
|V|,
P,Q
Slack bus
Quantities to be obtained
|V|,
3.3. NECESSITY OF SLACK BUS

Basically the power system has only two types of buses and they are load and generator
buses. In these buses only power injected by generators and power drawn by loads are
specified, but the power loss in transmission lines are not accounted.
In a power system the total power generated will be equal to sum of power consumed by
loads and losses
i.e
Sum of complex
Sum of complex
Total power loss in
power of generators
power of loads
transmission lines
or
Total power loss in
Sum of complex
Sum of complex
transmission lines
power of generators
power of loads
Page 2
The transmission line losses can be estimated only if the real and reactive power of all buses
is known. The power in the buses will be known only after solving the load flow equations.
For these reasons, the real and reactive power of one of the generator bus is not specified and
this bus is called slack bus. It is assumed that the slack bus generates the real and reactive
power required for transmission line losses. Hence for a slack bus, the magnitude and phase
of bus voltage are specified and real and reactive powers are obtained through the load flow
solution.
3.4. DATA FOR LOAD FLOW STUDIES
Irrespective of the method for the solution, the data required is common for any
load flow. These are presented below. All data is normally in p.u.
i) System Data: This should give information on
Number of buses n
Number of PV buses
Number of load buses
Number of loads
Slack bus number
Voltage magnitude of slack bus
Reactive power limits for the generator bus
Number of transmission lines
Number of transformers
Number of shunt elements
Base MVA
Tolerance limit
Maximum permissible number of iterations
ii) Generator bus Data: For every generator bus p the data required is
Bus number
Active power generation, PGp
Reactive power limits
Voltage magnitude, Vp,spec.
iii) Load Data: For all loads, the data required is
Bus number
Active power demand, PDp
Reactive power demand QDp
Page 3
iv) Transmission line Data: For every transmission line connected between buses p and q,
the data required is
Starting bus number, p
Ending bus number, q
Resistance and reactance of the line
Half line charging admittance
v) Transformer data: For every transformer connected between buses p and q, the data
required is
Starting bus number, p
Ending bus number, q
Resistance and reactance of the transformer
Off nominal turns ratio, a
vi) Shunt element data: The data needed for shunt element is
Bus number where element is connected
Shunt admittance (G sh+jBsh)
3.5. FORMULATION OF LOAD FLOW EQUATIONS USING YBUS MATRIX
The load flow equations can be formed using either the mesh or node basis equations of a
power system. However, from the view point of computer time and memory, the nodal
admittance formulation using the nodal voltages as the independent variables is the most
economic. As far as power system networks are concerned, the major advantages of the nodal
approach may be listed as follows:
Data preparation is easy.
The number of variables and equations is usually less than with the mesh
method for power networks.
Parallel branches do not increase the number of variables or equations.
Node voltages are available directly from the solution, and branch currents
are easily calculated.
Off-nominal transformer taps can easily be represented.
The load flow equations, using nodal admittance matrix formulation for a three bus system as
shown in fig.(3.1), are developed first and then they are generalized for n-bus system.
Page 4
Fig.(3.1): 3-bus system

By applying KCL at node 1
I1 = (V1-V2)y12+(V1-V3)y13
= V1(y12+y13) V2y23 V3y13
=> I1= V1Y11 +V2Y12 +V3Y13
Where
y12 =
1
,
z12
y 23 =
--- (3.1)
1
z 23
, y 31 =
1
z 31
Here Y11= (y12+y13) is shunt charging admittance at bus 1.

Y12=-y12 is the mutual admittance between the buses 1 and 2
Y13=-y13 is the mutual admittance between the buses 1 and 3
Similarly nodal current equations for the other nodes can be written as follows:
I2= V1Y21 +V2Y22 +V3Y23
--- (3.2)
I3= V1Y31 +V2Y32 +V3Y33
--- (3.3)
These equations can be written in matrix form as follows

I1
Y11
Y12
Y13 V1
I2
Y21 Y22
Y23 V2
I3
Y31 Y32
Y33 V3
--- (3.4)
In general the above equation can be written in matrix notation as

I =Y V
--- (3.5)
The elements Y11 , Y22, Y33 forming the diagonal terms are called self admittances. The self
admittance of a node n is equal to the sum of admittances of all the elements connected to
node n. In general the diagonal element Ypp of the bus admittance matrix is equal to the sum
of admittances of all the elements connected to bus p.
i.e. Ypp = yp1+yp2 +-----------+ypn
The elements Y12 , Y13, Y21 , Y23 , Y31 , Y32 forming the off-diagonal terms are called mutual
admittances.
Page 5
Y12 = Y21= - y12 , Y23 = Y32= - y 23 , Y13= Y31 = -y13

It is to be noted that all mutual admittance terms have a negative sign. In general, the offdiagonal term of the bus admittance matrix is equal to the negative of admittance connected
between nodes p and q
i.e Y pq
y pq
In compact form, the eq.(3.4) can be written as

3
Ip
Y pqVq , p=1,2,----n
--- (3.6)
q 1
From this we can write nodal current equation for an n bus system where each node is
connected to all other nodes.
n
Ip
Y pqVq , p=1,2,----n
--- (3.7)
q 1
Y ppVP
Y pqVq
q 1
q p
Vp
Ip
1
Y pp
Y pp
--- (3.8)
Y pqVq
q 1
q p
Ip has been substituted by the real and reactive powers, because normally in a power system
these quantities are specified.
Now, we know that
V p* I p
Ip
Pp
jQp
Pp
jQp
V
--- (3.9)
*
p
From equations (3.8) & (3.9)
Vp
1 Pp jQ p
Y pp
V p*
1
Y pp
Y pqVq , p 1,2
--- (3.10)
q 1
q p
If the power system elements have mutual coupling, the bus admittance matrix cannot be
found directly by inspection of the single line diagram. In presence of mutual coupling
between power system elements the inspection method fails. In such a case Y bus can be
formed from graph theory approach. However, the mutual coupling between power system
elements exist only in case of transmission lines running in parallel for a long distance. But
Page 6
this coupling is also weak. Therefore, for all practical purposes the mutual coupling can be
ignored and Ybus is formed by inspection method.
Properties of Ybus matrix
The nodal admittance matrix in (3.4) or (3.5) has a well-defined structure, which makes it
easy to construct. Its properties are as follows:
Square of order n n.
Symmetrical, since Ypq = Yqp
Complex.
Each off-diagonal element Ypq is the negative of the branch admittance between
nodes p and q, and is frequently of value zero.
Each diagonal element Ypp is the sum of the admittances of the branches which
terminate at node p including branches to ground.
Very few non-zero mutual admittances exist in practical networks. Therefore matrix
Y is generally highly sparse.
Problem-1: Determine the nodal admittance matrix for the power system represented by
the single line diagram as shown in the fig.
Solution:
Y11 = y12+y13 = 1+j2+3+j3 = 4+j5
Y12 = -1- j2
Y13 = -3- j3
Y22 = y21+y23 = 1+j2+2+j1 = 3+j3
Y23 = -2- j1= Y32
Y33 = y31+y32 = 3+j3+2+j1 = 5+j4
The nodal admittance matrix (
is
Page 7
Y11
YBUS
Y12
Y21 Y22
Y31 Y32
Y13
Y23 =
Y33
j5
j2
j3
1
3
j2
j3
j3
j1
j1
j4
3.6. ITERATIVE METHODS OF LOAD FLOW SOLUTION

Iterative methods can be used to solve the load flow equations which are non-linear.
The iterative methods are:
i.
Gauss- Seidal method

a. Without PV bus
b. With PV bus
ii.
Newton Rapshon method
iii.
Decoupled load flow method
iv.
Fast-Decoupled load flow method
The static load flow equations are of such complexity that it is not possible to obtain exact
analytical solution. We must use some approximate techniques that will give a sufficiently
accurate numerical solution.
The solution of the load flow problem is obtained in the following manner
1. Draw the single line diagram and write bus admittance matrix.
2. Identify the buses and branches by numbers.
3. Write the power flow equations for the given network in suitable form.
4. An initial solution is guessed for the given power system network.
5. This solution is used in conjunction with static load flow equations to compute a new
and better second estimation.
6. The second estimation is then used for finding the third estimation and so on.
7. The iterations are continued till the desired convergence is reached.
8. Calculate the desired quantities at the various buses.
3.6.1 GAUSS-SEIDAL (GS) METHOD
Case-1: Gauss-Seidal (GS) method without PV bus
Gauss-Seidal method without PV bus, is an iterative method can be chosen first
because of its inherent simplicity. We shall apply this method to solve our static load
flow equations of general n-bus system. Presently we shall consider the case, when
voltage control buses or PV buses are not present. This means that we have n-1 load
buses or PQ buses, the remaining one bus is the slack bus.
Thus in this method our unknown variables are
Vp = |Vp | e jp, p = 2,3..n
Page 8
Which are n-1 complex unknown variables V2, V3, V4 Vn and S1 = P1+jQ1
Where the complex power
at the slack bus can only be computed if the unknown
|Vp | and p at the (n-1) load buses are computed first.

The current entering at the
bus of an n-bus system is given by
Ip
Y pqVq , p=1,2,----n
q 1
Y ppVP
Y pqVq
q 1
q p
Vp
Ip
Ypp
1
Ypp
--- (3.11)
YpqVq
q 1
q p
Now, we know that for
bus
Vp*Ip= Pp- jQp
Ip
Pp
jQ p
V
---- (3.12)
*
p
By substituting Ip from eqn.(12) in eqn.(11) , we have

Vp
1 Pp jQ p
Y pp
V P*
--- (3.13)
Y pqV q
q 1
q p
For Gauss-Seidal iterative method without PV bus, we can write the above equation
as
V pK
1 Pp jQ p
Y pp
(V p* ) K
--- (3.14)
Y pqVqK
q 1
q p
Algorithm when PV buses are absent

1. Read the system data and formulate YBUS for the given power system network.
2.
Assume a flat voltage profile (1+j0) for all node voltages except the slack bus. Let
slack bus voltage be (a+j0) and it is not modified in any iteration.
3.
Assume a suitable value of convergence criterion
. If the absolute value of the
maximum change in voltage between any two consecutive iterations is less than a prespecified tolerance
, the convergence is achieved and the iterative procedure is
terminated.
4. Set iteration count k= 0
5. Set bus count p= 1
Page 9
6. Check for the slack bus. If it is a slack bus then go to step (8), since voltage at the
slack bus is fixed both in magnitude and phase, it does not vary during iterative
procedure. If it is not a slack bus then go to next step.
7. Calculate bus voltage
using equation VpK
V pk
difference in the bus voltage using
V pk
1 Pp jQp
Ypp (Vp* ) K
YpqVqK
and the
q 1
q p
V pk
8. Advance the bus count by 1 to evaluate other values of Vpk
and Vpk Check all the
buses have been taken into account or not. If yes, go to the next step, otherwise go
back to step (6).
9. Determine the largest absolute value of change in voltage
10. If
is less than the pre specified tolerance
, then evaluate line flows and
print the voltages and line flows. If not, advance the iteration count k= k+1 and go
back to step (5).
Case-2: Gauss -Seidal (GS) method including PV buses
The GS method is an iterative algorithm for solving a set of non-linear load flow
equations. The non-linear load flow equations are
given by eqn.(3.13) can be
represented for convenience as follows

Vp
1 Pp jQ p
Y pp
VP*
p 1
Y pqVq
q 1
Y pqVq
where p=1,2,3------n,
--- (3.15)
q p 1
The variables obtained from the above equation are node voltages V 1, V2, ------Vn.
In the GS method an initial value of voltages are assumed and they are denoted as
V10, V20-------, Vn0. On substituting these initial values in the above equation and by
taking p=1, the revised value of bus 1 voltage V11 is computed. The revised value of
bus voltage V11 is replaced for initial value V10 and the revised bus 2 voltage V21 is
computed. Now replace the value of V11 for V10 and V21 for V20 and again calculate
the voltage for bus 3 and so on.
The process of computing all the bus voltages as explained above is called one
iteration. The iterative procedure is repeated till the bus voltages converges within
prescribed accuracy.
Based on the above discussion the load flow eqn.(3.15) can be written in modified
form as given below
Page 10
V pK
1 Pp jQ p
Y pp (VP* ) K
Where V pk
V pk
p 1
Y pqVqK
q 1
Y pqVqK
--- (3.16)
q p 1
k th iteration value of bus voltage Vp

1
(k 1) th iteration value of bus voltage Vp
In eqn.(3.16), to compute the (k 1) th
iteration value of bus-p voltage,
the (k 1) th iteration values of voltages are used for all buses less than p and k th
iteration values off voltages are used for all buses greater than or equal to p.
The eqn.(3.16) is applicable for load bus, since in load bus, changes in both
magnitude and phase of voltages are allowed. But in generator bus the magnitude of
voltage remains constant and so the eqn.(3.16) is used to calculate the phase angle of
voltage.
In the load flow analysis one of the bus is taken as a slack bus or reference bus and so
its voltage will not change. Therefore in each iteration the slack bus voltage remains
constant and it is not modified.
For a generator bus, the reactive power is not specified. Therefore in order to calculate
the phase of bus voltage of a generator bus using eqn.(3.16), we have to estimate first
the reactive power, from the bus voltages and admittance matrix as shown below
Pp
jQ p
p 1
Y pqVq
VP*
Y pqVq
q 1
q p
p 1
Pp
jQ p
*
P
Y pqVq
q 1
Y pqVq
--- (3.17)
q p
From the above eqn.(3.17), the equation for complex power in bus-p during (k+1)th
iteration can be obtained as given below.
p 1
Ppk
jQ pk
(VPk ) *
Y pqVqk
Y pqVqk
q 1
--- (3.18)
q p
The reactive power of bus-p during (k 1) th iteration is given by imaginary part of

eqn.(18). So the reactive power of bus-p during (k 1) th iteration is given by
p 1
Q pk
( 1) I ..P.of (VPk ) *
Y pqVqk
q 1
Y pqVqk
--- (3.19)
q p
Also for a generator bus a lower and upper limits for reactive power will be specified.
In each iteration, the reactive power of generator bus is calculated using eqn.(3.19)
Page 11
and then checked with specified limits. If it violates the specified limits then the
reactive power of the bus is equated to the limit violated and it is treated as load bus.
If it does not violate the limits then the bus is treated as generator bus.
Computation of Slack bus power and Line flows
The slack bus power can be calculated after the voltages have converged. The eqn.(3.17)
can be used to calculate the slack bus power. Here, bus-p is slack bus.
n
Pp
jQ p
VP*
Y pqVqk
q 1
Consider a line connecting between buses p and q as shown in fig. (3.2).Usually the
transmission line is connected to buses using transformer at its ends. The -equivalent
model of a transmission line with transformer at its both ends is as shown in fig. (3.2).
Fig.(3.2)
From fig.(3.2)
I pq
I qp
(V p
(Vq
Vq ) y pq
V p ) y pq
V p y pq 0
Vq y pq 0
Complex power injected by bus-p in the line pq is

*
*
S pq V p I pq V p (V p Vq ) y pq V p y pq0
Complex power injected by bus-q in the line pq is
*
*
S qp Vq I pq Vq (Vq V p ) y pq Vq y pq 0
The complex power loss in the line pq is given by
S= Spq + Sqp
Algorithm when PV buses are present
1) Read the system data and formulate YBUS for the given power system network.
2) Assume a flat voltage profile (1+j0) for all the bus voltages except the slack bus. Let
3) Assume a suitable value of
called convergence criterion. Here
is a specified change
in the bus voltage that is used to compare the actual change in bus voltage between
and (k 1) th iteration.
4) Set iteration count k= 0
Page 12
5) Set bus count p=1.

6) Check for slack bus. If it is a slack bus then go to step (13), otherwise go to next step.
7) Check for generator bus. If it is a generator bus go to next step, otherwise go to step (9)
8) Replace the value of the voltage magnitude of generator bus in that iteration by the
specified value. Keep the phase angle same as in that iteration. Calculate Q for
generator bus.
The reactive power of the generator bus can be calculated by using the following
equation
p 1
1
Q pk ,cal
( 1) I .P.of (VPk ) *
Y pqVqk
q 1
Y pqVqk
q p
The calculated reactive power may be within specified limits or it may violate the
limits. If the calculated reactive power violates the specified limit for the reactive
power then treat this bus as the load bus. The magnitude of the reactive power at this
bus will correspond to the limit it has violated
1
i.e. if Q pk ,cal
1
or if Qpk ,cal
Q p,min
then Q p
Qp,max then Q p
Q p ,m in
Q p , m ax
Since the bus is treated as load bus, take actual value of V pk for (k+1)th iteration
i.e. | V pk | need not be replaced by | Vp |sep when the generator bus is treated as
load bus. Go to step (10).
9) For generator bus the magnitude of voltage does not change and so for all iterations the
magnitude of bus voltage is the specified value only. The phase of the bus voltage can
be calculated as shown below.
1
V pK,temp
k 1
p
1 Pp jQ p
Y pp
(VPk ) *
tan
p 1
Y pqVqK
q 1
Y pqVqK
q p
1
I .P. of V pk,temp
1
R.P. of V pk,temp
Now the (k+1)th iteration voltage of the generator bus is given by
V pk
| V p | spe
k 1
p
Where | Vp |spe is magnitude of specified voltage.

After calculating V pk
for generator bus go to step (12)
Page 13
10) For the load bus the (k+1)th iteration value of load bus-p voltage, V pk
can be
calculated with the following equation.
V pK
1 Pp jQ p
Y pp
(VPk ) *
11) An acceleration factor
p 1
Y pqVqK
q 1
Y pqVqK
q p 1
can be used for faster convergence. If acceleration factor is
specified then modify the (k+1)th iteration value of bus-p voltage using the following
equation.
1
V pk,acc
Then set
V pk
V pk
(V pk
V pk )
1
V pk,acc
12) Calculate the change in bus-p voltage, using the relation
V pk = V pk
V pk
Advance the bus count by 1 to evaluate other values of V pk
and
V pk
13) Check all the buses have been taken into account or not. If yes, go to the next step,
Otherwise go back to step (6).
14) Determine the largest absolute value of change in voltage | V| max
15) If | V| max is less than the pre specified tolerance , then evaluate line flows and print
the bus voltages and line flows. If not, advance the iteration count K= K+1 and go
back to step (5).
Important Note
For load bus the active and reactive powers are considered as negative, when
generation of active power(PG) and reactive power(QG) are not specified for
the given power system network. When both generation and demand of load
bus are given then the active power is P=PG-PD and reactive power is
Q=QG-QD.
For generator bus the active and reactive powers are always considered as
positive.
In a particular iteration if the calculated reactive power for the generator bus
violates the given limits, then in that iteration that bus is taken as load bus. But
the signs of active and reactive powers will remains positive, even if the bus is
treated as load bus.
Page 14
In buses having generators and loads connected to it, either the net power will
be specified or the generator and load power will be individually specified.
3.6.1.1 Acceleration factor ()
In the GS method, a large number of iterations are required to arrive at the specified
convergence. The rate of convergence can be increased by the use of acceleration
factor to the solution obtained after each iteration. The acceleration factor is a
multiplier that enhances correction between the values of voltage in two successive
iterations.
If the acceleration factor() is specified then modify the (k+1) th iteration value of the
bus-p voltage using the following equation
1
V pk,acc
V pk
Then set V pk
(V pk
1
V pk )
1
V pk,acc
The choice of a specific value of acceleration factor depends upon the system
parameters. The optimum value of
is 1.6
Problem 2: The system data for load flow solution are given in the following tables.
Determine the voltages at various buses at the end of the first iteration by using GS
method. Take =1.6
Bus specifications
Line admittances
Bus code
Admittance
Bus code
Remarks
1-2
2-j8
1.06
Slack bus
1-3
1-j4
0.5
0.2
PQ bus
2-3
0.666-j2.664
0.4
0.3
PQ bus
2-4
1-j4
0.3
0.1
PQ bus
3-4
2-j8
Solution :
From the table-1,the admittances of various branches are calculated as follows
y12= 2-j8, y13=1-j4, y23= 0.666-j2.664, y24=1-j4, y34= 2-j8
Y11= y12 + y13= 2-j8 +1-j4 = 3-j12
Y22= y12+ y23 + y24= 2- j8 +0.666- j2.664+ 1- j4=3.666- j 14.664
Y33= y31+y32+ y34= 1- j4 +0.666-j2.664 +2-j8 = 3.666-j14.664
Y44= y42+y43= 1- j4 +2- j8 =3-j12
Y12=Y21= -y12 =-2+j8
Page 15
Y13=Y31 = -y13 = -1+j4

Y14= Y41=0
Y23=Y32= -y23= -0.666+j2.664
Y24= Y42 = -y24 = -1+j4
The bus admittance matrix of the given power system is
Y11
YBUS
Y12
Y13
Y14
Y21 Y22
Y31 Y32
Y23 Y24
Y33 Y34
Y41 Y42
Y43 Y44
j12
j8
j4
3.666
j8
j14.664
0.666
1
0.666
j 2.664
j4
3.666
j4
j 2.664
j4
j14.664
j8
j8
j12
The initial values of the bus voltages are considered as 1p.u except the slack bus.
V 20 = V30 = V 40 = 1+j0
The bus 1 is slack bus and so its voltage remains at the specified value for all iterations.
i.e V10 = V11 = V1k = 1.06+ j0.0
Since the buses are PQ buses the specified real and reactive powers are considered as
load powers. Therefore negative sign is attached to the specified powers. For first
iteration k=0, the system has four buses and p will take values from 1 to 4. Here all the
buses are load buses except bus1.
The calculations of bus voltages for first iterations are shown below.
V10 = V11 = 1+j0 (slack bus)
V pK
V21
1 Pp jQ p
Y pp
(VPk ) *
1 P2
Y22 1
jQ2
j0
1
3.666 j14.664
p 1
Y pqVqK
q 1
Y pqVqK
q p 1
Y21V11 Y23V30 Y24V40

0.5 j 2
( 2 j8)1.06 ( 0.666
1 j0
( 1 j 4)(1 j 0)
j 2.664)(1
j 0)
Page 16
0.5 + j0.2 + 2.12 - j8.48 + 0.666 - j2.664 + 1 - j4

3.666 - j14.664
3.286 - j14.944
3.666 - j14.664
= (1.0119-j0.029) p.u
V21,acc
V20
(V21 V20 )
=1+1.6(1.0119-j0.029-1)
= (1.019-j0.0464) pu
Now V21
V21,acc = (1.019-j0.0464) pu
1 P3
Y33 1
V31
jQ3
j0
Y31V11 Y32V21 Y34V40
1
3.666 j14.664
- 0.4 + j0.3
(-1+ j4)1.06 ( 0.666
1 j0
( 2 j8)(1 j 0)
j 2.664)(1.019- j0.0464)
- 0.4 + j0.3 + 1.06 - j4.24 + 0.555 - j2.755 + 2 - j8

3.666 - j14.664
3.215 - j14.6855
3.666 - j14.664
= ( 0.9942-j0.0293) pu
V31,acc
V30
(V31 V30 )
= 1+1.6(0.9942-j0.0293-1)
= (0.9907-j0.0469) pu
Now V31
V41
V31,acc
1 P4
Y44 1
0.9907-j0.0469
jQ4
j0
Y41V11 Y42V21 Y43V31
1
- 0.3 + j0.1
- 0 *1.06 - (-1+ j4)(1.019- j0.0464) - (-2 + j8)(0.9907- j0.0469)
3 - j12
1 - j0
=
2.1396 - j12.0418
3 - j12
= (0.9864-j0.0683) pu
V41,acc
V40
(V41 V40 )
= 1+1.6(0.9864-j0.0683-1)
Page 17
= (0.9762-j0.1093) pu
Now V 41 = V41,acc = (0.9762-j0.1093)pu
The bus voltages at the end of first iteration are
V11 = (1+j0) pu
V 21 = (1.019-j0.0464) pu
V 31 = (0.9907-j0.0469) pu
V 41 =(0.9762-j0.1093) pu
Problem-3: Given the simultaneous equations

x1+x2=4, 2x1+x2=5
0
Using initial values x1 =2 and x20= 3, write down the values for x11 and x21 using GS method.
Solution: From the given equation, we can get
x1=4-x2, x2=5-2x1
By using GS method
x11 = 4 x20 =4 3= 1
x21 = 5 2x11=5 2x2 = 5-2 = 3
Problem-4: For the system shown in the figure, determine the voltages at the end of first
iteration by GS method. Take =1 and bus specifications are given in the following table.
Generator is connected at bus 2.
Bus specifications
Bus
Admittance data
Remarks
code
Element
Admittance
1-2
1-j5
1.06
Slack
1-3
1.2-j4
0.5
0.1Q21
1.04
PV
2-3
0.5-4j
0.4
0.3
PQ
2-4
1.1-j2
3-4
1.2-j3
0.2
0.1
PQ
Solution:
The bus admittance matrix is as given below
Page 18
2.2
j9
j5
2.6
1.2
j4
j5
1.2
j4
j11
0.5
j4
1.1
j2
j11
1.2
j3
0.5
j4
1.1
j2
2.9
1.2
j3
2.3
j5
In the given system bus-1 is slack, bus-2 is generator bus and bus-3 , bus-4 are load
buses. The initial voltages of load buses are assumed as (1+j0) pu. For slack and
generator buses the specified voltages are used as initial values.
V10= V11= .. = V1k= 1.06 (slack bus)
V20 = 1.4+j0 (generator bus) [initial phase is assumed as 0]
V30= (1+j0) pu (load bus)
V40= (1+j0) pu (load bus)
For the generator bus the specified powers are considered as positive powers but for
load buses the specified powers are considered as ve powers.
For first iteration, k=0, in each iteration the slack bus voltage need not be
recalculated. In each iteration the reactive power for generator bus as to be calculated
and checked for violation of the specified limits. If the limits are violated then it is
treated as load bus.
The calculation of bus voltage for first iteration is shown below.
V11= V10= (1.06+j0) pu
( bus 2 is slack bus)
The bus-2 is a generator bus and to calculate its reactive power Q2
Q p,cal = -I.P of
Here p=2, k=0, n=4
Q2,cal
Note: Here V2
I .Pof
0
is same as V2
spec
not same then we have to replace V2
Q2,cal
I .Pof 1.04
I .Pof 0.0458
j5 1.06
2.6
and so V2 is used for calculation as such if it is

0
with V2
j11 1.04
spec.
0.5
j4 1
1.1
j2 1
j 0.1456
0.1456 pu
Page 19
The specified range for Q2 is 0.1 Q2 1. The calculated value of Q2 is within this
range and so the reactive limit is not violated. Therefore the bus can be treated as
generator bus.
Now P2
0.5, Q2
0.1456 , V2
1.04
j0
Since the bus-2 is treated as generator bus, then V2

by the phase of V2
V p ,temp
V2
temp
temp
1 P2 jQ2
0 *
Y22
v2
p 1
Y pqVq k
2.6
3.1408
2 .6
Y pqVq
q 1
spec
and phase of V2 is given
Y21V1
q p 1
Y23V3
0.5 j 0.1456
1.04 j 0
j11
1.1 j 2 1
V2
1 Pp jQ p
k *
Y pp
Vp
k 1
Y24V4
j5 1.06
j0
0.5
j4 1
j0
j0
j11 .44
j11
11 .8633
74 .65 0
11 .3031
76 .7 0
1.0496 2.05 0
V2
1
temp
V2
V2
1.0496 2.05 0 pu
V2
2.05 0
temp
1
2
spec
1.04 2.05 pu 1.0393
j 0.0372 pu
The bus-3 and bus-4 are load buses. The voltages of load bus are calculated using the
following equation
Vp
k 1
V3
1 Pp jQ p
k *
Y pp
Vp
1 P3 jQ3
0 *
Y33
V3
1
2.9
j11
p 1
Y pqVq
k 1
Y pqVq
q 1
q p 1
Y31V1
0.4 j 0.3
1 j0
1.2 j3 1
Y32V2
1.2
Y34V4
j 4 1.06
0.5
j 4 1.0489
j 0.0375
j0
Page 20
2.7405 j11 .0786

2.9 j11
11 .4125
11 .3759
1.0032
76 .11 0
75 .23 0
0.880
1.0031 j 0.01545 pu
1
V4
1 P4 jQ4
1
1
1
Y41V1 Y42V2 Y43V3
*
0
Y44
V4
1
2.3
j5
0.2 j 0.1
0 1.06
1.1 j 2 1.0393 j 0.0372
1 j0
1.2 j 3 1.0031 j 0.01545
2.1751 j 4.9655
2.3 j5
5.4210
5.5036
66 .34 0
65 .30
0.9850
1.040 pu
0.9848
j 0.0178 pu
The bus voltages at the end of the first iteration are

1
1.06
1.0393
1.0031 j0.01545 1.0032
0.9848
V1
V2
V3
V4
j 0 1.06 0 0 pu
j 0.0371
1.04 2.05 0 pu
j 0.0179
0.9850
0.880 pu
1.04 0 pu.
3.6.1.2 Advantages and disadvantages of GS method

Advantages of GS method
i) Calculations are simple and so programming task is lesser.
ii) The memory requirement is less.
iii) Use full for small size systems.
Disadvantages of GS method
i) Requires large number of iterations to reach convergence.
ii) Not suitable for large systems.
iii) Convergence time increases with the size of the system.
Page 21
3.6.2 NEWTON-RAPSHON METHOD

3.6.2.1 Development of load flow equations
The NR method of load flow analysis is an iterative method which
approximates the set of non-linear simultaneous load flow equations to a set of linear
simultaneous equations using Taylors series expansion and the terms are limited to first
order approximation. The equations for NR method are derived as follows
Case 1: In rectangular form

We know that for an n-bus system
Pp jQp = Vp*Ip = Vp*
Y pqVq
q 1
Let Vp= ep+jfp,

Vq = eq+jfq
Ypq = Gpq-jBpq
Where ep &fp are real and imaginary parts of V p respectively.
eq & fq are real and imaginary parts of V q respectively.
Gpq & Bpq are conductance and susceptance of admittance Ypq
n
Pp-jQp = Vp*
Y pqV q
q 1
n
[G pq
= (ep-jfp)
jB pq )( eq
jf q )
q 1
n
= (ep-jfp)
[(G pq eq
f q B pq )
j ( f q G pq
f q B pq )
f p ( f q G pq
eq B pq )
eq B pq )]
q 1
n
e p (eq G pq
Pp-jQp =
q 1
f p (eq G pq f q B pq ) e p ( f q G pq eq B pq )
j
q 1
By comparing real and imaginary parts on both sides

n
[e p (eq G pq
Pp=
f q B pq )
f p ( f q G pq
eq B pq )]
--- (3.20)
f q B pq ) e p ( f q G pq
eq B pq )]
--- (3.21)
q 1
[ f p (eq G pq
Qp=
q 1
Vp
e 2p
f p2
--- (3.22)
Page 22
The above three set of equations are the load flow equations and it can be seen that
they are non-linear equations in terms of the real and imaginary components of nodal
voltages.
The voltages of a slack bus will be a known quantity in a power system and so it need
not be solved. For load buses Pp and Qp will be specified and we have to solve V p.
For a generator bus Pp and V p will be specified and we have to solve Q p and phase
angle of Vp i.e. p.
Case 2: In polar form
We can also formulate the load flow problem using NR method in polar
coordinates , say for any buses p & q we have
Vp
Vp e
Vq
Vq e
Vp
Vp e
, Ypq
Ypq e
pq
pq
For any bus p we have

Pp
jQ p
Vp
Y pqVq
q 1
Vp e
n
p
Y pq e
Vq e
j(
V pVq Y pq e
q 1
pq
q)
q 1
By comparing real and imaginary parts on both sides

n
Pp =
V p Y pqVq cos(
pq
q 1
n
= V p Y pp cos
V p Y pqVq cos(
pp
pq
--- (3.23)
q 1
q p
Qp =
V p Y pqVq sin(
pq
q 1
n
= V p Y pp sin
V p Y pqVq sin(
pp
pq
--- (3.24)
q 1
q p
The above two equations are the load flow equations in polar form.
3.6.2.2 Mathematical background for N-R method
Let (x1,x2,------------xn) be a set of unknown variables and (y1,y2,-----------yn) be set of
specified quantities. Now the specified quantities can be expressed as a non-linear
function of unknown variables as shown below.
Page 23
y1 = f1(x1, x2, --------------xn)

y2 = f2(x1, x2, --------------xn)
:
--- (3.25)
yn = fn(x1, x2, --------------xn)

Let us assume an approximate initial solution x10 , x 20 ,
x n0 for the above
equations. The prefix zero refers to zeroth iteration in the processing of solving the
above non-linear equations.
x n0 are the corrections required for x10 , x20
x10 , x 20 ,
Let
0
---------- xn
respectively for the next better solution.

Now the non-linear equations can be expressed as shown below, i.e. they can be
expressed as functions of modified variables x10
x10 , x 20
y1 = f1( x10
x10 , x 20
x 20 ,
x n0
x n0 )
y2 = f2( x10
x10 , x 20
x 20 ,
x n0
x n0 )
yn = fn ( x10
x 20 ,
x n0
x n0
--- (3.26)
x10 , x 20
x 20 ,
x n0
x n0 )
The above equations are linearized about the initial guess using Tailors expansion.
The linearized equations with second order and higher order derivatives neglected are given
below.
0
-------xn0)
-------xn0)+
y1 = f1(x1 , x2 ,
y2 = f2(x1 , x2 ,
f2
x1
0
1
f1
x1
0
1
x
0
0
2
f1
x2
0
2
f2
x2
0
2
fn
x2
0
n
f1
xn
0
n
f2
xn
0
n
fn
xn
x
0
---(3.27)
yn = fn(x1 , x2 ,
--------xn0)+
fn
x1
0
1
Note: Tailors series expansion for any function f(x) is given by

f( x
x ) = f (x )
f
x
( x0 )2
2!
x2
Page 24
The above equations can be written as

y1 f1( x , x ,
0
1
0
1
0
2
0
n
x )=
0
2
0
n
y2 - f2( x , x ,
x )=
f1
x1
0
1
x
0
f2
x1
0
1
0
1
f2
x2
0
2
fn
x2
0
n
0
0
n
f2
xn
0
n
fn
xn
f1
xn
---(3.28)
0
2
0
n
yn - fn( x , x ,
Let
0
2
f1
x2
0
2
x )=
fn
x1
0
1
y1- f1( x10 , x 20 ,
x n0 ) =
y1
y2- f2( x10 , x 20 ,
x n0 ) =
y2
x n0 ) =
yn
yn- f1( x10 , x 20 ,
Now the above equations can be written as

y1 =
y2 =
0
1
f1
x1
0
1
f2
x1
0
2
f1
x2
0
2
f2
x2
yn =
0
n
f1
xn
0
n
f2
xn
0
n
fn
xn
--- (3.29)
0
1
fn
x1
fn
x2
0
2
The above equations can be written in matrix form
y1
y2
=
y3
f1
x1
f2
x1
fn
x1
f1
x2
f2
x2
fn
x2
. .
0
f1
xn
f2
xn
fn
xn
x10
x 20
x n0
--- (3.30)
Page 25
B =JC
Where
B=
--- (3.31)
f1
x1
y1
y2
y3
; J=
f2
x1
fn
x1
f1
x2
f2
x2
fn
x2
. .
0
f1
xn
f2
xn
fn
xn
; C=
x10
x 20
x n0
Here J is the first derivative matrix and it is called Jacobian matrix. The elements of
Jacobian matrix are obtained by evaluating the first derivatives at the assumed solution.
The B matrix is called residual column vector. The elements of B are the difference
between the specified quantities and calculated quantities at the assumed solution. With
the elements of J and B are known, the elements of matrix C are obtained by
solving matrix eqn. (3.31)
The solution of matrix eqn.(3.31) gives
x10 , x 20 ,
x n0 . The next better
solution is obtained as follows.

x11
x10
x10
x 12
x 20
x 20
--- (3.32)
x 1n
x n0
x n0
With the new solution given by eqn. (3.32) the process is repeated to find next solution.
The iterative process is terminated if any one of the following condition is satisfied.
(i) The largest (magnitude of the) element in the B matrix is less than a pre-specified value.
(ii)The largest (magnitude of the ) element in the C matrix is less than pre-specified value.
3.6.2.3 Applying NR method to load flow problem
Consider a power system with n-buses. The bus-1 is usually selected as slack bus. The other
buses (i.e bus-2 to bus-n) can be either generator bus or load bus. The specified quantities for
load buses are Pp and Qp and for generator buses are Pp and V p .
Case (i) : When the power system has all the (n-1) buses are load buses.
In this case, bus-1 is slack bus and bus-2 to bus-n are load buses. Let P2,P3-------- Pn be
the specified real powers and Q2,Q3, ---------Qn be the specified reactive powers of (n-1) load
Page 26
buses. The unknown variables are real part of voltages e2,e3, --------en and imaginary part of
bus voltages f2,f3, --------fn. Now the matrix equation B=JC for this power system problem
will be in the form shown below.
P2
e2
Pn
e2
Q2
e2
Qn
e2
P2
Pn
Q2
Qn
P2
e3
Pn
e3
Q2
e3
Qn
e3
P2
en
Pn
en
Q2
en
Qn
en
P2
f2
Pn
f2
Q2
f2
Qn
f2
P2
f3
Pn
f3
Q2
f3
P2
fn
Pn
fn
P2
fn
Qn
fn
e2
en
--(3.33)
f2
fn
2(n-1)x2(n-1)
P
Q
J1
J3
J2
J4
--- (3.34)
Case (ii) : When the power system has both load and generator buses
In this case also bus-1 is slack bus and buses 2 to m be load buses and buses (m+1) to
n are generator buses. Let P2, P3,----------Pn be the specified real power of (n-1) buses. Let
Q2, Q3, -------Qm be the specified reactive powers of load buses. Let Vm
, Vm
,-------- Vn
be the specified magnitude of voltages of generator buses. The unknown variables are real
part of bus voltages e2,e3,----en and imaginary part of bus voltages f2,f3, --------fn . Now the
matrix equation B = JC for this case will be in the following form
P2
Pn
Q2
Qm
Vm
Vn
2
1
P2
e2
Pn
e2
Q2
e2
Qm
e2
2
Vm 1
P2
e3
Pn
e3
Q2
e3
Qm
e3
2
Vm 1
e2
2
Vn
e3
2
Vn
e2
e3
P2
en
Pn
en
Q2
en
Qm
en
2
Vm 1
P2
f2
Pn
f2
Q2
f2
Qm
f2
2
Vm 1
P2
f3
Pn
f3
Q2
f3
Qm
f3
2
Vm 1
en
2
Vn
f2
2
Vn
f3
2
Vn
en
f2
f3
P2
fn
Pn
fn
Q2
fn
Qm
fn
2
Vm 1
fn
2
Vn
fn
e2
en
f2
fm
fm 1
fn
--- (3.35 )
Page 27
P
Q
V
J1
J2
= J3
J5
J4
e
f
J6
--- (3.36)
Computing the elements of Jacobian matrix

The elements of Jacobian matrix (J) can be derived from the load flow equations as given
below.
Case 1: NR method in Rectangular form
The load flow equations can be written in rectangular form as given below
n
e p eq G pq
PP =
f q B pq
f p f q G pq
eq B pq
q 1
= e p e p G pp
f p B pp
e p eq G pq
f p f p G pp e p B pp +
f q B pq
f p f q G pq eq B pq
q 1
q p
n
f p eq G pq
QP =
f q B pq
e p f q G pq
eq B pq
q 1
= f p e p G pp
f p B pp
e p f p G pp
f p eq G pq
e p B pp +
f q B pq
e p f q G pq eq B pq
q 1
q p
VP
ep
fp
J1 : Off-diagonal elements are
Pp
eq
q 1
e p G pq
= e p G pq
f p B pq , q
f p B pq
---(3.37)
Diagonal elements are
Pp
eP
= 2e p G pp
f p B pp
f p B pp
eq G pq
f q B pq
q 1
q p
n
= 2e p G pp
eq G pq
---(3.38)
f q B pq
q 1
q p
Pp
fq
q 1
e p B pq
= e p B pq
f p G pq , q
f p G pq
---(3.39)
Page 28

Pp
= 2 f p G pp
fp
e p B pp
e p B pp
f q G pq
eq B pq
q 1
q p
n
= 2 f p G pp
f q G pq
eq B pq
---(3.40)
q 1
q p

Qp
eq
= e p B pq
f p G pq , q
---(3.41)

Qp
2e p B pp
ep
f q G pq
e q B pq
---(3.42)
q 1
q p

Qp
e p G pq
fq
f p B pq , q
--- (3.43)

Qp
= 2 f p B pp
fp
eq G pq
f q B pq
--- (3.44)
q 1
q p

2
Vp
=0, q
eq
--- (3.45)

Vp
ep
= 2ep
--- (3.46)

Vp
=0, q
fq
---(3.47)
diagonal elements are

Vp
fp
= 2fp
--- (3.48)
Page 29
The elements of the matrices are obtained by partially differentiating the load flow
equations w.r.t. a unknown variable and then evaluating the first derivatives using
the solution of previous iteration. For first iteration the initial assumed values are ep0
and fp0 for p = 2,3, -----,n
The elements of the residual column matrix B is the difference between the
specified value of the quantity and the calculated value of the quantity using the
solution of previous iteration.
Let Pp,spec, Qp,spec. and V p
solution the value of
spec
be the specified quantities at the bus-P. for the initial

can be calculated using the load flow equations.
Now for the first iteration

= Pp,spec = Qp.spec-
Vp
Vp
2
spec
Vp
0 2
After calculating the elements of Jacobian matrix J and residual column vector B
the elements of increment voltage vector C can be calculated by using any standard
technique.
Now the next better solution will be
These values of voltages will be used in the next iteration.

The process will be repeated and in general the new better estimation for bus voltages
will be
The process is repeated till the magnitude of the largest element in the residual column
vector B is less than a pre-specified value.
Algorithm for NR method in rectangular form
2) Assume a flat voltage profile (1+j0) for all nodal voltages except the slack bus. Let
Page 30
called convergence criterion i.e. if the largest of the
absolute value of the residues exceeds
, the process is repeated, otherwise
terminated.
6) Check for slack bus. If it is a slack bus go to step (11), otherwise go to next step.
7) Calculate the real and reactive power of bus-P using the following equations
n
Ppk
e kp eqk G pq
f qk B pq
f pk f qk G pq
eqk B pq
f pk eqk G pq
f qk B pq
e kp f qk G pq
eqk B pq
q 1
n
Q pk
q 1
8) Calculate the change in active power
Ppk
Pp,spec
Ppk
9) Check for generator bus. If yes, compare the calculated reactive power, Q pk with the
limits. The calculated reactive power may be within specified limits or it may violate
given the limits. If the calculated reactive power violates the specified limit, then fix
the reactive power generation to the corresponding violated limit and treat this bus as
load bus and go to the next step.
if Qpk
i.e
or if Q pk
Qp,min
then Q p , spec
Q p,max then Q p , spec
Q p ,m in
Q p ,m ax
If the reactive power limit is not violated then evaluate the voltage residue
VPK
VP , spec
VPK
And go to step (11).

10) Calculate the change in reactive power for load bus (or for the generator bus treated as
load bus)
Change in reactive power QPK
QP,Spec QPK
11) Advance the bus count by 1 i.e p=p+1 and check if all the buses have been taken into
account or not. If yes, go to the next step, Otherwise go back to step (6).
12) Determine the largest of the absolute value of the residue.
13) If the largest of the absolute value of the residue is less than
, go to step (18).
14) Evaluate the elements for Jacobian matrix.

Page 31
15) Calculate voltage increments
by using matrix inverse technique.
16) The new bus voltages can be calculated as follows

; p=1,2,----,n. except slack bus
; p=1,2,----,n. except slack bus
K 1
p
tan
VPK
f PK 1
e pK 1
VPK
K 1
P
17) Advance iteration count i.e k=k+1 and go back to step (5).
18) Evaluate the line flows and slack bus power.
Case 2: NR method in Polar form
The load flow equations can be written in polar form as given below
n
Pp =
V p Y pqVq cos(
pq
q 1
n
= V p Y pp cos
V p Y pqVq cos(
pp
pq
--- (3.49)
q 1
q p
Qp =
V p Y pqVq sin(
pq
q 1
n
= V p Y pp sin
V p Y pqVq sin(
pp
pq
--- (3.50)
q 1
q p
Now the linear equation in the polar form becomes
P
Q
J1
J3
J2
J4
--- (3.51)
Where J1,J2,J3,J4 are the elements of Jacobian matrix, which can be calculated in the
following manner
J1 : off-diagonal elements are
Pp
V pY pqVq sin(
pq
), q
--- (3.52)
Page 32
Pp
V p Y pqVq cos(
pq
--- (3.53)
--- (3.54)
q 1
q p

Pp
Vq
= V p Y pq cos(
pq
), q

Pp
Vp
= 2 V p Y pp cos
pp
VqY pq cos(
pq
--- (3.55)
q 1
q p

Qp
= - V p Y pqVq cos(
pq
), q
--- (3.56)

Qp
p
V p Y pqVq cos(
pq
--- (3.57)
q 1
q p

Qp
Vq
= V p Y pq sin(
pq
),
--- (3.58)

Pp
Vp
= 2 V pY pp sinpp +
V p Y pq sin(
pq
--- (3.59)
q 1
q p
Algorithm for NR method in polar form

called convergence criterion i.e. if the largest of the
absolute value of the residues exceeds
, the process is repeated, otherwise
terminated.
6) Check for slack bus. If it is a slack bus go to step (11), otherwise go to next step.
7) Calculate the real and reactive power of bus-p using the following equations
Page 33
Pp
V p Vq Y pq cos(
k
pq
q 1
Qp
V p Vq Y pq sin(
k
pq
q 1
8) Calculate the change in active power
Ppk
Ppk
Pp,spec
9) Check for generator bus. If yes, compare the calculated reactive power, Q p with the
given limits. The calculated reactive power may be within the specified limits or it
may violate given the limits. If the calculated reactive power violates the specified
limit, then fix the reactive power generation to the corresponding violated limit and
treat this bus as load bus and go to the next step.
if Qpk
i.e
Qp,min
or if Q pk
then Q p , spec
Q p ,m in
Q p,max then Q p , spec
Q p ,m ax
If the reactive power limit is not violated treat this bus as generator bus.
10) Calculate the change in reactive power for load bus or generator bus (or for the
generator bus treated as load bus)
Change in reactive power Qpk
Qp,Spec Qpk
11) Advance the bus count by 1 i.e p = p+1 and check if all the buses have been taken into
account or not. If yes, go to the next step, Otherwise go back to step (6).
12) Determine the largest of the absolute value of the residue.
13) If the largest of the absolute value of the residue is less than
, go to step (18).
14) Evaluate the elements for Jacobian matrix.

15) Calculate phase angle and voltage increments
and
V by using matrix inverse
technique.
k 1
Vp
p
k 1
V pk
Vp
V pk
p
k
Vp
k 1
p
Page 34
Symmetry property in Jacobian matrix

Using the rectangular coordinates (see equations3.37, 3.39, 3.41, 3.43,) a careful
examination for the off-diagonal elements of the sub-matrices [J1], [J2], [J3] and [J4]
would reveal certain interesting properties.
And
[J4]pq= -[J1]pq
---(3.60)
[J2]pq = [J3]pq
---(3.61)
This property reduces the computational efforts considerably as it is enough only to

compute the off-diagonal elements of any two sub-matrices.
It may be noted that we do not see the symmetry property in the Jacobian, if polar
coordinates are used. However, if we replace
|V |
in the eqn.(3.51), we
|V |
|v| by
have.
J1 J 2
J3J4
---(3.62)
|V |
|V |
Now the expressions for off-diagonal terms

Pp
J1
V pY pqVq sin(
pq
pq
), q
--- (3.63)
Qp
J4
pq
J2
J3
Vq = V pVqY pq sin(
Vq
Pp
pq
pq
Vq = V pVqY pq cos(
Vq
Qp
= - V p Y pqVq cos(
pq
pq
pq
), q
), q
), q
---(3.64)
---(3.65)
--- (3.66)
From the above equations (3.63), (3.64), (3.65) and (3.66), it is seen
J1
pq
J2
pq
J4
--- (3.67)
pq
J3
pq
---(3.68)
Thus with slight modification in equation (3.51), we get the symmetry property in the
Jacobian, which is observed in the case of expressing Jacobian in rectangular coordinates.
The elements of Jacobian (J) are calculated with the latest voltage estimate and calculated
power. However, the procedure (i.e algorithm) here, is the same as that of the rectangular
Page 35
coordinates. The formulation in the polar coordinates takes less computational effort and also
requires less memory space.
3.6.2.4. Advantages and Disadvantage of NR method
Advantages:
1. The NR method is faster, more reliable and the results are accurate.
2. Requires less number of iterations for convergence.
3. The number of iterations are independent of the size of the system (i.e. no. of buses)
4. Suitable for large size systems
Disadvantages:
1. The programming logic is more complex than GS method.
2. The memory requirement is more.
3. No. of calculations per iteration are higher than GS method
3.6.2.5 COMPARISON OF GS AND NR METHOD
1. For GS method the variables are expressed in rectangular coordinates where as in NR
method, they are expressed in polar coordinates. If rectangular coordinates are used for
NR method then memory requirement will be more.
2. The no. of mathematical operations per iteration will be lesser in GS method than NR
method. Hence computation time per iteration is less in GS method.
3. The GS method has linear convergence characteristics where as the NR method has
quadratic convergence characteristics. Hence NR method converges faster than GS
method.
4. In GS method no. of iteration increases with no. of buses but in NR method the no. of
iterations remains constant and it does not depend on the size of the system.
5.In GS method convergence is affected by the choice of slack bus and the presence of
series capacitors but the NR method is less sensitive to these factors.
6. The NR method needs only 3 to 5 iterations to reach an acceptable solution for a large
system. But GS method requires large no. of iteration (30 or more) for same level of
accuracy.
Problem-4: In the load flow data for the sample power system are given below. The voltage
magnitude of bus 2 is to be maintained at 1.04 p.u. The maximum and minimum reactive
power limits of the generator at bus 2 are 0.35 and 0.0 p.u respectively. Determine the set of
load flow equations at the end of first iteration by using NR method.
Impedance for sample system
Page 36
Bus code
Impedance
Line changing admittance
1-2
0.08+j0.24
0.0
1-3
0.02+j0.06
0.0
2-3
0.06+j0.18
0.0
Schedule of generation and loads
Solution:
y12 =
Bus
Assumed
Generation
code
voltages
MW
1.06+j0.0
1.0+j0.0
0.2
0.0
0.0
0.0
1.0+j0.0
0.6
0.25
1
z12
1
0.08
j 0.24
Load
MVAR
MW
MVAR
= 1.25-j3.75
y13 = 5-j15, y23 = 1.667-j5

6.25
j18.75
1.25
YBUS =
1.25
j 3.75
j 3.75 2.916
j8.75
1.666
j15
1.666
j5
j15
6.666
j5
j 20
From the nodal admittance matrix and assumed voltage solution

G11 = 6.25
B11 = 18.75
e1 = 1.06,
f1 = 0.0
G12 = -1.25 = G21
B12 = -3.75 = B21
e2 = 1.0,
f2 = 0.0
G13 = -5 = G31
B13 = -15 = B31
e3= 1.0,
f3 = 0.0
G22 = -2.916
B22 = 8.75
G23 = -1.666 = G32
B23 = -5
G33 = 6.66
B33 = 20
PP =
e p eq G pq
f q B pq
f p f q G pq
eq B pq
e2 eq G 2 q
f q B2 q
f 2 f q G2 q
eq B2 q
q 1
P2 =
q 1
= e2 (e1G21 + f1B21) + f2 (f1G21 e1B21) + e2 (e2G22 + f2B22)

+ f2 (f2G22 e2B22) + e2 (e3G23+ f3B23) + f2 (f3G23 e3B23)
= 1(1.06x-1.25 + 0x-3.5) + 0 1.6x-3.75 + 1(1x2.916) + 0+ 1x1x-1.666 + 0
= -1.325 + 2.916 + - 1.666
= -0.075
Page 37
Similarly P3 = -0.3
n
f p eq G pq
QP =
f q B pq
e p f q G pq
eq B pq
q 1
Q2 = -0.225
Q3 = -0.9
i.e P20 = 0.2 -(-0.075) = 0.275
P2 = P2 specified P2 calculated
P3 = -0.6 (-0.3) = -0.3
Since the lower limit on Q2 is 0.0 and the value of Q2 as calculated above violates this limit,
bus-2 is treated as a load bus
Q2, spec = 0.0
Q2 = 0.0 - (-0.225) = 0.225
Q3 = - 0.25 - (-0.9) = 0.65
Diagonal elements:
Pp
ep
= 2epGPP +
P2
= 2e2G22 +
e2
n
q 1
q P
3
q 1
q 2
(eqGpq + fqBpq)
(eqG2q + fqB2q)
= 2e2G22 + e1G21 + f1B21 + e3G23 + f3B23

= 2(1) (2.916) +1.06(-1.25) +0(-3.75) +1(-1.666) +0(-5)
= 2.848
P3
= 2e3G33 + e1G31 + f1B31 + e2G32 + f2B32
e3
= 2(1) (6.666) +1.06(-5) +0(-15) +1(-1.666) +0(-5)

= 6.3666
Pp
fp
= 2fpGPP +
P2
= 2f2G22 +
f2
n
q 1
q P
3
q 1
q 2
(fqGpq -eqBpq)
(fqG2q -eqB2q)
=2f2G22 + f1G21 e1B21+f3G23 - e3B23

= 2(0) (2.916)-1.06(-3.75)-1(-5)
= 8.975
P3
f3
=20.9
Page 38
Off-diagonal elements:
Pp
eq
= epGpq - fpBpq
P2
= e2G23 - f2B23 = - 1.666
e3
P3
= e3G32 - f3B32 = - 1.666
e2
Pp
fq
= epBpq + fpGpq
P2
= e2B23 + f2G23 = - 5.0
f3
P3
= e3B32 + f3G32 = - 5.0
f2
Similarly we find out the derivatives of the reactive power

Diagonal elements:
Qp
ep
Q2
e2
= 2epBPP -
n
q 1
q P
(fqGpq -eqBpq)
= 2e2B22 - f1G21 + e1B21 - f3G23 + e3B23

= 2(2)(8.75)+1.06(-3.75)+1(-5) = 8.525
Q3
e3
Qp
fp
19 .1
= 2fpBPP
Q2
f2
n
q 1
q P
(eq Gpq + fqBpq)
2.991
Q3
f3
6.966
Off-diagonal elements:
Qp
eq
= epBpq+ fpGpq
Q2
= e2B23+ f3G23 = 1(-5) + 0 = -5
e3
Q3
= e3B32+ f2G32 = 1(-5) + 0 = -5
e2
Page 39
Qp
= - epGpq + fpBpq
fq
Q2
= - e2G23 + f2B23 = -1(-1.666) +0 =1.666
f3
Q3
= - e3G32 + f3B32 = -1(-1.666) +0
f2
The set of linear equations at the end of first iteration are
0.275
2.846
0.3
1.666
0.225
0.65
3.525
5
8.975
e2
6.366
5.0
20.9
e3
5
19.1
2.991
1.666
1.666
6.966
f2
1.666
f3
Problem-5: In case the reactive power constraints at bus-2 in the previous problem is
-0.3 Q2
0.3 Determine the equations at the end of first iteration.
Solution:
Since Q2 = 0.225and the lower limit is -0.3, therefore the bus-2 behaves like a generator bus
P2 = 0.2 (-0.075) = 0.275
P3 = -0.6 - (-0.3) = -0.3
Q3 = -0.25 (-0.9) = 0.65
Since bus-2 behave like a generator bus therefore
|V2|2 = |V2|2 - |V2cal|2
= 1.042 -1.02 = 0.0816
The Jacobian elements corresponding to rows P2, P3, Q3 remains same as in previous
problem, those of Q2 will be change and they are calculated as follows
| V2 | 2
= 2e2 = 2
e2
| V2 | 2
= 2f2 = 0
f2
| V2 | 2
=0
e3
| V2 | 2
=0
f3
The set of equations will be as given below
Page 40
0.275
2.846
e2
20.9
e3
6.966
0.0
f2
1.666 8.975
0.3
1.666
6.366
5.0
0.65
0.0816
5.0
2.0
19.1
0.0
1.666
0.0
f3
3.6.3 DECOUPLED LOAD FLOW (DLF) METHOD

In NR method, the inverse of the Jacobian has to be computed at every iteration. When
solving large interconnected power system, alternative solution methods are possible, taking
into account certain observations made of practical systems. These are
. The real power changes (P) are less sensitive to changes in voltage magnitude and
are mainly sensitive to change in bus voltage angles. In other words the coupling
between active power P and the bus voltage magnitude |V| is relatively weak
The reactive power changes (Q) are less sensitive to change in bus voltage angles
and are mainly sensitive to change in voltage magnitude. In other words the coupling
between reactive power (Q) and bus voltage phase angle () is also weak.
Thus the weak coupling is utilized in the development of decoupled load flow method in
which P is decoupled from V and Q is decoupled from .
With these assumptions the equation
P
Q
J1
J3
J2
J4
(3.69)
|V |
is reduced to the following form
J1
J4
Therefore
(3.70)
|V |
P = [J1] [ ]
(3.71)
Q = [J4] [ |V|]
(3.72)
The load flow equations of NR method in polar form are

Pp = |Vp2 Ypp| cos
| V pVqY pq | cos(
pp
pq
q 1
q p
Qp = |Vp2 Ypp| sin
| V pVqY pq | sin(
pp
pq
q 1
q p
J1: Off-diagonal elements are

Pp
V pVqY pq sin(
pq
),
--- (3.73)
Page 41
Pp
p
V pVqY pq sin(
pq
),
--- (3.74)
),
q 1
q p
J4: Off-diagonal elements are
Qp
Vq
V p Y pq sin(
pq
--- (3.75)
Qp
Vp
2 V pY pp sin
VqY pq sin(
pp
pq
--- (3.76)
),
q 1
q p
Equations (3.71) and (3.72) can be constructed and solved simultaneous with each
other at each iteration, updating the
and
in each iteration using the
equations (3.73), (3.74), (3.75) and (3.76).

A better approach is to conduct each iteration by first solving equation (3.71) for ,
and use the updated in constructing and then solving equation (3.72) for
. This
will result in faster convergence then in the simultaneous mode.

The main advantage of the decoupled load flow (DLF) method as compared to the NR
method is its reduced memory requirements in storing the Jacobian elements.
However, the time required per iteration of the DLF method is practically the same as
that of NR method. In DLF method more no. of iterations are required for
convergence because of the approximations made in it.
3.6.4 FAST DECOUPLED LOAD FLOW (FLDF) METHOD
This FDLF method is an extension of NR method formulated in polar coordinates
with certain approximations which results a faster algorithm for load flow solution.
In this method both the speed as well as the sparsity are exploited (make good use)
The load flow equations for NR method in polar form can be written as
--- (3.77)
--- (3.78)
These equations after linearization can be written in matrix form as
--- (3.79)
Where H, N, M and L are the elements (viz J1, J2, J3, J4 ) of the Jacobian matrix.
Page 42
The first assumption under decoupled load flow method is that the real power changes
are less sensitive to change in voltage magnitude and mainly sensitive to
changes in phase angle. Similarly, the reactive power changes
are less sensitive
to changes in phase angle but mainly sensitive to change in voltage magnitude. With
these assumptions, the equation (3.79) is reduced to
--- (3.80)
Equation (3.80) is decoupled equation which can be expressed as
--- (3.81)
--- (3.82)
Using equations (3.77) and (3.78) the elements of the Jacobian matrices H and L are
obtained as follows:
Off-diagonal elements of H:
=
=
=
--- (3.83)
Off-diagonal elements of L:
=
--- (3.84)
From equations (3.83) and (3.84)
The diagonal elements of H:
(From eqn. (3.78))

--- (3.85)
Page 43
The diagonal elements of L:
--- (3.86)
In the case of fast decoupled load flow methods following approximations are further
made for evaluating Jacobian elements
i)
Bpq>>Gpq(since, the X/R ratio of transmission lines is high in well designed

system
ii)
The voltage angle difference
between two buses in the system is
very small. This means that
and sin
=0
iii)
With these assumption the Jacobian elements now become
With these Jacobian elements, equations (3.81) and (3.82) becomes

--- (3.87)
--- (3.88)
Where
are the elements of
matrix.
Further decoupling and logical simplification of the FLDF algorithm is achieved by

1. Omitting effect of phase shifting transformers.
2. Setting off-nominal turns ratio of transformers to 1.0
3. In forming
, omitting the effect of shunt reactors and capacitors which
mainly effect reactive power.

4. Ignoring series resistance of lines in forming the Ybus.
5. Dividing each of the equations (3.87) and (3.88) by
and setting
p.u, we get
Page 44
P
V
B1
Q
V
B11
Here both
--- (3.89)
--- (3.90)
and
are real and sparse and have structures of H and L
respectively. Since, they contain only network admittances, they are constant and
need to be triangularised only once at the beginning of the iteration. This algorithm,
which results in a very fast solution of and V, is known as fast decoupled load
flow formulation of load flow studies.
Algorithm for FDLF method
called convergence criterion.
4) Set iteration count K= 0.

6) Check for slack bus. If it is a slack bus goes to step (11), other wise go to next step.
7) Calculate active power and reactive power by using the following formula
n
Pp = PP
V p Y pqVq cos(
pq
q 1
QP
V pY pqVq sin(
pq
q 1
8) Calculate the mismatches (i.e changes) in active power Pk and reactive powers Qk.
If the mismatches are within the desirable tolerance then stop the iteration process.
9) Normalize the mismatches by dividing each entry by its respective bus voltage
magnitude.
Page 45
P2k
V2k
k
p
Q2k
V2k
P3k
V3k
Q3K
V3k
Pnk
Vnk
Qnk
Vnk
Vk and
10) Solve the following equations for the correction factors
k by using the
constant matrices BI and BII which are extracted from the bus admittance matrix.
P
V
B1
Q
V
B11

pk+1 = pk +
pk
Vpk+1= Vpk + Vpk

VPK
VPK
K 1
P
12) Check if all the buses are taken into account or not. If yes go to next step otherwise
increase bus count by 1 i.e set p=p+1 and go back to step (6).
3.7. COMPARISON OF DIFFERENT LOAD FLOW METHODS
S.No
GS Method
Rectangular
NR Method
coordinates Polar
are preferred for solution

2
coordinates
FDLF Method
are Polar coordinates are
preferred for solution
preferred for solution
More no.of iteration(i.e 30 Less no.of iteration(i.e 3 Less no.of iteration(i.e

or more) are required to to 5) are required to get 2 to 5) are required to
get the acceptable solution
the acceptable solution
get
the
acceptable
solution
Page 46
The computation time per The computation time The computation time
iteration will be less due per iteration is more i.e 8 per iteration is more i.e
to
less
no.
of times than GS method
2/3
mathematical operation
times
than
GS
method and 5 times

than NR method
Acceleration factor is used No such factor is used
No such factor is used
to get fast convergence

5
The number of iterations The number of iterations The
number
increases as the size of the independent of the size of iterations

system increases.
the system.
of
independent
of the size of the

system
Less computer memory is More computer memory Memory requirement is

required
is required
intermediate of GS and
NR method
High computation cost
Less computation cost
Suitable for small size of Suitable for large size of Suitable for large size
systems
systems
of systems
Convergence is effected Convergence

on the selection of slack independent
10
Less computation cost
is Convergence
on
the independent
is
on
the
bus
selection of slack bus
selection of slack bus
Convergence is uncertain
Convergence is certain
Convergence is certain
3.8. DC LOAD FLOW METHOD

In certain power system studies (e.g. reliability studies) a very large no. of load flow
runs may be needed. Therefore, a very fast (and not necessarily accurate, due to the
linear approximation involved) method can be used for such studies.
The method of calculating the real power flows by solving first for the bus angles is
known as the dc load flow method, in contrast with the exact non 0linear solution,
which is known as the ac solution.
Assume that bus p is connected to bus q over an impedance of Z pq. Therefore, the
active power flow can be expressed as
Ppq
V pVq
Z pq
sin(
---- (3.91)
Page 47
Where V p
Vp
Vq
Vq
The following simplifying approximations are made

Since
Now, the active load flow eqn. (3.91) can be expressed as

p
Ppq
X pq
B pq (
---- (3.92)
In matrix form
---- (3.93)
---- (3.94)
---- (3.95)
Where [B] matrix is an (n-1) (n-1) matrix dimensionally for an n-bus system. The
diagonal and off-diagonal elements of the [B] matrix can be found by adding the series
susceptances of the branches connected to bus and by setting then equal to negated series
susceptance of branch pq, respectively.
The linear equation (3.93) can be solved for by using matrix techniques.
It is possible with the dc load flow method to carry out the thousands of load flow
runs that are required for comprehensive contingency analysis on large scale systems.
In summary, the choice of a load flow method is a matter of choice between speed
and accuracy. For a given degree of accuracy, the speed depends on the size,
complexity, and configuration of the power system and on the numerical approach
chosen.
Problem-1: The load flow data for the system shown in figure is given below in the
following tables
Bus-Code
Impedance (Zpq)
1-2
j0.05pu
1-3
j0.1pu
2-3
j0.05pu
Table (1)
Page 48
Assumed
Generation
Bus Voltage
MW
MVar
MW
MVar
1.03 + j0 pu
1.0 + j0 pu
50
20
10
1.0 + j0 pu
20
20
Bus Code
Load
Table(2)
The Voltage magnitude at bus-2 to be held at 1.0p.u.The maximum and minimum reactive
power limits at bus-2 are 50 and -10 MVars respectively .With bus-1 as the slack bus, use GS
method and Ybus matrix to obtain a load flow solution up to one iteration
[JNTU, Regular, Nov - 2004]
Solution: From the table (1)
y12 =
1
1
=
=-j20
Z 12 j 0.05
y13 =
1
1
=
=-j10
Z 13 j 0.1
y23 =
1
1
=
=-j20
Z 23 j 0.05
Y11
Y12
Y13
Ybus = Y21 Y22

Y31 Y32
Y23
Y33
Y11 = y12 + y13 = -j20-j10= -j30

Y22 = y12 + y23 = -j20-j20= -j40
Y33 = y31 + y32 = -j10-j20= -j30
Y12 = Y21 = -y12 = j20
Y23 = Y32 = -y23 = j20
Y31 = Y13 = -y13= j10
j 30
Nodal admittance Matrix
Ybus =
j 20
j10
j 20
j 40
j 20
j10
j 20
j 30
The data in table-2 is to be converted into per unit value
Page 49
Per unit value =
actual value
base value
Let the base value =30 MVA

i.e
so the data in table-2 is changed accordingly
PG2 = 50 MW =
50
=1.0 p.u
50
PD2 = 20 MW =
20
=0.4 p.u
50
QD2 = 10 MVar =
PD3 = 20 MW =
10
=0.2 p.u
50
20
=0.4 p.u
50
QD3 = 20 MVar =
20
=0.4 p.u
50
Assume flat voltage profile for all the buses except slack bus i.e
V20
V30
1 p.u ,
1 p.u ,
|V2| spec =1.0 p.u
Since bus-1 is slack bus, its voltage remains constant at the specified value for all the
iterations.
V10
V11
V12
.......... .....
V1k
(1.03
j 0.0) p.u
For generator bus ,the reactive power limits are

Q2min=
10
= -0.2 p.u
50
Q2max=
50
= 1.0 p.u
50
The bus-2 is a generator bus and so calculate its reactive power, Q2

p 1
1
Q pk ,cal
( 1) I .P.of (VPk ) *
Y pqVqk
q 1
Y pqVqk
q p
Here p = 2 , k=0 ,n=3
Q21,cal
( 1) I .P.of (V20 )* Y21V11 Y22V20
Y23V30
= ( 1) I .P.of (1) ( j 20 )(1.03) ( j 40 )(1) ( j 20 )(1)

= -0.6 p.u
The specified range for Q2 is -0.2<Q<1.0. The calculated value of Q2 violates the lower
limit of the specified range for Q2 . Therefore the reactive power generation for bus-2 is
fixed at -0.2 (lower limit) and the bus-2 is treated as load bus for 1st iteration. Now
Page 50
V20
j 0 , similar to other load buses for first iteration .But P and Q are considered
as positive for bus-2 and P and Q are negative for other load buses.
V pK
1 Pp jQ p
Y pp
(VPk ) *
V21
1 P2 jQ2
Y22 (V20 )
p 1
1
q 1
Y pqVqK
q p 1
Y21V11 Y23V30
1
0.6 j 0.2
j 40
1
(0.6
Y pqVqK
j 40 .4)
j 40
( j 20 )(1.03) ( j 20 )(1)
40 .4 90 .9
40 90
1.01 0.9 p.u
Voltage at bus-3
V31
1 P3 jQ3
Y33 (V30 )
1
j 30
Y31V11 Y32V21
0.4
( 0.083
j 0.4
( j10 )(1.03) ( j 20 )(1.01 0.9)
j 30 .1)
j 30
30 .1 89 .842
30 90
1.003
0.158 p.u
The voltages at the end of first iteration are

V11
1.03 0 p.u
V31
1.003
V21
1.01 0.9 p.u
0.158 p.u
Problem-2: The load flow data for the power system as shown in the figure is given in the
following tables.
Bus-Code
Impedance (Zpq)
1-2
0.08+j0.24pu
1-3
0.02+j0.06pu
2-3
0.06+j0.18pu
Table (1)
Bus
Assumed
Generation
Load
Code
Bus Voltage
MW
MVar
MW
MVar
0
1.05 + j0 pu
20
20
1.0 + j0 pu
50
25
Page 51
1.0 + j0 pu
60
Table(2)
The voltage magnitude at bus-2 is to be maintained at 1.03 p.u .The maximum and minimum
reactive power limits of the generator at bus-2 are 35 and 0 Mvars respectively. With bus1 as
slack bus, obtain voltage at bus-3 using GS method after first iteration (assume base Mva =
50)
(JNTU , supplementary Feb-2007)
Solution: From the table (1)

y12 =
1
1
=
=1.25-j3.75
Z 12 0.08 j 0.24
y13 =
1
1
=
=5-j15
Z 13 0.02 j 0.06
y23 =
1
1
=
=1.667-j5
Z 23 0.06 j 0.18
Nodal admittance matrix

Y11
Y12
Y13
Ybus = Y21 Y22

Y31 Y32
Y23
Y33
6.25
Nodal admittance Matrix
Ybus =
j18.75
1.25
5
1.25
j 3.75
j 3.75 2.916
j8.75
1.666
j15
1.666
j5
6.666
j15
j5
j 20
The data in table-2 is to be converted into per unit value

Per unit value =
actualvalue
basevalue
Let the base value =50 MVA. So the data in table-2 is changed accordingly
i.e
PG2 = 20 MW =
20
=0.4 p.u
50
QG2 = 0, QG3 = 0, PG3 = 0

PD2 = 50 MW =
50
=1.0p.u
50
QD2 = 20 MVar =
20
=0.4 p.u
50
Page 52
PD3 = 60 MW =
60
=1.2 p.u
50
QD3 = 25 MVar =
25
=0.2 p.u
50
Assume flat voltage profile for all the buses except slack bus i.e
V20
V30
1.03 p.u ,
1 p.u ,
V2
spec
1.03 p.u |
Since bus-1 is slack bus, its voltage remains constant at the specified value for all the
iterations.
V10
V11
V12
.......... .....
V1k
(1.05
j 0.0) p.u
For generator bus ,the reactive power limits are

Q2min=
0
= 0 p.u
50
Q2max=
35
= 0.7 p.u
50
The bus-2 is a generator bus and so calculate its reactive power, Q 2

p 1
1
Q pk ,cal
( 1) I .P.of (VPk ) *
Y pqVqk
Y pqVqk
q 1
q p
Here p = 2 , k=0 ,n=3
Q21,cal
( 1) I .P.of (V20 )* Y21V11 Y22V20 Y23V30
= ( 1) I .P.of (1.03) ( 1.25
j3.75)(1.05) (2.916
j8.75)(1.03) ( 1.666
j5)(1)
= I.P of (0.2544 j0.07725)

=0.07725 pu
The specified range for Q2 is 0<Q2<0.7p.u. The calculated value of Q2 is lies with in the
given reactive power limits and so this bus can be treated as generator bus.
Now P2= -0.6, Q2 = 0.07725, V20
bus, the V21
1
V pK,temp
V21,temp
V2
spec
1.03
j 0 , since the bus-2 is treated as generator
1
1
and phase of V 2 is given by the phase of V 2 temp.
1 Pp jQ p
Y pp
(VPk ) *
1 P2 jQ2
Y22 (V20 ) *
p 1
Y pqVqK
q 1
Y pqVqK
q p 1
Y21V11 Y23V30
Page 53
1
2.916
0.6 j 0.077
1.03
j8.75
2.38
2.916
j8.75
9.29
9.20
75.17
71.56
V21
1.03
( 1.25
j 3.75 )(1.05 ) ( 1.666
j 5)(1)
j8.99
1.009
3.61
3.61 p.u
Voltage at bus-3
1 P3 jQ3
Y33 (V30 ) *
V31
1
6.666 j 20
1.2
j 0.5
1
Y31V11 Y32V21
( 5.02
(5.451 j 20 .543 )
6.666 j 20
j15 .06 )(1.05 ) ( 1.66
21 .25
21 .12
75 .13
71 .56
j 5)(1.03
1.006
3.61)
3.57 p.u
Problem-3: The data for 2-bus system is given below SG1 = unknown , SD1 = unknown
V1=1.00 p.u S1 = To be determined , SG2 = 0.25 + j QG2 p.u ,SD2 = 1+j0.5p.u The two
buses are connected by a transmission line p.u reactance of 0.5 p.u . Find Q 2 and
Neglect shunt susceptance of the tie line .Assume V2
GS Method.
V2 .
1.0 . Perform Two iterations Using
[JNTU ,Supplementary , Feb-2008]
sol) From the given data, the single line diagram can be drawn as
Z12 = j0.5
y12=
1
Z 12
j2
P2 =0.25-1 = -0.75
, V1 = 1.0 p.u
Let us assume the bus-1 as slack bus , so its voltage remains constant throughout all the
iterations
V10
V20
V11
V12
.......... .....
V1k
(1
j 0.0) p.u
1.0 p.u
The bus admittance matrix
Ybus
Y11
Y12
Y21 Y22
j2
j2
j2
j2
Page 54
p 1
1
Q pk ,cal
( 1) I .P.of (VPk ) *
Y pqVqk
Y pqVqk
q 1
q p
Here p = 2 , k=0 ,n=2
Q21,cal
( 1) I .P.of (V20 )* Y21V11 Y22V20

= ( 1) I .P.of ( j 2
= ( 1) I .P.of (1) ( j 2)(1) ( j 2)(1)
j 2) =0
Q2=0 p.u
V pK
1 Pp jQ p
Y pp
(VPk ) *
1 P2 jQ2
Y22 (V20 ) *
V21
1
j2
0.75
1
p 1
Y pqVqK
q 1
Y pqVqK
q p 1
Y21V11
j0
j2
0.75 j 2
j2
j 0.375
V21
1.068
1
2
20 .56
20 .56 p.u
V21
V2
1
2
V22
1 P2 jQ2
Y22 (V21 ) *
1
0.75 j 0
j 2 0.936 j 0.35
1
2.674
j2
j 0.2628
20.56 p.u
0.936
j0.35
Y21V12
j2 1
j2
=0.868-j1.337 = 1.594-57
2
2
V22
57
V2
2
2
57 p.u
Page 55
Page 56

Load Flow Analysis

Uploaded by

Copyright:

Available Formats

Load Flow Analysis

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Load Flow Analysis

Uploaded by

Copyright:

Available Formats

3.

LOAD FLOW ANALYSIS

Representation of the system by single line diagram

Determine the impedance (admittance) diagram using the information in

Formulation of network equations

Solution of network equations

3.2. BUS CLASSIFICATION

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

i) Load bus or PQ bus

3.3. NECESSITY OF SLACK BUS

Total power loss in

Total power loss in

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

Fig.(3.1): 3-bus system

Here Y11= (y12+y13) is shunt charging admittance at bus 1.

I3= V1Y31 +V2Y32 +V3Y33

These equations can be written in matrix form as follows

In general the above equation can be written in matrix notation as

Y12 = Y21= - y12 , Y23 = Y32= - y 23 , Y13= Y31 = -y13

In compact form, the eq.(3.4) can be written as

From equations (3.8) & (3.9)

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

3.6. ITERATIVE METHODS OF LOAD FLOW SOLUTION

Gauss- Seidal method

Newton Rapshon method

Decoupled load flow method

Fast-Decoupled load flow method

at the slack bus can only be computed if the unknown

|Vp | and p at the (n-1) load buses are computed first.

bus of an n-bus system is given by

Now, we know that for

Vp*Ip= Pp- jQp

By substituting Ip from eqn.(12) in eqn.(11) , we have

Algorithm when PV buses are absent

Assume a suitable value of convergence criterion

. If the absolute value of the

, the convergence is achieved and the iterative procedure is

using equation VpK

difference in the bus voltage using

8. Advance the bus count by 1 to evaluate other values of Vpk

and Vpk Check all the

is less than the pre specified tolerance

, then evaluate line flows and

given by eqn.(3.13) can be

represented for convenience as follows

Prepared by P.Ramana, Associate Professor, EEE, GMRIT-Rajam

k th iteration value of bus voltage Vp

(k 1) th iteration value of bus voltage Vp

In eqn.(3.16), to compute the (k 1) th

iteration value of bus-p voltage,

The reactive power of bus-p during (k 1) th iteration is given by imaginary part of

Complex power injected by bus-p in the line pq is

called convergence criterion. Here

5) Set bus count p=1.

Now the (k+1)th iteration voltage of the generator bus is given by

Where | Vp |spe is magnitude of specified voltage.