Towards parametric model order reduction for nonlinear PDE systems in networks
|
|
- Alfred Hopkins
- 6 years ago
- Views:
Transcription
1 Towards parametric model order reduction for nonlinear PDE systems in networks MoRePas II 2012 Michael Hinze Martin Kunkel Ulrich Matthes Morten Vierling Andreas Steinbrecher Tatjana Stykel Fachbereich Mathematik Universität Hamburg October 4, 2012
2 page 1 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
3 page 2 Motivation: Coupled circuit and semiconductor models Aim Accurate reduced order models for semiconductors in networks Validity over relevant parameter range Accurate physical reduced order model of the coupled system
4 page 3 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
5 Parametric MOR for PDEs in networks page 4 Coupled circuit and semiconductor models [M. Günther 01, C. Tischendorf 03] Kirchhoff s laws (no semiconductors) read +12 V Aj = 0, v = A e R2 R3 C2 A: incidence matrix. C1 Voltage-current relations of components: vin R1 R5 R j C = dq C dt (v C, t), j R = g(v R, t), v L = dφ L (j L, t) dt R4 C3 Modified Nodal Analysis: join all equations to DAE system A C dq C dt ( ) ( ) A C e(t), t + A R g A R e(t), t + A L j L (t) + A V j V (t) = A I i s(t), dφ L dt (j L (t), t) A L e(t) = 0, A V e(t) = v s(t).
6 page 5 Coupled circuit and semiconductor models [M. Günther 01, C. Tischendorf 03] How can semiconductors be introduced? replace semiconductor by a (possibly nonlinear) electrical network, stamp semiconductor network into surrounding network, apply Modified Nodal Analysis. Here: use PDE model for semiconductors DD equations.
7 Parametric MOR for PDEs in networks page 6 Coupled circuit and semiconductor models [M. Günther 01, C. Tischendorf 03] PDE-model (drift-diffusion equations) for semiconductors div (ε ψ) = q(n p C), q tn + div J n = qr(n, p), q tp + div J p = qr(n, p), J n = µ nq( U T n n ψ), J p = µ pq( U T p p ψ), on Ω [0, T ] with Ω R d (d = 1, 2, 3). Dirichlet boundary constraints at Γ O,k : Γ Γ Ο,1 Ο,2 ψ(t, x) = next slide, n(t, x) = ñ(x), p(t, x) = p(x) and Neumann boundary constraints at Γ I : Γ Ι + ψ(t, x) ν(x) = J n ν(x) = J p(t, x) ν(x) = 0 Γ Ο,3 Ω or mixed boundary conditions at MI contacts (MOSFETs).
8 page 7 Couple semiconductor to circuit [M. Günther 01, C. Tischendorf 03]
9 page 8 Couple semiconductor to circuit [M. Günther 01, C. Tischendorf 03] Coupling conditions: Γ Γ Ο,1 Ο,2 j S,k (t) = (J n + J p ε t ψ) ν dσ, Γ O,k Γ Ι + ψ(t, x) = ψ bi (x) + (A S e(t)) k for (t, x) [0, T ] Γ O,k, Γ Ο,3 Ω and add current j S to Kirchhoff s current law: A C dq C dt ( ) ( ) A C e, t + A R g A R e, t + A L j L + A V j V +A S j S = A I i s, dφ L dt (j L, t) A L e = 0, A V e = v s. Add DD-equations + coupling conditions for each semiconductor.
10 page 9 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
11 page 10 Mixed formulation The electric field E = ψ plays dominant role in DD-equations. Mixed formulation [Brezzi et al. 05] Provide additional variable g ψ and equation g ψ = ψ. Scaled DD equations then read: λ div g ψ = n p C, tn + ν n div J n = R(n, p), tp + ν p div J p = R(n, p), g ψ = ψ, J n = n ng ψ, J p = p pg ψ.
12 page 11 Finite Element approximation Finite elements piecewise constant ansatz functions for ψ, n and p. Basis functions: ϕ i, i = 1,..., N, N = T. Raviart-Thomas elements for g ψ, J n and J p. Basis functions: φ j, i = 1,..., M, M = E E N. RT 0 := {y : Ω R d : y T (x) = a T + b T x, a T R d, b T R, [y] E ν E = 0, for all inner edges E}. Galerkin ansatz: ψ h (t, x) = N ψ i (t)ϕ i (x), gψ(t, h x) = i=1 M g ψ,j (t)φ j (x), j=1 and analogously for n, p, J n, and J p.
13 Parametric MOR for PDEs in networks page 12 Full model +12 V C1 R2 R3 C2 A C dq C dt ) ) (A C e(t), t + A R g (A R e(t), t +A L j L (t) + A V j V (t)+a S j S (t) = A I i s(t), vin R1 R4 R5 C3 R dφ L dt (j L (t), t) A L e(t) = 0, A V e(t) = vs(t), j S (t) C 1 J n(t) C 2 J p(t) C 3 ġ ψ (t) = 0, 0 ψ(t) M L ṅ(t) n(t) M L ṗ(t) p(t) + A 0 FEM + F(n h, p h, g g ψ (t) ψ h ) b(a S e(t)) = 0. 0 J n(t) 0 J p(t)
14 page 13 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
15 page 14 Basic test circuit, simulation results e 1(t) j V (t) e 2(t) v s(t) ψ(t, x) n(t, x) p(t, x) R input voltage: v s(t) = 5[V ] sin(2πf t) similar results obtained by MECS [Selva Soto] 0.2 frequency = 1 MHz 0.2 frequency = 1 GHz 0.2 frequency = 5 GHz current j V (t) [ma] current j V (t) [ma] current j V (t) [ma] scaled time scaled time scaled time
16 page 15 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
17 page 16 Snapshot-POD (Proper Orthogonal Decomposition) [L. Sirovich 87] Full simulation yields snapshots (here: y = ψ, n, p,...) { y(ti, ) } i=1,...,m span{ ϕ j } j=1,...,n, with y(t i, x) = N y j (t i )ϕ j (x). j=1 Gather coefficients in matrix Y := ( y(t 1 ),..., y(t m) ) R N m. POD in Hilbert space X as eigenvalue problem: Kv k = σ 2 k v k, with K ij := y(t i, ), y(t j, ) X. Note that K = Y MY with M ij = ϕ i, ϕ j X. Write POD in terms of SVD: Then, the s-dimensional POD basis is { u i := N j=1 ŨΣṼ = L Y, with LL := M. } u j i ϕ j ( ), U := ( u 1,..., u s ) := L Ũ (:,1:s). i=1,...,s
18 page 17 Model Order Reduction Simulate the complete network at one or more reference parameters. Take snapshots of the state of each semiconductor at time points t i. Perform POD component wise on ψ, n, p, g ψ, J n and J p. Use the POD basis functions as (non local) FEM ansatz functions: ψ POD (t, x) = s γ ψ,i (t)uψ(x) i i=1 1D FEM ansatz functions for J n first 5 POD basis functions for J n
19 page 18 Reduced model +12 V C1 R2 R3 C2 A C dq C dt ) ) (A C e(t), t + A R g (A R e(t), t +A L j L (t) + A V j V (t)+a S j S (t) = A I i s(t), vin R1 R4 R5 C3 R dφ L dt (j L (t), t) A L e(t) = 0, A V e(t) = vs(t), j S (t) C 1 U Jn γ Jn (t) C 2 U Jp γ Jp (t) C 3 U gψ γ gψ (t) = 0, 0 γ ψ (t) γ n(t) γ n(t) γ p(t) γ p(t) + A 0 POD + U F(n POD, p POD, g γ gψ (t) ψ POD ) U b(a S e(t)) = 0. 0 γ Jn (t) 0 γ Jp (t)
20 page 19 Computational complexity Computational complexity of reduced model still depends on n FEM : U F(n POD, p POD, gψ POD ) = }{{} U }{{} F ( U n γ }{{} n, U pγ p, U gψ γ gψ ). n POD n FEM n FEM n FEM n POD With matrix-matrix multiplications in Jacobian computation: U }{{} n POD n FEM, block-dense F (...) }{{} n FEM n FEM, sparse }{{} U. n FEM n POD, block-dense
21 page 20 Discrete Empirical Interpolation Md. (DEIM) [S. Chaturantabut, D. Sorensen 09] DEIM Do POD on snapshots {F(n(t i ), p(t i ), g ψ(t i ))}, obtain basis W R n FEM n DEIM (block diagonal matrix). Ansatz is overdetermined. Select n DEIM useful rows: F(U nγ n(t), U pγ p(t), U gψ γ gψ (t)) Wc(t) P F(...) P Wc(t). If P W is regular: F(...) Wc(t) = W (P W ) 1 P F(...) The regularity of P W can be guaranteed, see [CS09]. Again we apply the method component-wise.
22 page 21 Discrete Empirical Interpolation Md. (DEIM) [S. Chaturantabut, D. Sorensen 09] Reduced model with DEIM: U F(U nγ n, U pγ p, U gψ γ gψ ) (U W (P W ) 1 ) }{{} n POD n DEIM, block-dense P }{{ F } (U nγ n, U pγ p, U gψ γ gψ ) }{{} n DEIM n FEM Results for 1D-diode: n FEM FEM n POD ROM n DEIM ROM + DEIM sec sec sec sec sec sec sec sec sec. order nfem 1.18 nfem 1.10 nfem 0.578
23 page 22 Discrete Empirical Interpolation Md. (DEIM) [S. Chaturantabut, D. Sorensen 09]
24 page 23 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
25 page 24 Reduced model depends on position of diode in network Bridge rectifier with 4 diodes: e 2(t) 6 Input: v s (t) Output: e 3 (t) e 1 (t) S 1 S 2 4 v s(t) e 1(t) R e 3(t) potential [V] 2 0 S 4 S time [sec] x 10 9
26 page 25 Reduced model depends on position of diode in network The distance between the spaces U 1 and U 2 which are spanned, e.g., by the POD-functions Uψ 1 of the diode S 1 and Uψ 2 of the diode S 2 respectively, is measured by d(u 1, U 2 ) := max u U 1 u 2 =1 min u v 2 = 2 2 λ, v U 2 v 2 =1 where λ is the smallest eigenvalue of the positive definite matrix SS with S ij = u 1 ψ,i, u 2 ψ,j 2. d(u 1, U 2 ) d(u 1, U 3 ) Table: Distances between reduced models in the rectifier network.
27 page 26 Modes MOR yields a similar but different model for the diodes S 1 and S 2 : J n POD basis functions (modes) 0 L mode 1 0 L mode 2 0 L mode 3 0 L mode 4
28 page 27 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
29 page 28 Problem setting MOR test problem Basic circuit with frequency f of the voltage source v s(t) = 5[V ] sin(2πf t) as model parameter. e1(t) jv (t) e2(t) vs(t) ψ(t, x) n(t, x) p(t, x) R Lack of information Select number of snapshots so that (s) = m i=s+1 σ2 i mi=1 σ 2 i tol.
30 page 29 Reduced model at a fixed frequency First test: Compare reduced and unreduced system at a fixed frequency. relative error simulation time [sec] reduced simulation full simulation (s) (s)
31 page 30 Reduced model over parameter space Construction of reduced model requires snapshots from full simulations at reference parameters. Is the model valid over a large parameter space? reference parameter: P 1 := {f 1 } := {10 10 [Hz]} parameter space P = [10 8, ] 10 2 error 10 1 error reference frequencies parameter (frequency)
32 page 31 Reduced model over parameter space - sampling Goal Find new sampling parameter f k+1 (reference frequency) without simulating the full, unreduced system. Set P k+1 := P k {f k+1 }. We do not consider the PDE discretization error. Rigorous upper bound for the error not available E(f ; P k ) = y h (f ) y POD (f ; P k )?(s) where y h := (ψ h, n h, p h, g h ψ, J h n, J h p ), y POD := (ψ POD, n POD,...). Rigorous RB methods, Greedy algorithm [see e.g. A. Patera, G. Rozza 07]: a-posteriori error estimates required. Linear ODEs [see e.g. B. Haasdonk, M. Ohlberger 09]: build difference between residual and unreduced equation to derive an ODE for the error.
33 page 32 Residual based sampling Define residual R(z POD (f ; P k )): insert z POD (f ; P k ) into unreduced equation, 0 ψ POD (t) M L ṅ POD (t) n POD (t) M L ṗ POD (t) p POD (t) R := + A 0 FEM gψ POD(t) + F(n POD, p POD, g ψ POD ) b(epod (t)). 0 Jn POD (t) 0 Jp POD (t) Residual admits different scales. Scale with block diagonal matrix-valued function D(f ) := diag( d ψ(f )I, d n(f )I, d p(f )I, d gψ (f )I, d Jn (f )I, d Jp (f )I ) and choose d ψ(f ) according to d ψ(f j ) R ψ(y POD (f j ; P k )) = ψh (f j ) ψ POD (f j ; P k ), f ψ h j P k. (f j )
34 page 33 Residual based sampling Algorithm: sampling 1. Select f 1 P, P test P, tol > 0, and set k := 1, P 1 := {f 1 }. 2. Simulate the unreduced model at f 1 and calculate the reduced model with POD basis functions U Calculate weight functions d ( ) (f ) > 0 for all f P k. 4. Calculate the scaled residual D(f )R(z POD (f, P k )) for all f P test. 5. Check termination conditions, e.g. max f Ptest D(f )R(z POD (f, P k )) < tol, no progress in weighted residual. 6. Calculate f k+1 := arg max f Ptest D(f )R(z POD (f, P k )). 7. Simulate the unreduced model at f k+1 and create a new reduced model with POD basis U k+1 using also the already available information at f 1,..., f k. 8. Set P k+1 := P k {f k+1 }, k := k + 1 and goto 3.
35 page 34 Numerical example - sampling step 1 Let f 1 := [Hz], P 1 := {10 10 [Hz]}, P = [10 8, ] sampling step error residual reference frequencies parameter (frequency) f 2 = arg max f Ptest D(f )R(z POD (f, P 1 )) = 10 8 [Hz] P 2 = {10 8 [Hz], [Hz]}
36 page 35 Numerical example - sampling step 2 P 2 = {10 8 [Hz], [Hz]} 10 2 sampling step error residual reference frequencies parameter (frequency) f 3 = arg max f Ptest D(f )R(z POD (f, P 2 )) = [Hz] P 3 = {10 8 [Hz], [Hz], [Hz]}
37 page 36 Numerical example - sampling step 3 P 3 = {10 8 [Hz], [Hz], [Hz]} 10 2 sampling step error residual reference frequencies parameter (frequency) Terminate with no progress in residual.
38 page 37 Outline Motivation PDAE-model Finite Element Method Simulation results Construction of the reduced model Location dependence of reduced model Residual based parameter sampling PABTEC and POD, joint work with A. Steinbrecher & Tatjana Stykel Next steps
39 Parametric MOR for PDEs in networks page 38 Combination of PABTEC (Reis & Stykel 2010) and POD; joint work with [A. Steinbrecher, T. Stykel] +12 V subproject V R3 C2 R3 C2 R2 R2 C1 C1 R5 R R5 R vin R1 vin R1 R4 C3 R4 C3 subproject V subproject V vin vin
40 page 39 Combination of PABTEC and POD; Int. J. Numer. Model. 2012
41 page 40 Substitution of nonlinear components for PABTEC and recoupling A. Steinbrecher, T. Stykel (Int. J. Circuits Theory Appl., 2012): Nonlinear inductor current source Nonlinear capacitor voltage source Nonlinear resistor linear circuit with 2 serial resistors and one voltage source parallel to one of the resistors
42 page 41 Combination of PABTEC and POD; Int. J. Numer. Model input in V / output in A time net full / diode full net red / diode full net full / diode red net red / diode red input u V (t)
43 page 42 Next steps Include QDD models. Include EM effects. Generalize approach to other equation networks containing simple and complex components. Thank you for attending! The work reported in this talk is supported by the German Federal Ministry of Education and Research (BMBF), grants 03HIPAE5 & 03MS613D.
SyreNe System Reduction for Nanoscale IC Design
System Reduction for Nanoscale Max Planck Institute for Dynamics of Complex Technical Systeme Computational Methods in Systems and Control Theory Group Magdeburg Technische Universität Chemnitz Fakultät
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 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 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 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 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 informationProper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints
Proper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints Technische Universität Berlin Martin Gubisch, Stefan Volkwein University of Konstanz March, 3 Martin Gubisch,
More informationREDUCED ORDER METHODS
REDUCED ORDER METHODS Elisa SCHENONE Legato Team Computational Mechanics Elisa SCHENONE RUES Seminar April, 6th 5 team at and of Numerical simulation of biological flows Respiration modelling Blood flow
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 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 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 informationPriority Programme The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models
Priority Programme 1962 The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models Carmen Gräßle, Michael Hinze Non-smooth and Complementarity-based
More informationModel reduction for multiscale problems
Model reduction for multiscale problems Mario Ohlberger Dec. 12-16, 2011 RICAM, Linz , > Outline Motivation: Multi-Scale and Multi-Physics Problems Model Reduction: The Reduced Basis Approach A new Reduced
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 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 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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationModel Order Reduction of Electrical Circuits with Nonlinear Elements
Model Order Reduction of Electrical Circuits with Nonlinear Elements Tatjana Stykel and Technische Universität Berlin July, 21 Model Order Reduction of Electrical Circuits with Nonlinear Elements Contents:
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 informationDEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) DEIM-BASED PGD FOR PARAMETRIC NONLINEAR MODEL ORDER REDUCTION JOSE
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 informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationSome recent developments on ROMs in computational fluid dynamics
Some recent developments on ROMs in computational fluid dynamics F. Ballarin 1, S. Ali 1, E. Delgado 2, D. Torlo 1,3, T. Chacón 2, M. Gómez 2, G. Rozza 1 1 mathlab, Mathematics Area, SISSA, Trieste, Italy
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 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 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 informationComparison of the Reduced-Basis and POD a-posteriori Error Estimators for an Elliptic Linear-Quadratic Optimal Control
SpezialForschungsBereich F 3 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz Comparison of the Reduced-Basis and POD a-posteriori Error Estimators for an Elliptic
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 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 informationInductance, RL Circuits, LC Circuits, RLC Circuits
Inductance, R Circuits, C Circuits, RC Circuits Inductance What happens when we close the switch? The current flows What does the current look like as a function of time? Does it look like this? I t Inductance
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 information2017 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media,
217 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising
More informationPDAE models of integrated circuits and perturbation analysis
PDAE models of integrated circuits and perturbation analysis MARTIN BODESTEDT 1 AND CAREN TISCHENDORF 2 SUMMARY A model for a linear electric circuit containing semiconductors is presented. The modified
More informationReduced order modelling for Crash Simulation
Reduced order modelling for Crash Simulation Fatima Daim ESI Group 28 th january 2016 Fatima Daim (ESI Group) ROM Crash Simulation 28 th january 2016 1 / 50 High Fidelity Model (HFM) and Model Order Reduction
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 informationEileen Kammann Fredi Tröltzsch Stefan Volkwein
Universität Konstanz A method of a-posteriori error estimation with application to proper orthogonal decomposition Eileen Kammann Fredi Tröltzsch Stefan Volkwein Konstanzer Schriften in Mathematik Nr.
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 informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationIntroduction to Compact Dynamical Modeling. II.1 Steady State Simulation. Luca Daniel Massachusetts Institute of Technology. dx dt.
Course Outline NS & NIH Introduction to Compact Dynamical Modeling II. Steady State Simulation uca Daniel Massachusetts Institute o Technology Quic Snea Preview I. ssembling Models rom Physical Problems
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationReduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling Toni Lassila, Andrea Manzoni, Gianluigi Rozza CMCS - MATHICSE - École Polytechnique Fédérale
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 informationStructure preserving Krylov-subspace methods for Lyapunov equations
Structure preserving Krylov-subspace methods for Lyapunov equations Matthias Bollhöfer, André Eppler Institute Computational Mathematics TU Braunschweig MoRePas Workshop, Münster September 17, 2009 System
More informationComputer Science Technical Report CSTR-19/2015 June 14, R. Ştefănescu, A. Sandu
Computer Science Technical Report CSTR-19/2015 June 14, 2015 arxiv:submit/1279278 [math.na] 14 Jun 2015 R. Ştefănescu, A. Sandu Efficient approximation of sparse Jacobians for time-implicit reduced order
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
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 informationECS 40, Fall 2008 Prof. Chang-Hasnain Test #3 Version A
ECS 40, Fall 2008 Prof. ChangHasnain Test #3 Version A 10:10 am 11:00 am, Wednesday December 3, 2008 Total Time Allotted: 50 minutes Total Points: 100 1. This is a closed book exam. However, you are allowed
More informationAn Algorithmic Comparison of the Hyper-Reduction and the Discrete Empirical Interpolation Method for a Nonlinear Thermal Problem
Mathematical and Computational Applications Article An Algorithmic Comparison of the Hyper-Reduction and the Discrete Empirical Interpolation Method for a Nonlinear Thermal Problem Felix Fritzen 1, ID,
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More informationMultiobjective PDE-constrained optimization using the reduced basis method
Multiobjective PDE-constrained optimization using the reduced basis method Laura Iapichino1, Stefan Ulbrich2, Stefan Volkwein1 1 Universita t 2 Konstanz - Fachbereich Mathematik und Statistik Technischen
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationModel Order Reduction for Coupled Problems
Max Planck Institute Magdeburg Preprints Peter Benner Lihong Feng Model Order Reduction for Coupled Problems MAX PLANCK INSTITUT FÜR DYNAMIK KOMPLEXER TECHNISCHER SYSTEME MAGDEBURG MPIMD/15-02 February
More informationarxiv: v1 [math.na] 16 Feb 2016
NONLINEAR MODEL ORDER REDUCTION VIA DYNAMIC MODE DECOMPOSITION ALESSANDRO ALLA, J. NATHAN KUTZ arxiv:62.58v [math.na] 6 Feb 26 Abstract. We propose a new technique for obtaining reduced order models for
More informationM. C. Escher: Waterfall. 18/9/2015 [tsl425 1/29]
M. C. Escher: Waterfall 18/9/2015 [tsl425 1/29] Direct Current Circuit Consider a wire with resistance R = ρl/a connected to a battery. Resistor rule: In the direction of I across a resistor with resistance
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationCIRCUIT ELEMENT: CAPACITOR
CIRCUIT ELEMENT: CAPACITOR PROF. SIRIPONG POTISUK ELEC 308 Types of Circuit Elements Two broad types of circuit elements Ati Active elements -capable of generating electric energy from nonelectric energy
More informationOutline. Week 5: Circuits. Course Notes: 3.5. Goals: Use linear algebra to determine voltage drops and branch currents.
Outline Week 5: Circuits Course Notes: 3.5 Goals: Use linear algebra to determine voltage drops and branch currents. Components in Resistor Networks voltage source current source resistor Components in
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 informationReduced Modeling in Data Assimilation
Reduced Modeling in Data Assimilation Peter Binev Department of Mathematics and Interdisciplinary Mathematics Institute University of South Carolina Challenges in high dimensional analysis and computation
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 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 informationModel Order Reduction Techniques
Model Order Reduction Techniques SVD & POD M. Grepl a & K. Veroy-Grepl b a Institut für Geometrie und Praktische Mathematik b Aachen Institute for Advanced Study in Computational Engineering Science (AICES)
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
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 informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationModel reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature. Masayuki Yano
Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony Patera;
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationLecture: Quadratic optimization
Lecture: Quadratic optimization 1. Positive definite och semidefinite matrices 2. LDL T factorization 3. Quadratic optimization without constraints 4. Quadratic optimization with constraints 5. Least-squares
More informationFinding normalized and modularity cuts by spectral clustering. Ljubjana 2010, October
Finding normalized and modularity cuts by spectral clustering Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu Ljubjana 2010, October Outline Find
More informationProject Components. MC34063 or equivalent. Bread Board. Energy Systems Research Laboratory, FIU
Project Components MC34063 or equivalent Bread Board PSpice Software OrCAD designer Lite version http://www.cadence.com/products/orcad/pages/downloads.aspx#pspice More Details on the Introduction CONVERTER
More informationAn optimal control problem for a parabolic PDE with control constraints
An optimal control problem for a parabolic PDE with control constraints PhD Summer School on Reduced Basis Methods, Ulm Martin Gubisch University of Konstanz October 7 Martin Gubisch (University of Konstanz)
More informationUNITEXT La Matematica per il 3+2
UNITEXT La Matematica per il 3+2 Volume 92 Editor-in-chief A. Quarteroni Series editors L. Ambrosio P. Biscari C. Ciliberto M. Ledoux W.J. Runggaldier More information about this series at http://www.springer.com/series/5418
More informationECE 497 JS Lecture - 11 Modeling Devices for SI
ECE 497 JS Lecture 11 Modeling Devices for SI Spring 2004 Jose E. SchuttAine Electrical & Computer Engineering University of Illinois jose@emlab.uiuc.edu 1 Announcements Thursday Feb 26 th NO CLASS Tuesday
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 informationApplication of POD-DEIM Approach on Dimension Reduction of a Diffusive Predator-Prey System with Allee effect
Application of POD-DEIM Approach on 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 informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationEM Simulations using the PEEC Method - Case Studies in Power Electronics
EM Simulations using the PEEC Method - Case Studies in Power Electronics Andreas Müsing Swiss Federal Institute of Technology (ETH) Zürich Power Electronic Systems www.pes.ee.ethz.ch 1 Outline Motivation:
More informationChapter 3. Steady-State Equivalent Circuit Modeling, Losses, and Efficiency
Chapter 3. Steady-State Equivalent Circuit Modeling, Losses, and Efficiency 3.1. The dc transformer model 3.2. Inclusion of inductor copper loss 3.3. Construction of equivalent circuit model 3.4. How to
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationManuel Matthias Baumann
Delft University of Technology Faculty of Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics Nonlinear Model Order Reduction using POD/DEIM for Optimal Control
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 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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationThe Harmonic Balance Method
For Nonlinear Microwave Circuits Hans-Dieter Lang, Xingqi Zhang Thursday, April 25, 2013 ECE 1254 Modeling of Multiphysics Systems Course Project Presentation University of Toronto Contents Balancing the
More informationParameter Selection Techniques and Surrogate Models
Parameter Selection Techniques and Surrogate Models Model Reduction: Will discuss two forms Parameter space reduction Surrogate models to reduce model complexity Input Representation Local Sensitivity
More informationIntroduction to Scientific Computing
(Lecture 5: Linear system of equations / Matrix Splitting) Bojana Rosić, Thilo Moshagen Institute of Scientific Computing Motivation Let us resolve the problem scheme by using Kirchhoff s laws: the algebraic
More informationSystem Reduction for Nanoscale IC Design (SyreNe)
System Reduction for Nanoscale IC Design (SyreNe) Peter Benner February 26, 2009 1 Introduction Since 1993, the German Federal Ministry of Education and Research (BMBF Bundesministerium füa Bildung und
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report June 2009
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 09-21 June 2009 Model order reduction for nonlinear IC models with POD by A. Verhoeven, M. Striebel, E.J.W.
More informationComputing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices
Computing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices Alper Demir David Long Jaijeet Roychowdhury Bell Laboratories Murray Hill New Jersey USA Abstract The main effort in oscillator
More informationHighly-efficient Reduced Order Modelling Techniques for Shallow Water Problems
Highly-efficient Reduced Order Modelling Techniques for Shallow Water Problems D.A. Bistrian and I.M. Navon Department of Electrical Engineering and Industrial Informatics, University Politehnica of Timisoara,
More informationPOD/DEIM Nonlinear model order reduction of an ADI implicit shallow water equations model
Manuscript Click here to view linked References POD/DEIM Nonlinear model order reduction of an ADI implicit shallow water equations model R. Ştefănescu, I.M. Navon The Florida State University, Department
More informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationINDUCTANCE Self Inductance
NDUTANE 3. Self nductance onsider the circuit shown in the Figure. When the switch is closed the current, and so the magnetic field, through the circuit increases from zero to a specific value. The increasing
More informationThe Behavioral Approach to Systems Theory
The Behavioral Approach to Systems Theory... and DAEs Patrick Kürschner July 21, 2010 1/31 Patrick Kürschner The Behavioral Approach to Systems Theory Outline 1 Introduction 2 Mathematical Models 3 Linear
More informationA-Posteriori Error Estimation for Parameterized Kernel-Based Systems
A-Posteriori Error Estimation for Parameterized Kernel-Based Systems Daniel Wirtz Bernard Haasdonk Institute for Applied Analysis and Numerical Simulation University of Stuttgart 70569 Stuttgart, Pfaffenwaldring
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More information