SUPPLEMENTARY INFORMATION
|
|
- Diana Cunningham
- 6 years ago
- Views:
Transcription
1 DOI: 1.138/NPHYS535 Spontaneous synchrony in power-grid networks Adilson E. Motter, Seth A. Myers, Marian Anghel and Takashi Nishikawa Supplementary Sections S1. Power-grid data. The data required for power flow calculations were obtained as follows. For the 1-generator system, known as the New England test system, the parameters were taken from ref. 47. For the 3- and 5-generator systems, the parameters were taken from refs. 3 and 48, respectively. The data for the Guatemala and Northern Italy systems were provided by F. Milano (University of Castilla La Mancha), and the data for the Poland system were provided as part of the MATPOWER software package 49. The dynamic data required to calculate the synchronous state, determine its stability, and simulate equation () in the main text, are not all available for the real power grids, and were obtained as follows. For all systems, we assumed that before any optimization the damping coefficient and droop parameter satisfy D i +1/R i = 5 per unit for all generators. The parameters x d,i and H i are available for the three test systems from the respective references mentioned above. For the Guatemala, Northern Italy, and Poland systems, we estimated x d,i and H i using the strong correlation observed in the test systems between each of these parameters and the power P i injected by generator i into the network, as shown in Supplementary Fig. S1. The estimated values are x d,i in per unit, and H i is in seconds. 9.8P 1.3 i and H i.4p i, where P i is in megawatts, x d,i is S. Derivation of the swing equation. Equation () in the main text can be derived by setting the rate of change of the angular momentum of the rotor equal to the net torque acting on the rotor: J d δ i dt = T mi T ei, (S1) where J is the moment of inertia in kg m, T mi is the mechanical torque in N m accelerating the rotor, and T ei is the typically decelerating torque in N m due to electrical load in the network. Multiplying both sides by ω i and using the fact that the torque in N m multiplied by the angular velocity in radians per second gives the power in watts, the equation can be written in terms of power: d δ i Jω i dt = P mi P ei. (S) To make P mi and P ei per unit quantities, we divide both sides of the equation by the rated power P R (used as a reference). The factor Jω i then becomes H i /ω i, where we defined the inertia constant 1 NATURE PHYSICS 1
2 DOI: 1.138/NPHYS535 H i = W i /P R (in seconds) and the kinetic energy of the rotor W i = 1 Jω i (in joules). Noting that ω i is approximately equal to the reference frequency ω R in systems close to synchronization, we obtain equation () in the main text. A more detailed description of this derivation can be found in ref. 3 (second edition, pp ). S3. Relation to the master stability formalism. Equations (4) and (5) in the main text, when written individually for each node i, can be expressed as d δ i dt = 1 δ i n + δ j P ij. (S3) β j=1 i 1 δ i δ i This is in the same general form as the variational equation for the class of coupled oscillators considered in ref. 34, which for a synchronous state s = s(t) is δ j d dt x i = DF(s) x i + σ j G ij DH(s) x j, (S4) where ẋ i = F(x i ) describes the node dynamics, H is the coupling function, and σg ij represents the strength of coupling from node j to i. The factors DF(s) and DH(s) in equation (S4) are both constant matrices in equation (S3), and P ij in equation (S3) corresponds to σg ij in equation (S4), which are relations that can be used to derive Λ β (α) from equations (4) and (5) based on the results in ref. 34. Therefore, while the formalism in ref. 34 cannot be directly applied to equation () in the main text, it can be applied near the synchronization manifold, which is a procedure previously proposed for a broad class of coupled oscillators 35. S4. Local analysis of λ max (β opt + εβ ) along all directions β. For a fixed direction β = (β 1,...,β n) T in the space of all β i values and β opt =(β opt,...,β opt ) T, we consider β = β opt + εβ and study the behavior of the maximum Lyapunov exponent λ max as a function of ε in the limit of small ε. The same transformation that leads to equation (7) in the main text, now applied to the general case of nonidentical β i values, results in Ż = MZ, where 1 M =, Z = Z 1, J B and B = Q 1 BQ. Here and 1 denote the null and the identity matrix, respectively, and we recall that J, Z 1, Z, B, and Q are defined in the main text. We denote the eigenvalues of M as Z NATURE PHYSICS
3 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION ν j± = ν j± (ε), j =1,...,n, and index them so that, as ε, each ν j± approaches continuously the corresponding Lyapunov exponent given by equation (8) in the main text for the case β i = β opt, i.e., ν j± (ε) ε ν j± () = λ j± (α j,β opt )= β opt ± 1 βopt 4α j = α ± α α j. Here, as in the main text, the eigenvalues of the matrix P are denoted α 1,...,α n, among which α 1 =and α is the smallest nonzero eigenvalue. Also note that ν 1± (S5) are not relevant for determining the stability of synchronous states, since these eigenvalues are associated with perturbation modes that do not affect synchronization. Thus, the Lyapunov exponent λ max determining the stability is the largest real component among the remaining eigenvalues: λ max (ε) = max j n max {±} Re ( ν j± (ε) ). (S6) For ε =, which corresponds to β i = β opt, it follows from equation (S5) and α α j, j, that λ max () = β opt = α λ max opt. (S7) Our goal here is to show that this is a local minimum of λ max (ε) for any given direction β. To help with our analysis, we introduce auxiliary variables by writing the characteristic polynomial of M as (ν + c 1 ν + d 1 )(ν + c ν + d ) (ν + c n ν + d n ) (S8) and indexing the coefficients c j and d j so that the jth quadratic factor corresponds to the pair ν j±, i.e., ν + c j ν + d j =(ν ν j+ )(ν ν j ). The coefficients c j (ε) and d j (ε) are thus related to ν j± (ε) by c j = (ν j+ + ν j ) and d j = ν j+ ν j, (S9) or equivalently, ν j± = c j ± 1 c j 4d j. (S1) Equations (S5) and (S9) can be used to show that for ε =, we have c j () = β opt and d j () = α j. From equations (S6) and (S1), it follows that the Lyapunov exponent λ max (ε) is bounded from below by a function of c and d : λ max (ε) Re ν + = Re ( c + 1 ) c 4d f(c,d ). (S11) Since f(c (ε),d (ε)) = λ max opt when ε =, it suffices to show that f(c (ε),d (ε)) >λ max opt for all sufficiently small ε = in order to establish that ε =is a local minimum of λ max (ε). The lower 3 NATURE PHYSICS 3
4 DOI: 1.138/NPHYS535 bound function f(c,d ), whose functional form is system independent, is shown in Supplementary Fig. S. Note that the set of points at which f(c,d ) is not smooth form the curve c =4d, and at these points the associated eigenvalues are degenerate, i.e., ν + = ν. It is evident from the figure (and also from equation (S11)) that the region in the (c,d )-plane characterized by f(c,d ) >λ max opt can be expressed as { (c,d ): c <β opt and/or d α < } α. (S1) c β opt We thus aim to show that the point (c (ε),d (ε)) remains is this region for small ε = for any given direction β. We now expand the functions c j (ε) and d j (ε) in terms of ε around the points c j () = β opt and d j () = α j, respectively: c j = β opt + c (1) j ε + c () j ε + O(ε 3 ), d j = α j + d (1) j ε + d () j ε + O(ε 3 ). (S13) To derive equations satisfied by the linear coefficients c (1) j and d (1) j, we first write the characteristic polynomial as det(m ν1), set it equal to expression (S8), and compare the coefficients of ν k to obtain an equation relating c j and d j to the components of the matrix B for each k. We then substitute β i = β opt + εβ i and expansions (S13) into this equation and equate the first-order terms in ε to obtain a system of linear equations. In matrix form, this system reads Cc (1) + Dd (1) = Cb, (S14) where c (1) =(c (1) 1,...,c (1) n ) T, d (1) =(d (1) 1,...,d (1) n ) T, b =(b 1,...,b n ) T, and b j = i u jiv ji β i. Here C and D are constant matrices that depend only on α j (recall that α j is the jth eigenvalue of the matrix P), and u ji and v ji are the ith component of the left and right eigenvectors of P corresponding to the eigenvalue α j, respectively. We see that a solution of equation (S14) is given by c (1) j = b j and d (1) j =for all j, which corresponds to c j = β opt + b j ε + O(ε ), d j = α j + O(ε ). (S15) This means that to first order in ε, the function c j (ε) is linear (b j ) or constant (b j =), and d j (ε) is constant. In particular, the behavior of the point (c,d ) as ε varies, which determines whether it remains in region (S1), depends on the value of b. 4 4 NATURE PHYSICS
5 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION If b, this leads to a first-order approximation formula for f(c,d ): f(c,d ) λ max opt b ε/, εb, λ max opt + β opt b ε/, εb >. (S16) The function f attains the minimum value λ max opt at ε =, and the approximation is valid for all directions β = (β 1,...,β n) T outside the (n 1)-dimensional hyperplane L defined by b = i u iv i β i =in the space of all β i values. Thus, along such a direction, the Lyapunov exponent λ max (ε) attains a local minimum λ max opt at ε =. This establishes that the parameter choice β i = β opt locally minimizes λ max in an arbitrary power-grid network along almost all directions in the β -space (i.e., all directions, except those corresponding to b =). Supplementary Fig. S3 (black dashed curves) shows the behavior of c (ε) and d (ε) in equation (S15), as well as of the corresponding λ max (ε), along a randomly chosen direction in the β i -space of the 3-generator test system in Table 1 of the main text. The same conclusion holds true even if c (1) b, as long as c (1) and d (1) =. In this general case, we define the hyperplane L by c (1) = i c(1) i β i =, where c (1) i are the linear coefficients that arise in the full n-dimensional expansion of c and d : c (εβ )=β opt + d (εβ )=α + n c (1) i=1 n d (1) i=1 i εβ i + i k i εβ i + i k c () ik (εβ i)(εβ k)+o ( (ε β ) 3), d () ik (εβ i)(εβ k)+o ( (ε β ) 3). (S17) Since we can write d (1) = i d(1) i β i using the linear coefficients of d in equation (S17), the local minimum of λ max at β opt is established for all directions outside L (for which c (1) ) if d (1) i = for all i. This condition is numerically validated below. If b =(or c (1) =in the more general case), which corresponds to the hyperplane L, the first-order approximation is not sufficient to describe the local behavior of the Lyapunov exponent. We thus use the second-order terms in ε to show that λ max is locally minimum at β opt along all directions β L. The condition for (c,d ) to stay in region (S1) for small ε can be expressed in terms of the second-order coefficients as d () c () < α, (S18) 5 NATURE PHYSICS 5
6 DOI: 1.138/NPHYS535 since d α = d() ε + O(ε 3 ) c β opt c () ε + O(ε 3 ) d () as ε. c () In order to verify condition (S18) for the power-grid systems in Table 1, we rewrite it in terms of β i as d () = < α c (), (S19) i k d() ik β iβ k i k c() ik β i β k where c () ik and d() ik are the quadratic coefficients in equation (S17). We have thus reduced the problem of establishing the local optimality of λ max (β opt + εβ ) along the directions with b = to the problem of computing the n(n + 1)/ quadratic coefficients in equation (S17) and using the computed coefficients to verify condition (S19) for all β L. As a first step toward verifying condition (S19), we determine the linear and quadratic coefficients in equation (S17) for a neighborhood of size ε (ε = 1 for the Guatemala and Italy networks and ε = 1 4 for all other networks) from the values of c and d computed at n(n+3)/ sample points in the β -space (the minimum number of points required to determine the n linear coefficients and n(n + 1)/ quadratic coefficients of each function). For a given ε, we place two sample points on each β i-axis (β i = ε or ε and β k =for all k =i) and one for each pair of coordinate axes (β i = ε and β k = ε for a pair i and k =i, and β l =for all other l). The values of c and d are computed from the estimated ν + and ν through relations (S9). The resulting quadratic approximations for c (ε) and d (ε), as well as for the maximum Lyapunov exponent λ max (ε), are shown in Supplementary Fig. S3 (black solid curves) for a typical direction β outside the hyperplane L. Before proceeding to verify condition (S19) using the computed coefficients, we note that for all systems we consider, the linear coefficients satisfy d (1) i =within numerical precision; this condition was used above to study the local behavior of λ max (β opt + εβ ) within the hyperplane L. Due to the large condition numbers of the matrices involved in the computation, the observed numerical errors are larger than the machine precision. We thus further validate d (1) i = by quadratically fitting many values of d (taken to be 1 in our simulations) computed along each β i -axis to estimate d (1) i more accurately. Assuming that numerical errors in computing d are Gaussian, such estimates of d (1) i follow a distribution with a standard deviation σ. Using the observed residuals in the quadratic fit to estimate σ, we find that for each system and for each i the standard deviation σ is at most of order 1 9 and zero is within σ of our estimate of d (1) i. This provides further validation of the local behavior of λ max established above for all power grids. 6 6 NATURE PHYSICS
7 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION Given the values of the quadratic coefficients c () ik and d() ik, along with eigenvector components u i and v i, we now seek to verify condition (S19) for all β in the hyperplane L. To this end, we first formulate an auxiliary constrained optimization problem: Find β =(β 1,...,β n) T that maximizes subject to i k i k d () ik β iβ k c () ik β iβ k = K and b = i u i v i β i = (S) for a given arbitrary constant K. If condition (S19) is satisfied for all solutions of this optimization problem for all K (i.e., the maximum value found is less than K α ), then the condition is also satisfied for all directions within the hyperplane L. Using the method of Lagrange multipliers, this optimization problem can be reduced to the generalized eigenvalue problem Dη = µ Cη, where the matrices D and C are defined by D = R T ˆDR, C = R T ĈR, ˆDik = d () ik + d() ii δ ik, and Ĉik = c () ik + c() ii δ ik. Here R is an n (n 1) matrix of orthonormal basis vectors for the hyperplane L. Specifically, a solution β of the optimization problem must satisfy β = Rη k, where η k is a generalized eigenvector satisfying Dη k = µ k Cηk. Note that the eigenvectors do not depend on the constant K except via the normalization factors, since the matrices D and C do not depend on K. Since the r.h.s. of inequality (S19) does not change upon scaling of β, the condition can be verified using normalized eigenvectors η k without loss of generality. It can also be shown that the verification of condition (S19) is independent of the choice of the matrix R. Thus, for a given system, the local minimum of λ max (β opt + εβ ) can be established for all directions β by computing the generalized eigenvectors η k with an arbitrary choice of K and R, and verifying condition (S19) for each β = Rη k. To proceed with the verification of condition (S19), we implement an additional numerical step in the procedure. This is necessary because the variations of c and d for ε ε ε along the direction corresponding to a generalized eigenvector η k can be smaller than the numerical error involved in the computation if (as in the procedure above) the sample points are not taken along this direction, and this could prevent the estimates of c () and d () from being sufficiently accurate. We thus compute c () and d () in equations (S13) directly through a quadratic fit using a number points much larger than 3 sampled along this particular direction on an interval of ε much wider than ε (to capture sufficiently large variations in c and d ). Carrying this out for β corresponding to each generalized eigenvector η k, we find that condition (S19) is indeed satisfied 7 NATURE PHYSICS 7
8 DOI: 1.138/NPHYS535 for each system with a significant margin (see Supplementary Table S4 for the smallest margins we found). This shows that the rate of increase in λ max (ε) is strictly positive, thus showing that ε = is locally a minimum, for all directions β in the hyperplane L. Together with the result from the first-order approximation outside the hyperplane L, this establishes the local optimality of λ max at β i = β opt along all directions β (within and outside L). Supplementary Fig. S4 illustrates the behavior of c (ε), d (ε), and λ max (ε) for the 3-generator test system in Table 1 along the direction corresponding to the maximum of d () /c (), which is also the direction β of least increase in λ max (ε). S5. Further implications. The optimization in equation (1) of the main text can be effective even when modeling load dynamics 5 (Supplementary Fig. S5). In that context, our approach for enhancing synchronization stability is complementary to an approach proposed recently 51 to mitigate saddle-node instability by adjusting power scheduling or line impedances (Supplementary Fig. S6). Conversely, the adjustment of power scheduling and/or line impedances can also be effective within the reduced network model used in our study (Supplementary Fig. S7). Supplementary References 47. Pai, M. Energy function analysis for power system stability (Kluwer Academic Publishers, Norwell, 1989). 48. Vittal, V. Transient stability test systems for direct stability methods. IEEE T. Power Syst. 7, (199). 49. Zimmerman, R. D., Murillo-Sánchez, C. E. & Thomas, R. J. MATPOWER: Steady-state operations, planning and analysis tools for power systems research and education. IEEE T. Power Syst. 6, 1 19 (11). 5. Bergen, A. R. & Hill, D. J. A structure preserving model for power system stability analysis. IEEE T. Power Ap. Syst. PAS-1, 5 35 (1981). 51. Mallada, E. & Tang, A. Improving damping of power networks: Power scheduling and impedance adaptation. IEEE Conference on Decision and Control and European Control Conference (11). 8 NATURE PHYSICS
9 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION Supplementary Figures Supplementary Fig. S1: Relationship between the power produced by individual generators and their parameters in the three test systems of Table 1. a, The transient reactance x d,i versus the power P i injected into the network. b, The inertia constant H i versus P i. We used the approximate functional relations revealed by these data to estimate the transient reactance and the inertia constant for each generator in the Guatemala, Northern Italy, and Poland systems. NATURE PHYSICS 9
10 DOI: 1.138/NPHYS535 Slope Supplementary Fig. S: Lower bound function f(c,d ) for the Lyapunov exponent λ max. Region (S1), in which f(c,d ) >λ max opt holds true, is shown by the shades of blue, while the other region is shown by the shades of green. The whole color scheme is used to indicate the values of f(c,d ). In Supplementary Section S4, we show that, given a direction β in the β i -space, the point (c (ε),d (ε)) remains in the blue region for any sufficiently small ε =, thus establishing λ max (ε) >λ max opt. 1 1 NATURE PHYSICS
11 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION 1 a b c d Supplementary Fig. S3: Synchronization stability near β i = β opt along a typical direction β. For a fixed direction β =(β 1,β,β 3) T in the β i -space, which we chose randomly and normalized to unit length, we consider the deviation β i = β opt + εβ i with ε for the 3-generator test system in Table 1. In all panels, the blue dots represent values computed directly from the eigenvalues of the matrix M defined in Supplementary Section S4, and the black curves represent the linear (dashed) and quadratic (solid) approximation from expansion (S13) in the same section. a, b, The characteristic polynomial coefficients c (panel a) and d (panel b) as a function of ε, relative to their nominal values at ε =(which are β opt and α, respectively). c, The path followed by the point (c,d ) as ε varies in the interval [, ]. The red curve indicates the points where the Lyapunov exponents are degenerate (ν + = ν ) and f(c,d ) is not smooth as a function of c and d (Supplementary Fig. S). It is evident that no blue dots (and hence no eigenvalues of M) fall in the gray shaded region, in which f(c,d ) <λ max opt = α. d, The maximum Lyapunov exponent λ max as a function of ε. NATURE PHYSICS 11
12 DOI: 1.138/NPHYS535.1 a b c.6 d.4 Slope. 1 1 Supplementary Fig. S4: Synchronization stability near β i = β opt along the direction β of least increase in λ max. The information presented is the same as in Supplementary Fig. S3 except that β is chosen to be the direction that gives the maximum value of d () /c() (which corresponds to the approximate slope of the curves followed by the blue dots in panel c). 1 NATURE PHYSICS
13 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION Improvement factor Supplementary Fig. S5: Enhancement of synchronization stability in power-grid model incorporating load dynamics. Colors indicate the factor by which the stability improves after adjusting all β i to the value β opt, where the darkest red corresponds to any factor > 7 and black to any factor < 1. The stability of a synchronous state is measured by the largest nonzero Lyapunov exponent λ max computed within the Bergen-Hill model 5. In this model, instead of eliminating non-generator nodes by the Kron reduction, load dynamics is modeled by assuming that the real power is a linear function of the voltage frequency for all power-consuming nodes, while the voltage magnitude is assumed to be constant for all nodes. For illustration, we use the 3-generator system studied in ref. 51, with the power injection modified to P i =5, 6, 7 per unit. The improvement factor is shown as a function of the parameters D gen and D load, where D gen is a coefficient given by (D i +1/R i )/ω R (assumed be the same for all generator nodes) and D load is the frequency coefficient (assumed to be the same for all the other nodes). In most cases, our method yields significant stability improvement, demonstrating its effectiveness beyond the reduced network model. NATURE PHYSICS 13
14 DOI: 1.138/NPHYS535 Supplementary Fig. S6: Combination of complementary approaches to mitigate instabilities associated with saddle-node bifurcations. As an example, we use the same 3-generator system described in the caption of Supplementary Fig. S5 with power injection P i =7.994, 3.6, 7., simulated using the Bergen-Hill model 5. For D gen = D load =1per unit, the system is near a saddle-node bifurcation, and thus the Lyapunov exponent λ max is negative real and close to zero. At each iteration of the gradient descent-like method described in ref. 51, which changes the power injected by the generators (independently of the values of β i and thus of D i ), we compute λ max both with and without adjusting β i to β opt by tuning the damping coefficients D i. The improvement factor, measured by the ratio between the smallest λ max obtained through the iterative process with and without the β i -adjustment, is shown as a function of D gen for different values of D load. In most cases, significant additional improvement results from the adjustment of β i, illustrating that near instabilities the two methods can be combined to achieve enhancement not possible by either method alone. This result is robust against heterogeneity in the network parameters, as illustrated in the inset histogram for the system with D gen =.5and D load =1(blue dot in the main plot), where each of these coefficients is independently perturbed according to the Gaussian distribution with mean zero and standard deviation NATURE PHYSICS
15 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION Iterations Iterations Supplementary Fig. S7: Enhancement of synchronization stability via tuning of power injection within the reduced network model. We use the same 3-generator system as in Supplementary Figs. S5 and S6, where all system parameters are the same as in ref. 51, except for the power injection and consumption, which are initially set to 6% of the values used in that reference. This makes the system unstable, as indicated by the positivity of the Lyapunov exponent λ max for the synchronous state. The power injections were then adjusted iteratively according to the gradient descent-like scheme of ref. 51 and λ max is plotted against the number of iterations for several values of the coefficient D gen. In all cases, the power adjustment stabilizes the synchronous state, which is a consequence of increasing α across zero, as shown in the inset. Note that α is independent of D gen because the synchronous state in the reduced network model is unaffected by changes in D gen. 15 NATURE PHYSICS 15
16 DOI: 1.138/NPHYS535 Supplementary Tables Supplementary Table S1: Phase difference in synchronous states and the real versus imaginary parts of the admittances in the physical and effective networks. Physical network Effective network System Mean δij /π Mean G ij Mean B ij Mean G ij Mean B ij 3-generator system generator system generator system Guatemala Northern Italy: original x d,i adjusted Poland: original x d,i adjusted The real and imaginary components of the admittance Y ij are denoted G ij and B ij, respectively. Supplementary Table S: -norm of the symmetric and antisymmetric parts of the matrix P. System Symmetric Antisymmetric 3-generator system generator system generator system Guatemala Northern Italy: original.45.3 x d,i adjusted Poland: original x d,i adjusted NATURE PHYSICS
17 DOI: 1.138/NPHYS535 SUPPLEMENTARY INFORMATION Supplementary Table S3: Fraction of positive off-diagonal elements in the matrices (B ij ), (cos δij ), and (B ij cos δij ). System B ij cos δij B ij cos δij 3-generator system generator system generator system Guatemala Northern Italy: original x d,i adjusted Poland: original x d,i adjusted Supplementary Table S4: Property of the estimated quadratic coefficients of c and d. ( System min β α d () ) /c() 3-generator system generator system generator system.87 Guatemala.9 Northern Italy (x d,i adjusted) 1.85 Poland (x d,i adjusted).61 The minimum is taken over all directions β within the hyperplane L. 17 NATURE PHYSICS 17
arxiv: v3 [physics.soc-ph] 21 Dec 2018
Spontaneous synchrony in power-grid networks Adilson E. Motter 1,2,, Seth A. Myers 3, Marian Anghel 4 and Takashi Nishikawa 1 Nature Physics 9, 191 197 (2013) doi:10.1038/nphys2535 1 Department of Physics
More informationStability and Control of Synchronization in Power Grids
Stability and Control of Synchronization in Power Grids Adilson E. Motter Department of Physics and Astronomy Northwestern Institute on Complex Systems Northwestern University, Evanston, USA Power Generators
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationREAL-TIME TRANSIENT STABILITY ASSESSMENT
CHAPTER 3 REAL-TIME TRANSIENT STABILITY ASSESSMENT In Chapter 2, the basic multimachine theory and the different dynamic, state and transient energy equations have been presented and the basic criterion
More informationSCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS. Transient Stability LECTURE NOTES SPRING SEMESTER, 2008
SCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS LECTURE NOTES Transient Stability SPRING SEMESTER, 008 October 7, 008 Transient Stability Transient stability refers to the ability of a synchronous
More informationCURENT Course Power System Coherency and Model Reduction
CURENT Course Power System Coherency and Model Reduction Prof. Joe H. Chow Rensselaer Polytechnic Institute ECSE Department November 1, 2017 Slow Coherency A large power system usually consists of tightly
More informationA Generalized Hamiltonian Model for Power System Dynamics with Relay Action. Raja Timihiri. Project Advisor: Christopher DeMarco
A Generalized Hamiltonian Model for Power System Dynamics with Relay Action by Raja Timihiri Project Advisor: Christopher DeMarco In partial fulfillment of the degree of Masters In Electrical Engineering
More informationExamples of Applications of Potential Functions in Problem Solving (Web Appendix to the Paper)
Examples of Applications of otential Functions in roblem Solving (Web Appendix to the aper) Ali Mehrizi-Sani and Reza Iravani May 5, 2010 1 Introduction otential functions may be exploited to formulate
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationPower Grid State Estimation after a Cyber-Physical Attack under the AC Power Flow Model
Power Grid State Estimation after a Cyber-Physical Attack under the AC Power Flow Model Saleh Soltan, Gil Zussman Department of Electrical Engineering Columbia University, New York, NY Email: {saleh,gil}@ee.columbia.edu
More informationPOWER SYSTEM DYNAMIC SECURITY ASSESSMENT CLASSICAL TO MODERN APPROACH
Abstract POWER SYSTEM DYNAMIC SECURITY ASSESSMENT CLASSICAL TO MODERN APPROACH A.H.M.A.Rahim S.K.Chakravarthy Department of Electrical Engineering K.F. University of Petroleum and Minerals Dhahran. Dynamic
More informationPower Grid Partitioning: Static and Dynamic Approaches
Power Grid Partitioning: Static and Dynamic Approaches Miao Zhang, Zhixin Miao, Lingling Fan Department of Electrical Engineering University of South Florida Tampa FL 3320 miaozhang@mail.usf.edu zmiao,
More informationSynchronization in power-grid networks
Technische Universiteit Delft Faculteit Elektrotechniek, Wiskunde en Informatica Delft Institute of Applied Mathematics Synchronization in power-grid networks Report on behalf of the Delft Institute of
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationCHAPTER 2 DYNAMIC STABILITY MODEL OF THE POWER SYSTEM
20 CHAPTER 2 DYNAMIC STABILITY MODEL OF THE POWER SYSTEM 2. GENERAL Dynamic stability of a power system is concerned with the dynamic behavior of the system under small perturbations around an operating
More informationThe N k Problem using AC Power Flows
The N k Problem using AC Power Flows Sean Harnett 5-19-2011 Outline Introduction AC power flow model The optimization problem Some results Goal: find a small set of lines whose removal will cause the power
More informationSupplementary Information for Enhancing synchronization stability in a multi-area power grid
Supplementary Information for Enhancing synchronization stability in a multi-area power grid Bing Wang, 1, 2, Hideyuki Suzuki, 3 and Kazuyuki Aihara 2 1 School of Computer Engineering and Science, Shanghai
More informationKernel Methods. Machine Learning A W VO
Kernel Methods Machine Learning A 708.063 07W VO Outline 1. Dual representation 2. The kernel concept 3. Properties of kernels 4. Examples of kernel machines Kernel PCA Support vector regression (Relevance
More informationSimulating a Power System
Simulating a Power System Presented by Prof. Tyrone Fernando School of Electrical and Electronic Engineering (EECE), University of Western Australia (UWA) 1. Motivations In an actual power system, it is
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More information254 IEEE TRANSACTIONS ON CONTROL OF NETWORK SYSTEMS, VOL. 2, NO. 3, SEPTEMBER 2015
254 IEEE TRANSACTIONS ON CONTROL OF NETWORK SYSTEMS, VOL 2, NO 3, SEPTEMER 215 The Price of Synchrony: Evaluating the Resistive Losses in Synchronizing Power Networks Emma Tegling, assam amieh, Fellow,
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationGradient Descent. Dr. Xiaowei Huang
Gradient Descent Dr. Xiaowei Huang https://cgi.csc.liv.ac.uk/~xiaowei/ Up to now, Three machine learning algorithms: decision tree learning k-nn linear regression only optimization objectives are discussed,
More informationMinimax Approximation Synthesis in PSS Design by Embedding
Minimax Approximation Synthesis in PSS Design by Embedding Gravitational Search Algorithm Dr. Akash Saxena Department of Electrical Engineering Swami Keshvanand Institute of Technology Jaipur, India Power
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationAlgorithmic Construction of Lyapunov Functions for Power System Stability Analysis
1 Algorithmic Construction of Lyapunov Functions for Power System Stability Analysis M. Anghel, F. Milano, Senior Member, IEEE, and A. Papachristodoulou, Member, IEEE Abstract We present a methodology
More informationConvex Optimization. Newton s method. ENSAE: Optimisation 1/44
Convex Optimization Newton s method ENSAE: Optimisation 1/44 Unconstrained minimization minimize f(x) f convex, twice continuously differentiable (hence dom f open) we assume optimal value p = inf x f(x)
More informationStatistical Machine Learning from Data
Samy Bengio Statistical Machine Learning from Data 1 Statistical Machine Learning from Data Support Vector Machines Samy Bengio IDIAP Research Institute, Martigny, Switzerland, and Ecole Polytechnique
More informationEVALUATION OF THE IMPACT OF POWER SECTOR REFORM ON THE NIGERIA POWER SYSTEM TRANSIENT STABILITY
EVALUATION OF THE IMPACT OF POWER SECTOR REFORM ON THE NIGERIA POWER SYSTEM TRANSIENT STABILITY F. I. Izuegbunam * Department of Electrical & Electronic Engineering, Federal University of Technology, Imo
More information1. Introduction. Keywords Transient Stability Analysis, Power System, Swing Equation, Three-Phase Fault, Fault Clearing Time
Energy and Power 17, 7(1): -36 DOI: 1.593/j.ep.1771.3 Numerical Simulations for Transient Stability Analysis of Two-Machine Power System Considering Three-Phase Fault under Different Fault Clearing Times
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationKinematics of fluid motion
Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx
More informationFinal Exam, Second Semester: 2015/2016 Electrical Engineering Department
Philadelphia University Faculty of Engineering Student Name Student No: Serial No Final Exam, Second Semester: 2015/2016 Electrical Engineering Department Course Title: Power II Date: 21 st June 2016 Course
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationControl of Chatter using Active Magnetic Bearings
Control of Chatter using Active Magnetic Bearings Carl R. Knospe University of Virginia Opportunity Chatter is a machining process instability that inhibits higher metal removal rates (MRR) and accelerates
More informationYou know for EE 303 that electrical speed for a generator equals the mechanical speed times the number of poles, per eq. (1).
Stability 1 1. Introduction We now begin Chapter 14.1 in your text. Our previous work in this course has focused on analysis of currents during faulted conditions in order to design protective systems
More informationAppendix A Taylor Approximations and Definite Matrices
Appendix A Taylor Approximations and Definite Matrices Taylor approximations provide an easy way to approximate a function as a polynomial, using the derivatives of the function. We know, from elementary
More informationLet x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call. r := b A x.
ESTIMATION OF ERROR Let x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call the residual for x. Then r := b A x r = b A x = Ax A x = A
More informationLeast Sparsity of p-norm based Optimization Problems with p > 1
Least Sparsity of p-norm based Optimization Problems with p > Jinglai Shen and Seyedahmad Mousavi Original version: July, 07; Revision: February, 08 Abstract Motivated by l p -optimization arising from
More informationAnalysis of Coupling Dynamics for Power Systems with Iterative Discrete Decision Making Architectures
Analysis of Coupling Dynamics for Power Systems with Iterative Discrete Decision Making Architectures Zhixin Miao Department of Electrical Engineering, University of South Florida, Tampa FL USA 3362. Email:
More information1 Unified Power Flow Controller (UPFC)
Power flow control with UPFC Rusejla Sadikovic Internal report 1 Unified Power Flow Controller (UPFC) The UPFC can provide simultaneous control of all basic power system parameters ( transmission voltage,
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationCoordinate Descent and Ascent Methods
Coordinate Descent and Ascent Methods Julie Nutini Machine Learning Reading Group November 3 rd, 2015 1 / 22 Projected-Gradient Methods Motivation Rewrite non-smooth problem as smooth constrained problem:
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More informationModel of Induction Machine to Transient Stability Programs
Model of Induction Machine to Transient Stability Programs Pascal Garcia Esteves Instituto Superior Técnico Lisbon, Portugal Abstract- this paper reports the work performed on the MSc dissertation Model
More informationRobust Tuning of Power System Stabilizers Using Coefficient Diagram Method
International Journal of Electrical Engineering. ISSN 0974-2158 Volume 7, Number 2 (2014), pp. 257-270 International Research Publication House http://www.irphouse.com Robust Tuning of Power System Stabilizers
More informationLecture Notes: Geometric Considerations in Unconstrained Optimization
Lecture Notes: Geometric Considerations in Unconstrained Optimization James T. Allison February 15, 2006 The primary objectives of this lecture on unconstrained optimization are to: Establish connections
More informationHopf bifurcations induced by SVC Controllers: A didactic example
Electric Power Systems Research 77 (2007) 234 240 Hopf bifurcations induced by SVC Controllers: A didactic example Wei Gu a,, Federico Milano b, Ping Jiang a, Guoqing Tang a a SouthEast University, Department
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 2014-15) Time Allowed: 2 Hours INSTRUCTIONS TO CANDIDATES 1. Please write your student number only. 2. This examination
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationPerformance of Flocking-Based Control Schemes in Smart Grid Applications
Performance of Flocking-Based Control Schemes in Smart Grid Applications Abdallah K. Farraj Eman M. Hammad Jin Wei Deepa Kundur Karen L. Butler-Purry Department of Electrical and Computer Engineering,
More informationScientific Computing
Scientific Computing Direct solution methods Martin van Gijzen Delft University of Technology October 3, 2018 1 Program October 3 Matrix norms LU decomposition Basic algorithm Cost Stability Pivoting Pivoting
More informationB553 Lecture 5: Matrix Algebra Review
B553 Lecture 5: Matrix Algebra Review Kris Hauser January 19, 2012 We have seen in prior lectures how vectors represent points in R n and gradients of functions. Matrices represent linear transformations
More informationarxiv: v1 [cs.sy] 24 Mar 2016
Nonlinear Analysis of an Improved Swing Equation Pooya Monshizadeh, Claudio De Persis, Nima Monshizadeh, and Arjan van der Schaft arxiv:603.07440v [cs.sy] 4 Mar 06 Abstract In this paper, we investigate
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationJACOBI S ITERATION METHOD
ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes
More informationPower system modelling under the phasor approximation
ELEC0047 - Power system dynamics, control and stability Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct October 2018 1 / 16 Electromagnetic transient vs. phasor-mode simulations
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationHere each term has degree 2 (the sum of exponents is 2 for all summands). A quadratic form of three variables looks as
Reading [SB], Ch. 16.1-16.3, p. 375-393 1 Quadratic Forms A quadratic function f : R R has the form f(x) = a x. Generalization of this notion to two variables is the quadratic form Q(x 1, x ) = a 11 x
More informationStable Limit Cycle Generation for Underactuated Mechanical Systems, Application: Inertia Wheel Inverted Pendulum
Stable Limit Cycle Generation for Underactuated Mechanical Systems, Application: Inertia Wheel Inverted Pendulum Sébastien Andary Ahmed Chemori Sébastien Krut LIRMM, Univ. Montpellier - CNRS, 6, rue Ada
More informationEconomic Operation of Power Systems
Economic Operation of Power Systems Section I: Economic Operation Of Power System Economic Distribution of Loads between the Units of a Plant Generating Limits Economic Sharing of Loads between Different
More informationSolving Linear Systems of Equations
November 6, 2013 Introduction The type of problems that we have to solve are: Solve the system: A x = B, where a 11 a 1N a 12 a 2N A =.. a 1N a NN x = x 1 x 2. x N B = b 1 b 2. b N To find A 1 (inverse
More informationReflections and Rotations in R 3
Reflections and Rotations in R 3 P. J. Ryan May 29, 21 Rotations as Compositions of Reflections Recall that the reflection in the hyperplane H through the origin in R n is given by f(x) = x 2 ξ, x ξ (1)
More informationLinear Algebra and Matrices
Linear Algebra and Matrices 4 Overview In this chapter we studying true matrix operations, not element operations as was done in earlier chapters. Working with MAT- LAB functions should now be fairly routine.
More informationCHAPTER 11. A Revision. 1. The Computers and Numbers therein
CHAPTER A Revision. The Computers and Numbers therein Traditional computer science begins with a finite alphabet. By stringing elements of the alphabet one after another, one obtains strings. A set of
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationECEN 667 Power System Stability Lecture 25: FFT, Energy Methods
ECEN 667 Power System Stability Lecture 5: FFT, Energy Methods Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements Read Chapter 9 Final
More informationLinear Algebra for Machine Learning. Sargur N. Srihari
Linear Algebra for Machine Learning Sargur N. srihari@cedar.buffalo.edu 1 Overview Linear Algebra is based on continuous math rather than discrete math Computer scientists have little experience with it
More informationSymmetric Matrices and Eigendecomposition
Symmetric Matrices and Eigendecomposition Robert M. Freund January, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 2 1 Symmetric Matrices and Convexity of Quadratic Functions
More informationPower Angle & Time Domain Plots 1.0 The power-angle curves
Power Angle & Time Domain Plots 1.0 The power-angle curves Consider a generator connected to an infinite bus through a network, as shown in Fig. 1 below (all in pu). E E NETWORK X V V 0 Fig. 1 We assume
More informationStability Analysis Through the Direct Method of Lyapunov in the Oscillation of a Synchronous Machine
Modern Applied Science; Vol. 12, No. 7; 2018 ISSN 1913-1844 E-ISSN 1913-1852 Published by Canadian Center of Science and Education Stability Analysis Through the Direct Method of Lyapunov in the Oscillation
More informationLeast Squares Optimization
Least Squares Optimization The following is a brief review of least squares optimization and constrained optimization techniques, which are widely used to analyze and visualize data. Least squares (LS)
More informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationPOWER SYSTEM STABILITY
LESSON SUMMARY-1:- POWER SYSTEM STABILITY 1. Introduction 2. Classification of Power System Stability 3. Dynamic Equation of Synchronous Machine Power system stability involves the study of the dynamics
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5
Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5 Instructor: Farid Alizadeh Scribe: Anton Riabov 10/08/2001 1 Overview We continue studying the maximum eigenvalue SDP, and generalize
More informationECEN 667 Power System Stability Lecture 15: PIDs, Governors, Transient Stability Solutions
ECEN 667 Power System Stability Lecture 15: PIDs, Governors, Transient Stability Solutions Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements
More informationarxiv: v1 [nlin.ao] 8 Dec 2017
Nonlinearity 31(1):R1-R23 (2018) Identical synchronization of nonidentical oscillators: when only birds of different feathers flock together arxiv:1712.03245v1 [nlin.ao] 8 Dec 2017 Yuanzhao Zhang 1 and
More informationPhase Boundary Computation for Fault Induced Delayed Voltage Recovery
IEEE th Annual Conference on Decision and Control (CDC) December -,. Osaka, Japan Phase Boundary Computation for Fault Induced Delayed Voltage Recovery Michael W. Fisher Ian A. Hiskens Abstract Distribution
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationCHAPTER 3 MATHEMATICAL MODELING OF HYDEL AND STEAM POWER SYSTEMS CONSIDERING GT DYNAMICS
28 CHAPTER 3 MATHEMATICAL MODELING OF HYDEL AND STEAM POWER SYSTEMS CONSIDERING GT DYNAMICS 3.1 INTRODUCTION This chapter focuses on the mathematical state space modeling of all configurations involved
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
More informationDistributed Stability Analysis and Control of Dynamic Networks
Distributed Stability Analysis and Control of Dynamic Networks M. Anghel and S. Kundu Los Alamos National Laboratory, USA August 4, 2015 Anghel & Kundu (LANL) Distributed Analysis & Control of Networks
More informationSynchronization Stability of Coupled Near-Identical Oscillator Network
Synchronization Stability of Coupled Near-Identical Oscillator Network Jie Sun, Erik M. Bollt, and Takashi Nishikawa Department of Mathematics and Computer Science Clarkson University Potsdam, NY 13699-5815,
More informationPredicting, controlling and damping inter-area mode oscillations in Power Systems including Wind Parks
3rd IASME/WSEAS Int. Conf. on Energy & Environment, University of Cambridge, UK, February 3-5, 008 Predicting, controlling and damping inter-area mode oscillations in Power Systems including Wind Parks
More informationEE2351 POWER SYSTEM ANALYSIS
EE351 POWER SYSTEM ANALYSIS A.Ahamed Riazudeen EEE DEPARTMENT 1 UNIT I INTRODUCTION Power system network 3 SINGLE LINE DIAGRAM It is a diagrammatic representation of a power system in which the components
More informationCS 246 Review of Linear Algebra 01/17/19
1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector
More informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More informationFOR REDUCE SUB-SYNCHRONOUS RESONANCE TORQUE BY USING TCSC
FOR REDUCE SUB-SYNCHRONOUS RESONANCE TORQUE BY USING TCSC Shrikant patel 1, N.K.Singh 2, Tushar kumar 3 Department of Electrical Engineering, Scope College of Engineering Bhopal,(M.P.) Emil:Shrikantcseb@gmail.com
More informationDistributed Frequency Control in Power Grids Under Limited Communication
Distributed Frequency Control in Power Grids Under Limited Communication Marzieh Parandehgheibi, Konstantin Turitsyn and Eytan Modiano Laboratory for Information and Decision Systems, Massachusetts Institute
More information