Clustered Model Reduction of Interconnected Second-Order Systems and Its Applications to Power Systems

Similar documents
Clustering-based State Aggregation of Dynamical Networks

Model Reduction of Multi-Input Dynamical Networks based on Clusterwise Controllability

Network Clustering for SISO Linear Dynamical Networks via Reaction-Diffusion Transformation

Glocal Control for Network Systems via Hierarchical State-Space Expansion

Structure preserving model reduction of network systems

Eigenstructure Analysis from Symmetrical Graph Motives with Application to Aggregated Controller Design

Retrofitting State Feedback Control of Networked Nonlinear Systems Based on Hierarchical Expansion

Network Structure Preserving Model Reduction with Weak A Priori Structural Information

Clustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents

Projective State Observers for Large-Scale Linear Systems

Average State Observers for Large-Scale Network Systems

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

QUANTIZED SYSTEMS AND CONTROL. Daniel Liberzon. DISC HS, June Dept. of Electrical & Computer Eng., Univ. of Illinois at Urbana-Champaign

Model reduction of large-scale dynamical systems

Hierarchical Decentralized Observer Design for Linearly Coupled Network Systems

Outline. Input to state Stability. Nonlinear Realization. Recall: _ Space. _ Space: Space of all piecewise continuous functions

Glocal Control for Hierarchical Dynamical Systems

POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS

Consider the following example of a linear system:

4.8 Arnoldi Iteration, Krylov Subspaces and GMRES

Retrofit Control: Localization of Controller Design and Implementation

Sparsity-Promoting Optimal Control of Distributed Systems

Hierarchical Decentralized Observers for Networked Linear Systems

Research Reports on Mathematical and Computing Sciences

Discrete Time Coupled Logistic Equations with Symmetric Dispersal

Approximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method

Order Reduction of a Distributed Parameter PEM Fuel Cell Model

Parallel implementation of primal-dual interior-point methods for semidefinite programs

Solving Large Nonlinear Sparse Systems

Paul Heckbert. Computer Science Department Carnegie Mellon University. 26 Sept B - Introduction to Scientific Computing 1

Planning of Optimal Daily Power Generation Tolerating Prediction Uncertainty of Demand and Photovoltaics

Encapsulating Urban Traffic Rhythms into Road Networks

Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS

Matrix Decomposition in Privacy-Preserving Data Mining JUN ZHANG DEPARTMENT OF COMPUTER SCIENCE UNIVERSITY OF KENTUCKY

Modeling with Itô Stochastic Differential Equations

Distributed Randomized Algorithms for the PageRank Computation Hideaki Ishii, Member, IEEE, and Roberto Tempo, Fellow, IEEE

EE5900 Spring Lecture 5 IC interconnect model order reduction Zhuo Feng

AMS526: Numerical Analysis I (Numerical Linear Algebra)

Graph-based codes for flash memory

Survey of Synchronization Part I: Kuramoto Oscillators

Scientific Computing with Case Studies SIAM Press, Lecture Notes for Unit VII Sparse Matrix

Feedback stabilisation with positive control of dissipative compartmental systems

Detecting Wormhole Attacks in Wireless Networks Using Local Neighborhood Information

Projection of state space realizations

Simulation of quantum computers with probabilistic models

Dissipativity. Outline. Motivation. Dissipative Systems. M. Sami Fadali EBME Dept., UNR

Model Reduction for Unstable Systems

Minimizing Total Delay in Fixed-Time Controlled Traffic Networks

Design of structured optimal feedback gains for interconnected systems

Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems

New Directions in Computer Science

Distributed Optimization over Networks Gossip-Based Algorithms

Iterative Methods for Linear Systems

Bisimilar Finite Abstractions of Interconnected Systems

Control and synchronization in systems coupled via a complex network

AN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL MOTIVATING EXAMPLE INVERTED PENDULUM

Control Systems Design

Discrete-state Abstractions of Nonlinear Systems Using Multi-resolution Quantizer

Math 405: Numerical Methods for Differential Equations 2016 W1 Topics 10: Matrix Eigenvalues and the Symmetric QR Algorithm

NOWADAYS, many control applications have some control

Synchronization of a General Delayed Complex Dynamical Network via Adaptive Feedback

Facebook Friends! and Matrix Functions

Kronecker Product of Networked Systems and their Approximates

Model reduction of interconnected systems

Manifold Coarse Graining for Online Semi-supervised Learning

Ateneo de Manila, Philippines

Copyrighted Material. 1.1 Large-Scale Interconnected Dynamical Systems

Boundary Value Problems and Iterative Methods for Linear Systems

Introduction to Model Order Reduction

QR FACTORIZATIONS USING A RESTRICTED SET OF ROTATIONS

This paper presents the

We showed that adding a vector to a basis produces a linearly dependent set of vectors; more is true.

Model reduction of coupled systems

A Control-Theoretic Perspective on the Design of Distributed Agreement Protocols, Part

Centralized Supplementary Controller to Stabilize an Islanded AC Microgrid

SOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS. Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA

Lyapunov small-gain theorems for not necessarily ISS hybrid systems

Robust solution of Poisson-like problems with aggregation-based AMG

