Bang bang control of elliptic and parabolic PDEs
|
|
- Catherine Henderson
- 5 years ago
- Views:
Transcription
1 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 2/26 Model problem (P) min J(y, u) = 1 u U ad,y Y ad 2 y z 2 + α D 2 u 2 U subject to y = G(Bu). y = G(Bu) iff Ay = Bu in D, plus b.c. (plus i.c.) Examples Ay = d i,j=1 x j (a ij y xi ) + d i =1 b i y xi + cy uniformly elliptic operator, Y = H 1 (Ω) Ay = y t d i,j=1 x j (a ij y xi ) + d i =1 b i y xi + cy parabolic operator, Y = W (0, T ). Constraint sets: U ad = {a u b} and Y ad = {y c}, or { y δ}. Bang bang control: α = 0.
3 3/26 Existence and uniqueness of solutions, optimality conditions For every u we have a unique y(u) = G(Bu). So we may minimize the reduced functional J(u) = J(G(Bu), u) instead. Problem (P) admits a unique solution u with corresponding state y = G(Bu). J (u), v u U,U 0 for all v F ad = {v U ad, G(Bv) Y ad }.
4 4/26 Necessary optimality conditions (α > 0), Y ad = {y c} Constraint qualification (Slater condition) ũ U ad such that G(Bũ) < c in D Then there exist µ M( D) and some p such that with y and u there holds Ay = Bu A p = y z + µ u = P Uad ( 1 α RB p), µ 0, y c in D and D(c y)dµ = 0, Regularity elliptic case: p W 1,s (Ω) for all s < d/(d 1). parabolic case: p L s (W 1,σ ) for all s, σ [1, 2) with 2/s + d/σ > d + 1. Low regularity of adjoint states introduces difficulties in the numerical analysis
5 5/26 Only state constraints (α > 0), U ad U Optimality conditions: Ay = Bu A p = y z + µ u = 1 α RB p, µ 0, y c in D and y)dµ = 0, D(c The state y in general is smoother than the adjoint p. The optimal control and the adjoint state p are coupled via an algebraic relation a discretization of p induces a discretization of u. A discretization of y ideally should deliver feasible discrete states.
6 6/26 Only control constraints (α > 0), Y ad Y Optimality conditions: Ay = Bu A p = y z u = P Uad ( 1 α RB p), The adjoint p in general is smoother than the state y. The optimal control and the adjoint state p are coupled via an algebraic relation a discretization of p induces a discretization of u. Reduction to p, say delivers AA p BP Uad ( 1 α RB p) = Az, which in the parabolic case is a bvp for p in space-time.
7 7/26 Control constraints, α = 0 The function u U ad is a solution of the optimal control problem iff there exists an adjoint state p such that y = G(u), p = G(y z) and (αu + RB p, v u) 0 for all v U ad. Recall: u = P Uad ( 1 α RB p) for α > 0, i.e. u = a, αu + RB p > 0, 1 α RB p, αu + RB p = 0, b, αu + RB p < 0, If α = 0: u = a, RB p > 0, [a, b], RB p = 0, = b, RB p < 0. Control is determined through the sign of RB p. Control in the generic case is either at upper or lower bound. Use homothopy α 0 to compute bang bang control.
8 8/26 Tailored discretization Discrete counterpart to (P): (P h ) min J h (y h, u h ) s.t. PDE h (y h ) = B h (u h ) u h U h ad,y h Y h ad Questions: Appropriate choice of U h ad and Ansatz for u h? Appropriate choice of Y h ad and Ansatz for y h? Aim: Capture as much structure as possible of (P) on the discrete level.
9 9/26 Discretization a variational concept 1 Discrete optimal control problem: (P h ) min J h (u) = 1 u U ad 2 y h z 2 + α D 2 u 2 U subject to y h = G h (Bu) and y h I h c. Here, y h (u) = G h (Bu) denotes e.g. the This problem is still dimensional. p.l. and continuous fe approximation to y(u) (elliptic case), dg(0) in time and p.l. and continuous fe in space approximation to y(u) (parabolic case), i.e. a(y h, v h ) = Bu, v h for all v h X h. The control is not discretized explicitely! 1 M. Hinze: A variational discretization concept in control constrained optimization: the linear-quadratic case. J. Comput. Optim. Appl. 30, (2005)
10 10/26 Discrete necessary optimality conditions (α > 0) Problem (P h ) admits a unique solution u h U ad. Furthermore, uniform convergence of discrete states implies discrete Slater condition G h (Bũ) < I h b in D for h small enough Then there exist µ h M( D) and some p h X h such that with y h and u h there holds a(y h, v h ) = Bu h, v h v h X h, a(w h, p h ) = D (y h z)w h + D w h dµ h w h Y h, u h = P Uad ( 1 α RB p h ) µ h 0, y h I h c, and D(I h c y h )dµ h = 0.
11 11/26 Discrete necessary optimality conditions (α > 0) Problem (P h ) admits a unique solution u h U ad. Furthermore, uniform convergence of discrete states implies discrete Slater condition G h (Bũ) < I h b in D for h small enough Then there exist µ h M( D) and some p h X h such that with y h and u h there holds a(y h, v h ) = Bu h, v h v h X h, a(w h, p h ) = D (y h z)w h + D w h dµ h w h Y h, u h = P Uad ( 1 α RB p h ) µ h 0, y h I h c, and D(I h c y h )dµ h = 0.
12 12/26 Discrete necessary optimality conditions (α > 0) Problem (P h ) admits a unique solution u h U ad. Furthermore, uniform convergence of discrete states implies discrete Slater condition G h (Bũ) < I h b in D for h small enough Then there exist µ h M( D) and some p h X h such that with y h and u h there holds a(y h, v h ) = Bu h, v h v h X h, a(w h, p h ) = D (y h z)w h + D w h dµ h w h Y h, u h = P Uad ( 1 α RB p h ) µ h 0, y h I h c, and D(I h c y h )dµ h = 0. Post processing 2 exploits the projection formula with a discrete adjoint state associated to a fully discrete control. 2 C. Meyer, A. Rösch. Superconvergence properties of optimal control problems. SIAM J. Control Optim. 43: (2004)
13 13/26 α = 0: optimality conditions for discrete problem The variational-discrete optimal control problems admits a solution u h U ad, which is unique in the case α > 0. The state y h is unique (also in the case α = 0). Let u h U ad be a solution of the optimal control problem. Then there exists a unique adjoint state p h such that y h = G h (u h ), p h = G h (y h z) and (αu h + RB p h, v u h ) 0 for all v U ad. Recall u h = P Uad ( 1 α RB p h ) for α > 0, i.e. u h = a, αu h + RB p h > 0, 1 α RB p h, αu h + RB p h = 0, b, αu h + RB p h < 0. If α = 0: u h = a, RB p h > 0, [a, b] RB p h = 0, = b, RB p h < 0.
14 Sketch of concept in 1d Bang bang control problems 14/26
15 15/26 Variational versus conventional discretization (α > 0)
16 16/26 Variational discretization in the bang bang case α = 0 u(x) = sign(p(x)), p(x) = sin(8πx 1)sin(8πx 2 )
17 17/26 Parabolic bang bang with α homothopy p(t, x) = T 2πa 2πat sin ( T ) sin 2πx 1 sin 2πx 2, a = 2, a 1 = 0.2, b 1 = 0.4 bounds. where B p(t) = sin 2πx 1 sin 2πx 2 p(t, x)dx. Ω u(t) = { a 1, B p > 0, b 1, B p < 0, Here, α = h 2 = k 4/3 with h = 2 l. Furthermore, with this coupling 3 u u k,h,α L 1 ch 2. 3 N. von Daniels, M. Hinze: Variational discretization of a control-constrained parabolic bang-bang optimal control problem. J. Comput. Math., to appear
18 18/26 Algorithms for P α h Define G h (u) = u P Uad ( 1 α p h(y h (u))). The optimality condition reads G h (u) = 0 and motivates the fix point iteration u given, do until convergence u + = P Uad ( 1 α p h(y h (u))), u = u Is this algorithm numerically implementable? Yes, whenever for given u it is possible to numerically evaluate the expression P Uad ( 1 α p h(y h (u))) in the i th iteration, with an numerical overhead which is independent of the iteration counter of the algorithm.
19 19/26 Semi smooth Newton algorithm for α > 0 2. Does the fix point algorithm converge? Yes, if α > RB S h S hb L(U), since P Uad is non expansive. Condition too restrictive for our purpose semi smooth Newton method applied to G h (u) = 0: u given, solve until convergence G h (u)u+ = G h (u) + G h (u)u, u = u+. 1. This algorithm is implementable whenever the fix point iteration is, since G h (u) + G h (u)u = P U ad ( 1 α p h(u)) 1 α P U ad ( 1 α p h(u)) S h S hu. 2. For every α > 0 this algorithm is locally fast convergent 4. 4 M. Hinze, M. Vierling: Variational discretization and semi-smooth Newton methods; implementation, convergence and globalization in pde constrained optimzation with control constraints. Optim. Meth. Software 27: (2012)
20 20/26 Error analysis, sketch for elliptic case 5 It is well known that y y h + α u u h y y h (u) + p p h (y(u)) So one expects estimates for y y h also in the case α = 0. Estimates for u u h? Estimate for the states (S = {x Ω p(x) 0} Ω) y y h C (h 2 + (b a) p R h p L 1 (Ω S) + p R h p L u u h L 1 (S) ), p p h L C y y h + p R h p L, follow from 0 (p p h, u h u) = (R h p p h, u h u) + (p R h p, u h u) I + II. I 1 2 y y h y R hy 2 II = Ω S (p R h p)(u h u) + S (p R h p)(u h u). 5 K. Deckelnick, M. Hinze: A note on the approximation of elliptic control problems with bang-bang controls. Comput. Optim. Appl. 51: (2012)
21 21/26 Error estimates If the adjoint state p in the solution with some β (0, 1] satisfies the structural assumption C > 0 ɛ > 0 L({x Ω; p(x) ɛ}) Cɛ β, one for h small enough gets a unique variational discrete control u h, which together with the discrete state y h and the discrete adjoint p h satisfies y y h + p p h L C 1 h2 2 β + p R h p L, u u h L 1 C β h2β 2 β + p R h p L.
22 22/26 Sketch of proof for β = 1 u u h L 1, y y h, p p h L C {h 2 + p R h p L } Sketch of proof: 0 (p p h, u h u) = (R h p p h, u h u) + (p R h p, u h u) I + II. I 1 2 y y h y R hy 2 II = S (p R h p)(u h u). Combine now u u h L 1 (b a)l({p > 0, p h 0} {p < 0, p h 0}) {p > 0, p h 0} {p < 0, p h 0} { p(x) p p h } L({ p(x) p p h }) C p p h u u h L 1 C p p h p p h p R h p + R h p p h R h p p h C y y h to estimate II.
23 23/26 Numerical example with 2 switching points, fix-point iteration u(x) = sign( sin(mπx)), y(x) = sin(mπx), and p(x) = sin(mπx), m= exact discrete FE grid 0.2 exact discrete FE grid Experimental order of convergence: u u h L 1 : Function values p p h L : y y h L : p p h L 2 :
24 24/26 Special cases 1. If the desired state z is reachable with a feasible control, then y y h + p p h L Ch If p C 1 ( Ω) satisfies min p(x) > 0, where K = {x Ω p(x) = 0}. x K Then, the structural assumption is satisfied with β = If p W 2, (Ω) and satisfies the structural assumption, then y y h + p p h L + u u h L 1 Ch 2 log h γ(d).
25 25/26 Further reading M. Hinze: A variational discretization concept in control constrained optimization: the linear-quadratic case. J. Comput. Optim. Appl. 30:45-63 (2005). C. Meyer, A. Rösch. Superconvergence properties of optimal control problems. SIAM J. Control Optim. 43: (2004). K. Deckelnick, M. Hinze: A note on the approximation of elliptic control problems with bang-bang controls. Comput. Optim. Appl. 51: (2012). E. Casas, D. Wachsmuth, and G. Wachsmuth. Sufficient Second-Order Conditions for Bang-Bang Control Problems. SIAM J. Control Optim. 55: (2017). D. Wachsmuth and G. Wachsmuth. Regularization error estimates and discrepancy principle for optimal control problems with inequality constraints. Control and Cybernetics 40: (2011). N. von Daniels, M. Hinze: Variational discretization of a control-constrained parabolic bang-bang optimal control problem. J. Comput. Math., accepted for publication. N. von Daniels. Tikhonov regularization of control-constrained optimal control problems. Comput. Optim. Appl. 70: (2018). U. Felgenhauer. On Stability of Bang-Bang Type Controls. SIAM J. Control Optim. 41: (2003). M. Seydenschwanz. Convergence results for the discrete regularization of linear-quadratic control problems with bang-bang solutions. Comput. Optim. Appl. 61: (2015). D. Wachsmuth. Adaptive regularization and discretization of bang-bang optimal control problems. ETNA 40: (2013). D. Wachsmuth. Robust error estimates for regularization and discretization of bang-bang control problems. Comp. Opt. Appl. 62: (2014).
26 26/26 End Thank you very much for your attention!
Robust 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 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 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 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 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 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 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 informationNecessary conditions for convergence rates of regularizations of optimal control problems
Necessary conditions for convergence rates of regularizations of optimal control problems Daniel Wachsmuth and Gerd Wachsmuth Johann Radon Institute for Computational and Applied Mathematics RICAM), Austrian
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 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 informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
More informationNumerical Solution of Elliptic Optimal Control Problems
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,
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 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 informationPriority Program 1253
Deutsce Forscungsgemeinscaft Priority Program 53 Optimization wit Partial Differential Equations K. Deckelnick, A. Günter and M. Hinze Finite element approximation of elliptic control problems wit constraints
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
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 informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationMidterm Solution
18303 Midterm Solution Problem 1: SLP with mixed boundary conditions Consider the following regular) Sturm-Liouville eigenvalue problem consisting in finding scalars λ and functions v : [0, b] R b > 0),
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
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 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 informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
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 informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
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 informationON CONVERGENCE OF A RECEDING HORIZON METHOD FOR PARABOLIC BOUNDARY CONTROL. fredi tröltzsch and daniel wachsmuth 1
ON CONVERGENCE OF A RECEDING HORIZON METHOD FOR PARABOLIC BOUNDARY CONTROL fredi tröltzsch and daniel wachsmuth Abstract. A method of receding horizon type is considered for a simplified linear-quadratic
More informationAn Introduction to Optimal Control of Partial Differential Equations with Real-life Applications
An Introduction to Optimal Control of Partial Differential Equations with Real-life Applications Hans Josef Pesch Chair of Mathematics in Engineering Sciences University of Bayreuth, Bayreuth, Germany
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
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 informationModel order reduction & PDE constrained optimization
1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell,
More 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 informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationA Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems
A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems D. Meidner and B. Vexler Abstract In this article we discuss a priori error estimates for Galerkin
More informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
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 informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More informationG: Uniform Convergence of Fourier Series
G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) u t = Du xx 0 < x < l, t > 0 u(0, t) = 0 = u(l, t) t > 0 u(x, 0) = f(x) 0 < x < l () we developed
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
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 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 informationSolving the Heat Equation (Sect. 10.5).
Solving the Heat Equation Sect. 1.5. Review: The Stationary Heat Equation. The Heat Equation. The Initial-Boundary Value Problem. The separation of variables method. An example of separation of variables.
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 informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
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 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 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 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 informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More 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 informationSimulation based optimization
SimBOpt p.1/52 Simulation based optimization Feb 2005 Eldad Haber haber@mathcs.emory.edu Emory University SimBOpt p.2/52 Outline Introduction A few words about discretization The unconstrained framework
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationOPTIMALITY CONDITIONS FOR STATE-CONSTRAINED PDE CONTROL PROBLEMS WITH TIME-DEPENDENT CONTROLS
OPTIMALITY CONDITIONS FOR STATE-CONSTRAINED PDE CONTROL PROBLEMS WITH TIME-DEPENDENT CONTROLS J.C. DE LOS REYES P. MERINO J. REHBERG F. TRÖLTZSCH Abstract. The paper deals with optimal control problems
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
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 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 informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationA PRIORI ERROR ESTIMATE AND SUPERCONVERGENCE ANALYSIS FOR AN OPTIMAL CONTROL PROBLEM OF BILINEAR TYPE *
Journal of Computational Mathematics, Vol.??, No.??, 2008, 1 17. A PRIORI ERROR ESTIMATE AND SUPERCONVERGENCE ANALYSIS FOR AN OPTIMAL CONTROL PROBLEM OF BILINEAR TYPE * Danping Yang Department of Mathematics,
More informationCoordinate Update Algorithm Short Course Operator Splitting
Coordinate Update Algorithm Short Course Operator Splitting Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 25 Operator splitting pipeline 1. Formulate a problem as 0 A(x) + B(x) with monotone operators
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 informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationAnalysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017
Analysis III (BAUG Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017 Question 1 et a 0,..., a n be constants. Consider the function. Show that a 0 = 1 0 φ(xdx. φ(x = a 0 + Since the integral
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
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 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 informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationMaximum principle for the fractional diusion equations and its applications
Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationPARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr
More informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
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 informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationThe Normal Equations. For A R m n with m > n, A T A is singular if and only if A is rank-deficient. 1 Proof:
Applied Math 205 Homework 1 now posted. Due 5 PM on September 26. Last time: piecewise polynomial interpolation, least-squares fitting Today: least-squares, nonlinear least-squares The Normal Equations
More informationMATH 4211/6211 Optimization Basics of Optimization Problems
MATH 4211/6211 Optimization Basics of Optimization Problems Xiaojing Ye Department of Mathematics & Statistics Georgia State University Xiaojing Ye, Math & Stat, Georgia State University 0 A standard minimization
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
More 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
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 informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationAdaptive 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 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 informationA Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow
A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,
More informationL -Estimates of Rectangular Mixed Methods for Nonlinear Constrained Optimal Control Problem
BULLETIN of the MALAYSIAN MATHEMATICAL SCIENCES SOCIETY http:/math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 37(1) (2014), 271 284 L -Estimates of Rectangular Mixed Methods for Nonlinear Constrained
More informationMultigrid methods for pde optimization
Multigrid methods for pde optimization Alfio Borzì Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Alfio Borzì (Università degli Studi del Sannio, Multigrid Italy)
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More information