Model reduction of large-scale systems
|
|
- Gary Wilson
- 5 years ago
- Views:
Transcription
1 Model reduction of large-scale systems An overview Thanos Antoulas joint work with Dan Sorensen URL: www-ece.rice.edu/ aca CITI Lecture, Rice, February 2005 Model reduction of large-scale systems p.1
2 Support: Collaborators: NSF + NSF ITR Grants Dan Sorensen (CAAM), TA (ECE) Paul van Dooren (Belgium) Ahmed Sameh (Purdue) Mete Sozen (Purdue) Ananth Grama (Purdue) Chris Hoffmann (Purdue) Kyle Gallivan (Florida State) Serkan Gugercin (Virginia Tech) Model reduction of large-scale systems p.2
3 Contents 1. Problem statement 2. Motivating examples 3. Approximation methods (very briefly) SVD Krylov Krylov/SVD 4. Open problems/concluding remarks Model reduction of large-scale systems p.3
4 The big picture Physical System + Data discretization Modeling ODEs Model reduction reduced # of ODEs PDEs Simulation Control Model reduction of large-scale systems p.4
5 Dynamical systems y 1 ( ) y 2 ( ). y p ( ) x 1 ( ) x 2 ( ). x n ( ) u 1 ( ) u 2 ( ). u m ( ) We consider explicit state equations Σ : ẋ(t) = f(x(t),u(t)), y(t) = h(x(t),u(t)) with state x( ) of dimension n m,p. Model reduction of large-scale systems p.5
6 Problem statement Given: dynamical system Σ = (f,h) with: u(t) R m, x(t) R n, y(t) R p. Problem: Approximate Σ with: ˆΣ = ( ˆf,ĥ), u(t) Rm, ˆx(t) R k, y(t) R p, k n : (1) Approximation error small - global error bound (2) Preservation of stability, passivity (3) Procedure computationally efficient Model reduction of large-scale systems p.6
7 Approximation by projection Unifying feature of approximation methods: projections. Let V, W R n k, such that W V = I k Π = VW is a projection. Define ˆx = W x. Then ˆΣ : d dt ˆx(t) = W f(vˆx(t), u(t)) y(t) = h(vˆx(t), u(t)) Thus ˆΣ is "good" approximation of Σ, if x Πx is "small". Model reduction of large-scale systems p.7
8 Special case: linear dynamical systems Σ: ẋ(t) = Ax(t) + Bu(t), y(t) = Cx(t) + Du(t) Σ = A B R (n+p) (n+m) C D Problem : Approximate Σ by projection: Π = V W, Π 2 = Π, ˆΣ = Â Ĉ ˆB ˆD = W AV CV W B D R (k+p) (k+m), k n Norms : H norm: worst output error y(t) ŷ(t) for u(t) = 1. H 2 -norm: h(t) ĥ(t) Model reduction of large-scale systems p.8
9 Motivating Examples: Simulation/Control 1. Passive devices VLSI circuits Thermal issues 2. Data assimilation North sea forecast Air quality forecast Sensor placement 3. Biological/Molecular systems Honeycomb vibrations MD simulations Heat capacity 4. CVD reactor Bifurcations 5. Mechanical systems: Windscreen vibrations Buildings 6. Optimal cooling Steel profile 7. MEMS: Micro Electro-Mechanical Systems Elf sensor Model reduction of large-scale systems p.9
10 Passive devices : VLSI circuits 1960 s: IC 1971: Intel : Intel Pentium IV 10µ details 2300 components 64KHz speed 0.18µ details 42M components 2GHz speed 2km interconnect 7 layers Model reduction of large-scale systems p.10
11 Passive devices : VLSI circuits Conclusion: Simulations are required to verify that internal electromagnetic fields do not significantly delay or distort circuit signals. Therefore interconnections must be modeled. Electromagnetic modeling of packages and interconnects resulting models very complex: using PEEC methods (discretization of Maxwell s equations): n SPICE: inadequate Source: van der Meijs (Delft) Model reduction of large-scale systems p.11
12 Passive devices : Thermal management of semiconductor device Semiconductor package ANSYS grid Temperature distribution Semiconductor manufacturers provide to customers models of heat dissipation. Procedure: from finite element models generated using ANSYS, a reduced order model is derived, which is provided to the customer. Model reduction of large-scale systems p.12
13 Weather: Storm surge forecast depth Wick n North Shields Huibertgat Harlingen Delfzijl Den Helder Lowestoft IJmuiden-II HoekvanHolland Sheerness Vlissingen Dover m Wave propagation in shallow sea described by shallow water equations. Computational complexity: 60, 000 equations more than 2 days computational time. Allowed: 6 hours! Actually, there are 8 measurement stations Data assimilation Kalman Filter Propagate low-rank factors of data covariance matrix. Source: Verlaan, Heemink (Delft University) 0 Model reduction of large-scale systems p.13
14 Weather: Air quality simulations Data assimilation for atmospheric chemistry models: model the formation and transport of air pollutants (CO). Chemistry transport model (CTM): governed by PDEs + boundary conditions. After spatial discretization ODEs. Grid: x, y direction 1km = 100 points; z direction 10km = 30 points. Terra satellite pictures of CO concentration at the pacific rim Source: Glowinski, He (University of Houston) Model reduction of large-scale systems p.14
15 CVD reactors Stability analysis of steady flows as a function of Grashof, Rayleigh, Reynolds numbers. 16 million variable FEM model of CVD reactor (full 3D Navier Stokes) Used P_ARPACK on Sandia-Intel Tflop, 1024 processors. Computed 6 eigenvalues of largest real part Source: Lehoucq, Salinger (Sandia) Model reduction of large-scale systems p.15
16 Mechanical systems: cars Car windscreen simulation subject to acceleration load. Problem: compute noise at points away from the window. PDE: describes deformation of a structure of a specific material; FE discretization: 7564 nodes (3 layers of 60 by 30 elements). Material: glass with Young modulus N/m 2 ; density 2490kg/m 3 ; Poisson ratio 0.23 coefficients of FE model determined experimentally. The discretized problem has dimension: 22,692. Notice: this problem yields 2nd order equations: Mẍ(t) + Cẋ(t) + Kx(t) = f(t). Source: Meerbergen (Free Field Technologies) Model reduction of large-scale systems p.16
17 Mechanical Systems: Buildings Earthquake prevention Yokohama Landmark Tower: two active tuned mass dampers damp the dominant vibration mode of 0.185Hz. Source: S. Williams Model reduction of large-scale systems p.17
18 Mechanical Systems: Buildings Building Height Control mechanism Damping frequency Damping mass CN Tower, Toronto 533 m Passive tuned mass damper Hancock building, Boston 244 m Two passive tuned dampers 0.14Hz, 2x300t Sydney tower 305 m Passive tuned pendulum 0.1,0.5z, 220t Rokko Island P&G, Kobe 117 m Passive tuned pendulum Hz, 270t Yokohama Landmark tower 296 m Active tuned mass dampers (2) 0.185Hz, 340t Shinjuku Park Tower 296 m Active tuned mass dampers (3) 330t TYG Building, Atsugi 159 m Tuned liquid dampers (720) 0.53Hz, 18.2t Source: S. Williams Model reduction of large-scale systems p.18
19 MEMS: Elk sensor Mercedes A Class semiconductor rotation sensor Source: Laur (Bremen) Model reduction of large-scale systems p.19
20 Approximation methods: Overview Model reduction of large-scale systems p.20
21 Approximation methods: Overview SVD Nonlinear systems POD methods Empirical Gramians Linear systems Balanced truncation Hankel approximation Model reduction of large-scale systems p.20
22 Approximation methods: Overview Krylov SVD Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians Linear systems Balanced truncation Hankel approximation Model reduction of large-scale systems p.20
23 Approximation methods: Overview Krylov SVD Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians Linear systems Balanced truncation Hankel approximation Krylov/SVD Methods Model reduction of large-scale systems p.20
24 SVD Approximation methods SVD: A = UΣV. A prototype problem. Supernova Clown Singular values of Clown and Supernova Supernova: original picture Clown: original picture Clown: rank 6 approximation green: clown red: supernova (log log scale) Supernova: rank 6 approximation Supernova: rank 20 approximation Clown: rank 12 approximation Clown: rank 20 approximation Singular values provide trade-off between accuracy and complexity Model reduction of large-scale systems p.21
25 POD: Proper Orthogonal Decomposition Consider: ẋ(t) = f(x(t), u(t)), y(t) = h(x(t), u(t)). Snapshots of state: X = [x(t 1 ) x(t 2 ) x(t N )] R n N Model reduction of large-scale systems p.22
26 POD: Proper Orthogonal Decomposition Consider: ẋ(t) = f(x(t), u(t)), y(t) = h(x(t), u(t)). Snapshots of state: X = [x(t 1 ) x(t 2 ) x(t N )] R n N SVD: X = UΣV U k Σ k Vk, k n. Approximate the state: ξ(t) = U kx(t) x(t) U k ξ(t), ξ(t) R k Model reduction of large-scale systems p.22
27 POD: Proper Orthogonal Decomposition Consider: ẋ(t) = f(x(t), u(t)), y(t) = h(x(t), u(t)). Snapshots of state: X = [x(t 1 ) x(t 2 ) x(t N )] R n N SVD: X = UΣV U k Σ k Vk, k n. Approximate the state: ξ(t) = U kx(t) x(t) U k ξ(t), ξ(t) R k Project state and output equations. Reduced order system: ξ(t) = U kf(u k ξ(t), u(t)), y(t) = h(u k ξ(t), u(t)) ξ(t) evolves in a low-dimensional space. Model reduction of large-scale systems p.22
28 POD: Proper Orthogonal Decomposition Consider: ẋ(t) = f(x(t), u(t)), y(t) = h(x(t), u(t)). Snapshots of state: X = [x(t 1 ) x(t 2 ) x(t N )] R n N SVD: X = UΣV U k Σ k Vk, k n. Approximate the state: ξ(t) = U kx(t) x(t) U k ξ(t), ξ(t) R k Project state and output equations. Reduced order system: ξ(t) = Ukf(U k ξ(t), u(t)), y(t) = h(u k ξ(t), u(t)) ξ(t) evolves in a low-dimensional space. Issues with POD: (a) Choice of snapshots, (b) singular values not I/O invariants. Model reduction of large-scale systems p.22
29 The Hankel singular values Trade-off between accuracy and complexity for linear dynamical systems is provided by the Hankel Singular Values, which are computed as follows: Gramians P = 0 e At BB e A t dt, Q = 0 e A t C Ce At dt To compute gramians solve 2 Lyapunov equations : AP + PA + BB = 0, P > 0 A Q + QA + C C = 0, Q > 0 σ i = λ i (PQ) σ i : Hankel singular values of system Σ Model reduction of large-scale systems p.23
30 Approximation by balanced truncation There exists basis where P = Q = S = diag (σ 1,, σ n ). This is the balanced basis of the system. In this basis partition: A = A 11 A 12 A 21 A 22, B = B 1 B 2, C = (C 1 C 2 ), S = Σ Σ 2 ROM obtained by balanced truncation ˆΣ = A 11 B 1 C 1 Σ 2 contains the small Hankel singular values. ( Projector: Π = V W I k where V = W = 0 ). Model reduction of large-scale systems p.24
31 Properties of balanced reduction 1. Stability is preserved 2. Global error bound : σ k+1 Σ ˆΣ 2(σ k σ n ) Model reduction of large-scale systems p.25
32 Properties of balanced reduction 1. Stability is preserved 2. Global error bound : σ k+1 Σ ˆΣ 2(σ k σ n ) Drawbacks 1. Dense computations, matrix factorizations and inversions may be ill-conditioned 2. Need whole transformed system in order to truncate number of operations O(n 3 ) Model reduction of large-scale systems p.25
33 Approximation methods: Krylov Krylov Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians Krylov/SVD Methods SVD Linear systems Balanced truncation Hankel approximation Model reduction of large-scale systems p.26
34 Krylov approximation methods Given Σ = ( A C B D ), expand transfer function around s 0 : G(s) = η 0 + η 1 (s s 0 ) + η 2 (s s 0 ) 2 + η 3 (s s 0 ) 3 + Moments at s 0 : η j. Model reduction of large-scale systems p.27
35 Krylov approximation methods Given Σ = ( A C B D ), expand transfer function around s 0 : G(s) = η 0 + η 1 (s s 0 ) + η 2 (s s 0 ) 2 + η 3 (s s 0 ) 3 + Moments at s 0 : η j. Find ˆΣ = ( Â ˆB Ĉ ˆD ), with Ĝ(s) = ˆη 0 + ˆη 1 (s s 0 ) + ˆη 2 (s s 0 ) 2 + ˆη 3 (s s 0 ) 3 + such that for appropriate l: η j = ˆη j, j = 1, 2,,l Model reduction of large-scale systems p.27
36 Key fact for numerical reliability But: computation of moments is numerically problematic ( ) A B Solution: given Σ = moment matching methods can be implemented in a numerically efficient way: C moments can be matched without moment computation iterative implementation : ARNOLDI LANCZOS Special case: Expansion around infinity: η j : Markov parameters. Problem: partial realization. Expansion around zero: η j : moments. Problem: Padé interpolation. In general: arbitrary s 0 C: Problem: rational interpolation. Model reduction of large-scale systems p.28
37 Properties of Krylov methods (a) Number of operations: O(kn 2 ) or O(k 2 n) vs. O(n 3 ) efficiency (b) Only matrix-vector multiplications are required. No matrix factorizations and/or inversions. No need to compute transformed model and then truncate. Model reduction of large-scale systems p.29
38 Properties of Krylov methods (a) Number of operations: O(kn 2 ) or O(k 2 n) vs. O(n 3 ) efficiency (b) Only matrix-vector multiplications are required. No matrix factorizations and/or inversions. No need to compute transformed model and then truncate. (c) Drawbacks global error bound? ˆΣ may not be stable. Remedy: implicit restart. Lanczos breaks down if det H i = 0. Remedy: look-ahead. ˆΣ tends to approximate the high frequency poles of Σ. Remedy: match expansions around other frequencies (rational Lanczos). Model reduction of large-scale systems p.29
39 Approximation methods: Krylov/SVD Krylov Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians SVD Krylov/SVD Methods Linear systems Balanced truncation Hankel approximation Model reduction of large-scale systems p.30
40 Application: LA Administration Building 2nd order model: Mẍ + Dẋ + Kx = Bu, y = Cx, x R n, n = Proportional damping: D = αm + βk Model reduction of large-scale systems p.31
41 300 Damped Evals alpha = 0.67, beta = Freq. Response alpha = 0.67, beta = reduced original Response G( jω) Frequency ω Projection matrices: W = V, spanv = spanned by eigenvectors of 200 dominant eigenvalues ˆM = V MV, ˆD = V DV, ˆK = V KV, ˆB = V B, Ĉ = CV. Model reduction of large-scale systems p.32
42 Bowen Large-Scale Laboratory (Purdue) Test facility: Model reduction of large-scale systems p.33
43 Next application: WTC more than 1,000,000 degrees of freedom Model reduction of large-scale systems p.34
44 Projectors and complexity Projectors are a unifying feature of model reduction methods Project with Π = V W, Π 2 = Π, to get reduced model: ˆx = (W AV )ˆx + (W B)u, y = (CV )ˆx + Du Quality is measured in terms of the frequency response G(jω) = C(jωI A) 1 B + D and in particular its peak, called the H -norm, or its 2-norm (energy) called the H 2 -norm. Computational complexity: n: dimension of original system; k: dimension of reduced system; α: average number of non-zero elements per row of A Dense operations SVD: solve for gramians n 3 αkn SVD: perform balancing n 3 αkn Krylov kn 2 αkn Approximate/ sparse operations Model reduction of large-scale systems p.35
45 Software Goal: reliable solution of large dimensional numerical problems encountered in complex control applications. Turning numerically reliable algorithms into high performance, portable and robust numerical software relies on numerical linear algebra, computer science, computer hardware, etc. General purpose software including model reduction: MATLAB, SCILAB Model reduction software: SLICOT Parallel model reduction software: PSLICOT For large-scale eigenvalue problems: ARPACK SUGAR: MEMS simulation MOZART: Ozone propagation Model reduction of large-scale systems p.36
46 Approximation methods: Properties Krylov Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians SVD Krylov/SVD Methods Linear systems Balanced truncation Hankel approximation Model reduction of large-scale systems p.37
47 Approximation methods: Properties Properties Krylov Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians SVD Linear systems Balanced truncation Hankel approximation numerical efficiency n 500 Krylov/SVD Methods Model reduction of large-scale systems p.37
48 Approximation methods: Properties Properties numerical efficiency n 500 Krylov Realization Interpolation Lanczos Arnoldi Nonlinear systems POD methods Empirical Gramians Krylov/SVD Methods SVD Linear systems Balanced truncation Hankel approximation Properties Stability Error bound n 500 Model reduction of large-scale systems p.37
49 Open problems Decay rate of the Hankel singular values Choice of expansion points for rational Krylov Choice of spectral zeros in passive model reduction Iterative SVD model reduction methods Model reduction of uncertain systems Model reduction of second-order systems Model reduction of descriptor (DAE) systems Model reduction for e.g. parabolic PDE systems Parallel algorithms for sparse computations in model reduction Development/validation of control algorithms based on reduced models In general: Model reduction of structured systems Model reduction of large-scale systems p.38
50 References A.C. Antoulas, Approximation of large-scale dynamical systems, Advances in Design and Control, SIAM, Philadelphia (2005). Dan Sorensen s web site: sorensen Thanos Antoulas web site: aca Model reduction of large-scale systems p.39
Model reduction of large-scale dynamical systems
Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview and some new results Thanos Antoulas Rice University email: aca@rice.edu URL: www.ece.rice.edu/ aca LinSys2008, Sde Boker, 15-19 September 2008 Thanos
More informationGramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csam.ucl.ac.be
More informationGramians of structured systems and an error bound for structure-preserving model reduction
Gramians of structured systems and an error bound for structure-preserving model reduction D.C. Sorensen and A.C. Antoulas e-mail:{sorensen@rice.edu, aca@rice.edu} September 3, 24 Abstract In this paper
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationFluid flow dynamical model approximation and control
Fluid flow dynamical model approximation and control... a case-study on an open cavity flow C. Poussot-Vassal & D. Sipp Journée conjointe GT Contrôle de Décollement & GT MOSAR Frequency response of an
More informationAN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL MOTIVATING EXAMPLE INVERTED PENDULUM
Controls Lab AN OVERVIEW OF MODEL REDUCTION TECHNIQUES APPLIED TO LARGE-SCALE STRUCTURAL DYNAMICS AND CONTROL Eduardo Gildin (UT ICES and Rice Univ.) with Thanos Antoulas (Rice ECE) Danny Sorensen (Rice
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview and some new results Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/ aca Electrical and Computer Engineering
More informationModel reduction of large-scale systems by least squares
Model reduction of large-scale systems by least squares Serkan Gugercin Department of Mathematics, Virginia Tech, Blacksburg, VA, USA gugercin@mathvtedu Athanasios C Antoulas Department of Electrical and
More informationModel Reduction for Linear Dynamical Systems
Summer School on Numerical Linear Algebra for Dynamical and High-Dimensional Problems Trogir, October 10 15, 2011 Model Reduction for Linear Dynamical Systems Peter Benner Max Planck Institute for Dynamics
More informationGramian Based Model Reduction of Large-scale Dynamical Systems
Gramian Based Model Reduction of Large-scale Dynamical Systems Paul Van Dooren Université catholique de Louvain Center for Systems Engineering and Applied Mechanics B-1348 Louvain-la-Neuve, Belgium vdooren@csamuclacbe
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationComputation and Application of Balanced Model Order Reduction
Computation and Application of Balanced Model Order Reduction D.C. Sorensen Collaborators: M. Heinkenschloss, A.C. Antoulas and K. Sun Support: NSF and AFOSR MIT Nov 2007 Projection Methods : Large Scale
More informationPrinciple Component Analysis and Model Reduction for Dynamical Systems
Principle Component Analysis and Model Reduction for Dynamical Systems D.C. Sorensen Virginia Tech 12 Nov 2004 Protein Substate Modeling and Identification PCA Dimension Reduction Using the SVD Tod Romo
More informationAutomatic Order Reduction for Finite Element Models
Automatic Order Reduction for Finite Element Models Andreas Greiner, Jan Lienemann, Evgenii Rudnyi, and Jan G. Korvink IMTEK, Albert-Ludwigs-Universität Freiburg Geroges-Köhler-Allee 103, D-79110 Freiburg,
More informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationOn Dominant Poles and Model Reduction of Second Order Time-Delay Systems
On Dominant Poles and Model Reduction of Second Order Time-Delay Systems Maryam Saadvandi Joint work with: Prof. Karl Meerbergen and Dr. Elias Jarlebring Department of Computer Science, KULeuven ModRed
More informationModel Reduction of Inhomogeneous Initial Conditions
Model Reduction of Inhomogeneous Initial Conditions Caleb Magruder Abstract Our goal is to develop model reduction processes for linear dynamical systems with non-zero initial conditions. Standard model
More informationMoving Frontiers in Model Reduction Using Numerical Linear Algebra
Using Numerical Linear Algebra Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory Group Magdeburg, Germany Technische Universität Chemnitz
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More informationChris Beattie and Serkan Gugercin
Krylov-based model reduction of second-order systems with proportional damping Chris Beattie and Serkan Gugercin Department of Mathematics, Virginia Tech, USA 44th IEEE Conference on Decision and Control
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationS. Gugercin and A.C. Antoulas Department of Electrical and Computer Engineering Rice University
Proceedings of the 39" IEEE Conference on Decision and Control Sydney, Australia December, 2000 A Comparative Study of 7 Algorithms for Model Reduct ion' S. Gugercin and A.C. Antoulas Department of Electrical
More informationBALANCING AS A MOMENT MATCHING PROBLEM
BALANCING AS A MOMENT MATCHING PROBLEM T.C. IONESCU, J.M.A. SCHERPEN, O.V. IFTIME, AND A. ASTOLFI Abstract. In this paper, we treat a time-domain moment matching problem associated to balanced truncation.
More informationModel reduction of large-scale systems
Model reduction of large-scale systems Lecture IV: Model reduction from measurements Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/ aca International School,
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 WeC10.4 An iterative SVD-Krylov based method for model reduction
More informationVirtual Prototyping for Power Electronics
Virtual Prototyping for Power Electronics Cross-Theme Project (Design Tools and Modelling) EPSRC Centre for Power Electronics Dr Xibo Yuan 5 th July 2016 Contents Background Challenges and approaches Accuracy
More informationModel reduction of interconnected systems
Model reduction of interconnected systems A Vandendorpe and P Van Dooren 1 Introduction Large scale linear systems are often composed of subsystems that interconnect to each other Instead of reducing the
More informationPassivity Preserving Model Reduction for Large-Scale Systems. Peter Benner.
Passivity Preserving Model Reduction for Large-Scale Systems Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik Sonderforschungsbereich 393 S N benner@mathematik.tu-chemnitz.de SIMULATION
More informationFourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics
Fourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics K. Willcox and A. Megretski A new method, Fourier model reduction (FMR), for obtaining stable, accurate, low-order models
More informationStability preserving post-processing methods applied in the Loewner framework
Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 1 / 20 Stability preserving post-processing methods applied in the Loewner framework Ion Victor
More informationPassive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form
Passive Interconnect Macromodeling Via Balanced Truncation of Linear Systems in Descriptor Form Boyuan Yan, Sheldon X.-D. Tan,PuLiu and Bruce McGaughy Department of Electrical Engineering, University of
More informationUsing Model Reduction techniques for simulating the heat transfer in electronic systems.
Using Model Reduction techniques for simulating the heat transfer in electronic systems. -Pramod Mathai, Dr. Benjamin Shapiro, University of Maryland, College Park Abstract: There is an increasing need
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
An iterative SVD-Krylov based method for model reduction of large-scale dynamical systems Serkan Gugercin Department of Mathematics, Virginia Tech., Blacksburg, VA, USA, 24061-0123 gugercin@math.vt.edu
More informationParallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation
Parallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation Peter Benner 1, Enrique S. Quintana-Ortí 2, Gregorio Quintana-Ortí 2 1 Fakultät für Mathematik Technische Universität
More informationAn iterative SVD-Krylov based method for model reduction of large-scale dynamical systems
Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 1964 1986 www.elsevier.com/locate/laa An iterative SVD-Krylov based method for model reduction of large-scale dynamical
More informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More informationH 2 -optimal model reduction of MIMO systems
H 2 -optimal model reduction of MIMO systems P. Van Dooren K. A. Gallivan P.-A. Absil Abstract We consider the problem of approximating a p m rational transfer function Hs of high degree by another p m
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More informationProblem set 5 solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION Problem set 5 solutions Problem 5. For each of the stetements below, state
More informationLarge-scale and infinite dimensional dynamical systems approximation
Large-scale and infinite dimensional dynamical systems approximation Igor PONTES DUFF PEREIRA Doctorant 3ème année - ONERA/DCSD Directeur de thèse: Charles POUSSOT-VASSAL Co-encadrant: Cédric SEREN 1 What
More informationKrylov Subspace Methods for Nonlinear Model Reduction
MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationModel Order Reduction using SPICE Simulation Traces. Technical Report
Model Order Reduction using SPICE Simulation Traces Paul Winkler, Henda Aridhi, and Sofiène Tahar Department of Electrical and Computer Engineering, Concordia University, Montreal, Canada pauwink@web.de,
More informationRational Krylov Methods for Model Reduction of Large-scale Dynamical Systems
Rational Krylov Methods for Model Reduction of Large-scale Dynamical Systems Serkan Güǧercin Department of Mathematics, Virginia Tech, USA Aerospace Computational Design Laboratory, Dept. of Aeronautics
More informationArnoldi Methods in SLEPc
Scalable Library for Eigenvalue Problem Computations SLEPc Technical Report STR-4 Available at http://slepc.upv.es Arnoldi Methods in SLEPc V. Hernández J. E. Román A. Tomás V. Vidal Last update: October,
More informationKrylov-based model reduction of second-order systems with proportional damping
Krylov-based model reduction of second-order systems with proportional damping Christopher A Beattie and Serkan Gugercin Abstract In this note, we examine Krylov-based model reduction of second order systems
More informationComputing Transfer Function Dominant Poles of Large Second-Order Systems
Computing Transfer Function Dominant Poles of Large Second-Order Systems Joost Rommes Mathematical Institute Utrecht University rommes@math.uu.nl http://www.math.uu.nl/people/rommes joint work with Nelson
More informationFast Frequency Response Analysis using Model Order Reduction
Fast Frequency Response Analysis using Model Order Reduction Peter Benner EU MORNET Exploratory Workshop Applications of Model Order Reduction Methods in Industrial Research and Development Luxembourg,
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More informationWavelet-Based Passivity Preserving Model Order Reduction for Wideband Interconnect Characterization
Wavelet-Based Passivity Preserving Model Order Reduction for Wideband Interconnect Characterization Mehboob Alam, Arthur Nieuwoudt, and Yehia Massoud Rice Automated Nanoscale Design Group Rice University,
More informationComparison of Model Reduction Methods with Applications to Circuit Simulation
Comparison of Model Reduction Methods with Applications to Circuit Simulation Roxana Ionutiu, Sanda Lefteriu, and Athanasios C. Antoulas Department of Electrical and Computer Engineering, Rice University,
More informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
More informationIntroduction to Model Order Reduction
Introduction to Model Order Reduction Lecture 1: Introduction and overview Henrik Sandberg Kin Cheong Sou Automatic Control Lab, KTH ACCESS Specialized Course Graduate level Ht 2010, period 1 1 Overview
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationA comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control
A comparison of model reduction techniques from structural dynamics, numerical mathematics and systems and control B. Besselink a, A. Lutowska b, U. Tabak c, N. van de Wouw a, H. Nijmeijer a, M.E. Hochstenbach
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationBalancing-Related Model Reduction for Large-Scale Systems
Balancing-Related Model Reduction for Large-Scale Systems Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz D-09107 Chemnitz benner@mathematik.tu-chemnitz.de
More informationModel reduction of nonlinear circuit equations
Model reduction of nonlinear circuit equations Tatjana Stykel Technische Universität Berlin Joint work with T. Reis and A. Steinbrecher BIRS Workshop, Banff, Canada, October 25-29, 2010 T. Stykel. Model
More informationModel Reduction (Approximation) of Large-Scale Systems. Introduction, motivating examples and problem formulation Lecture 1
Lecture outines Motivating examples Norm and linear algebra reminder Introduction Summary Model Reduction (Approximation) of Large-Scale Systems Introduction, motivating examples and problem formulation
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationApplication of Lanczos and Schur vectors in structural dynamics
Shock and Vibration 15 (2008) 459 466 459 IOS Press Application of Lanczos and Schur vectors in structural dynamics M. Radeş Universitatea Politehnica Bucureşti, Splaiul Independenţei 313, Bucureşti, Romania
More informationModel Reduction of Linear Systems, an Interpolation point of View
Université Catholique De Louvain Faculté des sciences appliquées Département d ingénierie mathématique Center for Systems Engineering and Applied Mechanics Model Reduction of Linear Systems, an Interpolation
More informationParametric Model Order Reduction for Linear Control Systems
Parametric Model Order Reduction for Linear Control Systems Peter Benner HRZZ Project Control of Dynamical Systems (ConDys) second project meeting Zagreb, 2 3 November 2017 Outline 1. Introduction 2. PMOR
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationKrylov-based model reduction of second-order systems with proportional damping
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 TuA05.6 Krylov-based model reduction of second-order systems
More informationA global Arnoldi method for the model reduction of second-order structural dynamical systems
A global Arnoldi method for the model reduction of second-order structural dynamical systems Thomas Bonin b, Heike Faßbender a, Andreas Soppa,a, Michael Zaeh,b a TU Braunschweig, Carl-Friedrich-Gauß Fakultät,
More informationMASS, STIFFNESS AND DAMPING IDENTIFICATION OF A TWO-STORY BUILDING MODEL
COMPDYN 2 3 rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 25-28 May 2 MASS,
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationModel reduction of coupled systems
Model reduction of coupled systems Tatjana Stykel Technische Universität Berlin ( joint work with Timo Reis, TU Kaiserslautern ) Model order reduction, coupled problems and optimization Lorentz Center,
More informationModel Order Reduction of Continuous LTI Large Descriptor System Using LRCF-ADI and Square Root Balanced Truncation
, July 1-3, 15, London, U.K. Model Order Reduction of Continuous LI Large Descriptor System Using LRCF-ADI and Square Root Balanced runcation Mehrab Hossain Likhon, Shamsil Arifeen, and Mohammad Sahadet
More informationNew Directions in the Application of Model Order Reduction
New Directions in the Application of Model Order Reduction D.C. Sorensen Collaborators: M. Heinkenschloss, K. Willcox S. Chaturantabut, R. Nong Support: AFOSR and NSF MTNS 08 Blacksburg, VA July 2008 Projection
More informationA Novel Approach for Generating Dynamic Compact Models of Thermal Networks Having Large Numbers of Power Sources
A Novel Approach for Generating Dynamic Compact Models of Thermal Networks Having Large Numbers of Power Sources Lorenzo Codecasa, Dario D Amore, Paolo Maffezzoni and Nadia Spennagallo Dipartimento di
More informationA Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay
A Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay Zheng-Jian Bai Mei-Xiang Chen Jin-Ku Yang April 14, 2012 Abstract A hybrid method was given by Ram, Mottershead,
More informationModel Order Reduction for Parameter-Varying Systems
Model Order Reduction for Parameter-Varying Systems Xingang Cao Promotors: Wil Schilders Siep Weiland Supervisor: Joseph Maubach Eindhoven University of Technology CASA Day November 1, 216 Xingang Cao
More informationEfficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers
Efficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers Jens Saak joint work with Peter Benner (MiIT) Professur Mathematik in Industrie und Technik (MiIT) Fakultät für Mathematik
More informationProgress in Parametric Model Reduction. J. Mohring, Paris,
Progress in Parametric Model Reduction J. Mohring, Paris, 15.12.2009 Overview Need for parametric model reduction Optimal H 2 controller for LMS concrete car benchmark The LBF acoustic aquarium Numerical
More informationDATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT
DATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT JULIUS H. SUMIHAR AND MARTIN VERLAAN,2 Faculty of Electrical Engineering, Mathematics and Computer Science, Department
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationMODEL REDUCTION FOR VEHICLE SYSTEMS MODELLING. Khanh V. Q. Nguyen
MODEL REDUCTION FOR VEHICLE SYSTEMS MODELLING by Khanh V. Q. Nguyen A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Applied Science
More informationComparison of methods for parametric model order reduction of instationary problems
Max Planck Institute Magdeburg Preprints Ulrike Baur Peter Benner Bernard Haasdonk Christian Himpe Immanuel Maier Mario Ohlberger Comparison of methods for parametric model order reduction of instationary
More informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationOptimal Approximation of Dynamical Systems with Rational Krylov Methods. Department of Mathematics, Virginia Tech, USA
Optimal Approximation of Dynamical Systems with Rational Krylov Methods Serkan Güǧercin Department of Mathematics, Virginia Tech, USA Chris Beattie and Virginia Tech. jointly with Thanos Antoulas Rice
More informationPseudospectra and Nonnormal Dynamical Systems
Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures
More informationS N. hochdimensionaler Lyapunov- und Sylvestergleichungen. Peter Benner. Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz
Ansätze zur numerischen Lösung hochdimensionaler Lyapunov- und Sylvestergleichungen Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz S N SIMULATION www.tu-chemnitz.de/~benner
More informationOn the solution of large Sylvester-observer equations
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 200; 8: 6 [Version: 2000/03/22 v.0] On the solution of large Sylvester-observer equations D. Calvetti, B. Lewis 2, and L. Reichel
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationModel Reduction and Controller Design for a Nonlinear Heat Conduction Problem Using Finite Element Method
International Journal of Automation and Computing 9(5), October 2012, 474-479 DOI: 10.1007/s11633-012-0669-6 Model Reduction and Controller Design for a Nonlinear Heat Conduction Problem Using Finite Element
More informationImplicit Volterra Series Interpolation for Model Reduction of Bilinear Systems
Max Planck Institute Magdeburg Preprints Mian Ilyas Ahmad Ulrike Baur Peter Benner Implicit Volterra Series Interpolation for Model Reduction of Bilinear Systems MAX PLANCK INSTITUT FÜR DYNAMIK KOMPLEXER
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More information