www.allsyllabus.

com

PART - A
UNIT - 1
NETWORK TOPOLOGY: Introduction, Elementary graph theory oriented graph, tree,
co-tree, basic cut-sets, basic loops; Incidence matrices Element-node, Bus incidence,
Tree-branch path, Basic cut-set, Augmented cut-set, Basic loop and Augmented loop;
Primitive network impedance form and admittance form.
6 Hours

yl
la
bu
s.

co

UNIT - 2
NETWORK MATRICES: Introduction, Formation of YBUS by method of inspection
(including transformer off-nominal tap setting), by method of singular transformation
(YBUS = AT yA); Formation of Bus Impedance Matrix by step by step building algorithm
(without mutual coupling elements).
6 Hours

PART - B

w
w

w
.a

lls

UNIT - 3 & 4
LOAD FLOW STUDIES: Introduction, Power flow equations, Classification of buses,
Operating constraints, Data for load flow; Gauss-Seidal Method Algorithm and flow
chart for PQ and PV buses (numerical problem for one iteration only), Acceleration of
convergence; Newton Raphson Method Algorithm and flow chart for NR method in
polar coordinates (numerical problem for one iteration only); Algorithm for Fast
Decoupled load flow method; Comparison of Load Flow Methods.
14 Hours

UNIT - 5 & 6
ECONOMIC OPERATION OF POWER SYSTEM: Introduction, Performance curves,
Economic generation scheduling neglecting losses and generator limits, Economic
generation scheduling including generator limits and neglecting losses; Iterative
techniques; Economic Dispatch including transmission losses approximate penalty
factor, iterative technique for solution of economic dispatch with losses; Derivation of
transmission loss formula; Optimal scheduling for Hydrothermal plants problem
formulation, solution procedure and algorithm.
12 Hours
UNIT - 7 & 8
TRANSIENT STABILITY STUDIES: Numerical solution of Swing Equation Pointby-point method, Modified Eulers method, Runge-Kutta method, Milnes predictor
corrector method. Representation of power system for transient stability studies load
representation, network performance equations. Solution techniques with flow charts.
14 Hours
Page 1
vtu.allsyllabus.com

www.allsyllabus.com

TEXT BOOKS:
1. Computer Methods in Power System Analysis- Stag, G. W., and EI-Abiad, A.
H.- McGraw Hill International Student Edition. 1968
2. Computer Techniques in Power System Analysis- Pai, M. A- TMH, 2nd
edition, 2006.

w
w

w
.a

lls

yl
la
bu
s.

co

REFERENCE BOOKS:
1. Modern Power System Analysis- Nagrath, I. J., and Kothari, D. P., -TMH,
2003.
2. Advanced Power System Analysis and Dynamics- Singh, L. P.,
New Age International (P) Ltd, New Delhi, 2001.
3. Computer
Aided
Power
System
Operations
and
Analysis- Dhar, R. N- TMH, New Delhi, 1984.
4. Power System Analysis- Haadi Sadat, -TMH, 2nd , 12th reprint, 2007

Page 2
vtu.allsyllabus.com

www.allsyllabus.com

TABLE OF CONTENTS
TOPIC

Sl. No.

1.

PAGE
NO.

UNIT-1
NETWORK TOPOLOGY-1.1 Introduction,

3-18

1.1 Elementary linear graph theory: important terms

1.2 Examples on Basics of LG Theory
1.3 1.3 Incidence matrices

6
7

1.4 Bus incidence matrix: A

NETWORK MATRICES
2.1 Rule of Inspection

yl
la
bu
s.

co

1.5 Primitive networks

UNIT-2

2.2 Bus impedance matrix

lls

2.3 Singular transformations

w
.a

2.4 Node elimination by matrix algebra

2.5 ZBUS building

13
19-51
19
21
22
26
30

w
w

2.6 Modification of Zbus for network changes

UNIT-3 &4
LOAD FLOW STUDIES
3.1 Review of numerical solution of equations
3.2 Load flow studies

41
52-81
52
53

3.3 Classification of buses for LFA

3.4 The load flow problem
3.4 The load flow problem
3.5 Algorithm for GS method
3.6 Examples on GS load flow analysis:
3.7 NR method

53
55
58
59
62
74
Page 3

vtu.allsyllabus.com

www.allsyllabus.com

3.8 Algorithm for NR method

UNIT-5 & 6
ECONOMIC OPEARATION OF POWER SYSTEM
5.1 Introduction

88-113

88

5.2 performance curves

89

5.3 economic generation scheduling neglecting losses and

generator limits

yl
la
bu
s.

UNIT-7 & 8

co

5.6 derivation of transmission loss formula

5.4 Economic schedule including limits on generator

(neglecting losses)
5.5 economic dispatch including transmission losses

TRANSIENT STABILITY STUDIES

7.1 Introduction

lls

7.2 Mechanics of rotatory motion

w
.a

7.3 Powerangle equation

7.4 Power angle equation for a salient pole machine:
7.5 Transient stability:

91

96
99
106

114-149

114
115
124
125
125

w
w

78

7.6 Modified Eulers method:

7.7 Runge-Kutta method

127
136

Page 4
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-1
NETWORK TOPOLOGY
2. INTRODUCTION

yl
la
bu
s.

co

The solution of a given linear network problem requires the formation of a set of equations
describing the response of the network. The mathematical model so derived, must describe
the characteristics of the individual network components, as well as the relationship which
governs the interconnection of the individual components. In the bus frame of reference
the variables are the node voltages and node currents. The independent variables in any
reference frame can be either currents or voltages. Correspondingly, the coefficient matrix
relating the dependent variables and the
independent variables will be either an impedance or admittance matrix. The formulation
of the appropriate relationships between the independent and dependent variables is an
integral part of a digital computer program for the solution of power system problems. The
formulation of the network equations in different frames of reference requires the
knowledge of graph theory. Elementary graph theory concepts are presented here,
followed by development of network equations in the bus frame of reference.
1.1 ELEMENTARY LINEAR GRAPH THEORY: IMPORTANT TERMS

w
w

w
.a

lls

The geometrical interconnection of the various branches of a network is called the

topology of the network. The connection of the network topology, shown by replacing all
its elements by lines is called a graph. A linear graph consists of a set of objects called
nodes and another set called elements such that each element is identified with an ordered
pair of nodes. An element is defined as any line segment of the graph irrespective of the
characteristics of the components involved. A graph in which a direction is assigned to
each element is called an oriented graph or a directed graph. It is to be noted that the
directions of currents in various elements are arbitrarily assigned and the network
equations are derived, consistent with the assigned directions. Elements are indicated by
numbers and the nodes by encircled numbers. The ground node is taken as the reference
node. In electric networks the convention is to use associated directions for the voltage
drops. This means the voltage drop in a branch is taken to be in the direction of the current
through the branch. Hence, we need not mark the voltage polarities in the oriented graph.
Connected Graph : This is a graph where at least one path (disregarding orientation)
exists between any two nodes of the graph. A representative power system and its oriented
graph are as shown in Fig 1, with:
e = number of elements = 6
n = number of nodes = 4
b = number of branches = n-1 = 3
l = number of links = e-b = 3
Tree = T(1,2,3) and
Co-tree = T(4,5,6)
Sub-graph : sG is a sub-graph of G if the following conditions are satisfied:
sG is itself a graph
Every node of sG is also a node of G
Every branch of sG is a branch of G
Page 5
vtu.allsyllabus.com

www.allsyllabus.com

For eg., sG(1,2,3), sG(1,4,6), sG(2), sG(4,5,6), sG(3,4),.. are all valid sub-graphs of
the oriented graph of Fig.1c.
Loop : A sub-graph L of a graph G is a loop if
L is a connected sub-graph of G
Precisely two and not more/less than two branches are incident on each node in L

yl
la
bu
s.

co

In Fig 1c, the set{1,2,4} forms a loop, while the set{1,2,3,4,5} is not a valid, although
the set(1,3,4,5) is a valid loop. The KVL (Kirchhoffs Voltage Law) for the loop is
stated as follows: In any lumped network, the algebraic sum of the branch voltages
around any of the loops is zero.

w
w

w
.a

lls

Fig 1a. Single line diagram of a power system

Fig 1b. Reactance diagram

Page 6
vtu.allsyllabus.com

www.allsyllabus.com

Fig 1c. Oriented Graph

lls

yl
la
bu
s.

co

Cutset : It is a set of branches of a connected graph G which satisfies the following

conditions :
The removal of all branches of the cutset causes the remaining graph to have two
separate unconnected sub-graphs.
The removal of all but one of the branches of the set, leaves the remaining graph
connected.
Referring to Fig 1c, the set {3,5,6} constitutes a cutset since removal of them isolates node
3 from rest of the network, thus dividing the graph into two unconnected subgraphs.
However, the set(2,4,6) is not a valid cutset! The KCL (Kirchhoffs Current Law) for the
cutset is stated as follows: In any lumped network, the algebraic sum of all the branch
currents traversing through the given cutset branches is zero.

w
w

w
.a

Tree: It is a connected sub-graph containing all the nodes of the graph G, but
without any closed paths (loops). There is one and only one path between every pair of
nodes in a tree. The elements of the tree are called twigs or branches. In a graph with n
nodes,
The number of branches: b = n-1
(1)
For the graph of Fig 1c, some of the possible trees could be T(1,2,3), T(1,4,6), T(2,4,5),
T(2,5,6), etc.
Co-Tree : The set of branches of the original graph G, not included in the tree is called the
co-tree. The co-tree could be connected or non-connected, closed or open. The branches of
the co-tree are called links. By convention, the tree elements are shown as solid lines while
the co-tree elements are shown by dotted lines as shown in Fig.1c for tree T(1,2,3). With e
as the total number of elements,
The number of links: l = e b = e n + 1
(2)
For the graph of Fig 1c, the co-tree graphs corresponding to the various tree graphs are as
shown in the table below:

Page 7
vtu.allsyllabus.com

www.allsyllabus.com

Basic loops: When a link is added to a tree it forms a closed path or a loop. Addition of
each subsequent link forms the corresponding loop. A loop containing only one link and
remaining branches is called a basic loop or a fundamental loop. These loops are defined
for a particular tree. Since each link is associated with a basic loop, the number of basic
loops is equal to the number of links.

