Model order reduction & PDE constrained optimization
|
|
- Barbra Lewis
- 5 years ago
- Views:
Transcription
1 1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01
2 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell, Glavaski, Graham, Gunzburger, Hinze, Iannou, Ito, King, Kunisch, Lall, Leibfritz, Le Letty, Marsden, Noack, Patera, Peraire, Prudhomme, Ravindran, Sachs, Schroeter, Serban, Sorensen, Tang, Volkwein, Willcox,...
3 3 Outline Problem formulation Model reduction Validation of reduced model Optimization using the reduced model Basic adaptive approach Numerical experiments Recent developments
4 4 Motivation We have a validated mathematical model for physical process (here a pde system) We intend to use this model to tailor and/or optimize the physical process.
5 5 Motivation: -dimensional optimization problem with pde constraints min (y,u) W U ad J(y, u) s.t. y + Ay + G(y) t = Bu in Z y(0) = y 0 in H. Central tasks: Develop solution strategies which obey the rule Effort of optimization = K Effort of simulation with K small, Propose surrogate models for the pde and quantify their errors, Present a complete (numerical) analysis.
6 6 Examples of pde systems 1 Heat equation: A :=. y + Ay + G(y) t = Bu in Z y(0) = y 0 in H. Burgers: A :=, G(y) := yy, 3 Ignition (Bratu): A :=, G(y) := δe y, δ > 0, 4 Navier Stokes: A := P, G(y) := P[(y )y], P Leray projector, 5 Boussinesq [ Approximation: P Gr g A := 0 ], G(y) = G(v, θ) := [ P[(v )v] (v )θ ].
7 7 DOF diagram Spatial Discretization DOF for Moving Horizon Approach DOF for full optimization DOF for Moving Horizon combined with Model Reduction DOF for Model Reduction Approach
8 8 System reduction for model dynamics Navier-Stokes equations y + (y )y ν y + p = f t in Q = (0, T ) Ω, div y = 0 in Q, y(t, ) = g y(0, ) = y 0 in Ω. on Σ = (0, T ) Ω, Aim: Reduced description of the Navier-Stokes equations α + Aα + n(α) = r in (0, T ) α(0) = a 0
9 1. Construction and validation of the reduced model MOR & PDE constrained optimization 9
10 10 System reduction: Expansions w.r.t. base flows Let ȳ denote a base flow and Φ i, i = 1,..., n Modes. Ansatz for the flow: n y = ȳ + α i Φ i i=1 Possibilities: ȳ stationary solution of Navier-Stokes system, Φ i eigenfunctions of the Navier-Stokes system linearized at ȳ. ȳ mean value of instationary Navier-Stokes solution, Φ i eigenfunctions of the Navier-Stokes system linearized at ȳ. ȳ mean value of instationary Navier-Stokes solution, Φ i normalized Modes obtained from snapshot form of Proper Orthogonal Decomposition.
11 11 Snapshot form of POD Let s take snapshots:
12 1 POD with Snapshots Let y 1,..., y n denote an ensemble of snapshots (of the flow or the dynamical system). Build mean ȳ and modes Φ i as follows: n 1 Compute mean ȳ = 1 y n i i=1 Build correlation matrix K = k ij, k ij = y i ȳ, y j ȳ 3 Compute eigenvalues λ 1,..., λ n and eigenvectors v 1,..., v n of K 4 Define modes Φ i := n j=1 5 Normalize modes Φ i = Φi Φ i Properties: v i j (y j ȳ) The modes are pairwise orthogonal w.r.t. inner product, No other basis can contain more information in fewer elements (Information w.r.t. the norm induced by, ).
13 MOR & PDE constrained optimization 13 First 10 Modes containing % of the information Vorticity Y Vorticity Y X Y Vorticity Vorticity Vorticity X Vorticity Y Y X Vorticity X - 6 Y Y X Vorticity X - 6 Y Y X Vorticity X - 6 Y Vorticity X X 6 7 8
14 14 Galerkin projection Ansatz for the flow n y = ȳ + α i Φ i i=1 Galerkin method with basis Φ 1,..., Φ n yields reduced system α + Aα + n(α) = r α(0) = a 0. Here,, denotes the L inner product. n A = ( ) n a i,j i=1, a i,j = ν Φ i Φ j dx, n(α) = (y y)φ i dx r = ν Ω Ω n ȳ Φ i + f Φ i dx Note that Φ 1,..., Φ n are solenoidal. i=1 Ω c n and a 0 = y 0 Φ i dx Ω i=1. i=1,
15 MOR & PDE constrained optimization 15 Long time behaviour of the POD model Amplitudes of the first mode in [34,44] when using N modes in the POD model N=4 N=8 N=10 N=16 N=0 amplitude time
16 MOR & PDE constrained optimization 16 Cylinder flow at Re = 100, reduced versus full model, 50 snapshots 4 Amplitudes of the first four modes (DS1 with 50 equidistant snapshots on [0,6,8]) projections predictions amplitude time 6 8
17 17 Error estimate (Kunisch and Volkwein) Let y(t 1 ),..., y(t n) denote snapshots taken on an equidistant time grid of [0, T ] with gridsize δt. Let λ 1 > > λ d > 0 denote the strictly positive eigenvalues of the correlation matrix K. For l d let V l = Φ 1,..., Φ l. Further set Then l Y k := α i (t k )Φ i. i=1 n d δt Y i y(t i ) H C y 0, Φ i V + 1 d δt λ i + δt. i=1 i=l+1 i=l+1 This result also extends to the case of distinguish time and snapshot grids. Improvements of reduced models and error estimate by different weighting of snapshots.
18 . Optimization with reduced model MOR & PDE constrained optimization 18
19 19 System reduction & control Model optimization problem: min J(y, u) := 1 y z γ dxdt + (y,u) W U Q o Q c u dxdt s.t. y + (y )y ν y + p t = Bu in Q = (0, T ) Ω, div y = 0 in Q, y(t, ) = g on Σ = (0, T ) Ω, y(0, ) = y 0 in Ω. Q o := (0, T ) Ω o Observation cylinder, Q c := (0, T ) Ω c control cylinder, B control operator
20 0 POD model as pde surrogate in the optimization problem Ansatz for state, control and desired state n n n y = ȳ + α i Φ i, u = β i Φ i, z = ȳ + i=1 i=1 i=1 α z i Φ i. Optimization problem with surrogate model (reduced order model) min (y,u) J(y, u) = J(α, β) = 1 T (α α z ) t M 1 (α α z ) dt + γ T 0 0 β t Mβ dt s.t. α + Aα + n(α) = r + Mβ, α(0) = a 0.
21 1 Validity of surrogate model Fact: Control changes system dynamics. Consequence: Mean and modes should be suitably modified during the optimization process. Idea: Adaptively modify the surrogate model and thus, the reduced optimization problems.
22 Adaptive POD control Afanasiev, Hinze Snapshots y 0 i, i = 1,..., N 0 given, u 0 given control, δ [0, 1] required relative information content, j=0. Compute M = argmin { I(M) := } M λ k / N λ k ; I(M) δ. k=1 k=1 3 Compute POD modes and solve min J(α, β) (ROM) s.t. α + Aα = n(α) + Mβ, α(0) = a 0. for u j = u j (β j ). 4 Compute y j corresponding to u j and new snapshots y j+1 i, i = N j + 1,..., N j+1 to the snapshot set y j i, i = 1,..., N j. 5 While u j+1 u j is large, j = j+1 and goto.
23 3 Numerical comparison Flow around a circular cylinder at Re=100. Control gain: Tracking of Stokes flow (or mean flow) ȳ in an observation volume Ω obs behind the cylinder by applying a volume force in the control volume Ω c. Cost functional: J(y, u) = γ T u 1 T dxdt + 0 Ω c 0 Ω obs y ȳ dxdt CPU time needed to compute the suboptimal controls 40 times smaller than that needed to compute the optimal open loop control. But the quality of the controls is very similar. Runtime(Optimization Problem) = 6 8 Runtime(PDE)
24 4 Uncontrolled flow, target flow = mean flow, controlled flow at t = 3.4.
25 MOR & PDE constrained optimization 5 Numerical results cont. Cost, Ω o = Ω = Ω c, control in space of fluctuations, tracking of Stokes flow 10 8 Evolution of cost 1 iteration iteration 3 iteration 4 iteration 5 iteration uncontrolled 1 grad. iteration (optimal) optimal control Time 3
26 MOR & PDE constrained optimization 6 Numerical results cont. Cost, Ω o = Ω = Ω c, control in the space of fluctuations, tracking of mean flow J 10 unkontrolliert 1 Iteration (ISO) Iteration 3 Iteration Iteration 5 Iteration 6 Iteration 7 Iteration (ISO) optimale Kontrolle (NM) optimale Kontrolle (GM) Zeit
27 MOR & PDE constrained optimization 7 Numerical results cont. Control cost, Ω o = Ω = Ω c, control in the space of fluctuations, tracking of mean flow Iteration (ISO) Iteration (ISO) 3 Iteration (ISO) 4 Iteration (ISO) 5 Iteration (ISO) 6 Iteration (ISO) 7 Iteration (ISO) optimale Kontrolle (NM) 10 0 J f Phase 1 0 Phase 4 6
28 8 Recent developments TRPOD by Arian, Fahl and Sachs 000 Idea: Use a POD surrogate model as model function in the Trust Region process. Let J(u) = J(y(u), u), Ĵ(u) = J(ŷ(u), u), with ŷ(u) the response of the POD surrogate model. Pseudo Algorithm: 1 Given u, compute POD model Compute s = argmin u s Ĵ(u + s) 3 ρ := J(u + s ) J(u) Ĵ(u + s ) Ĵ(u) large: u = u + s, increase moderate: u = u + s, decrease small: keep u, decrease Global convergence under standard TR assumptions plus J (u) Ĵ (u) Ĵ (u) sufficiently small POD of the adjoint equation required.
29 9 Recent developments OSPOD by Kunisch and Volkwein 006 Idea: Include choice of trajectory dependent POD modes as subsidiary condition into the optimization problem. This reads min α,φ,u Ĵ(α, Φ, u) s.t. M(Φ) α + A(Φ)α + n(φ)(α) = B(Φ)u, M(Φ)α(0) = α 0 (Φ), (POSPOD l ) y t + Ay + G(y) = Bu, y(0) = y 0, R(y)Φ i = λ i Φ i for i = 1,..., l, Φ i X = 1 for i = 1,..., l. Here Φ = [Φ 1,..., Φ l ], y l = l α i (t)φ, and i=1 T R(y)(z) := y(t), z X y(t)dt for z X. 0
30 30 Upcoming Numerical analysis of POD Galerkin approximations Error estimates for POD approximations to pde constrained optimal control problems Goal oriented concepts in POD Future perspectives Break
31 31 Outline Motivation: PDE constrained optimization Mathematical setting Construction of the POD spaces Tailored discrete concept for optimal control problem and numerical analysis Adaptive concepts - why?
32 3 Motivation: optimization problem with pde constraints min (y,u) W U ad J(y, u) s.t. y + Ay + G(y) t = Bu in Z y(0) = y 0 in H. Approach: Solve this problem by using a POD surrogate model; min J l (y l, u l ) s.t. (y l,u l ) W l U ad Tasks: Error estimation, y l t + Al y l + G l (y l ) = Bu l in (Z l ) y l (0) = y l 0 in Hl. adaption of the POD surrogate model during the optimization loop.
33 33 Mathematical setting, state equation V, H separable Hilbert spaces, (V, H = H, V ) Gelfand triple. a : V V R bounded, coercive and symmetric. Set, V := a(, ). U Hilbert space, B : U L (U, L (V )) linear control operator, y 0 H. State equation d dt (y(t), v) H + a(y(t), v) = (Bu)(t), v V,V, t [0, T ], v V, (y(0), v) H = (y 0, v) H, v V. For every u U the solution y = y(u) W := {w L (V ), w t L (V )} is unique.
34 34 Optimization problem Cost functional J(y, u) := 1 y z L (H) + α u U. Admissibility: u U ad U closed, convex, y y(u) unique solution of state equation associated to u, i.e. d dt (y(t), v) H + a(y(t), v) = (Bu)(t), v V,V, t [0, T ], v V, (y(0), v) H = (y 0, v) H, v V. Minimization problem: (P) min J(y, u) s.t. Admissibility. (y,u) W (0,T ) U ad (P) admits a unique solution (y, u) W U ad.
35 35 Optimality conditions With the reduced cost functional Ĵ(u) := J(y(u), u) there holds (Ĵ (u), v u) 0 for all v U ad. Here Ĵ (u) = αu + B p(y(u)). The function p solves the adjoint equation d dt (p(t), v) H + a(v, p(t)) = (y z, v) H, t [0, T ], v V, (p(t ), v) H = 0, v V. Variational inequality equivalent to nonsmooth operator equation u = P Uad ( 1 ) α B p(y(u)) with P Uad denoting the orthogonal projection onto U ad.
36 36 Construction of the POD Galerkin spaces Let (X = H or X = V ) T Y : L (0, T ) X, Yϕ = ϕ(t)y(t) dt for ϕ L (0, T ). 0 Then Y : X L (0, T ), ( Y z ) (t) = z, y(t) X for z X and t [0, T ], and T R = YY : X X, Rz = z, y(t) X y(t) dt for z X, 0 K = Y Y : L (0, T ) L (0, T ), ( Kϕ ) T (t) = y(s), y(t) X ϕ(s) ds for ϕ L (0, T ) 0 K continuous version of the correlation matrix, and K is compact, R is compact.
37 37 Construction of the POD spaces cont. Right-eigenfunctions: Kv i = λ i v i, λ 1 λ..., and λ i 0 as i Left eigenfunctions: Rψ i = λ i ψ i, i = 1,,... with {ψ i } i N a complete orthonormal basis for X and non negative eigenvalues {λ i } i N. Eigenvalues contain information: T y(t) X dt = λ i and y 0 X = y0, ψ i X. 0 i=1 i=1
38 38 Discrete concept for the state equation For l N choose a POD subspace V l := χ 1,..., χ l of V with the property l y(t) (y(t), χ k ) V χ k L (0,T ;V ) λ k. k=1 k=l+1 Galerkin semi-discretization y l of state y using subspace V l : ( d dt y l (t), v ) H + a(y l (t), v) = (Bu)(t), v V,V, t [0, T ], v V l, (y(0), v) H = (y 0, v) H, v V l. If needed, define similarly a Galerkin semi-discretization p l of p; d ( dt p l (t), v ) H + a(v, pl (t)) = ( y l z, v ) H, t [0, T ], v V l, ( p l (T ), v ) = 0, v V l. H
39 39 Optimization problem with POD surrogate model Discrete minimization problem: ( ˆP l ) min u U ad Ĵ l (u) := J(y l (u), u). ( ˆP l ) admits a unique solution u l U ad. Optimality condition: (Ĵl (u), v u l ) 0 for all v U ad. Here Ĵ l (u) = αu + B p l (y l (u)). The function p l solves the adjoint equation d ( dt p l (t), v ) H + a(v, pl (t)) = ( y l z, v ) H, t [0, T ], v V l, ( p l (T ), v ) = 0, v V l. H Variational inequality equivalent to nonsmooth operator equation u l = P Uad ( 1 ) α B p l (y l (u)).
40 40 Error estimate Theorem: Let u, u l denote the unique solutions of (P) and ( ˆP l ), respectively. Then u u l U 1 {( ) B (p(y(u)) p l (y(u))), u l u α + U T + 0 ( ) y l (u l ) y l (u), y(u) y l (u) dt H Using the analysis of Kunisch and Volkwein for POD approximations one gets u u l U y 0 P l y 0 H + λ k + k=l+1 + y t P l y t L (0,T ;V ) + p(y(u)) Pl (p(y(u))) W (0,T )
41 41 Conclusions from the analysis Get rid of (y P l y) t L (0,T ;V ) include derivative information into your snapshot set. Get rid of p P l p W (0,T ) include adjoint information into your snapshot set. Recipe: For l N choose a POD subspace V l := χ 1,..., χ l of V with the property l y(t) (y(t), χ k ) V χ k W (0,T ) λ k, k=1 k=l+1 and if one intends to solve optimization problems, also ensure l p(t) (p(t), χ k ) V χ k W (0,T ) λ k, k=1 k=l+1
42 41 Conclusions from the analysis Get rid of (y P l y) t L (0,T ;V ) include derivative information into your snapshot set. Get rid of p P l p W (0,T ) include adjoint information into your snapshot set. Recipe: For l N choose a POD subspace V l := χ 1,..., χ l of V with the property l y(t) (y(t), χ k ) V χ k W (0,T ) λ k, k=1 k=l+1 and if one intends to solve optimization problems, also ensure l p(t) (p(t), χ k ) V χ k W (0,T ) λ k, k=1 k=l+1
43 4 Numerical confirmation of approximation order of the states Linear heat equation with y 0 0 and inhomogeneous boundary data. FE-solution {y h (t j )} m j=0 computed on equi distant time grid. Snapshots: y h (t j 1 ) for 1 j m + 1, y j = y h (t j m 1 ) y h (t j m ) for m + j m + 1. t Correlation matrix (k ij ) m+1 i,j=1, k ij = y i, y j V Expected decay of its eigenvalues: Experimental order of decay: λ i = λ 1 e α(i 1) for i 1. Q(l) = ln y l y L (0,T ;X ) y l+1 y EOD := 1 l max Q(k) α. l L (0,T ;X ) max k=1
44 43 Decay of eigenvalues and of norms 10 0 Decay of the eigenvalues and estimated rates Decay of the norms L (H 1 (Ω)) error L (L (Ω)) error λ i λ 1 e α (i 1) mit L (0,T;H 1 (Ω)) λ 1 e α (i 1) mit L (0,T;L (Ω)) i axis Number of POD basis functions
45 44 Optimization problem with POD surrogate model Discrete minimization problem: ( ˆP l ) min u U ad Ĵ l (u) := J(y l (u), u). ( ˆP l ) admits a unique solution u l U ad. Optimality condition: (Ĵl (u), v u l ) 0 for all v U ad. Here Ĵ l (u) = αu + B p l (y l (u)). The function p l solves the adjoint equation d ( dt p l (t), v ) H + a(v, pl (t)) = ( y l z, v ) H, t [0, T ], v V l, ( p l (T ), v ) = 0, v V l. H Variational inequality equivalent to nonsmooth operator equation u l = P Uad ( 1 ) α B p l (y l (u)).
46 45 Error estimate Theorem: Let u, u l denote the unique solutions of (P) and ( ˆP l ), respectively. Then u u l U 1 {( ) B (p(y(u)) p l (y(u))), u l u α + U T + 0 ( ) y l (u l ) y l (u), y(u) y l (u) dt H Consequently, if u y l (u) is Lipschitz continuous w.r.t. the L (H)-norm, there holds u u l U y 0 P l y 0 H + λ k + k=l+1 + y t P l y t L (0,T ;V ) + p(y(u)) Pl (p(y(u))) W (0,T )
47 46 Proof of error estimate Add necessary optimality conditions and use u, u l as test functions: (αu l + B p l (y l (u)), v u l ) 0 for all v U ad. Set v = u (αu + B p(y(u)), v u) 0 for all v U ad. Set v = u l α u u l U (B (p p l (y(u)) + (p l (y(u)) p l (y l (u)))+ ) + (p l (y l (u)) p l )), u l u =: U 1. gives first term in estimate.. gives second term in estimate 3. 0.
48 47 Error between state and and its orthogonal projection 10 0 y P l y L (0,T;H 1 (Ω)) 10 0 y P l y L (0,T;L (Ω)) {ψ 1 } {ψ } {ψ 3 } Number l of POD basis functions {ψ 1 } {ψ } {ψ 3 } Number l of POD basis functions
49 48 Error between co-state and its orthogonal projection 10 0 p P l p L (0,T;H 1 (Ω)) 10 p P l p L (0,T;L (Ω)) {ψ 1 } {ψ } {ψ 3 } Number l of POD basis functions {ψ 1 } {ψ } {ψ 3 } Number l of POD basis functions
50 49 Neumann boundary control of the heat equation Comparison of the optimal controls Comparison of u h u u h u 15 for ψ i 1 u 15 for ψ i u 15 for ψ i t axis 0. u h u 15 for ψ i 1 u h u 15 for ψ i u h u 15 for ψ i t axis
51 50 Recent developments TRPOD by Arian, Fahl and Sachs 000 Idea: Use a POD surrogate model as model function in the Trust Region process. Let J(u) = J(y(u), u), Ĵ(u) = J(ŷ(u), u), with ŷ(u) the response of the POD surrogate model. Pseudo Algorithm: 1 Given u, compute POD model Compute s = argmin u s Ĵ(u + s) 3 ρ := J(u + s ) J(u) Ĵ(u + s ) Ĵ(u) large: u = u + s, increase moderate: u = u + s, decrease small: keep u, decrease Global convergence under standard TR assumptions plus J (u) Ĵ (u) Ĵ (u) sufficiently small.
52 51 Recent developments OSPOD by Kunisch and Volkwein 006 Idea: Include choice of trajectory dependent POD modes as subsidiary condition into the optimization problem. This reads min α,φ,u Ĵ(α, Φ, u) s.t. M(Φ) α + A(Φ)α + n(φ)(α) = B(Φ)u, M(Φ)α(0) = α 0 (Φ), (POSPOD l ) y t + Ay + G(y) = Bu, y(0) = y 0, R(y)Φ i = λ i Φ i for i = 1,..., l, Φ i X = 1 for i = 1,..., l. Here Φ = [Φ 1,..., Φ l ], y l = l α i (t)φ, and i=1 T R(y)(z) := y(t), z X y(t)dt for z X. 0 A very similar approach is proposed by Ghattas, van Bloemen Waanders and Willcox 005.
53 5 How many snapshots? Meyer, Matthies, Heuveline, H. How many snapshots? iterative goal oriented procedure. Goal: Resolve J(y) Start on coarse equi-distant time grid and compute snapshots Build POD model and compute y h and adjoint z h of reduced dynamics Becker and Rannacher: J(y) J(y h ) η(y h, z h ) η(y h, z h ) > tol: double number of snapshots (re-computation)
54 53 Where to take snapshots? Where to take snapshots? time-step adaption via sensitivity of POD model. Goal: Optimal time-grid for system dynamics Start on coarse (equi-distant) time grid and compute snapshots Build POD model and compute y h and adjoint z h of reduced dynamics Becker, Johnson, Rannacher: η(y h, z h ) = I j ρ loc j (y h )ωj loc (z h ) New time-grid: equi-distribute ρ loc j (y h )ωj loc (z h )
55 54 Further developments and improvements Efficient treatment of nonlinearities Chaturantabut, Sorensen (010) MOR for the input output map Heiland, Mehrmann A posteriori POD concept Tröltzsch and Volkwein (010) Which modes? DWR concepts (Matthies, Meyer 003) How many snapshots? iterative goal oriented DWR procedure Were take snapshots? time step adaption via sensitivity of the POD model POD in the context of space mapping Sampling of parameter ( control) space Greedy sampling by Patera and Rozza (007) Use of linear MOR techniques for nonlinear problems SQP context, semi linear time integration, domain decomposition Thank you for your attention
56 54 Further developments and improvements Efficient treatment of nonlinearities Chaturantabut, Sorensen (010) MOR for the input output map Heiland, Mehrmann A posteriori POD concept Tröltzsch and Volkwein (010) Which modes? DWR concepts (Matthies, Meyer 003) How many snapshots? iterative goal oriented DWR procedure Were take snapshots? time step adaption via sensitivity of the POD model POD in the context of space mapping Sampling of parameter ( control) space Greedy sampling by Patera and Rozza (007) Use of linear MOR techniques for nonlinear problems SQP context, semi linear time integration, domain decomposition Thank you for your attention
Proper 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 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 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 (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 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 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 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 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 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 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 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 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 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 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 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 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 informationModel Reduction via Proper Orthogonal Decomposition
Model Reduction via Proper Orthogonal Decomposition René Pinnau Fachbereich Mathematik, Technische Universität Kaiserslautern, D-67663 Kaiserslautern, Germany pinnau@mathematik.uni-kl.de 1 Introduction
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 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 (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 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 informationBang bang control of elliptic and parabolic PDEs
1/26 Bang bang control of elliptic and parabolic PDEs Michael Hinze (joint work with Nikolaus von Daniels & Klaus Deckelnick) Jackson, July 24, 2018 2/26 Model problem (P) min J(y, u) = 1 u U ad,y Y ad
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 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 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 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 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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik A finite element approximation to elliptic control problems in the presence of control and state constraints Klaus Deckelnick and Michael Hinze Nr. 2007-0
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 informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Klaus Deckelnick and Michael Hinze A note on the approximation of elliptic control problems with bang-bang
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 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 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 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 informationSTATE-CONSTRAINED OPTIMAL CONTROL OF THE THREE-DIMENSIONAL STATIONARY NAVIER-STOKES EQUATIONS. 1. Introduction
STATE-CONSTRAINED OPTIMAL CONTROL OF THE THREE-DIMENSIONAL STATIONARY NAVIER-STOKES EQUATIONS J. C. DE LOS REYES AND R. GRIESSE Abstract. In this paper, an optimal control problem for the stationary Navier-
More informationNECESSARY OPTIMALITY CONDITIONS FOR AN OPTIMAL CONTROL PROBLEM OF 2D
NECESSARY OPTIMALITY CONDITIONS FOR AN OPTIMAL CONTROL PROBLEM OF 2D g-navier-stokes EQUATIONS Nguyen Duc Loc 1 Abstract Considered here is the optimal control problem of 2D g-navier-stokes equations.
More informationarxiv: v1 [math.oc] 5 Jul 2017
Variational discretization of a control-constrained parabolic bang-bang optimal control problem Nikolaus von Daniels Michael Hinze arxiv:1707.01454v1 [math.oc] 5 Jul 2017 December 8, 2017 Abstract: We
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 informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationImprovement of Reduced Order Modeling based on Proper Orthogonal Decomposition
ICCFD5, Seoul, Korea, July 7-11, 28 p. 1 Improvement of Reduced Order Modeling based on Proper Orthogonal Decomposition Michel Bergmann, Charles-Henri Bruneau & Angelo Iollo Michel.Bergmann@inria.fr http://www.math.u-bordeaux.fr/
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 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 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 informationHessian-based model reduction for large-scale systems with initial condition inputs
Hessian-based model reduction for large-scale systems with initial condition inputs O. Bashir 1, K. Willcox 1,, O. Ghattas 2, B. van Bloemen Waanders 3, J. Hill 3 1 Massachusetts Institute of Technology,
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationSonderforschungsbereich 609 Elektromagnetische Strömungsbeeinflussung in Metallurgie, Kristallzüchtung und Elektrochemie
SFB 69 Sonderforschungsbereich 69 Elektromagnetische Strömungsbeeinflussung in Metallurgie, Kristallzüchtung und Elektrochemie M. Hinze, O. Pfeiffer Active Closed Loop Control Of Weakly Conductive Fluids.
More informationOptimal control of the cylinder wake in the laminar regime by Trust-Region methods and POD Reduced-Order Models
Optimal control of the cylinder wake in the laminar regime by Trust-Region methods and POD Reduced-Order Models M Bergmann a,, L Cordier b a MAB, UMR - Université Bordeaux, 3 cours de la Libération, 33
More informationA-posteriori error estimation of discrete POD models for PDE-constrained optimal control
Wegelerstraße 6 53115 Bonn Germany phone +49 228 73-3427 fax +49 228 73-7527 www.ins.uni-bonn.de M. Gubisch, I. Neitzel, S. Volkwein A-posteriori error estimation of discrete POD models for PDE-constrained
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 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 informationHybrid Model Reduction for Compressible Flow Controller Design
Hybrid Model Reduction for Compressible Flow Controller Design Xiaoqing Ge* and John T. Wen Department of Electrical, Computer and Systems Engineering Rensselaer Polytechnic Institute, Troy, NY, USA Abstract
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
More informationAlgorithms for PDE-Constrained Optimization
GAMM-Mitteilungen, 31 January 2014 Algorithms for PDE-Constrained Optimization Roland Herzog 1 and Karl Kunisch 2 1 Chemnitz University of Technology, Faculty of Mathematics, Reichenhainer Straße 41, D
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
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 informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
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 informationOptimal rotary control of the cylinder wake using POD reduced order model p. 1/30
Optimal rotary control of the cylinder wake using POD reduced order model Michel Bergmann, Laurent Cordier & Jean-Pierre Brancher Michel.Bergmann@ensem.inpl-nancy.fr Laboratoire d Énergétique et de Mécanique
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 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 informationRegularization and Inverse Problems
Regularization and Inverse Problems Caroline Sieger Host Institution: Universität Bremen Home Institution: Clemson University August 5, 2009 Caroline Sieger (Bremen and Clemson) Regularization and Inverse
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 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 informationOptimal control of partial differential equations with affine control constraints
Control and Cybernetics vol. 38 (2009) No. 4A Optimal control of partial differential equations with affine control constraints by Juan Carlos De Los Reyes 1 and Karl Kunisch 2 1 Departmento de Matemática,
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 informationConvergence of a finite element approximation to a state constrained elliptic control problem
Als Manuskript gedruckt Technische Universität Dresden Herausgeber: Der Rektor Convergence of a finite element approximation to a state constrained elliptic control problem Klaus Deckelnick & Michael Hinze
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 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 informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationSome recent results for two extended Navier-Stokes systems
Some recent results for two extended Navier-Stokes systems Jim Kelliher UC Riverside Joint work with Mihaela Ignatova (UCR), Gautam Iyer (Carnegie Mellon), Bob Pego (Carnegie Mellon), and Arghir Dani Zarnescu
More informationA Posteriori Error Estimation for Reduced Order Solutions of Parametrized Parabolic Optimal Control Problems
A Posteriori Error Estimation for Reduced Order Solutions of Parametrized Parabolic Optimal Control Problems Mark Kärcher 1 and Martin A. Grepl Bericht r. 361 April 013 Keywords: Optimal control, reduced
More informationError estimates for the discretization of the velocity tracking problem
Numerische Mathematik manuscript No. (will be inserted by the editor) Error estimates for the discretization of the velocity tracking problem Eduardo Casas 1, Konstantinos Chrysafinos 2 1 Departamento
More informationThe Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems
The Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems Mehdi Badra Laboratoire LMA, université de Pau et des Pays
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 informationOptimal control of semilinear elliptic PDEs with state constraints - numerical analysis and implementation
Optimal control of semilinear elliptic PDEs with state constraints - numerical analysis and implementation Dissertation with the aim of achieving a doctoral degree at the Faculty of Mathematics, Informatics
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 informationMaster Thesis. POD-Based Bicriterial Optimal Control of Time-Dependent Convection-Diffusion Equations with Basis Update
Master Thesis POD-Based Bicriterial Optimal Control of Time-Dependent Convection-Diffusion Equations with Basis Update submitted by Eugen Makarov at the Department of Mathematics and Statistics Konstanz,.6.18
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 informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
More informationMODEL ORDER REDUCTION
MODEL ORDER REDUCTION Approximate yet accurate surrogates for large-scale simulation KAREN E. WILLCOX Science at Extreme Scales: Where Big Data Meets Large-Scale Computing Tutorials Institute for Pure
More 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 informationON TWO OPTIMAL CONTROL PROBLEMS FOR MAGNETIC FIELDS
ON TWO OPTIMAL CONTROL PROBLEMS FOR MAGNETIC FIELDS SERGE NICAISE, SIMON STINGELIN, AND FREDI TRÖLTZSCH Abstract. Two optimal control problems for instationary magnetization processes are considered in
More informationPOD A-POSTERIORI ERROR BASED INEXACT SQP METHOD FOR BILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS
ESAIM: M2AN 46 22) 49 5 DOI:.5/m2an/26 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org POD A-POSTERIORI ERROR BASED INEXACT SQP METHOD FOR BILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS
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 informationSemismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations
Semismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations Axel Kröner 1 Karl Kunisch 2 and Boris Vexler 3 1 Lehrstuhl für Mathematische Optimierung Technische Universität München
More informationPOD/DEIM Strategies for reduced data assimilation systems
POD/DEIM Strategies for reduced data assimilation systems Răzvan Ştefănescu 1 Adrian Sandu 1 Ionel M. Navon 2 1 Computational Science Laboratory, Department of Computer Science, Virginia Polytechnic Institute
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationInexact Algorithms in Function Space and Adaptivity
Inexact Algorithms in Function Space and Adaptivity Anton Schiela adaptive discretization nonlinear solver Zuse Institute Berlin linear solver Matheon project A1 ZIB: P. Deuflhard, M. Weiser TUB: F. Tröltzsch,
More informationSnapshot location by error equilibration in proper orthogonal decomposition for linear and semilinear parabolic partial differential equations
Numerical Analysis and Scientific Computing Preprint Seria Snapshot location by error equilibration in proper orthogonal decomposition for linear and semilinear parabolic partial differential equations
More informationHessian-Based Model Reduction with Applications to Initial-Condition Inverse Problems. Omar Shahid Bashir
Hessian-Based Model Reduction with Applications to Initial-Condition Inverse Problems by Omar Shahid Bashir Submitted to the Department of Aeronautics and Astronautics in partial fulfillment of the requirements
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
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 informationReduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD
3rd AIAA Flow Control Conference, 5 8 June 26, San Francisco Reduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD M. Ilak and C. W. Rowley Dept. of Mechanical and Aerospace Engineering,
More informationCIMPA Summer School on Current Research on Finite Element Method Lecture 1. IIT Bombay. Introduction to feedback stabilization
1/51 CIMPA Summer School on Current Research on Finite Element Method Lecture 1 IIT Bombay July 6 July 17, 215 Introduction to feedback stabilization Jean-Pierre Raymond Institut de Mathématiques de Toulouse
More informationFEEDBACK STABILIZATION OF TWO DIMENSIONAL MAGNETOHYDRODYNAMIC EQUATIONS *
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LV, 2009, f.1 FEEDBACK STABILIZATION OF TWO DIMENSIONAL MAGNETOHYDRODYNAMIC EQUATIONS * BY CĂTĂLIN-GEORGE LEFTER Abstract.
More information