4.3 Z BUS formation considering mutual coupling between elements

Size: px
Start display at page:

Download "4.3 Z BUS formation considering mutual coupling between elements"

Transcription

1 4.3 Z BUS formation considering mutual coupling between elements Assume that the bus impedance matrix [ Z m BUS ] is known for a partial network of m nodes and a reference node. The bus voltage and bus current relation for the partial network, shown in Fig. 4.26, can be expressed as: Figure 4.26: Partial Network with m-buses [ V BUS ] = [ Z m BUS ] [Ī BUS ] (4.33) In equation (4.33), V BUS is (m 1) bus voltage vector Ī BUS is (m 1) bus current vector Z m BUS is (m m) bus impedance matrix The new added element p-q may be a branch or may be a link as discussed in the previous algorithm. 128

2 4.3.1 Addition of a branch to this partial network: The performance equation of the network with an added branch p-q is V 1 V 2 V p V m V q = Z 11 Z12 Z 1p Z 1m Z1q Z 21 Z22 Z 2p Z 2m Z2q Z p1 Zp2 Z pp Z pm Zpq Z m1 Zm2 Z mp Z mm Zmq Z q1 Zq2 Z qp Z qm Zqq Ī 1 Ī2 Ī p Ī mīq (4.34) The network is assumed to contain bilateral passive elements and hence, Zqi = Z iq for i = 1, 2, m, i q. The added branch p-q is assumed to be mutually coupled with one or more elements of the partial network. To determine element Z qi, inject a current at i th node and calculate the voltage at q th node with respect to reference, as shown in Fig Figure 4.27: Calculation of Z qi for the addition of branch Calculation of Z qi As all other bus currents are zero, bus voltages can be written as, 129

3 V 1 = Z 1i Ī i V 2 = Z 2i Ī i V p = Z pi Ī i (4.35) V m = Z mi I i V q = Z qi Ī i Also from Fig. 4.27, V p and V q can be related as, V q = V p v pq (4.36) Where v pq is the voltage across the added element p q. Also, the currents in the elements of the network can be related to the voltages across the elements as, [īpq ȳ ρσ,pq ī ρσ ] = [ ȳpq,pq ȳ pq,ρσ ] [ v pq ] (4.37) ȳ ρσ,ρσ v ρσ Where, ī pq = the current through the added element p q. ī ρσ = (m 1) current vector of the elements of the partial network. v ρσ = (m 1) voltage vector of the elements of the partial network. = Self-admittance of the added element. ȳ pq,ρσ = (m 1) vector of mutual admittances between the added element p q and the elements ρ σ of the partial network. ȳ ρσ,ρσ = (m m) primitive admittance matrix of the partial network. ȳ ρσ,pq = [ȳ pq,ρσ ] T The diagonal elements of primitive impedance matrix [ z] are the self impedance of the individual elements while the off-diagonal elements are mutual impedances between the elements. The inverse of primitive impedance matrix is the primitive admittance matrix [ȳ]. This can be explained with the help of an illustrative example. A single line diagram of a power system is shown in Fig The self-impedances of lines are written by the side of the line and are in p.u.. The two lines between nodes 1 and 3 are mutually coupled with a mutual impedance of j0.10. The primitive impedance matrix for the network can be written as, j j [ z] = 0 0 j j0.25 j j1.0 j

4 Figure 4.28: Sample Power System The inverse of [ z] is [ȳ], the primitive admittance matrix. j j [ȳ] = 0 0 j j5.0 j j2.5 j6.25 The current ī pq in the added branch p q equal to zero as node q is open. ī pq = 0 (4.38) The voltage v pq, however, is not zero as the added branch is mutually coupled to one or more elements of the partial network. Thus, the voltage across other elements of the network can be expressed as, v ρσ = V ρ V σ (4.39) where V ρ and V σ are the voltages of the nodes of the partial network. With ī pq = 0 from equation (4.37) one can write, v pq + ȳ pq,ρσ v ρσ = 0 Hence, v pq = ȳpq,ρσ v ρσ (4.40) 131

5 Substituting v ρσ from equation (4.39) and v pq equation (4.36) in equation (4.40) one gets, V q = V p + ȳpq,ρσ( V ρ V σ ) Substitution of Īi = 1 pu in equation (4.35) results in V p, V q, V ρ and V σ being replaced by their corresponding impedances and hence, Z qi = Z pi + ȳpq,ρσ( Z ρi Z σi ) (4.41) i = 1, 2,, m, i q For calculating the self impedance Z qq, a current Īq = 1 p.u. is injected into q th node with all other currents equal to zero as shown in Fig Then the voltages of the nodes are calculated from equation (4.35), as V 1 = Z 1q Ī q V 2 = Z 2q Ī q V p = Z pq Ī q (4.42) V m = Z mq I q V q = Z qq Ī q With Īq = 1 p.u., Zqq Can be calculated directly by calculating V q. also, V q = V p v pq (4.43) and ī pq = Īq = 1 (4.44) Hence, from equation (4.37) one gets ī pq = 1 = v pq + ȳ pq,ρσ v ρσ (4.45) And thus v pq can be written as, Substituting v pq and v ρσ, the above equation can be rewritten as, v pq = 1 + ȳ pq,ρσ v ρσ (4.46) V q = V p ȳ pq,ρσ( V ρ V σ ) y pq,pq (4.47) 132