Basic cut-sets: Cut-sets which contain only one branch and remaining links are called
basic cutsets or fundamental cut-sets. The basic cut-sets are defined for a particular tree.
Since each branch is associated with a basic cut-set, the number of basic cut-sets is equal
to the number of branches.

co

1.2 Examples on Basics of LG Theory:

w
.a

lls

yl
la
bu
s.

Example-1: Obtain the oriented graph for the system shown in Fig. E1. Select any four
possible trees. For a selected tree show the basic loops and basic cut-sets.

w
w

Fig. E1a. Single line diagram of Example System

Fig. E1b. Oriented Graph of Fig. E1a.

Page 8
vtu.allsyllabus.com

www.allsyllabus.com

yl
la
bu
s.

co

For the system given, the oriented graph is as shown in figure E1b. some of the valid Tree
graphs could be T(1,2,3,4), T(3,4,8,9), T(1,2,5,6), T(4,5,6,7), etc. The basic cutsets
(A,B,C,D) and basic loops (E,F,G,H,I) corresponding to the oriented graph of Fig.E1a and
tree, T(1,2,3,4) are as shown in Figure E1c and Fig.E1d respectively.

w
w

w
.a

lls

Fig. E1c. Basic Cutsets of Fig. E1a.

Fig. E1d. Basic Loops of Fig. E1a.

2.3 INCIDENCE MATRICES

Elementnode incidence matrix: A

Page 9
vtu.allsyllabus.com

www.allsyllabus.com

w
.a

lls

yl
la
bu
s.

co

The incidence of branches to nodes in a connected graph is given by the element-node

incidence matrix,A .
An element aij of A is defined as under:
aij = 1 if the branch-i is incident to and oriented away from the node-j.
= -1 if the branch-i is incident to and oriented towards the node-j.
= 0 if the branch-i is not at all incident on the node-j.
Thus the dimension of A is e n, where e is the number of elements and n is the number
of nodes in the network. For example, consider again the sample system with its oriented
graph as in fig. 1c. the corresponding element-node incidence matrix, is obtained as under:

w
w

It is to be noted that the first column and first row are not part of the actual matrix and
they only indicate the element number node number respectively as shown. Further, the
sum of every row is found to be equal to zero always. Hence, the rank of the matrix is less
than n. Thus in general, the matrix A satisfies the identity:

1.4 Bus incidence matrix: A

By selecting any one of the nodes of the connected graph as the reference node, the
corresponding column is deleted from A to obtain the bus incidence matrix, A. The
dimensions of A are e (n-1) and the rank is n-1. In the above example, selecting node-0 as
reference node, the matrix A is obtained by deleting the column corresponding to node-0,
as under:

Page 10
vtu.allsyllabus.com

co

www.allsyllabus.com

yl
la
bu
s.

It may be observed that for a selected tree, say, T(1,2,3), the bus incidence matrix can be
so arranged that the branch elements occupy the top portion of the A-matrix followed by
the link elements. Then, the matrix-A can be partitioned into two sub matrices Ab and Al
as shown, where,
(i) Ab is of dimension (bxb) corresponding to the branches and
(ii) Al is of dimension (lxb) corresponding to links.

w
w

w
.a

lls

A is a rectangular matrix, hence it is singular. Ab is a non-singular square matrix of

dimension-b. Since A gives the incidence of various elements on the nodes with their
direction of incidence, the KCL for the nodes can be written as
AT i = 0
(4)
where AT is the transpose of matrix A and i is the vector of branch currents. Similarly for
the branch voltages we can write,
v = A bus E
(5)

Examples on Bus Incidence Matrix:

Example-2: For the sample network-oriented graph shown in Fig. E2, by selecting a tree,
T(1,2,3,4), obtain the incidence matrices A and A . Also show the partitioned form of the
matrix-A.

Page 11
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Fig. E2. Sample Network-Oriented Graph

Corresponding to the Tree, T(1,2,3,4), matrix-A can be partitioned into two submatrices
as under:

Page 12
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

lls

Example-3: For the sample-system shown in Fig. E3, obtain an oriented graph. By
selecting a tree, T(1,2,3,4), obtain the incidence matrices A andA . Also show the
partitioned form of the matrix-A.

Consider the oriented graph of the given system as shown in figure E3b, below.

Page 13
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

Fig. E3b. Oriented Graph of system of Fig-E3a.

w
w

w
.a

lls

Corresponding to the oriented graph above and a Tree, T(1,2,3,4), the incidence matrices
and A can be obtained as follows:

Corresponding to the Tree, T(1,2,3,4), matrix-A can be partitioned into two submatrices
as under:

Page 14
vtu.allsyllabus.com

www.allsyllabus.com

yl
la
bu
s.

co

2.4 PRIMITIVE NETWORKS

So far, the matrices of the interconnected network have been defined. These matrices
contain complete information about the network connectivity, the orientation of current,
the loops and cutsets. However, these matrices contain no information on the nature of the
elements which form the interconnected network. The complete behaviour of the network
can be obtained from the knowledge of the behaviour of the individual elements which
make the network, along with the incidence matrices. An element in an electrical network
is completely characterized by the relationship between the current through the element
and the voltage across it.

w
w

w
.a

lls

General representation of a network element: In general, a network element may

contain active or passive components. Figure 2 represents the alternative impedance and
admittance forms of representation of a general network component.

Fig.2 Representation of a primitive network element (a) Impedance form (b)

Admittance form
The network performance can be represented by using either the impedance or the
admittance form of representation. With respect to the element, p-q, let,
vpq = voltage across the element p-q,
epq = source voltage in series with the element p-q,
ipq= current through the element p-q,
Page 15
vtu.allsyllabus.com

www.allsyllabus.com

jpq= source current in shunt with the element p-q,

zpq= self impedance of the element p-q and
ypq= self admittance of the element p-q.
Performance equation: Each element p-q has two variables, Vpq and ipq. The
performance of the given element p-q can be expressed by the performance equations as
under:
vpq + epq = zpqipq (in its impedance form)
ipq + jpq = ypqvpq (in its admittance form)
(6)

yl
la
bu
s.

co

Thus the parallel source current jpq in admittance form can be related to the series source
voltage, epq in impedance form as per the identity:
jpq = - ypq epq
(7)
A set of non-connected elements of a given system is defined as a primitive Network and
an element in it is a fundamental element that is not connected to any other element. In the
equations above, if the variables and parameters are replaced by the corresponding vectors
and matrices, referring to the complete set of elements present in a given system, then, we
get the performance equations of the primitive network in the form as under:
v + e = [z] i
i + j = [y] v
(8)

w
w

w
.a

lls

Primitive network matrices:

A diagonal element in the matrices, [z] or [y] is the self impedance zpq-pq or self
admittance, ypq-pq. An off-diagonal element is the mutual impedance, zpq-rs or mutual
admittance, ypq-rs, the value present as a mutual coupling between the elements p-q and rs. The primitive network admittance matrix, [y] can be obtained also by inverting the
primitive impedance matrix, [z]. Further, if there are no mutually coupled elements in the
given system, then both the matrices, [z] and [y] are diagonal. In such cases, the self
impedances are just equal to the reciprocal of the corresponding values of self
admittances, and vice-versa.
Examples on Primitive Networks:
Example-4: Given that the self impedances of the elements of a network referred by the
bus incidence matrix given below are equal to: Z1=Z2=0.2, Z3=0.25, Z4=Z5=0.1 and
Z6=0.4 units, draw the corresponding oriented graph, and find the primitive network
matrices. Neglect mutual values between the elements.

Page 16
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Solution:
The element node incidence matrix, A can be obtained from the given A matrix, by preaugmenting to it an extra column corresponding to the reference node, as under.

Based on the conventional definitions of the elements of A , the oriented graph can be
formed as under:

Page 17
vtu.allsyllabus.com

co

www.allsyllabus.com

yl
la
bu
s.

Fig. E4 Oriented Graph

w
w

w
.a

lls

Thus the primitive network matrices are square, symmetric and diagonal matrices of
order e=no. of elements = 6. They are obtained as follows.

Example-5: Consider three passive elements whose data is given in Table E5 below.
Form the primitive network impedance matrix.
Page 18
vtu.allsyllabus.com

www.allsyllabus.com

co

Table E5

w
w

w
.a

lls

yl
la
bu
s.

Solution:

Note:
The size of [z] is e e, where e= number of elements,
The diagonal elements are the self impedances of the elements
The off-diagonal elements are mutual impedances between the corresponding elements.
Matrices [z] and [y] are inter-invertible.

Page 19
vtu.allsyllabus.com

www.allsyllabus.com

Expected questions:

w
w

w
.a

lls

yl
la
bu
s.

co

1. Define the following and give an illustrative example: i) tree and cotree ii) Basic loops iii) Basic cut sets iv) primitive network v) Bus frame
of reference
2. Derive an expression for obtaining Y-bus using singular
transformations.
3. The Bus incidence matrix for a network of 8 elements and 5 nodes is
given below. Reconstruct the oriented graph, by forming the element
node incidence matrix.
4. Define the following terms with examples: i) Graph ii) branch-path
incidence matrix.
5. The bus incidence matrix is given below. Draw the oriented graph.
Obtain the augmented loop incidence matrix.
6. Define with examples: i) oriented connected graph ii) primitive
network
7. Show that for a power system Ybus=At[y]A Where A is bus incidence
and [y] is primitive admittance matrix
8. a. Define the following terms with the illustrative examples: i)
oriented graph ii) Tree iii) Co-tree
9. Derive the expression for Ybus using singular transformation
10. What is primitive network? Explain with circuit and equations the
significance of it in both impedance and admittance forms.
11. Derive an expression for obtaining Y-bus using singular
transformations.

Page 20
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-2
NETWORK MATRICES
2. FORMATION OF YBUS AND ZBUS
The bus admittance matrix, YBUS plays a very important role in computer aided power
system analysis. It can be formed in practice by either of the methods as under:
1. Rule of Inspection
2. Singular Transformation
3. Non-Singular Transformation
4. ZBUS Building Algorithms, etc.

The performance equations of a given power system can be considered in three different
frames of reference as discussed below:

yl
la
bu
s.

co

Frames of Reference:
Bus Frame of Reference: There are b independent equations (b = no. of buses) relating the
bus vectors of currents and voltages through the bus impedance matrix and bus admittance
matrix:

