Power System Fundamentals
Power System Fundamentals
Fundamentals
Power System
Fundamentals
Pedro Ponce • Arturo Molina • Omar Mata
Luis Ibarra • Brian MacCleery
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts have
been made to publish reliable data and information, but the author and publisher cannot assume responsibility
for the validity of all materials or the consequences of their use. The authors and publishers have attempted to
trace the copyright holders of all material reproduced in this publication and apologize to copyright holders if
permission to publish in this form has not been obtained. If any copyright material has not been acknowledged
please write and let us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmit-
ted, or utilized in any form by any electronic, mechanical, or other means, now known or hereafter invented,
including photocopying, microfilming, and recording, or in any information storage or retrieval system, without
written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.com
(http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Dan-
vers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and registration for a
variety of users. For organizations that have been granted a photocopy license by the CCC, a separate system of
payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used only
for identification and explanation without intent to infringe.
Preface xiii
Acknowledgments xv
vii
viii Contents
3 Hydroelectricity 123
3.1 Overall characteristics and operation . . . . . . . . . . . . . 124
3.1.1 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . 131
3.2 Integration with the infinite bus . . . . . . . . . . . . . . . . 135
3.2.1 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . 138
3.3 Theoretical problems . . . . . . . . . . . . . . . . . . . . . . 139
3.4 Homework problems . . . . . . . . . . . . . . . . . . . . . . . 142
3.5 Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
Bibliography 415
Index 421
Preface
Power systems are changing in the world because the conventional electric grid
systems cannot face ecological and technological demands from the end-user
of the electric grid. However, control systems, communication technologies,
digital systems, and power electronics are opening new electric grids called
smart grids. The Smart Grid has advanced engineering concepts that improve
the efficiency inside the smart grid but at present there are new challenges
that need to be solved; for instance, the real-time communication to meet
customer electrical demand in the smart grid system.
This book shows the fundamentals of power systems that are the basic
elements for understanding the smart grid technology, based on theoretical and
experimental challenges in which the conventional and advanced electric grids
have to be understood in undergraduate and graduate courses. Moreover, this
book also describes the technological knowledge that has to be understood by
the new generations of electrical engineers when they are dealing with power
systems and smart grid technologies.
Although conventional electrical energy is a well-known engineering area,
there are topics that have to be presented in a different manner because new
technologies have arisen such as power electronics and digital systems. On the
other hand, smart grid is a new concept for integrating several areas such as
renewable energy, transmission lines, conventional electric sources, communi-
cations, SCADA systems, and so on. Smart grid will bring new possibilities to
solve generation, distribution, and consumption problems of power systems.
On the other hand, one of the main challenges is to generate engineers who
understand power systems, so it is important to present in a clear and friendly
manner the principles that drive modern power systems. As a result, this book
illustrates step by step the topics that are needed to cover power systems, an
introduction to smart grid courses.
The content of the book moves from electric circuits to renewable electric
energy in order to connect the concepts that are needed; the book shows the-
oretical and experimental exercises so the reader can compare the simulation,
co-simulation, experimental, and theoretical results. For additional material,
visit https://www.crcpress.com/9781138554436
xiii
Acknowledgments
xv
1
Linear Electric Circuits
1
2 1 Linear Electric Circuits
In real problems, the initial voltage across the capacitor often cannot be
discerned. In this case it is convenient to set t0 = −∞ and v(−∞) = 0, so
that ˆ
1 t
v(t) = i(c) dt0 . (1.4)
C −∞
The last passive element is the inductor, and as the capacitor, its induc-
tance L is defined by a voltage-current equation:
di
v=L . (1.5)
dt
The unit in which inductance is measured is the henry [H]. The operating
principle of the inductor began with Michael Faraday, when he discovered that
a changing magnetic field could induce a voltage, proving that this voltage was
proportional to the time rate of change of current producing the magnetic field.
Equation (1.5) is the mathematical model of an ideal inductor. A physical
inductor can be constructed by winding a length of wire into a coil. Some
electrical characteristics can be obtained through this equation. First it is
shown that the voltage across an inductor is proportional to the time rate of
change of the current through it, therefore for a constant current there is no
voltage across the inductor, so it can be viewed as a “short circuit to DC.”
Another fact is that a sudden or discontinuous change in the current must
be associated with an infinite voltage across the inductor, but in reality an
infinite voltage does not exists.
The inductor current may be expressed in terms of voltage, solving Equa-
tion (1.5) for i(t):
1
di = vdt
L
ˆ
1 t
i(t) = vdt0 + t(t0 ). (1.6)
L −∞
1.2 Type of connection of passive elements 3
FIGURE 1.1
Series combination of N resistors
FIGURE 1.2
Series combination of N inductors
FIGURE 1.3
Series combination of N capacitors
In case the capacitors are connected in series (see Figure 1.3), the way to
obtain its equivalent value is not the same as for the previous passive elements.
The equivalent capacitor for a series connection is obtained as the inverse
sum of each capacitor inversed:
1 1 1 1
= + + ··· + . (1.10)
Ceq C1 C2 CN
FIGURE 1.4
Parallel combination of N resistors
1.3 Frequency-domain analysis 5
FIGURE 1.5
Parallel combination of N inductors
FIGURE 1.6
Parallel combination of N capacitors
the source, and substituting with ohm’s law, this can be obtained:
1 1 1 1
= + + ··· + . (1.11)
Req R1 R2 RN
For inductor elements, the equivalent circuit is obtained the same way as
for the resistors.
The equation of the equivalent inductor, represented in Figure 1.5, for a
parallel connection is
1 1 1 1
= + + ··· + . (1.12)
Leq L1 L2 LN
Finally, for N capacitors in a parallel combination as shown in Figure 1.6,
the equation of the equivalent capacitor is not the same as the last two, but
the algebraic sum of the capacitors in parallel.
The equation of the equivalent capacitor for a parallel connection is
Ceq = C1 + C2 + · · · + CN . (1.13)
v(t) = Re(Vmej(wt+θ))
Re
Vm
Vm
θ
to
Im t
to
-Vm
FIGURE 1.7
Frequency domain representation and time domain representation
voltage
v(t) = Vm cos (ωt + θ), (1.14)
where V m is the amplitude of the sinusoid, ω the angular frequency in [rad/s],
and θ is the phase in [rad].
Equation (1.14) can also be expressed as the real part of a complex quantity
accordingly to Euler’s identity
V = Vm ejθ . (1.16)
Since a complex number can also be represented in a polar form, the phasor
is generally represented as
V = Vm ∠θ. (1.17)
To differentiate a complex number from one that is not, it will be printed
as boldface type.
This process is called a phasor transformation, in which a function is taken
from a time-domain representation to a frequency-domain representation. It
is noted that the frequency-domain expression of a voltage or a current does
not implicitly include the frequency, as the time-domain expression does.
To obtain the sinusoid corresponding to a given phasor V , it has to be
multiplied by the time factor ejωt and just the real part is taken; see Figure 1.7.
The importance of the phasor in electrical engineering lies in the fact that
it is possible to define algebraic relationships between the voltage and current
for each of the three passive elements instead of using differential equations
[80]. The simplest case to begin with is the resistor, where
v(t) = Ri(t).
1.3 Frequency-domain analysis 7
The second passive element is the inductor. Its defining equation in a time-
domain is
di(t)
v(t) = L .
dt
It can be shown that the derivate of the voltage or the current is trans-
formed to the phasor domain as jωX, thus the complex expression for the
inductor is
Vm ej(ωt+θ) = jωLIm ej(ωt+θ) ,
and dividing by ejωt , the phasor relation for an inductor is found:
1.3.2 Impedances
An impedance can be defined as the ratio of the phasor voltage to the phasor
current that represents the opposition that the circuit exhibits to the flow
of sinusoidal current. Therefore an impedance can be treated, in a circuit
analysis, as a resistance with the exception that the impedance is a complex
quantity and its algebraic manipulation must be the one for a complex number
[27].
Previously, the phasor relation for the three passive elements, Equations
(1.5–1.7), were obtained,
I
V = RI, V = jωLI, V = .
jωC
Writing Equations (1.5–1.7) in terms of the ratio of the phasor voltage to
the phasor current, the impedance of each element is obtained:
V V V 1
= ZR = R, = ZL = jωL, = ZC = .
I I I jωC
8 1 Linear Electric Circuits
FIGURE 1.8
Phasor diagram for V1
V1
V1+V2
V2
FIGURE 1.9
Phasorial addition of two phasors
10 1 Linear Electric Circuits
Imaginary axis
w
V
Vm
Leading direction
Φ
Real axis
-θ
Lagging direction
Im
I
w
FIGURE 1.10
Leading or lagging direction of phasors
converter microprocessor
FIGURE 1.11
Phasor measurement unit (PMU)
a high sampling rate and a time synchronizing signal to determine the health
of the electricity distribution system. In Figure 1.11, PMU is represented in a
block diagram.
From the voltage and current measured with the PMU, typically at a rate
of 48 samples per cycle, the magnitudes and phase angles of the signals are
calculated in the phasor microprocessor. A phase-lock oscillator along with
the clock signal of the global positioning system (GPS) is used to provide
high-speed synchronized sampling with 1 microsecond accuracy for adding
a timestamp to every measured point. The measured phasors are known as
synchrophasors [44].
When data from multiple PMUs are combined together, the information
provides a global behavioral view of the interconnection. This made possible
some applications like the automation in power systems for smart grids, in-
1.3 Frequency-domain analysis 11
creasing the reliability of the power grid by detecting faults early, prevention
of power outages, increasing power quality, etc.
Calculate the current I for the three measured impedances if the voltage
is 226∠0◦ [Vrms ]. (V = IZ).
Feed the three modules individually with 220 [Vrms ] from the three-phase
supply unit and measure the current drawn with the maximum demand meter
(see Diagram 1.1). Write down your results.
Draw a phasor diagram and write R, C, L next to the corresponding phasor
current according to its nature and your previous calculations. Also write if
it is lagging, leading, or in phase with the voltage phasor.
The current phasor of a capacitive load leads its corresponding voltage
phasor while the current phasor of an inductive load lags the voltage phasor.
In this exercise the 90◦ shift of the capacitive load and the −90◦ shift of the
inductive load can be noticed.
When the resistive load phasors are in phase, it means that there is no
shift between the current phasor and voltage phasor.
12 1 Linear Electric Circuits
Since an AC power analysis is being made, let the voltage and current at
the terminals of the circuit be
v(t) = Vm cos(ωt + θv ),
i(t) = Im cos(ωt + θi ).
Diagram 1.1
Connections for different loads
FIGURE 1.12
Instant power function
14 1 Linear Electric Circuits
and depends on the phase difference between voltage and current. The second
term of equation (1.24) is a sinusoidal function that has a cyclic variation at
twice the applied frequency.
In Figure 1.12, it can be observed that p(t) has a period of T /2 due to
the second term and a vertical shift due to the first term of the instantaneous
power expression. Also noted is that the function is positive in some parts and
negative in others, which means that in the positive parts power is absorbed
by the circuit while in the negative parts power is absorbed by the source.
v(t) = Im cos(ωt + θv ),
The apparent power is the product of the effective values of the voltage
and current for sinusoidal waveforms. It is called apparent because of its re-
semblance to the average power if the voltage and current responses had been
DC quantities [27]. It is measured in volt-amperes [V A] to distinguish it from
other kinds of power already seen.
16 1 Linear Electric Circuits
The ratio of the average power to the apparent power is called the power
factor (PF), so it is the cosine of the phase difference between voltage and
current,
P
PF = = cos(θv − θi ). (1.31)
S
The PF may be seen as the factor that relates the apparent power to
the average power. Its importance lies in the effectiveness in which power is
transferred and therefore it affects the consumption cost of electricity. For
a low power factor (≈ 0) a load draws more current than for a high power
factor (≈ 1) with the same load, so the utility company may charge more to
a customer with a low power factor.
S = P + jQ, (1.34)
where P is the average power and Q is the new imaginary term, reactive
power. The dimensions of Q have to be the same for P and S, but in order to
avoid confusion, the unit of Q is defined as volt-ampere-reactive [VAr]. While
P is the actual power dissipated by the load, Q is a measure of the energy
exchange between the reactive components of the load and the source.
A physical interpretation can be made based on three cases for the reactive
power Q [27]; see Figure 1.13.
1. Q = 0 represents resistive loads (unity PF)
2. Q > 0 represents capacitive loads (leading PF)
3. Q < 0 represents inductive loads (lagging PF)
1.4 AC power analysis 17
Im
+Q Inductive loads
Re
-Q Capacitive loads
FIGURE 1.13
Power triangle with different reactive powers
It was said before that the complex power contains all the relevant power
information in a given load, which is shown in the next recap:
∗
Complex power = S = Vrms Irms = P + jQ = Vrms Irms ∠(θv − θi )
Average power = P = Re(s) = S cos(θv − θi )
Reactive power = Q = Im(S) = S sin(θv − θi )
p
Apparent power = S = |S| = P 2 + Q2 = Vrms Irms
P
Power factor = = cos(θv − θi )
S
θ1 V
V IL
IL
FIGURE 1.14
Circuit with an RL load
I
IC
θ2
V θ1 V
IL IC I
IL
FIGURE 1.15
Circuit with an RL load in parallel with a capacitor
If a capacitor C is added in parallel with the load into the circuit, a new
current IC is generated so that the total current drawn from the generator
becomes I = IC + IL . Again, making the phasor diagram, it can be seen that
the power factor is increased by the reduction of the phase angle between the
supplied voltage and the current from θ1 to θ2 ; see Figure 1.15.
Let us note that the magnitude of current I is less than the current IL with
the same supplied voltage. So, the reason to make a power factor correction
is to reduce the current drawn from the generators, representing savings for
both the power companies and the consumer.
To determine the value of the required capacitor for improving the PF of a
circuit, it is necessary to take into consideration that the supplied real power
has to remain the same.
From the power triangle for an inductive load with an apparent power S1
(see Figure 1.16),
Q2 = P tan(θ2 ).
Qc
S1
Q1
S2
Q2
θ1 θ2
FIGURE 1.16
Power triangle with Q reduction
determined as
Q P (tan(θ1 − θ2 ))
C= 2
= 2
, (1.35)
ωVrms ωVrms
where P = supplied real power, tan(θ1 ) = present power factor, tan(θ2 ) =
target power factor, ω = frequency of operation, Vrms = effective voltage.
1.4.7 Exercises
Most of the household appliances are a combination of inductive and resistive
loads. In this exercise a comparison between an inductive-resistive load and a
capacitive-resistive load will be made regarding their power consumption and
power factor. Finally, a correction in power factor for the inductive-resistive
load will be made to improve the power stats.
Recommended equipment and connections:
• Capacitive load
• Switchable capacitor battery
• LCR meter
20 1 Linear Electric Circuits
Connect the previously selected resistive load and capacitive load in series,
then feed them with 220 [Vrms ] from the three-phase supply unit. Measure
all the power parameters of the load with the maximum demand meter and
compare your results (see Diagram 1.2).
Now, for a series combination of the previously selected resistor and induc-
tor, calculate the current I of the circuit if it has a voltage of 226∠0◦ [Vrms ].
Then obtain the total complex power of the load and fill the next table:
Connect the previously selected resistive load and inductive load in series,
1.4 AC power analysis 21
Diagram 1.2
RC load connection
Diagram 1.3
RL load connection
then feed them with 220 [Vrms ] from the three-phase supply unit. Measure
all the power parameters of the load with the maximum demand meter and
compare your results (see Diagram 1.3).
For the RL load that is already connected, calculate the value of the ca-
pacitor that has to be added in parallel to improve the power factor to 0.97.
C = 1.36 [µF ].
Finally, connect the switchable capacitor battery module in parallel with
the RL load (see Diagram 1.4). Activate the first switch from left to right;
this will add a 2 µF capacitor to the circuit. Observe the new power factor
22 1 Linear Electric Circuits
Diagram 1.4
RL load with power factor correction
Vp Φ
Vp Φ
FIGURE 1.17
Single-phase three-wire system
where Van , Vbn , Vcn are respectively between lines a, b, and c, and the neu-
24 1 Linear Electric Circuits
a
Van Vbn Vcn
b
c c
b’ N a’
0
120° 240° wt
a S b
c’
FIGURE 1.18
A three-phase generator with generated voltages
a a
Vca Vab
Van
Vcn n b
Vbn
b Vbc
c c
FIGURE 1.19
Y-connected and ∆-connected source voltages
1.5 Polyphase circuits 25
a a
b ZY
ZY Z∆ Z∆
n
Z∆
ZY b
c
c
FIGURE 1.20
Y-connected and ∆-connected loads
tral line n known as phase voltages. Since the three voltages are 120◦ out of
phase and arbitrarily Van is chosen as the reference, there are two possible
combinations of sequence:
1. Van = Vp ∠0◦ , Vbn = Vp ∠ − 120◦ , Vcn = Vp ∠ − 240◦ ;
Za = Zb = Zc = ZY ; (1.38)
and similarly for a balanced ∆-connected load, the impedance per phase is
written Z∆ :
ZA = ZB = ZC = Z∆ . (1.39)
In this way any balanced load can be transformed from Y-connected to
∆-connected and vice versa using the next equations:
1
ZY = Z∆ , or Z∆ = 3ZY . (1.40)
3
26 1 Linear Electric Circuits
a IaA A
Van ZY
InN
Vcn n N
a Vbn
c b IbB B
ZY ZY
C
IcC
FIGURE 1.21
Balanced Y-Y circuit
Since the three-phase source and the three-phase load can be either Y-
connected or ∆-connected, four different connections can be made [27, p. 508]:
1. Y-Y connection (Y-connected source and Y-connected load)
2. Y-∆ connection (Y-connected source and ∆-connected load)
3. ∆-Y connection (∆-connected source and Y-connected load)
4. ∆-∆ connection (∆-connected source and ∆-connected load).
Vcn Vab
Vca
30°
Van
120°
Vbn
Vbc
FIGURE 1.22
Phasor diagram for a Y-Y circuit
√
Thus, it can be observed that the magnitude of the line voltages VL is 3
the magnitude of Vp : √
VL = 3Vp . (1.41)
It can also be noted that the line voltages are 120◦ out of phase between
each other and 30◦ out of phase between their corresponding phase voltages.
This is shown in the phasor diagram of Figure 1.22.
Now let’s apply KVL to each phase to obtain the line currents IaA , IbB ,
and IcC .
Van
IaA = ,
Zy
Vbn
IbB = = IaA ∠ − 120◦ ,
Zy
Vcn
IcC = = IaA ∠ − 240◦ .
Zy
Therefore,
InN = IaA + IbB + IcC = 0. (1.42)
Thus, no current flows through the neutral line if the source and the load
are both balanced. This means that as long as the system is balanced, no
matter the size of the impedance inserted in the neutral line, the neutral
current will remain zero.
In this kind of system, the magnitude
√ of the line voltage is different from
the phase voltage by a factor of 3 while the line current is the same as the
phase current.
28 1 Linear Electric Circuits
a IaA
A
Van
IAB
ICA
Vcn n Z∆ Z∆
a Vbn IBC
IbB
c Z∆
b IcC B C
FIGURE 1.23
Balanced Y-∆ circuit
It is noted by inspection that the line voltages are equal to the voltages
across the load impedances, making it possible to obtain the phase currents
as
Vab Vbc Vca
IAB = , IBC = , ICA = . (1.43)
Z∆ Z∆ Z∆
Applying KCL at nodes A, B, and C, the line currents are obtained from
the phase current previously calculated,
while the line current lags the corresponding phase currents by 30◦ , as shown
in the phasor diagram of Figure 1.24.
1.5 Polyphase circuits 29
IcC
Vcn
VCA VAB
ICA
IAB
Van
IaA
IBC
Vbn
IbB
VBC
FIGURE 1.24
Phasor diagram for a Y-∆ circuit
IaA
a A
c
IBC
b IbB
Vbc Z∆
IcC B C
FIGURE 1.25
Balanced ∆-∆ circuit
IaA
a A
ZY
Vca Vab N
c IbB B
ZY ZY
b
Vbc IcC C
FIGURE 1.26
Balanced ∆-Y circuit
1.5 Polyphase circuits 31
one way to obtain the line currents is to apply KVL to the loop that involves
IaA and IbB :
ZY IaA − ZY IbB = Vab .
Thus,
Vp ∠0◦
IaA − IbB = ,
Zy
Since it is a balanced system in positive sequence, the line current IbB
must be the same as the current IaA but 120◦ out of phase. Hence,
√
IaA − IbB = IaA 3∠30◦ .
Therefore, the other two line currents may be found by lagging 120◦ the
current IaA .
Another way to find the line currents is to either replace the ∆-connected
source with an equivalent Y-connected source or to replace the Y-connected
load with a ∆-connected load. Both cases were covered in the previous sec-
tions.
IaA A
VAN ZA
InN
N
VBN ZB ZC
IbB
B C
VCN
IcC
FIGURE 1.27
Unbalanced Y-connected loads
For a balanced system, it was said that the current in the neutral line
must be zero, hence it could be removed. In this case each of the line currents
produces a current though the neutral line, and applying KCL at node N, the
neutral line current is obtained:
InN = −(IaA + IbB + IcC ). (1.49)
To find the line currents for any of the other three-wire systems ∆-Y, Y-∆,
and ∆-∆, a mesh analysis has to be made, taking into account that √ phase and
line voltages and currents are no longer necessarily related by a 3 factor,
nor a difference in phase of 120◦ .
Sp
Sp Sp B
N OR = 3Sp
Sp
Sp Sp B C
C
B C
FIGURE 1.28
Power calculation for a balanced three-phase system
A A A
S1
S1 S2 B
N OR = S1+S2+S3
S3
S2 S3 B C
C
B C
FIGURE 1.29
Power calculation for an unbalanced three-phase system
In the case for an unbalanced system, the total complex power cannot be
obtained as three times the power of any of the loads because they are no
longer equal; rather, the total complex power is the sum of the three loads
(see Figure 1.29),
Stotal = S1 + S2 + S3 . (1.51)
It is important to notice that the relation of √
phase voltages to line voltages
and phase currents to line currents by the factor 3 is no longer valid, therefore
a mesh analysis is required to obtain the power of each phase.
An advantage of a three-phase system for power distributions is that uses
a lesser amount of wire than the single-phase system for the same line voltage
and the same absorbed power [27, p. 521].
1.5.9 Exercises
In some apartment buildings or condos, a three-phase motor is used to pump
water from the underground deposit to a water tank in the roof. This motor
can be connected in either a Y or ∆ connection. In this exercise a comparison
between these two connections can be made by measuring the power absorbed
by the three-phase load circuit.
Recommended equipment and connections:
34 1 Linear Electric Circuits
Connect a balanced three-phase Y-Y circuit (see Figure 1.30) where each
source’s phase is taken from the three phase supply unit and each load is an
RL combination of the previous values chosen. Note that each phase has to
have the same loads (see Diagram 1.5).
a IaA A
Van ZY
InN
Vcn n N
a Vbn
c b IbB B
ZY ZY
C
IcC
FIGURE 1.30
Three-phase Y-Y circuit electric diagram
Calculate the current IaA of the circuit if it has a voltage Van = 226∠0◦
[Vrms ]. Then obtain the power of the load ZY (S = V I ∗ ):
Diagram 1.5
Three-phase Y-Y circuit with RL load
This power corresponds to the absorbed power by the load of one of the
three phases. Since it is a balanced three-phase circuit, the total absorbed
power is obtained by multiplying the previous power by 3:
With the maximum demand meter, measure the real power absorbed for
each phase and the total power, and then compare your results.
SY STOTAL
Power (S) 16+j27 50+j80
Finally, modify the values of the loads for the second and third phase, now
the circuit is unbalanced. With the maximum demand meter, measure the
power being absorbed for each phase and the total.
Note that the total power absorbed now is obtained by the sum of each
phase power.
Connect a balanced three-phase Y-∆ circuit (see Figure 1.31) where each
source’s phase is taken from the three-phase supply unit and each load is an
RL combination of the previous values chosen. Note that each phase has to
have the same loads (see Diagram 1.6).
With the first maximum demand meter (the one on the left side in Dia-
gram 1.6), measure each of the line currents.
a IaA
A
Van
IAB
ICA
Vcn n Z∆ Z∆
a Vbn IBC
IbB
c Z∆
b IcC B C
FIGURE 1.31
Three-phase Y-∆ circuit electric diagram
Diagram 1.6
Three-phase Y-∆ circuit with RL load
Now, with the second module maximum demand meter (the one on the
right side of Diagram 1.6), measure each of the phase currents.
Note that in a balanced Y-∆ √ system, the phase current is related to its
corresponding line current by a 3 factor.
Finally, modify the values of the loads for the second and third phase
for the circuit to be unbalanced. Again, with the maximum demand meters,
measure the phase and line currents.
Note that in an unbalanced Y-∆ system, the√phase currents are no longer
related to their corresponding line current by a 3 factor.
1.6 Harmonics 37
1.6 Harmonics
So far it has been assumed that the voltage supplied to the loads is a perfectly
sinusoidal signal, but this is not necessarily true in reality. Mainly because of
the increasing use of nonlinear loads, the voltage and current waveforms suffer
from distortions, often expressed as harmonic distortion [50, p. 1].
Harmonics in AC power systems are components of a periodic current or
voltage waveform which have frequencies that are integer multiples of the
fundamental frequency of the system. A harmonic waveform is a perfectly
sinusoidal signal that is added to another sinusoidal waveform, which results
in a distorted waveform, as shown in Figure 1.32.
Figure 1.32 (left) shows the main waveform and its 3th, 5th, and 7th
harmonics, which can be expressed as
FIGURE 1.32
60 Hz signal and its 3th, 5th and 7th harmonics
Figure 1.32 (right) shows the resultant distorted waveform, which can be
38 1 Linear Electric Circuits
expressed as
vtotal = Vm1 sin (ωt) + Vm3 sin (3ωt + δ3 ) + Vm5 sin (5ωt + δ5 )
+ Vm7 sin (7ωt + δ7 ). (1.56)
The frequency of any harmonic can be easily found with the expression
V/3
V/5
V/7
1 2 3 4 5 6 7
Harmonic order
FIGURE 1.33
Harmonic spectrum
and is defined as the ratio between RMS values of signals, including harmonics
and signals considering only the fundamental frequency:
v
u∞
1 u X
THDF = t Fh2 , (1.64)
F1
h=2
P [kW]
pf = = 6= cos φ (1.65)
S [kVA]
p
S = P 2 + Q2 + H 2
q
[kV A] = [kW]2 + [kVAr]2 + [kVAr]2H
40 1 Linear Electric Circuits
S=kVA
H=kVARH
Q=kVAR
Φ P=kW
FIGURE 1.34
Power vector relation
3 90
Penalty = − 1 x 100 (1.66)
5 F.P.
1 90
Reward = 1− x 100 (1.67)
4 F.P.
1.6.2 Exercises
In a “healthy” electric installation, the AC current has a sinusoidal waveform,
but due to nonlinear loads, this waveform can be distorted by some harmonics
introduced in the signal. In this exercise, a sinusoidal waveform is going to
be observed in an oscilloscope with and without harmonics to appreciate the
distortion in the signal.
Recommended equipment and connections:
• Lamps for photovoltaic trainer
• Photovoltaic inclinable module
1.7 Theoretical problems 41
• Inverter grid
Diagram 1.7
Connection for harmonic distortion exercise
1.7 Theoretical problems 43
FIGURE 1.35
Current waveform in the oscilloscope
FIGURE 1.36
Current waveform with harmonics
FIGURE 1.37
Theoretical problem circuit
44 1 Linear Electric Circuits
FIGURE 1.38
Theoretical problem reduced circuit
Solution 1.1. Let’s begin reducing the circuit to one equivalent impedance
with a voltage source (see Figure 1.38). First note that the impedances of
4 [Ω] and –j6 [Ω] are in parallel with the impedances of 3 [Ω] and j4 [Ω]:
(4 − j6)(3 + j4)
Z1 = = 4.83 + j1.094Ω.
(4 − j6) + (3 + j4)
So, the equivalent impedance is zeq = 6.83 + j1.094. To obtain the current
Ix that corresponds to the source current, the source voltage is divided by the
equivalent impedance according to Ohm’s law:
(120∠10◦ )
Ix = = 17.34∠0.897◦ [A].
(6.83 + j1.094)
Problem 1.2. From the circuit of Problem 1.1, obtain the apparent, average,
and reactive power supplied to the impedances and the power factor of the
source.
Solution 1.2. The power supplied refers to the power delivered by the source.
To obtain the average power, Equation (1.26) is used:
But the easiest way to obtain the three parameters at the same time is by
getting the complex power. According to Equation (1.33):
∗ ◦ ◦ ◦
S = VRM S IRM S = (120∠10 )(17.34∠ − 0.897 ) = 2080.8∠9.103 [V A],
where apparent power is |S| = 2080.8 [V A], and average power is ReS =
2054.59 [W ], and reactive power is ImS = 329.20 [V Ar].
Now, the power factor can be easily obtained by the cosine of the angle of
the complex power or by Equation (1.31):
P
PF = = cos(θv − φi ) = 0.9874.
VRM S IRM S
1.7 Theoretical problems 45
Note that the average power and the power factors of the source are re-
quired. First, let’s obtain the complex power:
∗
220
S = 220 = 6643.13 + j5694.12 = 8749.52∠40.6◦ .
4.2 + j3.6
Hence, the power factor is PF = cos(40.6) = 0.7593.
The required capacitor connected in parallel to increase the power factor to
1 is:
6643.13(tan(40.6◦ − tan(0◦ )))
C= = 312µ[F ].
377(2202 )
Problem 1.4. A 4 wire three-phase Y-Y unbalanced system with negative
sequence has the following characteristics: Van = 220∠0◦ [V ], Z1 = 10+j8 [Ω],
Z2 = 8 + j6 [Ω], Z3 = 5 − j5 [Ω]. Obtain the three-phase current and the total
average power absorbed by the loads.
Solution 1.4. Since the sequence of the voltages is negative:
Van = 220∠0◦ ,
Vbn = 220∠120◦ ,
Vcn = 220∠240◦ .
Despite the system being unbalanced, it is easy to obtain the current of each
phase because of the fourth wire (the neutral one). The system can be analyzed
with three independent loops, each one having one of the source voltages and
its corresponding phase load:
220∠0◦
IaA =
10 + 8j
IaA = 17.18∠ − 38.66◦ [A];
220∠120◦
IbB =
8 + 6j
IbB = 22∠83.13◦ [A];
220∠240◦
IcC =
5 + 5j
IcC = 31.11∠ − 75◦ [A].
46 1 Linear Electric Circuits
To obtain the total average power, let’s get the complex power for all three
phases and then add just the real part of them:
∗
S1 = Van IaA
S1 = 2951.36 + j2361.10 [V A];
∗
S2 = Vbn IbB
S2 = 3872 + j2904 [V A];
∗
S3 = Vcn IcC
S3 = 4839.58 − j4839.58 [V A].
PT = 11, 662.94 [W ].
The first part of the function is known, so only the second part needs to be
computed. It can be written as
ˆ 2π
377 377
380 cos(754t + 0.75)dt,
2π 0
having
ˆ
1
cos(ωt + φ) = (cos(ωt) sin(φ) + cos(φ) sin(ωt)) , then,
ω
ˆ ˆ
sin(φ) sin(φ)
cos(ωt + φ)|t=0 = and cos(ωt + φ)|t=T = ,
ω ω
and implying
ˆ 2π
377 377
380 cos(754t + 0.75)dt = 0, so,
2π 0
ˆ
1 T 1
p(t)dt = Vm Im cos(θv − θi ) = 278.04 [W ].
T 0 2
Up to this point, we have already calculated the apparent power and active
power. The reactive power can then be obtained by any known method:
V (ω)
= i(ω),
jωL
and accordingly,
V (ω)
vωL
= ki(ω)k .
48 1 Linear Electric Circuits
FIGURE 1.39
Current plot
Simplifying:
Vrms
irms (ω) = ,
ωL
which can be plotted as shown in Figure 1.39.
As ω → 0, the inductance behaves as a short-circuit (its DC behavior).
Problem 1.7. Consider a capacitive load connected to a variable frequency
voltage source. If the capacitance magnitude is C and it is fed by Vm sin(ωt),
find the equation that describes ki(ω)k.
Solution 1.7. The capacitance can be described by
V 1
= ZC = .
i jωC
So, easily, the current can be derived from the above equation as
V (ω)jωC = i(ω),
and, accordingly,
kV (ω)jωCk = ki(ω)k .
Simplifying:
irms (ω) = Vrms ωC,
which can be plotted as shown in Figure 1.40. As ω → 0, the capacitance
behaves as an open-circuit (its DC behavior).
Problem 1.8. Power factor correction implies a current balancing between
the inductive part of the load and the added capacitor. If a load of ZL =
4.2 + j3.6 [Ω] is fed with 220 [Vrms ] at 60 [Hz], find the required capacitor to
obtain a PF= 1 without using Equation (1.22). Hint: see Figure 1.9.
Solution 1.8. We first need to find the current due to the inductive part of
the load:
220 [Vrms ]
= 30.20 − j25.88 [Arms ].
4.2 + j3.6[Ω]
1.8 Homework problems 49
FIGURE 1.40
Current plot
40
Current [A]
30
25.88
20
FIGURE 1.41
Inductance behavior
The inductance is of 9.55 mH, so around 377 [rad/s] its behavior is iL (ω)
(notice that the resistive part has been subtracted from the voltage as VR =
iR R, and iR is in phase with Vin ):
VL
220 − (30.2 × .2)
kıL (ω)k =
jωL
377 × 9.55 × 10−3 )
,
=
50 Inductor's current
C=312 µ F
C=250 µ F
C=370 µ F
40
Current [A]
30
25.88
20
FIGURE 1.42
Sweeping capacitance values
Problem 1.9. There are many buildings that take power from a three-phase
1.8 Homework problems 51
+
-
-
+
FIGURE 1.43 +
Electric circuit
1.9 Simulation
It is possible to analyze the effects of balanced or unbalanced loads in terms of
voltage and current (see Figure 1.44). Besides the fact that unbalanced loads
and harmonic effects will be a detriment to a grid’s performance, it is difficult
to find voltage variations due to the source’s stiffness. However, current is
directly affected, so it is more effective to measure current disturbances if
some faulty source or load is being evaluated.
Furthermore, an unbalanced load will account for a differentiated power
factor per phase. As the power factor is commonly measured as an average
value, it will not represent the grid’s reality for unbalanced networks. It is
important to acknowledge the possible variations regarding unbalanced loads
and harmonic distortion, and link their effects to common electric measure-
ments.
1.9 Simulation 53
FIGURE 1.44
Load unbalancing simulation
2
Power Flow and Electric Machinery Basics
As a first step into the electric grid, the fundamental topics about power
generation, transmission, distribution, and consumption must be established.
Each phase must then be first described and appropriately placed among
the others so the power flow is thoroughly understood. This chapter aims to
introduce the electric grid by presenting how the power flow occurs and by
offering the modeling basics around its major electromagnetic components.
An electric machine is a device that can convert either mechanical energy
to electrical energy or electrical energy to mechanical energy. When such a
device is used to convert mechanical energy to electrical energy, it is called
a generator. When it converts electrical energy to mechanical energy, it is
called a motor. Since any given electric machine can convert power in either
direction, any machine can be used as either a generator or a motor. Almost
all practical motors and generators convert energy through the interaction of
a magnetic field, and only machines using magnetic fields to perform such
conversions are considered in this book.
Many concepts must be established before performing any analysis of elec-
tric machines. The principle of electromechanical energy conversion is the
most important law of machine analysis. This theory allows us to establish
an expression of electromagnetic torque in terms of machine variables, like
the currents and the displacement of the mechanical system. In this chapter,
basic principles that will be mentioned are the equivalent circuit representa-
tions of magnetically coupled circuits, the concept of a sinusoidally distributed
winding, and the winding inductances, among others.
55
56 2 Power flow basics
water turbines, but diesel engines are sometimes used in remote locations.
Modern generators generate electrical power at voltages of 13.8 − 24 [kV ].
To construct an electric generator, the principle of “electromagnetic in-
duction” is used, as one of the fundamental elements of matter is precisely
the electromagnetic charge composed of a magnetic field and an associated
particle-motion electric field. A bosons generator uses a magnetic field to en-
ergize electrons kinetically and cause an interaction with other electrons, which
results in the generation of electric current and a voltage. Manipulating an
electromagnetic force may induce the displacement or movement of electrons,
and as a result, an electric current will flow.
As the rotor rotates at high speeds thanks to external mechanical energy
from the turbine, current flows occur in the stator copper wires (1 [A] is equal
to 6.25 × 1018 electrons moving through a wire per second). All power plants
have turbines and generators. Some turbines are powered by wind, water, and
steam from the earth or from the combustion of biomass, fossil fuels, and
other forms of energy. The electricity produced by a generator is transmitted
through the transmission lines, taking the power from plants to homes, schools,
and industry. Large-scale power generation requires huge power plants, com-
posed of many turbines.
A specific view of modern generation technologies is portrayed in some
chapters throughout this book: Hydroelectric generation is covered in Chap-
ter 3, Chapter 4 introduces wind power, and Chapter 5 deals with photo-
voltaic solar power generation. On the other hand, the basic tools to perform
electric analyses of transformers and electric machinery are provided in the
remainder of this chapter.
2.1.2 Transformers
Transformers convert AC electrical energy at one voltage magnitude into AC
electrical energy at another voltage magnitude. They are essential for the
operation of a modern power system, since transformers allow power to be
transmitted with minimal losses over long distances. Transformers are the
“connection” that holds the entire system together by increasing line voltages
and reducing currents for transmission over long distances, and then reducing
the voltages to the levels required by the end-users. They make modern power
systems possible, and furthermore, they do it with very high efficiency.
The invention of the transformer and the current development of AC power
sources forever eliminated range and power-level restrictions of power systems.
A transformer ideally changes one AC voltage level to another without affect-
ing the actual power supplied. If the transformer steps up the voltage level
of a circuit, it decreases the current due to energy conservation. Therefore,
AC electric power can be generated at one central location and stepped up
for transmission over long distances (exhibiting very low losses), to finally be
stepped down for its final use. Since the transmission losses in the lines of a
power system are proportional to the square of the current in the lines, raising
2.1 A glance into power flow 57
F = ma, (2.3)
where F is the net force applied to an object, m represents the mass of the
object, and a is the acceleration. A similar equation describes the relationship
between the torque applied to an object and its resulting angular acceleration.
This relationship is given by
τ = Jα, (2.4)
where τ is the net applied torque, α is the resulting angular acceleration, and
the term J serves the same purpose as an object’s mass in linear motion, and
is called the moment of inertia of the object.
2.2 Modeling basics 61
It is assumed that the force is collinear with the direction of motion. For
the special case of a constant force applied collinearly with the direction of
motion, this equation becomes
W = F dr. (2.6)
W = τ dθ. (2.8)
where B̂ is the magnetic field in Teslas [T], n̂ is the unit vector normal to the
surface, and da is a change in the area in [m2 ].
There are three cases for calculating the magnetic field. The first one is
when the magnetic field vector is uniform and perpendicular to S. Thus, the
magnetic flux is calculated by (2.13).
→
−
φB = B × A, (2.13)
where A is the area of the surface.
The second case is the magnetic field, uniform and at an angle to S ().
→
−
φB = B · n̂ × A. (2.14)
2.2 Modeling basics 63
Current - carrying Bar magnet
straight wire
I B N N
N
magnetic field lines
B S N N
FIGURE 2.1
Magnetic fields
The third case is the magnetic field that is non-uniform and at a variable
angle to S: ˛
→
−
B · n̂ da = 0 (2.15)
S
I I
I I S N
FIGURE 2.2
Lenz law defines the direction of the current
˛ ˆ
− →
→ − d →
−
E · dl = − B · n̂ da, (2.16)
C dt S
→
−
where the electric field is expressed by E that is a vector, and the increment
→
−
segment of the path C is dl .
The sign of the right-hand side of Faraday’s law is defined by the Lenz
law according to Figure 2.2. The EMF is produced when the magnetic flux is
changing, so a current is generated. This current generates a magnetic field
opposing the flux change that produced it. The Lenz law can be defined using
the energy conservation principle because the direction of the current opposes
the change. When the change appears, the current appears to limit it. Thus,
the current avoids an increment of the change that induces a bigger current
that increases the change.
An electrical charge in motion experiences a force when an electromagnetic
field is an area of influence of this electromagnetic field. The Lorentz equation
is defined by (2.17) in which the force direction and magnitude are described
by the cross-product [61].
F = qu × B, (2.17)
where F is the force [N], Q is the value of the electric charge [C], u is the
velocity of the charge [m/s], and B is the magnetic flux density [T].
If it is supposed that a current of magnitude I flows in a conductor of
length l within a uniform magnetic field, it is possible to rewrite Equation
(2.17) and get Equation (1.7)
F = lI × B. (2.18)
Motion
Field
Current
FIGURE 2.3
Fleming’s left-hand rule
which depicts Fleming’s left-hand rule [61] for determining the direction of
the force (motion).
Another basic law governing the generation of a magnetic field by some
current is Ampere’s law: ˛
H · dl = Inet , (2.19)
where H is the magnetic field intensity produced by the current Inet , and dl
is a differential element of length along the path of integration.
The magnetic field intensity H is, in a sense, a measure of the effort that
current is putting into the establishment of the magnetic field. The strength
of the magnetic field flux produced in the core also depends on the material
of the core. The relationship between the magnetic field intensity H and the
resulting magnetic flux density B produced within a material is given by
B = µH, (2.20)
FIGURE 2.4
Magnetizing B-H curve
where dA is the differential unit of area. If the flux density vector is perpen-
dicular to a plane of area A, and if the flux density is constant throughout
the area, then this equation reduces to
φ = BA. (2.24)
2.2 Modeling basics 67
Saturation
B (Ho, Bo)
a
2.0
Residual
Magnetism 1.5
1.0 b
0.5
c f
0 200 400 600 800 1000
H
d Coersive
(-Ho, -Bo) Force
FIGURE 2.5
Hysteresis
Thus, the total flux in the core in Figure 2.8 due to the current i in the
winding is
µN i
φ = BA = , (2.25)
lc
where A is the cross sectional area of the core.
An important effect called hysteresis is presented when the magnetic field
increases from positive to negative value. Hysteresis is characterized by the
residual magnetism and the coercive force that appear in the B-H curve (see
Figure 2.5). This magnetizing curve is determined according to the ferromag-
netic material.
To study magnetic circuits, analogies between magnetic and electric cir-
cuits can be created. Those analogies allow us to describe the magnetic circuit
and solve it using the same laws as for electric circuits. Table 2.1 shows a com-
parison of electric and magnetic circuits that was presented by [24]. In this
table, the reluctance is one of the most important quantities that is deter-
mined using the material permeability, the mean length, and the area. The
reluctance is an analog to the resistance.
Using the information of Table 2.1, an electric circuit can be created by
the magnetic circuit (see Figure 2.6). Electric circuits can represent complex
magnetic circuits, so the laws of current and voltages can be used for solving
such circuits.
Example 2.1. Consider Figure 2.7 to represent a quadrangular piece of ma-
terial with B = 1 [T], Lc =0.8 [m], Lg = 0.8 × 10−3 [m], µR =3000, and a
transversal area Ag = Ac = 16 × 10−4 [m2 ]. Calculate the current necessary
to hold B.
First we can calculate the reluctances of the metallic part of the circuit and
68 2 Power flow basics
TABLE 2.1
Analogies between electric and magnetic quantities [24]
Electrical Magnetic
Voltage V [V] mmf F [A-t]
Current I [A] Flux φ [Wb]
Resistance R [Ω] Reluctance R [H−1 ]
Electric field E [V/m] Magnetic field H [A-
intensity intensity t/m]
Current density J [A/m2 ] Flux density B [T]
Conductivity σ [S/m] Permeability µ [H/m]
Ohm’s law V = IR F = φR
Flow density J = σE B = µH
Admitance l
R = σA = 1
G R = µA
l
= 1
P
Electrical Magnetic
I
V R
FIGURE 2.6
Analogy between electric and magnetic circuits
FIGURE 2.7
Schematic for Example 2.1
2.2 Modeling basics 69
the air-gap:
Lg 0.8 × 10−3
Rg = = = 3.98 × 105 ,
µ0 Ag (4π × 10−7 ) (16 × 10−4 )
Lc 0.8
Rc = = = 1.33 × 105 .
µ0 µR Ac (4π × 10 ) (3000) (16 × 10−4 )
−7
and the mmf can be obtained as the flux was previously computed:
With the previous results, the current can be directly computed as:
F 84.88 [A − t]
I= = = 0.1415 [A]. (2.26)
N 600 [t]
FIGURE 2.8
Magnetically coupled circuit
The leakage flux Φl1 is produced by current flowing in coil 1 and it links
only the turns of coil 1. Likewise, the leakage flux Φl2 is produced by current
flowing in coil 2 and it links only the turns of coil 2. The magnetizing flux
Φm1 is produced by current flowing in coil 1 and it links all turns of coils 1
and 2. Similarly, the magnetizing flux Φm2 is produced by current flowing in
coil 2 and it also links all turns of coils 1 and 2.
With the selected positive direction of current flow and the way in which
the coils are wound, the magnetizing flux produced by positive current in one
coil adds to the magnetizing flux produced by positive current in the other
coil. In other words, if both currents are actually flowing in the same direction,
the magnetizing fluxes produced by each coil are in the same direction, making
the total magnetizing flux or the total core flux the sum of the instantaneous
magnitudes of the individual magnetizing fluxes. If the actual currents are in
opposite directions, the magnetizing fluxes are in opposite directions. In this
case, one coil is said to be magnetizing the core, the other demagnetizing.
All of the leakage flux may not link all the turns of the coil producing it.
Likewise, all of the magnetizing flux of one coil may not link all of the turns
of the other coil. The number of turns is considered an equivalent number
rather than the actual number. The inductances of the electric circuit resulting
from the magnetic coupling are generally determined from tests. The voltage
2.2 Modeling basics 71
FIGURE 2.9
Block diagram of an electromechanical system
λ1 = N1 Φ1 , (2.31)
λ2 = N2 Φ2 , (2.32)
magnetic fields while dielectric losses occur in all electric fields. Eddy currents
are electric currents produced by the changing magnetic field in the conductor
while hysteresis is the dependence of the output of a system on its current
input and past inputs.
If WE is the total energy supplied by the electrical source and WM the
total energy supplied by the mechanical source, then the energy distribution
could be expressed as
- - -
WE + We + + Wm + WM
Σ Σ Σ
- - -
WeS Wf WmS
Electrical System Coupling Field Mechanical System
FIGURE 2.10
Energy balance
• The energies stored in the electric or magnetic field that are not common
to both systems, denoted by WeS
• The energies stored in the mechanical system, shown as WmS . If the losses
of the coupling field are neglected, then the field is conservative and be-
comes
Wf = We + Wm . (2.38)
The examples of elementary electromechanical systems are shown in Figure
2.11, representing a magnetic coupling field, while the electromechanical sys-
tem shown in Figure 2.12, employs an electric field as a means of transferring
energy between the electrical and mechanical systems.
L r K
+ +
i ef M
v f
-
- fe
x0
FIGURE 2.11
Electromechanical system with magnetic field
L r M f
+ + fe
i
v ef
- -q +q
- D
x0
FIGURE 2.12
Electromechanical system with electric field
2.3 Machine parameters 75
a2’ a1’
a3’
bs axis a4’ b4
c1 b3
f1’ b2
c2 ωr
f2’
Θr
c3 f3’ b1
f4’ Φr
c4 f4
as axis
f3
f2 c4’
b1’
f1 c3’
b2’ c2’
b3’ c1’
b4’ a1 a2 a3 a4
cs axis fd axis
FIGURE 2.13
Elementary 2-pole, 3-phase, Y-connected salient-pole synchronous machine
It is assumed that each coil spans π radians of the stator for a 2-pole
machine. One side of the coil (coil side) is represented by ⊗, indicating that
the assumed positive direction of the current is down the length of the stator
(like if it were moving into this page). The point inside the circle, , indicates
that the assumed positive direction of the current is flowing from the inside
of this page. Each coil contains nc conductors, each of which makes up an
individual single conductor coil within the main coil. Thus the number of
turns of each winding is determined by the product of nc and the number of
coils’ sides carrying current in the same direction.
The winding configuration shown in Figure 2.13 is an oversimplification of
a practical machine. The coil sides of each phase winding are considered to
be distributed uniformly over 60◦ of the stator circumference. Generally, the
coil sides of each phase are distributed over a larger area, perhaps as much
as 120◦ , in which case it is necessary for some of the coil sides of two of the
phase windings to occupy the same slot.
In some cases, the coil sides may not be distributed uniformly over the
part of the circumference that the coil occupies. The windings span π radians
for the 2-pole machine. This is referred to as full-pitch winding. In order to
reduce voltage and current harmonics, the windings are often wound so that
they span slightly less than π radians for a 2-pole machine. This is referred to
as a fractional-pitch winding.
A salient-pole synchronous machine is selected for consideration since the
analysis of this type of machine may be easily modified to account for other
machine types. However, a salient-pole synchronous machine would seldom be
a 2-pole machine except in the case of small reluctance machines, which are
76 2 Power flow basics
FIGURE 2.14
A synchronous machine
of the synchronous class but do not have a field winding. Generally, 2- and 4-
pole machines are round-rotor machines with the field winding embedded in a
solid steel (nonlaminated) rotor. Salient-pole machines generally have a large
number of poles composed of laminated steel whereupon the field winding is
wound around the poles similar to that shown in Figure 2.13.
After directly measuring stator resistances and finding the coils’ distribu-
tion, the next step is to determine the self- and mutual inductances of the ma-
chine windings. It is advantageous to use the elementary 2-pole, 3-phase syn-
chronous machine to explain these inductance relationships. This procedure
can be readily modified to account for additional windings (damper windings)
placed on the rotor of the synchronous machine or for a synchronous machine
with a uniform air gap (round rotor). These inductance relationships are easily
altered to describe the winding inductances of an induction machine.
Figure 2.13 and Figure 2.15 are redrawn with the windings portrayed as
sinusoidally distributed windings. In a magnetically linear system, the self-
inductance of a winding is the ratio of the flux linked by a winding to the
current flowing in the winding with all other winding currents zero. A syn-
chronous machine is shown in Figure 2.14.
Mutual inductance is the ratio of flux linked by one winding due to current
flowing in a second winding with all other winding currents zero, including
the winding for which the flux linkages are being determined. For this analysis
2.3 Machine parameters 77
bs axis as’
bs
ar axis
cs fd’
ωr Φr Θr
as axis
fd
bs’ cs’
as
cs axis fd axis
FIGURE 2.15
Elementary 2-pole, 3phase, Y-connected salient-pole synchronous machine
FIGURE 2.16
Induction machine
bs axis as’
br axis
ar’ bs’
cs
br
Θr ωr ar a
cr
as a
Φr
cs’
br’
bs
ar cr’
cs axis as
cr axis
FIGURE 2.17
Two-pole, 3-phase, Y-connected symmetrical induction machine
2.3 Machine parameters 79
FIGURE 2.18
Wound rotor and squirrel cage rotor
in the short–circuited rotor windings. The frequency of the rotor currents cor-
responds to the difference in the speed of the rotating MMF due to the stator
currents and the speed of the rotor.
The induced rotor currents will in turn produce an air-gap MMF that
rotates relative to the rotor at a speed corresponding to the frequency of the
rotor currents. The speed of the rotor air-gap MMF superimposed upon the
rotor speed is the same speed as that of the air-gap MMF established by the
currents flowing in the stator windings. These two air-gap MMFs rotating
in unison may be thought of as two synchronously rotating sets of magnetic
poles. Torque is produced due to an interaction of these magnetic poles. It
is clear, however, that torque is not produced when the rotor is running in
synchronism with the air-gap MMF due to the stator currents since in this
case currents are not induced into the short-circuited rotor windings.
The winding inductance of the induction machine is expressed from the
inductance relationships given for the salient-pole synchronous machine. In
the case of the induction machine, the air gap is uniform. Thus, 2θr variations
in the self- and mutual inductances do not occur. These variations can be
eliminated by setting α2 = 0 in the inductance relationship given for a salient-
pole synchronous machine.
An induction machine can be a squirrel cage or wound rotor, as shown in
Figure 2.18.
The squirrel cage is the most-employed electric machine due to its robust-
ness, high efficiency, and low cost, while the wound rotor allows access to the
rotor windings to modify its torque and speed curve, or even feed the electric
machine to the rotor in order to increase its speed controllability. The disad-
vantage is that its maintenance is high due to the weathering of the brushes;
80 2 Power flow basics
FIGURE 2.19
Frontal view slip-rings induction machine
FIGURE 2.20
Induction machine assembly
IA
P
1 A
Three-phase
V
power source IB
Load
voltage and A Test
frequency
variables IC
P
2 A
FIGURE 2.21
Load test circuit
IA
a
A P
Three-phase
V
power source IB
b Locked
voltage and A Rotor
andjustable
frequency IC
c
A P
FIGURE 2.22
Locked rotor test circuit
82 2 Power flow basics
FIGURE 2.23
SCADA representation of the water governor to increase and decrease gener-
ator speed
FIGURE 2.24
SCADA representation of the switch that connects the generator to the loads
Close the correct bus bar R1 to connect the load R-L. The load R-L has
to be set to the position 0.
Increment the excitation current to reach the nominal armature voltage
Us .
Change the resistive load in a balanced way and observe the variations of
VL12, V23, and VL31 due to the variation of load. Increment the excitation
current Ie to level the output voltage Us .
2.3 Machine parameters 83
FIGURE 2.25
Excitation voltage
FIGURE 2.26
Loads
Change the inductive load in a balanced way and observe the variations of
VL12, V23, and VL31 due to load variation. Increment the excitation current
Ie to level the output voltage Us .
Register the data and draw a graph relating the excitation current Ie to
the variation of a purely resistive load and inductive load.
Table 2.2 has been plotted in Graphs 2.1 and 2.2.
Suggested material:
• Brushless driver for generating mechanical energy
• Power circuit breaker
• Feeder manager relay
• Module for measuring the electric power
• Three-phase asynchronous generator
• Programmable automatic power supply unit
• RS485 communication
• Resistive load
• Inductive load
84 2 Power flow basics
Graph 2.1
Resistive load with constant Us
Graph 2.2
Inductive load with constant Us
Graph 2.3
Resistive load with constant Ie
2.3 Machine parameters 85
Diagram 2.1
Suggested connection to load test
86 2 Power flow basics
TABLE 2.2
Data registered with U s[VVRMS ] constant, having Ie as the excitation current,
Iarm as the armature current, and US as the supplied voltage
Resistive Load Inductive Load
Balanced Ie [mA] Iarm [mA] U s[VRM S ] Ie [mA] Iarm [mA] U s[VRM S ]
Off 364 0 380 x x 380
R1 390 202 380 442 150 380
R2 403 285 380 478 218 380
R3 450 490 380 550 366 380
R4 516 704 380 640 536 380
R5 615 990 380 745 740 380
TABLE 2.3
Data registered with Ie [mA] constant, having Ie as the excitation current,
Iarm as the armature current, and US as the supplied voltage
Resistive Load Inductive Load
Balanced Iarm [mA] U s[VRM S ] Ie [mA] Iarm [mA] U s[VRM S ] Ie [mA]
Off 0 380 365 0 0 365
R1 192 360 365 130 326 365
R2 262 350 365 176 307 365
R3 410 318 365 260 269 365
R4 522 281 365 326 233 365
R5 618 236 365 390 206 365
R6 692 186 365 458 176 365
R7 722 161 365 492 156 365
Graph 2.4
Inductive load with constant Ie
FIGURE 2.27
SCADA representation of the water governor to increase and decrease gener-
ator speed
FIGURE 2.28
Excitation of the generator
Vexc ). Modify the excitation voltage and observe the armature voltage varia-
tion. Register the data in a table and draw the graph.
Sugested Material:
• Brushless control
• Power circuit breaker
88 2 Power flow basics
TABLE 2.4
Data registered with three different speeds, having Us , Ut , and Uw as the
voltage of each of the three phases
Speed (min-1) 1800 1000 700
Ie(mA) (Us,Ut,Uw) Vrms (Us,Ut,Uw) Vrms (Us,Ut,Uw) Vrms
100 117 63 46
150 172 95 69
200 225 124 88
250 275 154 107
300 323 181 126
350 366 202.6 143
400 408 225 158
450 440 243 170
500 467 259 181
550 487 274 192
Graph 2.5
Graph from Table 2.4
• Inductive load
2.3 Machine parameters 89
Diagram 2.2
Suggested connection to load test
90 2 Power flow basics
φm = φ1 + φ2
i1 φl 1 i2
φl 2
FIGURE 2.29
Transformer schematic view
Φ1 =Φl1 + Φm
Φ2 =Φl2 + Φm ,
where Φli stands for the leakage flux of winding i. It is called leakage flux as
it can be considered as a loss, because it is not used for core magnetization.
In this way, the linked flux λi of the ith winding is
λi = Ni Φi = Ni (Φl1 + Φm ),
The linkage flux can be given in terms of system’s inductances. The sub-
scripts related to each inductance Lij stand for the relation such inductance
keeps towards each coil. For instance, L11 focuses on those fluxes going ex-
clusively through winding 1, and is called self-inductance. On the other hand,
any Lij , i 6= j takes into account the flux going through two different windings
and thus, it is called mutual inductance. Our two-windings transformer can
be described by two equations as
λ1 =L11 i1 + L12 i2
λ2 =L21 i1 + L22 i2 .
Physically, the self-inductance of both windings is due to the leakage and
the actual magnetizing component Lmi (only due to current ii ). In order to
calculate both inductances, the opposite current is considered to be zero, so
λi |ij =0 N1 (Φli + Φmi )
L11 = = = N12 Pli + N12 Pm ) .
ii ii | {z } | {z }
Ll1 Llm
2 !
(1)
(1) di d N1 N1
e1 = Ll2 2 + Lm2 i1 i2
dt dt N2 N2
(1) 2
(1) di N1 d N2
= Ll2 2 + Lm2 i1 i2
dt N2 dt
N1
| {z } | {z }
Lm1 (1)
i2
(1)
(1) di2 d (1)
= Ll2 + Lm1 i1 i2 .
dt dt
The induced voltages are not the same as the terminal ones because the
windings themselves exhibit some resistance. Thus, the terminal voltage is the
algebraic sum of the induced voltage plus the voltage drop at the winding.
(i)
dii d i i + ij
vi = ii ri + ei = ii ri + Lli + Lmi .
dt dt
So, having the terminal voltage of winding i referred to winding j would
be
(j)
(j) d i + i
(j) di
j i
(j) (j) (j)
vi = ii ri + Lli i + Lmj .
dt dt
If both windings are seen from the primary side, the equations describing
the behavior of the transformer are
(1)
di1 d i1 + i 2
v1 = i1 r1 + Ll1 + Lm1 ,
dt dt
(1)
di2 d i1 + i2
v2 = i2 r2 + Ll2 + Lm1 .
dt dt
It is quite clear that the right-hand side of both equations is the same, phys-
ically implying that both voltages are influenced by a common inductance. As
2.4 Mathematical models of electric machinery 93
(1) (1)
e1 em e2 v2
r1 Ll 1 (1)
r2
(1)
Ll 2
N1 i2
v1 i1 (1)
i2 N2 v2
Lm1 Ideal
Transformer
FIGURE 2.30
Transfomer equivalent circuit
(1)
a result, it can be seen that the magnetization flux is due to a current i1 + i2 .
Both primary and secondary voltages share a common path to ground, pro-
vided by Lm1 . This voltage drop is later followed by two branches with a single
series resistive-inductive load, respectively. The preceding description points
to the equivalent circuit shown in Figure 2.30.
It is important to notice that the magnetization inductance is not con-
stant as the current rises. As the magnetic core gets saturated, the inductance
value will decrease as permeability does. It is also important to notice that
the equivalent circuit does not consider a load connected to the secondary
terminals. The output voltage can be associated to an electric load as it can
be defined as an impedance, so a current can be calculated and used to modify
transformer’s operating conditions. On the other hand, a full reformulation
of the transformer model can be performed for generic impedance consider-
ations as secondary loads; however, this process could lead to unnecessary
complexity.
It is noteworthy that a full grid simulation conformed by many subsystems
will impose unconsidered variations. However, a full descriptive model of the
transformer could consider the electric source connected to the primary to have
an associated series impedance (i.e., a Thevenin equivalent), while its output
at the secondary would have a generic load. A different approach could lead to
a step-by-step simulation where values at each stage are needed to recalculate
the system’s performance.
A three-phase transformer would require new considerations to be taken as
the flux linked by each winding depends on the other two phases. However, if
all three phases are balanced, the same monophasic approach used so far can
be used to estimate threefold electric-magnetic relations by neglecting new
mutual inductances.
It is clear that the application of arbitrary voltages to the primary can be
modeled by the simple equivalent circuit shown in Figure 2.30; however, trans-
formers are mostly connected to well-defined three-phase distribution systems,
so no arbitrary voltages are applied and the electric relation between phases
is determined in a standard manner. For instance, a three-phase transformer
94 2 Power flow basics
A a
VA
rn
Magnetic To Load
N Coupling
VB VC
B C b c
FIGURE 2.31
Y-Y connected transformer
rA Ll A (A)
ra
(A)
Ll a
NA ia
A iA (A)
ia Na va
Lm 1 Ideal
rn Transformer
N
FIGURE 2.32
Equivalent circuit of Y-Y connected transformer for phase A-a
To Load
N Magnetic
Coupling
C B
VB VC
iB b c
iC
FIGURE 2.33
∆-Y transformer connection
iCA
rA Ll A ra
(A) (A)
Ll a
VA NA ia
N iAB (A)
ia Na va
VB Lm 1 Ideal
Transformer
iBC N
FIGURE 2.34
Equivalent circuit of ∆-Y connected transformer for phase A-a
are classified depending on the number of voltage lines, rated power, NEMA
design, and speed (number of poles and operating frequency).
An induction motor can be understood as a transformer with rotating sec-
ondary; if the rotor is locked, the motor behaves as a short-circuited power
transformer. This implies that the transformer modeling presented in Sec-
tion 2.4.1 is useful to derive a steady-state model of the motor.
This machine can be represented by an equivalent electrical circuit to de-
scribe its behavior (see Figure 2.35). This circuit has a variable load dependent
on rotor’s slip. The sliding of an engine is the basis of an asynchronous machine
operation,
nsync − nm
s= , (2.41)
nsync
120fe
nsync = [rpm], (2.42)
p
where nsync is the maximum engine’s speed-dependent electrical frequency
and number of poles, and nm is the actual speed of the rotor.
The equivalent circuit of induction motor (shown in Figure 2.36) in Fig-
ure 2.35 can be reduced to the Thevenin equivalent circuit in Figure 2.37,
neglecting some electric interactions, where
Xm
VT H ≈ Vin , (2.43)
Xs + Xm
96 2 Power flow basics
FIGURE 2.35
Equivalent circuit of induction motor
Fs Fs
t
Va Vb Vc
Va Va
c+
Vb θr
Vb b- a-
Fs
Fr Fs Fr
a+ b+
Fsc
Fsb c-
Vc Vc
(b) Induced rotor’s magnetomotive force
FIGURE 2.36
Graphical representation of asynchronous motor auto-induction
2.4 Mathematical models of electric machinery 97
FIGURE 2.37
Thevenin equivalent circuit
2
Xm
RT H ≈ Rs , (2.44)
Xs + Xm
XT H ≈ Xs . (2.45)
Then, the torque motor with sliding is calculated as
3×VT2H ×R2
s
τind = . (2.46)
sωsync × [(RT H + Rs2 )2 + (XT H + X2 )2 ]
By knowing the machine’s parameters, it is possible to graph its speed-
torque relation, as shown in Figure 2.38. The most important point of the
graph is the maximum torque, because as its name indicates, the engine will
not be able to exceed this value.
The National Electrical Manufacturers Association (NEMA) has standard-
ized the design of electric motors to establish certain characteristics and re-
quirements that must be met by manufacturers. The main characteristics of
the engines as designed are:
Design A The standard design has a current of moderate start, and an
equally medium torque, and small slip in the range of nominal operation.
The maximum permissible slip is 5% and must be less than that of type B
design. The maximum torque is 200 to 300% nominal torque and slippage
occurs at less than 20%, while the starting torque for this engine is rated
at 60% for large motors and 200% or more for small engines. The main
drawback for this type of engine is the large amount of current required
to start the line (usually 5 − 7 times the rated current), and if there is
sensitive equipment connected to the same equipment, it can be damaged
by voltage drop. Very large motors should be started at reduced voltage to
mitigate the problem of voltage drop. Typical applications for this engine
are: pumps, machine tools, blowers, compressors, fans.
98 2 Power flow basics
Tem
max Tem
starting torque
− ωs 0 ωs ωr
2 1 smax 0 −1 s
FIGURE 2.38
Induction motor operating modes
Design B This type of engine has a medium starting torque, but with a lower
starting current (about 25% less than design A), little slip in the operating
range, and minimum torque of 200% nominal torque, although this is less
than design A because it has a higher resistance rotor (this also reduces
its efficiency). The applications are very similar to design A. Design B has
shifted to this design because of its lower starting current.
Design C Design C is for high starting torques (over 250% nominal torque)
and low starting current, slippage less than 5% in nominal operation, and
maximum torque is slightly lower than design A. As shown in Figure 2.39,
the pattern C has a design performance between A and B. The disad-
vantage of this engine is cost, because achieving this behavior requires
manufacturing of dual cage rotors. These engines respond very well for
jobs that require motor starting loaded, loaded as pumps, compressors,
and conveyors.
Design D These motors have a very large starting torque; some motors even
have maximum torque at start-up (275% over the nominal torque) and a
small current boot, but they have a huge slide in nominal operation (7 to
11 − 17%) and a very poor efficiency (because it requires a lot of resistance
in the rotor to get the boot with little current). These motors are used
to start moving large inertial mass; for example, panels that open the
floodgates of dams which require a lot of torque to reach full speed. After
performing a task, these engines are left off until their next operation is
required.
2.4 Mathematical models of electric machinery 99
450
300
Design D
Design C
150
Design B
0 10 20 30 40 50 60 70 80 90 100
Speed (% of synchronous speed)
FIGURE 2.39
Torque-speed curves of different NEMA standard motors
TABLE 2.5
Rotor stopped Rotor stopped
NEMA Code [kVA/hp] NEMA Code [kVA/hp]
A 0-3.15 L 9.00-10.00
B 3.15-3.55 M 10.00-11.00
C 3.55-4.00 N 11.20-12.50
D 4.00-4.50 P 12.50-14.00
E 4.50-5.00 R 14.00-16.00
F 5.00-5.60 S 16.00-18.00
G 5.60-6.30 T 18.00-20.00
H 6.30-7.10 U 20.00-22.40
J 7.20-8.00 V Más de 22.40
K 8.00-9.00
Given the above descriptions, the best design for the implementation of the
retrofit design is A, because it has the greatest efficiency in nominal operation
and weakness can be mitigated with start control.
To know more precisely what the starting current at the plate will be,
NEMA code determines which will be the starting power full motor voltage.
Table 2.5 shows the values of the constant that is multiplied by the engine
power design for the PTO:
With the design constant, we obtain the initial complex power S,
FIGURE 2.40
Brushless motor
N
Φr Rm Φl
Rs
S
Rl Φg
Φ
N 2Rg
Φr Rm
S Φr 2Rm 2KrRg
Rr
(a) (b)
FIGURE 2.41
a) Magnetic circuit model, b) magnetic circuit simplified
determine the magnetic flux linkage, then the generated electromotive force
(EMF), and finally the induced torque. The flow exits the north pole of the
rotor magnet, passes through the air gap by the stator, and travels through
the gap again to reach the south pole of a magnet. Finally it closes the cir-
cuit in the inner part of the rotor. A small part of the stream does not flow
through the stator, representing the leakage flow in parallel with reluctance
(Rl) coupled with the flow through the return gap. This is the reason why the
circuit of Figure 2.41(a) has 2Rg, the magnetic circuit, which is in series with
the stator. Figure 2.41(b) represents the first circuit model that describes the
direction of the magnetic flux (only one pole pair).
Permanent magnets are represented by sources with their respective reluc-
tance flow in parallel. The metallic inner part of the rotor, the air gap, and the
stator are represented as reluctances, and the elements that make the circuit
are: Rm, is the reluctance of the magnet, Rs, is the reluctance of the stator,
g, is the air gap reluctance, and Rr, is the reluctance of the metal part of the
rotor.
The airgap flux can be written in terms of the total flow Φg = kl Φr where
Kl is known as factor binding, usually close to 1. Thus, we can eliminate the
air-gap flux and simplify the magnetic circuit by removing Rl , then you can
add Rs , Rr , and 2RG , since they are in series and the magnets can be repre-
sented by one in parallel with a resistor of 2W D value, shown in Equation 2.50.
Considering also that the reluctance of the air gap is greater than the stator
and rotor, then we can simplify these two reluctances representing circuit (Rr
and Rs ) as factor Kr reluctances, multiplying the air gap reluctance Rg . Kr
is typically slightly larger than 1. With the simplified circuit Figure 2.41(b),
102 2 Power flow basics
λ = N Φg [W b − turns]. (2.54)
Equation 2.53 depends only on the angular position of the motor. In Fig-
ure 2.42(a), the magnetic flux of the magnet pole is opposite to that generated
by the stator winding sense, so that the flow is equal to λ = −N Φg . Moreover,
in Figure 2.42(b), the stator is 90 electrical degrees, so that the flow goes in
and out, resulting in a flow of λ = 0. Figure 2.42(c) shows that the flux of the
rotor magnet coincides with that generated by the stator winding current, so
that flow is λ = N Φg .
Considering that the flow varies linearly with respect to the angular posi-
tion, we get the graph shown in Figure 2.43; for the induced EMF in Equa-
tion (2.53), we apply Faraday’s law.
In Figure 2.43(a), the waveform of inductive coupling is shown when it is
pursuing a counter electromotive force at π. The electromotive force opposes
the passage of electrical current in an inductor, in this case the rotor, which will
tend to slow down. In Figure 2.43(b), the square waveform of the electromotive
force due to the triangular waveform applied to the inductive coupling is shown
at π. If we consider the law of conservation of energy,
eb i = T ωn , (2.55)
S
N
Rotor Rotor
S
θe = 0 N θe = π/2 N
N
S
S
Estator Estator
(a) (b)
S Rotor N
θe = π
N S
Estator
(c)
FIGURE 2.42
a) Motor with a fractional step, b) motor with 90◦ E, c) motor with 180◦ E
the energy necessary to re-align the magnetic domains, which are different
groups of magnetic currents within a magnetic field. A more realistic graph
of an engine is shown in Figure 2.44, where it can be seen that actually EMF
is trapezoid.
Traditional conductors have an electrical resistance that is proportional to
the length of the conductor and inversely proportional to the cross-sectional
area through which electric current flows. In the case of electric motors, the
windings contain a varnish layer that isolates the cables and prevents shorting
of windings in the stator. Thus, the resultant electric resistance is equal to
ρl
R= , (2.57)
A
where ρ is a constant known as electrical resistivity of the material, l is the
104 2 Power flow basics
λ eb
N Φg
Eb
π 2π θe
π/2 π 2π θe
-N Φg
(a) (b)
FIGURE 2.43
a) Waveform of the inductive coupling in π, b) square waveform of the elec-
tromotive force
λ
eb
N Φg
Θe Θe
π/2 π 2π π 2π
-N Φg
(a) (b)
FIGURE 2.44
a) Inductive coupling in its trapezoidal shape, b) electromotive force in its
trapezoidal shape
length of the conductor, and A is the cross sectional area of the conductor.
The electrical conductivity of the metals varies approximately linearly with
respect to temperature, thus
ρ(T ) = ρ(T0 ) [a + α(T − T0 )] , (2.58)
where α is the coefficient of thermal resistivity and T0 is the base temperature.
The magnetic flux is related to the current through a constant, modeled by
an inductance, so it should be included in the model brushless motor,
λ
L=
. (2.59)
i
If we consider the magnetic circuit of Figure 2.45, we can notice that Fr
contains a value point “potential” of 0, then the circuit is balanced. Moreover,
each flow Φ is equal to
Φ1 = Φ2 = Φ3 = Φ4 = Rg + Rm (2.60)
so that the total flow is
Nm Ni2
Φ1 + ΦN (2 + Φ3 + Φ4 ) = λ. (2.61)
Rg + Rm
2.4 Mathematical models of electric machinery 105
Ni
φ2
Rg
Rm
Rg Rm φ1 Ni
φ3
Fr
Ni Rm Rg
Rm
Rg
φ4
Ni
FIGURE 2.45
Magnetic circuit
Nm Ni2
L= λ. (2.62)
Rg + Rm
FIGURE 2.46
Digital control speed display
106 2 Power flow basics
Press the “Torque” and “Enable” buttons to enable the torque to the shaft.
Set the knob to the minimum and set the “Volt” switch to the ON position,
as shown in Figure 2.47
FIGURE 2.47
Analog control speed
Adjust the speed by varying the knob Volt. Check the actual speed on the
graphic. Reduce the speed to 0 [rpm] and switch OFF the level of the analog
control, as shown in Figure 2.47.
FIGURE 2.48
Software driven control speed
Press the combo box “control speed” and select SOFTWARE and then
press Apply to confirm (see Figure 2.48).
Press the “Control Type” combo box to select “Speed control”. Press the
“Set Ramp” combo box to select “Trapezoidal”. Set the acceleration/deceler-
ation profiles to 0, which implies a slow speed increase toward the speed set
point.
Increment the acceleration by using the buttons INC & DEC. Connect the
three-phase induction generator and check the mechanical parameters.
Measure the Mechanical power with the load to determine the power loss.
Type in the red block the speed at < 3000 and check the speed plot. Observe
the time delay of the actual speed with respect to the target speed. A similar
graph is shown in Figure 2.50.
2.5 Theoretical problems 107
FIGURE 2.49
Actual speed
FIGURE 2.50
Speed graph
eind = (v × B)I˙
= (vBsin90◦ )/cos30◦
= (10m/s)(0.5T )(1m)cos30◦
= 4.33 [V ].
108 2 Power flow basics
. . . . . -
.
. . . .
. . .
eind
v
.
i
. + . . .
vxB
. . . . . .
. . . . . .
FIGURE 2.51
Movement of a conductor
Problem 2.2. Figure 2.52 shows a coil of wire wrapped around an iron core.
If the flux in the core is given by the equation Φ = 0.05sin377t [W b], and if
there are 100 [turns] on the core, what voltage is produced at the terminals of
the coil? Of what polarity is the voltage during the time when flux is increasing
in the reference direction shown in the figure? Assume that all the magnetic
flux stays within the core (assume that the flux leakage is zero).
FIGURE 2.52
Wire wrapped around an iron core
2.5 Theoretical problems 109
Solution 2.2. The direction of the voltage while the flux is increasing in the
reference direction must be positive to negative. The magnitude of the voltage
is given by:
dΦ
eind = N
dt
d
= 377t0.05sin(100turns)
dt
= 1885cos377t,
or alternatively,
eind = 1885sin(377t + 90◦ ) [V ].
Problem 2.3. Find the relative permeability of the typical ferromagnetic ma-
terial at H = 50 [A − turns/m].
Solution 2.3. The permeability of a material is given by
µ = BH,
and the relative permeability is given by
µ
µr = .
µ0
Thus, it is easy to determine the permeability at any given magnetizing
intensity. At H = 50 [A − turn/m] and B = 0.25 [T ], so
0.25 [T ] H
µ = BH = = 0.005 ,
50 [A − turn/m] m
and H
µ 0.005 m
µr = = H = 3980.
µ0 4π × 10−7 m
Problem 2.4. A square magnetic core has a mean path length of 55 [cm] and a
cross-sectional area of 150 [cm2 ]. A 200 [turn] coil of wire is wrapped around
one leg of the core. The core is made of a material that has a magnetizing
intensity of H = 115 [A − turn/m]. How much current is required to produce
0.012 [W b] of flux in the core?
Solution 2.4. The required flux density in the core is
1.012 [W b]
B = ΦA = = 0.8 [T ].
0.015 [m2 ]
The magnetomotive force needed to produce this magnetizing intensity is
F = HLc = (115 [A − turn/m])(0.55 [m]) = 63.25 [A − turn].
So the required current is
[A − turn]
i = F/N = 63.25 = 0.316 [A].
200 [turn]
110 2 Power flow basics
Problem 2.5. From Problem 2.4, calculate the core’s relative permeability at
that current level and its reluctance.
Solution 2.5. The core’s permeability at this current is
0.8T H
µ = BH = = 0.00696 .
115A − turns/m m
Therefore, the relative permeability is
µ 0.00696H/m
µr = = = 5540.
µ0 4π × 10−7 H/m
The reluctance of the core is
63.25A − turns A − turns
F
R= = = 5270 .
Φ 0.012W b Wb
So, using the first and third equation together with short-circuit data:
Estimate the missing impedances and compare the results with those from
Problem 2.6, considering the transformer’s ratio is 1:2.
Solution 2.7. The open-circuit test can yield the value of X1 and Xm directly.
As the resistance is known, only the imaginary part of the equation provides
new information: Im( 120i1 ) = 710.059 [Ω] = X1 + Xm . The real part does
not match with the measured resistance. Such remainder is due to numeric
approximation issues.
For the short-circuit situation, the equivalent circuit can be solved through:
(1)
5 = i1 (R1 + jX1 ) + jXm (i1 − i2 )
(1) (1)
0 = i2 (R2 + jX2 ) + jXm (i2(1) − i1 )
X1 + Xm = 710.059
(1)
So, all the three inductances can be solved directly (remember to use R2
instead of R2 ) X1 = 659 [µΩ] + jω204.8 [µH], Xm = 659 [µΩ] + jω1.88 [H],
X2 (1) = 362 [µΩ] + jω204.9 [µH].
Regarding the first problem, there are some differences regarding the values
of the inductances. Notice that those differences are due to numerical preci-
sion, so conclusively, the common tests performed on transformers require high
measuring precision and calculation consistency.
Problem 2.8. Assuming a transformer with the following parameters, draw
the RMS value of the output voltage, from 1/3 to 3 times its nominal power
value.
Estimate the missing impedances and compare the results with those from
the first problem, considering the transformer’s ratio is 1:2.
112 2 Power flow basics
FIGURE 2.53
Voltage - Power curve
Solution 2.8. The circuit equivalent of the transformer can be used to find
its efficiency for some set of output loads. The equations that depict such
interaction are as follows:
(1)
V = i1 (R1 + jωL1 ) + jωLm (i1 − i2 ),
(1) (1) (1) (1) (1)
0 = i2 (R2 + RL + jωL2 ) + jωLm (i2 − i1 ),
(1)
where RL is the load resistance seen from the primary. By sweeping through
different resistance values, the following curve can be drawn:
Notice that there is a loss of about 5 [V ] when the transformer is consuming
1500 [W ]. Such loss is seen as a 10 [V ] loss if seen from the secondary winding.
Problem 2.9. Some self-made inductancehas a circular sectional area of
49 [mm2 ] and is made of iron (6.3 × 10−3 m H
), and a mean path length of
l = 5 [cm]. How many turns are required to generate a 3 [mH] inductance?
Φ
Solution 2.9. B = A , and L = NiΦ , so B = Li
NA .
In addition, F = N i = Hl, so
B Ll
H = N 2 A = µ, and finally
2 −3 −6 2
L(N ) = µNl A = (6.3×10 )(49×10
0.05
)N
= 6.17 × 10−6 N 2 .
So, for a 3 [mH] inductance, only 22 [turns] are needed.
Problem 2.10. How strong will the flux be if a current of 3 [A] is passed
through a 3 [mH] cylindrical
H inductance (r = 7 [mm], l = 5 [cm]), made with
electrical steel (5 × 10−3 m )? Do not consider saturation.
Φ NΦ Li
Solution 2.10. B = A ,L= i ∴B= NA .
F = N i = Hl ∴ H = Nl i .
B
µ= H = NLi Ll
AH = N 2 A . q
Li Ll lΦ2 µLi2 A
As N = Φ, µ= ( Li 2 = Li2 A ∴Φ= l .
q Φ ) A
(5×10−3 )(3×10−3 )(9)(154×10−6 )
Substituting: Φ = 0.05 = 644.8 × 10−6 [W b].
2.7 Simulation 113
FIGURE 2.54
Resulting hysteresis curve
2.7 Simulation
2.7.1 Transfomers programs
Problem 2.11. Hysteresis test
Plot the hysteresis phenomena of a monophasic transformer, with the fol-
lowing features: V = 120 [V ], V = 24 [V ], and I = 1.2 [A]. The measurement
of voltage and current of transformer are in the files “CurrentTransformerFil-
Pas.csv” and “VoltageTransformaderInteger.csv,” respectively. The MATLAB
program has been developed to get hysteresis phenomena in a plot and the code
is shown below.
Solution 2.11. After running the code, the resulting plot should be as shown
in Figures 2.54 and 2.55.
clc , clear ;
% Read Data
data_corriente = csvread ( ' C u r r e n t T r a n s f o r m e r F i l P a s .
csv ' , 9 , 0)
data_voltage = csvread ( ' V o l t a g e T r a n s f o r m a d e r I n t e g e r .
csv ' , 9 , 0)
114 2 Power flow basics
20
Voltage (V)
0
-20
-40
-0.02 -0.01 0 0.01 0.02
Time (s)
Current Input Transformer
-0.01
-0.015
Current (A)
-0.02
-0.025
-0.03
-0.02 -0.01 0 0.01 0.02
Time (s)
FIGURE 2.55
Resulting voltage and current
% Plot Hystesis
figure
ax1 = subplot (2 ,1 ,1) ;
ax2 = subplot (2 ,1 ,2) ;
% Voltage Output Transformer
plot ( ax1 , var1 , var2 )
title ( ax1 , ' Voltage Output Transformer ')
ylabel ( ax1 , ' Voltage ( V ) ')
xlabel ( ax1 , ' Time ( s ) ')
% Plot Results
figure
axy1 = subplot (2 ,2 ,[1 ,3]) ;
axy2 = subplot (2 ,2 ,2) ;
axy3 = subplot (2 ,2 ,4) ;
% Hysteresis Transformer
plot ( axy1 , var4 , var8 )
title ( axy1 , ' Hysteresis in Transformer ')
ylabel ( axy1 , ' Voltage ( V ) ')
xlabel ( axy1 , ' Current ( A ) ')
• The transformers operate at partial load, but this partial load is constant.
• Additional cost and inflation factors are not considered.
• Demand charges are based on 100% load
2.7 Simulation 117
The total cost of owning and operating the transformer for one year is thus
defined as follows:
• Capital cost (Cc ), taking into account the purchase price (Cp ), the interest
rate (p), and the depreciation period (n).
• Cost of no-load loss (CP 0 ), based on the no-load loss (P0 )and energy cost
(Ce ).
• Cost of load loss (CP k ), based on the load loss(Pk ), the equivalent annual
load factor (a) and energy cost (Ce ).
• Cost resulting from demand charges (Cd ), based on the amount set by the
utility and the total kW of connected load.
The next MATLAB program have been developed for analyzing a low-cost
transformer and a loss-optimized transformer.
clc ;
clear ;
% Transformer Loss Evaluation
% Distribution tranformer whith a low - cost
transformer
% Data :
% Depreciation period n = 20
years
% Interest rate p = 12 %
p.a.
% Energy charge Ce = 0.25
Euro / kWh
% Demand charge Cd = 350
Euro /( kW . year )
% Equivalent annual load factor alpha =
0.8
% Cc Capital Cost
% where :
% Cp Purcharse price
% r Depreciation factor
% q Interest factor
% p Interest rate
% n Depretation period
n = 20; % Years
p = 12; % 12% p . a .
Ce = 0.25; % Euro / kWh
Cd = 350; % Euro /( kW . year )
118 2 Power flow basics
q = p /(100) +1
r = ( p *( q ^ n ) ) /(( q ^ n ) -1)
clc ;
clear ;
% Transformer Loss Evaluation
% Distribution tranformer whith a loss - optimized
transformer
% Data :
% Depreciation period n = 20
years
% Interest rate p = 12 %
p.a.
2.7 Simulation 119
% Cc Capital Cost
% where :
% Cp Purcharse price
% r Depreciation factor
% q Interest factor
% p Interest rate
% n Depretation period
n = 20; % Years
p = 12; % 12% p . a .
Ce = 0.25; % Euro / kWh
Cd = 350; % Euro /( kW . year )
alpha = 0.8; % Equivalent annual load factor
Po = 16; % kW No - load loss
Pk = 124; % kW Load loss
Cp = 585000; % Purcharse price ($)
q = p /(100) +1
r = ( p *( q ^ n ) ) /(( q ^ n ) -1)
FIGURE 2.56
Power manipulation through a synchronous generator
3
Hydroelectricity
123
124 3 Hydroelectricity
Other alternative classifications are made about flow control, load condi-
tions, and head altitude [33, pp. 214–218].
Impoundment hydroelectric generators (Figure 3.1) use the most basic op-
eration principle, as they are based on converting potential energy of raised
accumulated water into a quick flow to move a generator’s propeller. There are
six fundamental components in this type of power plant: dam, reservoir, pen-
stock, turbine, generator and governor [32, p. 55]. Other secondary elements
are machine room, internal pipelines, and oscillation tanks [33, p. 183].
The dam is the concrete wall that holds the water in the reservoir and
allows potential energy to be accumulated. In order to move the turbine at-
tached to the generator, the governor (control valve) must allow the water
to flow through the penstock; once the water has completed its trajectory, it
is released to the outlet, which is an open water body. The head represents
the height difference between the reservoir surface and the water outlet. This
3.1 Overall characteristics and operation 125
Dam
Pe
ns
Reservoir
toc
Head Governor
k
Generator
Transformer
Outlet
Transmision lines
FIGURE 3.1
Impoundment hydroelectic power plant representation
TABLE 3.1
Example features of two large impoundment hydroelectric plants [32, p. 56]
Feature Grand Coulee Three Gorges
Length of dam [km] 1.6 2.34
Height of dam [m] 170 185
Width at river base [m] 150 115
Width at crest [m] 9 40
Volume of concrete [m3 ] 9.16 × 106 28 × 106
Surface area of reservoir 320 72, 128
behind dam [km2 ]
2000
1500
Water runoff [m3/s]
1000
500
0
J F M A M J J A S O N D
Months
FIGURE 3.2
Water runoff example graph
surpluses. The reservoir must at least hold 200 [m3 /s-month] for October and
300 [m3 /s-month] for November, which results in a total minimum capacity of
1.31 × 109 [m3 ]. However, power consumption is not constant and in this case
has been only estimated as a mean demand value; moreover, hydrographical
data is not the same every year and the reservoir must be prepared to hold a
bigger amount of water. Previous calculations were only demonstrative.
Considering the mean water load value of the previous example, a prelim-
inary nominal power value can be assigned to the impoundment hydroelectric
plant. Assume the plant’s head is 90 [m], with 20% in losses due to friction,
turbulence, and generator efficiency,
FIGURE 3.3
Water load through time example graph
the turbine axis. The momentum gained by each cup makes the runner turn
at constant speed if mechanical load equals the gained momentum through
time; regulation must be made by the governor valve in order to modify the
amount of water and its speed through turbine nozzles.
In order to prevent the water jet from arbitrarily splashing in the cup and
canceling the axial momentum of the runner, a different design, which consid-
ers two attached cups (Figure 3.5), can be used so the reflected water takes
a known path. The harnessed momentum equals mi (vi − vc ), where subscript
i represents the water jet and c, the cup. On the other hand, the remaining
momentum (lost) is mr (vr + vc ), where subscript r stands for reflected water.
Notice that both speeds are added as their direction is opposite; at last, the
total momentum acquired by the cup is the difference between the gained
momentum and the lost one. The mass of incoming water must be equal to
that of the reflected one as well as its speed [32, pp. 58–59].
Water
flow
Nozzle
Runner
Water flow
ets
uck
Volute Blades
Chamber
FIGURE 3.4
Turbine types schematics
Reflected Jet
Cup
Cup
Reflected Jet
FIGURE 3.5
Improved cup design and water jets indication
130 3 Hydroelectricity
TABLE 3.2
Turbine types summary
Type Max. Max. Power Max.
Head [m] [MW] Diameter [m]
Francis 2000 >250 5.5
Pelton 500 >750 10
Kaplan 70 >225 10
3.1.1 Exercises
Based on the nominal supply of 380 [VL-L ] 60 [Hz], achieve a steady no-
load operating point for the hydroelectric generator. Attain desired frequency
by modifying the governor position using the SCADA system and set the
output voltage to 380 [V ] manually (without an automatic controller). Use
Diagram 3.1 for connections.
FIGURE 3.6
Power flow of a hydroelectric plant
Diagram 3.1
Testing synchronous generator’s void operation
3.1 Overall characteristics and operation 133
• Resistive load
Make the hydraulic generator operate steadily at 380 [VL−L ] 60 [Hz] at full
load conditions (lowest resistive load value possible) and measure delivered
power. Use Diagram 3.2 for connections. If the didactic generator is scaled
1:1000 and represents a hydroelectric power plant, answer and justify the
type of plant it would be, the needed water flow at full-load conditions, and
the potential energy needed at the dam to attain an operation of 1 month (30
days) at maximum rating (estimate total losses to be of 20% and the head of
50 [m]).
Max. Power [W] × 1000 Type Water flow [m3 /s] Energy [GJ]
332 332 [kW ] Small 0.846 860
• Resistive load
• Maximum demand meter
Make the hydraulic generator operate steadily at 380 [VL−L ] 60 [Hz] at
average load conditions and measure delivered power. Use Diagram 3.3 for
connections. If the didactic generator is scaled 1:1000 and represents a hy-
droelectric power plant, answer and justify the amount of water needed at
the reservoir to operate 3 months (90 days) without water runoff, the type
of turbine needed for such a plant, and the average water flux. Considering
134 3 Hydroelectricity
Diagram 3.2
Synchronous generator with resistive loads
Diagram 3.3
Instrumented operation
3.2 Integration with the infinite bus 135
the average demand line, estimate the needed mean water load if a risk factor
of 20% is desired. Use a mean resistive value as the average consumption.
(Estimate total losses to be of 20% and the head of 50 [m].)
Without risk factor:
Transformer Transformer
Load Center 1 Load Center m
44 kV 44 kV
Transformer Transformer
transformer
transformer
Pole
Pole
4.16 kV 4.16 kV
Transformer Transformer
FIGURE 3.7
Infinite bus representation
Vcn Vbvn
120 ω 120 ω
° °
120°
120°
o o
Van Van
12 12
0° 0°
Vbn Vcn
FIGURE 3.8
Positive and negative phase sequences
Synch Lamps
Vb
Breaker
A B C
Prime Vg
Mover
Generator
FIGURE 3.9
Schematic of integration of a generator to the infinite bus
138 3 Hydroelectricity
erator parameters need to be adjusted: shaft speed, field current, and phase
sequence. The shaft speed shall be adjusted so the output voltage frequency
is 60/50 [Hz]; similarly, the field current must be set so both Vb and Vg read
the same; if the phase sequence or phase are not correct, the lamps will never
hold a brightness level [68, p. 301].
There are four possible scenarios considering Figure 3.9 as described by
Sen [68, pp. 301–303]. The following list assumes all other parameters are the
same; obviously, if all voltages are exactly the same all lamps would be totally
off.
1. Voltages are different
All lamps’ brightness will be steady and equal. Adjustments must be made
in the generator’s field current.
2. Frequencies are different
All lamps brighten and darken simultaneously (change speed depends on
how different the frequencies are). Shaft speed must be modified by using
the prime mover (Coaxial electric machine connected to the shaft of the
generator — Figure 3.9) of a different breaking/accelerating device. If shaft
speed is changed, the output voltage also will change, so adjustments in
the field current will also be necessary.
3. Phase sequences are different
Lamps will glow and fade independently. Phase connections must be in-
terchanged between two of them, e.g., A with B.
4. Phases are different
All lamps glow with the same intensity. Generator speed must be slightly
modified so the lamps slowly fade; the connection must be made when
all of them are totally dark. Once the generator is connected, its speed
will be held constant by the infinite bus; however, the delivered active and
reactive powers can be controlled by adjusting the prime mover power and
the field current.
As told before, lamps usage is suitable for laboratory test conditions. Mod-
ern development in electronics and digital systems has made it possible to pro-
duce automatic breakers and field control systems that enable the automatic
connection to the infinite bus.
3.2.1 Exercises
A synchronous generator is to be incorporated to the infinite bus. Achieve a
successful connection through the generator synchronizing relay and measure
the generator’s contribution to the grid under various R-L load conditions and
water flow. Use Diagram 3.4 for connections.
Start the generator and make its voltage through the power transformer
3.3 Theoretical problems 139
equal to the grid’s voltage, and attain a speed such that output frequency is
60 [Hz]. Varying the brushless servo speed in very little amounts will slowly
make the phase of both grid and generator match, and the automatic re-
lay will close the connection automatically. Be very careful with phase
sequence; make sure the rotation direction is consistent with the
grid’s sequence before starting.
L stands for line — G stands for generator
LOAD Speed PL vs PG QL vs QG SL vs SG
[RMS] [W] [VAr] [VA]
R2 1818 110.08 / 144.8 -556.1 / 632 566.88 / 648.3
R2 1845 92.2 / 157.2 -498.38 / 573.3 506.82 / 594.2
R4 1845 280.86 / 160.8 -350.56 / 552.4 449.64 / 575
..
.
Diagram 3.4
Hydroelectrical power plant synchronization connections
3.3 Theoretical problems 141
Graph 3.1
Water load vs time availability graph
J F M A M J J A S O N D
1200 1000 900 650 800 1200 1750 1750 900 750 650 1400
m3
day h min h s i
1000 × 2[month] 30 24 60 60 × 1.15 =
s month day h min
1.79 × 109 m3 .
Due to power capacity, this plant is large and needs a Pelton turbine due
to head restrictions. So results are:
Problem 3.2. Plot the water load vs time availability for the conditions de-
scribed below. Plot the mean water load and demand lines (550 [m3 /s]). All
values in [m3 /s]
J F M A M J J A S O N D
500 600 600 400 600 800 900 700 700 500 400 500
Problem 3.3. How tall does the head need to be so an average power of
300 [M W ] could be taken from a hydroelectrical plant, considering 350 [m3 /s]
mean water flow and an overall efficiency of 75%?
Solution 3.3.
300[M W ]
h i 3 = 116[m].
kg
× 350 ms × 9.81 sm2 × 0.75
1000 m3
Problem 3.4. Considering the data provided in Problem 3.2, what is the total
amount of deficiencies expected in a whole year?
Solution 3.4. Whole year = −450 [m3 /s] and the longest gap = −300 [m3 /s].
Problem 3.5. Considering the data provided in Problem 3.2, calculate the
reservoir capacity considering a demand line of 600 [m3 /s].
Solution 3.5. The longest gap deficiency is −500 [m3 /s], in a period of 4
months, so:
m3
day h min h s i
500 × 4 [month] × 30 24 60 60 =
s month day h min
5.18 × 109 m3 .
J F M A M J J A S O N D
700 750 900 800 900 1200 1300 1200 1000 800 700 650
Calculate the overall required efficiency of the plant if the average demand
is about 935 [M W ] and the head is 140 [m].
Solution 3.6. The mean can be easily computed to be 908.33 [m3 /s]. Con-
sidering a head of 140 [m], the power generation can be expressed as Wn =
δQghη, implying:
3
kg m h m i
935 M W = 100 908.33 9.81 (140 [m]) η.
m3 s s2
So, η = 0.749.
3.4 Homework problems 143
J F M A M J J A S O N D
700 750 900 800 900 1200 1300 1200 1000 800 700 650
620 500 440 490 580 720 800 700 650 600 600 710
If its head is 90 [m] and its overall efficiency is 0.8, build a table to show
how much water must be released during each month, considering the initial
reserved water to be 500 × 106 [m3 ], and the total dam capacity to be 1 ×
109 [m3 ]. Report whether a power shortage is possible.
Solution 3.8. It would be useful to calculate the available power per month,
considering each month separately. In addition, calculate the difference regard-
ing the available power and finally, the compensating runoff.
The generated power and equivalent runoff can be easily computed by fol-
lowing the equation W = δQghη.
144 3 Hydroelectricity
3.5 Simulation
Once a hydraulic power plant is connected to the infinite bus, the behavior
and characteristics of the generator will depend on clear parameters. As the
voltage is led by the bus because it dominates the generator (in terms of power
capacity), variations of input mechanical power (water power) or field voltage
will not exhibit a direct effect on output voltage. However, they will enable
the power flow to/from the infinite bus.
A synchronous machine can provide its own magnetization between stator
and rotor so it can provide reactive power to the network. As the relation
between active/reactive power is clearly dominated by field current and me-
chanical power, respectively, this generator can quickly manage both types of
grid variations. The diagram of such simulated system is shown in Figure 3.10.
146 3 Hydroelectricity
FIGURE 3.10
Generator interconnection to the infinite bus simulation
TABLE 3.3
3.5 Simulation
J F M A M J J A S O N D
Runoff 700.00 750.00 900.00 800.00 900.00 1200.00 1300.00 1200.00 1000.00 800.00 700.00 650.00
[m3 /s]
Gen Power 494.42 529.74 635.69 565.06 635.69 847.58 918.22 847.58 706.32 565.06 494.42 459.11
[M W ]
Req Power 620.00 500.00 440.00 490.00 580.00 720.00 800.00 700.00 650.00 600.00 600.00 710.00
[M W ]
Computations to solve Problem 3.8
OverPower −125.58 29.74 195.69 75.06 55.69 127.58 118.22 147.58 56.32 −34.94 −105.58 −250.89
[M W ]
OverRunoff −177.79 42.11 277.05 106.26 78.84 180.63 167.37 208.95 79.74 −49.47 −149.47 −355.21
[m3 /s]
Compensate −4.61 1.09 7.18 2.75 2.04 4.68 4.34 5.42 2.07 −1.28 −3.87 −9.21
×108 [m3 ]
Dam status 3.92 1.48 8.66 1.00 1.00 1.00 1.00 1.00 1.00 8.72 4.84 −4.36
×108 [m3 ]
Release 0.00 0.00 0.00 −1.42 −2.04 −4.68 −4.34 −5.42 −2.07 0.00 0.00 0.00
×108 [m3 ]
147
4
Wind Energy
Wind energy has been exploited for thousands of years. Its principal ancient
applications are extracting water from wells, cutting lumber, grinding grain
for making flour, and other mechanical power applications. Nowadays wind is
employed mostly for electric power generation. Wind generators are also called
wind turbines, because their function is similar to the gas and steam turbines
that are used for electric energy generation. Since the end of the nineteenth
century, attempts at electric generation from wind power have been made;
nevertheless, combustion-based turbines had much greater development and
applications due to their higher reliability and the possibility of being used
whenever necessary. Due to oil crises in the 1970s, wind energy has been
growing and has been the subject of major investment in Europe and North
America, and more recently, in China. In the 1980s, the power electronics
technology and aerodynamics of wind turbines were sufficiently mature to
begin large-scale manufacturing of wind turbines for commercial use. Today,
wind energy is one of the most cost-effective methods of electricity generation
and it is the most important in terms of power capacity installed as a renewable
energy source.
149
150 4 Wind Energy
9
11
2 8 10
5
6
12 7
1. Blades
1 2. Rotor
3. Pitch
4 13 14
4. Brake
5. Low-speed shaft
6. Gearbox
7. Generator
8. Controller
9. Anemometer
10. Wind Vane
15 11. Nacelle
12. High-speed shaft
13. Yaw drive
14. Yaw motor
15. Tower
FIGURE 4.1
Wind turbine structure
its elements. The wind speed increases with the height, and therefore a tower
is used to elevate the turbine. The energy conversion chain can be organized
into:
300 GW
250 GW
200 GW
150 GW
100 GW
50 GW
0 GW
1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013
FIGURE 4.2
Global cumulative installed wind capacity from 1996 to 2013, modified from
[29]
40 GW
30 GW
20 GW
10 GW
0 GW
FIGURE 4.3
Global new installed wind capacity from 1996 to 2013, modified from [29]
Patm
P2 P3
Patm
V1 V2 V3 V4
1 2 3 4
Streamline
FIGURE 4.4
Control volumes for the analysis of a wind turbine
154 4 Wind Energy
P1 V2 P2 V2
+ 1 + z1 = + 2 + z2 , (4.6)
ρg 2g ρg 2g
and
P3 V2 P4 V2
+ 3 + z3 = + 4 + z4 . (4.7)
ρg 2g ρg 2g
If we add Equations (4.6) and (4.7) considering z1 = z2 = z3 = z4 , V2 = V3 ,
and P1 = P4 = Patm , this is obtained:
V42 − V12 P3 − P2
= . (4.8)
2 ρ
If ṁ = ρAV2 = ρAV is substituted in Equation 4.8,
V1 + V4
V2 = . (4.9)
2
Now let’s consider the speed of point 2 as a fraction of the speed: V2 =
(1 − 2a)V1 ; if we substitute this in Equation 4.9 it yields
V4 = (1 − 2a)V1 , (4.10)
The power extracted is equal to the energy difference between point 1 and
point 4, that is
Upper surface
Upper surface
Chord line
Shoulder Chord line
Leading edge
Trailing edge
Trailing edge Leading edge
Lower surface Lower surface
FIGURE 4.5
Symmetrical and unsymmetrical airfoil terminology and shape
4.4.1 Airfoils
An airfoil is the cross-sectional area of a turbine’s blade. There are many
designs of airfoils, depending on the characteristics of the operating environ-
ment. Figure 4.5 shows a scheme of an airfoil. The chord line is an imaginary
line that connects the leading edge with the trailing edge. The shoulder is the
point of higher thickness.
On a symmetrical airfoil, when the airflow comes parallel to the chord line,
the lift force equals zero, because as the upper and lower surfaces are identical,
the air will flow at the same speed and the pressure gradients are equal. On
the other hand, unsymmetrical airfoils have a larger area and are designed to
have a greater velocity of air on the upper surface; this provokes a decrease of
pressure on the upper side that tends to lift the airfoil (see Figure 4.6).
The angle of the chord line with respect to the airflow is known as the angle
of attack (α). Lift and drag force are functions of the angle of attack; as a
convention, lift force is perpendicular to the direction of air and drag is parallel
to it. The physical explanation of lifting is the curvature of streamlines around
the curve, as indicated in Figure 4.7. From fluid mechanics a pressure gradient
PV 2
is necessary to curve the streamlines: δp δr = r , where r is the curvature of
streamlines, V is the velocity of the fluid, and ρ is the density. The pressure
difference between the upper surface and the lower surface gives a force acting
like a centripetal force from the circular motion of the particles. The dragging
force at low attack angles is mainly caused by the friction of the surface with
the blade.
156 4 Wind Energy
Lift Force
R
Drag Force
V
α
FIGURE 4.6
Lifting and dragging force in an airfoil
P0
Decreasing pressure
Low pressure
High pressure
Increasing pressure
P0
FIGURE 4.7
Generation of lift
4.4 Aerodynamics of wind turbines 157
1.20
1.05 Cl
0.90
0.75
Cd
0.60
0.45
0.30
0.15
0.0
0° 10° 20° 30° 40° 50° 60° 70° 80° 90° α
FIGURE 4.8
Lift and drag coefficient as functions of attack angle
Vrel
V
γ
Vtan
a) b)
Vrel
Vtan
γ
FIGURE 4.9
Velocity of the blade relative to the blade at two different locations from the
center
relative velocity; there is the necessity of twisting the turbine, which is the
reason for the form of large blades.
0.4
Φ=5°
0.3
Cp
Φ=10°
0.2
Φ=15°
0.1
Φ=20°
0.0
-0.1 3 6 9 12
λ
FIGURE 4.10
Cp − λ characteristic for different values of pitch angle
Variable-speed turbines can provide greater efficiency that fixed ones, and
in a large period of time, it is economically viable to implement a variable
speed control for obtaining maximum wind energy. The speed control strategy
consists of obtaining the maximum power point of the Cp − λ curve at any
time, as shown in Figure 4.11.
Pn
Excracted Power
Nominal Power
P4
Vn
P3
V5
V4
P2
V3
P1 V2
V1
0
ω1 ω2 . . . ωn
Rotor Speed
FIGURE 4.11
Wind energy conversion system maximum power extraction
wind speed during the day). The short time scale is very difficult to predict
due to the high variability of the resource. On the other hand, medium and
large scale are much more predictable; therefore, the wind generation projects
are studied usually in terms of monthly variations covering one year. This
variation is typically characterized in terms of statistical analysis; one of the
most common distributions employed is the Weibull distribution [19]. This
probability function expresses the probability of having certain average wind
speed vi during a year,
k vi k−1 vi k
p(vi ) = exp − , (4.18)
c c c
where c is the Weibull scale parameter and has the units of the wind speed,
k is the unitless shape parameter, and vi is a particular wind speed. The
cumulative distribution function is expressed as follows:
n h v x io
i
p(v < vi ) = 1 − exp − . (4.19)
c
The Weibull distribution parameters and the average wind speed are re-
lated with the complete gamma function
1
v =c·Γ 1+ , (4.20)
k
4.7 Power generation system 161
where v is the average wind speed, c is again the scale parameter, and k is
the shape parameter. When k = 2, the Weibull distribution becomes a special
case
p π known as Rayleigh distribution, and the gamma function has a value of
2 ≈ 0.8862.
The variation of the hourly mean speed around the annual mean speed is
small, as the shape parameter k is higher while the scale factor shows how
windy (how high the average annual wind speed) a location is. An optimal
scenario for an eolic park location is a place having high scale factor and a
reduced shape factor. The two parameters can be estimated by the maximum
likelihood method [12],
Pn k
Pn −1
i=1 vi ln(vi ) i=1 ln(vi )
k= Pn k
− , (4.21)
i=1 vi n
n
! k1
1X k
c= v , (4.22)
n i=1 i
where n is the number of non-zero data points, vi is the wind speed in the time
i. Equations 4.21 and 4.22 must be solved using an iterative procedure; some
authors [56] recommend starting with the guess of k = 2. Another important
definition that sometimes is employed is the wind power density (WPD), a
quantitative indicator of the wind energy potential of a particular zone. This
value is a combination of the effects of wind speeds, wind speed distribution,
and the density of the air.
The WPD is defined as
vX
i =25
where ρ is the air density, P (vi ) is the probability to have a wind speed vi
during the year. Pv =25
If we substitute the expression of the mean wind speed v = vii =0 vi P (vi );
W P D = 0.5ρv 3i . (4.24)
Figure 4.12 shows a mean wind speed map obtained from statistical studies
[43].
Wind classes at 80 m
< 9 m/s
8 m/s
7 m/s
6 m/s
5 m/s
4 m/s
> 3 m/s
FIGURE 4.12
Global wind speed at 80m
Softstarter
AC Grid
Gearbox
Induction
generator
Capacitive
compensator
FIGURE 4.13
General structure of a fixed-speed wind turbine
because it is the roughest machine and also the most economical. After the
turbine is rotating at the synchronous speed, the blades are positioned at
their maximum efficiency attack angle and are left for extracting some me-
chanical power from the wind stream. When the induction machine is working
like a generator, the rotational speed exceeds the synchronous speed and the
electromagnetic torque is negative, inverting the current flux. The mechanical
speed is very close to the synchronous speed, and the wind induces only small
variations on it.
As shown in Figure 4.11, the aerodynamic efficiency has only a maximum
at certain wind speeds, for the remaining wind speeds, the machine will op-
erate at a point corresponding to the synchronous speed on the x axis of the
curve. Another disadvantage is the incapability of the system to deal with the
mechanical stress resulting from large wind gusts.
On the other hand, fixed-speed wind turbines are mechanically simple and
have a low maintenance cost; there are no power electronic devices and their
efficiency is high despite the fixed speed.
In order to increase the speed range of the generator, old designs are used to
employ two winding sets: one for low wind speed, typically made of eight poles,
and the second one for medium and large wind speeds (typically six or four
poles). As seen in Chapter 2, the synchronous speed is inversely related to the
number of poles. This design made the stator very big and expensive, therefore
the wound rotor induction generator (WRIG) with access to the rotor windings
represented the evolution of the multi-pole induction generators. The rotor
of the WRIG was controlled by a variable resistance, controlled in turn by
power electronics, depending on the measured speed of the wind, changing
the torque-speed characteristic curve of the machine, and modifying the slip
for certain values of extracted torque. The general scheme of this configuration
is shown in Figure 4.14.
Obviously, the range of controllability remains very limited: about 10%
164 4 Wind Energy
Rotor-Side Grid-Side
Converter Converter
Gearbox AC Grid
PMSG
Pref Qref
Control System
Pitch angle
FIGURE 4.14
General structure of a full-converter variable-speed wind turbine
of the synchronous speed. The WRIG and the multi-pole design were the
first attempts at variable speed wind turbines. The next subsection presents
the modern techniques employed for increasing the speed range of the wind
generators.
Rotor-Side Grid-Side
Converter Converter
Gearbox AC Grid
DFIG
Qref
Pref
Control System
Pitch angle
FIGURE 4.15
General structure of a DFIG-scheme wind turbine
nected directly to the electric grid; the advantage of doing this is that only
a portion of the total power (about 30% according to [22]) flows through the
power electronics devices. The DFIG allows variable-speed operation over a
range typically about ±40% of the synchronous speed [13], and the size of
the converter is directly related to the range of speed. The power electronics
at the rotor side are responsible for controlling the active and reactive power
entering or leaving the electric machine, while the converter connector to the
grid side controls the DC-link voltage and ensures operation at large power
factor.
The DFIG wind turbines are highly controllable and allow maximum power
extraction at a wide range of wind velocity. In addition, the active and reactive
power control is fully decoupled because of the independent control of the rotor
current. When the turbine is operating at sub-synchronous speed, the current
flows into the rotor, consuming power from the grid, and is turned back with
the stator windings, while for super-synchronous speeds, both currents from
the rotor and stator flow in the grid direction. The speed of the rotor influence,
the fixed stator speed, and the rotational speed of the electric machine are
shown with the following equation:
nrotor × Npoles
fstator = ± frotor , (4.25)
120
where n is the rpm speed of the rotor, N is the number of poles of the machine,
and frotor is the electrical frequency of the current that flows through the rotor.
e1
e2 3
e3 ~
FIGURE 4.16
Back-to-back converter
R2 L2
L1 R1
C V2
V1
FIGURE 4.17
Back-to-back converter based on voltage source converters (VSC)
u1iL1 u2iL2
C Rp
L1 R1 R2 L2
u2Vc
+ IL1 u1Vc +
V1 + + IL2 V2
- - - -
FIGURE 4.18
Equivalent circuits of back-to-back converter
168 4 Wind Energy
4.8.3 DC bus
The left and right converters are replaced by the voltage source 1 and voltage
source 2, respectively. Then, applying Kirchhoff’s current law,
L1 ẋ1 + R1 x1 + u1 x3 = v1 ,
L2 ẋ2 + R2 x2 + u2 x3 = v2 , (4.31)
1
C ẋ3 + x3 − u1 x1 − u2 x2 = v1 .
Rp
The model of the B2B converter in (4.31) is a discontinuous model since the
4.8 Mathematical model of back-to-back converter 169
FIGURE 4.19
Simulation of tracking control of current i1
inputs u1 and u2 are discontinuous. In order to analyze the B2B converter with
tools for continuous systems, an average model is required. Average models of
B2B converters are obtained by replacing the discontinuous controls u1 and
u2 by pulse with modulated (PWM) signals.
The continuous model of the B2B converter is represented by
L1 ẋ1 + R1 x1 + µ1 x3 = v1 ,
L2 ẋ2 + R2 x2 + µ2 x3 = v2 , (4.32)
1
C ẋ3 + x3 − µ1 x1 − µ2 x2 = v1 ,
Rp
FIGURE 4.20
Simulation of tracking control for current i2
FIGURE 4.21
Simulation voltage in DC bus, Vc
FIGURE 4.22
Pitch angle vs power extracted exercise
4.8 Mathematical model of back-to-back converter 171
4.8.6 Exercises
The objective of this exercise is to obtain a pitch angle versus a power ex-
tracted characteristic curve for the simulated aero-generator and observe the
power extraction depending on the wind speed.
Press ON the buttons “Torque and Enable” on the CH2 AD1 modulus.
Press the Button START to simulate a force of wind of 18 [m/s]. Control
the pitch of the blades to regulate the shaft speed to 1800 [rpm]. Observe the
mechanical power with respect to the Pitch Angle.
With the obtained data, graph the mechanical power versus pitch angle to
obtain the optimal angle of the blades (see Graph 4.1). When the induction
Graph 4.1
Pitch angle vs power extracted graph
machine is rotating at 1800 rpm, close the relay R2 for mounting the electric
generator with the grid.
Once the electric generator and the grid are coupled, increment the speed
172 4 Wind Energy
of the wind and measure the mechanical power extracted from the wind and
the electrical power delivered to the grid. Keep the pitch angle at the optimal
value (31.5◦ ).
Use Diagram 4.1 for connections. Recommended equipment:
• Induction machine
• Brushless control
• Power circuit breaker
• Feeder manager relay
• Three-phase power meter
• Three-phase slip-ring asynchronous machine
• Three-phase supply unit
Diagram 4.1
Suggested electric connection for wind generation
174 4 Wind Energy
TABLE 4.1
Wind properties
Wind
Wind parameter Land based Offshore
Annual average speed at 80 [m] 7.75 [m/s] 8.9 [m/s]
Weibull K 2 2.1
Shear exponent 0.143 0.1
Air density 1.163 [kg/m3 ] 1.225 [kg/m3 ]
Turbine
Turbine parameter Land based Offshore
Turbine hub height 80 [m] 80 [m]
Rotor net Cp 0.39 0.39
Rotor diameter 82.5 [m] 107 [m]
Turbine capacity 1.5 [M W ] 3.6 [M W ]
With this data, let’s obtain the net energy production for year of each
power conversion system, using Equation 4.16:
Land based:
π(82.5)2
1
P = (0.39)(1.165) (7.75)3 = 0.564 [M W ].
2 × 106 4
Offshore:
π(107)2
1
P = (0.39)(1.225) (8.9)3 = 1.513 [M W ].
2 × 106 4
The net annual energy production is defined as the total energy produced
per installed capacity.
[M W ] [M W ](8760)
AEPnet = (365)(24) = (4.33)
Project capacity Project capacity
Now, applying Equation 4.33 to the previous power calculation:
Land based:
(0.564)(8760)
AEPnet = = 3293.8.
1.5
Offshore:
(1.513)(8760)
AEPnet = = 3684.68.
3.6
On the other hand, the cost of installed projects is summed up in the
following table. For land, project costs were calculated with data from an
old project, observed in Figure 4.23. Considering a market adjustment for
Land-installation based on a 1.5 [M W ] turbine Off-shore based on a 3.6 [M W ] turbine
Element cost Installed capacity cost Energy production Installed capacity cost Energy production
[$/kW] cost [$/MWh] [$/kW] cost [$/MWh]
Turbine cost
TABLE 4.2
Soft costs
Construction finance 57.333 1.667 167.437 5.812
Contingency 114.667 3.333 445.5 15.5
Insurance and surety 0 0 279.06 9.687
bond
Total 172 5 893 31
Market price 195 6 0 0
175
adjustment
176 4 Wind Energy
FIGURE 4.23
Installed capital cost for land-base and off-shore wind turbine projects
the price of typical payment in the 2011 market, the energy installation and
production costs estimated were as shown in Table 4.2.
The total installed capacity cost considering the market price adjustment
is:
Land based:
Offshore:
ICC = 1789 + 2918 + 893 = 5600.
For calculating the levelized cost of energy in [$/MWh], the formula
TABLE 4.3
Data for FCR calculation
Data Land based Offshore
Discount rate (d) 8% 10.5 %
Number of operational years (n) 20 20
Effective taxes (T) 38.9 % 38.9 %
Present value depreciation (PVdep) 81.1 % 77.8%
TABLE 4.4
U.S. levelized cost of energy expectancies for 2019
Plant type Levelized capital cost [$/MWh]
Conventional coal 60
Conventional gas combined cycle 14.3
Conventional combustion turbine 27.3
Advanced nuclear 71.4
Geothermal 34.2
Biomass 47.3
Wind 64.1
Wind offshore 175.4
Solar thermal 195
Solar PV 114
Hydroelectrical 72
However, only 59.3% of such power can be taken, as pointed out by physi-
cists Albert Betz, so:
W = 32.6 [kW ].
Problem 4.2. Suppose some given wind park is subject to a wind speed profile
that rises as much as 18 [m/s] and sinks as low as 9 [m/s] (evenly) throughout
a full operational day. If a mean power consumption of 3 [M W ] is desired,
how many 314 [m2 ] ideal wind turbines are needed to reach the mean power
requirement?
Solution 4.2. As both speeds are said to be distributed evenly, the mean wind
speed is 13.5 [m/s]. By using the kinetic energy equation, it is possible to
calculate the individual power of each turbine.
1 1 kg h m i3
Ė = ρAv 3 = 1.225 314 [m2 ] 13.5 =
2 2 m3 s
473.2 × 103 [W ].
Problem 4.3. Using the data from Problem 4.2 and Figure 4.8 (lift and drag
coefficients graph), determine the optimal dimensions of each blade (3 blades)
considering a 20◦ attack angle, an overall design limit of 72.63 [kN m] torque,
and 457 [N ] drag force.
Solution 4.3. The 20◦ attack angle exhibits a peak lift coefficient of about
1.22, while the drag coefficient is as small as 0.08. Both force equations can
be written as follows (the maximum wind gusts rise up to 18 [m/s]):
1 1 kg h m i2
L = ρSV02 Cl = 1.2225 S 18 (1.22) = 241.11S,
2 2 m3 s
h i
1 2 1 kg m 2
D = ρSV0 Cd = 1.2225 S 18 (0.08) = 15.88S.
2 2 m3 s
Knowing that S = c × l, and that the blades are 10 [m] long, both the above
equations can be constrained to determine the chord length (all three blades
contribute equally to such forces):
72.63 × 103 [N m]
= 241.11S = (242.11)(10 [m])c ∴ c = 1 [m],
3 [blade] × 10 [m]
457 [N ]
= 15.88S = (15.88)(10 [m])c ∴ c = 1 [m].
3 [blade]
Therefore, due to design constraints, the chord length must be 96 [cm].
Problem 4.4. For constant torque operation, an attack angle controller is
installed on a wind turbine that operates under the same conditions as in
Problem 4.2. Consider the airfoil area to be S = (10 [m])(0.96 [m]), the attack
angle to be constrained between 0◦ and 15◦ , and the lift coefficient to behave
linearly between (0◦ , 0) and (15◦ , 1). Find the attack angle interval in which
the generator must operate to compensate wind variations.
Solution 4.4.
1 1 kg kg
L(V0 , Cl ) = ρSV02 Cl = 1.225 (9.6 [m2 ])V02 Cl = 5.88 V 2 Cl .
2 2 m3 m 0
This leads to the maximum (18 [m/s]) and minimum (9 [m/s]) operating
points:
Problem 4.5. Measure the wind power density (WPD) for a place where the
wind speed follows the Weibull distribution with parameters k = 2 and c = 9.
Solution 4.5. The Weibull distribution can be solved for these parameters as
follows:
k vi k−1 −( vci )k v vi 2
e−( 9 ) .
i
P (vi ) = e = (0.22)
c c 9
4.11 Simulation
As wind exhibits different speed-rates, the incorporation of an asynchronous
machine is quite transparent, as shown in Figure 4.24 . Its rotor is not neces-
sarily moving at the synchronous speed, so every speed surplus can be taken as
active power generation. However, in the case of a squirrel-cage asynchronous
generator, the magnetization field is taken from the grid, so the machine will
always require reactive power, making necessary the incorporation of parallel
capacitors to compensate for its inductive nature. In addition, speed varia-
tions could lead to an improper power profile, which cannot be amended if
the access to the rotor’s currents is restricted.
This has led to more complex asynchronous machine designs (together
with increasingly complex power electronics boards), such as the wound rotor
double fed induction machine, WR DFIG, which is actually leading the wind
generation industry.
4.11 Simulation 181
FIGURE 4.24
Asynchronous machine integration to the bus simulation
5
Solar Energy
The sun is a huge source of energy that keeps most environmental processes
running. Its light and heat are also paramount for different human activities,
even economic ones like agriculture. However, most of the available energy
delivered by the sun (about 700×106 [TWh] daily) over the Earth’s surface
is not purposely transformed into useful forms of energy, so its enormous
potential is hindered by the inability to capture it.
Besides, it would be impossible to actually take most of the sun’s energy.
Yet there are lots of opportunities to take advantage of it. Its heat can be
directly used in water heaters, and its light can be used thanks to the photo-
voltaic effect of some semiconductors. This has led to a whole industry that
develops technology regarding the different ways in which the sun’s energy
can be converted, mostly to electrical energy.
There are two main ways to use solar energy. Firstly, the sun’s light can
be redirected to a single spot, greatly increasing the amount of energy per
unit area. This can boil or pre-heat water if the light is directed to a metallic
container. Secondly, and perhaps more importantly, the light can be captured
through photo-voltaic cells, directly transforming light’s energy into electric-
ity. Both techniques are referred to as thermal and photovoltaic solar energy
sources, respectively [76].
Although the sun can provide huge amounts of usable energy, there are
still some issues to address regarding its availability. The sun’s power is only
available during the daytime, and can only be used fully for a limited part
of the day and depending on the weather. This makes solar energy subject
to sometimes unknown factors, leading to an intermittent performance, which
make such systems depend on energy storage technologies.
Some of these issues are partially solved by forecasting the amount of
energy on a specific place on earth. Hence, solar power plants can be installed
h/m2
depending on the potential [ kWyear ] of that location, and assuming a 10%
overall efficiency [33, p. 300]. It is also noteworthy that solar panels’ efficiency
will drop to their 80% to 88% nominal performance in 25 years.
It is clear that besides the amount of available energy the sun provides, it
is nonsense to estimate its thorough usage as a justification for solar energy’s
incorporation to daily life. The sun provides (some hundreds of thousands
of times) more energy than what is needed by humans, but such energy is
not readily available, it is intermittent, and the means used to “catch” it are
expensive, polluting, and mostly inefficient. Thus, solar energy should be seen
183
184 5 Solar Energy
not as a flawless solution but a promising area in which all kinds of resources
must be invested so that solar power is taken advantage of effectively.
In the same spirit, this chapter shows data regarding solar energy avail-
ability, and deals specifically with photo-voltaic solar panels to introduce their
behavior, modeling techniques, and the areas in which they can be improved.
FIGURE 5.1
Solar reflectors as heaters — basic principle
the day goes by, so the incidence of the light rays is not always the same,
and to guarantee maximum conversion effectiveness, the reflectors or panels
should move optimally. In addition, the unpredictable conditions of the useful
rays make any control system devoted to solar power extraction operate over
always-changing operating points.
However, most solar panels’ installations are fixed (they do not orient to
face the sun directly), and clouds, leaves, dirt, etc., can shade them (even
partially). This implies that a solar power plant depends on many conditions
that make the system not as reliable as other traditional power sources. So,
their application is sometimes better as complementary systems to help the
energy transfer process of another (more dependable) plant.
For instance, solar thermal facilities are commonly used as a pre-heater for
the incoming water to a thermal power plant. Any amount of energy taken
from the sun will reduce the required combustible material needed to boil the
water. This process is not intended to work on its own as it depends on the
same solar power’s intermittent nature, which, anyway, reduces the required
fuel and pollutant gases.
The above description corresponds to that of many applications of alterna-
tive power sources. Some geothermal units and photo-voltaic systems are also
used in that way. The main benefit of such a cooperative interconnection is
that the required electrical reliability and dependability are still obtained, and
the intermittent drawbacks of the aforementioned technologies are eliminated.
A different way to build reliable solar power plants is to specifically forecast
“ideal” zones where the light incidence can be guaranteed for a specific amount
of minimum energy capacity. Installation of power plants in deserts (e.g., the
Atacama desert at Chile) and off-shore locations make a successful integration
186 5 Solar Energy
FIGURE 5.2
Different solar power technologies [76, pp. 238–240]
FIGURE 5.3
World zones where sun iradiance is higher than 2500 [kWh/m2 /yr]
FIGURE 5.4
Schematic representation of a solar cell [51, pp. 65–69]
188 5 Solar Energy
FIGURE 5.5
Simple setting of solar photovoltaic panel
FIGURE 5.6
Schematic view of a panel arrangement
current controlled, enabling the load controller operation and allowing differ-
ent sources to be merged on the same DC bus almost regardless of the original
electrical conditions of each device.
Once the output has been regulated, there are two options for available
power usage. For stand-alone operation, the solar panel acts like the only
power provider so that different loads can be connected directly to the load
controller if their required DC input matches that of the controller. If AC is
required, an inverter is needed. Stand-alone operation can be seen in Figure 5.5
if the PLL block is disregarded, and the power bus is considered to be only
locally available.
The second option is called grid-connected, and implies that the solar power
is directly injected into the infinite bus. As the public electric grid is AC
driven, the AC inverter is required. In addition, the electric grid operates at
fixed-frequency and fixed-amplitude voltage, so a phase-locked loop (PLL) is
required to drive the inverter. This means the inverter’s generated AC voltage
must match the voltage already present in the power bus, so current injection
can occur. Figure 5.5 shows this operation if the power bus is considered to
be the public electric grid.
FIGURE 5.7
Solar cells/panels arrangements — ideal behavior
FIGURE 5.8
Partial shading effects on electric output
FIGURE 5.9
Different topologies regarding inverters’ connections
the use of additional diodes to avoid panel damage; this is further explained
in Section 5.3.1.
In addition, solar arrangements intended for AC use require inverter mod-
ules to effectively provide or interact with an AC bus. There are different
topologies in which these modules are added to the panels’ arrangements.
Figure 5.9 shows three ways of connecting solar panels to an AC bus.
Notice that the different AC connections represented in Figure 5.9 mainly
differ in the amount of power the inverter module manipulates. Figure 5.9(a)
shows individual inverters associated to a single panel each. This implies that
the power through the inverter is minimal but the synchronization with the
power grid must be performed individually. Figure 5.9(b) presents a “string”
topology that takes series-connected panels into a single inverter module. In
this case, the inverter modules need to manipulate higher power but the syn-
chronization to the power grid is performed considering “blocks” of solar pan-
els. Finally, Figure 5.9(c) shows a single inverter which, regardless of the pan-
els’ arrangement behind it, performs a full-power conversion. Evidently, the
inverter required for “central” operation is the one that processes the maxi-
mum amount of power among topologies; however, it is the only module to
be synchronized to the power grid.
TABLE 5.1
Some characteristics of different cell types, adapted from [72]
Cell type Required space [m2 /kW] Efficiency [%]
Monocrystalline 7 to 9 17–22
Polycrystalline 9 to 11 13–18
Amorphous 16 to 20 7–13
electron donor and acceptor polymers are built and placed together so the
current flux is possible. Organic cells can be further classified depending on
the disposition of both different polymers and other characteristics that can
be parallel furnished to their finishing, such as transparency, light weight,
flexibility, etc.
Besides, organic cells are a promising technique, mostly due to their low
fabrication costs and cleaner deposition, but they also exhibit low efficiency
(about 1/3 toward silicon ones) and chemical degradation [52]. This has hin-
dered their usage in applications where inorganic cells are commonly used.
Organic cells are the subject of current research and are undoubtedly a
topic that will be developed in the years to come; however, inorganic cells
have also made themselves a place in power generation, so a more specific
classification around them is needed as they dominate the current solar power
panorama.
Inorganic solar cells are classified broadly as crystalline and thin-film. Crys-
talline cells are made from semiconductor “patches” sliced from bigger blocks.
Such blocks are normally required to be thick (hundreds of [µm]) due to man-
ufacturing obstacles. Thin-film cells are made through a chemical process that
ends with a semiconductor under 10 [µm] [72, p. 377].
The above classification clearly points to the thickness of the final semi-
conductor block to be used. However, additional features must be considered
as part of such differentiation. For instance, silicon can be arranged to exhibit
a mono-crystalline structure (continuous crystal lattice structure throughout
the block) only for crystalline cell types, while thin-film can only incorporate
amorphous arrangements. Some electrical characteristics in this regard are
shown in Table 5.1.
Table 5.1 and the above classification clearly show the main trade-off re-
lated to solar cell selection. A heavier, more expensive crystalline cell will
exhibit better electric characteristics, while a lighter, cheaper one will make
a poorer use of light incidence. It is worth mentioning that organic cells can
be found to exhibit efficiencies around 9%; however, they still cannot replace
thin-film cells due to some existing issues regarding their operational stability.
194 5 Solar Energy
Current dependent Parasitic Series
on sunlight capacitance resistance
Shunt Load
resistance
p-n junction
FIGURE 5.10
Solar cell equivalent circuit
where G and G∗ stand for present and STC irradiance, respectively, and ∆T =
T − T ∗ . If (5.8) is used together with (5.4) or (5.7), the behavior of a solar
cell can be fully computed for any given electric condition.
It is noteworthy that the previously mentioned ideality factor accounts for
manufacturing issues regarding the p-n junction or materials variability. This
variable is A = 1 for an ideal junction and can go as high as 2 or even more
depending on the manufacturer’s specifications. It can be seen mathematically
that a higher A would result in a “slower” rise of the diode’s current due to an
identical voltage. The solar cell’s model will be clearly affected by this variable
as the current will drop faster for higher As.
Another commonly neglected issue is related to the scale current I0 , some-
times called “reverse bias saturation current.” Any p-n junction will exhibit
a negative current flow if inversely biased. Alternatively, the exponential be-
havior presented in (5.2) uses −I0 to effectively account for such a negative
current, making the exponential equation reach −I0 instead of 0 at −∞ (that
is why it is also referred to as “scale” current). It is important to mention that
a p-n junction cannot hold too long before breaking if inversely biased; hence,
there is a negative threshold voltage after which (5.2) cannot hold anymore.
Solar cells (panels) commonly incorporate a blocking diode at their positive
terminal. This diode prevents the cell for being inversely biased, so the junction
will be protected and I0 will be a mere scale factor. Parallel diodes can also
be used to bypass solar arrangements subjected to partial or full shading,
allowing other connected solar cells/panels to deliver power continuously.
Finally, it is important to notice a feature present in Equations (5.4) and
(5.7). The load current IL can be found on both sides of the equation, mak-
ing it hard to fulfill a voltage-current separation. This difficulty is due to the
exponential nature of the equation, which hinders analytical solving possibili-
ties. For instance, the cell model is usually computed by taking IL (t− ) as the
current inside the exponential part of the equation, using software capable of
variables feedback.
It is possible to find numerical results for the above-presented models by
using a successive approximation method. Such a method can be as simple
as equation balancing (substituting IL until both sides become equal), or it
can be automated through a numeric solver. An analytic alternative is also
possible if the Lambert function W (z), z ∈ C (5.9) is used. In this case the
exponential can be directly solved and the voltage and current decoupled;
however, a way to calculate W (z) must be available.
will hold, having a total amount of N series connected cells. On the other
hand,
N
(i)
X
IL = IL (5.11)
i=1
q Vout
IL = IG − I0 exp + I L Rs −1 series, (5.12)
AKT N
q IL Rs
IL = N IG − N I0 exp Vout + −1 parallel. (5.13)
AKT N
However, cells are not ideal and partial shading may occur in cell arrange-
ments. This forces such solar modules to commonly incorporate bypass diodes,
as shown in Figure 5.11. Suppose the third cell is shaded as shown, so the cur-
rent flowing through it would be much less than the current the remaining
cells could provide. This behavior will make the shaded cell to be seen as a
load, actually exhibiting an inverse polarization. Whether or not a diode is
placed in parallel, the current is free to flow avoiding the shaded panel.
If no bypass diodes are used in series arrangements, a problem commonly
referred to as a hot spot can occur [51, p. 110]. The shaded cell will heat as a
result of its load behavior (power absorption). Panel manufacturers commonly
provide intermediate connections at the output to allow the addition of bypass
diodes for sets of series connected cells. Such sets of cells are commonly rated
under 24 [V] as a rule of thumb [51, p. 113].
As can be foreseen from the above comments, a similar issue regarding
parallel connection will also exist. If two cells or panels are connected in par-
allel and their behavior is different, it is possible to have one of the modules
absorb the energy the other one generates. Figure 5.12 shows this situation
assuming the second module has been somehow partially shaded.
The voltage-current behavior of both the aforementioned modules is also
shown in Figure 5.12. Notice that the difference between both modules is not
large; however, in open circuit operation, the parallel output voltage would
make one of the modules operate at I (1) and the other at I (2) . This clearly
shows that the second module will be taking I (2) = I (1) . A blocking diode can
198 5 Solar Energy
FIGURE 5.11
Bypass diode in series cell arrangements
FIGURE 5.12
Blocking diodes used for parallel connection
Current [A]
Power [W]
Voltage[V]
(a) Graphical view of maximum power point of a solar cell
Environmental
MPPT variables
algorithm
+
A
V
FIGURE 5.13
Maximum power extraction concept
FIGURE 5.14
Block diagram of the solar umbrella monitoring system
the umbrella, so the intensity of the LEDs was dependent on the amount of
ambient light detected. Finally, all the information collected was sent to a
remote computer wirelessly via X-Bee protocol. The remote computer uses
LabVIEW to show the information and analyze it, as shown in Figure 5.15.
This allowed the tracking of the electrical consumption of devices con-
nected to the inverter in real time. This information was used for didactic
purposes in courses of engineering. Also, for research purposes, information
about the efficiency of the solar panel in different weather conditions was
gathered.
5.5 Exercises
Assume that there is a factory connected to a three-phase generator, and the
possibility to add a source of solar energy to reduce the consumption of energy
from the generator. The solar panel will be tilted to recreate the movement
of the sun during the day to observe the differences in power.
Recommended equipment and connections:
• Lamps for photovoltaic trainer
• Photovoltaic inclinable module
• Inverter grid
• Electrical power digital measuring unit
202 5 Solar Energy
FIGURE 5.15
Information on remote computer
Diagram 5.1
Solar energy exercise
TABLE 5.2
Panel’s dependency on incidence angle
Position [deg] Voltage [V] Current [A] Power [W]
◦
90 18.4 1.06 19.504
80◦ 15.3 1.08 16.524
70◦ 14 0.9 12.6
60◦ 15.3 0.86 13.158
50◦ 17.9 0.728 13.0312
40◦ 17.9 0.656 11.7424
30◦ 15.8 0.55 8.69
20◦ 15.1 0.358 5.4058
5. Repeat steps 3 and 4 but modifying the light intensity of the lamps to
75%. Write down the results as shown in Table 5.4.
A solar panel is exposed directly to the sunlight to obtain the most energy
from the sun and convert it to electricity, but there is a downside for this. As
the solar panel is in direct contact with the sunlight, it begins to heat up and
its efficiency starts to degenerate. In this exercise the lamps will act as the
light source, and at the same time, they will heat the panels. Measurements
will be taken every five minutes to observe the efficiency.
Recommended equipment and connections:
TABLE 5.3
One-phase solar contribution
Position [deg] Line 1 [W] Line 2 [W] Line 3 [W]
◦
90 268 292 292
80◦ 232 288 292
70◦ 240 288 288
60◦ 244 288 292
50◦ 248 288 288
40◦ 256 292 288
30◦ 260 288 288
20◦ 272 288 288
TABLE 5.4
Results considering 75% of original light power
Position Voltage Current Power Line 1 Line 2 Line 3
[deg] [V] [A] [W] [W] [W] [W]
90◦ 16.9 0.552 9.3288 264 288 288
80◦ 17.8 0.544 9.6832 264 292 288
70◦ 16.8 0.57 9.576 264 288 288
60◦ 14.4 0.528 7.6032 268 288 288
50◦ 13.3 0.52 6.916 268 288 288
40◦ 14.8 0.344 5.0912 276 288 288
30◦ 14.1 0.304 4.2864 276 288 288
20◦ 14.7 <0.300 - 284 288 288
5.5 Exercises 205
Diagram 5.2
Solar energy exercise connection diagram
TABLE 5.5
Panel’s power dependency on temperature
Time [min] Voltage [V] Current [A] Power [W] Temp. [C]
0 18.4 1.134 20.8656 35◦
5 17.5 1.168 20.44 62.6◦
10 17 1.168 19.856 75.5◦
15 16.8 1.167 19.6056 79.4◦
IL − I0
AKT
VOC = 0.81 = ln =
q I0
2 1.38 × 10−24 K
J
(298 [K]) IL − 1 × 10−6 [A]
ln .
1.6022 × 10−19 [C] 1 × 10−6 [A]
h q i
I = IL − I0 exp (V + IRs ) − 1 .
AKT
−6 q
Ii = IL − 1 × 10 exp (Vi + Ii Rs ) − 1 .
3K(298)
Rs Eq1 Eq2
6.0000000 8.18E+151 3.57E+157
7.0000000 6.80E+168 1.03E+166
6.5000000 2.36E+160 6.06E+161
6.7500000 4.00E+164 7.90E+163
6.6250000 3.07E+162 6.92E+162
6.6875000 3.51E+163 2.34E+163
6.6562500 1.04E+163 1.27E+163
6.6718750 1.91E+163 1.72E+163
6.6640625 1.41E+163 1.48E+163
6.6679688 1.64E+163 1.60E+163
6.6660156 1.52E+163 1.54E+163
6.6669922 1.58E+163 1.57E+163
6.6665039 1.55E+163 1.55E+163
Problem 5.3. Draw the V-I curve of a panel composed of 20 series cells
with the following parameters (each) and at standard test conditions: A = 3,
Rs = 0.06 [Ω], Is c = 3 [A], I0 = 1 [µA].
Solution 5.3. The cell panel equation can be better processed (after substi-
tuting) in the following form:
v
i = 3 − 1 × 10−6 exp 778.83 × 10−3 i + 12.98 .
20
Notice that the voltage has been divided by the number of series cells to
depict series cells connection.
This can be solved through the Lambert equation. However, it can also be
solved numerically for different voltages. The table below shows some iterations
of a simple search method.
The V-I curve can be then plotted if multiple voltage values are solved as
shown in the above table.
208 5 Solar Energy
Problem 5.4. Using the same data from Problem 5.3, plot the effect of chang-
ing the form factor as follows: A = [1, 5] ∈ Z.
Solution 5.4. As the form factor intervenes in the weighing of both voltage
and current inside the diode equation, the limit at which the equivalent diode
will start to take current is modified.
The same considerations for the solution of Problem 5.3 hold. However,
it will be necessary to solve numerically for each form factor. Results can be
seen in the graph below.
Problem 5.5. Find the optimal operating point (maximum power point) of
the panel described in Problem 5.3.
Solution 5.5. The panel equation can be substituted to derive
v
i = 3 − 1 × 10−6 exp 778.83 × 10−3 i + 12.98 .
20
Power can be calculated directly as P = vi, so the preceding equation leads
to
P v
P = 3v − v(1 × 10−6 ) exp 778.83 × 10−3 + 12.98 .
v 20
5.7 Simulation 209
Again, this equation could be solved directly by using the Lambert equation;
however, a numerical approximation is also possible, as shown in the table
below. Hence, the maximum power point can be found at about 16.25 [V] and
2.69 [A]. An analytical solution leads to v = 16.29275 [V] and i = 2.684 [A].
5.7 Simulation
Solar panels can be seen as a direct current source, similar to a battery. This
implies they require an AC inverter so they can be integrated to an electric
network. In addition, their electric output varies depending on temperature
and light irradiance, so their power profile must be manipulated by voltage
converter circuits to guarantee the maximum power extraction and the proper
synchronization to the grid.
Even though automatic control and power electronics are currently inte-
grated in commercial solutions to install solar panels domestically, it is im-
portant to understand the effects of power tracking, cooling, and solar panel
variations with respect to reported nominal characteristics.
The block diagram shown in Figure 5.16 exhibits a solar panel as a simple
block. This model could be used to sweep the desired environmental parame-
ters to test panel’s performance.
210 5 Solar Energy
FIGURE 5.16
Simulation of a solar panel
6
Electric Power Transmission
Transmission of electric power is the key for the development of the modern
electric system in the world. The necessity of carrying power at high voltage
through tall steel towers that cross countries in all directions has enhanced
the analysis of the transmission lines with equivalent circuits (inductances,
resistance, and capacitance). There are underground transmission lines and
overhead ones; in this chapter, the latter will be analyzed due to their greater
popularity.
A network of transmission lines is referred as “power grids.” Traditional
power grids transmit three-phase alternating current at voltages from 69 [kV ]
to 765 kV from the generators to the load. Less impressive but not less impor-
tant are the city’s and railway’s power lines, which transmit power generated
to the loads.
In this chapter, series resistance, series inductance, and shunt capacitance
for three-phase power lines will be analyzed, after which, the ABCD model of
transmission lines is presented and its π equivalent circuit.
211
212 6 Electric Power Transmission
Overhead groud
wire
Beam
Insulator string
Conductor bundle
Tower window
Waist
Legs
Tower base
FIGURE 6.1
Waist-type power tower components
6.2 Transmission line parameters 213
Figure 6.1 shows a transmission line of four conductors per phase, decreasing
the strength of the electric fiend at the surface in order to reduce the losses
due to the corona effect.
Insulator strings of discs, made mainly of porcelain, are used to suspend
transmission lines and separate them from the tower. Figure 6.1 shows a V-
shaped arrangement in the middle bundle for avoiding swings due to air cur-
rents. The porcelain discs are designed to support mechanical loads of about
7500 [kgf ].
The transmission tower is the supporting unit for overhead transmission
lines carrying the transmission cables at a safe height from the ground. Many
designs are implemented according to economic factors, the ambient condi-
tions, and mechanical loads (wind loads, weight of cables, etc.).
Although high voltage towers are made of steel, the transmission lines of
cities can be made of wood and concrete because they are not designed for
discharging energy to the ground.
TABLE 6.1
Resistivity and temperature constant of some conductor metals
Material Resistivity at 20◦ C Temperature
[Ωm] × 108 constant (◦ C)
Annealed copper 1.72 234.5
Hard-drawn copper 1.78 241.5
Aluminum 2.85 228.3
Silver 1.59 243
Brass 6–9 480
Iron 10 180
Steel 12–88 180–980
6.2.2 Inductance
To model the three-phase inductance of a transmission line, let’s start with the
solid cylindrical conductor inductance model. Inductance is present due to the
permeability µ of the medium. The permeability of the air is very close to the
vacuum’s permeability, µ0 = 4π × 10−7 [H/m], and also, the permeability of
non-ferromagnetic materials is very close to µ0 , 0.999 and 1.000022 for copper
and aluminum, respectively.
The steps followed for determining the inductance due to both the internal
and external magnetic fluxes are:
1. Obtain the magnetic field intensity due to the current flowing through the
conductor (Ampere’s law)
˛
H dl = Ienclosed . (6.3)
B = µH. (6.4)
6.2 Transmission line parameters 215
φ = BA. (6.5)
L = λ/I. (6.7)
Assumptions:
• Current distributes uniformly through the cross-sectional area of the con-
ductor (no skin effect)
• Conductor material is non-magnetic(µ = µ0 )
• Electric conductors are made of a single solid bar (or bundle of bars)
Consider the electric conductor shown in Figure 6.2. The current is flowing
out to the page, r is the radius of the cross sectional area, and the lined internal
circle is the path of integration for obtaining the magnetic field. Also consider
that the conductor has a length l.
internal
l
x
r
external D1
D2
FIGURE 6.2
Solid cylindrical conductor and paths of integration
Once we have obtained the magnetic flux density, let’s consider the lined
area shown in Figure 6.2 for obtaining the differential magnetic flux dφ:
µ0 xl
dφx = Bx ldx = ldx. (6.12)
2πr2
For obtaining the differential flux linkage dλ it is important to remember
that only the enclosed fraction of the current proportional to (x/r)2 is linked.
So, the actual differential internal linkage flux is
x 2 µ0 x3 I
dλint = dφx = ldx. (6.13)
r 2πr4
Integrating the differential linkage flux from 0 to r,
ˆ r
µ0 I µ0 I
λint = l x3 dx = l. (6.14)
2πr4 0 8π
λint µ0 l 1 × 10−7
Lint = = = l. (6.15)
I 8π 2
6.2 Transmission line parameters 217
−7 D −7 D
λtotal = 2 × 10 Il ln 1 = 2 × 10 Il ln . (6.21)
re− 4 r0
1
r0 = re− 4 . (6.22)
D
Ltotal = 2 × 10−7 l ln . (6.23)
r0
218 6 Electric Power Transmission
Ik
D3n
Dnp
I3 D3p
p
D2n D2p
I2 D1p
D1n
I1
FIGURE 6.3
Array of N conductors
FIGURE 6.4
Transposition tower
• The conductors are identical (r’ is equal for the three phases)
Let’s consider only phase “a” of the system: in three-phase balanced sys-
tems usually analysis is reduced and represented as a single-phase-per-unit
system,
−7 1 1 1
λa = 2 × 10 l Ia ln + Ib ln + Ic ln . (6.36)
Daa Dab Dac
If distances between conductors are not equal, inductances per phase are
different and as this is not convenient for a power system, a technique called
transposition is implemented. Transposition consists of interchanging posi-
tions of the phases in order to maintain an average constant distance between
them.
In Figure 6.4 with a scheme of a transposition tower is shown, a double-
circuit line, which conductors of the same phase usually place opposite each
other, on the left side. Previous to the transposition the phases from top to
bottom could be A-B-C, while on the right, C’-B’-A’. Then, phases of the left
side are transposed to C-A-B, while the right side changes to B’-A’-C’. The
opposite phase issue is made in order to reduce the inductance because the
fluxes counteract each other.
Assume a complete transposed line. From Equation (6.36) let’s obtain the
6.2 Transmission line parameters 221
phase a
Dab
phase b
Dac
Dbc
phase c
FIGURE 6.5
Transposed three-phase line
2 × 10−7 l
1 1
λaverage = 3Ia ln 0 − Ia ln ,
3 ra Dab Dbc Dac
√
3
Dab Dbc Dac
λaverage = 2 × 10−7 l Ia ln . (6.41)
ra0
As all the phases have the same average distance, inductance of any of
them will be √
3
Dab Dbc Dac
L = 2 × 10−7 l ln . (6.42)
ra0
2’
. M’
. N . .
1
Da1 . . Da1’
1’
a
Conductor a Conductor b
FIGURE 6.6
Single-phase system with composite conductors
In this case, the linking flux λa is not equal to the total flux, because only
the fraction of current 1/N is linked by the magnetic flux, so the actual linkage
flux is
" N M
#
−7 1 1 X 1 1 X 1
λx = 2 × 10 l ln − ln =
N N n=1 Dxn M m=1 Dxm
" N M
#
−7 1 X 1 1 X 1
2 × 10 lI ln − ln . (6.44)
N 2 n=1 Dxn M N m=1 Dxm
Phase B B’
Phase C A’
FIGURE 6.7
Double circuit line
If
v
uN M
uY Y
MN
GMD = t Dxm , (6.46)
x=1 m=1
v
uN N
2
u YY
Nt
GMR = Dxn , (6.47)
x=1 n=1
then
GMD
Lx = 2 × 10−7 l ln . (6.48)
GMR
The same procedure can be made for a three-phase system. From the
previous section we know the inductance of a three-phase system. Replacing
ra0 by GMR in Equation (6.42) for a complete solution, we get:
√
3
Dab Dbc Dac
L = 2 × 10−7 l ln . (6.49)
GMR
Similarly, we can obtain the equivalent radius of each phase for the GMR:
p p
GMRA = 4 (ra0 Daa0 )2 = ra0 Da a0 , (6.54)
q
GMRB = rb0 Dbb0 , (6.55)
p
GMRC = rc0 Dcc0 . (6.56)
1
where ra0 is re− 4 for a solid conductor and the geometric mean radius of a
stranded conductor.
For a complete transposed line, the average GMR is
p
GMRaverage = 3 GMRA GMRB GMRC . (6.57)
2. Obtain the voltage between conductors (or ground) from the electric field:
ˆ D2
V = Ex dx. (6.60)
D1
6.2 Transmission line parameters 225
P1 Electric Field
+
+ + + + ++
+ + ++ + + + + +
D1
V12 +
– +++ + + + + ++
+
P2 r +++ + +
+ +
D2 +++
x
FIGURE 6.8
Solid cylindrical conductor with uniform charge distribution and Gaussian
surface of integration
q
C= . (6.61)
V
Assumptions:
• Conductor is uniformly charged, solid, and cylindrical
• Conductor is very long and end effects are negligible
• Resistivity of the conductor is zero, so there is no electrical field inside it
The electric field inside the conductors equals zero, because we have estab-
lished that we are analyzing a perfect conductor. The electric field outside can
226 6 Electric Power Transmission
FIGURE 6.9
Unsymmetrically spaced, complete transposed 3-phase transmission line
a’ c’
b’
FIGURE 6.10
Image method for a three phase transmission line
uniform distribution, the geometric mean distance and geometric mean radius
remain valid.
So, for bundled conductors with well-defined distance between phases, the
equation
qa 2π0 l
Can = = √
3 (6.73)
Va ln Dab Dbc Dac
GMR
is a good approximation, while the exact model used for not-well-defined dis-
tance between phases, as that of double-line transmission lines, is
qa 2π0 l
Can = = GMD
. (6.74)
Va ln GMR average
Equations (6.46), (6.47) and from (6.50) to (6.57) of the inductance section
are valid for obtaining the GMD and GMR in Equation (6.74). The only
difference is that the distance of a conductor with respect to itself is the
radius (Dxx = rx ) instead of rx0 = re(−1/4) . This comes from the assumption
that the electric field within the conductor is zero, while for the magnetic field
this is not valid because all the conductors have a magnetic permeability.
6.2.6 Exercises
Calculate the following transmission line parameters for the transmission line
modules.
Consider a 360 [km] transmission formed by bundles of conductors. The
bundles have the following parameters:
12 m 12 m
FIGURE 6.11
Length of transmission lines
Obtain the resistance in [Ω], the inductance in [mH], and the capacitance
in [µF ].
Obtain the π equivalent circuit placing half of the total capacitance at
each end of the transmission line and the total impedance and resistance in
the middle.
Insert all bridging plugins to the capacitance of the transmission line
model, and with an impedances meter, obtain the capacitance of one end
of the transmission line with respect to neutral and compare it with the cal-
culated value.
Also obtain the inductance of one line and compare your result with the
inductance calculated.
Calculate the transmission line parameters for the second module with the
configuration shown in Figure 6.12 and the following parameters:
Diagram 6.1
Connections for different loads
232 6 Electric Power Transmission
2.5 m 2.5 m
FIGURE 6.12
Length of transmission lines
Insert all bridging plugs connecting the capacitances to the line, then close
the Relay R1 and check the output voltage coming from the secondary trans-
former.
Close the Relay R2 and monitor the electrical parameters.
Set the load R and L, C to zero and take measurements of the power
consumed with the maximum demand meter by the shunt capacitance of the
line.
Measure the reactive power consumed.
Increment the capacitive and inductive load in a balanced way and observe
the increment of the current consumed and the reactive power.
Compare the measured charging reactive power with that which it requires
according to the calculation:
2
QC = ωCan VLL ,
QL = ωLL IL2 ,
QLoad = QL − QC .
Fill the following table from the previous table and the above calculations.
• Capacitive load
• Overhead line model — Long
• Overhead line model — Medium
• Power circuit breaker
• Three-phase power meter
• Three-phase transformer
• Three-phase supply unit
• Feeder manager relay
+ +
x+∆x x
FIGURE 6.13
Differential section of a transmission line
pacitance:
z = R + jωL [Ω/m],
y = jωC [S/m].
Rearranging:
V (x + ∆x) − V (x)
zI(x) = . (6.81)
∆x
If ∆x → 0,
dv(x)
zI(x) = . (6.82)
dx
On the other hand, for the current, the KCL equation gives
1 d2 V (x)
yV (x) = , (6.85)
z dx2
which can be expressed as
d2 V (x)
− yzV (x) = 0. (6.86)
dx2
6.3 Transmission line mathematical model 235
IS IR
+ +
VS ABCD VR
– –
FIGURE 6.14
ABCD representation of a transmission line
A1 eγx + A2 e−γx
I(x) = , (6.88)
ZC
q
z
where ZC = y is called the characteristic impedance. Now, solve A1 and
A2 with the following boundary conditions: x = 0 → V (0) = VR = A1 + A2 ,
I(0) = IR = A1Z−A
C
2
.
Substituting A1 and A2 leads to
where
√ z
r
γ= yz, Zc = . (6.92)
y
236 6 Electric Power Transmission
IS Z’ IR
+ +
VS Y’/2 Y’/2 VR
– –
FIGURE 6.15
Transmission line equivalent pi circuit
sinh (γl)
Z 0 = ZC sinh (γl) = ZF1 = Z , (6.93)
γl
Y0 tanh (γl/2) Y Y tanh (γl/2)
= = F2 = . (6.94)
2 ZC 2 2 γl/2
6.3.3 Exercises
Assemble Diagram 6.2 using the large transmission line model and connect
all the plugins in the transmission line model. Set the primary side of the
three-phase transformer in ∆ connection 380 [V ] and set the secondary-side
to Y UN − 5%.
Take measure of the voltage V, current I, and power consumed (both active
P and reactive Q) of the entire power system (including transmission lines) and
the voltage and current after the transmission lines. Fill the following table,
where 1 is the input measure of the transmission line and 2 is the measure at
the output of the transmission line.
Observe the increase in reactive power due to the increment of current
flowing through the transmission line’s inductor. Also, observe the voltage
drop due to the transmission line impedance and how it increments with the
current.
6.3 Transmission line mathematical model 237
Diagram 6.2
Connections for matched-load practice
238 6 Electric Power Transmission
The inductive effect of the transmission line is always much larger than
the capacitive effect, so the reactive power will be inductive even with a very
low current flowing through the transmission line.
Power companies must maintain the voltage supplied to the users within
a range typically ±5% of the nominal voltage, maintain the load impedance
constant (R4 = 312 [Ω]), and obtain the power and voltage drop of the trans-
mission line model.
Supply V1 I1 P1 V2 I2 P2 ∆P η Voltage
voltage [V ] [A] [W ] [V ] [A] [W ] [W ] drop [%]
UN+5% 409.4 0.73 519.5 398.2 0.73 508.1 11.4 0.978 97.264
UN 392.6 0.7 478.2 382.3 0.7 466 12.2 0.974 97.376
UN–5% 376.3 0.67 432.9 364 0.67 422 10.9 0.975 96.731
UN–10% 354.9 0.64 391.2 345.8 0.64 381.7 9.5 0.976 97.436
UN–15% 336.4 0.6 351.5 327.7 0.6 342.1 9.4 0.973 97.414
• Resistive load
• Inductive load
• Capacitive load
where the diagonal elements Ykk are the sum of admittances connected to the
bus k. The off-diagonal elements Ykn , where k 6= n, are the negative of the
sum of the admittances connected between bus k and bus n.
The negative signs are present because for convention, it is assumed that
all the currents are leaving the bus,
I = Ybus V. (6.97)
240 6 Electric Power Transmission
where Ik∗ is the conjugate of the current of node k. If we substitute Ik∗ with
the current of each bus, we get
"N #∗
X
Sk = Pk + jQk = Vk Tkn Vn . (6.100)
n=1
By using polar notation, Vn = Vn ejδn ; Ykn = Ykn ejθkn , and separating real
and imaginary parts, we get
N
X
Pk = Vk Ykn Vn cos (δk − δn − θkn ), (6.101)
n=1
N
X
Qk = Vk Ykn Vn sin (δk − δn − θkn ). (6.102)
n=1
There are several methods for solving the power flow problem. One of
the most popular because of its simplicity is the Gauss-Seidel method. This
iterative method can be computed as
" k−1 k−1
#
1 Pk − jQk X X
Vk (i + 1) = − Ykn Vn (i + 1) − Ykn Vn (i) , (6.103)
Ykk Vk∗ (i) n=1 n=k+1
N
X
Pk = Vk (i) Ykn Vn (i) cos (δk (i) − δn (i) − θkn ), (6.104)
n=1
N
X
Qk = Vk (i) Ykn Vn (i) sin (δk (i) − δn (i) − θkn ). (6.105)
n=1
6.4.1 Exercises
This set of exercises requires a connection like the one shown in Diagram 6.3.
Obtain the active and reactive power consumed by an RL load. Enter the
measured values into the following table:
6.4 Power flow analysis 241
Now remove the connection to the resistive load and repeat the measure-
ment for L4 = 1.27 [H].
The inductive load also consumes active power due to ohmic resistance and
iron losses of the inductor.
Now observe the power consumption due to RC loads. Enter the measured
242 6 Electric Power Transmission
Diagram 6.3
Connections for RC and RL power consumed practice
6.4 Power flow analysis 243
Remove the connections to the resistive load and repeat the measurements
for C3 = 5[µF ].
Naturally, the loads have an inductive and resistive nature. In the indus-
try, it is important to keep the power factor within a certain limit to avoid
economic penalties. Capacitor banks are used for this purpose.
Recommended equipment and connections:
• Resistive load
• Inductive load
• Capacitive load
• Overhead line model — Long
• Overhead line model — Medium
• Power circuit breaker
Diagram 6.4
Connections for capacitive compensation 1/2
246 6 Electric Power Transmission
Diagram 6.5
Connections for capacitive compensation 2/2
6.5 Theoretical problems 247
Electric Resistance
(Ω / Conductor / km)
DC 50 [Hz] 60 [Hz]
0.188 0.189 0.2
Solution 6.2.
0.2
% Increase = × 100 = 106.16%.
0.1884
Problem 6.3. Consider a complete transposed 60 [Hz] bundled transmis-
sion line conformed by two ACSR “Drake” conductors with the characteris-
tics below. Calculate the inductance and series reactance per [km] for a flat
configuration with horizontal spacing between phases of 10 [m] (ab = 10 [m],
bc = 10 [m]). Also consider a bundle spacing of 0.5 [m] (d = 0.5 [m]).
ab bc
A B C
d
Outside GMR Resistance 60 [Hz] Current
Diameter [cm] [m] 50◦ C ([Ω/km]) Capacity [A]
2.8143 0.01143 0.08 900
Solution 6.3. Let’s call the GMD of each single conductor Ds, then the geo-
metric mean distance of each phase bundle is:
v
uN N
2
u YY p p
Nt
GMR = Dxn = 4 (d × DS )2 = (d × DS ) =
x=1 n=1
√
0.5 × 0.01143 = 0.0756 [m].
The series reactance is formed by the resistance and the inductive reac-
tance:
Problem 6.4. Now consider a double circuit line formed by the same con-
ductors with the following distances between conductors. Calculate the induc-
tance and reactance. Calculate the inductance and series reactance per km
and obtain the inductive reactance percentage reduction compared with the flat
configuration.
ac'
A C'
ab'
ab
B B'
aa'
C A'
Solution 6.4. First, let’s obtain the GMD from Equations (6.50) to (6.57):
√
GMRA = 0.01143 × 22.36 = 0.506 [m],
√
GMRB = 0.01143 × 10 = 0.3381 [m],
√
GMRC = 0.01143 × 22.36 = 0.506 [m],
p3
GMRaverage = 0.5062 × 0.3381 = 0.4424 [m].
Problem 6.5. For Problem 6.4, obtain the capacitance with reference to neu-
tral in [F/km] and the shunt admittance in [Siemens/km].
Solution 6.5. For the flat configuration, the radius of the conductor is:
0.02814 m
r= = 0.01407 [m].
2
Then, the GMR for the bundle is:
√ √
GMR = d × r = 0.5 × 0.01407 = 0.08387 [m].
2π × 8.854 × 10−12 × l
Can = √
3
10×10×20
= 1.1099 × 10−11 × l [F ].
ln 0.08387
2π × 8.854 × 10−12 × l
Can = = 1.7139 × 10−11 × l [F ].
ln 12.5971
0.4905
has the following impedance and admittance (solutions of Section 6.2.6 for the
double circuit line series impedance and shunt admittance):
Ω
z = 0.04 + j0.2525
km
y = j6.4614 × 10−6 [S/km].
Calculate the exact ABCD model of the transmission line and present the
result in matrix form. Note: To obtain the square root of a complex number,
it is convenient to employ DeMoivre’s identity: (cos x + j sin x)n = cos nx +
j sin nx.
Solution 6.6. To use Equation (6.91),
p p
γ = (j6.4614 × 10−6 ) × (0.04 + j0.2525) = 1.6518 × 10−6 ∠171◦
= 1.2852 × 10−3 ∠85.5◦ ,
s
0.04 + j0.2525 √
Zc = = 39565.5252∠−9◦ = 198.9108∠4.5◦ .
j6.4614 × 10−6
γl = 0.0504 + j0.6406,
γl 0.0504 j0.6406
e =e ×e = 1.0517 × (0.8017 + j0.5977) = 0.8431 + j0.6286,
−γl −0.0504
e =e × e−j0.6406 = 0.7623 − j0.5683,
eγl + e−γl
cosh γl = = 0.8027 + j0.03015,
2
eγl − e−γl
sinh γl = = 0.0404 + j0.59845.
2
Then,
VS 0.8027 + j0.03015 −1.3284 + j199.3017 VR
= .
IS (4.3853 + j29.8342) × 10−4 0.8027 + j0.03015 IR
Problem 6.7. For Problem 6.6, obtain the nominal π circuit and then cal-
culate the exact π circuit model. Compare the nominal circuit with the exact
model.
Solution 6.7. The nominal model is obtained by directly substituting the
impedance and shunt admittance of the line.
Then,
γl
= 0.0504 + j0.6406 = 0.0252 + j0.3203,
2
eγl = e0 .0252 × ej0.3203 = 0.9733 + j0.3229,
e−γl = e−0.0252 × e−j0.3203 = 0.9255 − j0.3070.
So,
γl eγl − e−γl
tanh = γl = 0.02795 + j0.3315,
2 e + e−γl
Y0 Y tanh (γl/2)
= = 1.6726 × 10−3 ∠89.679◦ .
2 2 γl/2
In addition,
2
0.46 [in] 1 [m]
2
A = πr = π = 107.22 × 10−6 [m2 ],
2 39.37 [in]
1.99 × 10−8 [Ωm]
ρ60◦ C −6 Ω
Rdc = = −6
= 185.6 × 10 .
A 107.22 × 10 [m ] 2 m
DC 50 [Hz] 60 [Hz]
100% 100.532% 106.383%
Solution 6.10. The reactance can be calculated for each frequency as follows:
Xf = Rf + j2πf L.
Thus,
X50
= 186.59 × 10−6 + j98.65 × 10−6 [Ω/m],
m
X60
= 197.45 × 10−6 + j118.38 × 10−6 [Ω/m].
m
Problem 6.11. Consider a transmission line that exhibits the following pa-
rameters:
• Series impedance Z = 0.02 + j0.3 [Ω/km]
254 6 Electric Power Transmission
FIGURE 6.16
Power vs load plot
6.7 Simulation
Transmission lines can be seen as a load themselves due to the non-ideal ef-
fects they impose on the network see Figure 6.17). A pure resistive load will
consume pure active power; however, whenever it is placed after a transmis-
sion line, the generators will have demands for reactive power as well. The
transmission lines also permit the voltage levels to be different at their termi-
nals, so compensation must occur for the systems to reach desired electrical
magnitudes.
Transmission lines increase network complexity and make it more suscep-
tible to failure; they also make an automatic control system to be required.
However, electric distribution needs transmission lines despite their shortcom-
ings. Depending on their total length, a generator sees non-ideal loads while
the load sees a non-ideal power source. In order to improve network perfor-
mance, short transmission lines are implemented in micro-grid systems, which
generates and consumes energy in a “close” or “local” geographical area.
256 6 Electric Power Transmission
FIGURE 6.17
Transmission lines simulation
7
Power Distribution System, Unsymmetrical
Failures, and Power System Protections
257
258 7 Distribution, failures, and protections
time and location of energy demands, duration and frequency of outages, and
the always-important economic factors.
In Figure 7.1, the basic structure of a power distribution system is shown.
This is the simplest structure where the load is unidirectional from the gener-
ators to the loads. The most basic elements of a distribution scheme are also
presented: electric generators, transformers, circuit breakers, and disconnector
switches. The principal difference between circuit breakers and disconnectors
is that the former are used to protect the circuit from overloads and short cir-
cuit conditions, while the disconnection switches are employed when a circuit
is required to be disconnected from the power supply. This element does not
have any protection at high voltages and it requires a circuit breaker before
the switch blades are opened.
FIGURE 7.1
Basic structure of a power system
260 7 Distribution, failures, and protections
FIGURE 7.2
Basic radial distributions
7.1 Types of distribution systems 261
FIGURE 7.3
Aerial radial network
FIGURE 7.4
Subterranean radial networks
kind of network needs switches for interconnecting one branch with the oth-
ers for maintenance or emergency issues. This flexibility makes subterranean
networks more reliable than aerial networks.
The disadvantages of this configuration are the higher cost of installation
and maintenance, and overdesign, so that the lines can handle the overload
when a contingency occurs.
Nowadays the popularity of subterranean installations has increased for
residential suburban areas. Typically, the power supplies consist of ring net-
works taken from an aerial line. Step-down transformers are connected to the
ring and then the loads. In Figure 7.4 a diagram of a subterranean radial
network is shown.
7.1 Types of distribution systems 263
FIGURE 7.5
Ring network
FIGURE 7.6
Mesh network
flow of any of the energy sources can be incremented in order to supply the
entire system. This type of configuration is implemented when failures are
critical; for example, on a smelter or in a hospital, because all the loads are
connected to at least two power sources. The obvious disadvantages are the
high cost and overdesign of the power system. Figure 7.6 shows a scheme of a
mesh network.
7.1.4 Exercise
Consider the electric installation of a metallurgical plant, where the power
supply is critical and the power factor has to be regulated in order to avoid
penalty taxes. Build the circuit shown in Diagram 7.1 by connecting as shown
7.1 Types of distribution systems 265
Diagram 7.1
Mesh network
in Diagrams 7.2 and 7.3. Notice that the incoming high voltage has to be
reduced. Also observe the connection of the automatic power factor compen-
sator, and take care to have the same order for the phases on both ends of
the transmission lines.
Set all loads to zero position, close the relay, and connect the transmission
lines in parallel, increment the resistance in a balanced way, and observe the
portion of the power that flows for each transmission line:
Now activate the capacitive compensator for power factor correction and
increase the inductance load in a balanced way while keeping the resistance
fixed and observing the value of correction. Finally, imagine that a power line
is undergoing maintenance. Deactivate one power relay and observe the power
flow increment in one line.
Recommended equipment and connections:
• Resistive load
266 7 Distribution, failures, and protections
• Inductive load
• Capacitive load
• Overhead line model — Long
• Overhead line model — Medium
• Power circuit breaker
• Three-phase power meter
• Reactive power controller
• Switchable capacitor battery
• Three-phase transformer
• Three-phase supply unit
• Feeder manager relay
Diagram 7.2
Connections for mesh network 1/2
268 7 Distribution, failures, and protections
Diagram 7.3
Connections for mesh network 2/2
7.2 Faults on power systems 269
I1C
ω ω ω
I2C I2B I0A
I1A
+ + I0B
I0C = System
Real
I2A
I1B Positive Sequence Negative Sequence Zero Sequence
FIGURE 7.7
Symmetrical components explanation
• Phase to ground
• Phase to phase
• Two phases to ground
• Three phases to ground
• Three phases
• Phase to ground with electric arc
Whether or not the real system is balanced, it will present only a positive
270 7 Distribution, failures, and protections
I1C IB
I2C I2B
ω ω
I2A
I1A
+ ω
= IA
I0A IC
I1B
I0B
I0C
FIGURE 7.8
Equivalence between an unbalanced system and its components
where a = ej2π/3 and a2 = aj4π/3 represent phase shifts of 120◦ and 240◦ ,
respectively. Notice that every symmetrical component is equally phase shifted
(by 120◦ multiples) and has equal magnitude Ii ; i = 0, 1, 2; consequently,
construction of real phase characteristics can be made by using a with respect
to A phase. The equations shown above are also valid for voltage magnitudes.
In matrix form,
1 1 1
Ts = 1 a2 a , (7.2)
1 a a2
where Ts is called the symmetrical component transformation matrix [53,
p. 250]. The magnitudes of each component can be computed by inverting Ts
as follows:
adj(Ts )
Ts−1 = , (7.3)
|Ts |
7.2 Faults on power systems 271
where 4
a − a2 a − a2 a − a2
adj(Ts ) = a − a2 a2 − 1 1 − a , (7.4)
a − a2 1−a a2 − 1
|Ts | = a4 − a2 − a2 − a − a − a2 = a4 − 3a2 + 2a.
Before performing the calculations to attain the inverse, some properties
of a should be taken into account:
a4 − a2 = a − a2 , (7.7)
2 2 3 4
a (a − a ) = a − a = 1 − a, (7.8)
2 2 3 2
a(a − a ) = a − a = a − 1. (7.9)
The preceding equation is the basis for determining whether or not the
symmetrical components exist, taking as input the direct measurements of
phases’ electrical magnitudes [49, p. 79].
In three-wire systems, all instantaneous voltages and currents must sum to
zero, as there is no neutral connection. Whether or not the neutral connection
exists, single-phase currents can flow; otherwise, zero phase-sequence must be
substituted as zero in previous equations [53, p. 250]. As single-phase currents
flow through the neutral connection, it can be easily calculated as
IA + IB + IC IN
I0 = = . (7.11)
3 3
Another analysis permits the inclusion of impedance variations between
symmetrical systems,
VA ZA 0 0 I0
VB = 0 ZB 0 [Ts ] I1 . (7.12)
VC 0 0 Zc I2
The preceding equation allows the analysis of some faults like phase dis-
connection (open circuit) [53, p. 250].
272 7 Distribution, failures, and protections
Transformer type
of connection ZN ZN ZN
+ + + + + +
3ZN 3ZN 3ZN
Zero sequence
– – – – – –
+ + + + + +
Positive/negative
– – – – – –
e j30° : 1
e -j30° : 1
FIGURE 7.9
Sequence networks for Y-Y, Y-∆, and ∆-∆ transformers
X3 H3
H1
X2
X1 H2
FIGURE 7.10
Voltage phasor diagram for a ∆-Y step-up transformer with HV side leading
by 30◦
due to the magnetizing currents and reactance are ignored because they are
much larger than leakage reactance and resistance of windings.
In Figure 7.11 the sequence networks for transformers are presented.
Synchronous generators are usually grounded through an impedance Zn .
The impedances of electric machines are not equal for all the sequence net-
works as for transformers; an important problem emerges, since rotating ma-
chines are complex to describe mathematically and many assumptions are
made for delivering impedances’ values. To give the reader a sense, the phe-
nomena considered for obtaining the equivalent circuits are the degree of satu-
ration of the core, the linearity of the magnetic circuit, and the rotating speed,
and others must be considered.
For the positive-sequence network, the magnetomotive force produced by
the positive sequence current rotates at the synchronous speed, and a high
quantity of flux crosses the rotor, such that the saturation provokes a high
value positive-sequence impedance. For the negative-sequence network, the
magnetomotive force rotates opposite to the synchronous speed and then its
speed with respect to the rotor is twice the synchronous speed. In this case,
currents are induced for preventing the magnetic flux from penetrating the ro-
tor, and the negative sequence impedance is smaller than the positive sequence
impedance. The zero sequence impedance is generated because of leakage flux,
harmonic flux, and end turns, and its magnitude is smaller than positive and
negative sequence impedances. Exhaustive analysis can be found in many
classic books on the subject [46, 64, 34].
Also, it is important to know that the impedances needed for analyz-
ing failure conditions in rotating machines are the subtransient ones, which
are smaller than synchronous and transient impedances. The synchronous
impedances are used for normal, steady-state conditions analysis while the
transient impedances are used for stability studies. For transmission lines, it
is obvious that the negative and positive sequence networks are equal be-
cause those are passive, bilateral devices. As explained previously, for com-
plete transposed lines, the average impedance of the phases is equal. For
274 7 Distribution, failures, and protections
Negative sequence
– – –
FIGURE 7.11
Sequence networks for rotating machinery and transmission lines
the zero sequence network the currents in each phase are identical in mag-
nitude and phase. Those currents have to return by both the ground and the
ground wires grounded at towers. The magnetic field patterns caused by zero-
sequence currents is very different from those caused by positive or negative
sequences where the currents have a phase difference of 120◦ . The zero se-
quence impedance of transmission lines usually ranges from 2 to 3.4 times the
positive sequence impedance [30].
Figure 7.11 shows simplified models for network sequences, which do not
consider saturation effects, magnetic losses, machine saliency, etc. Neverthe-
less, these networks are accurate enough for power system studies [38].
Once we have transformed a power system into three-sequence equivalent
circuits, we can obtain the Thévenin equivalent on the failure point of the
power system. As we have studied, the zero and the negative sequence do not
have power sources, and the voltage sources of generators are only present in
the positive sequence. So the Thévenin equivalents of negative and positive
are only an impedance, while for the positive there is a voltage source in series
with an impedance, as shown in Figure 7.12.
Thévenin equivalents
Zth0
+
Zero sequence
–
Zth1
+ +
Positive sequence
Vth ~
– –
Zth2
+
Negative sequence
–
FIGURE 7.12
Thevénin equivalents of sequence networks
Thévenin equivalents of the sequence networks are used. The currents Ia0 ,
Ia1 , and Ia2 are calculated and then transformed to phase currents by Equa-
tion (7.1):
00
Ia 1 1 1 Ia0
00
I = 1 a2 a Ia1 v(t) = Vm cos(ωt + θ).
b
00
Ic 1 a a2 Ia2
In this section the most common symmetrical faults of coupled circuits are
presented.
ZA
A
IA IF
E
n
aE a2E ZB
IF B
IB ZC
C
IC
FIGURE 7.13
Phase-to-ground fault
ZA
A
IA
E
n
aE a2E ZB
IF B
IB ZC
C
IC
FIGURE 7.14
Phase-to-phase fault
leading to
E
I1 = . (7.15)
Z1 + Z2 + Z3
The fault current is that flowing through phase A, which in this case is
If = IA = 3I1 . (7.16)
ZA
A
IA
E
n
aE a2E ZB
B
IB ZC
C
IC
FIGURE 7.15
Phase-to-phase to ground fault
IA = 0, (7.17)
IB = −IC , (7.18)
6 0.
VB = VC = (7.19)
So,
I0 = 0, (7.20)
1
I1 = IB (a − a2 ), (7.21)
3
1
I2 = IB (a2 − a), (7.22)
3
I1 = −I2 . (7.23)
Leading to:
E
I1 = . (7.26)
Z1 + Z2
IA = I1 + I2 + I0 = 0, (7.27)
VB = VC = 0. (7.28)
a2 Z1 I1 + aZ2 I2 + Z0 I0 = a2 E, (7.29)
2
aZ1 I1 + a Z2 I2 + Z0 I0 = aEc, (7.30)
I1 + I2 + I3 = 0. (7.31)
7.2.4 Exercises
7.2.4.1 Line-to-earth fault protection
In this exercise, the feeder manager relay will be programmed in order to break
the circuit when a fault condition appears: in this case a solid fault-to-earth
fault in one phase of the transmission line.
Connect Diagram 7.4 proposed for the fault-to-earth and start the practice
shown in Figure 7.17.
Close contactor 1 and notice the measured voltages of the capacitors’ cur-
rents. Close contactor 2 and observe that the voltage in one phase becomes
zero with respect to earth, as will occur in a phase-to-earth fault condition.
The next step is to configure the feeder manager relay for opening the
7.2 Faults on power systems 279
Zth0
Zth2
c Vth ~ Ia2
–
Zf
Line to line
Zth1 Ia1 3Zf Ia2
a +
b
Zth0
Zth2
c Vth ~
Ia0=0
–
Diagram of a LL fault Sequence network connection for a LL fault
Zf
Zf
c
–
Zf
FIGURE 7.16
Diagram of the electric fault and interconnected sequence network for most
commons faults
280 7 Distribution, failures, and protections
FIGURE 7.17
Connection diagram
power supply when a fault occurs. Close the smart grid SCADA program for
avoiding communication problems with the MSCom configuration program.
Open the MSCom 2.exe and click on search − > on the serial port. In the
emerging window select COM 3 and a Baud rate of 9600 (see Figure 7.18).
Once the device is connected, open the function settings (see Figure 7.19).
Configure the 2Is current with the following info; this allows you to set up
the failure condition as a percentage of the In value (also configurable) and
the time of continuous failure state measured before the device does an action
(see Figure 7.20).
Then configure the action that the device will do when the fault state is
triggered. You can find the do configuration in the menu in Figure 7.21.
Now link the current fault condition with its corresponding time for trigger
with an AND logic function (see Figure 7.22).
Close the MSCOM 2 program and again open the smart grid’s SCADA
program, and connect relays 1 and 2. Observe how the device starts counting
and then open the circuit (see Figure 7.23).
7.2 Faults on power systems 281
FIGURE 7.18
Baud rate configuration
FIGURE 7.19
Function settings
282 7 Distribution, failures, and protections
FIGURE 7.20
2Is current configuration
Diagram 7.4
Connections for line-to-ground fault
284 7 Distribution, failures, and protections
FIGURE 7.21
Configuring the action
FIGURE 7.22
Linking the fault condition
of the transformer and connect it to the +10% voltage (shown in Figure 7.29).
Observe how the failure condition is triggered.
Recommended equipment and connections:
• Resistive load
• Inductive load
7.2 Faults on power systems 285
FIGURE 7.23
Device tripped
FIGURE 7.24
Configuring function settings
• Capacitive load
• Overhead line model — Long
• Overhead line model — Medium
• Power circuit breaker
286 7 Distribution, failures, and protections
FIGURE 7.25
Configuration of do menu
FIGURE 7.26
Example
FIGURE 7.27
Configuring function settings
FIGURE 7.28
DO configuration
• Three-phase transformer
7.3 Electric protections 287
FIGURE 7.29
Generating overvoltage
Diagram 7.5
Connections for undervoltage and overvoltage faults
7.3 Electric protections 289
Active: Implies current flowing from one phase to other phase or to ground.
Solid: It results from an immediate insulation breakdown. It is destructive
and is usually accompanied by explosions and flashing: e.g., a rockfall
crushes a cable.
Incipient: It starts mildly. If the fault develops further with time it can
become solid: corona effect.
Passive: They are not actually faults but conditions that will lead to them.
7.3 Electric protections 291
• Overload
• Overvoltage
Transient: They do not imply permanent damage to the system or its insu-
lation. The system can be used again after a transient fault.
Iin Iout
Protected system
or module
Ie’ Ie’’
Relay Ie’ - Ie’’
FIGURE 7.30
Differential principle
Keeper Disc
Contact
Spring
Lag
Movable
Coil
Contact
Stop
Input Electromagnet
FIGURE 7.31
Induction disc relay
on its input coil, which generates magnetic flux along the electromagnet and
the keeper. The lag short-circuited coil present in one of the electromagnet
legs lags the magnetic flux linked to it, making the movable disc rotate. The
disc is attached to a movable contact on its top, which moves radially with
disc rotation. If the fault condition takes long enough to allow both contacts to
join, the relay is triggered and further equipment can be enabled or disabled.
The disc includes a radial spring on its center that opposes contact movement
and returns to the stop condition if the fault is not present anymore [49,
pp. 186–189].
These relays were commonly known as time-overcurrent relays. Their op-
eration can be understood quite straightforwardly for overcurrent conditions,
where a high current would move the disc quickly while a low current would
take some time before the contact reaches the active state. The stop position
can also be used as a normally closed contact, so two actions can be performed
by the same relay at once.
There are cases in which it is important to know the direction of the
current or the relation towards a certain reference. Whenever the measured
electrical magnitude’s direction is desired as a means of protection, a cylinder
294 7 Distribution, failures, and protections
Magnetic
Core
1st Coil
Contacts
FIGURE 7.32
Cylinder relay
relay based on an induction cup (Figure 7.32) can be used. Imagine that for
a given system, the current is known to lag by 30◦ phase A voltage; whether
or not voltage is connected to the first coil and current to the second coil of
cylindrical relay, the magnetic core at the center will orientate axially to the
magnetic flux generated, moving a contact to a certain operation position.
Whenever the current changes in magnitude or phase, the moving contact
will also move itself towards a contact, announcing the fault condition [49,
pp. 191–192].
Cylinder relays were extensively used for directional sensing, which in-
cludes among its applications phase and polarity fault detection. Its use to-
gether with time-overcurrent relays covered many fault scenarios and was the
prime option in protective relaying during the first half of the twentieth cen-
tury.
Besides the two most common and useful relays mentioned above, there is
another that presents a way to compare magnitudes and calculate tolerable
impedances. Distance relays (Figure 7.33) compare two magnitudes between
its contacts through an electromagnetic bar that is influenced by both of
them, and can tilt depending on the magnitude of each one over a fixed pivot.
If voltage and current are connected to the relay and it is adjusted so it reaches
equilibrium at a certain threshold condition, a change in impedance can be
detected as Z = V /I [49, pp. 192–193].
7.3 Electric protections 295
Contacts
Pivot
I
V
FIGURE 7.33
Distance relay
Very
Inverse Extremely
Inverse
Inverse
Time [s]
Instantaneous
FIGURE 7.34
Distance relay
Time
Distance
3 2 2 1 1
Generator Fault
FIGURE 7.35
Time-delay in trigger characteristic of overcurrent relays
Example 7.1. Determine the CT ratio, pickup, and time-delay settings for
the transformer at circuit breaker 1 assuming non-pilot relaying. Assume that
nominal current is 60 [A], minimum fault current is 300 [A], and maximum
fault current is 600 [A]. Select the CT ratio in order to give 1 [A] at the
secondary current for the maximum load (see Figure 7.36).
The current ratio shall be 60:1, since we want to reduce the trigger signal
from 60 [A] to 1 [A]. Some authors recommend configuring the relay settings
as twice the maximum load and one third of the minimum fault divided by 60.
The resulting configuration is: [1.66, 20].
The operation current is 60 [A] and the minimum current of triggering is
set to be 1.66 × 60 = 100, giving 1.66 pu to prevent false operation. Assuming
no coordination, the relay shall trigger at its fastest configuration.
298 7 Distribution, failures, and protections
1
F1
R
FIGURE 7.36
Example diagram
Voltage
Current
Magnitude
Time
Current
interruption
FIGURE 7.37
Overvoltage caused by current interruption in a transmission line
Phase c
FIGURE 7.38
Faulted circuit breaker and resulting remaining circuit
+ +
- -
FIGURE 7.39
Synchronous generator equivalent circuit
TABLE 7.1
Electric machinery maximum operation time at overvoltage conditions
Generators Transformers
Percentage of Permissible time Percentage of Permissible time
overvoltage of operation overvoltage of operation
105 Continuous 110 Continuous
110 30 minutes 115 30 minutes
115 5 minutes 120 5 minutes
125 2 minutes 130 3 minutes
1. What does the “Time Delay” parameter imply? It makes the system wait
for the amount of time specified to trip after a fault is detected.
2. What is the use of a trip delay on a protection system? It permits the sys-
tem to continue working if the fault detected was momentary or a bad mea-
surement. If the system remains in a faulty condition for a given amount
of time, then there will be no doubt of a faulty scenario. It is like a digital
anti-bump device.
3. What does the “ResetRatio” parameter achieve at step 9? It works like
a hysteresis latch by permitting the system to release the trip only if the
measured magnitude is below a predefined percent from the trip condition.
7.3 Electric protections 301
18
16
14
12
10
Time [s]
0
120 140 160 180 200
Percentage of rated voltage [%]
FIGURE 7.40
Tripping characteristics of overvoltage relay
302 7 Distribution, failures, and protections
So, the system will trip beyond a given voltage but will release the fault
state only if the voltage is measured below the “ResetRatio” percent.
4. What are your conclusions about relay logic in this exercise? How would it
be changed if the relay used was R3? Applied relay logic allows the system
to hold the contactor open even if a malfunction occurs, and disconnects
the grid whenever the protection module is not able to continue working.
R3 has an opposite normal state, so it would let the grid be connected even
in fault scenarios.
Example 7.3. A synchronous generator with the following characteristics is
being used as power supply to an industrial plant:
• 480 [V ] (∆ connected)
• 60 [Hz]
• Synchronous impedance: 0.015 + j0.1 [Ω]
• Full load current: 1200 [A] at 0.8 inductive PF
• Friction losses: 40 [kW ]
• Core losses: 30 [kW ]
If the generator is supplying current at full load and then the load is sud-
denly disconnected, how much voltage will be at the generator terminals in
volts and in percentage of the nominal voltage? According to the recommen-
dations in the text, for how much time will the generator be able to continue
operating?
If the load is disconnected, the terminals voltage of the generator will be
the internal voltage of the electric machine, since the current will stop flowing
and no voltage will be dropt at the impedance of the generator:
EA = Vφ + IA (ZA ),
where:
• EA is the internal voltage generated at phase a, which will also be the
voltage at the terminals because the electric machine is connected in ∆
• Vφ is the external voltage
• IA is the phase a current
• ZA is the generator impedance of phase a
1200
IA = √ ∠ − arccos 0.8A = 692.82∠ − 36.87◦ ,
3
EA = 480 + (692.82∠ − 36.87◦ ) × (0.015 + j0.1) = 532.161∠5.304◦ .
7.3 Electric protections 303
FIGURE 7.41
Equivalent circuit
What will be the minimum voltage that can be applied to the motor before
pullout?
The pullout torque is 110 [N m] with the conditions imposed, applying the
torque reduction formula of the text:
V
τavg = τR − τL .
VR
7.3 Electric protections 305
FIGURE 7.42
Torque speed characteristic
FIGURE 7.43
Torque speed characteristic
conditions may lead to residual current flow, such as harmonic distortion; its
effects must be detected and restricted as representing a power leakage derived
from unwanted variations on energy consumption and generation, which may
damage the system.
This protection is used in equipment or a networks component in which
direction of flow could change as a result of a fault in the power system; for
example, when a short circuit occurs between phases or phases to ground, the
natural direction of the current may change since the location of the fault has
zero voltage.
Residual overcurrent protection is used for feeders with a capacitive current
of the same order of magnitude as the earth fault current. In those circuits, the
phase-to-earth capacitance is high enough for a zero-sequence current to be
detected by a non-directional protection when a short circuit occurs. There-
fore, directional protection and overcurrent protection are complementary.
In order to determine the direction of a fault, a reference (or polarizing
quantity) needs to be provided to the relay. The polarizing variable can be a
zero-sequence voltage or zero-sequence current. In this case we are studying
the current polarization.
Current polarization: The current in the neutral line of a power circuit
can be used as a polarization value. This type of polarization has the impor-
tant advantage of only requiring one current transformer, while the voltage
polarization requires 3 VT. Another advantage is the possibility of measuring
high currents with a much more economical device. The main disadvantages
are: the saturation of CT when a short circuit occurs, producing a false value
in the false residual current. When there are mutually coupled lines, current
reversal may occur.
The zero-sequence current can be defined as:
1
I0 = (Ia + Ib + Ic ),
3
where Ia , Ib , and Ic are the currents of phase a, phase b, and phase c, respec-
tively.
The direction of the fault can be determined with the phase displacement
between the current and the zero sequence current. When the angle of the
phase’s currents exceeds 90◦ of the characteristic angle of the directional relay,
the relay is tripped. The characteristic angle must be determined so that only
changes in current direction trip the relay. Figure 7.44 shows a characteristic
curve of a directional overcurrent protection device.
1. Why does a new current appear when the loads are unbalanced? The new
current is the one that flows through the neutral line due to load unbalanc-
ing.
7.3 Electric protections 307
Non-tripping
zone
Tripping
zone
Characteristic angle
Zero-sequence current
FIGURE 7.44
Tripping characteristic of directional relay
2. How can you know if the unbalancing has occurred on the generator/load
side? The phase (angle) of the current is opposed to that of its corre-
sponding phase whether or not the fault is found on the load side; on the
other hand, the neutral current can have the same direction (angle) as the
corresponding phase’s current if the fault is found on the generator’s side.
3. Why does the disconnection of a generator’s phase trigger the residual
overcurrent fault? The protection is configured to detect the magnitude of
neutral current; in this way, it does not matter if the fault is found at any
side. Briefly, this fault detection is not directed.
4. What happens at s 8 and 9? Why? The resulting neutral current is directed
to a composite angle conformed by both generator and load unbalancing.
This angle depends on which phase is actually faulty and on which side it
is. Knowledge about the neutral current’s angle may lead to a fully opera-
tive directed protection.
5. Why does a new current appear when the loads are unbalanced? The new
current is the one that flows through the neutral line due to load unbalanc-
ing.
6. How can you know if the unbalancing has occurred on the generator or
the load side? The phase (angle) of the current is opposed to that of its
corresponding phase whether or not the fault is found on the load side; on
the other hand, the neutral current can have the same direction (angle) as
the corresponding phase’s current if the fault is found on the generator’s
side.
308 7 Distribution, failures, and protections
Zone 1: This is also called the instantaneous zone because the relay opens
as fast as it can when a measured impedance is within the protection
zone of the R-X diagram. For modern digital distance relays, this zone is
configured for a distance of about 85% of the total length until the next
relay. This is made in order to avoid overreaching the protected line due
to possible errors in measuring of current and voltage.
Zone 2: The remaining 15% of the line must be covered by the second zone, as
seen in Figure 7.45. It is recommended that the second zone be configured
to at least 120% of the protected line. In the traditional relays, zone 2
7.3 Electric protections 309
Zone of protection
FIGURE 7.45
Zones of protection and time to trigger the circuit breaker
was activated by a timer signal trigger, but in the modern devices this is
implemented in software. Some authors recommend setting the zone to up
to 50% of the length of the adjacent transmission line, to ensure that zone
2 does not extend beyond zone 1 of the next relay.
Zone 3: Also known as remote back-up protection, usually set at least 1.2
times the impedance that will be measured for a fault at the remote end of
the second line section. When there are interconnected power systems, the
impedance measured will be much greater than the actual fault impedance;
therefore, more than three protection zones will be unpractical.
In polygonal relays it is possible to work with zero and negative sequence
voltages and currents, making fault resistance results less sensitive to the
positive sequence load. There is no advantage to extending the R2 zone to
the left as in MHO relays; although for convenience, many times it is defined
as equal to R1, the impedance of the line is set in the center of the relay’s
X-R characteristic, and the quadrilateral is aligned with the line angle. The
resistance lines are parallel to the line angle because the fault resistance in a
quadrilateral relay is computed independently of the line resistance.
Quadrilateral elements can detect phase and ground faults by applying
mathematical computation. They can be configured in all directions to pro-
vide the desired sensitivity and can be combined with conventional MHO
relays to maintain the advantages of the well-known MHO relay and provide
additional protection to the arc resistance zone. Differences in the source and
line impedances may cause the reactance element to overreach or underreach.
It is common to reduce the reactive reach for covering high fault resistance
and use the MHO characteristics for giving a nice protection along the line
and avoiding overreaching and underreaching (see Figure 7.46).
Example 7.7. Review questions:
310 7 Distribution, failures, and protections
X
ZL1
R2
R1
Line
angle
R
B
FIGURE 7.46
Quadrilateral relay R-X characteristic
Bus A
Fault
FIGURE 7.47
Electric diagram
FIGURE 7.48
Line diagram
Problem 7.1. Obtain the per-unit equivalent circuit for the zero, positive,
and negative sequence networks taking the generator’s base.
314 7 Distribution, failures, and protections
Solution 7.1. First, let’s obtain the impedances at a common base with the
following equation:
2
Vold Snew
Xnew = Xold × ,
Vnew Sold
where “new” refers to the desired base for the voltage and power, and “old” is
the original base.
2
18 [kV ] 1000 [M V A]
Xm1 = 0.162 × = 0.2025,
18 [kV ] 800 [M V A]
2
18 [kV ] 1000 [M V A]
Xn = 0.02835 × = 0.03544,
18 [kV ] 800 [M V A]
2
20 [kV ] 1000 [M V A]
XT 1 = 0.1 × = 0.1,
20 [kV ] 1000 [M V A]
2
18 [kV ] 1000 [M V A]
XT 2 = 0.081 × = 0.10125.
18 [kV ] 800 [M V A]
For the transmission line, the impedance base is at the high voltage side:
2
Vbase 400 [kV ]2
Zbase = = = 160 [Ω],
Sbase 1000 [M V A]
20
XT L = = 0.125.
160
Problem 7.2. Obtain the Thévenin equivalent for each of the sequences seen
from bus 2.
Solution 7.2. The sequence networks can be obtained from Figure 7.9 and
Figure 7.11. The Thévenin equivalent for the zero sequence is Zth0 = j0.081 +
j0.03544.
7.5 Theoretical problems 315
1 2 +
-
Zero sequence Zero sequence
1 2
+
+
+ -
-
-
Positive sequence Positive sequence
1 2
+
-
Negative sequence Negative sequence
FIGURE 7.49
Negative sequence impedance
Problem 7.3. If a single line-to-ground fault from phase “a” occurs at bus 2
of the system described in Problem 7.2, calculate the fault current in per unit
and in amperes, and the voltages at bus 2.
Solution 7.3. From the previous exercise we know the Thévenin equivalent
for each sequence network viewed from bus 2. From Figure 7.16 we know that:
1
Ia0 = Ia1 = Ia2 = = −j2.468.
(0.1516 + 0.1371 + 0.1164)j
From Equation 7.1, the subransient currents are:
Ia 1 1 1 −j2.468 −j7.406
Ib = 1 a2 a −j2.468 = 0 .
2
Ic 1 a a −j2.468 0
316 7 Distribution, failures, and protections
S 1000 [M V A]
Ibase = √ =√ = 32075.01 [A],
3 [V ] 3 × 18 [kV ]
Ia = j7.406 × 32075.01 [A] = 237.547 [kA].
Problem 7.4. Now consider a solid line-to-line fault at half the transmission
line length. Calculate the voltages of the fault at buses 1 and 2 and the current
of the fault in pu. Note: It is convenient to create a new bus halfway along the
transmission line.
Solution 7.4. From previous problems we know the sequence networks for
the power system. The first step is to obtain the Thévenin equivalent viewed
from half the transmission line:
j0.375
j0.375
Zth0 = j0.1 +
j0.10125 + 2
= j0.1441,
2
j0.125
j0.125
Zth1 = j0.2 +
j0.2020 + j0.10125 + 2
= j0.1528,
2
j0.125
j0.2226 + j0.10125 + j0.125 = j0.1728.
Zth2 = j0.25 +
2
2
Ia0 = 0,
1
Ia1 = −Ia2 = = −j3.0712.
(0.1528 + 0.1728)j
S 1000 [M V A]
Ibase = √ = √ × 400 [kV ] = 1443.37 [A],
3 [V ] 3
Ib = −5.32 × 1443.37 [A] = −7678.76 [A] = −Ic .
79 × 10−3 ∠−86◦
I0
I1 = 97 × 10−3 ∠−74◦ .
I2 3.44∠1.75◦
Problem 7.6. The symmetrical current components were acquired for some
loads and are provided below. If those currents were acquired for 28 [Ω] loads
and the substation voltage is known to be a stiff 120 [Vrms ] source, estimate
the distribution reactances responsible for such unbalancing.
318 7 Distribution, failures, and protections
79 × 10−3 ∠−86◦
I0
I1 = 97 × 10−3 ∠−74◦ .
I2 3.44∠1.75◦
Solution 7.6. This can be computed directly from the following equation:
3.6∠242◦
VA ZA 0 0 I0 ZA + 28 0 0
VB = 0 ZB 0 Ts I1 = 0 ZB + 28 0 3.3∠121.23◦ .
VC 0 0 ZC I2 0 0 ZC + 28 3.6∠242◦
It is necessary to solve as follows:
120∠0◦ 3.6∠242◦
ZA + 28 0 0
120∠120◦ = 0 ZB + 28 0 3.3∠121.23◦ .
120∠240◦ 0 0 ZC + 28 3.6∠242◦
The above system leads to:
ZA → 7.07 − j1.22[Ω]
ZB → 8.36 − j0.78[Ω]
ZC → 5.31 − j1.16[Ω].
• VB = 128∠123◦
• VC = 126∠241◦
If the neutral terminal is known to have a reactance of 0.03 + j0.12 [Ω],
calculate the power lost through it.
Solution 7.7.
V0 1 1 1 VA
1
V1 = 1 a a2 VB .
3
V2 1 a2 a VC
So, V0 = −2.27 − j0.95 [V ], and the power is
V02
P = = 48.84∠−58.5◦ [V A].
Xn∗
Solution 7.8. It is necessary to calculate the current of the neutral line from
the symmetrical components analysis equations:
I0 1 1 1 IA
1
I1 = 1 a a2 IB .
3
I2 1 a2 a IC
VA VB VC
IA + IB + IC 10+j3 [Ω] + (10+j3 [Ω])p + 10+j3 [Ω]
I0 = = =
3 3
1 1
0.93 1 − + j3.95 −1 .
p p
RMS current can be obtained as the magnitude of the above calculated
current, so:
s 2 2
p−1 1−p
kI0 k = 5 = 0.93 + 3.95 ,
p p
25 = 16.46(1 − 2p−1 + 1p−2 ).
The B phase load must be of p = 0.448 its original value to trigger the
protection.
Problem 7.9. Suppose that some symmetrical analysis triggers due to a fault.
After checking the system’s status, all symmetric components are the same.
What happened to the system before triggering?
Solution 7.9. The only possibility is a phase-to-earth fault, as:
I0 IA IA
1 1
I1 = T −1 0 = IA .
3 S 3
I2 0 IA
7.7 Simulation
Electric networks, besides being quite complex, are also subject to failures.
There are some common types of faults, mostly related to transmission lines,
which involve the short-circuit of a combination of conductors. Once a short-
circuit appears, the system must disconnect the faulty section of the grid in
order to prevent further damage to loads and generators.
It is important to notice that current protections and power analyzers are
not fully capable of identifying the fault’s origin, and sometimes, a distance
protection is needed. After all, the full model of the grid must be taken into
account in order to effectively detect the type of fault, its location, and the
correction that can be performed in order to keep the power flowing to users
(see Figure 7.50).
320 7 Distribution, failures, and protections
FIGURE 7.50
Simulation of transmission faults
8
Renewable Energy
321
322 8 Renewable Energy
- Conversion losses
- Distribution losses
Secondary Energy
(coke, briquettes, gas,
- Self-consumption
biodiesel, heating oil, power,
- Non-energetic
wood logs, district heating)
consumption
- Conversion losses
- Distribution losses
Final Energy
(briquettes, gas, heating oil,
- Self-consumption
power, wood chips, natural
- Non-energetic
gas, district heating)
consumption
FIGURE 8.1
Energy conversion chain
rectly consumed by the final user. They result from secondary and possibly
from primary energy carriers, or energies, minus conversion and distribu-
tion losses, self-consumption of the conversion system, and non-energetic
consumption. They are available for the conversion into useful energy.
• Useful energy refers to the energy available to the consumer after the last
conversion step to satisfy the respective requirements or energy demands.
It is produced from the final energy carrier, reduced by losses of this last
conversion.
resources are energy resources that are currently generated, for instance, by
biological processes.
Energy sources, by contrast, provide energy streams over a long period of
time; they are thus regarded as almost inexhaustible in terms of human time.
But these energy flows are released by natural and technically uncontrollable
processes from exhaustible fossil energy resources. Even if such processes take
place within very long time periods, and are unlimited in human time, they
are nevertheless exhaustible. The available energies can be further subdivided
into fossil biogenous, fossil mineral, and renewable energy carriers.
heat from the earth from decay of radioactive particles, and residual heat
from gravitation during formation of the earth. Volcanoes are fiery examples
of geothermal energy reaching the surface from the interior, which is hotter
than the surface. Tidal energy is primarily due to the gravitational interaction
of the earth and the moon.
Overall, 14% of the world’s energy comes from bioenergy, primarily wood
and charcoal, but also crop residue and even animal dung for cooking and some
heating. This contributes to deforestation and the loss of topsoil in developing
countries. Production of ethanol from biomass is now a contributor to liquid
fuels for transportation. In contrast, fossil fuels are stored solar energy from
past geological ages. Even though the quantities of oil, natural gas, and coal
are so large they are infinite and for the long term of hundreds of years, they
are not sustainable.
Due to the great variety of possibilities for using renewable energy sources
with the aim of fulfilling the demand for end or useful energy, it is very difficult
to present the different possibilities in a similar manner.
The possibilities and boundaries for converting renewable energies into
end or useful energy largely depend on the respective physical and technical
conditions. Whenever possible, also among others the theoretically maximum
achievable efficiencies and technical availabilities among others will also be
indicated. These technical key figures are defined below:
• Efficiency, which is defined as the ratio of useful power output to the power
input, depends on the respective operating conditions of the conversion
plant, as well as a series of other factors, which vary over time.
• Utilization ratio is defined as the ratio of the total output of useful energy
to the total energy input within a certain period of time. The observed
time periods may include part load periods and breaks as well as start-up
and shutdown times. Utilization ratios are thus usually smaller compared
to the efficiency of conversion plants indicated to the design point at full
load.
• Technical availability describes that portion of the time period under ob-
servation, within which a plant has actually been available for its intended
purpose and thus considers time periods during which the plant has been
unavailable due to malfunctions.
8.5 Storage
Energy on demand means stored energy, and the most common means of
storage are water in dams, batteries, and biomass. Fossil fuels are stored solar
energy from past geological ages. Storage is a billion dollar idea. Economical
storage would mean no new electric power plants would have to be constructed
for many years, as energy could be stored from existing power plants during
periods of low demand.
Energy cannot be created or destroyed, only transformed from one form
to another, so in reality there are only two forms for storage: as kinetic energy
and potential energy.
Storage as kinetic energy could be as flywheels, and thermal storage could
be as heat. For example, a passive solar home would use concrete or rock and
maybe water for 2–4 days of thermal storage, and a ground source heat pump
would use the earth as seasonal thermal storage. Compressed air is king of
a mixture; it is mechanical, but there is a thermal change. Super flywheels
with high revolutions per minute and composite materials for strength have
been designed and are used in uninterrupted power supplies and prototypes
on buses.
Potential energy is due to generalized interactions, of which we consider
gravitational and electromagnetic for storage systems. The gravitational po-
tential energy is primarily water: dams, tidal basins, and pumped storage.
The electromagnetic interactions include chemical, phase change, magnetic,
electric, and mechanical interactions. Chemical storage is by batteries, pho-
tosynthesis, production of methane and hydrogen, fertilizer, and other types.
The storage of gas requires high pressure, converting it to liquid or as a chemi-
cal compound, or, for example, storing hydrogen in metal hydrides. Of course,
lots of solar energy is stored in chemical compounds as food and fiber.
Thermal storage could be thermal mass, ice, molten salt, cryogenic, earth,
solar pound, and phase change. Molten salt for thermal storage is being used
in combined heat and power systems. The main components for considera-
8.5 Storage 327
tion of different storage systems are: energy density, efficiency and rate of
charge/discharge, duration or lifetime (number of cycles), and economics.
Different storage technologies can be utilized for power or energy. Energy
density and then size and weight for some applications are important factors.
Liquid fuels have large energy density, while hydrogen gas has low energy den-
sity. In general, storage efficiencies range from 50% to 80%, and the lifetimes
vary widely, from 100 years for dams, to 5 to 10 years for lead acid batteries
in a photovoltaic system, to less than an hour for nonrechargeable batteries.
The maximum rate and best rate of charging and discharging the storage
is related to type and use of storage. Whether energy storage is included in an
application and the type of storage are driven by economics and specific power
and energy. An electric car that takes 6 hours to charge its batteries could
not be used on a cross-country trip, even if charging stations were available.
Thus, it can be seen how versatile liquid fuels are for transportation. Also,
the storage requirements for high power over a short time is different from
energy storage for a few days. For utilities, the only large storage systems are
pumped storage and compressed air energy systems; however, battery systems
for power shaving, conditioning and reducing the variability from renewable
energy sources have been installed. Both pumped hydro and compressed sys-
tems have long life and a large number of cycles. For remote village power
systems and stand-alone systems, batteries are the most common means of
storage.
The electricity produced by a generator cannot be stored; energy minus
energy losses is the demand, so the generation supplies that amount of de-
mand, which varies by time of day and season. In addition, a utility system
must meet peak demands and have spinning reserve for unforeseen conditions.
If demand exceeds capacity, then users are taken off the grid. In some parts
of the world, there are rolling blackouts, or electricity is only available for
certain time periods. Finally, extreme events may force shutdown of the total
grid. A good source for information on storage for electricity is provided by
the Electricity Storage Association.
maximum efficiency through the use of pumped storage to absorb their out-
put at night during low demand. There are motor/generator, friction, and
evaporation losses, so the overall efficiency of a storage cycle is around 75%.
The advantages are: improved energy regulation and operation of the sup-
ply grid, ancillary services such as standby and reserve, black-station star,
frequency control, and flexible reactive loading. The disadvantages are the
high capital cost due to fairly large reservoirs, and, as for other dams, the
land area needed, and the fact that collapse could happen.
8.5.3 Flywheels
Flywheels store energy due to rotational kinetic energy, which is proportional
to the mass and the square of the rotational speed,
Increasing the revolutions per minute increases the energy density, so high-
speed flywheels have revolutions per minute in tens of thousands. Low-speed
flywheels are made from steel, and high-speed flywheels are made from carbon
fiber or fiberglass. High-speed flywheels are housed in a low vacuum and use
magnetic bearings to reduce or eliminate frictional losses. Advances in power
electronics, magnetic bearings, and materials have resulted in direct current
8.5 Storage 329
(DC) flywheels. Note that if there is material failure, the container has to
retain that energy inside. Cycle efficiency is around 80%. Flywheels have been
used in trains, cars, and buses; however, they were primarily experimental or
prototypes. In the past, there were flywheels on tractors to smooth out the
rotation of the crankshaft of two-cylinder engines. A hybrid vehicle could use
flywheels for acceleration in conjunction with a smaller internal combustion
engine, much like hybrid vehicles with batteries. In a car, the flywheel is a
gyroscope, which will change the handling.
8.5.4 Batteries
Batteries are common all over the world, as lead-acid batteries are used for ve-
hicles. Batteries are used for low-power and low-energy applications for lights,
radio, TV, and electronic devices. Photovoltaic (PV) with rechargeable bat-
teries has now replaced batteries for very low-power applications. Batteries
convert chemical energy into electrical energy using electrodes immersed in
a medium (liquids, gels, and even solids) that supports the transport of ions
or electrolyte reactions at the two electrodes. Individual cells are placed in
series for higher voltage and in parallel for higher current. Since there is an
internal resistance and because of other factors, there are losses in charging
and discharging a battery.
The power is the product of the voltage and current; however, batteries are
generally specified by volts and storage capacity (battery charge CB), which
is related to stored energy. CB is the amount of charge that a battery can
deliver to a load. It is not an exact number because it depends on the age
of the battery (number of cycles), temperature, state of charge, and rate of
discharge. If a lead-acid battery is discharged to essentially zero a few times,
it drastically reduced its lifetime. As a first approximation, the energy is
where V is the voltage and CB is the battery capacity [Ah]. Decreased tem-
perature results in less battery capacity, and for a lead-acid battery, storage
capacity decreases around 1% for every 1◦ C drop in temperature. Explosions
due to short circuits and generation of hydrogen plus disposal of used batteries
and toxic chemicals are problems.
sodium-sulfur (NaS) batteries for load leveling, enabling the storage of low-
cost off-peak wind power for sale and/or distribution during peak demand
times. The sodium sulfur battery sets are charged primarily at night when
demand for power is lower. The stored electricity is then supplied to the utility
grid, along with power directly generated by the wind turbines, during peak
demand times. This system ensures a steady supply of power to the grid even
during those periods when wind generation falls as a result of low wind speeds.
To control the transmission of power from the Rokkasho-Futamata Wind
Farm to the national power grid, a smart grid monitoring and control sys-
tem from Yokogawa Electric Corporation systems was installed. The system
includes STARDOM network-based controllers and FA-M3 range-free con-
trollers. The Yokogawa system factors in fluctuations in power production as
a result of varying wind conditions and develops a power generation operation
plan. With over 200 [MW] of NaS battery storage devices installed through-
out the world, the company has leveraged their expertise with ceramics and
electrical insulation systems to become a world leader in the development of
sodium sulfur batteries. NGK NaS battery technology is typically used to sta-
bilize intermittent renewable energy, provide emergency power, and as part of
stand-alone power systems for small islands and isolated grids. NaS batteries
have a high energy density, high efficiency of charge-discharge (89–92%), long
cycle life, and they are fabricated from inexpensive materials. However, be-
cause of the operating temperatures of 300 to 350◦ C and the highly corrosive
nature of sodium polysulfides, NaS battery technology is primarily suited for
large-scale, non-mobile applications such as grid energy storage.
8.7 Economics
Business entities always couch their concerns in terms of economics, such as
claiming that we cannot have a clean environment because it is uneconomical.
The thought here is that renewable energy is not economical in comparison
to coal, oil, and natural gas. The different types of economics to consider
are pecuniary, social, and physical. Pecuniary is what everybody thinks of as
economics, or money. On that note, we look at life-cycle costs rather than our
ordinary way of doing business, low initial costs. Life-cycle costs refer to all
costs over the lifetime of the system.
Social economics are those costs borne by everybody, and many businesses
want the general public to pay for their environmental costs. A good example
is the use of coal in China, where there are laws for clean air, but they are not
enforced. The cost will be paid in the future in terms of health problems. The
estimates of pollution costs for generation of electricity by coal range from
$0.005 to $0.1/[kW h].
Physical economics are the energy costs and the efficiency of the process.
There are fundamental limitations in nature due to physical laws. Energetics,
which is the energy input versus energy in the final product for any source,
should be positive. For example, production of ethanol from irrigated corn
has close to zero energetics. So, physical economics is the final arbitrator in
energy production and consumption. What each entity wants are subsidies for
itself and penalties for its competitors. Penalties come in the form of taxes
and environmental and other regulations, while incentives come in the form
of subsidies, breaks on taxes, lack of social costs to pay on the product, and
governmental funding of research and development. It is estimated that we use
energy sources in direct proportion to the incentives that source has received
in the past.
There are many examples of incentives and penalties for all types of energy
production and use.
8.8 Solutions
We do not have an energy crisis, since we know that energy cannot be created
or destroyed. We have an energy dilemma because of the finite amount of
readily available fossil fuels, which are our main energy source today. The
problem is the overpopulation and overconsumption. The world population is
so large that we are doing an uncontrolled experiment on the environment of
the earth. However, the developed countries were also major contributors to
this uncontrolled experiment in terms of consumption. The solution depends
on world, national, and local policies and what policies to implement, and
8.8 Solutions 333
KE = 0.5 mv 2 .
I = 0.5(M )(R)2 ,
I = 0.5(50)(3)2 [Kgm2 ],
and
E = 0.5(I)(ω)2 [J],
200 [rpm] ≈ 2094.395 [rad/s],
E = 0.5(I)(2094.395)2 [J].
Problem 8.4. Compare kinetic energy storage to batteries in terms of energy
density and installed costs. Consider both systems to have the same amount
of energy capacity.
60 [s]
E = 456 [kW ] × (30 [min] ) = 821 [M J],
[min]
E = 0.5(I)(ω)2 η = 0.5(0.5(1522.1)π(3.5)4 3)ω 2 (0.8),
E = 430.5 × 103 ω 2 .
So,
r
821 × 106
rad
ω= = 43.69 ,
430 × 103 s
ω = 417 [rpm].
Problem 8.10. Perform some brief research regarding the latest lithium-based
batteries. List their characteristics regarding energy density, power density,
charge/discharge regimes, and lifespan.
Solution 8.10. There is no precise solution to this problem.
8.11 Summary
The integration of different energy sources can be seen as complementary and
should commonly obey a local energy plan, as shown in Figures 8.2 and 8.3.
This means that not all the types of energy are suitable to being applied in
a particular place due to their natural restrictions and storage requirements.
Transmission lines and phase unbalancing also need to be considered so that
the whole integration fulfills energy usage objectives.
Some types of generators can be seen as complementary in very specific
senses. For example, an induction generator with a squirrel cage rotor will
consume reactive power that can be provided by a synchronous machine. A
capacitive bank can also be used to compensate the requested reactive power
but will add complexity to the transients of the generating device. Broadly,
the effects of unbalanced loads, faulty systems, and potential complementary
generators must be measured and controlled. In the same spirit, a smart grid
is required to manage complex networks, as it provides a means to deal with
full-scale integration of electric devices.
338 8 Renewable Energy
FIGURE 8.2
Integration of various sources 1/2
8.11 Summary 339
FIGURE 8.3
Integration of various sources 2/2
9
Smart Grid
The smart grid appeared as the response to human energy needs, and uses a
communication system, that allows improvement of electric energy conditions.
A new professional curriculum is a vital element for the future electrical
power industry. Different engineering topics need to be created to face the
smart grid’s technological obstacles. The development of educational programs
requires skills and technologies outside of power engineering.
The electric system curriculum requires theory and experimental work in
numerous areas, with a major focus on the following basic topics:
• Electric machines
• Power electronics
• Control systems
• Telecommunications
• Cybersecurity
• Automation
• Instrumentation
• Digital systems
• Economics
• Prediction
• Signal processing
Different authors have defined the smart grid system as [9]: Green Energy
Act (Canada): “A nickname for an ever widening palette of utility applications
that enhance and automate the monitoring and control of electrical distribu-
tion.” “An automated, widely distributed energy delivery network character-
ized by a two-way flow of electricity and information, capable of monitoring
and responding to changes in everything from power plants to customer pref-
erences to individual appliances.” “A smart grid is the electricity delivery
system (from point of generation to point of consumption) integrated with
communications and information technology.” In fact, the smart grid concept
341
342 9 Smart Grid
Utilities
What?
at? Receiver
How?
Why?
Electric utilities Renewable energy
Sender
Smart What?
Why?
Who? Grid
Renewable energy
Prosumers/consumers
Conventional electricity Who?
generation
When?
FIGURE 9.1
Supply chain for describing a smart grid
Computer
intelligence
Testbed Environment
and economics
FIGURE 9.2
Key groups of the smart grid technology
• Economic factors
– Optimized generator operation
– Deferred generation capacity investments
– Reduced ancillary service cost
– Reduced congestion cost
– Deferred transmission (T) capacity investments
– Deferred distribution (D) capacity investments
344 9 Smart Grid
• Reliability factors
– Reduced sustained outages
– Reduced major outages
– Reduced restoration cost
– Reduced momentary outages
– Reduced sags and swells
• Environmental
– Reduced CO2 emissions
– Reduced SOX, NOX, and PM-2.5 emissions
• Security
– Reduced oil usage
– Reduced wide-scale blackouts
Since the smart grid has shown its energy potential, it is important to
make a prediction about the cost of upgrading the traditional electric grid
presented in Figure 9.3 to a smart grid. As an example, the Electric Power
Research Institute forecast that the cost to upgrade the U.S. electric grid to a
smart grid could be between $338 billion and $476 billion, and could generate
benefits of between $1,294 billion and $2,028 billion; as result, the upgrade
of the electric grid is possible and beneficial for the utility system and end
user. Some of the key drawbacks when a traditional grid is used are: centralized
architecture and control, passive transmission and distribution, very extensive
network (long paths and many components), lack of diversity.
As it could be observed, the communication systems inside a smart grid are
very important and have to be restructured in order to get the communication
system working properly. Figure 9.4 shows how the smart grid communication
has to be adjusted.
The micro-grids are local and independent electrical grids that could in-
tegrate energy resources, transmission, distribution, and end-user. A local is-
land that can operate as a stand-alone or grid-connected system is defined as
a micro-grid. Normally, gas turbines or renewable resources power it and an
9.1 Benefits of the smart grid 345
Medium-voltage
Tower
Tow
Distribution line Veryy h
Ver high-voltage Power plant
Transm
Tra nsmission
nsm ion line
Transmission
Medium-voltage
Distribution line
Circuit
breaker
Distribution Transformer
Transformer
Low voltage
FIGURE 9.3
Conventional electric grid
FIGURE 9.4
Communication in a smart grid system
346 9 Smart Grid
Energy storage
Fuel cells
Source side energy
storage Power electronics
interface Loads
Fuel supply
Microturbines
Fuel supply
Power electronics
interface
FIGURE 9.5
Basic topology of a micro-grid
Smart Metering
Standards
Home Networking
Electric Vehicle
Standards
Standards
(Vehicular communications)
FIGURE 9.6
IEEE standards at home
energy usage in real-time, just as the old ones did. But, they also send infor-
mation back to the utility by wireless signals. Instead of using a utility meter
reader who goes to the property and manually does the periodical electric
service reading, the information is send automatically to the utility. To im-
prove the overall performance of the smart meters, a new group of appliances
has been developed. Smart appliances like the toaster, the microwave, and
even the refrigerator are based on micro-chips that allow your appliances to
communicate with the smart meters to control the energy usage.
Smart meters are part of a global plan to upgrade the world electrical grid
to a smart grid, which reduces energy consumption. The smart grid is being
implemented in both Canada and Mexico in order to communicate with the
United States. Designers are working on new standards that will integrate all
of North America into a single, unified smart grid system. There is an initiative
to generate a global smart grid to cover a continent, as presented in Figure 9.8.
The Global Energy Network Institute (GENI) shows this Dymaxion (tm) Map
of the world from the perspective of the North Pole, to reveal the global grid
currently under construction. The lines are high-voltage electrical transmission
lines that are capable of transferring large amounts of energy from continent
to continent.
348
FIGURE 9.7
Smart meter and smart home
Displays, interfaces and controller devices
Smart TV
Wireless display
WEB Enabled PC
Wireless Wireless
page Internet thermostat
water sensor
Meter/HAN 70°F
Services
Wireless
TDSP contact sensor
Meter/HAN
Transactions over utility
infrastructure
Wireless Wireless
Mobile Wireless
switch dimmer
Home theater controller devices remote
9 Smart Grid
9.1 Benefits of the smart grid 349
FIGURE 9.8
Global smart grid
350 9 Smart Grid
9.2 Cyber-security
The presence of information and communication technologies (ICT), e.g., com-
puters, networks, programs, and data, gives flexibility and efficiency to the
smart grid operation. However, the vulnerability of ICT to threats must be
considered [71].
Smart-grid cyber security focuses on protecting ICT from unintended or
unauthorized access, changes, or destruction. By default, all the components
in the smart grid have protection against intruders, but hackers are always
looking for the weakest and most vulnerable link that designers overlooked.
For example, in 2015 a couple of hackers attacked during a self-drive car test,
and were able to remotely modify the a/c and radio stereo of the car. In
fact, they could modify the acceleration, brakes, or steering, but these hackers
didn’t. The same attacks could be done on the smart grid elements, resulting
in a big risk for both power distribution and people.
Cyberattacks can be classified as [25]:
Malware: Malicious software designed by hackers to run in a specific system.
Viruses, worms, Trojan horses, ransomware, spyware, rootkit, etc. are ex-
amples of malware. Stuxnet is an example of malware affecting SCADA
and PLC systems.
Denial-of-Service or Distributed Denial-of-Service: Makes a networked
system busy and unavailable to its users by overloading the system with
a huge amount of requests.
Zero-Day: This is an unpatched computer system vulnerability that is not
yet disclosed to the public. The name refers to the vendor of the affected
computer systems having zero days to fix it after the vulnerability has
been disclosed.
Advance Persistent Threats: This is a generalized term for a group of
professional hackers targeting a system continuously over a long period of
time.
Cyberattacks are presented all the time, so a secure smart grid should have
enough protection to limit adversary access. Furthermore, if the cyber-security
is broken, the smart grid should contain elements that minimize or eliminate
the effect of the attacks.
Cybersecurity is one of the 20 metrics for measuring progress toward im-
plementation of smart-grid technologies, practices, and services, proposed by
the U.S Department of Energy in [8]. Six major characteristics of a smart grid
are proposed in [8]:
1. Enables informed participation by customers
2. Accommodates all generation and store options
9.3 Exercises 351
9.3 Exercises
9.3.1 Contribution of solar energy
Assume that there is demand for energy from a distant point and there is solar
energy to be exploited. In this exercise, the student will intervene by reducing
the consumption of energy from an older generation plant, using the surplus
energy produced by solar photovoltaic systems.
The reduction of even a minimum of absorbed energy will certainly have
an impact on environmental pollution produced by a plant of the older gen-
eration.
Set the resistive load to the second value and close the R2 relay to supply
energy coming from the coal plan (see Figure 9.9).
Close R4 relay to transfer energy coming from the plant to the load and
observe the power consumption in the maximum demand meter.
In this situation you can see active power required from the resistive load
352 9 Smart Grid
FIGURE 9.9
Solar energy global panel
FIGURE 9.10
Transformer panel
and a minimum of reactive power required from the primary of the step-down
transformer (see Figure 9.10).
Insert the hydroelectrical energy to the network using the hydroelectrical
plant simulator.
9.3 Exercises 353
Observe the brushless motor starting, and after that, enter in the procedure
for the automatic synchronization connection (see Figure 9.11).
FIGURE 9.11
Automatic synchronization connection
Increment the speed to 1800 [rpm] and then increment the excitation field
until the delta output of the generator goes up to 230 [V ].
Look at the AC display phase-to-phase voltage in the electrical power
digital measuring unit. Check, at the same time, either the secondary voltage
or frequency on the software; the level has to be equal to that of the network
400 [V ].
Look to the display of the generator synchronizing relay and reduce the
speed angle shift and the angle at under 20 degrees (see Figure 9.12).
When the hydroelectrical part is connected to the grid, you have to pay
attention to stabilizing the system, incrementing the active power consump-
tion from the load, and stabilize the system regulating the speed and exci-
tation. Check the LED near the turbine; it has to be green and stable (see
Figure 9.13).
Observe the reducing of power coming from the coal plant when you in-
troduce active power to the network.
354 9 Smart Grid
FIGURE 9.12
Display panel
FIGURE 9.13
Stabilizing the system
FIGURE 9.14
Transfer power panel
FIGURE 9.15
Starting the motor
the reduction of energy coming from old plant generators. Observe the meter
near the coal plat simulator.
9.4 Homework problems 357
FIGURE 9.16
Speed increment
FIGURE 9.17
Wind generator increment
Solution 9.1. There is no specific solution to this problem. The teacher would
ask the student to write about resources usage.
Problem 9.2. Visit the Smartrid.gov site and find at least 5 reported smart-
grid or micro-grid projects throughout the available documents. Report them,
taking into account the funding required and the expected outcomes.
Solution 9.2. There is no specific solution to this problem. However, students
can be asked to be as specific as desired. Projects at the proposed site are
thoroughly explained and presented, so their technical issues, economic factors,
and follow-up inquiries are also available.
Problem 9.3. The smart grid is mostly supported by one key concept: ad-
vanced metering infrastructure (AMI). Perform some brief research about such
systems and conclude about their pertinence toward actual smart-grid integra-
tion.
Solution 9.3. AMI is a framework on which the smart grid is to be supported.
It is composed of electronic devices that perform precise and real-time mea-
surements so the electric grid can be monitored and key decisions can become
possible during operation. In addition, the digital systems required to make
their information available at a supervisory level are also part of an AMI.
AMI goes beyond smart meters. Briefly, the whole infrastructure is mostly
oriented to data analysis and its translation into information.
Problem 9.4. Another big concept behind smart-grid integration is the energy
management system (EMS). Perform some brief research about such systems
and conclude about their pertinence toward actual smart-grid development.
Solution 9.4. EMS is not a new concept and is commonly related to SCADA
systems. As the smart grid can easily move beyond a simple human-machine
interface conceptually, the EMS is a concept that makes smart-grid definition
easier, as it states the foundations where the smart-grid is to be implemented.
Problem 9.5. The smart grid approach to energy management has not been
thoroughly developed. There are many issues precluding its final realization.
Perform some brief research and determine the top five issues that hinder
smart grid realization. Provide an explanation for each concept.
9.5 Simulation 359
• The electric grid is not available for testing, so R&D efforts have been
obscured
9.5 Simulation
It is important to consider that most assumptions made regarding system
interconnections depend on the infinite bus stiffness, which is considered to
be infinite indeed. This implies that it can give/take any amount of active
or reactive power while holding the voltage characteristics of the network.
Although this is known to be unrealistic, it will be accurate as long as the bus
surpasses other sources in terms of installed capacity.
In addition, the incorporation of new sources should normally follow from
local needs and must be calculated so their power capacity is as required. If
a generator is not fully used, its efficiency will drop. On the other hand, the
incorporation of many generators will decimate the efficiency of the infinite
bus as its generators will require less energy (see Figure 9.18).
Energy projects must then be preceded by exhaustive analysis before their
application. Undoubtedly, the incorporation of new “domestic power” sources
to the grid should be regulated and measured, so the already installed systems
and the new ones fit a known consumption profile.
360 9 Smart Grid
FIGURE 9.18
Power contribution of different power sources
10
Power Electronics in Power Systems Using
LabView, LabView-FPGA, and Multisim
10.1 Introduction
“Today, approximately 30 percent of all power generation utilizes power elec-
tronics between the point of generation and consumption. By 2030, it is ex-
pected that up to 80 percent of all generated electricity will utilize power elec-
tronics” [59]. Moreover, the power electronic stages have reduced their losses
(see Figure 10.1), so they are very attractive for controlling power systems.
Hence, it is very important to design power electronic stages using software
tools such as LabVIEW.
Not only are National Instruments customers completing more projects
on time (http://forums.ni.com/t5/Power-Electronics-Development/g
p-p/grp-1891), but they’re doing it in less time with fewer resources than
those who are following the traditional design approaches. This is made possi-
ble by using a platform-based approach that can fusion the best of build and
the best of buy, as well as leveraging high-level software with flexible off-the-
shelf hardware. These are integrated together, so that you get the opportunity
to focus on innovation, not your implementation. Graphical system design and
the LabVIEW reconfigurable inputs and outputs (RIO) architecture provide
quality products and the confidence to innovate, while getting your project
done in less time. The RIO relationship between price-performance and num-
ber of systems deployed is presented in Figure 10.2.
The two main goals of the power electronics design using simulation tools
are the following. A traditional power electronics design is done using a V
process that is shown in Figure 10.3. This process starts with co-simulation
that allows us to move to prototype in a fast manner because the design V
process is highly improved, such that co-simulation provides real results that
are much closer to the experimental design (see Figure 10.4).
Embedded systems are used in the V process for designing power electron-
ics systems. Thus, the availability of heterogeneous field programmable gate
array (FPGA) hardware contains an array of integrated digital signal proces-
sor (DSP) cores capable of efficiently executing floating-point math operations,
so that the co-simulation and the deployment process can be improved. On
the other hand, these modern FPGAs are actually hybrid devices containing
361
362 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.1
Generations of power electronics and power losses
10.1 Introduction 363
FIGURE 10.2
NI-RIO deployment curve
364 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.3
Design V process
10.2 Co-simulation LabVIEW and Multisim 365
FIGURE 10.4
Co-simulation, interface board, and commercial deployment
FIGURE 10.5
Automatic Multisim to LabVIEW FPGA conversion enables this power sys-
tem model to be simulated at 2 MHz speeds using heterogeneous FPGA hard-
ware
10.2 Co-simulation LabVIEW and Multisim 367
FIGURE 10.6
Power system with power electronics elements
etc.). Those components can either conduct or block the electric current in
the power system when a control signal is provided. They are very effective
for controlling power systems since the time to change their conductivity con-
dition is in the range of microseconds. Thus, a power electric system can be
controlled using power electronic stages. When the power electronic stages
are designed for power systems, it is important to use simulation tools that
help to design the power electronic stage in a precise manner and to save time
and economic resources. As a result, the faster way to move from simulations
to commercial deployments, and product generation, is the co-simulation in
which the subsystems are linked as a coupled problem that is simulated in a
distributed manner. One of the most important industrial and academic tools
is the co-simulation that is a closer representation of the real system in which
the power electronic stage can be analyzed in a comprehensive manner (see
Figure 10.6). This chapter deals with the LabVIEW FPGA and Multisim co-
simulation tools to understand the power electronics design and fundamentals
for power systems. In this chapter, you will gain a basic understanding of
how the LabVIEW graphical system design tool chain can be used to allow
for power electronics control. Figure 10.7 illustrates a block diagram about
designing power electronics using LabVIEW FPGA and Multisim as the first
step in the power electronic design [54]. This diagram can be determined by
the design V process (see Figure 10.3).
One of the most important advantages of using Multisim (Multisim Man-
ual) and LabVIEW in co-simulation is designing and analyzing a complete
power electronic system, including applications such as wind turbines. Thus,
co-simulation helps to effectively design across both the analog and digital
domains. Traditional platforms are unable to accurately simulate both analog
and digital together, so design errors cascade to the physical prototype and
lead to an ineffective and lengthy design approach. With the new co-simulation
368 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.7
Design block diagram for power electronics based on LabVIEW FPGA, Mul-
tisim, and LabVIEW
10.2 Co-simulation LabVIEW and Multisim 369
Include the
Defining the
LabVIEW co-
electric circuit in
simulation
Multisim
terminals
Include the
LabVIEW external Create the
model for LabVIEW
connecting program
LabVIEW
Run the
co-simulation
FIGURE 10.8
Basic block diagram for running a co-simulation
Step 1: Place the R-L-C electric circuit and the voltage source in the topology in
Multisim, as shown in Figure 10.9.
Step 2: This circuit will be able to send and receive data to and from the Lab-
VIEW simulation engine when the communication terminals are placed.
These terminals in Multisim are hierarchical block or sub-circuit (HB/SC)
terminals. Right-click and select Place on the schematic HB/SC connec-
tor from the shortcut menu, or simply type <Ctrl-I>. Place one HB/SC
connector above and to the left of the schematic and another connector
370 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.9
Components placed in Multisim
above and to the right of the schematic. Rotate the second connector 180
degrees by holding Ctrl and tapping R twice. Connect the HB connectors
to the circuit by selecting them from the context menu (see Figure 10.10.)
To configure the HB/SC connector as an input or output to LabVIEW, the
LabVIEW Co-simulation Terminals window must be opened. Nav-
igate to ViewLabVIEW Co-simulation Terminals. Notice that the
HB/SC connectors placed earlier appear in this window. To configure each
connector as an input or output, select the desired choice from the Mode
settings. Then configure each connector as a voltage or current input/out-
put by selecting the Type. Finally, if it has placed two input or output
terminals that are needed to function as a differential pair, a Negative
Connection can be selected. Configure IO1 as an Input and IO2 as an
Output. You can immediately see the shape of the connector changing
on the schematic based on whether it’s an input or an output connector
(see Figure 10.11). The final electric circuit in Multisim can be seen in
Figure 10.12.
Step 3: Creating a Digital Controller in LabVIEW is the last step for connecting
LabVIEW and Multisim. To pass data back and forth between LabVIEW
and Multisim, you must use a Control & Simulation Loop in LabVIEW.
Navigate to the block diagram in LabVIEW (white window), right-click
to open the Functions Palette, and navigate to Control Design &
SimulationSimulationControl & Simulation Loop. Left-click
and drag the loop onto the block diagram (see Figure 10.13). Now add
the Halt Simulation function to your VI to stop the Control & Simu-
lation Loop. Right-click the block diagram and navigate to Control
Design & Simulation Simulation Utilities Halt Simulation.
Left-click to place the Halt Simulation function on the block diagram and
then right-click the Boolean input of the Halt Simulation VI and select
Create Control. This creates a Boolean control on the front panel
that you can use to halt, or stop, the simulation while it is running. Next,
place the Multisim Design VI that facilitates the communication between
the LabVIEW and Multisim simulation engines. Right-click the block di-
10.2 Co-simulation LabVIEW and Multisim 371
FIGURE 10.10
HB connectors in Multisim
FIGURE 10.11
Input or output section
372 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.12
Electric circuit (R-L-C)
FIGURE 10.13
Selecting the simulation loop
10.2 Co-simulation LabVIEW and Multisim 373
FIGURE 10.14
Selecting the simulation block diagram in Multisim
(a)
(b)
FIGURE 10.15
Block diagram LabVIEW simulation (a) and frontal panel (b)
10.2 Co-simulation LabVIEW and Multisim 375
FIGURE 10.16
Buck converter designed in Multisim
FIGURE 10.17
Co-simulation diagram between LabVIEW and Multisim
the results when the controller is adjusting the PWM command and the
output voltage reaches the set point of voltage.
A Boos DC-DC converter can also be analyzed using co-simulation. The
topology of a Boost converter is shown in Figure 10.19, and the co-
simulation is presented in Figure 10.20. The input voltage Vin is lower
than the output voltage Vout . When Q1 is on, the diode (D2A) is reversed-
biased and the output stage is isolated. Hence, the inductor receives energy
from Vin . On the other hand, when Q1 is off, the output stage receives
energy from the inductor as well as Vin ; thus, the output voltage (Vout )
is increased. The general relationship between the duty ratio (d), Vin ,
and Vout is defined by Equation 10.1, when the DC-DC converter is in
continuous-conduction mode. When the PWM command is changed, the
output voltage is directly affected (see Figure 10.21).
Vin
Vout = . (10.1)
1−d
376 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.18
Co-simulation results using a conventional PID controller
FIGURE 10.19
DC-DC boost converter in Multisim
FIGURE 10.20
Co-simulation block diagram
10.2 Co-simulation LabVIEW and Multisim 377
(a)
(b)
FIGURE 10.21
Results of Boost converter for a PWM command of 0.89 (a) and 0.3 (b)
378 10 LabVIEW-FPGA-Multisim co-Simulation
Classical wind turbines are characterized by lower inertia than classical power
plants. In the case of back-to-back converter/inverter topologies, common with
permanent magnet synchronous machine (PMSM) type turbines, the interme-
diate DC-voltage bus creates an electrical decoupling between machine and
grid. This decoupling can lead to a reduced amount of stored kinetic energy
in the turbine, but the effect can be compensated for through suitable imple-
mentations of machine control. In the case of variable-speed wind generators,
[42] suggested using the blade and machine inertia as a kinetic energy storage
system to participate in regulation of primary frequency control. However,
releasing or storing kinetic energy can only be considered as a portion of the
primary control scheme. Since wind persistence is limited, the power reserve
cannot be guaranteed for the long term.
It is also suggested in [17] to maintain a power reserve with the help of the
pitch control when the wind generator works close to rated power. It is also
possible to obtain a power reserve with the help of generator torque control
but at a cost with regards to overall turbine power efficiency. In the case of
the doubly-fed induction machine (DFIG), the power captured by the wind
turbine is converted into electrical power by the induction generator and it
is transmitted to the grid by both the stator and the rotor windings. The
control system generates the pitch angle command and the voltage command
signals at the rotor in order to control the power of the wind turbine, the
DC bus voltage, and the voltage and frequency at the grid terminals. This
section shows an intelligent control system that is proposed for DFIG wind
turbines, which is intended to be cost effective and provide pitch control while
reducing the mechanical stresses. In the near future the power rating of the
wind turbines will increase further, especially in offshore applications. Doubly
fed induction generator machines can be designed to work over a reasonable
speed range for improving overall power quality and system efficiency.
An important factor in the design of variable-speed wind turbine controls
is the random excitation of the system caused by wind turbulence. Obtaining
an optimal control design that balances all of these requirements, constraints,
and challenges requires high-performance digital control systems.
Intelligent control (IC) systems have the ability to work well when the
operation point of the system is changing, because the IC system can imple-
ment its derived control laws from a number of different sources, including
human knowledge. Fuzzy logic systems, for example, are able to synthesize
information from human experts. Also, artificial neural networks can imple-
ment complex non-linear control laws that are automatically derived from the
desired input-to-output response mapping of the overall wind turbine system.
the stator windings. Since the generator is an asynchronous machine, the value
of the slip, s, defined as the ratio ωωrs between the slip and synchronous angular
frequency, is lower than 1 and, consequently, the rotor power, Pr , is only a
fraction of the stator power, Ps . The sign of Pr is a function of the slip sign
when the mechanical torque applied to the rotor is positive. The synchronous
speed, ωs , is normally a positive constant for generation with a constant grid
voltage. Thus, the slip power can flow in both directions: to the rotor from
the supply, and from the supply to the rotor. The speed of the machine can
be controlled from either the rotor- or stator-side converter in both super-
and sub-synchronous speed ranges. Thus, the machine can be controlled as a
generator or a motor in both super- and sub-synchronous operating modes,
yielding all four quadrants of operation [17]. Below the synchronous speed in
the motoring mode and above the synchronous speed in the generating mode,
a rotor-side converter operates as a rectifier and a stator-side converter as an
inverter, whereby slip power is returned to the stator. Below the synchronous
speed in the generating mode and above the synchronous speed in the mo-
toring mode, the rotor-side converter operates as an inverter and stator-side
converter as a rectifier, whereby slip power is supplied to the rotor. At the
synchronous speed, slip power is taken from the supply to excite the rotor
windings, and in this case the machine behaves as a synchronous machine.
The following conditions can occur in the induction generator:
1. Pr is positive for negative slip (speed greater than synchronous speed).
2. Pr is negative for positive slip (speed lower than synchronous speed).
3. Pr is transmitted to the DC bus capacitor and raises the DC voltage for
super-synchronous speed operation.
4. Pr is taken out of the DC bus capacitor and decreases the DC voltage for
sub-synchronous speed operation.
The grid-tied power converter, Cgrid , is used to generate or absorb the
power flowing to the grid, Pgc , in order to keep the DC voltage constant. In
steady-state for a theoretical lossless AC/DC/AC converter, Pgc is equal to Pr
and the speed of the wind turbine is determined by the power, Pr , absorbed or
generated by the rotor-tied converter, Crotor . The phase-sequence of the AC
voltage generated by Crotor is positive for sub-synchronous speed and negative
for super-synchronous speed. The frequency of this voltage is equal to the
product of the grid frequency and the absolute value of the slip. A beneficial
feature for utility grid operation is that Crotor and Cgrid have the capability
for generating or absorbing reactive power and can be used to control the
reactive power or the voltage at the grid terminals [37].
(a)
(b)
(c)
FIGURE 10.22
Wind turbine simulation results (a and b) and wind turbine model in Lab-
VIEW (d)
382 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.23
Wind turbine complete simulation model using LabVIEW block diagram
Instruments LabVIEW control design and simulation module [62]. This model
does not include the following effects:
• Magnetic saturation
• Airgap harmonics
• Skin effects
The dynamic equations for the induction machine in reference frame k are
defined below. The reference frame, k, can be stationary or rotating.
d k
Vsk ejθk = Isk ejθk Rs + jωk Ψsk ejθk + ejθk Ψ
di s
d k
Vrk = Irk Rr + j(ωk − ωr )Ψrk + Ψ
di r
Ψsk = Ls Isk + Lm Irk
Ψrk = Lr Irk + Lm Isk
dωr
Tm − TL = Jm
dt
Tm = (Ψsα isβ − Ψsβ isα )
faults locally within the turbine control system itself. Taking, for example, the
power converter of Figure 10.1, a possible fault condition is the disconnection
of one converter bridge arm. In this scenario, the fault status of the converter
can be divided into five basic types.
1. Primary fault condition 001: There are no faults detected in the electric
power electronics components. The converter is operating normally.
2. Primary fault condition 010: A fault is detected in only one converter.
This fault can be sub-classified into six types:
a. Fault detected in component 1 (001)
b. Fault detected in component 2 (010)
c. Fault detected in component 3 (011)
d. Fault detected in component 4 (100)
e. Fault detected in component 5 (101)
f. Fault detected in component 6 (110)
3. Primary fault condition 011: A fault is detected in the same phase of two
converters. These faults can be sub-classified into three types:
a. Fault detected in components 1 and 4 (001)
b. Fault detected in components 2 and 5 (010)
c. Fault detected in components 3 and 6 (011)
4. Primary fault condition 100: A fault is detected in two components of a
half-bridge. These faults can be sub-classified into six types:
a. Fault detected in components 1 and 3 (001)
b. Fault detected in components 2 and 4 (010)
c. Fault detected in components 3 and 5 (011)
d. Fault detected in components 4 and 6 (100)
e. Fault detected in components 5 and 1 (101)
f. Fault detected in components 6 and 2 (110)
5. Primary fault condition 101: A fault is detected in two crossing compo-
nents. These faults can be sub-classified into six types:
a. Fault detected in components 1 and 2 (001)
b. Fault detected in components 2 and 3 (010)
c. Fault detected in components 3 and 4 (011)
d. Fault detected in components 4 and 5 (100)
e. Fault detected in components 5 and 6 (101)
384 10 LabVIEW-FPGA-Multisim co-Simulation
TABLE 10.1
Tasks commonly solved by artificial neural networks (ANNs)
Task Description
Function approximation Linear and nonlinear functions can
be approximated by neural
networks. The ANN performs a
model or curve fitting function.
1) Data classification: The ANN
Classification
assigns data to a specific class or
subclass. Useful for identifying
patterns.
2) Signal classification: Time-series
data is classified into subsets or
classes. Useful for identifying the
signal source.
Unsupervised clustering Identifies order in data. Creates
clusters of data of unknown classes.
Forecasting Predict the next values in a
time-series waveform.
Control systems ANN commonly performs function
approximation, classification,
unsupervised clustering, and
forecasting. Also, ANNs are used in
model identification and analysis of
open and closed loop systems.
An 8-bit coding scheme is used to describe the fault whereby the eighth bit
distinguishes converter 1 from converter 2 (a value of 0 indicates converter 1),
the primary fault type is encoded from the seventh bit to the fifth bit, the
sub-classification type code is indicated with the fourth bit to the second bit,
and an odd or even checksum is located in the first bit.
TABLE 10.2
Hebbian learning procedure
Algorithm Hebbian learning procedure
Determine the input space.
Step 1 Specify the number of iterations
iterN um and initialize t = 0
Generate small random values of
weights wi
Step 2 Evaluate the Hebbian neural
network and obtain the outputs xi
Step 3 Apply the updating rule (4)
If t = iterN um then STOP
Step 4
Otherwise, go to Step 2
fault detection information. The weights stored in the ANN can only be rein-
forced if the input stimulus provides sufficient output response values. In this
way, weights change proportionally to just the output signals. By this fact,
neurons compete to be a dedicated reaction to a part of the input. Hebbian
neural networks are thus considered the first self-organizing nets. The learn-
ing procedure is based on the following statement pronounced by Hebb: As
A becomes more efficient at stimulating B during training, A sensitizes B to
its stimulus, and the weight on the connection from A to B increases during
training as B becomes sensitized to A.
Then, Steven Grossberg developed a mathematical model for this sentence
[63]:
new old
wAB = wAB + βxB xA ,
where wAB is the weight between the interaction of two neurons A and B, xi
is the output signal of the i neuron. xB xA is so called the Hebbian learning
term. The next algorithm gives the Hebbian learning procedure (Table 10.2).
Figure 10.24 shows a Hebbian network for a classification problem of 8-bit
input signal based on LabVIEW. In the same way, a key point signal can be
supplied in order to get its classification process about the fault diagnosis.
In fact, one of the main signals to observe in fault diagnosis is the direct
current pulse voltage Vd of the three phases commute circuit and look at the
output port of inductance load that contains the information about whether
the transistor is fault or not, according to the fault model. If the fault type is
different, the wave of Vd is different, and the fault type can only be analyzed by
the six spectrum quantity, containing the direct current a0 quantity of voltage,
the 1,2,3-harmonic amplitude value (A1 A2,A3), and the 1,2-harmonic phase
(P1, P2); as a result, Vd is the main point signal in the wave method. Since
the phase information is needed in fault diagnosis, the time π6 with phase A
voltage passing the zero acts as a reference point at every diagnosis. For the
386 10 LabVIEW-FPGA-Multisim co-Simulation
(a)
(b)
FIGURE 10.24
(a) Hebbian network 8-bit classification and classification using generated fault
signals (b)
10.3 Case study: Intelligent wind turbine control with LabVIEW 387
Get data of Vd and Va
Get RMS of
phase voltage
FIGURE 10.25
Fault diagnosis process using neural networks
FIGURE 10.26
Results using Hebbian artificial neural network
wave method, the N pieces of unitary data collected in a period are directly
input into the neural networks. A simple flow chart of the fault diagnostic
process is presented in Figure 10.25.
The results obtained using Hebbian artificial neural networks when a fault
signal was generated are presented in Figure 10.26. The angle β acts like a
sampler at fixed intervals (8◦ ), the result presented was getting the angle β =
45◦ for a big-type fault. The results show a good classification and recognition
system. So, the use of neural networks for detecting all the fault signals of the
converter presented below gets zero error.
Using a control topology based on intelligent control systems, the perfor-
mance of the wind turbine can be improved. The proposed control system [62]
includes a fuzzy logic controller for the pitch angle, a grey system for predict-
ing the wind profile, and two ANFIS for controlling the generator currents.
Figure 10.27 shows the proposed topology in which the fuzzy logic controller
is a PI fuzzy logic controller [63]). The ANFIS was tuned with constant out-
puts; for tuning the ANFIS controller, the hybrid training was applied. The
simulation results are presented in Figures 10.27, 10.28, and 10.29.
The above results were obtained without the power electronics stage (in-
verter) that is one of the main stages in the simulation of a wind generator.
388 10 LabVIEW-FPGA-Multisim co-Simulation
(a)
(b)
(c)
FIGURE 10.27
(a) Main control and simulation loop, (b) model estimation process, and (c)
detailed view of the fuzzy controller
10.3 Case study: Intelligent wind turbine control with LabVIEW 389
FIGURE 10.28
Wind speed controller
390 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.29
ANFIS incorporation to rotor voltage computation
10.3 Case study: Intelligent wind turbine control with LabVIEW 391
FIGURE 10.30
Frontal panel of grey predictor results (real wind, grey predictor, and predic-
tion error)
Thus, the following test was done, including the power electronic stage, using
a space vector modulation and ideal switches. The parameters of the controller
were kept constant, using the previous tuning (without the inverter). As the
results show, the controllers can deal with the pulse distortion generated for
the power electronic stage. The selection of SPWM helps to decrease the har-
monics and also increase the DC bus energy that can be used. Figure 10.31
presents the rotor voltage using the space pulse width modulation (SPWM),
and Figure 10.32 shows the performance of the wind turbine using intelligent
control systems.
(a) (b)
FIGURE 10.31
Rotor voltage: (a) reference signals and rotor voltage and (b) three-rotor volt-
age using the SPWM
392 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.32
Controller performance including the power electronic stage using space pulse
width modulation
FIGURE 10.33
Reducing the simulation to deployment cycle from weeks to hours
394 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.34
Co-simulation between LabVIEW FPGA and Multisim
10.4 Co-simulation LabVIEW FPGA and Multisim 395
FIGURE 10.35
Co-Simulation flow diagram
396 10 LabVIEW-FPGA-Multisim co-Simulation
(a)
(b)
FIGURE 10.36
Induction motor drive (inverter A) and grid synchronized inverter (inverter
B), frontal panel (a), and Multisim diagram (b)
10.5 Case study: A half-bridge inverter control 397
front panel and block diagram, highlighting the differences from LabVIEW
programming on Windows (see Figure 10.37). Follow the next basic steps in
order to create and connect the co-simulation. Install the required software
using the Measurements and Automation Explorer (see Figure 10.38) and
connect the hardware (see Figure 10.39). Finally, the co-simulation file using
LabVIEW FPGA can be created (see Figure 10.40).
1. Disable WiFi
2. Launch NI MAX
3. Refresh remote systems
4. Format disk
5. Install firmware: LabVIEW RT Add-Ons
6. LabVIEW PID and Fuzzy Logic Toolkit
7. NI system configuration
8. System State Publisher / Network I/O
9. Network Streams/ Protocols and Buses
10. NI-Watchdog
11. Set LV Project IP address
12. Connect
10.5.1 Background
The field-programmable gate array (FPGA) allows a user to perform a variety
of tasks, including signal processing, at higher speeds than most processors.
The FPGA desktop execution node allows a user to run an FPGA appli-
cation on the development computer with simulated I/O. By specifying the
desired FPGA clock speed, as well as the number of ticks before inputs are
read and outputs are written, a user can determine the functionality of their
FPGA application without needing to compile and deploy the application to
398 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.37
Project explorer
FIGURE 10.38
Measurements and automation explorer
10.5 Case study: A half-bridge inverter control 399
FIGURE 10.39
Connect the GPIC inverter board
400 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.40
FPGA — Multisim co-simulation for manual IGBT control using LabVIEW
FPGA and Multisim
10.5 Case study: A half-bridge inverter control 401
the FPGA hardware [54]. Figure 10.41 shows the FPGA co-simulation. To ef-
fectively simulate an FPGA design using the FPGA desktop execution node,
you must understand two time paradigms, real-world time and simulated time.
Real-world time is the physical amount of time that elapses while something
occurs. Because an FPGA is a programmable circuit, it takes a fixed amount
of real-world time to execute such a circuit. To simulate a design, you can
create a model that reflects the functional behavior of the FPGA and execute
this model on a computer processor. However, such a model does not have the
same real-world timing as the FPGA. To create a model that has the same
real-world timing as the FPGA, you must use simulated time. Simulated time
is an event-driven model of real-world time. When LabVIEW executes block
diagram nodes, the execution causes simulated time to advance a certain num-
ber of steps. Clock ticks are the unit of simulated time where one clock tick
is representative of a single cycle of the referenced FPGA clock [6].
For this exercise, you will need to use the LabVIEW FPGA in order to
control the power inverter. A power inverter consists of a voltage source, two
insulated-gate bipolar transistors (IGBTS), an inductor, a capacitor, and a
resistor (see Figure 10.42). The IGBTs switch in order to change the current
flowing through the resistive load while the inductor and capacitor filter and
smooth the output signal. The power inverter performs DC-to-AC conversion.
10.5.2 Procedure
To run the half-bridge inverter control, a Multisim file and two TestBench
VIs are needed. As shown in Figure 10.43, a project was designed. On the
other hand, the Multisim file is used for studying the thermal conditions in
the semiconductors, as shown in Figure 10.44. This Multisim file has to be
constructed in order to create the co-simulation. In the Multisim editor envi-
ronment, if you do not see the LabVIEW Co-Simulation Terminals pane
(red box in the screenshot above), go to the View menu, select it, and then
position it in the desired location (see Figure 10.45).
It is necessary to define a LabVIEW project. This file is the front panel of
the half-bridge IGBT controller. In the LabVIEW project (see Figure 10.46),
double-click on [TestBench] Half-Bridge IGBT Control - Starting
Point.vi.
The goal of this example is to control the IGBTs (switches) in such a
way as to avoid damaging the hardware. In addition to the switches, other
controls include a slider for the modulation index and a button for stopping
the program. Important indicators include the IGBT temperatures, a dial for
the output voltage, both instant and RMS, as well as a chart that plots various
waveforms from the circuit. Open the block diagram by pressing CTRL+E or
by navigating to WindowShow Block Diagram and observe the layout.
Press CTRL+H to open the Context Help window and place the mouse
over the blue box that reads Desktop Execution Node. This will cause the
Context Help window to display information about the desktop execution
402 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.41
To enable high fidelity co-simulation, the LabVIEW FPGA desktop execution
node interfaces to the LabVIEW FPGA application while the Multisim co-
simulation interface connects to the power electronics circuit model. A desktop
testbench application is shown.
10.5 Case study: A half-bridge inverter control 403
(a)
(b)
FIGURE 10.42
The Multisim power electronics circuit (a) and LabVIEW FPGA control ap-
plication (b) are co-simulated with a common sense of simulated time
404 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.43
The project design for a half-bridge IGBT
FIGURE 10.44
Thermal analysis of IGBTs using Multisim
10.5 Case study: A half-bridge inverter control 405
FIGURE 10.45
LabVIEW co-simulation terminals
406 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.46
LabVIEW half-bridge IGBT control file
node (see Figure 10.47). The desktop execution node allows users to run FPGA
VIs on the development computer with simulation I/O.
To run the program and experiment with the IGBT switches, navigate to
the front panel of [TestBench] Half-Bridge IGBT Control - Starting
Point.vi. Run the application by pressing the Run button on the toolbar.
Check the box Close on successful completion. Select Close once the
deployment is complete (see Figure 10.48.)
The application will take a moment to begin running. Once running, the
waveforms will begin to be plotted on the graph. While the program is run-
ning, open the front panel to [FPGA] Half-Bridge Inverter Control.vi.
Left-click the top switch to close it. Observe the graph on the front panel of
the [TestBench] Half-Bridge IGBT Control - Starting Point.vi. (See
Figure 10.49.)
It is very important to understand that co-simulation using FPGA can
gives results closer to the experimental ones, as shown in Figure 10.50, in
which experimental and simulation results are compared.
FIGURE 10.47
Desktop execution node
FIGURE 10.48
Deployment process
408 10 LabVIEW-FPGA-Multisim co-Simulation
(a)
(b)
FIGURE 10.49
Frontal panel (a) and block diagram (b) of the co-simulation
10.6 Case study: the phase-lock loop (PLL) 409
FIGURE 10.50
Comparing simulated against experimental results
10.6.2 Background
Phase-locked loops are a form of feedback control. A PLL produces an oscil-
lating output, which begins operating with a set frequency. An input is then
passed to the PLL to serve as a reference point. The PLL adjusts the output to
then closely match the input phase. The longer the PLL runs, the more closely
the PLL output is able to track the input, given that the input is stable. PLLs
are useful in a variety of synchronous applications. Such applications include
clock multipliers, demodulation, and power-grid phase locking. In order to
connect an inverter to the grid, the inverter output’s phase must be within a
few degrees of the grid signal phase. This is where the PLL will be applied.
By using the grid signal as a reference, the PLL is capable of generating a
waveform that will be phase-locked to the grid signal. This output can then
be used to drive PWM in order to generate the inverter output, as shown in
Figure 10.51.
For creating a three-phase PLL, a project has to be created. There are a
Multisim file and three VIs in the project, as shown in Figure 10.52.
The VI includes a front panel of the PLL application in which a vertical
410 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.51
Frontal panel of a three-phase PLL
FIGURE 10.52
PLL files in the project
10.6 Case study: the phase-lock loop (PLL) 411
slider on the left to control frequency of the simulated power grid waveforms
is added, as shown in Figure 10.53.
The Desktop Execution Node to open up the FPGA application
[FPGA] PLL Testbench.vi is presented in Figure 10.54. For this appli-
cation, the Clock Ticks setting is set to 325 ticks. Thus, for every 325 ticks
of the 40 MHz FPGA clock, every 8.125 microseconds of simulated time, the
Desktop Execution Node will return the register values for the items chosen
in the Selected Resources dialogue (see Figure 10.55). Note that the FPGA
resources you can choose from include both front panel controls and indica-
tors as well as I/O resources. In this case, no changes to the configuration are
necessary, so click Cancel to exit.
When the TestBench application is run, the PLL frequency is initially 30
Hz and takes a half cycle to begin tracking the noisy simulated grid waveforms.
You can see the frequency of PLL in the indicator labeled PLL Frequency.
Due to the noisy nature of the 3-phase grid waveforms and the aggressive tun-
ing of the PLL proportional-integral (PI) control system, the PLL Frequency
oscillates around the correct grid frequency (see Figure 10.56).
As it was presented, the cos-simulation is a powerful tool for designing
power electronics in power systems which has to be considered during the
creating process in power electronics. As a result, extremely complex systems
can be created, such as the micro-grid experimental laboratory shown in Fig-
ure 10.57.
412 10 LabVIEW-FPGA-Multisim co-Simulation
(a)
(b)
FIGURE 10.53
Front panel (a) and block diagram (b) for the PLL co-simulation
10.6 Case study: the phase-lock loop (PLL) 413
FIGURE 10.54
PLL control loop
FIGURE 10.55
Configure desktop execution node
414 10 LabVIEW-FPGA-Multisim co-Simulation
FIGURE 10.56
PLL Co-simulation results
FIGURE 10.57
Micro-grid experimental laboratory [54]
Bibliography
[1] GE: GE wind turbine technology selected for new wind power station
in Japan. http://site.ge-energy.com/about/press/en/2003_press/
100103.htm. (Accessed 05/19/2017).
[2] NI multisim user manual — National Instruments. http://www.ni.com
/pdf/manuals/374483d.pdf. (Accessed 05/19/2017).
[3] Simulation fundamentals: Cosimulation in NI multisim — National In-
struments. http://www.ni.com/white-paper/6369/en/. (Accessed
05/19/2017).
[4] Stardom and fa-m3 ensure smooth supply of power to grid by wind farm
equipped with large-capacity NAS batteries — Yokogawa Electric Cor-
poration. https://www.yokogawa.com/library/resources/referen
ces/stardom-and-fa-m3-ensure-smooth-supply-of-power-to-gri
d-by-wind-farm-equipped-with-large-capacity-nas-batteries/.
(Accessed 05/19/2017).
[5] Storage boosts the power of renewable energy — Renewable Energy
World. http://www.renewableenergyworld.com/articles/2008/
06/storage-boosts-the-power-of-renewable-energy-52716.html.
(Accessed 05/19/2017).
[6] Using the LabVIEW FPGA desktop execution node — National In-
struments. http://www.ni.com/white-paper/51859/en/. (Accessed
05/19/2017).
415
416 Bibliography
[53] Weedy B. M., Cory B. J., N. Jenkins, Ekanayake J. B., and G. Strbac.
Electric power systems. John Wiley & Sons, 2012.
[67] Maurice George Say. Alternating current machines. Halsted Press, 1976.
[68] P. C. Sen. Principles of electric machines and power electronics. John
Wiley & Sons, New York, 2nd edition, 1996.
[69] Marcelo Godoy Simoes, Robin Roche, Elias Kyriakides, Sid Surya-
narayanan, Benjamin Blunier, Kerry D. McBee, Phuong H. Nguyen,
Paulo F. Ribeiro, and Abdellatif Miraoui. A comparison of smart grid
technologies and progresses in Europe and the U.S. IEEE transacitions
on industrial electronics, 48(4):1154 – 1162, 2012.
[70] Susan Solomon. Climate change 2007—the physical science basis: Work-
ing group I contribution to the fourth assessment report of the IPCC,
volume 4. Cambridge University Press, 2007.
[71] Siddharth Sridhar, Adam Hahn, and Manimaran Govindarasu. Cy-
ber–physical system security for the electric power grid. In Proceedings
of the IEEE, volume 100. IEEE, 2012.
[72] R. Strzelecki and G. Benysek. Power Electronics in Smart Electrical
Energy Networks, volume 1 of Power Systems. Springer, London, 1st
edition, 2008.
[73] Aman Abdulla Tanvir, Adel Merabet, and Rachid Beguenane. Real-time
control of active and reactive power for doubly fed induction genera-
tor (DFIG)-based wind energy conversion system. Energies, 8(9):10389–
10408, 2015.
[74] S Tegen, E Lantz, M Hand, B Maples, A Smith, and P Schwabe. Cost of
wind energy review, National Renewable Energy Laboratory. Technical
report, Technical Report, 2011.
[75] DM Triezenberg. Electric power systems. Perdue University: West
Lafayette, IN, 1978.
[76] Juan Carlos Vega de Kuyper and Santiago Ramı́rez Morales. Fuentes de
energı́a, renovables y no renovables. Aplicaciones. 2014.
[77] Blume S. W. Electric power system basics: for the nonelectrical profes-
sional. Wiley-Interscience; IEEE Press, 2007.
[78] C. L. Wadhwa. Electrical power systems. New Age International, 2005.
[79] G.J. Wakileh. Power systems harmonics: fundamentals, analysis, and
filter design. Springer, New York, 2001.
[80] Hayt W.H. Engineering circuit analysis. McGraw-Hill Higher Education,
Boston, 2006.
[81] David C White and Herbert H Woodson. Electromechanical energy con-
version. Wiley, 1959.
Index
A Apparent power, 15
AAAC, see All aluminum-alloy Artificial neural networks (ANNs),
conductor 384
AAC, see All-aluminum conductor Asynchronous machine, 94–99
ABCD model (transmission line),
233–235 B
ACAR, see Aluminum conductor Back-to-back (B2B) converters, 165,
alloy reinforced 166
AC power analysis, 12–22 Batteries, 329
apparent power and power Blocking diode (solar energy), 196,
factor, 15–16 197
average power, 14 Boost DC-DC converter, 375
complex power, 16–17 Brushless motor, 100–104
effective values or RMS values,
14–15 C
exercises, 19–22 Cable lines, 58
instantaneous power, 12–14 Compressed air energy storage
power factor correction, 17–19 (CAES), 328
ACSR, see Aluminum-conductor Counterclockwise (CCW) angle
steel-reinforced cable rotation, 60
Advanced metering infrastructure Current transformers (CT), 59
(AMI), 358 Cybersecurity, 350–351
Advance persistent threats, 350
Airfoils, 155–158 D
All aluminum-alloy conductor DC flywheels, see Direct current
(AAAC), 211 flywheels
All-aluminum conductor (AAC), 211 Denial-of-service attack, 350
Aluminum conductor alloy reinforced DFIG, see Doubly-fed induction
(ACAR), 211 generator
Aluminum-conductor steel-reinforced Digital signal processor (DSP), 361
cable (ACSR), 57, 211 Direct current (DC) flywheels,
AMI, see Advanced metering 328–329
infrastructure Displacement power factor (dPF), 39
Ampere’s law, 214 Distortion power factor (hPF), 39
Amplitude modulation (AM), 312 Diversion hydroelectric power plants,
Angle of attack, 115 124
Animal products, 324 Double line to ground fault, 277–278
421
422 Index