6 Figure 4.29: Calculation of Z qq for the addition of a branch With Īq = 1 p.u., from equation (4.35) V p, V q, V ρ and V σ can be replaced by respective transfer impedances, Z qq = Z pq ȳ pq,ρσ( Z ρq Z σq ) y pq,pq (4.48) Addition of a link to this partial network: If the added element p q is a link, then a fictitious node l is created by connecting an ideal voltage in series with the added element, as shown in Fig The value of the source voltage is selected such that the current Īl through the added link is zero. If ē l is the voltage of node l with respect to node q and Īl is the current injected into node l from node q. The performance equation of the partial network with the added link p l and ideal series voltage source e l is, V 1 Z Z1p... Z1m Z1l V p = Z p1... Zpp... Zpm Zpl V m Z m1... Zmp... Zmm Zml ē l Z l1... Zlp... Zlm Zll Ī 1 Ī p Ī mīl (4.49) 133

7 Figure 4.30: Calculation of Z li for the addition of a link Here, Zli represents the transfer impedance relating the current, I i, injected into the i th bus and the voltage of the added source ē l, connected between nodes l and q. They are conceptually different from the elements Z ij of Z BUS matrix which relate the current injected into the j th bus and the voltage of the i th bus with repect to the reference. As ē l = V l V q, The elements Z li, 1 = 1, 2, m, i l, of the added row and column, can be calculated by injecting a current Īi into the i th node and determining the voltage of l th node with repect to q th node. Hence, Z li = ēl Ī i, Īk = 0, k = 1, 2, m, k l Also, ē l = V p V q v pl (4.50) The current ī pl through the link can be written as, ī pl = ȳ pl,pl v pl + ȳ pl,pσ v ρσ (4.51) Since the current through the link is zero, ī pl = ī pq = 0 Hence, v pl = ȳpl,ρσ v ρσ ȳ pl,pl (4.52) 134

8 Since the voltage source is ideal source, one can write, ȳ pl,ρσ = ȳ pq,ρσ and So, ȳ pl,pl = (4.53) v pl = ȳpq,ρσ v ρσ = ȳpq,ρσ( V ρ V σ ) (4.54) With I i = 1 p.u., substituting V p, Vq, Vp and V σ from equation (4.35) and v pl from equation (4.52) in equation (4.50) one gets Z li = Z pi Z qi + ȳpq,ρσ( Z pi Z σi ) (4.55) i = 1, 2, m, i l To calculate Z ll, a current is injected at the l th node with respect to node q, as shown in Fig As all other node currents are zero, the node voltages can be written as, Figure 4.31: Calculation of Z ll for the addition of a link V k = Z kl Ī l, k = 1, 2, m 135

9 ē l = Z ll Ī l (4.56) With Īl = 1 p.u., Z ll can be directly computed by calculating ē l. The current in the element p l is ī pl = Īl = 1 p.u. (4.57) From equation (4.37) one gets ī pl = ȳ pl,pl v pl + ȳ pl,pσ v pσ = 1 Further as ȳ pl,pl = and ȳ pl,pσ = ȳ pq,pσ, hence v pl can be expressed as Substituting v pl from equation (4.50), one can write v pl = 1 + ȳ pq,pσ v pσ (4.58) ē l = V p V q ȳ pq,pσ( V p V σ ) With Īl = 1 p.u., substituting V p, V q, V ρ, V σ and ē l from equation (4.56), Zll is obtained as Z ll = Z pi Z qi ȳ pq,ρσ( Z pl Z σl ) (4.59) i = 1, 2, m, i l In the case of link addition the additional row and column corresponding to fictitious node l are to be eliminated. For this the fictitious series voltage source ē l is short circuited. From equation (4.49) the bus voltages can be written in compact from as [ V m Bus] = [ Z m Bus][Ī m Bus] + [ Z][Ī l ] (4.60) Where, [ Z] = [ Z1l Z2l Z ml ] T is an (m 1) vector comprising of the entries of the column added to the Z m Bus matrix [Ī m ], [ Bus V Bus m ] = (m 1) bus current and voltage vectors respectively, of the partial network before the addition of element p l. [ Z m ] = (m m) [ Z Bus Bus ] matrix of the partial network before the addition of element p l. Ī l = current injected in the link. Also, [ē l ] = [ Z] T [Ī Bus ] + Z ll Ī l = 0 (4.61) 136

10 On substituting Īl from equation (4.61) into equation (4.60) [ V Bus ] can be written as, [ V Bus ] = [ Z m Bus] [ Z. Z T Z ll ][ĪBus] (4.62) Hence, the final [ Z Bus ] which is (m m) in size can be written as [ Z Final Bus ] = [ Z m Bus] [ Z. ZT ] (4.63) Z ll An example illustrating the [ Z Bus ] building procedure will be discussed in the next lecture. 137

4.1.3 Addition of a link between an existing node and the reference node (Case 3):

4.1.3 Addition of a link between an existing node and the reference node (Case 3): 4.1.3 Addition of a link between an existing node and the reference node (Case 3): When an element is connected between an existing node and the reference, it creates a loop and thus, the addition of this

More information

INSTITUTE OF AERONAUTICAL ENGINEERING

INSTITUTE OF AERONAUTICAL ENGINEERING INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 043 Department of Electrical and Electronics Engineering COMPUTER METHODS IN POWER SYSTEMS B.Tech III-II semester LECTURE NOTES

More information

NETWORK MATRICES. voltages through the branch impedance matrix and branch admittance matrix: ELOOP = ZLOOP ILOOP ILOOP = YLOOP ELOOP (11)

NETWORK MATRICES. voltages through the branch impedance matrix and branch admittance matrix: ELOOP = ZLOOP ILOOP ILOOP = YLOOP ELOOP (11) NETWORK MATRICES 2. FORMATION OF Y BUS AND Z BUS The bus admittance matrix, YBUS plays a very important role in computer aided power system analysis. It can be formed in practice by either of the methods

