New Directions in the Application of Model Order Reduction
|
|
- Ashlyn Todd
- 6 years ago
- Views:
Transcription
1 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
2 Projection Methods for MOR Impossible Calculations Made Possible with ROM Experiments with many instances of same reduced model Brief Intro to Gramian Based Model Reduction Proper Orthogonal Decomposition (POD) Balanced Reduction Neural Modeling: Nonlinear MOR: Process/Design Variation: Local Reduction Many Interactions Application of Empirical Interpolation (EIM) Monte-Carlo via ROM
3 LTI Model Reduction by Projection ẋ = Ax + Bu y = Cx Approximate x S V = Range(V), a k-diml. subspace i.e. Put x = Vˆx, and then force W T [V ˆx (AVˆx + Bu)] = 0 ŷ = CVˆx If W T V = I k, then the k dimensional reduced model is ˆx = ˆx + ˆBu ŷ = Ĉˆx where  = W T AV, ˆB = W T B, Ĉ = CV.
4 Moment Matching Krylov Subspace Projection Based on Lanczos, Arnoldi, Rational Krylov methods Padé via Lanczos (PVL) Freund, Feldmann Bai Multipoint Rational Interpolation Grimme Gallivan, Grimme, Van Dooren Recent: Optimal H 2 approximation via interpolation Gugercin, Antoulas, Beattie
5 Gramian Based Model Reduction Proper Orthogonal Decomposition (POD) Principal Component Analysis (PCA) ẋ(t) = f(x(t), u(t)), y = g(x(t), u(t)) The Gramian Eigenvectors of P P = o x(τ)x(τ) T dτ P = VS 2 V T Orthogonal Basis x(t) = VSw(t)
6 PCA or POD Reduced Basis Low Rank Approximation x V kˆx k (t) Galerkin condition Global Basis ˆx k = Vk T f(v kˆx k (t), u(t)) Global Approximation Error? (H 2 bound for LTI) Snapshot Approximation to P x V kˆx k 2 σ k+1 P 1 m m x(t j )x(t j ) T = XX T j=1 Truncate SVD : X = VSU T V k S k U T k
7 SVD Compression n 100 SVD of Clown m k ( m + n) v.s. m x n Storage k 20 Advantage of SVD Compression
8 Image Compression - Feature Detection original rank = 10 rank = 30 rank = 50
9 POD in CFD Extensive Literature Karhunen-Loéve, L. Sirovich Burns, King Kunisch and Volkwein Gunzburger Many, many others Incorporating Observations Balancing Lall, Marsden and Glavaski K. Willcox and J. Peraire
10 POD vs. FEM Both are Galerkin Projection POD - Global Basis fns. vs FEM - Local Basis fns. FEM - Complex Behavior via Mesh Refinement/ Higher Order High Dimension - Sparse Matrices POD - Complex Behavior contained in Global Basis Low Dimension - Dense Matrices Caveat: POD is ad hoc: Must sample rich set of inputs and parameter settings Qx: How to automate parameter/input sampling for POD
11 POD for LTI systems Impulse Response: H(s) = C(sI A) 1 B, s 0 Input to State Map: Controllability Gramian: x(t) = e At B P = o x(τ)x(τ) T dτ = o e Aτ BB T e AT τ dτ State to Output Map: Observability Gramian: y(t) = Ce At x(0) Q = o e AT τ C T Ce Aτ dτ
12 Balanced Reduction (Moore 81) Lyapunov Equations for system Gramians AP + PA T + BB T = 0 A T Q + QA + C T C = 0 With P = Q = S : Want Gramians Diagonal and Equal States Difficult to Reach are also Difficult to Observe Reduced Model A k = W T k AV k, B k = W T k B, C k = C k V k PVk = W k S k QW k = V k S k Reduced Model Gramians Pk = S k and Q k = S k.
13 Hankel Norm Error estimate (Glover 84) Why Balanced Truncation? Hankel singular values = λ(pq) Model reduction H error (Glover) y ŷ 2 2 (sum neglected singular values) u 2 Extends to MIMO Preserves Stability Key Challenge Approximately solve large scale Lyapunov Equations in Low Rank Factored Form
14 CD Player Frequency Response y ŷ 2 2 (sum neglected singular values) u 2 Freq Response CD Player : τ = 0.001, n = 120, k = Original Reduced G(jω) Frequency ω
15 CD Player Frequency Response y ŷ 2 2 (sum neglected singular values) u 2 Freq Response CD Player : τ = 1e 005, n = 120, k = Original Reduced G(jω) Frequency ω
16 CD Player - Hankel Singular Values λ(pq) y ŷ 2 2 (sum neglected singular values) u Hankel Singular Values
17 Approximate Balancing AP + PA T + BB T = 0 A T Q + QA + C T C = 0 Sparse Case: Iteratively Solve in Low Rank Factored Form, P U k U T k, Q L kl T k [X, S, Y] = svd(u T k L k) W k = LY k S 1/2 k and V k = UX k S 1/2 k. Now: PW k V k S k and QV k W k S k
18 Recent Progress LTI MOR Low Rank Approximate Solutions to Lyapunov Eqns n = 1M Now Possible Large Scale BTMOR Possible Optimal H2 reduction IRKA Promising for Large no. Inputs Descriptor Systems e.g. Stykel (LAA 06) general theory and approach
19 Balanced Truncation MOR of Oseen Eqn. Semi-Discrete Oseen Equations: A Descriptor System E d v(t) = Av(t) + Bu(t), dt y(t) = Cv(t) + Du(t) 150 full reduced full reduced Time Frequency Figure: Time response (left) and frequency response (right) for the full order model (circles) and for the reduced order model (solid line).
20 Velocity Profile full 22K dof reduced 15 dof Figure: Velocities generated with the full order model (left column) and with the reduced order model (right column) at t = , , , , (top to bottom). Heinkenschloss, Sun, S., SISSC 08
21 NOTICE State Variables will be y for remainder
22 Reduced Order Neural Modeling Steve Cox Tony Kellems LINEAR MODELS Balanced Truncation Optimal H 2 Nan Xiao and Derrick Roos Ryan Nong Complex Model (Dim 160 K) 20 variable ROM NONLINEAR MODELS Empirical Interpolation (EIM) - Patera Saifon Chaturantabut T. Kellems Nonlinear H-H Neuron Model (Dim 1198) 30 variable ROM Complex nonlinear behavior well approximated
23 Neuron Image AR A Image from Neuromart J.O. Martinez, Rice-Baylor Archive of Neuronal Morphology, cox/neuromart/, accessed 29 July 08 Developed using Neurolucida and NeuroExplorer
24 90µm Neuron Cell
25 Hodgkin-Huxley Neuron Model Full Non-Linear Model I j,syn is the synaptic input into branch j a j 2R i xx v j = C m t v j + G Na m 3 j h j (v j E Na ) + G K n 4 j (v j E K ) + G l (v j E l ) + I j,syn (x, t) Kinetics of the potassium (n) and sodium (h, m) channels t m j = α m (v j )(1 m j ) β m (v j )m j t h j = α h (v j )(1 h j ) β h (v j )h j t n j = α n (v j )(1 n j ) β n (v j )n j.
26 Ultimate Goal for Neural Modeling Ultimate goal is to simulate a few-million neuron system over a minute of brain-time: Feasibility Demonstrated Currently limited to a 10K neuron system over a few brain-seconds without new technology Single Cell Simulation Time Reduction: times ROM computation time: seconds - few minutes neuron types sufficient for Cortex Parallel computing required - under development(kellems)
27 v Cell Response - Lin and Non-Lin Morphology: RC B.asc h n Lin HH BT Full HH m t (ms)
28 v Cell Response - Linear Morphology: RC B.asc h n Lin HH BT Full HH m t (ms)
29 v Cell Response - Near Threshold 2 h n 1.5 Morphology: RC B.asc Lin HH BT Full HH m t (ms)
30 Optimal H 2 Methods: IRKA PROBLEM: Many inputs Controllability Gramian Not Low Rank Kellems and Nong Using Variant of IRKA Gugercin, Antoulas, Beattie (2008) Reduced 160K Neuron Model to a ROM of order 20. Solve times to construct ROM are under 2 mins High End Workstation(Sun Ultra 20) in MATLAB Parameter Study Experiment 158 hrs (full) 3.4 hrs (ROM)
31 ROM Results on Realistic Neuron AR A (Rosenkranz lab) Full system size 6726 Reduced system size 25 B 60µm Soma Potential w.r.t. rest (mv) Comparison of Soma Voltages Nonlinear Linearized Reduced (k = 25) Soma Voltage Error (mv) 10 2 Errors in Soma Voltage (Lin. vs. Reduced) Abs. Error Rel. Error Time (ms) Time (ms)
32 Error vs HS-Values: AR A Error follows Hankel singular value decay: AR A, 35 Stimuli A 0 Max. Abs. Error 25 75% Error Bars Normalized HSV s log 10 of Error (v lin vs. v BT ) Number of Singular Values used (10 runs/value)
33 Highly Branched Neuron n408 (Pyapali et al., 1998) Full system size 41, , 330 Reduced system C size µm A Seconds Time to compute reduced model using IRKA 1000 h = 2 µm (N = 41364) 800 h = 0.5 µ m (N = ) B log 10 of Max. Abs. Error Output error vs. full system k (reduced system size) k (reduced system size)
34 Model Reduction of Nonlinear Terms Saifon Chaturantabut Implemenation of EIM method of Patera et.al. (2004) Test Case: FitzHugh-Nagumo equations εv t (x, t) = ε 2 v xx (x, t) + f (v(x, t)) w(x, t) w t (x, t) = βv(x, t) γw(x, t), f (v) = v(v 0.1)(1 v) IC s and BC s v(x, 0) = 0, w(x, 0) = 0 x [0, L] v x (0, t) = i 0 (t), v x (L, t) = 0 t 0 After FEM discretization, Ey t = Ay + g(t) + N (y), y(0) = 0
35 Problem with Direct POD If we apply the POD basis directly to construct a discretized system, the original system of order N: become a system of order k N: where the nonlinear term : E d dt y(t) = Ay(t) + g(t) + N (y(t)) Ẽ d dt ỹ(t) = Ãỹ(t) + g(t) + Ñ (ỹ(t)), Ñ (ỹ(t)) = }{{} U T N (Uỹ(t)) }{{} k N N 1 Computational Complexity still depends on N!!
36 Nonlinear Approximation via EIM WANT: Ñ (ỹ(t)) }{{} C ˆN (ỹ(t)) }{{} k n m n m 1 k, n m N Complexity k Independent of N
37 EIM Steps Ey t = Ay + g(t) + N (y), y(0) = 0 1) Run trajectory and collect snapshots Y = [y 1, y 2,... y m ] perhaps with several different inputs and parameter values. 2) Truncate SVD of snapshots to get a POD basis for Trajectory 3) Collect nonlinear snapshots s j = f (y j ) = [f (y 1,j ), f (y 2,j ),..., f (y N,j )] T 4) Truncate SVD of nonlinear snapshots to get another POD basis for the nonlinear term 5) Select EIM interpolation points and approximate nonlinear term via collocation in the non-linear POD basis 6) Construct the nonlinear ROM from the reduced linear and nonlinear terms
38 Nonlinear Approximation via EIM Contd. HOW: C ˆN (ỹ(t)) = (U T ) [Ψ(x)] T Q(x)dx } Ω {{ } C ( Q 1 z f z (t) ) } {{ } ˆN (ỹ(t)) where fz (t) = [f (Φ(z 1 )ỹ(t)), f (Φ(z 2 )ỹ(t)),... f (Φ(z nm )ỹ(t))] T {zj } Empirical Interpolation Points and where Ψ(x) = [ψ1 (x), ψ 2 (x),... ψ N (x)] - FEM basis Φ(x) = [φ1 (x), φ 2 (x),... φ k (x)] - POD basis - From Snapshots Ψ(x)y(t l ) Q(x) = [q1 (x), q 2 (x),... q nm (x)] - Nonlinear POD basis - From Snapshots s l = f (Ψ(x)y(t l ))
39 EIM: Numerical Example f (y; µ) = (1 y)cos(3πµ(y + 1))e (1+y)µ, 6 POD basis fns from 50 snapshots µ [1, π] uniform Plot of Approximate Functions (dim = 10) with Exact Functions (in black solid line) 2 µ= µ= µ= µ= s(x;µ) x
40 EIM Reduction of FitzHugh-Nagumo Fiber click figure for movie
41 EIM Reduction of FitzHugh-Nagumo Fiber Nonlinear Reduction N = 1024 k = 40 click figure for movie
42 EIM for Hodgkin-Huxley Equations Kellems After FEM or FD discretization HH-equations yield ODE system a(x)c mi I I I v t m t h t n t = H 2R i succinctly written as Ey t = Ay + g(t) + N (y) Implementation Issues v m h n g 0 (t) ) Choose input stimuli to generate snapshots over full wave 2) NL-snapshots generated for nonlinear term N v (v, m, h, n) N v (v, m, h, n) N m (v, m) N h (v, h) N n (v, n) 3) Terms N m (v, m), N h (v, h), N n (v, n) evaluated only at EIM points 3 7 5
43 EIM Reduction of Hodgkin-Huxley Fiber Three Inputs Forked neuron: voltages at node 10 Full (N = 1198) Reduced (k,nm = 30) Voltage (mv) Time (ms)
44 EIM Reduction of HH : 3 inputs Voltage Profile at Various Times (Full vs ROM)
45 Blade Variation Turbine Reliability/Performance Siemens press picture Small variations among blades large impact on forced response Determines high-cycle fatigue properties
46 Parametric ROM for Geometric Shape Variations T. Bui-Thanh,, K. Willcox, and O. Ghattas 2007 Efficient reduced order models for probabilistic analysis of the effects of blade geometry 2D problem governed by the Euler equations Three orders of magnitude reduction in number of states Accurately reproduce CFD Monte Carlo simulation results at fraction of computational cost. Results impossible without ROM
47 Euler Equations DG Discretization Mathematical Model governed by 2D Euler Equations DG + BC s E(w) dy + R(y, u; w) = 0 dt R nonlinear function of y, u; w u R m m-external forcing inputs via BC s w - random vector associated with geometric variability
48 Parametrically Dependent Linear Model Linearize about steady state y ss (g(w)) corresponding to geometry g(w) y = y ss + ỹ E(w) dỹ dt = A(w)ỹ + B(w)u z = C(ỹ) Linearize E, A, B, C : e.g. A(w) A 0 + A 1 w 1 + A 2 w A m w m One time off-line cost: nominal geometry base
49 Construct POD Basis ROM Adaptive Sampling Method Starting with existing POD basis Φ, construct ROM Then iterate: 1) Select new parameter value w j via w j = argmax w z(w) z r (w) s.t. PDE Constraints 2) Collect additional snapshots y(t; w j ) for w = w j 3) Update POD basis Φ 5) Project A i, B i, C i via e.g. Â i = Φ T A i Φ
50 ROM vs Full Linear Response Outputs: z(1) = C L Lift Coefficient z(2) = C M Moment Coefficient
51 ROM Preservation of Statistics Linearized CFD (left) ROM predictions of work per cycle (WPC) Monte Carlo simulation results for 10,000 blade geometries
52 Prediction and Cost CFD (Full) Reduced Model Size 103, Number nonzeros 2,846,056 84,100 Offline cost hrs Online cost hrs 1.10 hrs Blade 1 WPC mean Blade 1 WPC variance Blade 2 WPC mean Blade 2 WPC variance random parameters Computations on 64 bit PC with 3.2GH Pentium Processor WPC = integral of the blade motion times the lift force over one unsteady cycle
53 Summary Gramian Based Model Reduction: POD, Balanced Reduction Optimal H 2 Reduction via IRKA Neural Modeling - Single Cell ROM Many Interactions Example of important class of problems (including Monte Carlo of Stochastic Systems) Nonlinear MOR via EIM Demonstrated effective reduction and feasibility Process Variation - Blade Geometry Effects Example of important future application of MOR
Computation 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 informationApplications of DEIM in Nonlinear Model Reduction
Applications of DEIM in Nonlinear Model Reduction S. Chaturantabut and D.C. Sorensen Collaborators: M. Heinkenschloss, K. Willcox, H. Antil Support: AFOSR and NSF U. Houston CSE Seminar Houston, April
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 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 informationSolving Large Scale Lyapunov Equations
Solving Large Scale Lyapunov Equations D.C. Sorensen Thanks: Mark Embree John Sabino Kai Sun NLA Seminar (CAAM 651) IWASEP VI, Penn State U., 22 May 2006 The Lyapunov Equation AP + PA T + BB T = 0 where
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 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 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 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 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 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 informationDiscrete Empirical Interpolation for Nonlinear Model Reduction
Discrete Empirical Interpolation for Nonlinear Model Reduction Saifon Chaturantabut and Danny C. Sorensen Department of Computational and Applied Mathematics, MS-34, Rice University, 6 Main Street, Houston,
More informationModel 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 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 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 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 informationRICE UNIVERSITY. Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut
RICE UNIVERSITY Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy
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 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 informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationPOD/DEIM 4D-Var for a Finite-Difference Shallow Water Equations Model
.. POD/DEIM 4D-Var for a Finite-Difference Shallow Water Equations Model North Carolina State University rstefan@ncsu.edu POD/DEIM 4D-Var for a Finite-Difference Shallow Water Equations Model 1/64 Part1
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 informationParametric Reduced-Order Models for Probabilistic Analysis of Unsteady Aerodynamic Applications
Parametric Reduced-Order Models for Probabilistic Analysis of Unsteady Aerodynamic Applications T. Bui-Thanh K. Willcox O. Ghattas Methodology is presented to derive reduced-order models for large-scale
More informationPOD/DEIM 4DVAR Data Assimilation of the Shallow Water Equation Model
nonlinear 4DVAR 4DVAR Data Assimilation of the Shallow Water Equation Model R. Ştefănescu and Ionel M. Department of Scientific Computing Florida State University Tallahassee, Florida May 23, 2013 (Florida
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 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 informationParameterized Partial Differential Equations and the Proper Orthogonal D
Parameterized Partial Differential Equations and the Proper Orthogonal Decomposition Stanford University February 04, 2014 Outline Parameterized PDEs The steady case Dimensionality reduction Proper orthogonal
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 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 informationProper Orthogonal Decomposition (POD)
Intro Results Problem etras Proper Orthogonal Decomposition () Advisor: Dr. Sorensen CAAM699 Department of Computational and Applied Mathematics Rice University September 5, 28 Outline Intro Results Problem
More informationMODEL ORDER REDUCTION
MODEL ORDER REDUCTION Approximate yet accurate surrogates for large-scale simulation KAREN E. WILLCOX Science at Extreme Scales: Where Big Data Meets Large-Scale Computing Tutorials Institute for Pure
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 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 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 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 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 informationNONLINEAR MODEL REDUCTION VIA DISCRETE EMPIRICAL INTERPOLATION
SIAM J. SCI. COMPUT. Vol. 3, No. 5, pp. 737 764 c Society for Industrial and Applied Mathematics NONLINEAR MODEL REDUCTION VIA DISCRETE EMPIRICAL INTERPOLATION SAIFON CHATURANTABUT AND DANNY C. SORENSEN
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Methods for Nonlinear Systems Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 65 Outline 1 Nested Approximations 2 Trajectory PieceWise Linear (TPWL)
More informationModel reduction of large-scale systems
Model reduction of large-scale systems An overview Thanos Antoulas joint work with Dan Sorensen aca@rice.edu URL: www-ece.rice.edu/ aca CITI Lecture, Rice, February 2005 Model reduction of large-scale
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 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 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 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 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 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 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 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 informationModel order reduction & PDE constrained optimization
1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell,
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 informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
More informationPOD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM
POD-DEIM APPROACH ON DIMENSION REDUCTION OF A MULTI-SPECIES HOST-PARASITOID SYSTEM G. Dimitriu Razvan Stefanescu Ionel M. Navon Abstract In this study, we implement the DEIM algorithm (Discrete Empirical
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 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 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 informationInexact Solves in Krylov-based Model Reduction
Inexact Solves in Krylov-based Model Reduction Christopher A. Beattie and Serkan Gugercin Abstract We investigate the use of inexact solves in a Krylov-based model reduction setting and present the resulting
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
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 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 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 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 informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Proper Orthogonal Decomposition (POD) Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 43 Outline 1 Time-continuous Formulation 2 Method of Snapshots
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 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 informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationTechnische Universität Berlin
; Technische Universität Berlin Institut für Mathematik arxiv:1610.03262v2 [cs.sy] 2 Jan 2017 Model Reduction for Large-scale Dynamical Systems with Inhomogeneous Initial Conditions Christopher A. Beattie
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 informationBALANCING-RELATED MODEL REDUCTION FOR LARGE-SCALE SYSTEMS
FOR LARGE-SCALE SYSTEMS Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Recent Advances in Model Order Reduction Workshop TU Eindhoven, November 23,
More informationProper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein
Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition
More informationSnapshot location by error equilibration in proper orthogonal decomposition for linear and semilinear parabolic partial differential equations
Numerical Analysis and Scientific Computing Preprint Seria Snapshot location by error equilibration in proper orthogonal decomposition for linear and semilinear parabolic partial differential equations
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationCME 345: MODEL REDUCTION - Projection-Based Model Order Reduction
CME 345: MODEL REDUCTION - Projection-Based Model Order Reduction Projection-Based Model Order Reduction Charbel Farhat and David Amsallem Stanford University cfarhat@stanford.edu 1 / 38 Outline 1 Solution
More informationRapid Modeling of Nonlinearities in Heat Transfer. Jillian Chodak Free. Doctor of Philosophy In Mechanical Engineering
Rapid Modeling of Nonlinearities in Heat Transfer Jillian Chodak Free Dissertation submitted to the faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
More informationComplexity Reduction for Parametrized Catalytic Reaction Model
Proceedings of the World Congress on Engineering 5 Vol I WCE 5, July - 3, 5, London, UK Complexity Reduction for Parametrized Catalytic Reaction Model Saifon Chaturantabut, Member, IAENG Abstract This
More informationTransient Sensitivity Analysis CASA Day 13th Nov 2007 Zoran Ilievski. Zoran Ilievski Transient Sensitivity Analysis
CASA Day 13th Nov 2007 Talk Structure Talk Structure Introduction Talk Structure Introduction Recap Sensitivity Talk Structure Introduction Recap Sensitivity Examples and Results Talk Structure Introduction
More informationAdaptive rational Krylov subspaces for large-scale dynamical systems. V. Simoncini
Adaptive rational Krylov subspaces for large-scale dynamical systems V. Simoncini Dipartimento di Matematica, Università di Bologna valeria@dm.unibo.it joint work with Vladimir Druskin, Schlumberger Doll
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
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 informationLow Rank Approximation Lecture 3. Daniel Kressner Chair for Numerical Algorithms and HPC Institute of Mathematics, EPFL
Low Rank Approximation Lecture 3 Daniel Kressner Chair for Numerical Algorithms and HPC Institute of Mathematics, EPFL daniel.kressner@epfl.ch 1 Sampling based approximation Aim: Obtain rank-r approximation
More informationReduced Basis Method for Parametric
1/24 Reduced Basis Method for Parametric H 2 -Optimal Control Problems NUMOC2017 Andreas Schmidt, Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation, University of Stuttgart andreas.schmidt@mathematik.uni-stuttgart.de
More informationModel Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University
Model Order Reduction via Matlab Parallel Computing Toolbox E. Fatih Yetkin & Hasan Dağ Istanbul Technical University Computational Science & Engineering Department September 21, 2009 E. Fatih Yetkin (Istanbul
More informationModel Order Reduction
Model Order Reduction Wil Schilders NXP Semiconductors & TU Eindhoven November 26, 2009 Utrecht University Mathematics Staff Colloquium Outline Introduction and motivation Preliminaries Model order reduction
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
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 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 informationBalancing-Related Model Reduction for Parabolic Control Systems
Balancing-Related Model Reduction for Parabolic Control Systems Peter Benner Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany (e-mail: benner@mpi-magdeburg.mpg.de). Abstract:
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 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 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 informationFourier Series for Accurate, Stable, Reduced-Order Models for Linear CFD Applications.
Fourier Series for Accurate, Stable, Reduced-Order Models for Linear Applications. K. Willcox and A. Megretski Massachusetts Institute of Technology Cambridge, MA 39 Abstract A new method, Fourier model
More informationApplication of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect
Application of POD-DEIM Approach for Dimension Reduction of a Diffusive Predator-Prey System with Allee Effect Gabriel Dimitriu 1, Ionel M. Navon 2 and Răzvan Ştefănescu 2 1 The Grigore T. Popa University
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
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 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 informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationOrder Reduction of a Distributed Parameter PEM Fuel Cell Model
Order Reduction of a Distributed Parameter PEM Fuel Cell Model María Sarmiento Carnevali 1 Carles Batlle Arnau 2 Maria Serra Prat 2,3 Immaculada Massana Hugas 2 (1) Intelligent Energy Limited, (2) Universitat
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 informationMethods for eigenvalue problems with applications in model order reduction
Methods for eigenvalue problems with applications in model order reduction Methoden voor eigenwaardeproblemen met toepassingen in model orde reductie (met een samenvatting in het Nederlands) Proefschrift
More information