lls

EBUS = ZBUS IBUS

IBUS = YBUS EBUS
(9)
Branch Frame of Reference: There are b independent equations (b = no. of branches of a
selected Tree sub-graph of the system Graph) relating the branch vectors of currents and
voltages through the branch impedance matrix and branch admittance matrix:

w
w

w
.a

EBR = ZBR IBR

IBR = YBR EBR
(10)
Loop Frame of Reference: There are b independent equations (b = no. of branches of a
selected Tree sub-graph of the system Graph) relating the branch vectors of currents and
voltages through the branch impedance matrix and branch admittance matrix:
ELOOP = ZLOOP ILOOP
ILOOP = YLOOP ELOOP

(11)

Of the various network matrices refered above, the bus admittance matrix (YBUS) and the
bus impedance matrix (ZBUS) are determined for a given power system by the rule of
inspection as explained next.
2.1 Rule of Inspection
Consider the 3-node admittance network as shown in figure5. Using the basic branch
relation: I = (YV), for all the elemental currents and applying Kirchhoffs Current
Law principle at the nodal points, we get the relations as under:
At node 1: I1 =Y1V1 + Y3 (V1-V3) + Y6 (V1 V2)
At node 2: I2 =Y2V2 + Y5 (V2-V3) + Y6 (V2 V1)
At node 3: 0 = Y3 (V3-V1) + Y4V3 + Y5 (V3 V2)

(12)
Page 21

vtu.allsyllabus.com

www.allsyllabus.com

Fig. 3 Example System for finding YBUS

yl
la
bu
s.

co

These are the performance equations of the given network in admittance form and
they can be represented in matrix form as:

w
w

w
.a

lls

In other words, the relation of equation (9) can be represented in the form
IBUS = YBUS EBUS
(14)
Where, YBUS is the bus admittance matrix, IBUS & EBUS are the bus current and bus
voltage vectors respectively. By observing the elements of the bus admittance matrix,
YBUS of equation (13), it is observed that the matrix elements can as well be obtained by
a simple inspection of the given system diagram:
Diagonal elements: A diagonal element (Yii) of the bus admittance matrix, YBUS, is
equal to the sum total of the admittance values of all the elements incident at the bus/node
i,
Off Diagonal elements: An off-diagonal element (Yij) of the bus admittance matrix,
YBUS, is equal to the negative of the admittance value of the connecting element present
between the buses I and j, if any. This is the principle of the rule of inspection. Thus the
algorithmic equations for the rule of inspection are obtained as:
Yii = S yij (j = 1,2,.n)
Yij = - yij (j = 1,2,.n)

(15)

For i = 1,2,.n, n = no. of buses of the given system, yij is the admittance of element
connected between buses i and j and yii is the admittance of element connected between
bus i and ground (reference bus).

Page 22
vtu.allsyllabus.com

www.allsyllabus.com

2.2 Bus impedance matrix

In cases where, the bus impedance matrix is also required, it cannot be formed by direct
inspection of the given system diagram. However, the bus admittance matrix determined
by the rule of inspection following the steps explained above, can be inverted to obtain the
bus impedance matrix, since the two matrices are interinvertible.
Note: It is to be noted that the rule of inspection can be applied only to those power
systems that do not have any mutually coupled elements.
Examples on Rule of Inspection:

lls

yl
la
bu
s.

co

Example 6: Obtain the bus admittance matrix for the admittance network shown aside by
the rule of inspection

w
w

w
.a

Example 7: Obtain YBUS for the impedance network shown aside by the rule of
inspection. Also, determine YBUS for the reduced network after eliminating the eligible
unwanted node. Draw the resulting reduced system diagram.

Page 23
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

2.3 SINGULAR TRANSFORMATIONS

The primitive network matrices are the most basic matrices and depend purely on the
impedance or admittance of the individual elements. However, they do not contain any
information about the behaviour of the interconnected network variables. Hence, it is
necessary to transform the primitive matrices into more meaningful matrices which can
relate variables of the interconnected network.

Page 24
vtu.allsyllabus.com

www.allsyllabus.com

Bus admittance matrix, YBUS and Bus impedance matrix, ZBUS

In the bus frame of reference, the performance of the interconnected network is described
by n independent nodal equations, where n is the total number of buses (n+1nodes are
present, out of which one of them is designated as the reference node).
For example a 5-bus system will have 5 external buses and 1 ground/ ref. bus). The
performance equation relating the bus voltages to bus current injections in bus frame of
reference in admittance form is given by

IBUS = YBUS EBUS

(17)
Where EBUS = vector of bus voltages measured with respect to reference bus
IBUS = Vector of currents injected into the bus
YBUS = bus admittance matrix
The performance equation of the primitive network in admittance form is given by
i + j = [y] v
Pre-multiplying by At (transpose of A), we obtain
(18)

co

At i +At j = At [y] v
However, as per equation (4),

yl
la
bu
s.

At i =0,
since it indicates a vector whose elements are the algebraic sum of element currents
incident at a bus, which by Kirchhoffs law is zero. Similarly, At j gives the algebraic sum
of all source currents incident at each bus and this is nothing but the total current injected
at the bus. Hence,
(19)
(20)

IBUS = At [y] A EBUS

Comparing (21) with (17) we obtain,

(21)

w
w

w
.a

lls

At j = IBUS
Thus from (18) we have, IBUS = At [y] v
However, from (5), we have
v =A EBUS
And hence substituting in (20) we get,

YBUS = At [y] A
(22)
The bus incidence matrix is rectangular and hence singular. Hence, (22) gives a singular
transformation of the primitive admittance matrix [y]. The bus impedance matrix is given
by ,
ZBUS = YBUS-1
(23)
Note: This transformation can be derived using the concept of power invariance, however,
since the transformations are based purely on KCL and KVL, the transformation will
obviously be power invariant.

Page 25
vtu.allsyllabus.com

www.allsyllabus.com

Examples on Singular Transformation:

yl
la
bu
s.

co

Example 8: For the network of Fig E8, form the primitive matrices [z] & [y] and obtain
the bus admittance matrix by singular transformation. Choose a Tree T(1,2,3). The data is
given in Table E8.

Fig E8 System for Example-8

w
w

w
.a

lls

Table E8: Data for Example-8

Page 26
vtu.allsyllabus.com

www.allsyllabus.com

Solution:
The bus incidence matrix is formed taking node 1 as the reference bus.

yl
la
bu
s.

co

The primitive incidence matrix is given by

w
w

w
.a

lls

The primitive admittance matrix [y] = [z]-1 and given by,

The bus admittance matrix by singular transformation is obtained as

Page 27
vtu.allsyllabus.com

www.allsyllabus.com

SUMMARY
The formulation of the mathematical model is the first step in obtaining the solution of any
electrical network. The independent variables can be either currents or voltages.
Correspondingly, the elements of the coefficient matrix will be impedances or
admittances.
Network equations can be formulated for solution of the network using graph theory,
independent of the nature of elements. In the graph of a network, the tree-branches and
links are distinctly identified. The complete information about the interconnection of the
network, with the directions of the currents is contained in the bus incidence matrix.

yl
la
bu
s.

co

The information on the nature of the elements which form the interconnected network is
contained in the primitive impedance matrix. A primitive element can be represented in
impedance form or admittance form. In the bus frame of reference, the performance of the
interconnected system is described by (n-1) nodal equations, where n is the number of
nodes. The bus admittance matrix and the bus impedance matrix relate the bus voltages
and currents. These matrices can be obtained from the primitive impedance and
admittance matrices.

FORMATION OF BUS IMPEDANCE MATRIX

lls

2.4 NODE ELIMINATION BY MATRIX ALGEBRA

w
.a

Nodes can be eliminated by the matrix manipulation of the standard node equations.
However, only those nodes at which current does not enter or leave the network can be
considered for such elimination. Such nodes can be eliminated either in one group or by
taking the eligible nodes one after the other for elimination, as discussed next.

w
w

CASE-A: Simultaneous Elimination of Nodes:

Consider the performance equation of the given network in bus frame of reference in
admittance form for a n-bus system, given by:
IBUS = YBUS EBUS

(1)

Where IBUS and EBUS are n-vectors of injected bus current and bus voltages and YBUS
is the square, symmetric, coefficient bus admittance matrix of order n. Now, of the n buses
present in the system, let p buses be considered for node elimination so that the reduced
system after elimination of p nodes would be retained with m (= n-p) nodes only. Hence
the corresponding performance equation would be similar to (1) except that the coefficient
matrix would be of order m now, i.e.,
IBUS = YBUSnew EBUS

(2)

Where YBUSnew is the bus admittance matrix of the reduced network and the vectors
Page 28
vtu.allsyllabus.com

www.allsyllabus.com

IBUS and EBUS are of order m. It is assumed in (1) that IBUS and EBUS are obtained
with their elements arranged such that the elements associated with p nodes to be
eliminated are in the lower portion of the vectors. Then the elements of YBUS also get
located accordingly so that (1) after matrix partitioning yields,

yl
la
bu
s.

IBUS-m = YA EBUS-m + YB EBUS-p

IBUS-p = YC EBUS-m + YD EBUS-p = 0
(4)

co

Where the self and mutual values of YA and YD are those identified only with the nodes
to be retained and removed respectively and YC=YBt is composed of only the
corresponding mutual admittance values, that are common to the nodes m and p.
Now, for the p nodes to be eliminated, it is necessary that, each element of the vector
IBUS-p should be zero. Thus we have from (3):

w
w

w
.a

lls

Solving,
EBUS-p = - YD-1YC EBUS-m
(5)
Thus, by simplification, we obtain an expression similar to (2) as,
IBUS-m = {YA - YBYD-1YC} EBUS-m
(6)
Thus by comparing (2) and (6), we get an expression for the new bus admittance matrix in
terms of the sub-matrices of the original bus admittance matrix as:
YBUSnew = {YA YBYD -1YC}
(7)
This expression enables us to construct the given network with only the necessary nodes
retained and all the unwanted nodes/buses eliminated. However, it can be observed from
(7) that the expression involves finding the inverse of the sub-matrix YD (of order p). This
would be computationally very tedious if p, the nodes to be eliminated is very large,
especially for real practical systems. In such cases, it is more advantageous to eliminate
the unwanted nodes from the given network by considering one node only at a time for
elimination, as discussed next.
CASE-B: Separate Elimination of Nodes:
Here again, the system buses are to be renumbered, if necessary, such that the node to be
removed always happens to be the last numbered one. The sub-matrix YD then would be a
single element matrix and hence it inverse would be just equal to its own reciprocal value.
Thus the generalized algorithmic equation for finding the elements of the new bus
admittance matrix can be obtained from (6) as,
Yij new = Yij old Yin Ynj / Ynn " i,j = 1,2, n.