More information

Power System Analysis

Power System Analysis Power System Analysis BY A. P U R N A C H A N D E R A S S I S T A N T P R O F E S S O R D E P A R T M E N T O F E E E A C E E N G I N E E R I N G C O L L E G E Course Objectives: 1. To understand and develop

More information

4.10 Unbalanced fault analysis using Z BUS matrix:

4.10 Unbalanced fault analysis using Z BUS matrix: 4.10 Unbalanced fault analysis using Z BUS matrix: In the previous section, it is observed that, for fault calculations the Thevenin s equivalent networs, at the fault point, are needed for the three sequence

More information

2 NETWORK FORMULATION

2 NETWORK FORMULATION NTWRK FRMUATN NTRDUCTRY RMARKS For performing any power system studies on the digital computer, the first step is to construct a suitable mathematical model of the power system network The mathematical

More information

ELEC4612 Power System Analysis Power Flow Analysis

ELEC4612 Power System Analysis Power Flow Analysis ELEC462 Power Sstem Analsis Power Flow Analsis Dr Jaashri Ravishankar jaashri.ravishankar@unsw.edu.au Busbars The meeting point of various components of a PS is called bus. The bus or busbar is a conductor

More information

NETWORK CALCULATIONS updated 11/5/13 1:02 PM

NETWORK CALCULATIONS updated 11/5/13 1:02 PM NETWORK CALCULATIONS updated 11/5/13 1:02 PM 11/5/13 Network Calcula2ons (c) 2013 H. Zmuda 1 Introductory Comments The typical power transmission network span a large geographic area and involve a large

More information

2 Power System Network Matrices I

2 Power System Network Matrices I Power System Analysis Power System Network Matrices I. INRODUCION he various terms used in Graph heory are presented in this chapter. Formulation of different network matrices are discussed. Primitive

More information

= V I = Bus Admittance Matrix. Chapter 6: Power Flow. Constructing Ybus. Example. Network Solution. Triangular factorization. Let

= V I = Bus Admittance Matrix. Chapter 6: Power Flow. Constructing Ybus. Example. Network Solution. Triangular factorization. Let Chapter 6: Power Flow Network Matrices Network Solutions Newton-Raphson Method Fast Decoupled Method Bus Admittance Matri Let I = vector of currents injected into nodes V = vector of node voltages Y bus

More information

EE 581 Power Systems. Admittance Matrix: Development, Direct and Iterative solutions

EE 581 Power Systems. Admittance Matrix: Development, Direct and Iterative solutions EE 581 Power Systems Admittance Matrix: Development, Direct and Iterative solutions Overview and HW # 8 Chapter 2.4 Chapter 6.4 Chapter 6.1-6.3 Homework: Special Problem 1 and 2 (see handout) Overview

More information

Chapter 9 Balanced Faults, Part II. 9.4 Systematic Fault Analysis Using Bus Impedance Matrix

Chapter 9 Balanced Faults, Part II. 9.4 Systematic Fault Analysis Using Bus Impedance Matrix Chapter 9 Balanced Faults, Part II 9.4 Systematic Fault Analysis Using Bus Impedance Matrix In the previous analysis we employed the Thevenin model and ound the Thevenin voltage and impedance by means

More information

CHAPTER 2 LOAD FLOW ANALYSIS FOR RADIAL DISTRIBUTION SYSTEM

CHAPTER 2 LOAD FLOW ANALYSIS FOR RADIAL DISTRIBUTION SYSTEM 16 CHAPTER 2 LOAD FLOW ANALYSIS FOR RADIAL DISTRIBUTION SYSTEM 2.1 INTRODUCTION Load flow analysis of power system network is used to determine the steady state solution for a given set of bus loading

More information

TWO PORT NETWORKS Introduction: A port is normally referred to a pair of terminals of a network through which we can have access to network either for a source for measuring an output We have already seen

More information

PowerApps Optimal Power Flow Formulation

PowerApps Optimal Power Flow Formulation PowerApps Optimal Power Flow Formulation Page1 Table of Contents 1 OPF Problem Statement... 3 1.1 Vector u... 3 1.1.1 Costs Associated with Vector [u] for Economic Dispatch... 4 1.1.2 Costs Associated

More information

Improvement of Sparse Computation Application in Power System Short Circuit Study

Improvement of Sparse Computation Application in Power System Short Circuit Study Volume 44, Number 1, 2003 3 Improvement o Sparse Computation Application in Power System Short Circuit Study A. MEGA *, M. BELKACEMI * and J.M. KAUFFMANN ** * Research Laboratory LEB, L2ES Department o

More information

EE2351 POWER SYSTEM ANALYSIS

EE2351 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 information

(1) The open and short circuit required for the Z and Y parameters cannot usually be

(1) The open and short circuit required for the Z and Y parameters cannot usually be ECE 580 Network Theory Scattering Matrix 76 The Scattering Matrix Motivation for introducing the SM: () The open and short circuit required for the Z and Y parameters cannot usually be implemented in actual

More information

CHAPTER 2 DYNAMIC STABILITY MODEL OF THE POWER SYSTEM

CHAPTER 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 information

SECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems

SECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems SECTION 5: POWER FLOW ESE 470 Energy Distribution Systems 2 Introduction Nodal Analysis 3 Consider the following circuit Three voltage sources VV sss, VV sss, VV sss Generic branch impedances Could be

More information

