Martin Gubisch Stefan Volkwein
|
|
- Nicholas Lester
- 5 years ago
- Views:
Transcription
1 Universität Konstanz POD a-posteriori error analysis for optimal control problems with mixed control-state constraints Martin Gubisch Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 34, Juni 3 ISSN Fachbereich Mathematik und Statistik Universität Konstanz Fach D 97, Konstanz, Germany
2
3 Noname manuscript No. (will be inserted by the editor) POD a-posteriori error analysis for optimal control problems with mixed control-state constraints Martin Gubisch Stefan Volkwein Received: April, 3/ Accepted: date Abstract In this work linear-quadratic optimal control problems for parabolic equations with mixed control-state constraints are considered. These problems arise when a Lavrentiev regularization is utilized for state constrained linearquadratic optimal control problems. For the numerical solution a Galerkin discretization is applied utilizing proper orthogonal decomposition (POD). Based on a perturbation method it is determined how far the suboptimal control, computed on the basis of the POD method, is from the (unknown) exact one. Numerical examples illustrate the theoretical results. In particular, the POD Galerkin scheme is applied to a problem with state constraints. Introduction In this paper we consider a certain class of linear-quadratic optimal control problems governed by a linear evolution problem and mixed control-state constraints. Due to the following reasons mixed control-state constraints are of The authors gratefully acknowledge support by the German Science Fund DFG grant VO 658/- A-posteriori-POD Error Estimators for Nonlinear Optimal Control Problems governed by Partial Differential Equations. The first author is further supported by the Landesgraduiertenförderung of Baden-Wurttemberg. Martin Gubisch University of Konstanz Department of Mathematics and Statistics Universitätsstraße D Konstanz, Germany Martin.Gubisch@uni-konstanz.de Stefan Volkwein University of Konstanz Department of Mathematics and Statistics Universitätsstraße D Konstanz, Germany Stefan.Volkwein@uni-konstanz.de
4 Martin Gubisch, Stefan Volkwein interest: () They arise in Lavrentiev-type regularizations of state constrained problems and () they may appear in their own rights (e.g., if the control is resticted by a multiple of the state). For the numerical solution we apply a Galerkin approximation, which is based on proper orthogonal decomposition (POD). Recall that POD is a method for deriving reduced-order models of dynamical systems; see [], for instance. In order to ensure that the POD suboptimal solutions are sufficiently accurate, we derive an a-posteriori error estimate for the difference between the exact (unknown) optimal control and its POD suboptimal approximations. Moreover, it is shown that this error tends to zero if the number of POD basis functions in the Galerkin ansatz is increased. The proof relies on a perturbation argument [4] and a convergence analysis for the POD Galerkin scheme, where we make use of a modified POD approximation [6]. Although we can transform the optimal control problem with mixed control-state constraints into a purely control constrained optimal control problem, we can not directly apply the results from [], because the transformation itself depends on the POD discretization as well. Furthermore, we propose a new POD Galerkin ansatz for state and adjoint equations which avoids discretization errors coming from approxiamtions of the initial values. In the numerical examples we combine the a-posteriori error estimator with an adaptive basis update strategy; see []. Of course, other strategies can be applied as well; see [], for instance. Although linear-quadratic optimal control problems with mixed-control constraints can be cast into linear-quadratic optimal control problems with purely control constraints, we can not simply apply the results from []. This is due to the fact, that the transformation itself depends on the POD discretization as well. Let us mention that the a- posteriori analysis can also be utilized for nonlinear problems in an inexact sequential quadratic programming (SQP) approach, where in each level of the SQP iteration a linear-quadratic optimal control problem has to be solved. For instance, this is done in [] utilizing the a-posteriori analysis from []. The paper is organized as follows: In Section we introduce our linearquadratic optimal control problems and review first-order optimality conditions. The a-posteriori error analysis is carried out in Section 3. Section 4 is devoted to the POD approximation and the POD convergence analysis. Finally, numerical test examples are studied in Section 5. The optimal control problem In this section we introduce a class of linear-quadratic optimal control problems. We recall the associated first-order optimality conditions and formulate the optimization problem as a reduced problem for the control variable only. For the solution of the optimal control problem we apply a primal-dual active set strategy which is equivalent to a semismooth Newton method [7]. To utilize the error analysis presented in [, ] we transform the reduced problem to an optimal control problem which is governed by bilateral control constraints for a transformed control variable.
5 POD a-posteriori error analysis for mixed control-state constraints 3. Problem formulation Let V and H be real, separable Hilbert spaces and suppose that V is dense in H with compact embedding. In particular, there exists a constant C V > such that ϕ H C V ϕ V ϕ V. (.) By, H and, V we denote the inner products in H and V, respectively. Let T > be the fixed final time. For t [,T] we define a time-dependent symmetric bilinear form a(t;, ) : V V R satisfying a(t;ϕ,ψ) α ϕ V ψ V ϕ V a.e. in [,T], (.a) a(t;ϕ,ϕ) α ϕ V α ϕ H ϕ V a.e. in [,T] (.b) for constants α, α > and α which do not depend on t. In (.), the abbreviation a.e. stands for almost everywhere. By identifying H with its dual H it follows that V H = H V, each embedding being continuous and dense. Recall that the space W(,T) W(,T) = { ϕ L (,T;V) ϕt L (,T;V ) } is a Hilbert space endowed with the common inner product [3, pp ]. The control space is given by U = L (,T;R m ) with m N. In particular, we identify U with its dual space U. For u U, y H and f L (,T;V ) we consider the linear evolution problem d dt y(t),ϕ H +a(t;y(t),ϕ) = (f +Bu)(t),ϕ V,V ϕ V a.e. in [,T], y() = y in H, (.3) where, V,V stands for the dual pairing between V and its dual space V and B : U L (,T;V ) is a continuous, linear operator. It is known that for every f L (,T;V ), u U and y H there is a unique weak solution y W(,T) satisfying (.3) and ( ) y W(,T) C y H + f L (,T;V ) + u U (.4) for a constant C > which is independent of y, f and u. For a proof of the existence of a unique solution we refer to [3, pp. 5-5]. The a-priori error estimate follows from standard variational techniques and energy estimates. If f+bu L (,T;H), a(t;, ) = a(, ) (independent of t) and y V hold, we have y L (,T;V) H (,T;H); see [3, pp ] and [5, pp ].
6 4 Martin Gubisch, Stefan Volkwein Remark. We split the solution to (.3) in one part, which depends on the fixed initial condition y and right-hand side f, and another part depending linearly on the control variable. Let ŷ W(,T) be the unique solution to the problem d dt ŷ(t),ϕ H +a(t;ŷ(t),ϕ) = f(t),ϕ V,V We define the subspace ŷ() = y in H. W (,T) = { ϕ W(,T) ϕ() = in H } ϕ V a.e. in [,T], endowedwiththetopologyofw(,t).letusnowintroducethelinearsolution operator S : U W (,T): for u U the function y = Su W (,T) is the unique solution to d dt y(t),ϕ H +a(t;y(t),ϕ) = (Bu)(t),ϕ V,V ϕ V a.e. in [,T]. From y W (,T) it follows that y() = in H. The boundedness of S follows from (.4). Now, the solution to (.3) can be expressed as y = ŷ +Su. We introduce the Hilbert space X = W(,T) U endowed with the natural product topology. Let I : L (,T;V) U be a bounded linear operator. By X ad we denote the closed, convex and bounded subset X ad = { (ŷ +Su,u) X ua εu+i(ŷ +Su) u b in R m a.e. in [,T] }, where u a,u b U satisfy u a u b componentwise in R m a.e. and ε > holds. The cost function J : X R is given by J(x) = σ Q T y(t) y Q (t) H dt+ σ Ω y(t) y Ω H + σ u U (.5) for x = (y,u) X, where (y Q,y Ω ) L (,T;H) H are desired states. Furthermore, σ Q,σ Ω and σ >. The optimal control problem is given by minj(x) subject to (s.t.) x X ad. (P) Applying standard arguments [5] one can prove that there exists a unique optimal solution x = (ȳ,ū) to (P). The uniqueness follows from the strict convexity properties of the objective functional on X ad. Throughout this paper, a bar indicates optimality. Next we formulate the first-order sufficient optimality conditions of (P) (see [], for instance):
7 POD a-posteriori error analysis for mixed control-state constraints 5 Theorem. Suppose that the feasible set X ad is nonempty and that (ȳ,ū) X ad is the solution to (P). Then there exists a unique Lagrange multiplier pair ( p, λ) X satisfying together with (ȳ,ū) the primal-dual system ȳ = ŷ +Sū, (.6a) d dt p(t),ϕ H +a(t; p(t),ϕ)+ (I λ)(t),ϕ V,V +σ Q ȳ(t) y Q (t),ϕ H = ϕ V a.e., (.6b) p(t)+σ Ω (ȳ(t) yω ) = in H, σū(t) (B p)(t)+ε λ(t) = in R m a.e., (.6c) λ(t) = max (, λ(t)+γ((iȳ)(t)+εū(t) u b (t)) ) (.6d) +min (, λ(t)+γ((iȳ)(t)+εū(t) u a (t)) ) in R m a.e., where B : L (,T;V) U and I : U L (,T;V ) denote the adjoints of B and I, respectively. Furthermore, γ is an arbitrary real number. In (.6d) the max- and min-operations are interpreted componentwise in the pointwise everywhere sense. Remark. Analogous to Remark. we split the adjoint variable into one part depending on the fixed desired states and into two other parts, which depend linearly on the control variable and on the multiplier λ. Recall that we have defined ŷ as well as the operator S in Remark.. For given y Q L (,T;H) and y Ω H let ˆp = ˆp W(,T) denote the unique solution to the adjoint equation d dt ˆp(t),ϕ H +a(t; ˆp(t),ϕ) = σ Q (y Q ŷ)(t),ϕ H ϕ V a.e., ˆp(T) = σ Ω ( yω ŷ(t) ) in H. Further, we define the linear, bounded operators A, A : U W(,T) as follows: for given u U the function p = A u is the unique solution to d dt p(t),ϕ H +a(t;p(t),ϕ) = σ Q (Su)(t),ϕ H p(t) = σ Ω (Su)(T) ϕ V a.e., in H and for given λ U the function p = A λ uniquely solves d dt p(t),ϕ H +a(t;p(t),ϕ)+ (I λ)(t),ϕ V,V = p(t) = in H. ϕ V a.e., Then, the solution to (.6b) can be expressed by p = ˆp+A ū+a λ.
8 6 Martin Gubisch, Stefan Volkwein. The reduced problem In Remark. we have introduced the solution operator S. The reduced cost functional Ĵ : U R is defined as Ĵ(u) = J(ŷ +Su,u), u U. We define the set of admissible controls by U ad = { u U ua εu+i ( ŷ +Su ) u b in R m a.e. in [,T] }, which is convex, closed and bounded in U; see [8, Prop..]. Now we consider the reduced optimal control problem: minĵ(u) s.t. u U ad. Clearly, if ū is the optimal solution to (ˆP), then x = (ŷ+sū,ū) is the optimal solution to (P). On the other hand, if x = (ȳ,ū) is the solution to (P), then ū solves ( ˆP). A first-order sufficient optimality condition for the convex linearquadratic problem ( ˆP) is given by the variational inequality (ˆP) Ĵ(ū),u ū U for all u U ad. (.7) Combining (.6c) with (.6d) and choosing γ = σε to prevent the dependency of the min- and max-terms on u, this is equivalent with u = ( B p εn(y,p) ) /σ in U, (.8) where (y, p) solves the coupled nonlinear primal-dual system y = ŷ + σ S( B p εn(y,p) ), p = ˆp+ σ A ( B p εn(y,p) ) +A N(y,p), (.9) and the nonlinearity N : W(,T) W(,T) U, N(y,p) = max(,ε B p+ε σ(iy u b )) +min(,ε B p+ε σ(iy u a )), coincides with the Lagrange multiplier λ..3 The primal-dual active set method To solve (P) numerically, a primal-dual active set strategy is utilized. This method is equivalent to a locally superlinearly convergent semi-smooth Newton algorithm applied to the first-order necessary optimality conditions; see [7, 8]. For given iterates (y k,p k ) W(,T) W(,T), k, and for i =,...,m we introduce the active and inactive sets { A k ai = t [,T] ( εσ B p k +Iy k) } (t) < u i ai(t), A k i = A k ai A k bi, { A k bi = t [,T] ( εσ B p k +Iy k) } (.) i (t) > u bi(t), I k i = [,T]\A k i. In Algorithm we formulate the semismooth Newton method for our problem.
9 POD a-posteriori error analysis for mixed control-state constraints 7 Algorithm (Primal-dual active set strategy) Require: Starting value (y,p ) and maximal iteration number k max. : Set k =, determine the active and inactive sets according to (.). : repeat 3: Compute the solution (y k+,p k+ ) to (.9) and set k = k +. 4: Compute the active and inactive sets according to (.). 5: until (A k ai = Ak ai and A k bi = Ak bi ) or k = k max. 6: Set u = (B p k ελ k )/σ..4 Restatement of (P) as a control constrained problem It is known [6] that (P) can be cast into a purely control constrained optimal control problem. Recall that I and S are linear and continuous operators, i.e. the operator F = ε + IS : U U is linear and bounded, too. The proof of Lemma. can be found in the Appendix. Lemma. Suppose that for every ǫ > there is a constant C ǫ > satisfying (BIϕ)(t) V C ǫ ϕ(t) H +ǫ ϕ(t) V (.) for all ϕ W(,T) a.e. in [,T]. Then, the linear operator F = ε+is has a bounded inverse. Next we introduce the transformed control variable v = Fu U for u U. From Remark. we conclude that the solution y to (.3) satisfies Su = y ŷ. Consequently, v has the representation v = εu+i(y ŷ). Using (y,u) X ad, it follows that v satisfies the bilateral box constraints where we set v a = u a Iŷ and v b = u b Iŷ. Let v a v v b in [,T] a.e., (.) V ad = { v U va v v b a.e. in [,T] }. Let us assume that F has a bounded inverse. We replace the state and control variable in J as follows: J(v) = Ĵ(F v) = J(ŷ +SF v,f v) for v U. Then, we consider the optimal control problem min J(v) s.t. v V ad. ( P) If v is the solution to ( P), then the pair x = (ȳ,ū) with ȳ = ŷ +SF v and ū = F v is the solution to (P).
10 8 Martin Gubisch, Stefan Volkwein Remark.3 Let u U, then v = Fu is given by computing y = Su and choosing v = εu+iy. On the other hand, suppose that we know v U, then u = F v can be calculated by solving d dt y(t),ϕ H +a(t;y(t),ϕ) + ε (BIy)(t),ϕ V,V = ε Bv(t),ϕ V,V ϕ V a.e., y() = in H and choosing u = (v Iy)/ε. Hence, the transformations v u and u v both require to solve a state equation..5 First-order sufficient optimality conditions for ( P) In this subsection we present the first-order optimality conditions for ( P). We suppose that F has a bounded inverse. Assume that v denotes the unique optimal control for ( P). Using a Lagrangian framework [] the optimal control v satisfies together with the corresponding state variable ȳ and Lagrange multiplier q the following first-order optimality conditions for ( P) : ȳ = ŷ +SF v, q = ˆp+A F v σf F v F B q,v v U v V ad, (.3a) (.3b) (.3c) where F : U U stands for the dual operator of F. The operator A can be expressed in terms of the adjoint S : W (,T) L (,T;V). This follows from an adaption of Lemma 4. in [] and Lemma.4 in []: Lemma. Define the linear and bounded operator Θ : W(,T) W (,T) Θy,ϕ W(,T),W(,T) = σ Q T y(t),ϕ(t) H dt+σ Ω y(t),ϕ(t) H for y W(,T) and ϕ W (,T). Then, we have B A = S ΘS..6 A distributed optimal control problem In this subsection we introduce an example for (P) and discuss the presented theory for this application. Let Ω R d, d {,,3}, be an open and bounded domain with Lipschitz-continuous boundary Γ = Ω. For T > we set Q = (,T) Ω and Σ = (,T) Γ. We choose H = L (Ω) and V = H(Ω) endowed with the usual inner products ϕ,ψ H = ϕψdx, ϕ,ψ V = ϕψ + ϕ ψdx Ω Ω
11 POD a-posteriori error analysis for mixed control-state constraints 9 and their induced norms, respectively. In (.) we have C V =. Let χ i H, i m,denotegivencontrolshapefunctions.then,forgivencontrolu U, initial condition y H and inhomogeneity f L (,T;H) we consider the linear heat equation m y t (t,x) y(t,x) = f(t,x)+ u i (t)χ i (x), a.e. in Q, i= y(t,x) =, a.e. in Σ, y(,x) = y (x), a.e. in Ω. (.4) We introduce the time-independent, symmetric bilinear form a(ϕ,ψ) = ϕ ψdx for ϕ,ψ V Ω and the bounded, linear operator B : U L (,T;H) L (,T;V ) as (Bu)(t,x) = m u i (t)χ i (x) for (t,x) Q a.e. and u U. i= Hence, we have α = α = α = in (.). It follows that the weak formulation of (.4) can be expressed in the form (.3). Moreover, the unique weak solution to (.4) belongs to the space L (,T;V) provided y V holds. We choose certain shape functions π,...,π m H and introduce the operator I : L (,T;V) U by (Iϕ)(t) = (I ϕ)(t). (I m ϕ)(t) with (I i ϕ)(t) = Ω π i (x)ϕ(t,x)dx for ϕ L (,T;V) a.e. in [,T]. Then, the mixed control-state constraints have the form u ai (t) εu i (t)+ π i (x)y(t,x)dx u bi (t) a.e. in [,T] Ω for (y,u) X and i,...,m with u ai,u bi L (,T). Notice that I is even a bounded operator from L (,T;H) to U. The adjoint operators B : L (,T;V) U and I : U L (,T;V ) have the explicit representations (Bip)(t) = Ω χ i (x)p(t,x)dx ( i m), (I λ)(t,x) = m λ i (t)π i (x) for (t,x) Q a.e., p L (,T;V) and λ U. In particular, if χ i = π i holds for i m, then we have B = I as well as I = B. i=
12 Martin Gubisch, Stefan Volkwein Let us discuss condition (.). For ϕ W(,T) we have (BIϕ)(t) H for t [,T] a.e. Using (.) with C V = we obtain (BIϕ)(t) V m π i,ϕ(t) H χ i H c ϕ(t) / H ϕ(t) / i= wherewesetc = m i= π i H χ i H.Now (.)followsdirectlyfromyoung s inequality (A.) with a = ϕ(t) / V and b = c ϕ(t) / N /. In particular, we find C ǫ = c /(8ǫ) for every ǫ >. V, 3 A-posteriori error analysis The goal of this section is to derive an a-posteriori error estimate for (ˆP). For that purpose we utilize the technique in [, Section 3] for the control constrained problem ( P). 3. Derivation of the a-posteriori error estimate Suppose that F has a bounded inverse and that u p is an arbitrary control in U ad. Our goal is to estimate the norm ū u p U without the knowledge of the optimal solution ū = F v. We set v p = Fu p, i.e., v p = εu p +I(y p ŷ) with y p = ŷ+su p. If u p ū holds, then v p v. Thus, v p does not satisfy the necessary (and by convexity sufficient) optimality condition (.3c). However, there exists a function ζ U such that σf F v p F B q p +ζ,v v p U v V ad, (3.) with q p = ˆq + Au p. Therefore, v p satisfies the optimality condition of a perturbed parabolic optimal control problem with perturbation ζ. The smaller ζ is, the closer v p is to v. Next we estimate ū u p U in terms of ζ U. By Lemma. we have B ( q q p ) = B A (ū ūp ) = B S ΘS ( u p ū ) = B S Θ ( y p ȳ ), (3.) with ȳ = ŷ +Sū. Choosing v = v p in (.3c), v = v in (3.) and using (3.) we obtain σf (ū u p )+F B ( q q p )+ζ,f(ū u p ) U = σ ū u p U B S Θ(ȳ y p ),ū u p U + ζ,f(ū u p ) U = σ ū u p U Θ(ȳ y p),ȳ y p W(,T),W(,T) + F ζ,ū u p U = σ ū u p U + F ζ,ū u p U σ ū u p U + F ζ U ū u p U. Hence, we get the a-posteriori error estimation ū u p U σ F ζ U (3.3)
13 POD a-posteriori error analysis for mixed control-state constraints If the evaluation of F is much more expensive compared to the calculation of ζ, it may be advisable to deduce a-priori estimates for the operator norm of F by the application of energy estimates for the state equations. However, these bounds are not rigorous in general, i.e. the true value of F ζ U might be overestimated significantly. 3. Computation of a perturbation ζ Now we want to determine an appropriate perturbation ζ U satisfying (3.). Supposethatthesuboptimalcontrolu p = F v p U isknown,wederivefrom (3.) the variational inequality ξ +ζ,v v p U v V ad, (3.4) where we set ξ = F (σu p B q p ) U. Hence, the element ξ solves ( ε+s I ) ξ = σu p B q p in U. (3.5) To determine an appropriate perturbation ζ we have to compute ξ. This is formulated in the next theorem. Its proof follows from Proposition 3. in []. Theorem 3. (A-posteriori error estimate for (ˆP)) Suppose that u p is an arbitrary control in U ad and q p = ˆq+A u p. Moreover, the function ξ solves (3.5). Define ζ U as follows: min(,ξ i (t)) a.e. in A p ai = { t [,T] v pi (t) = v ai (t) }, ζ i (t) = max(,ξ i (t)) a.e. in A p bi = { t [,T] v pi (t) = v bi (t) }, ξ i (t) a.e. in [,T]\ ( (3.6) A p ) ai Ap bi for i m. Then, we have ū u p U ǫ ape (3.7) with the a-posteriori error estimator ǫ ape = F ζ U /σ. We call (3.7) an a-posteriori error estimate, since, in the next section, we will apply it to suboptimal controls u p that have already been computed from a POD model. 4 Galerkin approximation for (P) and ( P) The goal of this section is to apply the a-posteriori error analysis for a suboptimal control u p = ū l, which is computed by a POD Galerkin scheme for (ˆP). In particular, we prove that the a-posteriori error estimator tends to zero if the number of POD ansatz functions tends to infinity. Let us mention that we avoid a discretization of the control space U; see also [9].
14 Martin Gubisch, Stefan Volkwein 4. Proper orthogonal decomposition (POD) Suppose that w,...,w k L (,T;V) are k given trajectories. Moreover, we introduce the linear subspace { } V = span w k (t) t [,T] and k {,...,k} with dimension d. The method of POD consists in choosing an orthonormalbasisinv suchthatforeveryl {,...,d}themeansquareerrorbetween theelementsw k (t),t [,T]and k k,andthecorrespondingl-thpartial sum is minimized on average: with w k min ψ,..., ψ l V l w k, ψ i V ψi i= k w k k= l w k, ψ i V ψi i= s.t. ψ i, ψ j V = δ ij, i,j l L (,T;V) = T w k (t) i= L (,T;V) V (4.) l w k (t), ψ i V ψi dt. A solution to (4.) is called a POD basis of rank l. We introduce the linear, bounded, nonnegative operator R : V V by Rψ = k k= T w k (t),ψ V w k (t)dt for ψ V. The solution of (4.) can be found in [,], for instance. Proposition 4. For w,...,w k L (,T;V) the linear operator R is nonnegative, self-adjoint and compact. Let {λ i } i N and {ψ i } i N denote the nonnegative eigenvalues and associated orthonormal eigenfunctions of R satisfying Rψ i = λ i ψ i, λ λ..., and λ i as i. Then a POD basis of rank l d is given by {ψ i } l i=, and we have k w k k= l w k,ψ i V ψ i i= L (,T;V) = i=l+ λ i. (4.) Remark 4. From the Hilbert-Schmidt theorem [8, p. 3] it follows that {ψ i } i N form a complete orthonormal basis for the separable space V.
15 POD a-posteriori error analysis for mixed control-state constraints 3 4. The POD Galerkin approximation In this subsection we introduce the POD schemes for the first-order optimality system using a POD Galerkin approximation for the primal and dual variables. Moreover, we study the convergence of the POD discretizations, where we makeuseoftheanalysisin[,4,].wemakeuseofthefollowinghypothesis. Assumption Let u,λ U be chosen such that the functions Su, A u and A λ are nonzero and belong to H (,T;V). In the context of Section 4. we choose k = 6, w = Su, w = (Su) t, w 3 = A u, w 4 = (A u) t, w 5 = A λ and w 6 = (A λ) t. For l we denote by {ψ i } l i= a POD basis of rank l and set V l = span{ψ,...,ψ l }. Let us define the linear and bounded projection operator P l ψ = l ψ,ψ i V ψ i V l for ψ V. i= Suppose that (ȳ,ū) X is the solution to ( P) and ( p, λ) the associated unique Lagrange mulitplier pair. Analogous to the operator S introduced in Remark. we define the operator S l : U W (,T) as follows: for u U the function y l = S l u is the unique solution to d dt yl (t),ψ H +a(t;y l (t),ψ) = (Bu)(t),ψ V,V ψ V l a.e. Notice that for any u U the element S l u belongs even to H (,T;V l ) which is continuously embedded into W(,T). Similar to Remark. we introduce the linear and bounded operators A l, A l : U W(,T): for given u U the function p l = A l u uniquely solves d dt pl (t),ψ H +a(t;p l (t),ψ) = σ Q (S l u)(t),ψ V,V p l (T) = σ Ω (S l u)(t) ψ V l a.e., in H and for given λ U the function p l = A λ is the unique solution to d dt pl (t),ψ H +a(t;p l (t),ψ) = (I λ)(t),ψ V,V ψ V l a.e., p l (T) = in H. Now, the POD Galerkin scheme for (.9) is given by y l = ŷ + σ Sl( B p l εn(y l,p l ) ), p l = ˆp+ σ Al ( B p l εn(y l,p l ) ) +A l N(y l,p l ). (4.3a)
16 4 Martin Gubisch, Stefan Volkwein Then we set λ l = N(y l,p l ), u l = ( B p l ελ l ) ) /σ in U. (4.3b) For a proof of the following proposition we refer the reader to [6, Theorems 3. and 3.]. Proposition 4. Let Assumption hold. If Sũ even belongs to H (,T;V) for every ũ U, then we have lim S l Sl L(U,W(,T)) = lim A A l l L(U,W(,T)) =, where L(U,W(,T)) stands for the Banach space of all linear and bounded operators from U to W(,T) equipped with its natural operator norm. To introduce a POD Galerkin scheme for (.3) let us introduce the POD aproximation F l : U U of the operator F by F l = ε+is l : U U. The proof of the following proposition is given in the Appendix. Proposition 4.3 Suppose that Assumption is valid. Then, we have: ) F l is linear, bounded (uniformly with respect to l), lim l F F l L(U) =. ) Suppose that F has a bounded inverse. Then, F l has a bounded inverse as well. We write F l, = (F l ). If L N be chosen arbitrarily large so that F L(U) F F l L(U) < holds for all l > L, we obtain F l, L(U) F L(U) F L(U) F F l L(U) for all l L. 3) We have lim l FF l, L(U) =, lim l F l, F L(U) = and lim l F F l, L(U) =, where F l, denotes the adjoint operator of F l,. Let us formulate a discrete version of the primal-dual active set method which is utilized in our numerical tests to solve (4.3a) and (4.3b). For given iterates (y lk,p lk ) W(,T) W(,T), k, and for i =,...,m we introduce the active and inactive sets { } A lk ai = A lk bi = t [,T] (εσ B p lk +Iy lk ) i (t) < u ai (t) { } t [,T] (εσ B p lk +Iy lk ) i (t) > u bi (t), A lk i = A lk ai A lk bi,, I lk i = [,T]\A lk i ; (4.4) compare (.). In Algorithm we state the semismooth Newton method for the POD discretized problem.
17 POD a-posteriori error analysis for mixed control-state constraints 5 Algorithm (POD discretized primal-dual active set strategy) Require: POD basis ψ l, starting value (y l,p l ) and maximal iteration number k max. : Set k =, determine the active and inactive sets according to (4.4). : repeat 3: Determine the solution (y l,k+,p l,k+ ) to (4.3a) and set k = k +. 4: Compute the active and inactive sets according to (4.4). 5: until (A lk ai = Al,k ai and A lk bi = Al,k bi ) or k = k max. 6: Set u l = (B p lk ελ lk )/σ. Utilizing the operator F l we introduce the POD Galerkin scheme for (.3) as follows ȳ l = ŷ +S l F l, v l, q l = ˆp+A l F l, v l σf l, F l, v F l, B q l,v v l U v V ad, (4.5a) (4.5b) (4.5c) In contrast to the error analysis in [, ] the discretized variational inequality involves the operators F l. This reflects the fact that the POD discretization U l ad = { u U u a εu+i ( ŷ +S l u ) u b in R m a.e. in [,T] } oftheadmissible setu ad dependsonthepodgalerkinscheme,whereasin[, ] the admissible set for the controls is independent of the POD discretization. On the other hand, V ad is independent of l. In particular, we have v l = F l ū l. The proof of the next theorem is given in the Appendix. Theorem 4. Suppose that Assumption holds. ) For q and q l we have lim l q q l W(,T) = ) Let v and v l be the solutions to (.3c) and (4.5c), respectively. Then, lim v l vl U =. 3) If ū and ū l are the solutions to (.7) and (4.3), respectively, we obtain lim l ū ūl U =. 4) Define, according to (3.6), the function ζ l U by min(,ξi(t)) l a.e. in A l ai = { t [,T] v i(t) l = v ai (t) }, ζi(t) l = max(,ξi(t)) l a.e. in A l bi = { t [,T] vi(t) l = v bi (t) }, ξi(t) l a.e. in [,T]\ ( A l ai A l bi), where ξ l U solves F ξ l = σū l B q l in U; compare (3.5). Then, and lim l ζ l U =. ū ū l U ǫ ape with ǫ ape = F ζ l U σ (4.6)
18 6 Martin Gubisch, Stefan Volkwein Remark 4. ) In addition to the a-posteriori error analysis in [], we have to solve the linear system F ξ l = σū l B q l. ) Part 3) of Theorem 4. shows that ζ l U can be expected smaller than any ǫ > provided that l is taken sufficiently large. Motivated by this result we set up Algorithm 3. Algorithm 3 (POD reduced-order method with a-posteriori estimator) Require: Initial control u l U, initial number l for the POD ansatz functions, a maximal number l max > l of POD ansatz functions, and a stopping tolerance ǫ >. : Determine ŷ, ˆq, w = Su l, w = Au l. : Compute a POD basis {ψ i } lmax i= choosing w and w. Set l =. 3: repeat 4: Establish the POD Galerkin discretization using {ψ i } l i=. 5: Call Algorithm to compute suboptimal control ū l. 6: Determine ǫ ape from (4.6). 7: if ǫ ape < ǫ or l = l max then 8: Return l and suboptimal control ū l and STOP. 9: end if : Set l = l+. : until l > l max 5 Numerical experiments In this section we carry out numerical test examples for the presented theoretical findings. Run 5. In the context of Section.6 we choose d =, Ω = (,), Ω i = m [i,i] for i m, χ i(x) = π i (x) = χ Ωi (x) the characteristic functions on the subdomains Ω i and T = 3. Let σ = e-3, σ Q =, σ Ω =, ε = e-5 and f(t,x) = t x 3, y Q (t,x) = t( (x ) ), y Ω (x) = y Q (T,x), y (x) = (χ [.4,.] χ [.,.6] )(x).
19 POD a-posteriori error analysis for mixed control-state constraints 7 optimal control optimal state Bū(t, x) 5 5 ȳ(t, x) direction x.5 time t 3.5 direction x.5 time t 3 Fig. 5. Run 5.: The optimal control Bū and the optimal state ȳ for u a =.5, u b =.5 and m = calculated by solving the full order model. Notice that for large m and small regularization ε, the mixed control-state constraints can be interpreted as pointwise state constraints: As one sees in Figure 5., m u a y m u b, i.e..5 y(t,x).5, holds approximatively. 5 Rom error and Aposti error bound A posteriori bound ROM error Singular values 5 Rom error and Aposti error bound A posteriori bound ROM error Singular values 5 Error(l) Error(l) Pod basis rank l Pod basis rank l Fig. 5. Run 5.: The ROM errors of the control ū l for different POD basis ranks l with initial control guesses u l (left) and u l = ū (right). In Figure 5. we investigate the decay of the ROM errors for increasing POD basis rank l. For this purpose, we apply Algorithm for the (fixed) POD bases corresponding to the snapshots Su l with u l and u l = ū. One sees that increasing the basis rank of the arbitrarily chosen POD basis does not lead to a satisfying accuracy of the reduced order model solution since the state solution Su l does not cover enough of the dynamics of the optimal state ȳ, so POD basis updates which exploite the information gained
20 8 Martin Gubisch, Stefan Volkwein from the actual suboptimal control are required to get snapshots which fit to ū. Algorithm 4 (POD Primal-dual active set strategy with basis adaptivity) Require: Initial control guess u l, tolerance ǫ and maximal iterations k max. : Set k =. : repeat 3: Compute the solution y k = ŷ +Su kl W(,T) to the system (.3). 4: Compute a rank-l POD basis ψ l V by solving (4.) with w = y k. 5: Execute Algorithm and set u kl = ū l, k = k +. 6: until ǫ ape(u kl ) < ǫ or k = k max. 7: Return ū l = u kl. Since in practice we do not have enough knowledge about ū to choose a- priori an admissible POD basis, we present an adaptive strategy in Algorithm 4 which actualizes the snapshots allocations iteratively with the information gained from the ROM solutions. Now, one single basis update is sufficient to achieve the same decay order of the error caused by the model reduction as we have for the (in general unknown) optimal POD basis. In contrast to our experiences with pure control constrained optimal control problems, the a posteriori error bounds are not rigorous, i.e. the true ROM error has a smaller order than the error estimation. This is due to the fact that we have chosen the regularization parameter ε of the state constraints very small which leads to a badly scaled discretization of the differential equations required to determine the perturbation vector ζ. For larger values of ε, the error bounds are rigorous if the POD basis rank l is not chosen too small, see Figure 5.3. In our example, the primal-dual active set strategy does not converge within the claimed maximal number of active set actualizations if l <. 5 Rom error and Aposti error bound A posteriori bound ROM error Singular values 5 Rom error and Aposti error bound A posteriori bound ROM error Singular values Error(l) 5 Error(l) Pod basis rank l Pod basis rank l Fig. 5.3 Run 5.: The ROM errors of the control ū l for different POD basis ranks l with adaptive basis selections and regularization parameters ε = e-5 (left) and ε = 5e-3 (right).
21 POD a-posteriori error analysis for mixed control-state constraints 9 Process Time # Total Assemble full system (.9).66 sec sec Solve full system (.9).7 sec 9.43 sec Total 6.4 sec Solve full snapshots equations (.3a). sec. sec Solve eigenvalue problem (4.).4 sec.84 sec Assemble ROM system (4.3a).53 sec 7 9. sec Solve ROM system (4.3a).45 sec sec Evaluate error estimator (4.6). sec.3 sec Total 8. sec Table Run 5.: The calculation times for solving the optimization problem with the primal-dual active set strategy with and without model reduction. With 5 POD elements, the reduced-order problem has to be solved two times; solvings of two eigenvalue problems are required in addition for updating the POD basis. Nevertheless, 9.7% of the calculation time is spared in total where the true error is 9.7e-5 and the a-posteriori error bound totals up to 5.68e-4. We finish the first example with a look at Table where the effort of the model reduction on the calculation times is illustrated. Since the most expensive part of the optimization process is the simultaneous solving of the primal and dual equations, ROM is very effective here even for POD basis ranks which are chosen so large that the accuracy of the full order model is reached. Run 5. Choosing ε = and I =, the techniques presented also cover the case of pure control constraints. The optimal control and state solutions induced by the same data functions as we used before for the state constrained problem are shown in Figure 5.4: optimal control optimal state.8.5 Bū(t, x).6 ȳ(t, x) direction x.5 time t 3.5 direction x.5 time t 3 Fig. 5.4 Run 5.: The optimal control Bū and the optimal state ȳ for u a =.5, u b =.75 and m = calculated by solving the full order model.
22 Martin Gubisch, Stefan Volkwein In this case, the a-posteriori error bounds are rigorous. Figure 5.5 illustrates that even the arbitrary POD basis corresponding to the initial control guess u l coveres enough dynamics of ū to decrease the ROM errors below the FEM accuracy which is of the order.e-4 without basis updates. 5 Rom error and Aposti error bound A posteriori bound ROM error Singular values Rom error and Aposti error bound A posteriori bound ROM error Singular values 5 4 Error(l) Error(l) Pod basis rank l Pod basis rank l Fig. 5.5 Run 5.: The ROM errors of the control ū l for different POD basis ranks l with a fixed POD basis (left) and with basis adaptivity (right). For practical application it does not make sense to enlarge the reduced order model, i.e. the POD basis rank, as soon as the a-posteriori error estimator indicates that the exactness of the full order model is reached since the FEM error is part of the snapshots and hence of the reduced order model anyway. Process Time # Total Assemble full system (.9).79 sec sec Solve full system (.9) 8.53 sec sec Total 77.5 sec Solve full snapshots equations (.3a). sec. sec Solve eigenvalue problem (4.).7 sec.7 sec Assemble ROM system (4.3a).4 sec 4.6 sec Solve ROM system (4.3a). sec 4.47 sec Evaluate error estimator (4.6).3 sec.3 sec Total.48 sec Table Run 5.: The calculation times for solving the optimization problem with the primal-dual active set strategy with and without model reduction. With 5 POD elements, the reduced-order problem has to be solved two times; solvings of two eigenvalue problems are required in addition to update the POD basis. Nevertheless, 96.79% of the calculation time is spared in total and both the true error as well as the a-posteriori error bound amount to 5.6e-6.
23 POD a-posteriori error analysis for mixed control-state constraints Run 5.3 We consider now a two-dimensional setting: Let Ω = (,π) (,π), Ω ij = π m [i,i] [j,j] for i,j m, χ ij(x) = π ij (x) = χ Ωij (x) be the characteristic functions on the subdomains Ω ij and T = π. Let σ = e-3, σ Q =, σ Ω =, ε = e-5. We choose f, y and consistent desired states y Q (t,x) = sin(t)sin(x )sin(x ), y Ω = y Q (T, ). Obviously, y = y Q is an optimal state solution to the optimization problem without state or control constraints if σ =. Figure 5.6 shows that the state solution which respects the weakly regularized pointwise state constraints resembles the projection of y Q on the admissible range. optimal control optimal state Bū(T, x, x) 3 ȳ(t, x, x) direction x direction x 3 3 direction x direction x 3 Fig. 5.6 Run 5.3: The optimal control (Bū)(T) and the optimal state ȳ(t) for u a,ij =.5 Ω ij, u b,ij = +.5 Ω ij and m = 4 calculated by solving the full order model. An arbitrary POD basis is not able to establish the area of the graph where the upper constraints are active in this case, see Figure 5.7: Basis updates are required instead to build up an accurate reduced-order model.
24 Martin Gubisch, Stefan Volkwein Rom error and Aposti error bound A posteriori bound ROM error Singular values Rom error and Aposti error bound A posteriori bound ROM error Singular values Error(l) Error(l) Pod basis rank l Pod basis rank l Fig. 5.7 Run 5.3: The ROM errors of the control ū l for different POD basis ranks l with a fixed POD basis (left) and with basis adaptivity (right). Process Time # Total Assemble full system (4.3a).5 sec sec Solve full system (4.3a) 3.6 sec sec Total 64. sec Solve full snapshots equations (.3a). sec.4 sec Solve eigenvalue problem (4.). sec.9 sec Assemble ROM system (4.3a). sec sec Solve ROM system (4.3a).3 sec 5.69 sec Evaluate error estimator (4.6).54 sec.8 sec Total 9.56 sec Table 3 Run 5.3: The calculation times for solving the optimization problem with the primal-dual active set strategy with and without model reduction. With 5 POD elements, the reduced-order problem has to be solved two times; solvings of two eigenvalue problems are required in addition to update the POD basis. Nevertheless, 99.4% of the calculation time is spared in total. The a posteriori error bound is of the same order as the ROM error which amount to 8.74e-. Appendix Proof of Lemma. We utilize arguments from the proof of Lemma 4. in [6]. First we show that F is injective. From Fu = we infer (ε+is)u = in U. Since ε > holds, we derive u = ISu in U. ε By Remark., we conclude that y = Su W(,T) satisfies y t(t),ϕ V,V +a(t;y(t),ϕ) = ε (BIy)(t),ϕ V,V, ϕ V a.e. in [,T], y(),ϕ H = ϕ H.
25 POD a-posteriori error analysis for mixed control-state constraints 3 Choosing ϕ = y(t) and using (.), (.) and Young s inequality [5, p. 6], ab ǫa + b ǫ for all a, b R, ǫ > (A.) with a = y(t) V, b = y(t) H /(C ǫε), we deduce that d dt y(t) H +α y(t) V α y(t) H d dt y(t) H +a(t;y(t),y(t)) ε (BIy)(t) V y(t) V ε (Cǫ y(t) H +ǫ y(t) V ) y(t) V ( ǫ + ǫ ) y(t) V ε + Cǫ ǫ ε y(t) H a.e. in [,T]. Taking ǫ = εα /4 and ǫ = α /, we get the a-priori energy estimate d dt y(t) H + α ( ) y(t) V α + C ǫ ε y(t) H a.e. in [,T]. α (A.) Applying Gronwall s lemma [5, pp ] and y() = in H, we obtain by standard arguments for linear evolution problems that y holds in W(,T). Thus, the operator F is injective. Next we prove that F is surjective. Let v U be chosen arbitrarily. Then, F is surjective if there exists an u U satisfying εu+isu = v in U. Deriving again an a-priori energy estimate analogously to (A.), we conclude that d dt y(t),ϕ H +a(t;y(t),ϕ)+ ε BIy(t),φ V,V = ε Bv(t),ϕ V,V y(),ϕ H = ϕ V a.e., ϕ H admits a unique solution y W(,T). Defining u = (v Iy)/ε U, y = Su holds, i.e. Fu = (ε + IS)u = v in U which implies the surjectivity of F. Now the claim follows from the bounded inverse theorem. Proof of Proposition 4.3 ) Clearly, F l is linearand bounded by a constant which is independent of l. Since F F l L(U) = I(S S l ) L(U) I L(L (,T;V ),U) S Sl L(U,W(,T)) is satisfied, the convergence follows directly from Proposition 4.. ) By part ) there exists a constant L N so that F F l L(U) < F L(U) for all l L. Then, the claim follows from the perturbation lemma [7, p. 45]. 3) Using FF l, L(U) = (F l F)F l, L(U) F l F L(U) F l, L(U), F l, F L(U) = F l, (F l F) L(U) F l, L(U) F l F L(U) and parts ), 3) the limits hold.
26 4 Martin Gubisch, Stefan Volkwein Proof of Theorem 4. ) The claim follows from [6, Theorem 3.]. ) Choosing v = v l in (.3c) and v = v in (4.5c) we get the variational inequality σ ( F F v F l, F l, v l) F l, B q l +F B q, v l v. U (A.3) From v l V ad we infer that v l is bounded. Thus, there exists a constant C > with v l U C. Moreover, F l, L(U) C for l sufficiently large with C = F L(U). Thus, σ F F v F l, F l, v l, v l v U = σ F F ( v v l )+F F v l F l, F l, v l, v l v U = σ F ( v v l ) U +σ (F F l, ) v l,f ( v l v) U + ( F F l, )F l, v l,f ( v l v) U σ ( F ( v v l ( ) U +C 3 FF l, L(U) + F F l, ) ) L(U) F ( v l v) U with C 3 = C C. We set C 4 = C max(c, B A F L(U) and C 5 = B L(U). Recall that the operator S l is bounded independently of l. Hence, the constant C 6 = C C S l, ΘS l L(U) does not depend on l. Moreover, B A l = Sl, ΘS l is uniformly bounded. Hence, there exists a constant C 7 > which does not depend on l so that B A l L(U) C 7. Consequently, Lemma. implies that F B q F l, B q l, v l v U = F B A F ( v v l) +F B ( A F F l, ) v l, v l v U + F B ( A A l ) F l, v l + ( F F l, B A l F l, v l, v l v U ΘSF ( v v l),sf ( v l v ) U +C 4 ( FF l, L(U) + B (A A l ) L(U) ) F ( v v l) U +C 6 F F l, L(U) F ( v v l) U C 7 ( FF l, L(U) + F F l, L(U) + B (A A l ) L(U) ) F ( v v l) U with C 7 = max(c 4,C 6 ). Hence, we have F ( v v l ) U C 8 σ + C 8 σ ( FF l, L(U) + F F l, L(U) ) B (A A l ) L(U) with C 8 = (C 3 +C 7 ). From Proposition 4. we infer that lim B (A A l ) l L(U) = so that F ( v v l ) U for l. Since F is invertible, we conclude that lim v l vl U =.
27 POD a-posteriori error analysis for mixed control-state constraints 5 3) From ū = F v, ū l = F l, v l and Theorem 4. we infer that ū ū l U = F v F l, ū l U F L(U) ( v v l U + FF l, L(U) v l U ) l. 4) The first part of the claim follows directly from Theorem 3.. We infer from part ) and part 3) that {B q l } l N and {ū l } l N converge to ū respectively B q. Hence, {σū l B q l } l N tends to σū B q. Since F is bounded, we conclude that {ξ l } l N converge to ξ = F (σū B q). Now, the proof follows by the same arguments as the proof of Theorem 4., part (), in []. References. K. Afanasiev and M. Hinze. Adaptive control of a wake flow using proper orthogonal decomposition. Lecture Notes in Pure and Applied Mathematics, 6: 37-33,.. E. Arian, M. Fahl, and E.W. Sachs. Trust-region proper orthogonal decomposition for flow control. Technical Report -5, ICASE,. 3. R. Dautray and J.-L. Lions. Mathematical Analysis and Numerical Methods for Science and Technology. Volume 5: Evolution Problems I. Springer, Berlin,. 4. A. L. Dontchev, W. W. Hager, A. B. Poore, and B. Yang. Optimality, stability, and convergence in nonlinear control. Applied Math. and Optimization, 3:97 36, L.C. Evans. Partial Differential Equations. American Math. Society, Providence, Rhode Island, M. Gubisch and S. Volkwein. Modified POD Galerkin schemes for linear-quadratic optimal control problems. In preparation, 3 7. M. Hintermüller, K. Ito, and K. Kunisch. The primal-dual active set strategy as a semismooth Newton method. SIAM Journal on Optimization, 3: , M. Hintermüller, I. Kopacka, and S. Volkwein. Mesh-independence and preconditioning for solving control problems with mixed control-state constraints. ESAIM: Control, Optimisation and Calculus of Variations, 5:66-65, M. Hinze and F. Tröltzsch. Discrete concepts versus error analysis in pde constrained optimization. GAMM-Mitteilungen, 33:48-6,.. M. Hinze and S. Volkwein. Error estimates for abstract linear-quadratic optimal control problems using proper orthogonal decomposition. Computational Optimization and Applications, 39:39-345, 8.. P. Holmes, J.L. Lumley, G. Berkooz, and C.W. Rowley. Turbulence, Coherent Structures, Dynamical Systems and Symmetry. Cambridge Monographs on Mechanics, Cambridge University Press, second edition,.. M. Kahlbacher and S. Volkwein. POD a-posteriori error based inexact SQP method for bilinear elliptic optimal control problems. ESAIM: Mathematical Modelling and Numerical Analysis, 46:49-5,. 3. E. Kammann, F. Tröltzsch and S. Volkwein. A method of a-posteriori error estimation with application to proper orthogonal decomposition. ESAIM: Mathematical Modelling and Numerical Analysis, to appear, K. Kunisch and S. Volkwein. Galerkin proper orthogonal decomposition methods for parabolic problems. Numerische Mathematik, 9:7-48,. 5. J.L. Lions. Optimal Control of Systems Governed by Partial Differential Equations. Springer, Berlin, C. Meyer, U. Prüfert, and F. Tröltzsch. On two numerical methods for state-constrained elliptic control problems. Optimization Methods and Software, :87-899, J.M. Ortega and W.C. Rheinboldt Iterative Solution of Nonlinear Equations in Several Variables. Computer Science and Applied Mathematics. Academic Press, M. Reed and B. Simon. Methods of Modern Mathematical Physics I: Functional Analysis. Academic Press, New York, 98.
28 6 Martin Gubisch, Stefan Volkwein 9. A. Studinger and S. Volkwein. Numerical analysis of POD a-posteriori error estimation for optimal control To appear in Control and Optimization with PDE Constraints, International Series of Numerical Mathematics, Kunisch, Clason (eds), Birkhäuser, Basel, 3.. F. Tröltzsch. Optimal Control of Partial Differential Equations. Theory, Methods and Applications. American Math. Society, Providence, volume,.. F. Tröltzsch and S. Volkwein. POD a-posteriori error estimates for linear-quadratic optimal control problems. Computational Optimization and Applications, 44:83-5, 9.. S. Volkwein. Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling. Lecture Notes, University of Konstanz,.
Proper 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 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 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 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 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 informationPOD-Based Economic Optimal Control of Heat-Convection Phenomena. Luca Mechelli Stefan Volkwein
Universität Konstanz POD-Based Economic Optimal Control of Heat-Convection Phenomena Luca Mechelli Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 365, Oktober 2017 ISSN 1430-3558 Konstanzer Online-Publikations-System
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 informationEileen Kammann Fredi Tröltzsch Stefan Volkwein
Universität Konstanz A method of a-posteriori error estimation with application to proper orthogonal decomposition Eileen Kammann Fredi Tröltzsch Stefan Volkwein Konstanzer Schriften in Mathematik Nr.
More informationReduced-Order Methods in PDE Constrained Optimization
Reduced-Order Methods in PDE Constrained Optimization M. Hinze (Hamburg), K. Kunisch (Graz), F. Tro ltzsch (Berlin) and E. Grimm, M. Gubisch, S. Rogg, A. Studinger, S. Trenz, S. Volkwein University of
More informationModel order reduction & PDE constrained optimization
1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell,
More informationProper Orthogonal Decomposition. 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 informationOPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL
OPTIMALITY CONDITIONS AND POD A-POSTERIORI ERROR ESTIMATES FOR A SEMILINEAR PARABOLIC OPTIMAL CONTROL O. LASS, S. TRENZ, AND S. VOLKWEIN Abstract. In the present paper the authors consider an optimal control
More informationReduced-Order 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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
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 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 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 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 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 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 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 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 informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
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 information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
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 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 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 informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
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 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 informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
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 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
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 informationSome SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen
Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text
More informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
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 informationEquations paraboliques: comportement qualitatif
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents
More informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
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 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 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 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 informationConvexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls
1 1 Convexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls B.T.Polyak Institute for Control Science, Moscow, Russia e-mail boris@ipu.rssi.ru Abstract Recently [1, 2] the new convexity
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 informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationUniversité de Metz. Master 2 Recherche de Mathématiques 2ème semestre. par Ralph Chill Laboratoire de Mathématiques et Applications de Metz
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Systèmes gradients par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 26/7 1 Contents Chapter 1. Introduction
More informationL p MAXIMAL REGULARITY FOR SECOND ORDER CAUCHY PROBLEMS IS INDEPENDENT OF p
L p MAXIMAL REGULARITY FOR SECOND ORDER CAUCHY PROBLEMS IS INDEPENDENT OF p RALPH CHILL AND SACHI SRIVASTAVA ABSTRACT. If the second order problem ü + B u + Au = f has L p maximal regularity for some p
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
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 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 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 informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
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 informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More informationRESEARCH ARTICLE. A strategy of finding an initial active set for inequality constrained quadratic programming problems
Optimization Methods and Software Vol. 00, No. 00, July 200, 8 RESEARCH ARTICLE A strategy of finding an initial active set for inequality constrained quadratic programming problems Jungho Lee Computer
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
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 informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 15, No. 1, pp. 39 62 c 2004 Society for Industrial and Applied Mathematics SEMISMOOTH NEWTON AND AUGMENTED LAGRANGIAN METHODS FOR A SIMPLIFIED FRICTION PROBLEM GEORG STADLER Abstract.
More informationKey words. optimal control, heat equation, control constraints, state constraints, finite elements, a priori error estimates
A PRIORI ERROR ESTIMATES FOR FINITE ELEMENT DISCRETIZATIONS OF PARABOLIC OPTIMIZATION PROBLEMS WITH POINTWISE STATE CONSTRAINTS IN TIME DOMINIK MEIDNER, ROLF RANNACHER, AND BORIS VEXLER Abstract. In this
More informationOn m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry
On m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry Ognjen Milatovic Department of Mathematics and Statistics University of North Florida Jacksonville, FL 32224 USA. Abstract
More informationGALERKIN APPROXIMATIONS OF NONLINEAR OPTIMAL CONTROL PROBLEMS IN HILBERT SPACES
Electronic Journal of Differential Equations, Vol. 217 (217), No. 189, pp. 1 4. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GALERKIN APPROXIMATIONS OF NONLINEAR OPTIMAL
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More 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 informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More information