(8)

Page 29
vtu.allsyllabus.com

www.allsyllabus.com

Each element of the original matrix must therefore be modified as per (7). Further, this
procedure of eliminating the last numbered node from the given system of n nodes is to be
iteratively repeated p times, so as to eliminate all the unnecessary p nodes from the
original system.

Examples on Node elimination:

yl
la
bu
s.

co

Example-1: Obtain YBUS for the impedance network shown below by the rule of
inspection. Also, determine YBUS for the reduced network after eliminating the eligible
unwanted node. Draw the resulting reduced system diagram.

w
w

w
.a

lls

The admittance equivalent network is as follows:

The bus admittance matrix is obtained by RoI as:

The reduced matrix after elimination of node 3 from the given system is determined as per the
equation:

Page 30
vtu.allsyllabus.com

www.allsyllabus.com

yl
la
bu
s.

co

Alternatively,

w
w

w
.a

lls

Thus the reduced network can be obtained again by the rule of inspection as shown below.

Example-2: Obtain YBUS for the admittance network shown below by the rule of
inspection. Also, determine YBUS for the reduced network after eliminating the eligible
unwanted node. Draw the resulting reduced system diagram.
Page 31
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Thus the reduced system of two nodes can be drawn by the rule of inspection as under:
Page 32
vtu.allsyllabus.com

2.5 ZBUS building

yl
la
bu
s.

FORMATION OF BUS IMPEDANCE MATRIX

co

www.allsyllabus.com

w
.a

lls

The bus impedance matrix is the inverse of the bus admittance matrix. An alternative
method is possible, based on an algorithm to form the bus impedance matrix directly from
system parameters and the coded bus numbers. The bus impedance matrix is formed
adding one element at a time to a partial network of the given system. The performance
equation of the network in bus frame of reference in impedance form using the currents as
independent variables is given in matrix form by

w
w

When expanded so as to refer to a n bus system, (9) will be of the form

Now assume that the bus impedance matrix Zbus is known for a partial network of m
buses and a known reference bus. Thus, Zbus of the partial network is of dimension mxm.
If now a new element is added between buses p and q we have the following two
possibilities:
(i)
p is an existing bus in the partial network and q is a new bus; in this case p-q is a
branch added to the p-network as shown in Fig 1a, and
Page 33
vtu.allsyllabus.com

www.allsyllabus.com

w
.a

lls

yl
la
bu
s.

co

both p and q are buses existing in the partial network; in this case p-q is a link
added to the p-network as shown in Fig 1b.

w
w

(ii)

Page 34
vtu.allsyllabus.com

www.allsyllabus.com

If the added element ia a branch, p-q, then the new bus impedance matrix would be of
order m+1, and the analysis is confined to finding only the elements of the new row and
column (corresponding to bus-q) introduced into the original matrix. If the added element
ia a link, p-q, then the new bus impedance matrix will remain unaltered with regard to its
order. However, all the elements of the original matrix are updated to take account of the
effect of the link added.

yl
la
bu
s.

co

ADDITION OF A BRANCH
Consider now the performance equation of the network in impedance form with the added
branch p-q, given by

It is assumed that the added branch p-q is mutually coupled with some elements of the
partial network and since the network has bilateral passive elements only, we have

w
.a

lls

Vector ypq-rs is not equal to zero and Zij= Zji " i,j=1,2,m,q
(12)

w
w

To find Zqi:
The elements of last row-q and last column-q are determined by injecting a current of 1.0
pu at the bus-i and measuring the voltage of the bus-q with respect to the reference bus-0,
as shown in Fig.2. Since all other bus currents are zero, we have from (11) that
Ek = Zki Ii = Zki " k = 1, 2,i....p,.m, q
(13)
Hence, Eq = Zqi ; Ep = Zpi
Also, Eq=Ep -vpq ; so that Zqi = Zpi - vpq " i =1, 2,i....p,.m, _q
(14)
To find vpq:
In terms of the primitive admittances and voltages across the elements, the current through
the elements is given by

Page 35
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Fig.2 Calculation for Zqi

Page 36
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 37
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Special Cases
The following special cases of analysis concerning ZBUS building can be considered with
respect to the addition of branch to a p-network.

Page 38
vtu.allsyllabus.com

co

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

ADDITION OF A LINK
Consider now the performance equation of the network in impedance form with the added
link p-l, (p-l being a fictitious branch and l being a fictitious node) given by

It is assumed that the added branch p-q is mutually coupled with some elements of the
partial network and since the network has bilateral passive elements only, we have

To find Zli:
The elements of last row-l and last column-l are determined by injecting a current of 1.0
pu at the bus-i and measuring the voltage of the bus-q with respect to the reference bus-0,
as shown in Fig.3. Further, the current in the added element is made zero by connecting a
voltage source, el in series with element p-q, as shown. Since all other bus currents are
zero, we have from (25) that

Page 39
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

To find vpq:
In terms of the primitive admittances and voltages across the elements, the current through
the elements is given by

Page 40
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 41
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

From (39), it is thus observed that, when a link is added to a ref. bus, then the situation is
similar to adding a branch to a fictitious bus and hence the following steps are followed:
Page 42
vtu.allsyllabus.com

www.allsyllabus.com

1. The element is added similar to addition of a branch (case-b) to obtain the new matrix
of order m+1.
2. The extra fictitious node, l is eliminated using the node elimination algorithm.
Case (d): If there is no mutual coupling, then elements of pq rs y , are zero. Further, if p is
not the reference node, then

2.6 MODIFICATION OF ZBUS FOR NETWORK CHANGES

lls

yl
la
bu
s.

co

An element which is not coupled to any other element can be removed easily. The Zbus is
modified as explained in sections above, by adding in parallel with the element (to be
removed), a link whose impedance is equal to the negative of the impedance of the
element to be removed. Similarly, the impedance value of an element which is not coupled
to any other element can be changed easily. The Zbus is modified again as explained in
sections above, by adding in parallel with the element (whose impedance is to be
changed), a link element of impedance value chosen such that the parallel equivalent
impedance is equal to the desired value of impedance. When mutually coupled elements
are removed, the Zbus is modified by introducing appropriate changes in the bus currents
of the original network to reflect the changes introduced due to the removal of the
elements.

w
.a

Examples on ZBUS building

w
w

Example 1: For the positive sequence network data shown in table below, obtain ZBUS
by building procedure.

Solution:
The given network is as shown below with the data marked on it. Assume the elements to
be added as per the given sequence: 0-1, 0-3, 1-2, and 2-3.

Page 43
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Fig. E1: Example System

Consider building ZBUS as per the various stages of building through the consideration of
the corresponding partial networks as under:
Step-1: Add element1 of impedance 0.25 pu from the external node-1 (q=1) to internal
ref. node-0 (p=0). (Case-a), as shown in the partial network;

Step-2: Add element2 of impedance 0.2 pu from the external node-3 (q=3) to internal ref.
node-0 (p=0). (Case-a), as shown in the partial network;

Page 44
vtu.allsyllabus.com

co

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

Step-3: Add element3 of impedance 0.08 pu from the external node-2 (q=2) to internal node1 (p=1). (Case-b), as shown in the partial network;

Step-4: Add element4 of impedance 0.06 pu between the two internal nodes, node-2
(p=2) to node-3 (q=3). (Case-d), as shown in the partial network;

Page 45
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

lls

The fictitious node l is eliminated further to arrive at the final impedance matrix as under:

Page 46
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 47
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 48
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Solution: The specified system is considered with the reference node denoted by node-0.
By its inspection, we can obtain the bus impedance matrix by building procedure by
following the steps through the p-networks as under:
Step1: Add branch 1 between node 1 and reference node. (q =1, p = 0)

Step2: Add branch 2, between node 2 and reference node. (q = 2, p = 0).

Page 49
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Step3: Add branch 3, between node 1 and node 3 (p = 1, q = 3)

Step 4: Add element 4, which is a link between node 1 and node 2. (p = 1, q = 2)

Page 50
vtu.allsyllabus.com

www.allsyllabus.com

yl
la
bu
s.

co

Now the extra node-l has to be eliminated to obtain the new matrix of step-4, using the
algorithmic relation:

w
w

w
.a

lls

Step 5: Add link between node 2 and node 3 (p = 2, q=3)

Page 51
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 52
vtu.allsyllabus.com

www.allsyllabus.com

Expected questions:

w
w

w
.a

lls

yl
la
bu
s.

co

1. Obtain the general expressions for Zbus building algorithm when a

branch is added to the partial network.
2. For the network shown. Obtain the Ybus by singular
transformation analysis. The line data is as follows
3. Obtain the general expressions for Zbus building algorithm when a
link is added to the partial network.
4. Prepare the Zbus for the system shown using Zbus building algorithm
5. Prepare the Zbus for the system shown using Zbus building algorithm
6. Explain the formation of Zbus using Zbus building algorithm
7. What is a primitive network? Give the representation of a typical
component and arrive at the performance equations both in
impedance and admittance forms.

Page 53
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-3 &4
LOAD FLOW STUDIES
3.1 REVIEW OF NUMERICAL SOLUTION OF EQUATIONS

co

yl
la
bu
s.

1. Solution Linear equations:

* Direct methods:
- Cramers (Determinant) Method,
- Gauss Elimination Method (only for smaller systems),
- LU Factorization (more preferred method), etc.

The numerical analysis involving the solution of algebraic simultaneous equations forms
the basis for solution of the performance equations in computer aided electrical power
system analyses, such as during linear graph analysis, load flow analysis (nonlinear
equations), transient stability studies (differential equations), etc. Hence, it is necessary to
review the general forms of the various solution methods with respect to all forms of
equations, as under:

* Iterative methods:
- Gauss Method
- Gauss-Siedel Method (for diagonally dominant systems)

w
.a