EE 6501 POWER SYSTEMS UNIT I INTRODUCTION

EE 6501 POWER SYSTEMS UNIT I INTRODUCTION EE 6501 POWER SYSTEMS UNIT I INTRODUCTION PART A (2 MARKS) 1. What is single line diagram? A Single line diagram is diagrammatic representation of power system in which the components are represented by

More information

Chapter 8: Unsymmetrical Faults

Chapter 8: Unsymmetrical Faults Chapter 8: Unsymmetrical Faults Introduction The sequence circuits and the sequence networks developed in the previous chapter will now be used for finding out fault current during unsymmetrical faults.

More information

IMPEDANCE and NETWORKS. Transformers. Networks. A method of analysing complex networks. Y-parameters and S-parameters

IMPEDANCE and NETWORKS. Transformers. Networks. A method of analysing complex networks. Y-parameters and S-parameters IMPEDANCE and NETWORKS Transformers Networks A method of analysing complex networks Y-parameters and S-parameters 1 ENGN4545/ENGN6545: Radiofrequency Engineering L#7 Transformers Combining the effects

More information

Weighted Least Squares Topology Error Detection And Identification

Weighted Least Squares Topology Error Detection And Identification Weighted Least Squares Topology Error Detection And Identification A THESIS SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Jason Glen Lindquist IN PARTIAL FULFILLMENT

More information

Last Comments on Short-Circuit Analysis

Last Comments on Short-Circuit Analysis Last Comments on Short-Circuit Analysis.0 Exam summary Work all HW problems and know them; Read notes Read related book sections Topics:. Symmetrical faults (no book section) 2. Zbus building (9.3-9.6)

More information

Power Grid Partitioning: Static and Dynamic Approaches

Power 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 information

Power system model. Olof Samuelsson. EIEN15 Electric Power Systems L2

Power system model. Olof Samuelsson. EIEN15 Electric Power Systems L2 Power system model Olof Samuelsson EIEN15 Electric Power Systems L2 1 Outline Previously: Models for lines, generator, power electronic converter, transformer Single line diagram Per unit Bus admittance

More information

Adjoint networks and other elements of circuit theory. E416 4.Adjoint networks

Adjoint networks and other elements of circuit theory. E416 4.Adjoint networks djoint networks and other elements of circuit theory One-port reciprocal networks one-port network is reciprocal if: V I I V = Where and are two different tests on the element Example: a linear impedance

More information

Power system model. Olof Samuelsson. EIEN15 Electric Power Systems L2 1

Power system model. Olof Samuelsson. EIEN15 Electric Power Systems L2 1 Power system model Olof Samuelsson 1 Outline Previously: Models for lines, generator, power electronic converter, transformer Single line diagram Per unit Bus admittance matrix Bus impedance matrix Thévenin

More information

Module 3 : Sequence Components and Fault Analysis

Module 3 : Sequence Components and Fault Analysis Module 3 : Sequence Components and Fault Analysis Lecture 12 : Sequence Modeling of Power Apparatus Objectives In this lecture we will discuss Per unit calculation and its advantages. Modeling aspects

More information

J. Electrical Systems x-x (2010): x-xx. Regular paper

J. Electrical Systems x-x (2010): x-xx. Regular paper JBV Subrahmanyam Radhakrishna.C J. Electrical Systems x-x (2010): x-xx Regular paper A novel approach for Optimal Capacitor location and sizing in Unbalanced Radial Distribution Network for loss minimization

More information

ECEN 615 Methods of Electric Power Systems Analysis Lecture 19: State Estimation

ECEN 615 Methods of Electric Power Systems Analysis Lecture 19: State Estimation ECEN 615 Methods of Electric Power Systems Analysis Lecture 19: State Estimation Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University overbye@tamu.edu Announcements Homework

More information

ANALYSIS OF TRANSMISSION SYSTEM FAULTS IN THE PHASE DOMAIN. A Thesis JUN ZHU

ANALYSIS OF TRANSMISSION SYSTEM FAULTS IN THE PHASE DOMAIN. A Thesis JUN ZHU ANALSS OF TRANSMSSON SSTEM FAULTS N THE PHASE DOMAN A Thesis by JUN ZHU Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of

More information

State Estimation and Power Flow Analysis of Power Systems

State Estimation and Power Flow Analysis of Power Systems JOURNAL OF COMPUTERS, VOL. 7, NO. 3, MARCH 01 685 State Estimation and Power Flow Analysis of Power Systems Jiaxiong Chen University of Kentucky, Lexington, Kentucky 40508 U.S.A. Email: jch@g.uky.edu Yuan

More information

Analysis of Low-Frequency Electromagnetic Devices using Finite Elements

Analysis of Low-Frequency Electromagnetic Devices using Finite Elements Analysis of Low-Frequency Electromagnetic Devices using Finite Elements Ingeniería Energética y Electromagnética Escuela de Verano de Potencia 2014 Rafael Escarela (UAM-A) Coupling Circuit and Field Systems

More information

. ffflffluary 7, 1855.

