ECE 530 - Analysis Techniques For Large-Scale Electrical Systems
ECE 530 - Analysis Techniques For Large-Scale Electrical Systems
1
Numerical Integration
• Power system application: transient stability
• Power systems are dynamic state-space systems,
represented by DAEs and ODEs
• Typically the power flow solution represents an
equilibrium point
• At some point a contingency occurs, perturbing the
state away from the steady-state equilibrium point
• Time domain simulation is used to determine whether
the system returns to the equilibrium point
2
Methods
• One-step methods
– Forward Euler’s Method
– Runge-Kutta Methods
• Multistep methods
– Adams-Bashforth Method
• Multi-Rate Methods
• All these methods have numerical stability concerns
depending on the stepsize
• Implicit Methods
3
Multistep Methods
• Euler's and Runge-Kutta methods are single step
approaches, in that they only use information at x(t) to
determine its value at the next time step
• Multistep methods take advantage of the fact that using
we have information about previous time steps x(t-Dt),
x(t-2Dt), etc
• These methods can be explicit or implicit [dependent on
x(t+Dt) values]; we'll just consider the explicit Adams-
Bashforth approach
4
Multistep Motivation
• In determining x(t+Dt) we could use a Taylor series
expansion about x(t)
Dt 2
x(t Dt ) x(t ) Dt x(t ) x(t ) O(Dt 3 )
2
Dt 2 f (x(t )) f (x(t Dt ))
x(t Dt ) x(t ) Dt f (t ) O(Dt )
2 Dt
x(t Dt ) x(t ) Dt f x(t ) f x(t Dt ) O(Dt 3 )
3 1
2 2
6
Adams-Bashforth Versus
Runge-Kutta
• The key Adams-Bashforth advantage is the approach
only requires one function evaluation per time step
while the RK methods require multiple evaluations
• A key disadvantage is when discontinuities are
encountered, such as with limit violations;
• Another method needs to be used until there are
sufficient past solutions
• They also have difficulties if variable time steps are
used
7
Numerical Instability
• All explicit methods can suffer from numerical
instability if the time step is not correctly chosen for the
problem eigenvalues
10
Multi-Rate Methods
• At each macro step the slow variables are integrated, at
each micro step the fast variables are integrated
• Macro variables can be interpolated during the micro
steps
Source: J. Chen and M. L. Crow, "A Variable Partitioning Strategy for the Multirate Method in Power
Systems," IEEE Trans. Power Systems, vol. 23, pp. 259-266, 2008
11
Multi-Rate Example: Transient Stability
• The power system transient stability problem is usually
solved with a time step of ¼ or ½ cycle
• Some subsystems can have much faster time constants
– When starting induction machines can exhibit very fast
(relative to the time step) transients
– Some types of exciters can have very fast time constants, in
which the dynamics only come into play during close by
faults
12
Implicit Methods
• Implicit solution methods have the advantage of being
numerically stable over the entire left half plane
• Only methods considered here are the is the Backward
Euler and Trapezoidal
Euler
x(t Dt ) x(t ) Dt f (x(t ))
Backward Euler
x(t Dt ) x(t ) Dt f (x(t ) Dt )
Trapezoidal
Dt
x(t Dt ) x(t ) f (x(t )) f (x(t ) Dt )
2 13
Implicit Methods
• The obvious difficulty associated with these methods is
x(t) appears on both sides of the equation
• Easiest to show the solution for the linear case:
x f (x(t )) Ax(t ))
Then using backward Euler
x(t Dt ) x(t ) Dt A(x(t Dt ))
I Dt A x(t Dt ) x(t )
1
x(t Dt ) I Dt A x(t )
14
Backward Euler Example
• Returning to the cart example
0 1
x x(t ))
1 0
Then using backward Euler with Dt 0.25
1
1 1 0.25
x(t Dt ) I Dt A x(t ) x(t )
0.25 1
15
Backward Euler Example
• Results with Dt = 0.25 and 0.05
time actual x1(t) with x1(t) with
x1(t) Dt=0.25 Dt=0.05
0 1 1 1
0.25 0.9689 0.9411 0.9629
0.50 0.8776 0.8304 0.8700
0.75 0.7317 0.6774 0.7185
1.00 0.5403 0.4935 0.5277
2.00 -0.416 -0.298 -0.3944
Note: Just because the method is numerically stable
doesn't mean it is doesn't have errors! RK2 is more
accurate than backward Euler for this case.
16
Trapezoidal Linear Case
• For the trapezoidal with a linear system we have
x f (x(t )) Ax(t ))
Dt
x(t Dt ) x(t ) A(x(t )) A(x(t Dt ))
2
I Dt A x(t Dt ) I Dt A x(t )
2 2
1 Dt
x(t Dt ) I Dt A I A x(t )
2
17
Trapezoidal Example
• Results with Dt = 0.25, comparing between backward
Euler and trapezoidal
time actual Backward Trapezoidal
x1(t) Euler
0 1 1 1
0.25 0.9689 0.9411 0.9692
0.50 0.8776 0.8304 0.8788
0.75 0.7317 0.6774 0.7343
1.00 0.5403 0.4935 0.5446
2.00 -0.416 -0.298 -0.4067
18
Overview of Commercial Transient
Stability Algorithm
• Most commercial packages use an explicit approach,
such as a second order Runge-Kutta because of 1) the
large number of nonlinear models and 2) the number of
limit violations that are encountered
• Power flow solution provides the initial starting point
• Initial state variables are determined by back solving
from the power flow conditions
– Example: TGOV1 governor
initialization
19
Overview of Commercial Transient
Stability Algorithm
• For simplicity showing Euler's, which is not used
• While (t <= tend) Do Begin
– Any events? If so, apply event and solve algebraic
equations, g(x,y) = 0
– Solve differential equations: x f (x, y )
– Solve algebraic equations, g(x,y) = 0
– Output results
– End while
• Many complications to consider, including events that
occur in the middle of a time step
20
Power System Equivalents
• For many power system applications it is not necessary
to study the entire interconnected network
– Usually we are only concerned with a portion of the network
– For real-time operations, real-time information is only
available for a portion of the network
• System is partitioned into study system, for which a
detailed model is desired, and an external system, for
which an equivalent model is used
– Boundary buses (within the study system) connect the two
21
Power System Equivalents
• For decades power system network models have been
equivalenced using the approach originally presented by
J.B. Ward in 1949 AIEE paper “Equivalent Circuits for
Power-Flow Studies”
– Paper’s single reference is to 1939 book by Gabriel
Kron, so this also known as Kron’s reduction
• Additional classical techniques are discussed in S.
Deckmann, A. Pizzolante, A. Monticelli, B. Stott, and
O. Alsac,“Studies on power system load flow
equivalencing,” IEEE Trans. Power App. Syst., vol.
PAS-99, no. 6, pp. 2301–2310, Nov./Dec. 1980.
22
Ward Equivalents (Kron Reduction)
• Equivalent is performed by doing a reduction of the bus
admittance matrix
I e Yee Yes Ve
I Y Yss Vs
s se
I s Y Y 1
se ee eI Yss Yse ee Yes Vs
Y 1
0 MW
MVA
MVA
A
Case Hourly Cost
3454 $/h
A
MVA AGC OFF system is
0 MW
0 MW
1.04 pu
0 MW
0 MW
A
MVA
0 MW 1.04 pu
modeled
Bus 2 Top Area Cost
0 MW
0 MW
2122 $/h
Bus 5
0 MW 0 MW
unloaded
A
A 0 MW AGC OFF
MVA
MVA 0 MW
0 MW 0 MW 0 MW
A
26
PowerWorld B7Flat_Eqv Example
• Original Ybus
20.83 16.67 4.17 0 0 0 0
16.67 52.78 5.56 5.56 8 .33 16 .67 0
4.17 5.56 43.1 33.3 0 0 0
Ybus j 0 5.56 33.3 43.1 4.17 0 0
0 8.33 0 4.17 29.17 0 16.67
0 16.67 0 0 0 25 8.33
0 25
0 0 0 16.67 8.33
20.833 4.167 0
Yee j 4.167 43.056 33.333
0 33.333 43.056
27
PowerWorld B7Flat_Eqv Example
16.667 5.556 5.556
16.667 0 0 0 0
Yes j 5.556 0 0
0 4 .167
0 Yse j
0 0 0
5.556 4.167 0 0
0 0 0
31