Model Reduction for Dynamical Systems
|
|
- Jean Beasley
- 5 years ago
- Views:
Transcription
1 MAX PLANCK INSTITUTE Otto-von-Guericke Universitaet Magdeburg Faculty of Mathematics Summer term 2015 Model Reduction for Dynamical Systems Lecture 10 Peter Benner Lihong Feng Max Planck Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory Magdeburg, Germany FOR DYNAMICS OF COMPLEX TECHNICAL SYSTEMS MAGDEBURG Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 1/31
2 Contents 1 Linear parametric systems 2 PMOR based on Multi-moment matching 3 A Robust Algorithm 4 IRKA based PMOR 5 Steady systems 6 Extension to nonlinearities Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 2/31
3 Example 1 A microthruster Upper-left 1 : the structure of an array of pyrotechnical thrusters. Lower-right: the structure of a 2D-axisymmetric model. A model of the microthruster unit. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 3/31
4 Example 1 When the PolySilicon (green) in the middle is excited by a current, the fuel below is ignited and the explosion will occur through the nozzle. The thermal process can be modeled by a heat transfer partial differential equation, while the heat exchange through device interfaces is modeled by convection boundary conditions with different film coefficients h t, h s, h b. The film coefficients h t, h s, h b respectively describe the heat exchange on the top, side, and bottom of the microthruster with the outside surroundings. The values of the film coefficients can change from 1 to 10 9 Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 4/31
5 Example 1 After finite element discretization of the 2D-axisymmetric model, a parameterized system is derived, Eẋ = (A h t A t h s A s h b A b )x + B y = Cx. (1) Here, h t, h s, h b are the parameters and the dimension of the system is n = 4, 257. We observe the temperature at the center of the PolySilicon heater changing with time and the film coefficient, which defines the output of the system 2. 2 Detailed description of the parameterized system can be find at Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 5/31
6 Example 2 The second example is a butterfly gyroscope. The parameterized system is obtained by finite element discretization of the model for the gyroscope (The details of the model can be found in [Moosmann07]). Scheme of the butterfly gyroscope [Moosmann07]. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 6/31
7 Example 2 The paddles of the device are excited to a vibration z(t), where all paddles vibrate in phase. With the external rotation φ, the Coriolis force acts upon the paddles, which causes an out-of-phase movement measured as the z-displacement difference δz between the two red dotted nodes. The interesting output of the system is δz, the difference of the displacement z(t) between the two end nodes depicted as red dots on the same side of the bearing. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 7/31
8 Example 2 The system is of the following form: M(d)ẍ + D(θ, α, β, d)ẋ + T (d)x = Bu(t) y = Cx. (2) M(d) = (M 1 + dm 2 ), D(θ, α, β, d) = θ(d 1 + dd 2 ) + αm(d) + βt (d), and T (d) = (T d T 2 + dt 3 ). Parameters d, θ, α, β. d is the width of the bearing, and θ is the rotation velocity along the x axis. α, β are used to form the Rayleigh damping matrices αm(d), βt (d) in D(θ, α, β, d). The dimension of the system is n = Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 8/31
9 Example 3 The third example is a silicon-nitride membrane 3. This structure resembles a micro-hotplate similar to other micro-fabricated devices such as gas sensors [GrafBT04] and infrared sources [SpannSH05]. Temperature distribution over the silicon-nitride membrane. 3 Picture courtesy of T. Bechtold, IMTEK, University of Freiburg, Germany. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 9/31
10 Example 3 The model of the silicon-nitride membrane is a system with four parameters [BechtoldHRG10]. (E 0 + ρc p E 1 )ẋ + (K 0 + κk 1 + hk 2 )x = Bu(t) y = Cx. (3) The mass density ρ in kg/m 3, the specific heat capacity c p in J/kg/K, the thermal conductivity in W/m/K, and the heat transfer coefficient h in W/m 2 /K. The dimension of the system is n = Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 10/31
11 PMOR based on Multi-moment matching In frequency domain Using Laplace transform, the system in time domain is transformed into E(s 1,..., s p )x = Bu(s p ), y = L T x, (4) where the matrix E R n n is parametrized. The new parameter s p is in fact the frequency parameter s, which corresponds to time t. In case of a nonlinear and/or non-affine dependence of the matrix E on the parameters, the system in (4) is first transformed to an affine form (E 0 + s 1 E 1 + s 2 E s p E p )x = Bu(s p ), y = L T x. (5) Here the newly defined parameters s i, i = 1,..., p, might be some functions (rational, polynomial) of the original parameters s i in (4). Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 11/31
12 PMOR based on multi-moment matching To obtain the projection matrix V for the reduced model, the state x in (5) is expanded into a Taylor series at an expansion point s 0 = ( s 0 1,..., s0 p) T as below, x = [I (σ 1 M σ p M p )] 1 Ẽ 1 Bu(s p ) = [σ 1 M σ p M p ] m Ẽ 1 Bu(s p ) m=0 = m=0 m (k k p) k 2=0... m k p k p 1=0 F m k 2,...,k p (M 1,..., M p ) (6) where σ i = s i s 0 i, Ẽ = E 0 + s 0 1 E s 0 pe p, M i = Ẽ 1 E i, i = 1, 2,... p, and B M = Ẽ 1 B. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 12/31
13 PMOR based on multi-moment-matching σ 0 : L T B M : the 0th order multi-moment; the columns in B M : the 0th order moment vectors. σ 1 : L T M i B M, i = 1, 2,..., p: the first order multi-moments; the columns in M i B M, i = 1, 2,..., p: the first order moment vectors. σ 2 :... ; the columns in M 2 i B M, i = 1, 2,..., p, (M 1 M i + M i M 1 )B M, i = 2,..., p, (M 2 M i + M i M 2 )B M, i = 3,..., p,..., (M p 1 M p + M p M p 1 )B M : the second order moment vectors..... Since the coefficients corresponding not only to s = s p, but also to those associated with the other parameters s i, i = 1,..., p 1 are, we call them as multi-moments of the transfer function. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 13/31
14 PMOR based on multi-moment-matching For the general case, the projection matrix V is constructed as range {V } = colspan{ m q m=0 m (k p+...+k 3) k 2=0... m k p m k p 1=0 k p=0 F m k 2,...,k p (M 1,..., M p )B M } = colspan{b M, M 1 B M, M 2 B M,..., M p B M, (M 1 ) 2 B M, (M 1 M 2 + M 2 M 1 )B M,..., (M 1 M p + M p M 1 )B M, (M 2 ) 2 B M, (M 2 M 3 + M 3 M 2 )B M,...}. (7) Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 14/31
15 A Robust Algorithm Taking a closer look at the power series expansion of x in (6), we get the following equivalent, but different formulation, x = [I (σ 1 M σ p M p )] 1 Ẽ 1 Bu = [σ 1 M σ p M p ] m B M u m=0 = B M u + [σ 1 M σ p M p ]B M u +[σ 1 M σ p M p ] 2 B M u [σ 1 M σ p M p ] j B M u +... (8) By defining x 0 = B M, x 1 = [σ 1 M σ p M p ]B M, x 2 = [σ 1 M σ p M p ] 2 B M,..., x j = [σ 1 M σ p M p ] j B M,..., we have x = (x 0 + x 1 + x x j + )u and obtain the recursive relations Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 15/31
16 A Robust Algorithm x 0 = B M, x 1 = [σ 1 M σ p M p ]x 0, x 2 = [σ 1 M σ p M p ]x 1,... x j = [σ 1 M σ p M p ]x j 1,.... If we define a vector sequence based on the coefficient matrices of x j, j = 0, 1,... as below, R 0 = B M, R 1 = [M 1 R 0, M 2 R 0,..., M p R 0 ], R 2 = [M 1 R 1, M 2 R 1,..., M p R 1 ],. R j = [M 1 R j 1, M 2 R j 1,..., M p R j 1 ], (9). Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 16/31
17 A Robust Algorithm and let R be the subspace spanned by the vectors in R j, j = 0, 1,, m: R = colspan{r 0,..., R j,..., R m }, then there exists z R q, such that x Vz. Here the columns in V R n q is a basis of R. We see that the terms in R j, j = 0, 1,..., m are the coefficients of the parameters in the series expansion (8). They are also the j-th order moment vectors. How to compute an orthonormal basis V? Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 17/31
18 Algorithm 1: Compute V = [v 1, v 2,..., v q ] [Benner, Feng 14] Initialize a 1 = 0, a 2 = 0, sum = 0. Compute R 0 = Ẽ 1 B. if (multiple input) then Orthogonalize the columns in R 0 using MGS: [v 1, v 2,..., v q1 ] = orth{r 0 } with respect to a user given tolerance ε > 0 specifying the deflation criterion for numerically linearly dependent vectors. sum = q 1 % q 1 is the number of columns remaining after deflation w.r.t. ε.) else Compute the first column in V: v 1 = R 0 / R 0 2 sum = 1 end if Compute the orthonormal columns in R 1, R 2,..., R m iteratively as below: Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 18/31
19 continued for i = 1, 2,..., m do a 2 = sum; for t = 1, 2,..., p do if a 1 = a 2 then stop else for j = a 1 + 1,... a 2 do w = Ẽ 1 E t v j ; col = sum + 1; for k = 1, 2,..., col 1 do h = vk T w; w = w hv k end for if w 2 > ε then v col = w w 2 ; sum = col; end if end for end if end for a 1 = a 2 ; end for Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 19/31
20 Adaptively select expansion points Let µ = ( s 1,..., s p ), (µ) is an error estimation, or error bound for ˆx/ŷ, the state/output of the system computed from ROM. Greedy algorithm: Adaptive selection of the expansion points µ i V = []; ɛ = 1; Initial expansion point: µ 0 ; i = 1; Ξ train : a large set of the samples of µ WHILE ɛ > ɛ tol i=i+1; µ i = ˆµ; Use Algorithm 1 to compute V i = span{r 0,..., R q } µ i ; V = [V, V i ]; ˆµ = arg max µ Ξ train (µ); ɛ = (ˆµ); END WHILE. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 20/31
21 Experimental results Example 1: A MEMS model with 4 parameters (benchmark available at M(d)ẍ + D(θ, α, β, d)ẋ + T (d)x = Bu(t), y = Cx. Here, M(d) = (M 1 + dm 2 ), T (d) = (T d T 2 + dt 3 ), D(θ, α, β, d) = θ(d 1 + dd 2 ) + αm(d) + βt (d) R n n, n=17,913. Parameters, d, θ, α, β. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 21/31
22 θ [10 7, 10 5 ], s 2π 1 [0.05, 0.25], d [1, 2]. Ξ train : 3 random θ, 10 random s, 5 random d, α = 0, β = 0 [Salimbahrami et al. 06]. Totally 150 samples of µ (µ i ) ε max true ε tol = (µ i ) ε max true ith iteration step ith iteration step V µ i = span{b M, R 1, R 2} µ i, i = 1,..., 33. ɛ tol = 10 7, ε max true = max H(µ) Ĥ(µ), ROM size=804. µ Ξ train Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 22/31
23 (µ i ) ε max true ε tol = (µ i ) ε max true ith iteration step ith iteration step V µ i = span{b M, R 1} µ i, i = 1,..., 36. ɛ tol = 10 7, ROM size=210. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 23/31
24 Example 2: a silicon nitride membrane (E 0 + ρc p E 1 )dx/dt + (K 0 + κk 1 + hk 2 )x = bu(t) y = Cx. Here, the parameters ρ [3000, 3200], c p [400, 750], κ [2.5, 4], h [10, 12], f [0, 25]Hz Ξ train : 2250 random samples have been taken for the four parameters and the frequency. ε re true = max H(µ) Ĥ(µ) / H(µ), ˆ re (µ) = ˆ (µ)/ Ĥ(µ) µ Ξ train V µ i =span{b M,R 1 }, ɛ re tol = 10 2, n = 60, 020, r = 8, iteration ε re true ˆ re (µ i ) Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 24/31
25 Ξtrain : 3 samples for κ, 10 samples for the frequency. Ξvar : 16 samples for κ, 51 samples for the frequency κ Frequency (Hz) Relative error of the final ROM over Ξvar. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 25/31
26 IRKA based PMOR Consider a linear parametric system C(p 1, p 2,, p l ) dx dt = G(p 1, p 2,, p l )x + B(p 1, p 2,, p l )u(t), y(t) = L(p 1, p 2,, p l ) T x, (10) where the system matrices C(p 1, p 2,, p l ), G(p 1, p 2,, p l ), B(p 1, p 2,, p l ), LT (p 1, p 2,, p l ), are (maybe, nonlinear, non-affine) functions of the parameters p 1, p 2, p l. A straight forward way is [Baur, et.al 11]: Set a group of samples of µ = (p 1,..., p l ): µ 0,..., µ l. For each sample µ i = (p1 i,..., pi l ), i = 1..., l, implement IRKA to get a projection matrix, W i, V i. The final projection matrix: range(v ) = orth(v 1,..., V l ), range(w ) = orth(w 1,..., W l ), W = W (V T W ) 1. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 26/31
27 IRKA based PMOR The reduced parametric model is: Parametric ROM W T C(p 1, p 2,, p l )V dx dt = W T G(p 1, p 2,, p l )Vx +W T B(p 1, p 2,, p l )u(t), y(t) = L(p 1, p 2,, p l ) T Vx, Question: How to select the samples of µ? Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 27/31
28 How to deal with nonaffine matrices? Nonafine matrices are those matrices that cannot be written as: E(p 1,..., p l ) = E 0 + p 1 E , p l E l. PMOR based on multi-moment-matching cannot directly deal with nonaffine case. We must first approximate with affine matrices. IRKA can deal with nonaffine matrices directly. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 28/31
29 Why and How MOR for Steady systems? Steady parametric systems E(p 1,..., p l )x = B(p 1,..., p l ) Solving steady systems for multi-query tasks is also time-consuming. Application of PMOR based on multi-moment-matching to steady systems is straight forward. IRKA cannot. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 29/31
30 Applicable to nonlinear parametric systems? Nonlinear parametric systems: or f (µ, x) = b(µ), E(µ) dx dt = A(µ)x + f (µ, x) = B(µ)u(t), y(t) = L(µ) T x, µ = (p 1,..., p m ), x = x(µ, t). PMOR could be extended to solve weakly nonlinear parametric systems. IRKA can only deal with linear parametric systems. Good candidates for MOR of general nonlinear parametric systems are POD and reduced basis methods. To be introduced: POD and reduced basis method for linear and nonlinear parametric systems. Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 30/31
31 References 1 Daniel, L., Siong, O., Chay, L., Lee, K., White, J.: A multiparameter moment-matching model-reduction approach for generating geometrically parameterized interconnect performance models. IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst. 22(5), (2004). 2 Peter Benner and Lihong Feng: A robust algorithm for parametric model order reduction based on implicit moment-matching. In Reduced Order Methods for Modeling and Computational Reduction. A. Quarteroni, G. Rozza (editors), 9: , Series MS &A, 2014, Springer. 3 U. Baur and C. Beattie and P. Benner and S. Gugercin: Interpolatory Projection Methods for Parameterized Model Reduction. SIAM Journal on Scientific Computing, 33: , And many more... Max Planck Institute Magdeburg Peter Benner Lihong Feng, Model Reduction for Dynamical Systems 31/31
Model Order Reduction for Systems with Non-Rational Transfer Function Arising in Computational Electromagnetics
Model Order Reduction for Systems with Non-Rational Transfer Function Arising in Computational Electromagnetics Lihong Feng and Peter Benner Abstract We consider model order reduction of a system described
More informationParametrische Modellreduktion mit dünnen Gittern
Parametrische Modellreduktion mit dünnen Gittern (Parametric model reduction with sparse grids) Ulrike Baur Peter Benner Mathematik in Industrie und Technik, Fakultät für Mathematik Technische Universität
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 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 informationSome a posteriori error bounds for reduced order modelling of (non-)parametrized linear systems
Max Planck Institute Magdeburg Preprints Lihong Feng Athanasios C. Antoulas Peter Benner Some a posteriori error bounds for reduced order modelling of (non-)parametrized linear systems MAX PLANCK INSTITUT
More informationALBERT-LUDWIGS- UNIVERSITÄT FREIBURG
Compact Electro-Thermal Model for a Microthruster: A Case Study Evgenii B. Rudnyi IMTEK Institute for Microsystem Technology Albert Ludwig University Freiburg Introduction FW5 EU MicroPyros project Developing
More informationA reduced basis method and ROM-based optimization for batch chromatography
MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max
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 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 informationOptimal sensor placement based on model order reduction
P. Benner a, R. Herzog b, N. Lang b, I. Riedel b, J. Saak a a Max Planck Institute for Dynamics of Complex Technical Systems, Computational Methods in Systems and Control Theory, 39106 Magdeburg, Germany
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 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 informationTwo-Sided Arnoldi in Order Reduction of Large Scale MIMO Systems
Two-Sided Arnoldi in Order Reduction of Large Scale MIMO Systems B. Salimbahrami and B. Lohmann Institute of Automation, University of Bremen Otto-Hahn-Allee, NW1 28359 Bremen, Germany e-mail: salimbahrami@iat.uni-bremen.de
More informationParametric Model Order Reduction of a Thermoelectric Generator for Electrically Active Implants
Parametric Model Order Reduction of a Thermoelectric Generator for Electrically Active Implants Cheng Dong Yuan 1,2, Onkar Sandip Jadhav 2, Evgenii B. Rudnyi 3, Dennis Hohlfeld 2, Tamara Bechtold 1,2 1
More informationUsing Model Order Reduction to Accelerate Optimization of Multi-Stage Linear Dynamical Systems
Using Model Order Reduction to Accelerate Optimization of Multi-Stage Linear Dynamical Systems Yao Yue, Suzhou Li, Lihong Feng, Andreas Seidel-Morgenstern, Peter Benner Max Planck Institute for Dynamics
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 informationAutomatic Generation of Compact Electro-Thermal Models for Semiconductor Devices
Automatic Generation of Compact Electro-Thermal Models for Semiconductor Devices Tamara Bechtold, Evgenii B. Rudnyi and Jan G. Korvink IMTEK, University of Freiburg, Georges-Köhler-Allee 103, 79110 Freiburg,
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 informationParameter Preserving Model Reduction for MEMS System-level Simulation and Design
Parameter Preserving Model Reduction for MEMS System-level Simulation and Design E. B. Rudnyi, C. Moosmann, A. Greiner, T. Bechtold, J. G. Korvink IMTEK, Freiburg University, Germany Corresponding Author:
More informationClustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents
Max Planck Institute Magdeburg Preprints Petar Mlinarić Sara Grundel Peter Benner Clustering-Based Model Order Reduction for Multi-Agent Systems with General Linear Time-Invariant Agents MAX PLANCK INSTITUT
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 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 informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationSOME A POSTERIORI ERROR BOUNDS FOR REDUCED-ORDER MODELLING OF (NON-)PARAMETRIZED LINEAR SYSTEMS
ESAIM: MAN 51 (017) 17 158 DOI: 10.1051/man/01701 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-man.org SOME A POSTERIORI ERROR BOUNDS FOR REDUCED-ORDER MODELLING OF (NON-)PARAMETRIZED
More informationThermal Sensors and Actuators
Thermal Sensors and Actuators Part I Fundamentals of heat transfer Heat transfer occurs where there is a temperature gradient until an equilibrium is reached. Four major mechanism Thermal conduction Natural
More informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationOrder reduction numerical methods for the algebraic Riccati equation. V. Simoncini
Order reduction numerical methods for the algebraic Riccati equation V. Simoncini Dipartimento di Matematica Alma Mater Studiorum - Università di Bologna valeria.simoncini@unibo.it 1 The problem Find X
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 informationMethods for Computing Periodic Steady-State Jacob White
Introduction to Simulation - Lecture 15 Methods for Computing Periodic Steady-State Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Periodic Steady-state problems Application
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 information13. Nonlinear least squares
L. Vandenberghe ECE133A (Fall 2018) 13. Nonlinear least squares definition and examples derivatives and optimality condition Gauss Newton method Levenberg Marquardt method 13.1 Nonlinear least squares
More informationTHE AUTOMATIC extraction of parameterized macromodels
2116 IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 26, NO. 12, DECEMBER 2007 A Piecewise-Linear Moment-Matching Approach to Parameterized Model-Order Reduction for
More informationModel Order Reduction of MEMS for Efficient Computer Aided Design and System Simulation
Model Order Reduction of MEMS for Efficient Computer Aided Design and System Simulation Abstract E. B. Rudnyi, J. G. Korvink IMTEK, Freiburg University, Germany, D-79110, korvink@imtek.uni-freiburg.de
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationData-to-Born transform for inversion and imaging with waves
Data-to-Born transform for inversion and imaging with waves Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann Arbor,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
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 informationReduced Order Modeling Enables System Level Simulation of a MEMS Piezoelectric Energy Harvester with a Self-Supplied SSHI-Scheme
Reduced Order Modeling Enables System Level Simulation of a MEMS Piezoelectric Energy Harvester with a Self-Supplied SSHI-Scheme F. Sayed 1, D. Hohlfeld², T. Bechtold 1 1 Institute for Microsystems Engineering,
More informationEE C245 / ME C218 INTRODUCTION TO MEMS DESIGN FALL 2009 PROBLEM SET #7. Due (at 7 p.m.): Thursday, Dec. 10, 2009, in the EE C245 HW box in 240 Cory.
Issued: Thursday, Nov. 24, 2009 PROBLEM SET #7 Due (at 7 p.m.): Thursday, Dec. 10, 2009, in the EE C245 HW box in 240 Cory. 1. Gyroscopes are inertial sensors that measure rotation rate, which is an extremely
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3. 8. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Nonlinear Controllability and Observability Hanz Richter Mechanical Engineering Department Cleveland State University Definition of
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 informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationSTRUCTURED SPATIAL DISCRETIZATION OF DYNAMICAL SYSTEMS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds. Crete Island, Greece, 5 10 June
More informationAccepted to appear in SIAM Review, 2015
Accepted to appear in SIAM Review, 2015 A SURVEY OF PROJECTION-BASED MODEL REDUCTION METHODS FOR PARAMETRIC DYNAMICAL SYSTEMS PETER BENNER, SERKAN GUGERCIN, AND KAREN WILLCOX Abstract. Numerical simulation
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 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 informationInterconnection of Subsystem Reduced-Order Models in the Electro-Thermal Analysis of Large Systems
IEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY 1 Interconnection of Subsystem Reduced-Order Models in the Electro-Thermal Analysis of Large Systems Pramod Mathai, Graduate Research
More informationSGN Advanced Signal Processing Project bonus: Sparse model estimation
SGN 21006 Advanced Signal Processing Project bonus: Sparse model estimation Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 12 Sparse models Initial problem: solve
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationTowards high performance IRKA on hybrid CPU-GPU systems
Towards high performance IRKA on hybrid CPU-GPU systems Jens Saak in collaboration with Georg Pauer (OVGU/MPI Magdeburg) Kapil Ahuja, Ruchir Garg (IIT Indore) Hartwig Anzt, Jack Dongarra (ICL Uni Tennessee
More informationMEMS Gyroscope Control Systems for Direct Angle Measurements
MEMS Gyroscope Control Systems for Direct Angle Measurements Chien-Yu Chi Mechanical Engineering National Chiao Tung University Hsin-Chu, Taiwan (R.O.C.) 3 Email: chienyu.me93g@nctu.edu.tw Tsung-Lin Chen
More informationIncremental identification of transport phenomena in wavy films
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Incremental identification of transport phenomena in
More informationSystem Level Modeling of Microsystems using Order Reduction Methods
System Level Modeling of Microsystems using Order Reduction Methods Sven Reitz, Jens Bastian, Joachim Haase, Peter Schneider, Peter Schwarz Fraunhofer Institute for Integrated Circuits, Branch Lab Design
More informationModel order reduction of large-scale dynamical systems with Jacobi-Davidson style eigensolvers
MAX PLANCK INSTITUTE International Conference on Communications, Computing and Control Applications March 3-5, 2011, Hammamet, Tunisia. Model order reduction of large-scale dynamical systems with Jacobi-Davidson
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationNonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
More informationA Black-Box Method for Parametric Model Order Reduction based on Matrix Interpolation with Application to Simulation and Control
A Black-Box Method for Parametric Model Order Reduction based on Matrix Interpolation with Application to Simulation and Control Matthias Martin Georg Geuß Vollständiger Abdruck der von der Fakultät für
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 informationDefinition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition
6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method The minimization problem We are given a symmetric positive definite matrix R n n and a right hand side vector b R n We want to solve the linear system Find u R n such that
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 informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationOrder reduction of large scale second-order systems using Krylov subspace methods
Linear Algebra and its Applications 415 (2006) 385 405 www.elsevier.com/locate/laa Order reduction of large scale second-order systems using Krylov subspace methods Behnam Salimbahrami, Boris Lohmann Lehrstuhl
More informationModel Reduction for Dynamical Systems
Otto-von-Guericke Universität Magdeburg Faculty of Mathematics Summer term 2015 Model Reduction for Dynamical Systems Lecture 8 Peter Benner Lihong Feng Max Planck Institute for Dynamics of Complex Technical
More informationEmpirical Interpolation of Nonlinear Parametrized Evolution Operators
MÜNSTER Empirical Interpolation of Nonlinear Parametrized Evolution Operators Martin Drohmann, Bernard Haasdonk, Mario Ohlberger 03/12/2010 MÜNSTER 2/20 > Motivation: Reduced Basis Method RB Scenario:
More informationarxiv: v3 [math.na] 6 Jul 2018
A Connection Between Time Domain Model Order Reduction and Moment Matching for LTI Systems arxiv:181.785v3 [math.na] 6 Jul 218 Manuela Hund a and Jens Saak a a Max Planck Institute for Dynamics of Complex
More informationA hybrid control framework for impulsive control of satellite rendezvous
A hybrid control framework for impulsive control of satellite rendezvous Luca Zaccarian Joint work with Mirko Brentari, Sofia Urbina, Denis Arzelier, Christophe Louembet LAAS-CNRS and University of Trento
More informationHessian-based model reduction for large-scale systems with initial condition inputs
Hessian-based model reduction for large-scale systems with initial condition inputs O. Bashir 1, K. Willcox 1,, O. Ghattas 2, B. van Bloemen Waanders 3, J. Hill 3 1 Massachusetts Institute of Technology,
More informationA Trajectory Piecewise-Linear Approach to Model Order Reduction and Fast Simulation of Nonlinear Circuits and Micromachined Devices
A Trajectory Piecewise-Linear Approach to Model Order Reduction and Fast Simulation of Nonlinear Circuits and Micromachined Devices Michał Rewieński, Jacob White Department of Electrical Engineering and
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationLecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 54 Lecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable Peter J. Hammond latest revision 2017
More informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
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 informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationMulti-Robotic Systems
CHAPTER 9 Multi-Robotic Systems The topic of multi-robotic systems is quite popular now. It is believed that such systems can have the following benefits: Improved performance ( winning by numbers ) Distributed
More informationOne dimensional steady state diffusion, with and without source. Effective transfer coefficients
One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport
More informationKernel. Prop-Defn. Let T : V W be a linear map.
Kernel Aim lecture: We examine the kernel which measures the failure of uniqueness of solns to linear eqns. The concept of linear independence naturally arises. This in turn gives the concept of a basis
More informationDownloaded 06/21/12 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No. 5, pp. 2489 2518 c 2011 Society for Industrial and Applied Mathematics INTERPOLATORY PROJECTION METHODS FOR PARAMETERIZED MODEL REDUCTION ULRIKE BAUR, CHRISTOPHER BEATTIE,
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 informationChapter 6. Orthogonality
6.4 The Projection Matrix 1 Chapter 6. Orthogonality 6.4 The Projection Matrix Note. In Section 6.1 (Projections), we projected a vector b R n onto a subspace W of R n. We did so by finding a basis for
More informationBayesian Inference for DSGE Models. Lawrence J. Christiano
Bayesian Inference for DSGE Models Lawrence J. Christiano Outline State space-observer form. convenient for model estimation and many other things. Bayesian inference Bayes rule. Monte Carlo integation.
More informationCS491/691: Introduction to Aerial Robotics
CS491/691: Introduction to Aerial Robotics Topic: Midterm Preparation Dr. Kostas Alexis (CSE) Areas of Focus Coordinate system transformations (CST) MAV Dynamics (MAVD) Navigation Sensors (NS) State Estimation
More information1 Introduction. Minho Lee 1 Jihoon Lee 1 Gunhee Jang 1
DOI 10.1007/s005-015-5-5 TECHNICAL PAPER Stability analysis of a whirling rigid rotor supported by stationary grooved FDBs considering the five degrees of freedom of a general rotor bearing system Minho
More informationModel order reduction of mechanical systems subjected to moving loads by the approximation of the input
Model order reduction of mechanical systems subjected to moving loads by the approximation of the input Alexander Vasilyev, Tatjana Stykel Institut für Mathematik, Universität Augsburg ModRed 2013 MPI
More informationHeat Transfer Analysis
Heat Transfer 2011 Alex Grishin MAE 323 Chapter 8: Grishin 1 In engineering applications, heat is generally transferred from one location to another and between bodies. This transfer is driven by differences
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationCS Tutorial 5 - Differential Geometry I - Surfaces
236861 Numerical Geometry of Images Tutorial 5 Differential Geometry II Surfaces c 2009 Parameterized surfaces A parameterized surface X : U R 2 R 3 a differentiable map 1 X from an open set U R 2 to R
More informationAn Optimum Fitting Algorithm for Generation of Reduced-Order Models
An Optimum Fitting Algorithm for Generation of Reduced-Order Models M.M. Gourary 1, S.G. Rusakov 1, S.L. Ulyanov 1, M.M. Zharov 1, B.J. Mulvaney 2 1 IPPM, Russian Academy of Sciences, Moscow 1523, e-mail:
More informationCOMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT OF A NONLINEAR BISTABLE ELEMENT
11 th International Conference on Vibration Problems Z. Dimitrovová et.al. (eds.) Lisbon, Portugal, 9 12 September 2013 COMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationMODELLING OF RECIPROCAL TRANSDUCER SYSTEM ACCOUNTING FOR NONLINEAR CONSTITUTIVE RELATIONS
MODELLING OF RECIPROCAL TRANSDUCER SYSTEM ACCOUNTING FOR NONLINEAR CONSTITUTIVE RELATIONS L. X. Wang 1 M. Willatzen 1 R. V. N. Melnik 1,2 Abstract The dynamics of reciprocal transducer systems is modelled
More informationParallel Plate Heat Exchanger
Parallel Plate Heat Exchanger Parallel Plate Heat Exchangers are use in a number of thermal processing applications. The characteristics are that the fluids flow in the narrow gap, between two parallel
More informationSeismic imaging and multiple removal via model order reduction
Seismic imaging and multiple removal via model order reduction Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann
More informationMathematical Optimisation, Chpt 2: Linear Equations and inequalities
Mathematical Optimisation, Chpt 2: Linear Equations and inequalities Peter J.C. Dickinson p.j.c.dickinson@utwente.nl http://dickinson.website version: 12/02/18 Monday 5th February 2018 Peter J.C. Dickinson
More informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More information