Numerical Solution of Elliptic Optimal Control Problems
|
|
- Lucinda Haynes
- 5 years ago
- Views:
Transcription
1 Department of Mathematics, University of Houston Numerical Solution of Elliptic Optimal Control Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg Summer School Optimal Control of PDEs Cortona, July 12-17, 2010
2 Department of Mathematics, University of Houston Recent Books on Optimal Control of PDEs R. Glowinski, J.L. Lions, and J. He; Exact and Approximate Controllability for Distributed Parameter Systems: A Numerical Approach. Cambridge University Press, Cambridge, M. Hinze, R. Pinnau, and M. Ulbrich; Optimization with PDE Constraints. Springer, Berlin-Heidelberg-New York, F. Tröltzsch; Optimal Control of Partial Differential Equations. Theory, Methods, and Applications. American Mathematical Society, Providence, 2010.
3 Numerical Solution of Elliptic Optimal Control Problems Unconstrained Elliptic Optimal Control Problems Control Constrained Elliptic Optimal Control Problems State Constrained Elliptic Optimal Control Problems
4 Department of Mathematics, University of Houston Elliptic Optimal Control Problems Unconstrained Case
5 Elliptic Optimal Control Problems: Unconstrained Case Consider the unconstrained distributed elliptic optimal control problem inf J(y,u) := 1 y y d 2 dx + α u 2 dx, y,u 2 2 where u L 2 (Ω). Ω (EUC) 1 y = u in Ω, (EUC) 2 y = 0 on Γ, Ω
6 Existence and Uniqueness of a Solution We introduce the control-to-state map G : L 2 (Ω) H 1 0(Ω), y = Gu solves y = u. Substituting y = Gu allows to consider the reduced control problem inf u J red(u) := 1 2 Gu yd 2 L 2 (Ω) + α 2 u 2 L 2 (Ω). Theorem (Existence and uniqueness) The unconstrained elliptic control problem admits a unique solution y H 1 0(Ω). Proof. Direct method of the calculus of variations (minimizing sequence).
7 Optimality Conditions for the Distributed Control Problem There exists an adjoint state p H 1 0(Ω) such that the triple (y,p,u) satisfies a(y,v) = (u,v) L 2 (Ω), v H 1 0(Ω), a(p,v) = (y y d,v) L 2 (Ω), v H 1 0(Ω), p αu = 0, where a(, ) stands for the bilinear form a(w,z) = Ω w z dx, w,z H 1 0(Ω).
8 Optimality Conditions for the FE Discretized Control Problem There exists an adjoint state p h V h such that the triple (y h,p h,u h ) satisfies a(y h,v h ) = (u h,v h ) L 2 (Ω), v h V h, a(p h,v h ) = (y h y d,v h ) L 2 (Ω), v h V h, p h αu h = 0.
9 Finite Element Approximation of the Distributed Control Problem Let T h (Ω) be a shape regular, simplicial triangulation of Ω and let V h := { v h C(Ω) v h T P 1 (T), T T h (Ω), v h Ω = 0 } be the FE space of continuous, piecewise linear finite elements. Consider the following FE Approximation of the distributed control problem Minimize J(y h,u h ) := 1 2 y h y d 2 L 2 (Ω) + α 2 u h 2 L 2 (Ω), over (y h,u h ) V h V h, subject to a(y h,v h ) = (u h,v h ) L 2 (Ω), v h V h.
10 Algebraic Formulation of the Discrete Optimality System Let V l := span{ϕ l (1),,ϕ(N l ) l }, where ϕ (i) l,1 i N l, are the nodal basis functions. We refer to M l lr N l N l and A l lr N l N l as the associated mass matrix and stiffness matrix: (M l ) ij := (ϕ (i) l,ϕ(j) l ) L 2 (Ω), (A l ) ij := a(ϕ (i) l, ϕ(j) l ), 1 i,j N l. The solution of the discrete optimality system requires the computation of (y l,p l ) lr N l lr N l as the solution of A l α 1 M l M l A l y l p l = 0 M l y d l.
11 Left and Right Transforming Iterations For given K lr N N,d lr N, we want to solve the linear algebraic system ( ) Kx = d. Let K L,K R lr N N be regular. Then, ( ) can be equivalently wriiten as ( ) K L KK R K 1 R x = K L d. Assuming K to be a suitable preconditioner for K L KK R, we consider the transforming iteration K 1 R x (ν+1) = K 1 R x (ν) + K 1 (K L d K L Kx (ν) ). Backtransformation yields so that ˆK := K 1 L x (ν+1) = x (ν) + K R K 1 K }{{ L (d Kx (ν) ), } = (K 1 L KK 1 R ) 1 1 KK R is an appropriate preconditioner for the original system.
12 Department of Mathematics, University of Houston Left and Right Transforming Iterations II We apply the idea of left/right transformations to K = A l α 1 M l M l A l, d = 0 M l y d l We choose the left and right transformations according to K L = α 1/2 I 0 0 I, K R = α 1/2 I 0 0 I.. This leads to the transformed matrix K L KK R = A l α 1/2 M l α 1/2 M l A l.
13 Left and Right Transforming Iterations III Denoting by S l the Schur complement S l = A l + α 1 M l A 1 l M l, we have A l α 1/2 M l A α 1/2 = l α 1/2 M l M l A l α 1/2 M l S l + α 1 M l A 1 l M. l Choosing  l as a preconditioner for A l and Ŝ l := τ 1 l diag(a l + α 1 M l  1 l M l ), we obtain the symmetric Uzawa preconditioner K l =  l α 1/2 M l α 1/2 M l Ŝ l + α 1 M l  1 l M l Backtransformation results in the following preconditioner ˆK l = K 1 L K l K 1 R =.  l α 1 M l M l Ŝ l α 1 M l  l 1 M l.
14 Multigrid Solvers for Elliptic Optimal Control Problems A. Borzi, K. Kunisch, and D. Y. Kwak; Accuracy and convergence properties of the finite difference multigrid solution of an optimal control optimality system. SIAM J. Control Optimization 41, , A. Borzi and V. Schulz; Multigrid methods for PDE optimization. SIAM Rev. 51, , J. Schöberl, R. Simon, and W. Zulehner; A robust multigrid method for elliptic optimal control problems. Preprint, Inst. of Comput. Math., University of Linz, 2010.
15 Department of Mathematics, University of Houston Step 1: Given iterates y (ν) l,p(ν) l Preconditioned Richardson Iteration I Res (1) l Res (2) l lr N l, compute the residuals = 0 M l y d l Step 2: Solve the linear system  l α 1 M l M l Ŝ l α 1 M l  1 l M l Step 3: Compute the new iterates y (ν+1) l p (ν+1) l = y (ν) l p (ν) l K l y (ν) l p (ν) l + y (ν) l p (ν) l = y (ν) l p (ν) l.. Res (1) l Res (2) l.
16 Department of Mathematics, University of Houston Consider the linear system ( ) Preconditioned Richardson Iteration II  l α 1 M l M l Ŝ l α 1 M l  1 l M l y (ν) l p (ν) l = Res (1) l Res (2) l. Elimination of y l (ν) results in the Schur complement system Ŝ l p l (ν) = Res l (2) Ml 1 l Res(1) l. Hence, the solution of ( ) can be reduced to the successive solution of the three linear systems  l y l (ν) = Res l (1), Ŝ l p (ν) l = Res l (2) M (ν) l y l,  l y (ν) l = Res (1) l + α 1 M l p (ν) l.
17 Preconditioned Richardson Iteration III Theorem [Schöberl/Simon/Zulehner] If the preconditioners  l and Ŝ l are chosen according to then there holds  l = σ 1 l diag(a l ), Ŝ l = τ 1 l diag(s l ), σ l λ max ((diag(a l )) 1 A l ) 1, τ l λ max ((diag(s l )) 1 S l ) 1, where (y l (ν),p(ν) l ) (y l,p l ) C (1 + α 1/2 h 2 l ) η(ν) (y(0) l,p(0) l ) (y l,p l ), η(ν) 2 π(ν 1) 2 πν, if ν is even,, if ν is odd and C is a positive constant independent of h l and α.
18 Multigrid Solution of the Optimality System The previous theorem suggests to use the preconditioned iterative solver as a smoother within a multigrid solution of the optimality system. Theorem. In case of two nested grids T l 1 (Ω) and T l (Ω) let (y l (0),p(0) l ) be a level l startiterate and let (y l (1),p(1) l ) be the result of a two-grid cycle with ν > 0 smoothing steps. Then there holds (y l (1),p(1) l ) (y l,p l ) C 1 + α 1/2 h 2 l (y (0) ν l,p(0) l ) (y l,p l ), where C > 0 is a constant independent of h l and α.
19 Department of Mathematics, University of Houston Elliptic Optimal Control Problems Control Constraints
20 Elliptic Optimal Control Problems: Control Constraints Consider the control constrained distributed elliptic optimal control problem inf J(y,u) := 1 y y d 2 dx + α u 2 dx, y,u 2 2 and where u L 2 (Ω). Ω (ECC) 1 y = u in Ω, (ECC) 2 y = 0 on Γ, u K C := {v L 2 (Ω) v(x) ψ(x) f.a.a. x Ω}, Ω
21 Control Constraints: Necessary Optimality Conditions Theorem. If (y,u) H 1 0(Ω) K C is the optimal solution of the optimal control problem, then there exists (p,λ) H 1 0(Ω)) L 2 (Ω) such that p satisfies the adjoint state equation and there holds p = (y y d ) in Ω, p = 0 on Γ, p λ = αu in Ω, λ I Kc (u).
22 Department of Mathematics, University of Houston Elliptic Optimal Control Problems Primal-Dual Active Set Strategy
23 Moreau-Yosida Approximation of Multivalued Maps I Weighted Duality Mapping: Assume that V is a Banach space with dual V and let h : R + R + be a continuous and non-decreasing function such that h(0) = 0 and h(t) as t. Then the mapping J h : V 2 V J h (u) := {u V u,u = u u, u = h( u )} is called the duality mapping with weight h. Example: For V = L p (Ω),V = L q (Ω),1 < p,q <,1/p +1/q = 1, and h(t) = t p 1 we have J h (u)(x) = { u(x) p 1 sgn(u(x)), u(x) 0 0, u(x) = 0.
24 Moreau-Yosida Approximation of Multivalued Maps II Moreau-Yosida proximal map: Let f : V R be a lower semi-continuous proper convex function with subdifferential f. For c > 0, the Moreau-Yosida proximal map P f c : V 2 V is defined such that P f c (w), w V, is the set of minimizers of where j h = J h. inf f(v) + cj h( v w ), v V c Moreau-Yosida approximation: If J h is single-valued, then for c > 0 the Moreau- Yosida approximation ( f) c of f is given by ( f) c (w) := J h (c 1 w c 1 P f c (w)).
25 Department of Mathematics, University of Houston Moreau-Yosida Approximation of I KC Idea: Approximate I KC by its Moreau-Yosida approximation ( I KC ) c. Theorem. For any c > 0, we have λ ( I KC ) c, if and only if there holds ( ) λ = c u + c 1 λ Π KC (u + c 1 λ) = c max(0,u + c 1 λ ψ), and this is equivalent to u = Π KC (u + c 1 λ), where Π KC denotes the L 2 -projection onto K C.
26 Step 1 (Initialization): Primal-Dual Active Set Strategy I Choose c > 0, start-iterates y (0) h,u(0) h,λ(0) h and set n = 1. Step 2 (Specification of active/inactive sets): Compute the active/inactive sets A n and I n according to A n := {1 i N (u (n 1) h Step 3 (Termination criterion): + c 1 λ (n 1) h ) i > (ψ h ) i }, I n := {1,,N} \ A n. If n 2 and A n = A n 1, stop the algorithm. Otherwise, go to Step 4.
27 Department of Mathematics, University of Houston Primal-Dual Active Set Strategy II Step 4 (Update of the state, adjoint state, and control): Compute y (n) h,p(n) h and set (A h y (n) h ) i = as the solution of { (ψ h ) i, if i A n α 1 (M 1 h p, A h) i, if i I h p (n) h n { (u (n) h ) i := Step 5 (Update of the multiplier): Set λ (n) h (ψ h ) i, if i A n α 1 (M 1 h p(n) h ). i, if i I n := p(n) h αm hu (n) h, n := n + 1, and go to Step 2. = M h(y h y d h),
28 PDAS As A Semi-Smooth Newton Method I Newton differentiability: A function f : S X Y (X, Y Banach spaces) is called Newton differentiable at x S, if there exist a neighborhood U(x) S and a family of mappings g : U(x) L(X, Y) such that f(x + h) f(x) g(x + h)(h) Y = o( h X ) g(x) is said to be the generalized derivative of f at x. (as h X ). Example: The max-function max : L p (Ω) L q (Ω) is Newton differentiable for 1 q < p <. If f : L r (Ω) L p (Ω) is Fréchet-differentiable for some 1 < r, then the function x χ A (x) f(u(x)) is the generalized derivative of max(0,f( )). Here, χ A denotes the characteristic function of the set A where f(u( )) is nonnegative.
29 PDAS As A Semi-Smooth Newton Method II Setting z h := (y h,u h,λ h ), the discrete optimality conditions require the solution of the nonlinear system A h y h M h u h 0 F h (z h ) := αa h u h + M h y h yh d + λ h = 0. λ h + c max(0,u h + c 1 λ h ψ h ) 0 The primal-dual active set strategy (PDAS) is equivalent to the semi-smooth Newton method F h(z (n) h ) z(n) h z (n+1) h = F h (z (n) h ), = z (n) h + z(n) h, where F h (z(n) h ) is the generalized derivative of F h at z (n) h.
30 Department of Mathematics, University of Houston Elliptic Optimal Control Problems Interior Point Methods
31 Interior-Point Methods (General) Gonzaga [1992], Wright [1992], Zhang/Tapia/Dennis [1992] Conn/Gould/Toint [1996], Heinkenschloss [2001] Forsgren/Gill/Wright [2002], M. Ulbrich [2002], S. Ulbrich [2003] Interior-Point Methods (PDE Constrained Optimization) Ulbrich/Ulbrich/Heinkenschloss [1999], Ulbrich/Ulbrich [2000] H./Petrova/Schulz [2002], H./Petrova [2004], Weiser [2005] Ulbrich/Ulbrich [2007], Antil/H./Linsenmann [2007]
32 Interior Point Method I Idea: Coupling of control constraints u h ψ h by, e.g., logarithmic barrier functions with barrier parameter β > 0: N inf (y h,u h ) := J h (y h,u h ) + β ln((ψ h u h ) i ), B (β) y h,u h h which represents a parameter-dependent family of minimization subproblems. Coupling of the discrete state equation by a Lagrange multiplier: i=1 inf sup y h,u h L (β) h p h (y h,u h,p h ) := B (β) h (y h,u h ) + p h,a h y h M h u h.
33 Interior Point Method II The barrier path β x h (β) := (y h (β),u h (β),p h (β)) is the solution of the parameterdependent nonlinear system: A h p h + M h y h y d h A h y h M h u 0 h β F h (x h (β),β) = α(m h u h ) (ψ h u h ) 1 = 0. β α(m h u h ) N + 0 (ψ h u h ) N Solution by path-following predictor-corrector continuation methods in the barrier parameter.
34 Continuation Strategies in Interior-Point Methods (i) Static Strategies Fiacco/McCormick [1990], Byrd/Hribar/Nocedal [1999] Wächter/Biegler [2005] (ii) Dynamic Update Strategies El-Bakry/Tapia/Tsuchiya/Zhang [1996], Gay/Overton/Wright [1998] Vanderbei/Shanno [1999], H./Petrova/Schulz [2002] Tits/Wächter/Bakhtiari/Urban/Lawrence [2003], H./Petrova [2004] Ulbrich/Ulbrich/Vicente [2004], H./Linsenmann/Petrova [2006] Nocedal/Wächter/Waltz [2006], Armand/Benoist/Orban [2007]
35 Department of Mathematics, University of Houston The Barrier Path x( 0 ) Central Path x* "contraction tube"
36 Department of Mathematics, University of Houston Barrier Path: The Predicted Continuation Step x??? x k-1 k x( ) k-1 x( 0 ) x* may leave the Kantorovich neighborhood.
37 Department of Mathematics, University of Houston Barrier Path: The Newton Correction new x k-1 k x x( k ) x( k-1) x( 0 ) x* reduces the continuation steplength adaptively.
38 Numerical Results: Distributed Control Problem with Control Constraints I Minimize J(y,u) := 1 2 y yd 2 0,Ω + α 2 u 2 0,Ω over (y,u) H 1 0(Ω) L 2 (Ω) subject to y = u in Ω, u K := {v L 2 (Ω) v ψ a.e. in Ω} Data: Ω := (0,1) 2, { y d 200 x := 1 x 2 (x 1 0.5) 2 (1 x 2 ), 0 x (x 1 1) (x 2 (x 1 0.5) 2 (1 x 2 ), 0.5 < x 1 1, α = 0.01, ψ = 1.
39 Numerical Results: Distributed Control Problem with Control Constraints I Optimal state (left) and optimal control (right)
40 Numerical Results: Distributed Control Problem with Control Constraints I Grid after 6 (left) and 10 (right) refinement steps
41 Numerical Results: Distributed Control Problem with Control Constraints II Minimize J(y,u) := 1 2 y yd 2 0,Ω + α 2 u ud 2 0,Ω over (y,u) H 1 0(Ω) L 2 (Ω) subject to y = f + u in Ω, u K := {v L 2 (Ω) v ψ a.e. in Ω} Data: Ω := (0,1) 2, y d := 0, u d := û + α 1 (ˆσ 2 û), { (x1 0.5) ψ := 8, (x 1,x 2 ) Ω 1, (x 1 0.5) 2, α := 0.1, f := 0, otherwise { { ψ, (x û := 1,x 2 ) Ω 1 Ω 2, 2.25 (x1 0.75) 10, ˆσ := 4, (x 1,x 2 ) Ω 2, 1.01 ψ, otherwise 0, otherwise ( ) 1/2 Ω 1 := {(x 1,x 2 ) Ω (x 1 0.5) 2 + (x 2 0.5) }, Ω2 := {(x 1,x 2 ) Ω x }.,
42 Numerical Results: Distributed Control Problem with Control Constraints II Optimal state (left) and optimal control (right)
43 Numerical Results: Distributed Control Problem with Control Constraints II Grid after 6 (left) and 10 (right) refinement steps
44 Department of Mathematics, University of Houston Elliptic Optimal Control Problems State Constraints
45 Literature on State-Constrained Optimal Control Problems M. Bergounioux, K. Ito, and K. Kunisch (1999) M. Bergounioux, M. Haddou, M. Hintermüller, and K. Kunisch (2000) M. Bergounioux and K. Kunisch (2002) E. Casas (1986) E. Casas and M. Mateos (2002) J.-P. Raymond and F. Tröltzsch (2000) K. Deckelnick and M. Hinze (2006) M. Hintermüller and K. Kunisch (2007) K. Kunisch and A. Rösch (2002) C. Meyer and F. Tröltzsch (2006) C. Meyer, U. Prüfert, and F. Tröltzsch (2005) U. Prüfert, F. Tröltzsch, and M. Weiser (2004) H./M. Kieweg (2007) A. Günther, M. Hinze (2007) O. Benedix, B. Vexler (2008) M. Hintermüller/H. (2008) W. Liu, W. Gong and N. Yan (2008)
46 Elliptic Optimal Control Problems: State Constraints Consider the state constrained distributed elliptic optimal control problem inf J(y,u) := 1 y y d 2 dx + α u 2 dx, y,u 2 2 and where u L 2 (Ω). Ω (ESC) 1 y = u in Ω, (ESC) 2 y = 0 on Γ, y K S := {y H 1 0(Ω) y(x) ψ(x) f.a.a. x Ω}, Ω
47 State Constraints: Necessary Optimality Conditions Theorem. The state constrained optimal control problem admits a unique solution (y,u) with y W 1,r 0 (Ω),r > 2, and u L2 (Ω). Theorem. Let us assume that there exists u 0 L 2 (Ω) such that the associated solution y 0 W 1,r 0 (Ω) of the state equation satisfies y 0 int(k S ) (Slater condition). If (y,u) is the unique solution of the state constrained optimal control problem, there exist p W 1,s 0 (Ω),1/r + 1/s = 1, and λ M +(Ω) such that ( p, v) L 2 (Ω) = (y y d,v) L 2 (Ω) + λ,v, v W 1,r 0 (Ω), p = αu, λ,y ψ = 0.
48 Department of Mathematics, University of Houston Finite Element Approximation Let T h (Ω) be a simplicial triangulation of Ω and let V h := { v h C(Ω) v h T P 1 (T), T T h (Ω), v h Γ = 0 } be the FE space of continuous, piecewise linear functions. Let ψ h be the V h -interpoland of ψ. Consider the following FE Approximation of the state constrained control problem Minimize J h (y h,u h ) := 1 2 y h y d 2 0,Ω + α 2 u h 2 0,Ω, over (y h,u h ) V h V h, subject to ( y h, v h ) 0,Ω = (u h,v h ) 0,Ω, v h V h, y h K h,s := {v h V h v h (x) ψ h (x), x Ω}. Since the constraints are point constraints associated with the nodal points, the discrete multipliers are chosen from M h := {λ h M(Ω) λ h = a N h (Ω) κ a δ a, κ a lr}.
49 Step 1 (Initialization): Primal-Dual Active Set Strategy I Choose c > 0, start-iterates y (0) h,u(0) h,λ(0) h and set n = 1. Step 2 (Specification of active/inactive sets): Compute the active/inactive sets A n and I n according to A n := {1 i N (y (n 1) h Step 3 (Termination criterion): + c 1 λ (n 1) h ) i > (ψ h ) i }, I n := {1,,N} \ A n. If n 2 and A n = A n 1, stop the algorithm. Otherwise, go to Step 4.
50 Department of Mathematics, University of Houston Primal-Dual Active Set Strategy II Step 4 (Update of the state, adjoint state, control, and multiplier): Compute (y (n) h,u(n) h,p(n) h, λ(n) h ) as the solution of A h y (n) h A h p (n) h (M hy (n) h Set n := n + 1 and go to Step 2. M hu (n) h = 0, yd h) λ (n) h = 0, p (n) h αm hu (n) h = 0, (y (n) h ) i = (ψ h ) i for i A n, (λ (n) h ) i = 0 for i A n.
51 Lavrentiev Regularization: Mixed Control-State Constraints Introduce a regularization parameter ε > 0 and consider the mixed controlstate constrained optimal control problem Minimize J(y ε,u ε ) := 1 2 yε y d 2 L 2 (Ω) + α 2 uε 2 L 2 (Ω), subject to y ε = u ε in Ω, y ε = 0 on Γ D, εu ε + y ε K := {v L 2 (Ω) v(x) ψ(x), f.a.a. x Ω}. Theorem (Optimality conditions). The optimal solution (y ε,u ε ) V L 2 (Ω) is characterized by the existence of an adjoint state p ε V and a multiplier λ ε L 2 +(Ω) such that ( y ε, v) L 2 (Ω) = (u ε,v) L 2 (Ω), v V, ( p ε, v) L 2 (Ω) = (y ε y d,v) L 2 (Ω) + (λ ε,v) L 2 (Ω), v V, p ε α u ε + ελ ε = 0, (λ ε, εu ε + y ε ψ) L 2 (Ω) = 0.
52 Numerical Results: Distributed Control Problem with State Constraints I Minimize J(y,u) := 1 2 y yd 2 0,Ω + α 2 u ud 2 0,Ω over (y,u) H 1 0(Ω) L 2 (Ω) subject to y = u in Ω, y K := {v H 1 0(Ω) v ψ a.e. in Ω} Data: Ω := ( 2, +2) 2, y d (r) := y(r) + p(r) + σ(r), u d (r) := u(r) + α 1 p(r), ψ := 0, α := 0.1, where y(r),u(r),p(r), σ(r) is the solution of the problem: y(r) := r 4/3 + γ 1 (r), u(r) = y(r), p(r) = γ 2 (r) + r r r2, σ(r) :=, γ 1 := γ 2 := 1, r < (r 0.25) (r 0.25) 4 80(r 0.25) 3 + 1, 0.25 < r < , otherwise { 1, r < , otherwise. { 0.0, r < , otherwise,
53 Numerical Results: Distributed Control Problem with State Constraints I Optimal state (left) and optimal control (right)
54 Department of Mathematics, University of Houston Numerical Results: Distributed Control Problem with State Constraints II Minimize J(y,u) := 1 2 y yd 2 0,Ω + α 2 u ud 2 0,Ω over (y,u) H 1 (Ω) L 2 (Ω) subject to y + cy = u in Ω, y K := {v H 1 (Ω) v ψ a.e. in Ω} Data: Ω = B(0,1) := {x = (x 1,x 2 ) T x x 2 2 < 1}, y d (r) := π 1 4π r π ln(r), u d (r) := π r2 1 ln(r), ψ := 4 + r, α := 1. 2π The solution y(r),u(r),p(r), σ(r) of the problem is given by y(r) 4, u(r) 4, p(r) = 1 4π r2 1 2π ln(r), σ(r) = δ 0.
55 Numerical Results: Distributed Control Problem with State Constraints II Optimal state (left) and optimal control (right)
56 Numerical Results: Distributed Control Problem with State Constraints II Optimal adjoint state (left) and mesh after 16 adaptive loops (right)
Adaptive Finite Element Methods for Elliptic Optimal Control Problems
Adaptive Finite Element Methods for Elliptic Optimal Control Problems Ronald H.W. Hoppe 1,2 1 2 Institute of Mathematics, University of Augsburg Summer School Optimal Control of PDEs Cortona, July 12-17,
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Department of Mathematics, University of Houston Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of
More informationA POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS RONALD H.W. HOPPE AND MICHAEL KIEWEG Abstract. We provide an a posteriori error
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 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 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 informationA-posteriori error estimates for optimal control problems with state and control constraints
www.oeaw.ac.at A-posteriori error estimates for optimal control problems with state and control constraints A. Rösch, D. Wachsmuth RICAM-Report 2010-08 www.ricam.oeaw.ac.at A-POSTERIORI ERROR ESTIMATES
More informationNumerical Analysis of State-Constrained Optimal Control Problems for PDEs
Numerical Analysis of State-Constrained Optimal Control Problems for PDEs Ira Neitzel and Fredi Tröltzsch Abstract. We survey the results of SPP 1253 project Numerical Analysis of State-Constrained Optimal
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 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 informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
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 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 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 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 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 informationADAPTIVE MULTILEVEL CORRECTION METHOD FOR FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OPTIMAL CONTROL PROBLEMS
ADAPTIVE MULTILEVEL CORRECTION METHOD FOR FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OPTIMAL CONTROL PROBLEMS WEI GONG, HEHU XIE, AND NINGNING YAN Abstract: In this paper we propose an adaptive multilevel
More informationOn preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints
On preconditioned Uzawa-type iterations for a saddle point problem with inequality constraints Carsten Gräser and Ralf Kornhuber FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationLectures Notes Algorithms and Preconditioning in PDE-Constrained Optimization. Prof. Dr. R. Herzog
Lectures Notes Algorithms and Preconditioning in PDE-Constrained Optimization Prof. Dr. R. Herzog held in July 2010 at the Summer School on Analysis and Numerics of PDE Constrained Optimization, Lambrecht
More informationREGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS
REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,
More informationPreconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization problems
Zeng et al. Journal of Inequalities and Applications 205 205:355 DOI 0.86/s3660-05-0879-x RESEARCH Open Access Preconditioners for reduced saddle point systems arising in elliptic PDE-constrained optimization
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
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 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 informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
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 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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
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 informationRobust error estimates for regularization and discretization of bang-bang control problems
Robust error estimates for regularization and discretization of bang-bang control problems Daniel Wachsmuth September 2, 205 Abstract We investigate the simultaneous regularization and discretization of
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationOPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL
OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL EDUARDO CASAS, ROLAND HERZOG, AND GERD WACHSMUTH Abstract. Semilinear elliptic optimal control
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 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 informationEXISTENCE OF REGULAR LAGRANGE MULTIPLIERS FOR A NONLINEAR ELLIPTIC OPTIMAL CONTROL PROBLEM WITH POINTWISE CONTROL-STATE CONSTRAINTS
EXISTENCE OF REGULAR LAGRANGE MULTIPLIERS FOR A NONLINEAR ELLIPTIC OPTIMAL CONTROL PROBLEM WITH POINTWISE CONTROL-STATE CONSTRAINTS A. RÖSCH AND F. TRÖLTZSCH Abstract. A class of optimal control problems
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationAdaptive Space-Time Finite Element Approximations of Parabolic Optimal Control Problems. Dissertation
Adaptive Space-Time Finite Element Approximations of Parabolic Optimal Control Problems Dissertation zur Erlangung des akademischen Titels eines Doktors der Naturwissenschaften der Mathematisch-Naturwissenschaftlichen
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
More informationPontryagin s Maximum Principle of Optimal Control Governed by A Convection Diffusion Equation
Communications in Mathematics and Applications Volume 4 (213), Number 2, pp. 119 125 RGN Publications http://www.rgnpublications.com Pontryagin s Maximum Principle of Optimal Control Governed by A Convection
More informationPRECONDITIONED SOLUTION OF STATE GRADIENT CONSTRAINED ELLIPTIC OPTIMAL CONTROL PROBLEMS
PRECONDITIONED SOLUTION OF STATE GRADIENT CONSTRAINED ELLIPTIC OPTIMAL CONTROL PROBLEMS Roland Herzog Susann Mach January 30, 206 Elliptic optimal control problems with pointwise state gradient constraints
More informationA Mesh-Independence Principle for Quadratic Penalties Applied to Semilinear Elliptic Boundary Control
Schedae Informaticae Vol. 21 (2012): 9 26 doi: 10.4467/20838476SI.12.001.0811 A Mesh-Independence Principle for Quadratic Penalties Applied to Semilinear Elliptic Boundary Control Christian Grossmann 1,
More informationAdaptive Discontinuous Galerkin Approximation of Optimal Control Problems Governed by Transient Convection-Diffusion Equations
Max Planck Institute Magdeburg Preprints Hamdullah Yücel Martin Stoll Peter Benner Adaptive Discontinuous Galerkin Approximation of Optimal Control Problems Governed by Transient Convection-Diffusion Equations
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
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 informationApplications of Linear Programming
Applications of Linear Programming lecturer: András London University of Szeged Institute of Informatics Department of Computational Optimization Lecture 9 Non-linear programming In case of LP, the goal
More informationA stabilized finite element method for the convection dominated diffusion optimal control problem
Applicable Analysis An International Journal ISSN: 0003-6811 (Print) 1563-504X (Online) Journal homepage: http://www.tandfonline.com/loi/gapa20 A stabilized finite element method for the convection dominated
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 informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
More informationElliptic optimal control problems with L 1 -control cost and applications for the placement of control devices
Comput Optim Appl (2009) 44: 159 181 DOI 10.1007/s10589-007-9150-9 Elliptic optimal control problems with L 1 -control cost and applications for the placement of control devices Georg Stadler Received:
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 informationThe Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations
The Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations Iain Smears INRIA Paris Linz, November 2016 joint work with Endre Süli, University of Oxford Overview
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 informationA. RÖSCH AND F. TRÖLTZSCH
ON REGULARITY OF SOLUTIONS AND LAGRANGE MULTIPLIERS OF OPTIMAL CONTROL PROBLEMS FOR SEMILINEAR EQUATIONS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS A. RÖSCH AND F. TRÖLTZSCH Abstract. A class of nonlinear
More informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationOptimal Control of Partial Differential Equations I+II
About the lecture: Optimal Control of Partial Differential Equations I+II Prof. Dr. H. J. Pesch Winter semester 2011/12, summer semester 2012 (seminar), winter semester 2012/13 University of Bayreuth Preliminary
More informationSEVERAL PATH-FOLLOWING METHODS FOR A CLASS OF GRADIENT CONSTRAINED VARIATIONAL INEQUALITIES
SEVERAL PATH-FOLLOWING METHODS FOR A CLASS OF GRADIENT CONSTRAINED VARIATIONAL INEQUALITIES M. HINTERMÜLLER AND J. RASCH Abstract. Path-following splitting and semismooth Newton methods for solving a class
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
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 informationParallel-in-Time for Parabolic Optimal Control Problems Using PFASST
Zuse Institute Berlin Takustr. 7 14195 Berlin Germany SEBASTIAN GÖTSCHEL, MICHAEL L. MINION Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST ZIB Report 17-51 (May 2017) Zuse Institute
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationA Control Reduced Primal Interior Point Method for PDE Constrained Optimization 1
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-495 Berlin-Dahlem Germany MARTIN WEISER TOBIAS GÄNZLER ANTON SCHIELA A Control Reduced Primal Interior Point Method for PDE Constrained
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationDesign of optimal RF pulses for NMR as a discrete-valued control problem
Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationConferences Fictitious domain methods for numerical realization of unilateral problems. J.H. (jointly with T. Kozubek and R. Kučera) p.
p. /2 Conferences 28 Fictitious domain methods for numerical realization of unilateral problems J.H. (jointly with T. Kozubek and R. Kučera) Lyon 23-24.6. 28 p. 2/2 Outline Fictitious domain methods -
More informationConvex Optimization. Ofer Meshi. Lecture 6: Lower Bounds Constrained Optimization
Convex Optimization Ofer Meshi Lecture 6: Lower Bounds Constrained Optimization Lower Bounds Some upper bounds: #iter μ 2 M #iter 2 M #iter L L μ 2 Oracle/ops GD κ log 1/ε M x # ε L # x # L # ε # με f
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationNEW REGULARITY RESULTS AND IMPROVED ERROR ESTIMATES FOR OPTIMAL CONTROL PROBLEMS WITH STATE CONSTRAINTS.
ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/ Will be set by the publisher NEW REGULARITY RESULTS AND IMPROVED ERROR ESTIMATES FOR OPTIMAL CONTROL PROBLEMS WITH
More informationJOHANNES KEPLER UNIVERSITY LINZ. A Robust Preconditioner for Distributed Optimal Control for Stokes Flow with Control Constraints
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust Preconditioner for Distributed Optimal Control for Stoes Flow with Control Constraints Marus Kollmann Doctoral Program "Computational
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
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 informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
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 A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationLecture 15 Newton Method and Self-Concordance. October 23, 2008
Newton Method and Self-Concordance October 23, 2008 Outline Lecture 15 Self-concordance Notion Self-concordant Functions Operations Preserving Self-concordance Properties of Self-concordant Functions Implications
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Dietrich Braess, Carsten Carstensen, and Ronald H.W. Hoppe Convergence Analysis of a Conforming Adaptive Finite Element Method for an Obstacle
More informationMatrix square root and interpolation spaces
Matrix square root and interpolation spaces Mario Arioli and Daniel Loghin m.arioli@rl.ac.uk STFC-Rutherford Appleton Laboratory, and University of Birmingham Sparse Days,CERFACS, Toulouse, 2008 p.1/30
More informationSupport Vector Machine via Nonlinear Rescaling Method
Manuscript Click here to download Manuscript: svm-nrm_3.tex Support Vector Machine via Nonlinear Rescaling Method Roman Polyak Department of SEOR and Department of Mathematical Sciences George Mason University
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 informationDual methods for the minimization of the total variation
1 / 30 Dual methods for the minimization of the total variation Rémy Abergel supervisor Lionel Moisan MAP5 - CNRS UMR 8145 Different Learning Seminar, LTCI Thursday 21st April 2016 2 / 30 Plan 1 Introduction
More informationError estimates for the finite-element approximation of a semilinear elliptic control problem
Error estimates for the finite-element approximation of a semilinear elliptic control problem by Eduardo Casas 1 and Fredi röltzsch 2 1 Departamento de Matemática Aplicada y Ciencias de la Computación
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Matematik Micael Hintermüller, Ronald H.W. Hoppe Goal-Oriented Adaptivity in Pointwise State Constrained Optimal Control of Partial Differential Equations
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 informationPreconditioning Iterative Methods for PDE Constrained Optimization
Preconditioning Iterative Methods for PDE Constrained Optimization Tyrone Rees New College University of Oxford A thesis submitted for the degree of Doctor of Philosophy Hilary 2 This thesis is dedicated
More informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationEvolution Under Constraints: Fate of Methane in Subsurface
Evolution Under Constraints: Fate of Methane in Subsurface M. Peszyńska 1 Department of Mathematics, Oregon State University SIAM PDE. Nov.2011 1 Research supported by DOE 98089 Modeling, Analysis, and
More informationVariational Image Restoration
Variational Image Restoration Yuling Jiao yljiaostatistics@znufe.edu.cn School of and Statistics and Mathematics ZNUFE Dec 30, 2014 Outline 1 1 Classical Variational Restoration Models and Algorithms 1.1
More informationA robust optimal preconditioner for the mixed finite element discretization of elliptic optimal control problems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2010; 00:1 19 Published online in Wiley InterScience (www.interscience.wiley.com). A robust optimal preconditioner for the mixed finite
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More information