POD for Parametric PDEs and for Optimality Systems
|
|
- Jonathan O’Brien’
- 5 years ago
- Views:
Transcription
1 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, Bonn
2 Motivation : Parameter identification Model equations: ( ) div ( c u) + β u + au = f c u n + qu = g N u = g D in Ω R d on Γ N Γ on Γ D = Γ \ Γ N Problem: estimate parameters (e.g., β or a) in ( ) from given (perturbed) measurements u d for the solution u on (parts of) Γ Mathematical formulation: ( -dim.) optimization problem min α u u d 2 ds + κ p 2 s.t. (p, u) solves ( ) and p P ad Γ Numerical strategy: combine optimization methods with fast (local) rate of convergence and POD model reduction for the PDEs s.t. subject to
3 Motivation 2: Optimal control of time-dependent problems Model problem: min 2 s.t. Adjoint system: Ω y(t) y T 2 dx + κ 2 y t y + f (y) = y Γ = u y() = y T Γ u 2 dxdt p t p + f (y) p =, p Γ =, p(t) = y T y(t) Optimizer: second-order algorithms like SQP or Newton methods Challenge: large-scale fast/real-time optimizer
4 Motivation 3: Closed-loop control for time-dependent PDEs Open-loop control: input u(t) ẏ(t)=f (t,y(t),u(t)) y()=y R l (after spatial discretization) output y(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 equation v t (t, y ) + H(v y (t, y ), y ) = in (, T) R l Strategy: l-dim. spatial approximation by POD model reduction
5 Outline of the talk POD in Hilbert spaces Parameter estimation in elliptic systems POD for optimality systems (OS-POD) Conclusions
6 POD in Hilbert spaces Topology: Hilbert space X with inner product, Snapshots: y,..., y n X Snapshot ensemble: V = span {y,..., y n } X, d = dim V n POD basis of any rank l {,..., d}: with weights α j n yj l 2 min α j y j, ψ i ψ i j= i= s.t. ψ i, ψ j = δ ij Constrained optimization: min J(ψ,..., ψ l ) s.t. ψ i, ψ j = δ ij = { if i = j otherwise s.t. subject to
7 Optimality conditions and computation of POD basis EVP for linear, symmetric R n in X: n R n u i = α j u i, y j y j = λ i u i j= and set ψ i = u i EVP for linear, symmetric K n = (( y j, y i )) in R n : and set ψ i = λi K n v i = λ i v i n α j (v i ) j y j (methods of snapshots) j= Error for the POD basis of rank l: n yj l 2 d α j y j, ψ i ψ i = j= i= i=l+ λ i
8 Properties of the POD basis Uncorrelated POD coefficients: n α j y j, ψ i y j, ψ k = δ ik λ i j= Optimality of the POD basis: l n α j y j, ψ i l n 2 α j y j, χ i 2 i= j= i= j= where {χ i } l i= orthonormal in X
9 POD for λ-ω systems [Müller/V.] PDEs: s = u 2 + v 2, λ(s) = s, ω(s) = βs ( ) ( ) ( ut λ(s) ω(s) u = v t ω(s) λ(s) v Homogeneous boundary conditions: ) + ( σ u σ v ) u = v = or u n = v n = Initial conditions: u (x, x 2 ) = x 2.5, v (x, x 2 ) = (x.5)/2 u for β= and t= u for β=2 and t= u for β=3 and t= u for β= and t=
10 POD basis for λ-ω systems Offsets: ū(x) = n n u(t j, x) or ū j= Snapshots: û j (x) = u(t j, x) ū(x) for j n POD eigenvalue problem: X = L 2 (Ω) Kv i = λv i, i l, with K ij = POD basis computation: ψ i = λi n j= α j(v i ) j û j Ω û j (x)û i (x)dx Decay of the first eigenvalues Decay of the first eigenvalues β= β=.5 β=2 β= β= β=.5 β=2 β=
11 ROM for λ-ω systems POD Galerkin ansatz: u l (t, x) = ū(x) + l u j l (t)ψ j(x), j= Reduced-order model (ROM): insert ansatz into PDEs v l (t, x) = v(x) + l v j l (t)φ j(x) multiply by POD basis functions ψi respectively φ i integrate over Ω j= Numerical results: relative error of u for β=.5 t u l(t) u(t) 2 L 2 (Ω) u(t) 2 L 2 (Ω) 5 l= l=5 l=25 l= t axis
12 Relative POD errors for λ-ω systems Offsets: u m (x) = n Relative POD errors: n u(t j, x) or ū j= ū = ū = u m l = l = l = E rel (u) = np α j u l (t j) u(t j) 2 L j= 2 (Ω) np j= α j u(t j) 2 L 2 (Ω) ū = ū = u m l = l = l = for β =.5 (left) and β = 2 (right)
13 Continuous POD in Hilbert spaces [Henri/Yvon, Kunisch/V.,...] Snapshots: y(µ) X for all µ I Snapshot ensemble: V = {y(µ) µ I} X, d = dim V POD basis of rank l < d: l 2 min y(µ) y(µ), ψ i ψ i dµ s.t. ψi, ψ j = δ ij I i= EVP for linear, symmetric R in X: Rψi = ψi, y(µ) y(µ)dµ = λ i ψi Error for the POD basis of rank l: l y(µ) y(µ), ψi, ψi I I i= 2 dµ = i=l+ λ i
14 Relationship between discrete and continuous POD Operators R n and R: n R n ψ = α j ψ, y(µ j ) y(µ j ) j= Rψ = ψ, y(µ) y(µ)dµ I for ψ X = L 2 (Ω) for ψ X = L 2 (Ω) Operator convergence of R n R: y smooth and appropriate α j s Perturbation theory [Kato]: (λ i, ψ i ) n (λ i, ψ i ) for i l Choice of the weights α j?: ensure convergence R n n R
15 Parameter estimation [Kahlbacher/V.] Model equations: β(x) = ( x x 2 ) 2 u + β u + au = in Ω = (, ) (, ) 2 u n u = on Γ Snapshots: (FE) solutions {u j } 2 j= for a j = j POD basis of rank l: 2 (P l ) min u j j= l u j, ψ i ψ i s.t. ψi, ψ j = δ ij i= with ϕ, φ = Ω ϕφdx and ϕ = ϕ, ϕ Solution to (P l ): correlation matrix K ij = u i, u j Kv i = λ i v i, λ λ 2... λ l, ψ i = λi 2 (v i ) j u j j=
16 λ i Reduced-order modelling (ROM) Ansatz: u l = i l ul i ψ i and Galerkin projection Error estimate: u l (a) u(a) 2 da i>l λ i Exponential decay of the eigenvalues: λ i = λ e η(i ) Experimental order of decay (EOD) [Hinze/V.]: EOD := l max l max Q(l) with l= Q(l) = ln R u l (a) u(a) 2 da R u l+ (a) u(a) 2 da η Decay of Eigenvalues λ i 2 Eigenvalues exponential decay Relative errors for different numbers of POD functions upod(k) ufe(k) / ufe(k) i axis 5 npod= npod=2 npod=3 npod=4 npod= q axis
17 Parameter identification Model equations: g(x) = x ( ) 3 4 u + ( ) u + au = g in Ω = (, ) (, ) 3 u 4 n u = on Γ Data: choose a id and compute (FE) solution u(a id ) to ( ) Reconstruction: estimate a from u d = ( + εδ)u(a id ) Γ with random ε and factor δ = 5% Constrained optimization: min J(a, u) = α u u d 2 ds+κ a 2 s.t. (a, u) solves ( ) and a Relaxation of the inequality: Γ min J λ (a, u) = J(a, u)+ max{, λ + ( a) } 2 s.t. (a, u) solves ( )
18 Global convergent optimization method Outer loop: augmented Lagrangian method control of k and λ k Inner loop: globalized SQP algorithm with fixed ( k, λ k ) for min J k λ k (a, u) s.t. { 3 4 u + ( ) u + au = g in Ω 3 u 4 n u = on Γ Numerical results: α = 5, κ =.5, a id = 25, l = 7 Decay of the first eigenvalues 2 4 relative errors: u l u(a id ) u(a id ) a l a id a id
19 OS-POD [Kunisch/V.] Minimize: J(y, u) = β 2 T s.t. (e.g., Navier-Stokes) y(t) z(t) 2 H dt + 2 T u(t) T Ru(t)dt d dt y(t) + Ay(t) + N(y(t)) = m u k (t)b k in [, T] and k= y() = y H, V Hilbert spaces, V H = H V (e.g., H = L 2, V = H ) R R m m with R, z L 2 (, T;H), β > a : V V R bounded, symmetric, coercive A : V V with Aφ, ϕ V,V = a(φ, ϕ) for all φ, ϕ V N : V V, u L 2 (, T; R m ), y H, b k H
20 POD modelling Choose snapshots, e.g., V = {y(t) t [, T]} for fixed u Compute POD basis ψ,...,ψ l Galerkin ansatz: x(t) = l x i (t)ψ i i= Model reduction low dimension Problems: Quality of the basis for unknown optimal control? Can we avoid the computation of y for various inputs? Can we avoid the computation of p for various observations? Compare: Balanced truncation for ẋ(t) = Ax(t) + Bu(t), t > and x() = x y(t) = Cx(t), t >
21 Optimality System-POD Minimize for fixed l > J l (x, ψ, u) = β 2 s.t. T T x(t) T( E(ψ)x(t) 2z l (t, ψ) ) dt+ 2 u(t) T Ru(t)dt E(ψ)ẋ(t) + A(ψ)x(t) + N(x(t), ψ) = B(ψ)u(t) t, E(ψ)x() = x D E ij = ψ j, ψ i H, A ij = a(ψ j, ψ i ), B ij = b j, ψ i H, N i = N lp j= x j ψ j, ψ i E, z l i = z(t), ψ i H
22 Optimality System-POD Minimize for fixed l > J l (x, ψ, u) = β 2 s.t. T T x(t) T( E(ψ)x(t) 2z l (t, ψ) ) dt+ 2 u(t) T Ru(t)dt E(ψ)ẋ(t) + A(ψ)x(t) + N(x(t), ψ) = B(ψ)u(t) t, E(ψ)x() = x ẏ(t) + Ay(t) + N(y(t)) = (Bu)(t) = m u k (t)b k t, y() = y R(y)ψ i = T k= y(t), ψ i y(t)dt = λ i ψ i, ψ i, ψ j = δ ij D E ij = ψ j, ψ i H, A ij = a(ψ j, ψ i ), B ij = b j, ψ i H, N i = N lp j= x j ψ j, ψ i E, z l i = z(t), ψ i H
23 OS-POD 2 Constraints: E(ψ)ẋ(t) + A(ψ)x(t) + N(x(t), ψ) = B(ψ)u(t) t, E(ψ)x() = x ẏ(t) + Ay(t) + N(y(t)) = m u k (t)b k t, y() = y k= R(y)ψ i = λ i ψ i, ψ i, ψ j = δ ij State variables: w = (x, ψ,...,ψ l, y, λ,...,λ l ) POD basis: ψ i = ψ i (u) computed from trajectory y = y(u) Existence of optimal solutions under assumptions for N and A Existence of Lagrange multipliers provided λ >... > λ l >
24 Optimality system for OS-POD Optimal solution: w = (x, ψ,..., ψ l, y, λ,..., λ l ) Dual equation for x (POD) dynamics: q : [, T] R l, q(t) = E(ψ) q(t) + ( A(ψ) + N x (x(t), ψ) T) q(t) = β ( z l (t, ψ) E(ψ)x(t) ) Dual equations for y dynamics: p : [, T] V, p(t) = ṗ(t) + (A + N (y(t)) )p(t) = l ) ( y(t), µ i ψ i + y(t), ψ i µ i i= Optimality condition: Ru(t) = B(ψ) T q(t) + B p(t) Right-hand side: (R λ i I)µ i = G i (w) ker (R λ i I), i l
25 Boundary control of the Burgers equation Consider s.t. min J(y, u) = 2 T ( y z 2 dx + β ( u(t) 2 + v(t) 2) ) dt Ω y t νy xx + yy x = f in Q = (, T) Ω νy x (, ) = u in (, T) νy x (, ) + σy(, ) = v in (, T) y(, ) = y in Ω = (, ) Desired state z=z(t,x).5 T =, β =., ν =.5, f (t, x) = e 3t sin(2πx), σ =., y = sin(2πx) x axis.2.4 t axis.6.8
26 Numerical strategy: Update of POD basis Idea: combine fast POD solver and OS-POD information ) Choose u and set n = 2) Determine PODs {ψi n}l i= for yn = y n (u n ) 3) Solve inexactly [ũ n, q n ] = SQP l (m) (POD solver) 4) Solve dual quation for p n (OS-POD information) 5) Compute u n+ = ũ n τ ( Rũ n (t) B(ψ n ) T q n (t) B p n (t) ), τ > 6) Set n = n + and go back to 2) if n n max Stopping criterium: OS-POD gradient small, POD solver exact SQP l (m): m SQP steps for the control problem with fixed POD basis
27 Numerical test (figures) Finite element solution y=y(t,x) POD basis ψ POD basis ψ Initial condition y x axis 2.5 x axis.5 POD basis ψ 3 4 POD basis ψ x axis x axis.2.4 t axis x axis x axis Desired state z=z(t,x) POD optimal state y * =y * (t,x) POD basis ψ POD basis ψ x axis 2.5 x axis.5 2 POD basis ψ 3 4 POD basis ψ x axis t axis x axis.2.4 t axis x axis x axis
28 Numerical test (tables) J(y, u, v) Solve for u = v =.2234 OS-POD.383 FE-SQP.3765 Steps Generate snapshots Determine POD Determine ROM SQP solver Compute µ i s Dual FE solver Backtracking in 5) CPU time.62 s.9 s.3 s 4.6 s 2.7 s.7 s.32 s n = n = 4 with FE controls λ /tr (K h ) λ 2 /tr (K h ) λ 3 /tr (K h ) λ 4 /tr (K h )
29 Conclusions POD in Hilbert spaces: choice of weights and norms convergence estimates [Kunisch/V., Hinze/V., Kahlbacher/V.] Parameter estimation in elliptic systems Helmholtz equation [ACC Graz] OS-POD: update of the PODs within optimization optimal POD basis at optimal solution more complex problems Preprints: POD scriptum:
Proper 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 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) 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationHamburger Beiträge zur Angewandten Mathematik
Hamurger Beiträge zur Angewandten Mathematik POD asis updates for nonlinear PDE control Carmen Gräßle, Martin Guisch, Simone Metzdorf, Sarina Rogg, Stefan Volkwein Nr. 16-18 August 16 C. Gräßle, M. Guisch,
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 informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 informationAn Inexact Newton Method for Nonlinear Constrained Optimization
An Inexact Newton Method for Nonlinear Constrained Optimization Frank E. Curtis Numerical Analysis Seminar, January 23, 2009 Outline Motivation and background Algorithm development and theoretical results
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
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 informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationOn feedback stabilizability of time-delay systems in Banach spaces
On feedback stabilizability of time-delay systems in Banach spaces S. Hadd and Q.-C. Zhong q.zhong@liv.ac.uk Dept. of Electrical Eng. & Electronics The University of Liverpool United Kingdom Outline Background
More informationSECOND-ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR THE OPTIMAL CONTROL OF NAVIER-STOKES EQUATIONS
1 SECOND-ORDER SUFFICIENT OPTIMALITY CONDITIONS FOR THE OPTIMAL CONTROL OF NAVIER-STOKES EQUATIONS fredi tröltzsch and daniel wachsmuth 1 Abstract. In this paper sufficient optimality conditions are established
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationPDE Constrained Optimization selected Proofs
PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5
More informationQuestion 9: PDEs Given the function f(x, y), consider the problem: = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1. x 2 u. u(x, 0) = u(x, 1) = 0 for 0 x 1
Question 9: PDEs Given the function f(x, y), consider the problem: 2 u x 2 u = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1 u(x, 0) = u(x, 1) = 0 for 0 x 1 u(0, y) = u(1, y) = 0 for 0 y 1. a. Discuss how you
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 informationNumerics and Control of PDEs Lecture 7. IFCAM IISc Bangalore. Feedback stabilization of a 1D nonlinear model
1/3 Numerics and Control of PDEs Lecture 7 IFCAM IISc Bangalore July August, 13 Feedback stabilization of a 1D nonlinear model Mythily R., Praveen C., Jean-Pierre R. /3 Plan of lecture 7 1. The nonlinear
More informationExponential stability of abstract evolution equations with time delay feedback
Exponential stability of abstract evolution equations with time delay feedback Cristina Pignotti University of L Aquila Cortona, June 22, 2016 Cristina Pignotti (L Aquila) Abstract evolutions equations
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
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 information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
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 informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationDRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS
DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
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 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 informationAsymptotic distribution of eigenvalues of Laplace operator
Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour
More informationK. Krumbiegel I. Neitzel A. Rösch
SUFFICIENT OPTIMALITY CONDITIONS FOR THE MOREAU-YOSIDA-TYPE REGULARIZATION CONCEPT APPLIED TO SEMILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS WITH POINTWISE STATE CONSTRAINTS K. Krumbiegel I. Neitzel A. Rösch
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More information1. INTRODUCTION 2. PROBLEM FORMULATION ROMAI J., 6, 2(2010), 1 13
Contents 1 A product formula approach to an inverse problem governed by nonlinear phase-field transition system. Case 1D Tommaso Benincasa, Costică Moroşanu 1 v ROMAI J., 6, 2(21), 1 13 A PRODUCT FORMULA
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 informationStability, the Maslov Index, and Spatial Dynamics
Stability, the Maslov Index, and Spatial Dynamics Margaret Beck Boston University joint work with Graham Cox (MUN), Chris Jones (UNC), Yuri Latushkin (Missouri), Kelly McQuighan (Google), Alim Sukhtayev
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 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 informationCriterions on periodic feedback stabilization for some evolution equations
Criterions on periodic feedback stabilization for some evolution equations School of Mathematics and Statistics, Wuhan University, P. R. China (Joint work with Yashan Xu, Fudan University) Toulouse, June,
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationInvestigating advection control in parabolic PDE systems for competitive populations
Investigating advection control in parabolic PDE systems for competitive populations Kokum De Silva Department of Mathematics University of Peradeniya, Sri Lanka Tuoc Phan and Suzanne Lenhart Department
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
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 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 informationSpace time finite element methods in biomedical applications
Space time finite element methods in biomedical applications Olaf Steinbach Institut für Angewandte Mathematik, TU Graz http://www.applied.math.tugraz.at SFB Mathematical Optimization and Applications
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 informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
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 informationReduced basis method for efficient model order reduction of optimal control problems
Reduced basis method for efficient model order reduction of optimal control problems Laura Iapichino 1 joint work with Giulia Fabrini 2 and Stefan Volkwein 2 1 Department of Mathematics and Computer Science,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationNon-Differentiable Embedding of Lagrangian structures
Non-Differentiable Embedding of Lagrangian structures Isabelle Greff Joint work with J. Cresson Université de Pau et des Pays de l Adour CNAM, Paris, April, 22nd 2010 Position of the problem 1. Example
More informationSystems of Linear Differential Equations Chapter 7
Systems of Linear Differential Equations Chapter 7 Doreen De Leon Department of Mathematics, California State University, Fresno June 22, 25 Motivating Examples: Applications of Systems of First Order
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 informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationTechnische Universität Berlin
Technische Universität Berlin Institut für Mathematik Regularity of the adjoint state for the instationary Navier-Stokes equations Arnd Rösch, Daniel Wachsmuth Preprint 1-4 Preprint-Reihe des Instituts
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
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 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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationALADIN An Algorithm for Distributed Non-Convex Optimization and Control
ALADIN An Algorithm for Distributed Non-Convex Optimization and Control Boris Houska, Yuning Jiang, Janick Frasch, Rien Quirynen, Dimitris Kouzoupis, Moritz Diehl ShanghaiTech University, University of
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 informationMATH4406 (Control Theory) Unit 1: Introduction Prepared by Yoni Nazarathy, July 21, 2012
MATH4406 (Control Theory) Unit 1: Introduction Prepared by Yoni Nazarathy, July 21, 2012 Unit Outline Introduction to the course: Course goals, assessment, etc... What is Control Theory A bit of jargon,
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationOn Controllability of Linear Systems 1
On Controllability of Linear Systems 1 M.T.Nair Department of Mathematics, IIT Madras Abstract In this article we discuss some issues related to the observability and controllability of linear systems.
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 informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More information