lls

3. Solution of Nonlinear equations:

Iterative methods only:
- Gauss-Siedel Method (for smaller systems)
- Newton-Raphson Method (if corrections for variables are small)

w
w

4. Solution of differential equations:

Iterative methods only:
- Euler and Modified Euler method,
- RK IV-order method,
- Milnes predictor-corrector method, etc.
It is to be observed that the nonlinear and differential equations can be solved only by the
iterative methods. The iterative methods are characterized by the various performance
features as under:
_ Selection of initial solution/ estimates
_ Determination of fresh/ new estimates during each iteration
_ Selection of number of iterations as per tolerance limit
_ Time per iteration and total time of solution as per the solution method selected
_ Convergence and divergence criteria of the iterative solution
_ Choice of the Acceleration factor of convergence, etc.

Page 54
vtu.allsyllabus.com

www.allsyllabus.com

A comparison of the above solution methods is as under:

In general, the direct methods yield exact or accurate solutions. However, they are suited
for only the smaller systems, since otherwise, in large systems, the possible round-off
errors make the solution process inaccurate. The iterative methods are more useful when
the diagonal elements of the coefficient matrix are large in comparison with the off
diagonal elements. The round-off errors in these methods are corrected at the successive
steps of the iterative process.The Newton-Raphson method is very much useful for
solution of non linear equations, if all the values of the corrections for the unknowns are
very small in magnitude and the initial values of unknowns are selected to be reasonably
closer to the exact solution.
3.2 LOAD FLOW STUDIES

w
w

w
.a

lls

yl
la
bu
s.

co

Introduction: Load flow studies are important in planning and designing future expansion
of power systems. The study gives steady state solutions of the voltages at all the buses,
for a particular load condition. Different steady state solutions can be obtained, for
different operating conditions, to help in planning, design and operation of the power
system. Generally, load flow studies are limited to the transmission system, which
involves bulk power transmission. The load at the buses is assumed to be known. Load
flow studies throw light on some of the important aspects of the system operation, such as:
violation of voltage magnitudes at the buses, overloading of lines, overloading of
generators, stability margin reduction, indicated by power angle differences between buses
linked by a line, effect of contingencies like line voltages, emergency shutdown of
generators, etc. Load flow studies are required for deciding the economic operation of the
power system. They are also required in transient stability studies. Hence, load flow
studies play a vital role in power system studies. Thus the load flow problem consists of
finding the power flows (real and reactive) and voltages of a network for given bus
conditions. At each bus, there are four quantities of interest to be known for further
analysis: the real and reactive power, the voltage magnitude and its phase angle. Because
of the nonlinearity of the algebraic equations, describing the given power system, their
solutions are obviously, based on the iterative methods only. The constraints placed on the
load flow solutions could be:
_ The Kirchhoffs relations holding good,
_ Capability limits of reactive power sources,
_ Tap-setting range of tap-changing transformers,
_ Specified power interchange between interconnected systems,
_ Selection of initial values, acceleration factor, convergence limit, etc.

3.3 Classification of buses for LFA: Different types of buses are present based on
the specified and unspecified variables at a given bus as presented in the table below:

Page 55
vtu.allsyllabus.com

www.allsyllabus.com

Table 1. Classification of buses for LFA

lls

yl
la
bu
s.

co

Importance of swing bus: The slack or swing bus is usually a PV-bus with the largest
capacity generator of the given system connected to it. The generator at the swing bus
supplies the power difference between the specified power into the system at the other
buses and the total system output plus losses. Thus swing bus is needed to supply the
additional real and reactive power to meet the losses. Both the magnitude and phase angle
of voltage are specified at the swing bus, or otherwise, they are assumed to be equal to 1.0
p.u. and 00 , as per flat-start procedure of iterative
solutions. The real and reactive powers at the swing bus are found by the computer routine
as part of the load flow solution process. It is to be noted that the source at the swing bus is
a perfect one, called the swing machine, or slack machine. It is voltage regulated, i.e., the
magnitude of voltage fixed. The phase angle is the system reference phase and hence is
fixed. The generator at the swing bus has a torque angle and excitation which vary or
swing as the demand changes. This variation is such as to produce fixed voltage.

Importance of YBUS based LFA:

w
w

w
.a

The majority of load flow programs employ methods using the bus admittance matrix, as
this method is found to be more economical. The bus admittance matrix plays a very
important role in load flow analysis. It is a complex, square and symmetric matrix and
hence only n(n+1)/2 elements of YBUS need to be stored for a n-bus system. Further, in
the YBUS matrix, Yij = 0, if an incident element is not present in the system connecting
the buses i and j. since in a large power system, each bus is connected only to a fewer
buses through an incident element, (about 6-8), the coefficient matrix, YBUS of such
systems would be highly sparse, i.e., it will have many zero valued elements in it. This is
defined by the sparsity of the matrix, as under:

The percentage sparsity of YBUS, in practice, could be as high as 80-90%, especially

for very large, practical power systems. This sparsity feature of YBUS is extensively used
in reducing the load flow calculations and in minimizing the memory required to store the

Page 56
vtu.allsyllabus.com

www.allsyllabus.com

coefficient matrices. This is due to the fact that only the non-zero elements YBUS can be
stored during the computer based implementation of the schemes, by adopting the suitable
optimal storage schemes. While YBUS is thus highly sparse, its inverse, ZBUS, the bus
impedance matrix is not so. It is a FULL matrix, unless the optimal bus ordering schemes
are followed before proceeding for load flow analysis.
3.4 THE LOAD FLOW PROBLEM

w
w

w
.a

lls

yl
la
bu
s.

co

Here, the analysis is restricted to a balanced three-phase power system, so that the analysis
can be carried out on a single phase basis. The per unit quantities are used for all
quantities. The first step in the analysis is the formulation of suitable equations for the
power flows in the system. The power system is a large interconnected system, where
various buses are connected by transmission lines. At any bus, complex power is injected
into the bus by the generators and complex power is drawn by the loads. Of course at any
bus, either one of them may not be present. The power is transported from one bus to other
via the transmission lines. At any bus i, the complex power Si (injected), shown in figure
1, is defined as

where Si = net complex power injected into bus i, SGi = complex power injected by the
generator at bus i, and SDi = complex power drawn by the load at bus i. According to
conservation of complex power, at any bus i, the complex power injected into the bus must
be equal to the sum of complex power flows out of the bus via the transmission lines.
Hence,
Si = _Sij " i = 1, 2, ..n
(3)
where Sij is the sum over all lines connected to the bus and n is the number of buses in the
system (excluding the ground). The bus current injected at the bus-i is defined as
Ii = IGi IDi " i = 1, 2, ..n

(4)

Page 57
vtu.allsyllabus.com

www.allsyllabus.com

where IGi is the current injected by the generator at the bus and IDi is the current drawn
by the load (demand) at that bus. In the bus frame of reference

w
w

w
.a

lls

yl
la
bu
s.

co

IBUS = YBUS VBUS

(5)

Page 58
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Equations (9)-(10) and (13)-(14) are the power flow equations or the load flow
equations in two alternative forms, corresponding to the n-bus system, where each bus-i is
characterized by four variables, Pi, Qi, |Vi|, and di. Thus a total of 4n variables are
Page 59
vtu.allsyllabus.com

www.allsyllabus.com

involved in these equations. The load flow equations can be solved for any 2n unknowns,
if the other 2n variables are specified. This establishes the need for classification of buses
of the system for load flow analysis into: PV bus, PQ bus, etc.
3.4 DATA FOR LOAD FLOW
Irrespective of the method used for the solution, the data required is common for any load
flow. All data is normally in pu. The bus admittance matrix is formulated from these data.
The various data required are as under:
System data: It includes: number of buses-n, number of PV buses, number of loads,
number of transmission lines, number of transformers, number of shunt elements, the slack
bus number, voltage magnitude of slack bus (angle is generally taken as 0o), tolerance
limit, base MVA, and maximum permissible number of iterations.

co

Generator bus data: For every PV bus i, the data required includes the bus number,
active power generation PGi, the specified voltage magnitude i sp V , , minimum reactive
power limit Qi,min, and maximum reactive power limit Qi,max.

yl
la
bu
s.

Load data: For all loads the data required includes the the bus number, active power
demand PDi, and the reactive power demand QDi.
Transmission line data: For every transmission line connected between buses i and k the
data includes the starting bus number i, ending bus number k,.resistance of the line,
reactance of the line and the half line charging admittance.

w
.a

lls

Transformer data:
For every transformer connected between buses i and k the data to be given includes: the
starting bus number i, ending bus number k, resistance of the transformer, reactance of the
transformer, and the off nominal turns-ratio a.

w
w

Shunt element data: The data needed for the shunt element includes the bus number
where element is connected, and the shunt admittance (Gsh + j Bsh).
GAUSS SEIDEL (GS) METHOD
The GS method is an iterative algorithm for solving non linear algebraic equations. An
initial solution vector is assumed, chosen from past experiences, statistical data or from
practical considerations. At every subsequent iteration, the solution is updated till
convergence is reached. The GS method applied to power flow problem is as discussed
below.
Case (a): Systems with PQ buses only:
Initially assume all buses to be PQ type buses, except the slack bus. This means that (n1)
complex bus voltages have to be determined. For ease of programming, the slack bus is
generally numbered as bus-1. PV buses are numbered in sequence and PQ buses are
ordered next in sequence. This makes programming easier, compared to random ordering
of buses. Consider the expression for the complex power at bus-i, given from (7), as:

Page 60
vtu.allsyllabus.com

lls

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

Equation (17) is an implicit equation since the unknown variable, appears on both sides of
the equation. Hence, it needs to be solved by an iterative technique. Starting from an initial
estimate of all bus voltages, in the RHS of (17) the most recent values of the bus voltages
is substituted. One iteration of the method involves computation of all the bus voltages. In
GaussSeidel method, the value of the updated voltages are used in the computation of
subsequent voltages in the same iteration, thus speeding up convergence. Iterations are
carried out till the magnitudes of all bus voltages do not change by more than the tolerance
value. Thus the algorithm for GS method is as under:
3.5 Algorithm for GS method
1. Prepare data for the given system as required.
2. Formulate the bus admittance matrix YBUS. This is generally done by the rule of
inspection.
3. Assume initial voltages for all buses, 2,3,n. In practical power systems, the magnitude
of the bus voltages is close to 1.0 p.u. Hence, the complex bus voltages at all (n-1) buses
(except slack bus) are taken to be 1.0 00. This is normally refered as the flat start solution.
4. Update the voltages. In any (k +1)st iteration, from (17) the voltages are given by