Black Box Linear Algebra

Sieving for Shortest Vectors in Ideal Lattices:

Using Model Reduction techniques for simulating the heat transfer in electronic systems.

Lecture 11 FIR Filters

Graph-Theoretic Analysis of Power Systems

AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning

Power System Reliability Monitoring and Control. for Transient Stability

Observations on the Stability Properties of Cooperative Systems

Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White

Spatial Autocorrelation (2) Spatial Weights

EE 381V: Large Scale Learning Spring Lecture 16 March 7

Consensus Protocols for Networks of Dynamic Agents

Pseudospectra and Nonnormal Dynamical Systems

OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU

4 Second-Order Systems

Zeros and zero dynamics

March 27 Math 3260 sec. 56 Spring 2018

Preface to the Second Edition. Preface to the First Edition

Model reduction for linear systems by balancing

Large-scale eigenvalue problems

Power Grid State Estimation after a Cyber-Physical Attack under the AC Power Flow Model

Transcription:

Clustered Model Reduction of Interconnected Second-Order Systems and Its Applications to Power Systems Takayuki Ishizaki, Jun-ichi Imura (Tokyo Inst. of Tech.)

Outline Introduction: Why clustered model reduction? Clustered Model Reduction Theory Interconnected first-order systems Extension to second-order networks Application to power networks Conclusion 2/21

3/21 Control of Large-Scale Networks Power Networks Tokyo area: 20 million houses Instability may be caused by renewables such as PV, wind Micro Grid Power line homes PV Thermal power Model reduction is one prospective approach Traffic Networks Center of Tokyo area: 5 million cars Heavy traffic jam homes Hydroelectric power How to manage?

Standard Model Reduction Framework Main goal: Find such that is small enough + stability of error systems, error analysis, low computation cost Standard methods: Balanced truncation, Hankel norm approximation error bound, stability preservation high computational cost Krylov projection lower computation cost possibly unstable model, no error bound 4/21

Application of Standard Methods to Network Systems Drawback: Network structure is lost through reduction Network system Reduced model Dense matrix Sparse Dense 5/21

Clustered Model Reduction Network system Aggregated model Cluster state Aggregated state : row vector Sparse Sparse Preservation of network structure among clusters 6/21

7/21 Why Clustered Model Reduction? Gene Network [Mochizuki et al., J. Theoretical Biology (2010)] Clustered model reduction Extract essential structure to study mechanism of functions

Outline Introduction: Why clustered model reduction? Clustered Model Reduction Theory Interconnected first-order systems Extension to second-order networks Application to power networks Conclusion 8/21

9/21 System Description (First-Order Subsystems) [Definition] Bidirectional Network with is said to be bidirectional network if and is symmetric and stable. Reaction-diffusion systems:

Clustered Model Reduction Problem [Problem] Given, find a cluster set such that where and Bidirectional network Aggregated model Cluster state Aggregated state Sparse Sparse 10/21

11/21 How to Formulate Reducibility? Bidirectional network [State trajectory under random ] 7 clusters 50 trajectories 50 nodes, nonzero is randomly chosen from [Definition] Reducible cluster A cluster can be aggregated into Local 7-dim. uncontrollability! variable is said to be reducible if under any input signal

Positive Tridiagonalization [Lemma] For every bidirectional network such that, there exists a unitary has the following structure. Bidirectional network Positive tridiagonal realization (not necessarily positive) Metzler 12/21

13/21 Reducibility Characterization Bidirectional network reducible reducible -(graph Laplacian+diagonal)

Reducibility Characterization Bidirectional network : positive tridiagonal realization : transformation matrix Index matrix reducible reducible Characterization in frequency domain identical identical Equivalent characterization of cluster reducibility 14/21

-Reducible Cluster Aggregation [Definition] -Reducibility of Clusters A cluster is said to be -reducible if Similar behavior [Theorem] If all clusters are -reducible, then where : coarseness parameter 1000 nodes 47 clusters Extract essential aggregation cluster structure! about 5% error 15/21

Outline Introduction: Why clustered model reduction? Clustered Model Reduction Theory Interconnected first-order systems Extension to second-order networks Application to power networks Conclusion 16/21

17/21 Formulation by Second-Order Systems Second-order networks Aggregated model [Problem] Given, find a cluster set such that where and

Extension to Second-Order Networks First-order representation (2n-dim. system) where Index matrix w.r.t. position velocity [Definition] A cluster is said to be -reducible if [Theorem] If all clusters are -reducible, then where 18/21

Outline Introduction: Why clustered model reduction? Clustered Model Reduction Theory Interconnected first-order systems Extension to second-order networks Application to power networks Conclusion 19/21

20/21 Numerical Example Power network modeled by swing equation Original network (300 nodes) Agg. model (42 clusters) Zoom up Position trajectory of a mass smaller error Relative error when

Concluding Remarks Clustered model reduction extract essential information on input-to-state mapping application to power networks model by swing equation Future works application to more realistic power networks extension to nonlinear systems application to control system design [T. Ishizaki et al. IEEE TAC (2014)], [T. Ishizaki et al. NOLTA (2015)], My website, etc. Thank you for your attention! 21/21