ECE 422 Power System Operations & Planning 7 Transient Stability
|
|
- Elisabeth Eaton
- 6 years ago
- Views:
Transcription
1 ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun
2 References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3
3 Transient Stability The ability of the power system to maintain synchronism when subjecte to a severe transient isturbance such as a fault on transmission facilities, loss of generation, or loss of a large loa. The system response to such isturbances involves large excursions of generator rotor angles, power flows, bus voltages, an other system variables. Stability is influence by the nonlinear characteristics of the system If the resulting angular separation between the machines in the system remains within certain bouns, the system maintains synchronism. If loss of synchronism ue to transient instability occurs, it will usually be evient within 3 secons of the initial isturbances 3
4 Single Machine Infinite Bus (SMIB) System P e = =P max sin P max = E E B /X T Swing Equation: H w t = P - P = P -P m e m max sin P a is calle Accelerating Power P m = P e The rotor will accelerate if P m increases, or P e ecreases 4
5 Power Angle Relationship P e = =P max sin Question: What if both circuits are out of service? 5
6 Response to a step change in P m H = P a = Pm- Pe = Pm-P w t max sin Consier a suen increase in P m : P m P m. New equilibrium point b satisfying P e ( )=P m a b: Due to the rotor s inertia, cannot jump from to, so P a =P m -P e ( )> an r increases from. When b is reache, P a = but r >, so continues to increase. b c: > an P a <, so r ecreases until c. At c, r = an reaches the peak value max. P a r a > =, b = r, max, c < =, = max, c: At c, the rotor starts to ecelerate (since P a <) with r < an ecreases. With all resistances (amping) neglecte, an r oscillate aroun the new equilibrium point b with a constant amplitue. 6
7 Equal Area Criterion (EAC) If max exists where /t=: w D wr = ( Pm - Pe) H ò H w wr ( Pm Pe) w D = ò - æ w ö D = ( m- e) ç w ò çè ø max r J P P = max At max, r = an the integral= Moment of inertia in p.u. = ò +ò max ( - ) ( P - P ) P P m e m e (Note: all losses are neglecte) = area A - area A = 7
8 Equal Area Criterion (EAC): The stability is maintaine only if a ecelerating area A the accelerating area A can be locate above P m (from b to, i.e. the Unstable Equilibrium Point or UEP). SEP UEP If A < A, will continue increasing at UEP (since r > ), so it will lose stability. For the case with a step change in P m, the new P m oes matter for transient stability. 8
9 Transient stability limit for a step change of P m Following a step change P m P m, solve the transient stability limit of P m : Assume A = A in orer to solve the limit of P m P P P P m ò ò max ( - )- sin = sin - ( - ) max max m max ( - ) P = P (cos -cos ) max m max max - P + P (cos - cos ) =-P (cos -cos )-P m max max max m max At the UEP, Pm = Pmax sinmax ( - )sin + cos = cos max max max SEP UEP Solve max to calculate the transient stability limit for a step change of P m : P = P sin m max max = p-max Questions: Can we increase P m beyon that transient limit? If yes, how much further can we increase P m? 9
10 Solve max by the Newton Raphson metho ( - )sin + cos = cos max max max The nonlinear function form: Select an initial estimate: f ( ) = cos = c max p/< < p ( k ) max Calculate iterative solutions by the N R algorithm: + = +D ( k ) ( k) ( k) max max max where c- f( ) c- f( ) D = = ( k) ( k) ( k ) max max max ( k) ( k) f ( max - )cosmax max ( k ) max Give a solution when a specific accuracy is reache, i.e. + - e ( k ) ( k) max max
11 Response to a three phase fault P e = =P max sin P e,uring fault << P e, post-fault P e,post-fault <P e,pre-fault for a permanent fault (cleare by tripping the fault circuit) or P e, post-fault =P e,pre-fault for a temporary fault
12 Stable Unstable
13 Critical Clearing Angle (CCA) Saaat s Sec..6 (Example.5) Consier a simple case A three phase fault at the sening en P e, uring fault = if all resistances are neglecte Critical Clearing Angle c ò ò c max P = ( P sin -P) m max m A A Integrating both sies: ( ) c P - = P (cos -cos )-P ( - ) m c max c max m max c Pm cos c = ( max - c) + cos P max max Pm = ( p- - )-cos P max c = - 3
14 Critical Clearing Time (CCT) Solve the CCT from the CCA: Since P e, uring fault = for this case, uring the fault: H w t = P m t = w Pm t Pmt t H ò = w H w 4H = Pt m + t c = 4 H ( c -) w P m 4
15 For a more general case: P e (uring fault)> P 3max (pre-fault) A 3 A (post-fault) P 3max P max A (fault-on) P max s ( u ) ( ) c max P - P sn i = P sin - P ( - ) m c -ò ò max 3max m max c c cos c = P ( - ) + P cos -P cos P - P m max c 3max max max 3max max A = V ke ( c ), the kinetic energy at c A + A 3 = V( c )=V ke ( c )+V pe ( c ), total energy at c A + A 3 = V pe ( u )=V cr, i.e. the largest potential energy If an only if V( c ) V cr (i.e. A A ), the generator is stable 5
16 Factors influencing transient stability How heavily the generator is loae. The generator output uring the fault. This epens on the fault location an type The fault clearing time The post fault transmission system reactance The generator reactance. A lower reactance increases peak power an reuces initial rotor angle. The generator inertia. The higher the inertia, the slower the rate of change in angle. This reuces the kinetic energy gaine uring fault; i.e. the accelerating area A is reuce. The generator internal voltage magnitue (E ). This epens on the fiel excitation The infinite bus voltage magnitue E B 6
17 EAC for a Two Machine System Two interconnecte machines respectively with H an H The system can be reuce to an equivalent SMIB system w w t H H = ( P m- Pe ) = Pa w w t H H = ( P m- Pe) = Pa w P P ( ) t t t H H a a = - = - w HH H P - HP H P - HP a a m m e e = = - H + H t H+ H H+ H H + H H P - HP H w t = P -P m, e, P e ò,max, w H ( P - P ) = m, e, P m,,max 7
18 Methos for Transient Stability Analysis Analyzing a system s transient stability following a given contingency Time omain simulation: At present, the most practical available metho of transient stability analysis is time omain simulation in which the nonlinear ifferential equations are solve by using step by step numerical integration techniques. Direct methos: Those methos etermine stability without explicitly solving the system ifferential equations. The methos are base on Lyapunov s secon metho, efine a Transient Energy Function (TEF) as a possible Lyapunov function, an compare the TEF to a critical energy, enote by V cr, to juge stability EAC is a irect metho for a SMIB or two machine system 8
19 Numerical Integration Methos The ifferential equations to be solve in power system stability analysis are nonlinear ODEs (orinary ifferential equations) with known initial values x=x an t=t x t f ( x, t) where x is the state vector of n epenent variables an t is the inepenent variable (time). Our objective is to solve x as a function of t Explicit Methos In these methos, the value of x at any value of t is compute from the knowlege of the values of x from only the previous time steps, e.g. Euler metho an R-K methos Implicit Methos These methos use interpolation functions (involving future time steps) for the expression uner the integral, e.g. the Trapezoial Rule 9
20 Euler Metho The Euler metho is equivalent to using the first two terms of the Taylor series expansion for x aroun the point (x, t ), referre to as a first-orer metho (error is on the orer of t ) Approximate the curve at x=x an t=t by its tangent x(t ) x x t f ( xt, ) x t x f ( x, t ) x t x x t t x x x x x x t t x t t At step i+ x x x t i i t x The stanar Euler metho results in inaccuracies because it uses the erivative only at the beginning of the interval as though it applie throughout the interval i
21 Moifie Euler (ME) Metho Moifie Euler metho consists of two steps: (a) Preictor step: p x x t x t x Slope at the beginning of t x (t ) x c x p x x t f ( xt, ) t The erivative at the en of the t is estimate using x p x t x p p f( x, t ) Estimate slope at the en of t t t (b) Corrector step: x x x x t p x t t p x x t c c i x i x x t x x i i t It is a secon-orer metho (error is on the orer of t 3 ) Step size t must be small enough to obtain a reasonably accurate solution, but at the same time, large enough to avoi the numerical instability with the computer, e.g. increasing roun-off errors.
22 Runge Kutta (R K) Methos General formula of the n orer R-K metho: (error is on the orer of t 3 ) k f( x, t ) t k f( x k, t t) t x x ak ak At Step i+: x(t ) x + k x x f ( xt, ) t t t t k f ( x, t ) t i i x i x ak a k i k f( x k, t t) t i i The ME metho is a special case with a =a =/, = = General formula of the 4 th orer R-K metho: (error is on the orer of t 5 ) x i x i ( k k k3 k4 ) 6 k f( xi, ti) t k t k f( xi, ti ) t k t k3 f( xi, ti ) t k f( x k, t t) t 4 i 3 i
23 Numerical Stability of Explicit Integration Methos Numerical stability is relate to the stiffness of the set of ifferential equations representing the system The stiffness is measure by the ratio of the largest to smallest time constant, or more precisely by max / min of the linearize system. Stiffness in a transient stability simulation increases with moeling more etails (more small time constants are concerne). Explicit integration methos have weak stability numerically; with stiff systems, the solution blows up unless a small step size is use. Even after the fast moes ie out, small time steps continue to be require to maintain numerical stability 3
24 Implicit Methos Implicit methos use interpolation functions for the expression uner the integral. Interpolation implies the function must pass through the yet unknown points at t The simplest implicit integration metho is the Trapezoial Rule metho. It uses linear interpolation. The stiffness of the system being analyze affects accuracy but not numerical stability. With larger time steps, high frequency moes an fast transients are filtere out, an the solutions for the slower moes is still accurate. For example, for the Trapezoial rule, only ynamic moes faster than f(x n,t n ) an f(x n+,t n+ ) are neglecte. f x(t )=x(t )+ A + B B f(x,t ) f(x,t ) Δt x=x+ f x,t Δt x = x + f x,t + +f x,t f x, t n+ n n n n+ n+ A t Compare to ME metho: t x x f x,t + f x,t n p n n n n n+ t t 4
25 Overall System Equations The overall system equations are expresse in the general form comprising a set of first-orer ifferential equations (ynamic evices) an a set of algebraic equations (evices an network) where x V I Y N state vector of the system bus voltage vector x=f(x,v) I(x,V) = YN V DE AE current injection vector noe amittance matrix. It is constant except for changes introuce by network-switching operations; symmetrical except for issymmetry introuce by phase-shifting transformers 5
26 Solution of the Equations Schemes for the solution of equations DE an AE are characterize by the following factors The manner of interface between the DE an AE. Either a partitione approach or a simultaneous approach may be use The integration metho, i.e. an implicit metho or explicit metho, use to solve the DE. The metho use to solve the AE (power flow analysis), e.g. the Newton-Raphson metho. Most commercialize power system simulation programs provie the Moifie Euler, n orer R-K, 4 th orer R-K an Trapezoial Rule methos 6
27 A Simplifie Moel for Multi Machine Systems Consier these classic simplifying assumptions: Each synchronous machine is represente by a voltage source E with constant magnitue E behin X (neglecting armature resistances, the effect of saliency an the changes in flux linkages) The mechanical rotor angle of each machine coincies with the angle of E The governor s actions are neglecte an the input powers P mi are assume to remain constant uring the entire perio of simulation Using the pre-fault bus voltages, all loas are converte to equivalent amittances to groun. Those amittances are assume to remain constant (constant impeance loa moels) Damping or asynchronous powers are ignore. Machines belonging to the same station swing together an are sai to be coherent. A group of coherent machines is represente by one equivalent machine 7
28 Solve the initial power flow an etermine the initial bus voltage phasors V i. Terminal currents I i of m generators prior to isturbance are calculate by their terminal voltages V i an power outputs S i, an then use to calculate E i All loas are converte to equivalent amittances: To inclue voltages behin X i, a m internal generator buses to the n-bus power system network to form a n+m bus network (groun as the reference for voltages): E X E X X m E m Y bus nxn Y reuce bus m m 8
29 Noe voltage equation with groun as reference I bus is the vector of the injecte bus currents V bus is the vector of bus voltages measure from the reference noe Y bus is the bus amittance matrix : Y ii (iagonal element) is the sum of amittances connecte to bus i Y ij (off-iagonal element) equals the negative of the amittance between buses i an j Compare to the Y bus for power flow analysis, aitional m internal generator noes are ae an Y ii (i n) is moifie to inclue the loa amittance at noe i 9
30 To simplify the analysis, all noes other than the generator internal noes are eliminate as follows I where I I m where ij is the angle of Y ij nees to be upate whenever the network is change. 3
ECE 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 informationDynamics of the synchronous machine
ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course
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 informationA Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique
J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationState Space Analysis of Power System Stability Enhancement with Used the STATCOM
tate pace Analysis of Power ystem tability Enhancement with Use the ACOM M. Mahavian () - G. hahgholian () () Department of Electrical Engineering, Islamic Aza University, Naein Branch, Esfahan, Iran ()
More informationLecture 6: Control of Three-Phase Inverters
Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,
More informationDesign A Robust Power System Stabilizer on SMIB Using Lyapunov Theory
Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)
More informationPower System Stability and Control. Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India
Power System Stability and Control Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India Contents Chapter 1 Introduction to Power System Stability
More informationDetermine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Validation
Determine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Valiation Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper extens
More informationInvestigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model
ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationOutcome of this lecture
Outcome of this lecture At the en of this lecture you will be able to: List the ifferent parts of a synchronous machine Explain the operation principles of the machine Use the equivalent circuit moel of
More informationSimulink model for examining dynamic interactions involving electro-mechanical oscillations in distribution systems
University of Wollongong Research Online Faculty of Engineering an Information Sciences - Papers: Part A Faculty of Engineering an Information Sciences 205 Simulink moel for examining ynamic interactions
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 informationModelling of Three Phase Short Circuit and Measuring Parameters of a Turbo Generator for Improved Performance
Moelling of Three Phase Short Circuit an Measuring Parameters of a Turbo Generator for Improve Performance M. Olubiwe, S. O. E. Ogbogu, D. O. Dike, L. Uzoechi Dept of Electrical an Electronic Engineering,
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationECE 692 Advanced Topics on Power System Stability 2 Power System Modeling
ECE 692 Avance Topics on Power System Stability 2 Power System Moeling Spring 2016 Instructor: Kai Sun 1 Outline Moeling of synchronous generators for Stability Stuies Moeling of loas Moeling of frequency
More informationDeriving ARX Models for Synchronous Generators
Deriving AR Moels for Synchronous Generators Yangkun u, Stuent Member, IEEE, Zhixin Miao, Senior Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract Parameter ientification of a synchronous generator
More informationB.E. / B.Tech. Degree Examination, April / May 2010 Sixth Semester. Electrical and Electronics Engineering. EE 1352 Power System Analysis
B.E. / B.Tech. Degree Examination, April / May 2010 Sixth Semester Electrical and Electronics Engineering EE 1352 Power System Analysis (Regulation 2008) Time: Three hours Answer all questions Part A (10
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 informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationThe Energy Flow Approach for Oscillation Source Location and Damping Evaluation
1 The Energy Flow Approach for Oscillation Source Location an Damping Evaluation Lei CHEN PhD, Associate Professor Department of Electrical Engineering Tsinghua University chenlei08@tsinghua.eu.cn Backgroun
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationIPMSM Inductances Calculation Using FEA
X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationDynamic Modeling and Analysis of Large-scale Power Systems in the DQ0 Reference Frame
Dynamic Moeling an Analysis of Large-scale Power Systems in the DQ0 Reference Frame Juri Belikov Tallinn University of Technology juri.belikov@ttu.ee December 12, 2017 Juri Belikov (TUT) Moeling an Ientification
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationMATH2231-Differentiation (2)
-Differentiation () The Beginnings of Calculus The prime occasion from which arose my iscovery of the metho of the Characteristic Triangle, an other things of the same sort, happene at a time when I ha
More informationCritical clearing time evaluation of Nigerian 330kV transmission system
American Journal of Electrical Power and Energy Systems 2013; 2(6): 123-128 Published online October 20, 2013 (http://www.sciencepublishinggroup.com/j/epes) doi: 10.11648/j.epes.20130206.11 Critical clearing
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationModeling and analysis of parallel connected permanent magnet synchronous generators in a small hydropower plant
Proceeings of the 2006 IASME/WSEAS International Conference on Energy & Environmental Systems, Chalkia, Greece, May 8-10, 2006 (pp83-88) Moeling an analysis of parallel connecte permanent magnet synchronous
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff
More informationChapter 31: RLC Circuits. PHY2049: Chapter 31 1
Chapter 31: RLC Circuits PHY049: Chapter 31 1 LC Oscillations Conservation of energy Topics Dampe oscillations in RLC circuits Energy loss AC current RMS quantities Force oscillations Resistance, reactance,
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationChapter 6. Electromagnetic Oscillations and Alternating Current
hapter 6 Electromagnetic Oscillations an Alternating urrent hapter 6: Electromagnetic Oscillations an Alternating urrent (hapter 31, 3 in textbook) 6.1. Oscillations 6.. The Electrical Mechanical Analogy
More informationPredictive control of synchronous generator: a multiciterial optimization approach
Preictive control of synchronous generator: a multiciterial optimization approach Marián Mrosko, Eva Miklovičová, Ján Murgaš Abstract The paper eals with the preictive control esign for nonlinear systems.
More informationImpact of DFIG based Wind Energy Conversion System on Fault Studies and Power Swings
Impact of DFIG base Win Energy Conversion System on Fault Stuies an Power Swings Likin Simon Electrical Engineering Department Inian Institute of Technology, Maras Email: ee133@ee.iitm.ac.in K Shanti Swarup
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
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 informationEnhancement of transient stability analysis of multimachine power system
WWJMRD 2016; 2(6): 41-45 www.wwjmrd.com Impact Factor MJIF: 4.25 e-issn: 2454-6615 Oyediran Oyebode Olumide Department of Computer Engineering, Osun State Polytechnic Iree, Osun State, Nigeria Ogunwuyi
More informationPMU-Based System Identification for a Modified Classic Generator Model
PMU-Base System Ientification for a Moifie Classic Generator Moel Yasser Wehbe, Lingling Fan, Senior Member, IEEE Abstract The paper proposes to use PMU measurements (voltage phasor, real an reactive powers)
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationDYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES
Annals of the University of Craiova, Electrical Engineering series, No 33, 9; ISSN 184-485 7 TH INTERNATIONAL CONFERENCE ON ELECTROMECHANICAL AN POWER SYSTEMS October 8-9, 9 - Iaşi, Romania YNAMIC PERFORMANCE
More informationTwo Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices
Journal of Physics: Conference Series PAPER OPEN ACCESS Two Dimensional Numerical Simulator for Moeling NDC Region in SNDC Devices To cite this article: Dheeraj Kumar Sinha et al 2016 J. Phys.: Conf. Ser.
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationSummary: Differentiation
Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin
More informationModule 2. DC Circuit. Version 2 EE IIT, Kharagpur
Moule 2 DC Circuit Lesson 9 Analysis of c resistive network in presence of one non-linear element Objectives To unerstan the volt (V ) ampere ( A ) characteristics of linear an nonlinear elements. Concept
More informationmodel considered before, but the prey obey logistic growth in the absence of predators. In
5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an
More informationDifferentiability, Computing Derivatives, Trig Review. Goals:
Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an
More informationDefine each term or concept.
Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationMATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208
MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:
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 informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationKINGS COLLEGE OF ENGINEERING Punalkulam
KINGS COLLEGE OF ENGINEERING Punalkulam 613 303 DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING POWER SYSTEM ANALYSIS QUESTION BANK UNIT I THE POWER SYSTEM AN OVERVIEW AND MODELLING PART A (TWO MARK
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationSliding mode approach to congestion control in connection-oriented communication networks
JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University
More informationChapter 2 Derivatives
Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationTransient Stability Analysis of Single Machine Infinite Bus System by Numerical Methods
International Journal of Electrical and Electronics Research ISSN 348-6988 (online) Vol., Issue 3, pp: (58-66), Month: July - September 04, Available at: www.researchpublish.com Transient Stability Analysis
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More information6 Wave equation in spherical polar coordinates
6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationECE 585 Power System Stability
Homework 1, Due on January 29 ECE 585 Power System Stability Consider the power system below. The network frequency is 60 Hz. At the pre-fault steady state (a) the power generated by the machine is 400
More informationx f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.
Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationWITH the recent widescale deployment of Phasor Measurement. Using Effective Generator Impedance for Forced Oscillation Source Location
1 Using Effective Generator Impeance for Force Oscillation Source Location Samuel Chevalier, Stuent Member, IEEE, Petr Vorobev, Member, IEEE, Konstantin Turitsyn, Member, IEEE arxiv:178.1893v cs.sy 4 May
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationDirect Computation of Generator Internal Dynamic States from Terminal Measurements
Direct Computation of Generator nternal Dynamic States from Terminal Measurements aithianathan enkatasubramanian Rajesh G. Kavasseri School of Electrical En. an Computer Science Dept. of Electrical an
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationGyroscopic matrices of the right beams and the discs
Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More informationSection 2.1 The Derivative and the Tangent Line Problem
Chapter 2 Differentiation Course Number Section 2.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationOrdinary Differential Equations: Homework 1
Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationExam 2 Review Solutions
Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon
More informationFebruary 21 Math 1190 sec. 63 Spring 2017
February 21 Math 1190 sec. 63 Spring 2017 Chapter 2: Derivatives Let s recall the efinitions an erivative rules we have so far: Let s assume that y = f (x) is a function with c in it s omain. The erivative
More informationElectric Power Systems Research
Electric Power Systems Research 84 (22 35 43 Contents lists available at SciVerse ScienceDirect Electric Power Systems Research jou rn al h om epage: www.elsevier.com/locate/epsr Observer-base nonlinear
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 informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationmax clear Figure 1 below illustrates a stable case. P post (f) (g) (e) A 2 (d) (c)
Stability 4. Introduction In the previous notes (Stability 3, we developed the equal area iterion, which says that For stability, A =A, which means the decelerating energy (A must equal the accelerating
More informationThe dynamics of the simple pendulum
.,, 9 G. Voyatzis, ept. of Physics, University of hessaloniki he ynamics of the simple penulum Analytic methos of Mechanics + Computations with Mathematica Outline. he mathematical escription of the moel.
More informationSELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc.
SELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc. 1.0 SYSTEM DESCRIPTION The sel-erecting rotary motion inverte penulum consists of a rotary servo motor system (SRV-02) which rives
More information