Page 61
vtu.allsyllabus.com

www.allsyllabus.com

Here note that when computation is carried out for bus-i, updated values are already
available for buses 2,3.(i-1) in the current (k+1)st iteration. Hence these values are used.
For buses (i+1)..n, values from previous, kth iteration are used.

yl
la
bu
s.

co

Where,e is the tolerance value. Generally it is customary to use a value of 0.0001 pu.
Compute slack bus power after voltages have converged using (15) [assuming bus 1 is
slack bus].

lls

7. Compute all line flows.

8. The complex power loss in the line is given by Sik + Ski. The total loss in the system is
calculated by summing the loss over all the lines.

w
w

w
.a

Case (b): Systems with PV buses also present:

At PV buses, the magnitude of voltage and not the reactive power is specified. Hence it is
needed to first make an estimate of Qi to be used in (18). From (15) we have

Page 62
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

lls

Case (c): Systems with PV buses with reactive power generation limits specified:
In the previous algorithm if the Q limit at the voltage controlled bus is violated during any
iteration, i.e (k +1) i Q computed using (21) is either less than Qi, min or greater than
Qi,max, it means that the voltage cannot be maintained at the specified value due to lack
of reactive power support. This bus is then treated as a PQ bus in the (k+1)st iteration and
the voltage is calculated with the value of Qi set as follows:

If in the subsequent iteration, if Qi falls within the limits, then the bus can be switched
back to PV status.
Acceleration of convergence
It is found that in GS method of load flow, the number of iterations increase with increase
in the size of the system. The number of iterations required can be reduced if the
correction in voltage at each bus is accelerated, by multiplying with a constant , called
the acceleration factor. In the (k+1)st iteration we can let

where is a real number. When

=1, the value of (k +1) is the computed value. If 1<
<2 then the value computed is extrapolated. Generally _ is taken between 1.2 to 1.6, for
GS load flow procedure. At PQ buses (pure load buses) if the voltage magnitude violates

Page 63
vtu.allsyllabus.com

www.allsyllabus.com

the limit, it simply means that the specified reactive power demand cannot be supplied,
with the voltage maintained within acceptable limits.

w
w

w
.a

lls

yl
la
bu
s.

co

3.6 Examples on GS load flow analysis:

Page 64
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Since the difference in the voltage magnitudes is less than 10-6 pu, the iterations can be
stopped. To compute line flow

Page 65
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

lls

The total loss in the line is given by S12 + S21 = j 0.133329 pu Obviously, it is observed
that there is no real power loss, since the line has no resistance.

Page 66
vtu.allsyllabus.com

www.allsyllabus.com

Example-2:

w
w

w
.a

lls

yl
la
bu
s.

co

For the power system shown in fig. below, with the data as given in tables below, obtain
the bus voltages at the end of first iteration, by applying GS method.

Page 67
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 68
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 69
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Example-3:
Obtain the load flow solution at the end of first iteration of the system with data as given
below. The solution is to be obtained for the following cases
(i) All buses except bus 1 are PQ Buses
(ii) Bus 2 is a PV bus whose voltage magnitude is specified as 1.04 pu
Page 70
vtu.allsyllabus.com

www.allsyllabus.com

w
.a

lls

yl
la
bu
s.

co

Bus 2 is PV bus, with voltage magnitude specified as 1.04 and 0.25_Q2_1.0 pu.

w
w

(iii)

Page 71
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 72
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 73
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 74
vtu.allsyllabus.com

www.allsyllabus.com

Limitations of GS load flow analysis

yl
la
bu
s.

co

GS method is very useful for very small systems. It is easily adoptable, it can be
generalized and it is very efficient for systems having less number of buses. However, GS
LFA fails to converge in systems with one or more of the features as under:
Systems having large number of radial lines
Systems with short and long lines terminating on the same bus
Systems having negative values of transfer admittances
Systems with heavily loaded lines, etc.

w
w

w
.a

lls

GS method successfully converges in the absence of the above problems. However,

convergence also depends on various other set of factors such as: selection of slack bus,
initial solution, acceleration factor, tolerance limit, level of accuracy of results needed,
type and quality of computer/ software used, etc.

Page 75
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

3.7

Page 76
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

Here, the matrix [J] is called the Jacobian matrix. The vector of unknown variables is
updated using (30). The process is continued till the difference between two successive
iterations is less than the tolerance value.

Page 77
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 78
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 79
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

3. 8

Page 80
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 81
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 82
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 83
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 84
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 85
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 86
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 87
vtu.allsyllabus.com

FINAL WORD

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

w
w

In this chapter, the load flow problem, also called as the power flow problem, has been
considered in detail. The load flow solution gives the complex voltages at all the buses and
the complex power flows in the lines. Though, algorithms are available using the
impedance form of the equations, the sparsity of the bus admittance matrix and the ease of
building the bus admittance matrix, have made algorithms using the admittance form of
equations more popular. The most popular methods are the Gauss-Seidel method, the
Newton-Raphson method and the Fast Decoupled Load Flow method. These methods
have been discussed in detail with illustrative examples. In smaller systems, the ease of
programming and the memory requirements, make GS method attractive. However, the
computation time increases with increase in the size of the system. Hence, in large systems
NR and FDLF methods are more popular. There is a trade off between various
requirements like speed, storage, reliability, computation time, convergence characteristics
etc. No single method has all the desirable features. However, NR method is most popular
because of its versatility, reliability and accuracy.

Page 88
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-3&4

w
w

w
.a

lls

yl
la
bu
s.

co

1. Using generalized algorithm expressions for each case of analysis,

explain the load flow studies procedure, as per the G-S method for
power system having PQ and PV buses, with reactive power
generations constraints.
2. Derive the expression in polar form for the typical diagonal
elements of the sub matrices of the Jacobian in NR method of load
flow analysis.
3. Compare NR and GS method for load flow analysis procedure in
respect of the following i) Time per iteration ii) total solution time
iii) acceleration factor iv)number of iterations
4. Explain briefly fast decoupled load flow solution method for
solving the non linear load flow equations.
5. Draw the flow chart of NR method for load flow analysis.
6. Explain the representation of transformer with fixed tap changing
during the load flow studies
7. What are the assumptions made in fast decoupled load flow
method? Explain the algorithm briefly, through a flow chart.
8. Explain the NR load flow method? Explain the algorithm briefly,
through a flow chart.
9. What is load flow analysis? What is the data required to conduct
load flow analysis? Explain how buses are classified to carry out
load flow analysis in power system. What is the significance of
slack bus.
10. The following is the system data for a load flow analysis. The line
data are shown and also real and reactive powers data is given.
Determine the voltage at the end of first iteration using GS
method. Take =1.6

Page 89
vtu.allsyllabus.com

www.allsyllabus.com

PART-B
UNIT-5 & 6
ECONOMIC OPEARATION OF POWER SYSTEM
5.1 INTRODUCTION

w
.a

lls

yl
la
bu
s.

co

One of the earliest applications of on-line centralized control was to provide a central
facility, to operate economically, several generating plants supplying the loads of the
system. Modern integrated systems have different types of generating plants, such as coal
fired thermal plants, hydel plants, nuclear plants, oil and natural gas units etc. The capital
investment, operation and maintenance costs are different for different types of plants. The
operation economics can again be subdivided into two parts.
i) Problem of economic dispatch, which deals with determining the power output of each
plant to meet the specified load, such that the overall fuel cost is minimized.
ii) Problem of optimal power flow, which deals with minimum loss delivery, where in
the power flow, is optimized to minimize losses in the system. In this chapter we consider
the problem of economic dispatch.
During operation of the plant, a generator may be in one of the following states:
i) Base supply without regulation: the output is a constant.
ii) Base supply with regulation: output power is regulated based on system load.
iii) Automatic non-economic regulation: output level changes around a base setting as area
control error changes.
iv) Automatic economic regulation: output level is adjusted, with the area load and area
control error, while tracking an economic setting.

w
w

Regardless of the units operating state, it has a contribution to the economic operation,
even though its output is changed for different reasons. The factors influencing the cost of
generation are the generator efficiency, fuel cost and transmission losses. The most
efficient generator may not give minimum cost, since it may be located in a place where
fuel cost is high. Further, if the plant is located far from the load centers, transmission
losses may be high and running the plant may become uneconomical. The economic
dispatch problem basically determines the generation of different plants to minimize total
operating cost.
Modern generating plants like nuclear plants, geo-thermal plants etc, may require capital
investment of millions of rupees. The economic dispatch is however determined in terms
of fuel cost per unit power generated and does not include capital investment,
maintenance, depreciation, start-up and shut down costs etc.

Page 90
vtu.allsyllabus.com

www.allsyllabus.com

5.2 PERFORMANCE CURVES

INPUT-OUTPUT CURVE

yl
la
bu
s.

co

This is the fundamental curve for a thermal plant and is a plot of the input in British
thermal units (Btu) per hour versus the power output of the plant in MW as shown in Fig1.

w
w

w
.a

lls

HEAT RATE CURVE

The heat rate is the ratio of fuel input in Btu to energy output in KWh. It is the slope of the
input output curve at any point. The reciprocal of heat rate is called fuel efficiency.
The heat rate curve is a plot of heat rate versus output in MW. A typical plot is shown in
Fig .2

Fig .2 Heat rate curve.

Page 91
vtu.allsyllabus.com

www.allsyllabus.com

INCREMENTAL FUEL RATE CURVE

The incremental fuel rate is equal to a small change in input divided by the
corresponding change in output.

w
.a

lls

yl
la
bu
s.

co

The unit is again Btu / KWh. A plot of incremental fuel rate versus the output is shown in
Fig 3

w
w

Incremental cost curve

The incremental cost is the product of incremental fuel rate and fuel cost (Rs / Btu or $ /
Btu). The curve in shown in Fig. 4. The unit of the incremental fuel cost is Rs / MWh or
$ /MWh.

Page 92
vtu.allsyllabus.com

co

www.allsyllabus.com

w
.a