. ffflffluary 7, 1855. x B B - Y 8 B > ) - ( vv B ( v v v (B/ x< / Y 8 8 > [ x v 6 ) > ( - ) - x ( < v x { > v v q < 8 - - - 4 B ( v - / v x [ - - B v B --------- v v ( v < v v v q B v B B v?8 Y X $ v x B ( B B B B ) ( - v -

More information

CURENT Course Power System Coherency and Model Reduction

CURENT 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 information

STEADY STATE NETWORK EQUIVALENTS FOR LARGE ELECTRICAL POWER SYSTEMS THESIS. Submitted in Partial Fulfillment. of the REQUIREMENTS for the.

STEADY STATE NETWORK EQUIVALENTS FOR LARGE ELECTRICAL POWER SYSTEMS THESIS. Submitted in Partial Fulfillment. of the REQUIREMENTS for the. STEADY STATE NETWORK EQUIVALENTS FOR LARGE ELECTRICAL POWER SYSTEMS THESIS Submitted in Partial Fulfillment of the REQUIREMENTS for the Degree of MASTER OF SCIENCE (Electrical Engineering) at the POLYTECHNIC

More information

6. MESH ANALYSIS 6.1 INTRODUCTION

6. MESH ANALYSIS 6.1 INTRODUCTION 6. MESH ANALYSIS INTRODUCTION PASSIVE SIGN CONVENTION PLANAR CIRCUITS FORMATION OF MESHES ANALYSIS OF A SIMPLE CIRCUIT DETERMINANT OF A MATRIX CRAMER S RULE GAUSSIAN ELIMINATION METHOD EXAMPLES FOR MESH

More information

SHORT QUESTIONS AND ANSWERS. Year/ Semester/ Class : III/ V/ EEE Academic Year: Subject Code/ Name: EE6501/ Power System Analysis

SHORT QUESTIONS AND ANSWERS. Year/ Semester/ Class : III/ V/ EEE Academic Year: Subject Code/ Name: EE6501/ Power System Analysis Srividya colllege of Engg & Tech,Virudhunagar Sri Vidya College of Engineering And Technology Virudhunagar 626 005 Department of Electrical and Electronics Engineering QUESTION BANK SHORT QUESTIONS AND

More information

THE INVERSE BOUNDARY PROBLEM FOR GENERAL PLANAR ELECTRICAL NETWORKS

THE INVERSE BOUNDARY PROBLEM FOR GENERAL PLANAR ELECTRICAL NETWORKS THE INVERSE BOUNDARY PROBLEM FOR GENERAL PLANAR ELECTRICAL NETWORKS DMITRIY LEYKEKHMAN Abstract. We consider an electrical network where each edge is consists of resistor, inductor, and capacitor joined

More information

Math 180B Problem Set 3

Math 180B Problem Set 3 Math 180B Problem Set 3 Problem 1. (Exercise 3.1.2) Solution. By the definition of conditional probabilities we have Pr{X 2 = 1, X 3 = 1 X 1 = 0} = Pr{X 3 = 1 X 2 = 1, X 1 = 0} Pr{X 2 = 1 X 1 = 0} = P

More information

Lupe Latin squares of order odd, 3-odd, A 2 + 3B 2 with gcm(a, B) = 1

Lupe Latin squares of order odd, 3-odd, A 2 + 3B 2 with gcm(a, B) = 1 Toyama Math. J. Vol. 34(2011), 1-22 Lupe Latin squares of order odd, 3-odd, A 2 + 3B 2 with gcm(a, B) = 1 Kazuo Azukawa Abstract. Let p be an odd integer which is written as p = A 2 + 3B 2 with gcm(a,

More information

3- BASICS. YTransformation. for balanced load. \V ab 120 = \V bc. \V ab 240 = \V ca \I a 120 = \I b \I a 240 = \I c V ab I a

3- BASICS. YTransformation. for balanced load. \V ab 120 = \V bc. \V ab 240 = \V ca \I a 120 = \I b \I a 240 = \I c V ab I a 3- BASICS YTransformation for balanced load Z =3Z Y Balanced 3- Systems \V ab 10 = \V bc \V ab 40 = \V ca \I a 10 = \I b \I a 40 = \I c V ab I a = p 3 V an = p 3 I ab \V ab 30 = \V an \I ab 30 = \I a S

More information

Chapter 8 VOLTAGE STABILITY

Chapter 8 VOLTAGE STABILITY Chapter 8 VOTAGE STABIITY The small signal and transient angle stability was discussed in Chapter 6 and 7. Another stability issue which is important, other than angle stability, is voltage stability.

More information

Incorporation of Asynchronous Generators as PQ Model in Load Flow Analysis for Power Systems with Wind Generation

Incorporation of Asynchronous Generators as PQ Model in Load Flow Analysis for Power Systems with Wind Generation Incorporation of Asynchronous Generators as PQ Model in Load Flow Analysis for Power Systems with Wind Generation James Ranjith Kumar. R, Member, IEEE, Amit Jain, Member, IEEE, Power Systems Division,

More information

KINGS COLLEGE OF ENGINEERING Punalkulam

KINGS 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 information

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous)

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 043 ELECTRICAL AND ELECTRONICS ENGINEERING QUESTION BANK Course Name : Computer Methods in Power Systems Course Code : A60222

More information

EE 205 Dr. A. Zidouri. Electric Circuits II. Two-Port Circuits Two-Port Parameters. Lecture #42

EE 205 Dr. A. Zidouri. Electric Circuits II. Two-Port Circuits Two-Port Parameters. Lecture #42 EE 05 Dr. A. Zidouri Electric Circuits Two-Port Circuits Two-Port Parameters Lecture #4-1 - EE 05 Dr. A. Zidouri The material to be covered in this lecture is as follows: o ntroduction to two-port circuits

More information

EE2351 POWER SYSTEM ANALYSIS UNIT I: INTRODUCTION

EE2351 POWER SYSTEM ANALYSIS UNIT I: INTRODUCTION EE2351 POWER SYSTEM ANALYSIS UNIT I: INTRODUCTION PART: A 1. Define per unit value of an electrical quantity. Write equation for base impedance with respect to 3-phase system. 2. What is bus admittance

More information

P E R E N C O - C H R I S T M A S P A R T Y

P E R E N C O - C H R I S T M A S P A R T Y L E T T I C E L E T T I C E I S A F A M I L Y R U N C O M P A N Y S P A N N I N G T W O G E N E R A T I O N S A N D T H R E E D E C A D E S. B A S E D I N L O N D O N, W E H A V E T H E P E R F E C T R

More information

CHAPTER 6 : LITERATURE REVIEW

CHAPTER 6 : LITERATURE REVIEW CHAPTER 6 : LITERATURE REVIEW Chapter : LITERATURE REVIEW 77 M E A S U R I N G T H E E F F I C I E N C Y O F D E C I S I O N M A K I N G U N I T S A B S T R A C T A n o n l i n e a r ( n o n c o n v e

More information

Automatic Formulation of Circuit Equations

Automatic Formulation of Circuit Equations ECE 570 Session 3 IC 752-E Computer Aided Engineering for Integrated Circuits Automatic Formulation of Circuit Equations Objective: Basics of computer aided analysis/simulation Outline:. Discussion of

More information

Extensive Form Abstract Economies and Generalized Perfect Recall

Extensive Form Abstract Economies and Generalized Perfect Recall Extensive Form Abstract Economies and Generalized Perfect Recall Nicholas Butler Princeton University July 30, 2015 Nicholas Butler (Princeton) EFAE and Generalized Perfect Recall July 30, 2015 1 / 1 Motivation

More information

1 Unified Power Flow Controller (UPFC)

1 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 information

Operations Research. Duality in linear programming.

Operations Research. Duality in linear programming. Operations Research Duality in linear programming Duality in linear programming As we have seen in past lessons, linear programming are either maximization or minimization type, containing m conditions

More information

arxiv: v1 [cs.sy] 19 Dec 2018

arxiv: v1 [cs.sy] 19 Dec 2018 Saturated control without velocity measurements for planar robots with flexible joints P Borja, T Wesselink and JMA Scherpen Faculty of Science and Engineering, University of Groningen Nijenborgh 4, 9747

More information

Equivalent relaxations of optimal power flow

Equivalent relaxations of optimal power flow Equivalent relaxations of optimal power flow 1 Subhonmesh Bose 1, Steven H. Low 2,1, Thanchanok Teeraratkul 1, Babak Hassibi 1 1 Electrical Engineering, 2 Computational and Mathematical Sciences California

More information

Chapter 3. Loop and Cut-set Analysis

Chapter 3. Loop and Cut-set Analysis Chapter 3. Loop and Cut-set Analysis By: FARHAD FARADJI, Ph.D. Assistant Professor, Electrical Engineering, K.N. Toosi University of Technology http://wp.kntu.ac.ir/faradji/electriccircuits2.htm References:

More information

GENERAL FAULT ADMITTANCE METHOD LINE- TO- LINE-TO-LINE UNSYMMETRICAL FAULT

GENERAL FAULT ADMITTANCE METHOD LINE- TO- LINE-TO-LINE UNSYMMETRICAL FAULT Journal of Engineering Science and Technology ol., No. 9 (5) 3-43 School of Engineering, Taylor s University GENERAL FAULT ADMTTANCE METHOD LNE- TO- LNE-TO-LNE UNSMMETRCAL FAULT J. D. SAKALA*, J. S. J.

More information

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability

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 information

Two-Layer Network Equivalent for Electromagnetic Transients

Two-Layer Network Equivalent for Electromagnetic Transients 1328 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 18, NO. 4, OCTOBER 2003 Two-Layer Network Equivalent for Electromagnetic Transients Mohamed Abdel-Rahman, Member, IEEE, Adam Semlyen, Life Fellow, IEEE, and

More information

POWER flow studies are the cornerstone of power system

POWER flow studies are the cornerstone of power system University of Wisconsin-Madison Department of Electrical and Computer Engineering. Technical Report ECE-2-. A Sufficient Condition for Power Flow Insolvability with Applications to Voltage Stability Margins

More information

Australian Journal of Basic and Applied Sciences. General Fault Admittance Method Solution of a Balanced Line-to-Line-to-Line Fault

Australian Journal of Basic and Applied Sciences. General Fault Admittance Method Solution of a Balanced Line-to-Line-to-Line Fault Australian Journal of Basic and Applied Sciences, 8() January 4, Pages: 8-47 AENSI Journals Australian Journal of Basic and Applied Sciences Journal home page: www.ajbasweb.com General Fault Admittance

More information

Branch Flow Model. Computing + Math Sciences Electrical Engineering

Branch Flow Model. Computing + Math Sciences Electrical Engineering Branch Flow Model Masoud Farivar Steven Low Computing + Math Sciences Electrical Engineering Arpil 014 TPS paper Farivar and Low Branch flow model: relaxations and convexification IEEE Trans. Power Systems,

More information

, and ignoring all load currents, determine

, and ignoring all load currents, determine ECE43 Test 3 Dec 8, 5 Q. (33 pts.) The Zbus for the above 3-bus network with bus as reference, in per unit, is given to be 3.87 j.798 j.8 j Z.798 j.87 j.8 j bus.8 j.8 j j Assuming that the prefault values

More information

MA 201: Partial Differential Equations Lecture - 2

MA 201: Partial Differential Equations Lecture - 2 MA 201: Partial Differential Equations Lecture - 2 Linear First-Order PDEs For a PDE f(x,y,z,p,q) = 0, a solution of the type F(x,y,z,a,b) = 0 (1) which contains two arbitrary constants a and b is said

More information

Draft. Lecture 12 Gaussian Elimination and LU Factorization. MATH 562 Numerical Analysis II. Songting Luo

Draft. Lecture 12 Gaussian Elimination and LU Factorization. MATH 562 Numerical Analysis II. Songting Luo Lecture 12 Gaussian Elimination and LU Factorization Songting Luo Department of Mathematics Iowa State University MATH 562 Numerical Analysis II ongting Luo ( Department of Mathematics Iowa State University[0.5in]

More information

MATH 387 ASSIGNMENT 2

MATH 387 ASSIGNMENT 2 MATH 387 ASSIGMET 2 SAMPLE SOLUTIOS BY IBRAHIM AL BALUSHI Problem 4 A matrix A ra ik s P R nˆn is called symmetric if a ik a ki for all i, k, and is called positive definite if x T Ax ě 0 for all x P R

More information

Study Notes on Network Theorems for GATE 2017

Study Notes on Network Theorems for GATE 2017 Study Notes on Network Theorems for GATE 2017 Network Theorems is a highly important and scoring topic in GATE. This topic carries a substantial weight age in GATE. Although the Theorems might appear to

More information

Computational Intelligence Lecture 20:Neuro-Fuzzy Systems

Computational Intelligence Lecture 20:Neuro-Fuzzy Systems Computational Intelligence Lecture 20:Neuro-Fuzzy Systems Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Computational Intelligence

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 12: Gaussian Elimination and LU Factorization Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 10 Gaussian Elimination

More information

Application of Monte Carlo Simulation to Multi-Area Reliability Calculations. The NARP Model

Application of Monte Carlo Simulation to Multi-Area Reliability Calculations. The NARP Model Application of Monte Carlo Simulation to Multi-Area Reliability Calculations The NARP Model Any power system reliability model using Monte Carlo simulation consists of at least the following steps: 1.

More information

c Z ... Ln :E UJ UJ ...J. (f) (f) (f) UJ. (f) (f)

c Z ... Ln :E UJ UJ ...J. (f) (f) (f) UJ. (f) (f) - r1 N......J... co (f) Ln :E UJ UJ u UJ tn > (f) UJ S a. \ r1. c Z... (f) (f) UJ (f) -,... f.. U'....1i I /, IL I.., l..c: :! ta Q, - J "'C r.. c ::::> -u \("'t,; l I o I .\ ( X R+ IX = Zlli --...------------R

More information

Networks and Systems Prof. V. G. K. Murti Department of Electrical Engineering Indian Institution of Technology, Madras

Networks and Systems Prof. V. G. K. Murti Department of Electrical Engineering Indian Institution of Technology, Madras Networks and Systems Prof. V. G. K. Murti Department of Electrical Engineering Indian Institution of Technology, Madras Lecture - 32 Network Function (3) 2-port networks: Symmetry Equivalent networks Examples

More information

EE 742 Chapter 3: Power System in the Steady State. Y. Baghzouz

EE 742 Chapter 3: Power System in the Steady State. Y. Baghzouz EE 742 Chapter 3: Power System in the Steady State Y. Baghzouz Transmission Line Model Distributed Parameter Model: Terminal Voltage/Current Relations: Characteristic impedance: Propagation constant: π

More information

B.7 Uncertainty principle (supplementary)

B.7 Uncertainty principle (supplementary) B. BASIC CONCEPTS FROM QUANTUM THEORY 93 B.7 Uncertainty principle (supplementary) You might be surprised that the famous Heisenberg uncertainty principle is not among the postulates of quantum mechanics.

More information

A Novel Technique to Improve the Online Calculation Performance of Nonlinear Problems in DC Power Systems

A Novel Technique to Improve the Online Calculation Performance of Nonlinear Problems in DC Power Systems electronics Article A Novel Technique to Improve the Online Calculation Performance of Nonlinear Problems in DC Power Systems Qingshan Xu 1, Yuqi Wang 1, * ID, Minjian Cao 1 and Jiaqi Zheng 2 1 School

More information

EN221 - Fall HW # 1 Solutions

EN221 - Fall HW # 1 Solutions EN221 - Fall2008 - HW # 1 Solutions Prof. Vivek Shenoy 1. Let A be an arbitrary tensor. i). Show that II A = 1 2 {(tr A)2 tr A 2 } ii). Using the Cayley-Hamilton theorem, deduce that Soln. i). Method-1

More information

Distributed Optimization: Analysis and Synthesis via Circuits

Distributed Optimization: Analysis and Synthesis via Circuits Distributed Optimization: Analysis and Synthesis via Circuits Stephen Boyd Prof. S. Boyd, EE364b, Stanford University Outline canonical form for distributed convex optimization circuit intepretation primal

More information

Matrices: 2.1 Operations with Matrices

Matrices: 2.1 Operations with Matrices Goals In this chapter and section we study matrix operations: Define matrix addition Define multiplication of matrix by a scalar, to be called scalar multiplication. Define multiplication of two matrices,

More information

Optimal Placement & sizing of Distributed Generator (DG)

Optimal Placement & sizing of Distributed Generator (DG) Chapter - 5 Optimal Placement & sizing of Distributed Generator (DG) - A Single Objective Approach CHAPTER - 5 Distributed Generation (DG) for Power Loss Minimization 5. Introduction Distributed generators

More information

B.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 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 information

Analysis of Bifurcations in a Power System Model with Excitation Limits

Analysis of Bifurcations in a Power System Model with Excitation Limits Analysis of Bifurcations in a Power System Model with Excitation Limits Rajesh G. Kavasseri and K. R. Padiyar Department of Electrical Engineering Indian Institute of Science, Bangalore, India Abstract

More information

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages:

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages: CS100: DISCRETE STRUCTURES Lecture 3 Matrices Ch 3 Pages: 246-262 Matrices 2 Introduction DEFINITION 1: A matrix is a rectangular array of numbers. A matrix with m rows and n columns is called an m x n

More information

Chapter 5 continued. Chapter 5 sections

Chapter 5 continued. Chapter 5 sections Chapter 5 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions

More information

ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation

ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation Prof. om Overbye Dept. of Electrical and Computer Engineering exas A&M University overbye@tamu.edu Announcements

More information

Power system modelling under the phasor approximation

Power 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 information

An Electrical Interpretation of Mechanical Systems via the Pseudo-inductor in the Brayton-Moser Equations

An Electrical Interpretation of Mechanical Systems via the Pseudo-inductor in the Brayton-Moser Equations Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 2-5, 2005 WeC3.2 An Electrical Interpretation of Mechanical Systems via

More information

ECEN 667 Power System Stability Lecture 17: Transient Stability Solutions, Load Models

ECEN 667 Power System Stability Lecture 17: Transient Stability Solutions, Load Models ECEN 667 Power System Stability Lecture 17: Transient Stability Solutions, Load Models Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements

More information

DIRECTED GRAPH BASED POWER-FLOW ALGORITHM FOR 1-PHASE RDN

DIRECTED GRAPH BASED POWER-FLOW ALGORITHM FOR 1-PHASE RDN method is used to formulate element-currents during the Forward-Sweep and the inclusion of power variable in the algorithm during Backward Sweep leads to more accurate results. The proposed TT based power-flow

More information

Lossy DC Power Flow. John W. Simpson-Porco. Universidad de Chile Santiago, Chile. March 2, 2017

Lossy DC Power Flow. John W. Simpson-Porco. Universidad de Chile Santiago, Chile. March 2, 2017 Lossy DC Power Flow John W. Simpson-Porco Universidad de Chile Santiago, Chile March 2, 2017 Our 20th Century Bulk Power System A large-scale, nonlinear, hybrid, stochastic, distributed, cyber-physical...

More information

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for. 4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX Some reliminries: Let A be rel symmetric mtrix. Let Cos θ ; (where we choose θ π for Cos θ 4 purposes of convergence of the scheme)

More information

The Fundamental Insight

The Fundamental Insight The Fundamental Insight We will begin with a review of matrix multiplication which is needed for the development of the fundamental insight A matrix is simply an array of numbers If a given array has m

More information

11 a 12 a 13 a 21 a 22 a b 12 b 13 b 21 b 22 b b 11 a 12 + b 12 a 13 + b 13 a 21 + b 21 a 22 + b 22 a 23 + b 23

11 a 12 a 13 a 21 a 22 a b 12 b 13 b 21 b 22 b b 11 a 12 + b 12 a 13 + b 13 a 21 + b 21 a 22 + b 22 a 23 + b 23 Chapter 2 (3 3) Matrices The methods used described in the previous chapter for solving sets of linear equations are equally applicable to 3 3 matrices. The algebra becomes more drawn out for larger matrices,

More information

A DC Power Flow Extension

A DC Power Flow Extension 2013 4th IEEE PES Innovative Smart Grid Technologies Europe (ISGT Europe), October 6-9, Copenhagen 1 A DC Power Flow Extension Theodoros Kyriakidis, Rachid Cherkaoui, Maher Kayal Electronics Laboratory

More information

Lecture 16 Methods for System of Linear Equations (Linear Systems) Songting Luo. Department of Mathematics Iowa State University

Lecture 16 Methods for System of Linear Equations (Linear Systems) Songting Luo. Department of Mathematics Iowa State University Lecture 16 Methods for System of Linear Equations (Linear Systems) Songting Luo Department of Mathematics Iowa State University MATH 481 Numerical Methods for Differential Equations Songting Luo ( Department

More information

A COMPUTER PROGRAM FOR SHORT CIRCUIT ANALYSIS OF ELECTRIC POWER SYSTEMS

A COMPUTER PROGRAM FOR SHORT CIRCUIT ANALYSIS OF ELECTRIC POWER SYSTEMS NIJOTECH VOL. 5 NO. 1 MARCH 1981 EJEBE 46 A COMPUTER PROGRAM FOR SHORT CIRCUIT ANALYSIS OF ELECTRIC POWER SYSTEMS BY G.C. EJEBE DEPARTMENT OF ELECTRICAL/ELECTRONIC ENGINEERING UNIVERSITY OF NIGERIA, NSUKKA.

More information

Lecture 7: Two-Ports (2)

Lecture 7: Two-Ports (2) Lecture 7: Two-Ports (2) Amin Arbabian Jan M. Rabaey EE142 Fall 2010 Sept. 16 th, 2010 Announcements HW3 Posted, due Thursday September 23 rd, 330pm EE142 drop box. Office hours Two-Port notes by Professor

More information

Effect of Various Holomorphic Embeddings on Convergence Rate and Condition Number as Applied to the Power Flow Problem

Effect of Various Holomorphic Embeddings on Convergence Rate and Condition Number as Applied to the Power Flow Problem Effect of Various Holomorphic Embeddings on Convergence Rate and Condition Number as Applied to the Power Flow Problem Yuting Li Committee members: Dr. Daniel J. Tylavsky, Chair Dr. John Undrill Dr. Vijay

More information