Proper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
|
|
- Bertha Sutton
- 5 years ago
- Views:
Transcription
1 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, K. Kunisch, O. Lass, F. Leibfritz, H. Müller, T. Tonn, F. Tröltzsch, K. Urban AICES, 3 November 011
2 Motivation 1: Parameter identification Model equations (linear, for simplicity): div ( c u) + β u + au = f in Ω R d c u n + qu = g N on Γ N Γ = Ω ( ) u = g D on Γ D = Γ \ Γ N Problem: estimate parameters (e.g., c, β, a or q) in ( ) from given (perturbed) measurements u d for the solution u on (parts of) Γ Mathematical formulation: ( -dimensional) optimization problem min α u u d ds + κ µ s.t. (u,µ) solves ( ) and µ M ad Γ s.t. subject to Numerical strategy: combine optimization methods with fast (local) rate of convergence and POD model reduction for the PDEs
3 Motivation : Optimal control of time-dependent problems Model problem: min 1 y(t) y T dx + κ T u dxdt Ω 0 Γ y t y + f (y) = 0 in Q = (0,T) Ω s.t. y Γ = u on Σ = (0,T) Γ y(0) = y on Ω R d Adjoint system (for gradient computation): p t p + f (y) p = 0, p Γ = 0, p(t) = y T y(t) Optimizer: second-order methods like SQP or (semismooth) Newton Challenge: large-scale fast/real-time optimization
4 Motivation 3: Closed-loop control for time-dependent PDEs Open-loop control: input u(t) ẋ(t)=f (t,x(t),u(t)) x(0)=x R l (after spatial discretization) output y(t) = Cx(t) + Du(t) Closed-loop control: determine F with u(t) = F(t, y(t)) (feedback law) Linear case: LQR and LQG design Nonlinear case: Hamilton-Jacobi-Bellman eq. (v value function) v t (t,y ) + H(v y (t,y ),y ) = 0 in (0,T) R l Strategy: l-dim. spatial approximation by, e.g., POD basis
5 Model reduction: POD basis containing characteristic information 0,5% der Matrixbasis > 45% Information
6 Model reduction: POD basis containing characteristic information 0,5% der Matrixbasis > 45% Information 1% der Matrixbasis > 56% Information
7 Model reduction: POD basis containing characteristic information 0,5% der Matrixbasis > 45% Information 1% der Matrixbasis > 56% Information 5% der Matrixbasis > 76% Information
8 Model reduction: POD basis containing characteristic information 0,5% der Matrixbasis > 45% Information 1% der Matrixbasis > 56% Information 5% der Matrixbasis > 76% Information 10% der Matrixbasis > 85% Information
9 Model reduction: POD basis containing characteristic information 0,5% der Matrixbasis > 45% Information 1% der Matrixbasis > 56% Information 5% der Matrixbasis > 76% Information 10% der Matrixbasis > 85% Information 0% der Matrixbasis > 9% Information Originalbild
10 General outline of the lecture 1.) The Proper Orthogonal Decomposition (POD) method: What is a POD basis? How can we compute the basis numerically? Which theory is behind?.) Reduced-order modelling (ROM) with POD: What is a POD reduced-order model? How can we derive a-priori error estimates? Can we determine an improved ROM in an adaptive way? 3.) POD suboptimal control: How do we apply the POD in PDE constrained optimization? Can we control the error? What can be done for feedback control?
11 Part 1: The POD method What is a POD basis? How can we compute the basis numerically? Which theory is behind?
12 Outline of the first part: The POD method POD and singular value decomposition (SVD) POD method for ordinary differential equations (ODEs) Continuous POD method for ODEs POD method for partial differential equations (PDEs) References
13 POD method & SVD [Kunisch/V. 99, V. 01, V. 09] Given: y 1,...,y n R m ; set V = span {y 1,...,y n } R m Goal: Find l dim V orthonormal vectors {ψ i } l i=1 in Rm minimizing J(ψ 1,...,ψ l ) = n y j j=1 with the Euclidean norm y = y T y Constrained optimization: l i=1 min J(ψ 1,...,ψ l ) subject to ψ T i ψ j = ( y T j ψ i ) ψi { 1 if i = j 0 otherwise Equivalent problem: Find orthonormal ψ 1,...,ψ l R m maximizing J (ψ 1,...,ψ l ) = n j=1 i=1 l y T j ψ i since y j = m i=1 ( y T j ψ i ) ψi
14 Necessary optimality conditions (Part 1) Lagrange functional: l ( ) L(ψ 1,...,ψ l,λ 11,...,λ ll ) = J(ψ 1,...,ψ l ) + λ ij ψ T i ψ j δ ij i,j=1 with the Kronecker symbol δ ij = 1 for i = j and δ ij = 0 otherwise Optimality conditions: L (ψ 1,...,ψ l,λ 11,...,λ ll ) = 0 R m ψ i L (ψ 1,...,ψ l,λ 11,...,λ ll ) = 0 R λ ij for i = 1,...,l for i,j = 1,...,l
15 Necessary optimality conditions (Part ) L(ψ 1,...,ψ l,λ 11,...,λ ll ) = J(ψ 1,...,ψ l ) + l L ψ i = 0 n j=1 L λ ij = 0 ψ T i ψ j = δ ij i,j=1 λ ij ( ψ T i ψ j δ ij ) y j (y T j ψ i ) = λ ii ψ i and λ ij = 0 for i j Setting λ i = λ ii and Y = [y 1,...,y n ] R m n we have YY T ψ i = λ i ψ i for i = 1,...,l i.e., necessary optimality conditions are given by a symmetric m m eigenvalue problem Here: necessary optimality conditions are already sufficient
16 Computation of the POD basis (Part 1) Optimality conditions: YY T ψ i = λ i ψ i for i = 1,...,l Solution by SVD for Y R m n : d = rank Y, σ 1... σ d > 0, U = [u 1,...,u m ] R m m und V = [v 1,...,v n ] R n n orthogonal with ( ) U T D 0 YV = = Σ R m n 0 0 where D = diag (σ 1,...,σ d ) R d d. Moreover, for 1 i d Yv i = σ i u i, Y T u i = σ i v i, YY T u i = σ i u i, Y T Yv i = σ i v i POD basis: ψ i = u i and λ i = σ i > 0 for i = 1,...,l d = dim V with V = span {y 1,...,y n }
17 Computation of the POD basis (Part ) Data ensemble: V = span {y 1,...,y n } R m and d = dim V POD basis of rank l: ψ i = u i and λ i = σ i > 0 for i = 1,...,l d Three choices to compute the ψ i s SVD for Y R m n : Yv i = σ i u i EVD for YY T R m m : YY T u i = σi u i (if m n) EVD for Y T Y R n n : Y T Yv i = σi v i and u i = 1 σ i Yv i (if m n) Essential error formula for the POD basis of rank l: J(ψ 1,...,ψ l ) = n y j j=1 l i=1 ( y T j ψ i ) ψi = d i=l+1 λ i
18 Computation of the POD basis (Part 3) Essential error formula for the POD basis of rank l: n l ( ) J(ψ 1,...,ψ l ) = y j y T j ψ i ψi = j=1 i=1 YY T ψ i = λ i ψ i, 1 i l, and YY T ψ i = n j=1 j=1 d i=l+1 λ i ( y T j ψ i ) yj give λ i = λ i ψi T ψ i = ( ( YY T ) n ) Tψi ( ) T ψ i = y T j ψ i yj ψ i = y j = d i=1 ( y T j ψ i ) ψi, j = 1,...,m, and ψ T i ψ j = δ ij imply n y j j=1 l i=1 ( y T j ψ i ) ψi = n d j=1 i=l+1 y T j ψ i = n y T j ψ i j=1 d i=l+1 λ i
19 Outline of the first part: The POD method POD and SVD POD method for ODEs Continuous POD method for ODEs POD method for PDEs References
20 POD method for ODEs Nonlinear dynamical system in R m : ẏ(t) = f (t,y(t)) for t (0,T) and y(0) = y with given y 0 R m and f : [0,T] R m R m Time grid: 0 t 1 < t <...t n T, δt j = t j t j 1 for j n Available or known snapshots: y j = y(t j ), 1 j n Snapshot ensemble: V = span {y 1,...,y n }, d = dim V n POD basis of rank l < d: with weights α j 0 min n yj α j j=1 l i=1 ( y T j ψ i ) ψi s.t. ψ T i ψ j = δ ij
21 Computation of the POD basis EVD for linear and symmetric R n in ODE space R m : n R n ( ) u i = α j y j y T j u i = σ i u i (YY T u i = σi u i) j=1 and set λ i = σ i, ψ i = u i EVD for linear and symmetric K n = (( αi α j y T j y i )) in R n : and set λ i = σ i, ψ i = 1 λi K n v i = σ i v i (Y T Yv i = σ i v i) n αj (v i ) j y j j=1 Error formula for the POD basis of rank l: n yj l ( ) α j y T j ψ i ψi = j=1 i=1 d i=l+1 λ i
22 Outline of the first part: The POD method POD and SVD POD method for ODEs Continuous POD method for ODEs POD method for PDEs References
23 Continuous POD method [Kunisch/V. 0, Henri 04, V. 09, Chapelle et al. 11] Snapshots: y(t) for all t [0,T] Snapshot ensemble: V = {y(t) t [0,T]}, d = dim V POD basis of rank l < d: T l ( min y(t) y(t) T ) ψ i ψi dt s.t. ψ T i ψ j = δ ij 0 i=1 Optimality conditions: EVP for linear, symmetric, compact R Rψ i = T 0 ( ψ T i y(t) ) y(t)dt = λ i ψ i for i N Error for the POD basis of rank l: T l ( y(t) ψ T i y(t) ) d ψ i dt = 0 i=1 i=l+1 λ i
24 Relationship between discrete and continuous POD Operators R n and R: R n ψ = Rψ = n ( α j ψ T y(t j ) ) y(t j ) j=1 T 0 ( ψ T y(t) ) y(t)dt Operator convergence of R n R: y smooth and appropriate α j s Perturbation theory [Kato 80]: (λ n i,ψn i ) n (λ i,ψ i ) for 1 i l Choice of the weights α j?: ensure convergence R n n R
25 Outline of the first part: The POD method POD and SVD POD method for ODEs Continuous POD method for ODEs POD method for PDEs References
26 POD method for PDEs Heat equation (for instance): y t y = f y in Q = (0,T) Ω n = g on Σ = (0,T) Γ y(0) = y in Ω R l Variational formulation: for V = H 1 (Ω), f.a.a. t [0,T] y t (t)ϕ + y(t) ϕdx = f (t)ϕdx + g(t)ϕds Ω Ω Γ ϕ V FE Galerkin: y m (t) V m = span {ϕ 1,...,ϕ m }, f.a.a. t [0,T] yt m (t)ϕ + y m (t) ϕdx = f (t)ϕdx + g(t)ϕds ϕ V m Ω Ω Γ
27 POD basis computation Time grid: 0 t 1 < t <...t n T, δt = t j t j 1 for j n FE snapshots: y j = y m (t j ) V m, 1 j n Inner product: u,v = Ω uv dx or u,v = uv + u v dx Ω Sizes: # FE s # time instances, i.e., m n Computation of the correlation K n : α j = O(δt) (1 < j < n) αi α j y m j,y m i = α i α j n Y ik Y jl ϕ l,ϕ k = ( DY T MYD ) ij =: Kn ij k,l=1 with M ij = ϕ j,ϕ i (mass or stiffness matrix) and D = diag ( α i ) POD basis: K n v i = λ i v i and ψ i = 1 λi YDv i
28 Motivation Outline POD and SVD POD for ODEs Continuous POD POD for PDEs References Numerical example: Burgers equation Eigenvalues of the correlation matrix K= u(t ),u(t ) yt νyxx + yyx = f y (, 0 = y (, 1) = 0 y (0, ) = y j L (0,π) i 0 10 in Q = (0, T ) Ω λ % on (0, T ) λ1+λ % in Ω = (0, π) R 10 λ1+λ+λ % y (x) = sin(x) and ν = finite elements 6 10 Time integration with Matlab s ode15s Snapshots V = span {y (t1 ),..., y (t100 )} i axis ψ1 ψ Solution of the Burgers equation x axis x axis ψ ψ 3 6 t axis x axis x axis x axis 4 5 6
29 Numerical example: Energy transport (Boussinesq) u t + uu x + vu y + p x = ν u v t + uv x + vv y + p y = ν v + βθ u x + v y = 0 θ t + uθ x + vθ y = α θ in Q in Q in Q in Q α = 10 5, β = 10, ν = finite elements (Femlab) Time integration with Matlab s ode15s Snapshots at t 1,..., t 1 for u, v and θ 10 0 Eigenvalues of K= u(t i ),u(t j ) L (Ω) und K= v(t i ),v(t j ) L (Ω) u: 97.83% v: 95.98% Eigenvalues of the correlation matrix K= T(t i ),T(t j ) L (Ω) λ % 10 λ 1 +λ 99.84% λ 1 +λ +λ % i axis i axis
30 POD for λ-ω systems [Müller/V. 06] PDEs: s = u + v, λ(s) = 1 s, ω(s) = βs ( ) ( )( ut λ(s) ω(s) u = ω(s) λ(s) v v t ) + ( σ u σ v ) Homogeneous boundary conditions: u = v = 0 or u n = v n = 0 Initial conditions: u (x 1,x ) = x 0.5, v (x 1,x ) = (x 1 0.5)/ u for β=1 and t=100 u for β= and t=100 u for β=3 and t=100 u for β=1 and t=100
31 POD basis for λ-ω systems Offsets: ū(x) = 1 n n u(t j,x) or ū 0 j=1 Snapshots: û j (x) = u(t j,x) ū(x) for 1 j n POD eigenvalue problem: u,v = uv dx Ω K n v i = λv i, 1 i l, with Kij n = α i α j POD basis computation: ψ i = 1 λi 10 0 Decay of the first eigenvalues n αj (v i ) j û j j=1 Ω û j (x)û i (x)dx 10 0 Decay of the first eigenvalues β=1 β=1.5 β= β= β=1 β=1.5 β= β=
32 Elliptic-parabolic systems [Lass/V. 11] Elliptic-parabolic systems: T = 1, Ω = (a, b) y t (c 1 y) N(y,p,q;µ) = 0 (c p) N(y,p,q;µ) = 0 (c 3 q) + N(y,p,q;µ) = 0 in Q = (0,T) Ω in Q in Q
33 Elliptic-parabolic systems [Lass/V. 11] Elliptic-parabolic systems: T = 1, Ω = (a, b) y t (c 1 y) N(y,p,q;µ) = 0 (c p) N(y,p,q;µ) = 0 (c 3 q) + N(y,p,q;µ) = 0 in Q = (0,T) Ω in Q in Q Parameter-dependent nonlinearity: µ = (µ 1,µ ) 0 N(y,p,q;µ) = µ y sinh(µ1 (q p ln y)) Boundary conditions: y x (t,a) = y x (t,b) = p(t,a) = p x (t,b) = 0, q x (t,a) = q(t,b) = 0 Discretization: FE (nd order) and implicit Euler method Numerical solution method: (damped) Newton algorithm
34 POD basis computation POD criterium: l dim(span {y(t) t [0,T]}) T min y(t) 0 l y(t),ψ i ψ i dt s.t. ψi,ψ j = δ ij i=1 Inner product: L (Ω) or H 1 (Ω) (+b.c.) Solution to optimization problem: Rψ i = T 0 y(t),ψ i y(t)dt = λ i ψ i, i = 1,...,l (Kv i )(t) = T 0 y(t),y( ) v i ds = λ i v i (t), i = 1,...,l Relation via SVD: ψ i = T 0 v i(t)y(t)dt/ λ i Discrete variant: α j = O(N 1 t ) N t y(tj l min α j ) y(t j ),ψ i ψ i j=1 i=1 s.t. ψ i,ψ j = δ ij
35 Compare: POD and SVD Singular Values for POD σ Y σ P σ Q Singular and Eigen Values for Y σ λ σ σ and λ # #
36 References Antoulas, Chapelle, Heinkenschloss, Hinze, Petzold, Sachs... Karhunen-Loéve Decomp., Principal Component Analysis, SVD,... Kunisch/V. 99: Control of Burgers equation by a reduced order approach using POD Kunisch/V. 0: Galerkin POD methods for a general equation in fluid dynamics Lass/V. 11: POD Galerkin schemes for nonlinear elliptic-parabolic systems Müller/V. 06: Model reduction by POD for lambda-omega systems V. 01 : Optimal control of a phase-field model using POD V. 09: Model Reduction using POD. Lecture notes
Proper 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 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 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 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 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 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 informationINTERPRETATION OF PROPER ORTHOGONAL DECOMPOSITION AS SINGULAR VALUE DECOMPOSITION AND HJB-BASED FEEDBACK DESIGN PAPER MSC-506
INTERPRETATION OF PROPER ORTHOGONAL DECOMPOSITION AS SINGULAR VALUE DECOMPOSITION AND HJB-BASED FEEDBACK DESIGN PAPER MSC-506 SIXTEENTH INTERNATIONAL SYMPOSIUM ON MATHEMATICAL THEORY OF NETWORKS AND SYSTEMS
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 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 informationProper Orthogonal Decomposition: Theory and Reduced-Order Modelling
Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling S. Volkwein Lecture Notes (August 7, 13) University of Konstanz Department of Mathematics and Statistics Contents Chapter 1. The POD
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 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 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 informationReduced-Order Methods in PDE Constrained Optimization
Reduced-Order Methods in PDE Constrained Optimization M. Hinze (Hamburg), K. Kunisch (Graz), F. Tro ltzsch (Berlin) and E. Grimm, M. Gubisch, S. Rogg, A. Studinger, S. Trenz, S. Volkwein University of
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationProper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control
Proper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control Michael Hinze 1 and Stefan Volkwein 1 Institut für Numerische Mathematik, TU Dresden,
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
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 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 informationarxiv: v1 [math.na] 13 Sep 2014
Model Order Reduction for Nonlinear Schrödinger Equation B. Karasözen, a,, C. Akkoyunlu b, M. Uzunca c arxiv:49.3995v [math.na] 3 Sep 4 a Department of Mathematics and Institute of Applied Mathematics,
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationLinear Algebra and Dirac Notation, Pt. 3
Linear Algebra and Dirac Notation, Pt. 3 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 3 February 1, 2017 1 / 16
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 - 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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationAsymptotic Stability of POD based Model Predictive Control for a semilinear parabolic PDE
Asymptotic Stability of POD based Model Predictive Control for a semilinear parabolic PDE Alessandro Alla Dipartimento di Matematica, Sapienza Università di Roma, P.le Aldo Moro, 85 Rome, Italy Stefan
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 informationKernel Principal Component Analysis
Kernel Principal Component Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationDedicated to the 65th birthday of Prof Jürgen Sprekels
OPTIMAL CONTROL OF MAGNETIC FIELDS IN FLOW MEASUREMENT SERGE NICAISE, SIMON STINGELIN, AND FREDI TRÖLTZSCH Dedicated to the 65th birthday of Prof Jürgen Sprekels Abstract. Optimal control problems are
More informationSurface x(u, v) and curve α(t) on it given by u(t) & v(t). Math 4140/5530: Differential Geometry
Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt ( ds dt )2 Surface x(u, v)
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
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 informationOPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL
OPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL O. LASS, S. TRENZ, AND S. VOLKWEIN Abstract. In the present paper the authors consider an optimal control
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
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 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 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 informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationOptimal Control. Lecture 18. Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline Hamilton-Jacobi-Bellman (HJB) Equation Iterative solution of HJB Equation
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 informationEfficient Observation of Random Phenomena
Lecture 9 Efficient Observation of Random Phenomena Tokyo Polytechnic University The 1st Century Center of Excellence Program Yukio Tamura POD Proper Orthogonal Decomposition Stochastic Representation
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 informationOn the Karhunen-Loève basis for continuous mechanical systems
for continuous mechanical systems 1 1 Department of Mechanical Engineering Pontifícia Universidade Católica do Rio de Janeiro, Brazil e-mail: rsampaio@mec.puc-rio.br Some perspective to introduce the talk
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
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 information1. The Polar Decomposition
A PERSONAL INTERVIEW WITH THE SINGULAR VALUE DECOMPOSITION MATAN GAVISH Part. Theory. The Polar Decomposition In what follows, F denotes either R or C. The vector space F n is an inner product space with
More informationPOD Model Reduction of Large Scale Geophysical Models. Ionel M. Navon
POD Model Reduction of Large Scale Geophysical Models Ionel M. Navon School of Computational Science Florida State University Tallahassee, Florida Thanks to Prof. D. N. Daescu Dept. of Mathematics and
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More information1 Singular Value Decomposition and Principal Component
Singular Value Decomposition and Principal Component Analysis In these lectures we discuss the SVD and the PCA, two of the most widely used tools in machine learning. Principal Component Analysis (PCA)
More informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Numerical Linear Algebra Background Cho-Jui Hsieh UC Davis May 15, 2018 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationTowards parametric model order reduction for nonlinear PDE systems in networks
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
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
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 informationTopics. CS Advanced Computer Graphics. Differential Geometry Basics. James F. O Brien. Vector and Tensor Fields.
CS 94-3 Advanced Computer Graphics Differential Geometry Basics James F. O Brien Associate Professor U.C. Berkeley Topics Vector and Tensor Fields Divergence, curl, etc. Parametric Curves Tangents, curvature,
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 informationGeneral Relativity I
General Relativity I presented by John T. Whelan The University of Texas at Brownsville whelan@phys.utb.edu LIGO Livingston SURF Lecture 2002 July 5 General Relativity Lectures I. Today (JTW): Special
More informationThe Newton-ADI Method for Large-Scale Algebraic Riccati Equations. Peter Benner.
The Newton-ADI Method for Large-Scale Algebraic Riccati Equations Mathematik in Industrie und Technik Fakultät für Mathematik Peter Benner benner@mathematik.tu-chemnitz.de Sonderforschungsbereich 393 S
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationSTA141C: Big Data & High Performance Statistical Computing
STA141C: Big Data & High Performance Statistical Computing Lecture 5: Numerical Linear Algebra Cho-Jui Hsieh UC Davis April 20, 2017 Linear Algebra Background Vectors A vector has a direction and a magnitude
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationChapter 7: Symmetric Matrices and Quadratic Forms
Chapter 7: Symmetric Matrices and Quadratic Forms (Last Updated: December, 06) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). A few theorems have been moved
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters, transforms,
More informationReduced-Order Model Development and Control Design
Reduced-Order Model Development and Control Design K. Kunisch Department of Mathematics University of Graz, Austria Rome, January 6 Laser-Oberflächenhärtung eines Stahl-Werkstückes [WIAS-Berlin]: workpiece
More informationStabilization of Heat Equation
Stabilization of Heat Equation Mythily Ramaswamy TIFR Centre for Applicable Mathematics, Bangalore, India CIMPA Pre-School, I.I.T Bombay 22 June - July 4, 215 Mythily Ramaswamy Stabilization of Heat Equation
More informationLecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:
tifica Lecture 1: Computation of Karhunen-Loeve Expansion: Alexander Litvinenko http://sri-uq.kaust.edu.sa/ Stochastic PDEs We consider div(κ(x, ω) u) = f (x, ω) in G, u = 0 on G, with stochastic coefficients
More informationChapter 4 Euclid Space
Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationModel Reduction, Centering, and the Karhunen-Loeve Expansion
Model Reduction, Centering, and the Karhunen-Loeve Expansion Sonja Glavaški, Jerrold E. Marsden, and Richard M. Murray 1 Control and Dynamical Systems, 17-81 California Institute of Technology Pasadena,
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationVariational Assimilation of Discrete Navier-Stokes Equations
Variational Assimilation of Discrete Navier-Stokes Equations Souleymane.Kadri-Harouna FLUMINANCE, INRIA Rennes-Bretagne Atlantique Campus universitaire de Beaulieu, 35042 Rennes, France Outline Discretization
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationFull-State Feedback Design for a Multi-Input System
Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C
More informationON CONVERGENCE OF A RECEDING HORIZON METHOD FOR PARABOLIC BOUNDARY CONTROL. fredi tröltzsch and daniel wachsmuth 1
ON CONVERGENCE OF A RECEDING HORIZON METHOD FOR PARABOLIC BOUNDARY CONTROL fredi tröltzsch and daniel wachsmuth Abstract. A method of receding horizon type is considered for a simplified linear-quadratic
More informationMultigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients
Multigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Random fields
More informationTaylor expansions for the HJB equation associated with a bilinear control problem
Taylor expansions for the HJB equation associated with a bilinear control problem Tobias Breiten, Karl Kunisch and Laurent Pfeiffer University of Graz, Austria Rome, June 217 Motivation dragged Brownian
More information