lls

The incremental fuel cost is given by

yl
la
bu
s.

In general, the fuel cost Fi for a plant, is approximated as a quadratic function of the
generated output PGi.

w
w

The incremental fuel cost is a measure of how costly it will be produce an increment of
power. The incremental production cost, is made up of incremental fuel cost plus the
incremental cost of labour, water, maintenance etc. which can be taken to be some
percentage of the incremental fuel cost, instead of resorting to a rigorous mathematical
model. The cost curve can be approximated by a linear curve. While there is negligible
operating cost for a hydel plant, there is a limitation on the power output possible. In any
plant, all units normally operate between PGmin, the minimum loading limit, below which
it is technically infeasible to operate a unit and PGmax, which is the maximum output
limit.
5.3 ECONOMIC GENERATION SCHEDULING NEGLECTING LOSSES AND
GENERATOR LIMITS
The simplest case of economic dispatch is the case when transmission losses are
neglected. The model does not consider the system configuration or line impedances.
Since losses are neglected, the total generation is equal to the total demand PD. Consider a
system with ng number of generating plants supplying the total demand PD. If Fi is the

Page 93
vtu.allsyllabus.com

www.allsyllabus.com

cost of plant i in Rs/h, the mathematical formulation of the problem of economic

scheduling can be stated as follows:

yl
la
bu
s.

co

This is a constrained optimization problem, which can be solved by Lagranges method.

w
w

w
.a

lls

LAGRANGE METHOD FOR SOLUTION OF ECONOMIC SCHEDULE

The problem is restated below:

The second equation is simply the original constraint of the problem. The cost of a plant
Page 94
vtu.allsyllabus.com

www.allsyllabus.com

co

Fi depends only on its own output PGi, hence

w
w

w
.a

lls

yl
la
bu
s.

The above equation is called the co-ordination equation. Simply stated, for economic
generation scheduling to meet a particular load demand, when transmission losses are
neglected and generation limits are not imposed, all plants must operate at equal
incremental production costs, subject to the constraint that the total generation be equal to
the demand. From we have

Page 95
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

It can be seen that l is dependent on the demand and the coefficients of the cost function.

Page 96
vtu.allsyllabus.com

www.allsyllabus.com

Example 1.

w
w

w
.a

lls

yl
la
bu
s.

co

The fuel costs of two units are given by

Page 97
vtu.allsyllabus.com

yl
la
bu
s.

co

www.allsyllabus.com

w
w

w
.a

lls

5.4
ECONOMIC SCHEDULE INCLUDING LIMITS ON GENERATOR
(NEGLECTING LOSSES)
The power output of any generator has a maximum value dependent on the rating of the
generator. It also has a minimum limit set by stable boiler operation. The economic
dispatch problem now is to schedule generation to minimize cost, subject to the equality
constraint.

The procedure followed is same as before i.e. the plants are operated with equal
incremental fuel costs, till their limits are not violated. As soon as a plant reaches the limit
(maximum or minimum) its output is fixed at that point and is maintained a constant. The
other plants are operated at equal incremental costs.
Page 98
vtu.allsyllabus.com

www.allsyllabus.com

Example 3

w
.a

lls

yl
la
bu
s.

co

Incremental fuel costs in $ / MWh for two units are given below:

w
w

Now at light loads unit 1 has a higher incremental cost and hence will operate at its lower
limit of 20 MW. Initially, additional load is taken up by unit 2, till such time its
incremental fuel cost becomes equal to 2.2$ / MWh at PG2 = 50 MW. Beyond this, the
two units are operated with equal incremental fuel costs. The contribution of each unit to
meet the demand is obtained by assuming different values of l; When l = 3.1 $ / MWh,
unit 2 operates at its upper limit. Further loads are taken up by unit 1. The computations
are show in Table

Page 99
vtu.allsyllabus.com

co

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

For a particular value of l, PG1 and PG2 are calculated using (8.16). Fig 8.5 Shows plot of
each unit output versus the total plant output.

For any particular load, the schedule for each unit for economic dispatch can be obtained.
Example 4.

Page 100
vtu.allsyllabus.com

www.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

In example 3, what is the saving in fuel cost for the economic schedule compared to the
case where the load is shared equally. The load is 180 MW.
Solution:
From Table it is seen that for a load of 180 MW, the economic schedule is PG1 = 80 MW
and PG2 = 100 MW. When load is shared equally PG1 = PG2 = 90 MW. Hence, the
generation of unit 1 increases from 80 MW to 90 MW and that of unit 2 decreases from
100 MW to 90 MW, when the load is shared equally. There is an increase in cost of unit 1
since PG1 increases and decrease in cost of unit 2 since PG2 decreases.

5.5 ECONOMIC DISPATCH INCLUDING TRANSMISSION LOSSES

When transmission distances are large, the transmission losses are a significant part of the
generation and have to be considered in the generation schedule for economic operation.
The mathematical formulation is now stated as

Page 101
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 102
vtu.allsyllabus.com

www.allsyllabus.com

The minimum operation cost is obtained when the product of the incremental fuel cost and
the penalty factor of all units is the same, when losses are considered. A rigorous general
expression for the loss PL is given by

where Bmn, Bno , Boo called loss coefficients , depend on the load composition. The
assumption here is that the load varies linearly between maximum and minimum values. A
simpler expression is

w
w

w
.a

lls

yl
la
bu
s.

co

The expression assumes that all load currents vary together as a constant complex fraction
of the total load current. Experiences with large systems has shown that the loss of
accuracy is not significant if this approximation is used. An average set of loss coefficients
may be used over the complete daily cycle in the coordination of incremental production
costs and incremental transmission losses. In general, Bmn = Bnm and can be expanded
for a two plant system as

Page 103
vtu.allsyllabus.com

www.allsyllabus.com

Example 5
A generator is supplying a load. An incremental change in load of 4 MW requires
generation to be increased by 6 MW. The incremental cost at the plant bus is Rs 30 /MWh.
What is the incremental cost at the receiving end?

w
w

w
.a

lls

yl
la
bu
s.

co

Solution:

Page 104
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 105
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 106
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 107
vtu.allsyllabus.com

www.allsyllabus.com

5.6 DERIVATION OF TRANSMISSION LOSS FORMULA

w
w

w
.a

lls

yl
la
bu
s.

co

An accurate method of obtaining general loss coefficients has been presented by Kron.
The method is elaborate and a simpler approach is possible by making the following
assumptions:
(i) All load currents have same phase angle with respect to a common reference
(ii) The ratio X / R is the same for all the network branches.
Consider the simple case of two generating plants connected to an arbitrary number of
loads through a transmission network as shown in Fig a

Page 108
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 109
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 110
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 111
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 112
vtu.allsyllabus.com

w
w

w
.a

lls

yl
la
bu
s.

co

www.allsyllabus.com

Page 113
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 114
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-5&6
1. Derive the necessary condition for optimal operation of thermal
power plants with the transmission losses considered.
2. What are B- coefficients? Derive the matrix form of transmission loss
equation.
3. Explain the method of equal incremental cost for the economic
operation of generators with transmission loss considered.

co

4. Explain problem formation and solution procedure of optimal

scheduling for hydro thermal plants.

.a
lls

yl
la

bu

s.

5. What is the basic criterion for economical division of load between

units within a plant?

Page 115
vtu.allsyllabus.com

www.allsyllabus.com

UNIT-7 & 8
TRANSIENT STABILITY STUDIES
7.1 INTRODUCTION

bu

s.

co

Power system stability of modern large inter-connected systems is a major problem for
secure operation of the system. Recent major black-outs across the globe caused by system
instability, even in very sophisticated and secure systems, illustrate the problems facing
secure operation of power systems. Earlier, stability was defined as the ability of a system to
return to normal or stable operation after having been subjected to some form of disturbance.
This fundamentally refers to the ability of the system to remain in synchronism. However,
modern power systems operate under complex interconnections, controls and extremely
stressed conditions. Further, with increased automation and use of electronic equipment, the
quality of power has gained utmost importance, shifting focus on to concepts of voltage
stability, frequency stability,
inter-area oscillations etc.

.a
lls

yl
la

The IEEE/CIGRE Joint Task Force on stability terms and conditions have proposed the
following definition in 2004: Power System stability is the ability of an electric power
system, for a given initial operating condition, to regain a state of operating equilibrium
after being subjected to a physical disturbance, with most system variables bounded, so that
practically the entire system remains intact. The Power System is an extremely non-linear
and dynamic system, with operating parameters continuously varying. Stability is hence, a
function of the initial operating condition and the nature of the disturbance. Power systems
are continually subjected to small disturbances in the form of load changes. The system must
be in a position to be able to adjust to the changing conditions and operate satisfactorily. The
system must also withstand large disturbances, which may even cause structural changes due
to isolation of some faulted elements. A power system may be stable for a particular (large)
disturbance and unstable for another disturbance. It is impossible to design a system which
is stable under all disturbances. The power system is generally designed to be stable under
those disturbances which have a high degree of occurrence. The response to a disturbance is
extremely complex and involves practically all the equipment of the power system. For
example, a short circuit leading to a line isolation by circuit breakers will cause variations in
the power flows, network bus voltages and generators rotor speeds. The voltage variations
will actuate the voltage regulators in the system and generator speed variations will actuate
the prime mover governors; voltage and frequency variations will affect the system loads. In
stable systems, practically all generators and loads remain connected, even though parts of
the system may be isolated to preserve bulk operations. On the other hand, an unstable
system condition could lead to cascading outages and a shutdown of a major portion of the
power system.

Page 116
vtu.allsyllabus.com

www.allsyllabus.com

ROTOR ANGLE STABILITY

co

Small single (or small disturbance) rotor angle stability

Rotor angle stability refers to the ability of the synchronous machines of an interconnected
power system to remain in synchronism after being subjected to a disturbance. Instability
results in some generators accelerating (decelerating) and losing synchronism with other
generators. Rotor angle stability depends on the ability of each synchronous machine to
maintain equilibrium between electromagnetic torque and mechanical torque. Under steady
state, there is equilibrium between the input mechanical torque and output electromagnetic
torque of each generator, and its speed remains a constant. Under a disturbance, this
equilibrium is upset and the generators accelerate/decelerate according to the mechanics of a
rotating body. Rotor angle stability is further categorized as follows:

.a
lls

yl
la

bu

s.

It is the ability of the power system to maintain synchronism under small disturbances. In
this case, the system equation can be linearized around the initial operating point and the
stability depends only on the operating point and not on the disturbance. Instability may
result in
(i)
A non oscillatory or a periodic increase of rotor angle
(ii)
(ii) Increasing amplitude of rotor oscillations due to insufficient damping.
The first form of instability is largely eliminated by modern fast acting voltage regulators
and the second form of instability is more common. The time frame of small signal stability
is of the order of 10-20 seconds after a disturbance.
Large-signal rotor angle stability or transient stability

This refers to the ability of the power system to maintain synchronism under large
disturbances, such as short circuit, line outages etc. The system response involves large
excursions of the generator rotor angles. Transient stability depends on both the initial
operating point and the disturbance parameters like location, type, magnitude etc. Instability
is normally in the form of a periodic angular separation. The time frame of interest is 3-5
seconds after disturbance. The term dynamic stability was earlier used to denote the steadystate stability in the presence of automatic controls (especially excitation controls) as
opposed to manual controls. Since all generators are equipped with automatic controllers
today, dynamic stability has lost relevance and the Task Force has recommended against its
usage.
7.2 MECHANICS OF ROTATORY MOTION
Since a synchronous machine is a rotating body, the laws of mechanics of rotating bodies
are applicable to it. In rotation we first define the fundamental quantities. The angle _m is
defined, with respect to a circular arc with its center at the vertex of the angle, as the ratio of
the arc length s to radius r.

Page 117
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 118
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 119
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 120
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 121
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 122
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 123
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 124
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 125
vtu.allsyllabus.com

www.allsyllabus.com

This is the equivalent inertia constant on a base of 100 MVA and can be used when the two
machines swing coherently.
7.3 POWERANGLE EQUATION:
In solving the swing equation, certain assumptions are normally made (i) Mechanical power
nput Pm is a constant during the period of interest, immediately after the disturbance (ii)
Rotor speed changes are insignificant. (iii) Effect of voltage regulating loop during the
transient is neglected i.e the excitation is assumed to be a constant. As discussed in section
9.4, the powerangle relationship plays a vital role in the solution of the swing equation.
POWERANGLE EQUATION FOR A NONSALIENT POLE MACHINE:

.a
lls

yl
la

bu

s.

co

The simplest model for the synchronous generator is that of a constant voltage behind an
impedance. This model is called the classical model and can be used for cylindrical rotor
(nonsalient pole) machines. Practically all highspeed turbo alternators are of cylindrical
rotor construction, where the physical air gap around the periphery of the rotor is uniform.
This type of generator has approximately equal magnetic reluctance, regardless of the
angular position of the rotor, with respect to the armature mmf.

Page 126
vtu.allsyllabus.com

www.allsyllabus.com

7.4 POWER ANGLE EQUATION FOR A SALIENT POLE MACHINE:

.a
lls

yl
la

bu

s.

co

Here because of the salient poles, the reluctance of the magnetic circuit in which flows the
flux produced by an armature mmf in line with the quadrature axis is higher than that of the
magnetic circuit in which flows the flux produced by the armature mmf in line with the
direct axis. These two components of armature mmf are proportional to the corresponding
components of armature current. The component of armature current producing an mmf
acting in line with direct axis is called the direct component, Id. The component of armature
current producing an mmf acting in line with the quadrature axis is called the quadrature
axis component, Iq.

(9.42) gives the steady state power angle relationship for a salient pole machine. The second
term does not depend on the excitation and is called the reluctance power component. This
component makes the maximum power greater than in the classical model. However, the
angle at which the maximum power occurs is less than 90 o.
7.5 TRANSIENT STABILITY:
As defined earlier, transient stability is the ability of the system to remain stable under large
disturbances like short circuits, line outages, generation or load loss etc. The evaluation of
the transient stability is required offline for planning, design etc. and online for load
management, emergency control and security assessment. Transient stability analysis deals
with actual solution of the nonlinear differential equations describing the dynamics of the
machines and their controls and interfacing it with the algebraic equations describing the
Page 127
vtu.allsyllabus.com

www.allsyllabus.com

co

interconnections through the transmission network. Since the disturbance is large, linearized
analysis of the swing equation (which describes the rotor dynamics) is not possible. Further,
the fault may cause structural changes in the network, because of which the power angle
curve prior to fault, during the fault and post fault may be different. Due to these reasons, a
general stability criteria for transient stability cannot be established, as was done in the case
of steady state stability (namely PS > 0). Stability can be established, for a given fault, by
actual solution of the swing equation. The time taken for the fault to be cleared (by the
circuit breakers) is called the clearing time. If the fault is cleared fast enough, the probability
of the system remaining stable after the clearance is more. If the fault persists for a longer
time, likelihood of instability is increased. Critical clearing time is the maximum time
available for clearing the fault, before the system loses stability. Modern circuit breakers are
equipped with auto reclosure facility, wherein the breaker automatically recloses after two
sequential openings. If the fault still persists, the breakers open permanently. Since most
faults are transient, the first reclosure is in general successful. Hence, transient stability has
been greatly enhanced by auto closure breakers.

.a
lls

yl
la

bu

s.

Some common assumptions made during transient stability studies are as follows:
1. Transmission line and synchronous machine resistances are neglected. Since resistance
introduces a damping term in the swing equation, this gives pessimistic results.
2. Effect of damper windings is neglected which again gives pessimistic results.
3. Variations in rotor speed are neglected.
4. Mechanical input to the generator is assumed constant. The governor control loop is
neglected. This also leads to pessimistic results.
5. The generator is modeled as a constant voltage source behind a transient reactance,
neglecting the voltage regulator action.
6. Loads are modeled as constant admittances and absorbed into the bus admittance matrix.

The above assumptions, vastly simplify the equations. A digital computer program for
transient stability analysis can easily include more detailed generator models and effect of
controls, the discussion of which is beyond the scope of present treatment. Studies on the
transient stability of an SMIB system, can shed light on some important aspects of stability
of larger systems. The figure below shows an example of how the clearing time has an effect
on the swing curve of the machine.

Page 128
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

7.6 Modified Eulers method:

Eulers method is one of the easiest methods to program for solution of differential
equations using a digital computer . It uses the Taylors series expansion, discarding all
secondorder and higherorder terms. Modified Eulers algorithm uses the derivatives at the
beginning of a time step, to predict the values of the dependent variables at the end of the
step (t1 = t0 +t). Using the predicted values, the derivatives at the end of the interval are
computed. The average of the two derivatives is used in updating the variables.
Consider two simultaneous differential equations:

Page 129
vtu.allsyllabus.com

www.allsyllabus.com

Starting from initial values x0, y0, t0 at the beginning of a time step and a step size h we
solve as follows:

.a
lls

yl
la

bu

s.

co

Let

x 1 and y1 are used in the next iteration. To solve the swing equation by Modified Eulers
method, it is written as two first order differential equations:

Page 130
vtu.allsyllabus.com

www.allsyllabus.com

.a
lls

yl
la

bu

s.

co

Starting from an initial value _o, _o at the beginning of any time step, and choosing a step
size _t s, the equations to be solved in modified Eulers are as follows:

Page 131
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 132
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 133
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 134
vtu.allsyllabus.com

www.allsyllabus.com

.a
lls

yl
la

bu

s.

co

7.7

Page 135
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 136
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 137
vtu.allsyllabus.com

www.allsyllabus.com

.a
lls

yl
la

bu

s.

co

7.8

Page 138
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 139
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 140
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 141
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 142
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 143
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 144
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 145
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 146
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 147
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 148
vtu.allsyllabus.com

.a
lls

yl
la

bu

s.

co

www.allsyllabus.com

Page 149
vtu.allsyllabus.com

s.

co

www.allsyllabus.com

yl
la

bu

Generating units of present day have lower values of H, due to advanced cooling techniques,
which have made it possible to increase the rating of the machines without significant
increase in the size. Modern control schemes like generator excitation control, Turbine valve
control, single-pole operation of circuit breakers and fast-acting circuit breakers with auto
re-closure facility have helped in enhancing overall system stability.

.a
lls

Factors which can improve transient stability are

(i) Reduction of transfer reactance by using parallel lines.
(ii) Reducing transmission line reactance by reducing conductor spacing and increasing
conductor diameter, by using hollow cores.
(iii) Use of bundled conductors.
(iv) Series compensation of the transmission lines with series capacitors. This also increases
the steady state stability limit. However it can lead to problem of sub-synchronous
resonance.
(v) Since most faults are transient, fast acting circuit breakers with rapid re-closure facility
can aid stability.
(vi) The most common type of fault being the single-line-to-ground fault, selective single
pole opening and re-closing can improve stability.
(vii) Use of braking resistors at generator buses. During a fault, there is a sudden decrease in
electric power output of generator. A large resistor, connected at the generator bus, would
partially compensate for the load loss and help in decreasing the acceleration of the
generator. The braking resistors are switched during a fault through circuit breakers and
remain for a few cycles after fault is cleared till system voltage is restored.
(viii) Short circuit current limiters, which can be used to increase transfer impedance during
fault, there by reducing short circuit currents.
(ix) A recent method is fast valving of the turbine where in the mechanical power is lowered
quickly during the fault, and restored once fault is cleared.
Page 150
vtu.allsyllabus.com

www.allsyllabus.com

.a
lls

yl
la

bu

s.

co

UNIT-7 & 8
1. With the help of a flowchart, explain the method of finding the
transient stability of a given power system, based on Runge-Kutta
method.
2. Explain the representation of the following for power system
stability studies i) exciters ii) Governors
3. With the help of a flowchart, explain the method of finding the
transient stability of a given power system, based on modified
Eulers method.
4. Explain the solution of swing equation by point-by-point method.
Dec09/Jan10
5. With the help of block diagram explain the representation of
excitation control system and the speed governing system in
stability studies.
6. Explain clearly the representation of synchronous machine and load
for transient stability studies.
7. Derive the swing equation for a two machine system.
8. With the help of block diagram explain simplified representation of a
speed governor.
9. Write a node on i) Area control error ii)Automatic economic load
dispatch iii)Automatic voltage regulators

--------------------------------------------------*****------------------------------------------------------

Page 151
vtu.allsyllabus.com