Hamburger Beiträge zur Angewandten Mathematik
|
|
- Kelley Newton
- 6 years ago
- Views:
Transcription
1 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 Januar 2007
2
3 A finite element approximation to elliptic control problems in the presence of control and state constraints Klaus Deckelnick & Michael Hinze Abstract: We consider an elliptic optimal control problem with control and pointwise state constraints. The cost functional is approximated by a sequence of functionals which are obtained by discretizing the state equation with the help of linear finite elements and enforcing the state constraints in the nodes of the triangulation. The control variable is not discretized. Error bounds for control and state are obtained both in two and three space dimensions. Finally, we discuss some implementation issues of a generalized Newton method applied to the numerical solution of the problem class under consideration. Mathematics Subject Classification (2000): 49J20, 49K20, 35B37 Keywords: Elliptic optimal control problem, state constraints, control constraints, error estimates Introduction Let Ê d (d = 2,3) be a bounded domain with a smooth boundary and consider the differential operator d ( ) d Ay := xj aij y xi + b i y xi + cy, i,j= along with its formal adjoint operator A y = i= d ( d xi a ij y xj + b i y ) + cy i= j= where for simplicity the coefficients a ij,b i and c are assumed to be smooth functions on. We associate with A the bilinear form a(y,z) := ( d a ij (x)y xi z xj + i,j= d b i (x)y xi z + c(x)yz ) dx, and subsequently assume that there exists c 0 > 0 such that i= d a ij (x)ξ i ξ j c 0 ξ 2 for all ξ Ê d and all x. i,j= y,z H () Institut für Analysis und Numerik, Otto von Guericke Universität Magdeburg, Universitätsplatz 2, 3906 Magdeburg, Germany Schwerpunkt Optimierung und Approximation, Universität Hamburg, Bundesstraße 55, 2046 Hamburg, Germany.
4 Furthermore we suppose that the form a is coercive on H (), i.e. there exists c > 0 such that a(v,v) c v 2 H for all v H (). (.) From the above assumptions it follows that for a given f (H ()) the elliptic boundary value problem Ay = f in d i,j= a iju xi ν j = 0 on (.2) has a unique weak solution y H () which we denote by y = G(f). Here, ν is the unit outward normal to. Furthermore, if f L 2 (), then the solution y belongs to H 2 () and satisfies y H 2 C f, where have used to denote the L 2 norm. Next, let (U,(, ) U ) be a Hilbert space and B : U L 2 () a linear, continuous operator. We are interested in the following control problem min J(u) = y y α u U ad 2 2 u u 0 2 U (.3) subject to y = G(Bu) and y(x) b(x) in. Here, U ad U denotes the set of admissible controls which is assumed to be closed and convex. Furthermore, we suppose that α > 0 and that y 0 H (), u 0 U and b W 2, () are given. In the special case U L 2 () without control constraints, i.e. U ad L 2 () the finite element analysis of problem (.3) is carried out in [6]. In the present work we extend the analysis to the case of control and pointwise state constraints. Here we use techniques which are applicable to a wider class of control problems. In particular the results of [6] are contained as a special case. From here onwards we impose the following assumption which is frequently referred to as Slater condition or interior point condition. Assumption.. ũ U ad G(Bũ) < b in. Since the state constraints form a convex set and the set of admissible controls is closed and convex it is not difficult to establish the existence of a unique solution u U ad to this problem. In order to characterize this solution we introduce the space M( ) of Radon measures which is defined as the dual space of C 0 ( ) and endowed with the norm µ M( ) = sup f C 0 ( ), f fdµ. Using [3, Theorem 5.2] we then infer (compare also [2, Theorem 2]) Theorem.2. Let u U ad denote the unique solution to (.3). Then there exist µ M( ) and p L 2 () such that with y = G(Bu) there holds pav = (y y 0 )v + vdµ v H 2 () with d a ij v xi ν j = 0 on (.4) i,j= (B p + α(u u 0 ),v u 0 ) U 0 v U ad (.5) µ 0, y(x) b(x) in and (b y)dµ = 0. (.6) 2
5 Our aim is to develop and analyze a finite element approximation of problem (.3). We start by approximating the cost functional J by a sequence of functionals J h where h is a mesh parameter related to a sequence of triangulations. The definition of J h involves only the approximation of the state equation by linear finite elements and enforces constraints on the state in the nodes of the triangulation, whereas the controls are still sought in U ad. We shall prove that the minima of J h converge in L 2 to the minimum of J as h 0 and that the states convergence strongly in H with corresponding error bounds. We thereby extend the semi discrete approach to purely control constrained control problems presented by the second author in [9] to problems with control and state constraints. To the authors knowledge only few attempts have been made to develop a finite element analysis for elliptic control problems in the presence of control and state constraints. In [4] Casas proves convergence of finite element approximations to optimal control problems for semi-linear elliptic equations with finitely many state constraints. Casas and Mateos extend these results in [5] to a less regular setting for the states and prove convergence of finite element approximations to semi-linear distributed and boundary control problems. In [] Meyer considers a fully discrete strategy to approximate an elliptic control problem with pointwise state and control constraints. He obtains the approximation order O(h 2 d/2 ǫ ) for the state in H and for the control in L 2, where d denotes the spatial dimension and ǫ > 0 can be chosen arbitrarily. His results confirm those obtained by the authors in [6] for the purely state constrained case. Let us comment on further approaches that tackle optimization problems for pdes with control and state constraints. A Lavrentiev-type regularization of problem (.3) is investigated in [3]. In this approach the state constraint y b in (.3) is replaced by the mixed constraint ǫu + y b, with ǫ > 0 denoting a regularization parameter. It turns out that the associated Lagrange multiplier µ ǫ belongs to L 2 (). Numerical analysis for this approach with emphasis on the coupling of gridsize and regularization parameter ǫ is presented by the second author and Meyer in [0]. The resulting optimization problems are solved either by interior-point methods or primal-dual active set strategies, compare [2]. The development of numerical approaches to tackle (.3) is ongoing. An excellent overview can be found in [7, 8], where also further references are given, for the latter see also [6]. The paper is organized as follows: in 2 we describe our discretization and establish bounds on the relevant discrete quantities which are uniform in the discretization parameter. These bounds are used in 3 in order to prove the following error bounds O(h 2), if d = 2, u u h U, y y h H = O(h 4), if d = 3, where u h and y h are the discrete control and state respectively. If in addition Bu W,s () we obtain u u h U, y y h H Ch 3 2 d 2s log h. Roughly speaking, the idea is to test (.5) with u h and (2.8), the discrete counterpart of (.5), with the continuous solution u. This is feasible since controls are not discretized explicitly. An important tool in the analysis is the use of L error estimates for finite element approximations of the Neumann problem developed in [4]. The need for uniform estimates is due to the presence of the measure µ in (.4). 3
6 2 Finite element discretization Let T h be a triangulation of with maximum mesh size h := max T Th diam(t) and vertices x,...,x m. We suppose that is the union of the elements of T h so that element edges lying on the boundary are curved. In addition, we assume that the triangulation is quasi-uniform in the sense that there exists a constant κ > 0 (independent of h) such that each T T h is contained in a ball of radius κ h and contains a ball of radius κh. Let us define the space of linear finite elements, X h := {v h C 0 ( ) v h is a linear polynomial on each T T h } with the appropriate modification for boundary elements. In what follows it is convenient to introduce a discrete approximation of the operator G. For a given function v L 2 () we denote by z h = G h (v) X h the solution of the discrete Neumann problem a(z h,v h ) = vv h for all v h X h. It is well known that for all v L 2 () G(v) G h (v) Ch 2 v, (2.) G(v) G h (v) L Ch 2 d 2 v. (2.2) The estimate (2.2) can be improved provided one strengthens the assumption on v. Lemma 2.. Suppose that v W,s () for some < s < d d. Then G(v) G h (v) L Ch 3 d s log h v W,s. Proof. Let z = G(v), z h = G h (v). Elliptic regularity theory implies that z W 3,s () from which we infer that z W 2,q () with q = ds d s using a well known embedding theorem. Furthermore, we have z W 2,q c z W 3,s c v W,s. (2.3) Using Theorem 2.2 and the following Remark in [4] we have z z h L c log h inf χ X h z χ L, (2.4) which, combined with a well known interpolation estimate, yields z z h L ch 2 d q log h z W 2,q ch 3 d s log h v W,s in view (2.3) and the relation between s and q. Problem (.3) is now approximated by the following sequence of control problems depending on the mesh parameter h: min J h (u) := y h y α u U ad 2 2 u u 0,h 2 U (2.5) subject to y h = G h (Bu) and y h (x j ) b(x j ) for j =,...,m. Here, u 0,h denotes an approximation to u 0 which is assumed to satisfy u 0 u 0,h Ch. (2.6) Problem (2.5) represents a convex infinite-dimensional optimization problem of similar structure as problem (.3), but with only finitely many equality and inequality constraints for the state, which form a convex admissible set. Again we can apply [3, Theorem 5.2] which together with [2, Corollary ] yields (compare also the analysis of problem (P) in [4]) 4
7 Lemma 2.2. Problem (2.5) has a unique solution u h U ad. There exist µ,...,µ m Ê and p h X h such that with y h = G h (Bu h ) and µ h = m j= µ jδ xj we have a(v h,p h ) = (y h y 0 )v h + v h dµ h v h X h, (2.7) (B p h + α(u h u 0,h ),v u h ) U 0 v U ad, (2.8) ( ) µ j 0, y h (x j ) b(x j ),j =,...,m and Ih b y h dµh = 0. (2.9) Here, δ x denotes the Dirac measure concentrated at x and I h is the usual Lagrange interpolation operator. Remark 2.3. Problem (2.5) is still an infinite dimensional optimization problem, but with finitely many state constraints. This is reflected by the well known fact that the variational inequalities (.5) and (2.8) can be rewritten in the form u = Π Uad ( α B p + u 0 ) and u h = Π Uad ( α B p h + u 0,h ) (2.0) respectively, where Π Uad : U U ad denotes the orthogonal projection onto U ad. Due to the presence of Π Uad the function u h will in general not belong to X h even in the case U = L 2 (),B = Id. This is different for the purely state constrained problem, for which Π Uad Id, so that u h = α p h + u 0,h X h by (2.0). In that case the space U = L 2 () in (2.5) may be replaced by X h to obtain the same discrete solution u h, which results in a finite dimensional discrete optimization problem instead. However, we emphasize, that the infinite dimensional formulation of (2.5) is crucial for our numerical analysis in 3. As a first result for (2.5) we prove that the sequence of optimal controls, states and the measures µ h are uniformly bounded. Lemma 2.4. Let u h U ad be the optimal solution of (2.5) with corresponding state y h X h and adjoint variables p h X h and µ h M( ). Then there exists h > 0 so that y h, u h U, µ h M( ) C for all 0 < h h. Proof. Since G(Bũ) is continuous, Assumption. implies that there exists δ > 0 such that It follows from (2.2) that there is h 0 > 0 with G(Bũ) b δ in. (2.) G h (Bũ) b in for all 0 < h h 0 so that J h (u h ) J h (ũ) C uniformly in h giving u h U, y h C for all h h 0. (2.2) Next, let u denote the unique solution to problem (.3). We infer from (2.) and (2.2) that v := 2 u + 2ũ satisfies G h (Bv) 2 G(Bu) + 2 G(Bũ) + Ch2 d 2( Bu + Bũ ) (2.3) b δ 2 + Ch2 d 2 ( u U + ũ U ) b δ 4 in 5
8 provided that h h, h h 0. Since v U ad, (2.8), (2.7), (2.2) and (2.3) imply 0 (B p h + α(u h u 0,h ),v u h ) U = B(v u h )p h + α(u h u 0,h,v u h ) U = a(g h (Bv) y h,p h ) + α(u h u 0,h,v u h ) U = (G h (Bv) y h )(y h y 0 ) + (G h (Bv) y h )dµ h + α(u h u 0,h,v u h ) U m ( C + µ j b(xj ) δ 4 y h(x j ) ) = C δ m µ j 4 j= where the last equality is a consequence of (2.9). It follows that and the lemma is proved. 3 Error analysis µ h M( ) C An important ingredient in our analysis is an error bound for a solution of a Neumann problem with a measure valued right hand side. Let A as above and consider j= d i= A q = µ in ( d j= a ijq xj + b i q ) ν i = µ on. (3.4) Theorem 3.. Let µ M( ). Then there exists a unique weak solution q L 2 () of (3.4), i.e. d qav = vd µ v H 2 () with a ij v xi ν j = 0 on. Furthermore, q belongs to W,s d () for all s (, d ). For the finite element approximation q h X h of q defined by a(v h,q h ) = v h d µ for all v h X h. the following error estimate holds: i,j= q q h Ch 2 d 2 µ M( ). (3.5) Proof. A corresponding result is proved in [] for the case of an operator A without transport term subject to Dirichlet conditions, but the arguments can be adapted to our situation. We omit the details. We are now prepared to prove our main theorem for the optimal controls. Theorem 3.2. Let u and u h be the solutions of (.3) and (2.5) respectively. Then u u h U + y y h H Ch d 4. If in addition Bu W,s () for some s (, d d ) then u u h U + y y h H Ch 3 2 d 2s log h. 6
9 Proof. We test (.5) with u h, (2.8) with u and add the resulting inequalities. This gives which in turn yields α u u h 2 U (B (p p h ) α(u 0 u 0,h ) + α(u u h ),u h u) U 0, B(u h u)(p p h ) α (u 0 u 0,h,u h u) U. (3.6) Let y h := G h (Bu) X h and denote by p h X h the unique solution of a(w h,p h ) = (y y 0 )w h + w h dµ for all w h X h. Applying Theorem 3. with µ = (y y 0 )dx + µ we infer p p h Ch 2 d 2( y y0 + µ M( )). (3.7) Recalling that y h = G h (Bu h ),y h = G h (Bu) and observing (2.7) as well as the definition of p h we can rewrite the first term in (3.6) B(u h u)(p p h ) = B(u h u)(p p h ) + B(u h u)(p h p h ) = B(u h u)(p p h ) + a(y h y h,p h p h ) (3.8) = B(u h u)(p p h ) + (y y h )(y h y h ) + (y h y h )dµ (y h y h )dµ h = B(u h u)(p p h ) y y h 2 + (y y h )(y y h ) + (y h y h )dµ + (y h y h )dµ h. After inserting (3.8) into (3.6) and using Young s inequality we obtain in view of (3.7), (2.) and (2.6) α 2 u u h 2 U + 2 y y h 2 (3.9) C ( p p h 2 + y y h 2 + u 0 u 0,h 2) + (y h y h )dµ + (y h y h )dµ h Ch 4 d + (y h y h )dµ + (y h y h )dµ h. It remains to estimate the integrals involving the measures µ and µ h. Since y h y h (I h b b) + (b y) + (y y h ) in we deduce with the help of (.6) (y h y h )dµ µ M( ) ( I h b b + y y h ). Similarly, (2.9) implies ( ) (y h y h )dµ h µ h M( ) b I h b + y y h. 7
10 Inserting the above estimates into (3.9) and using Lemma 2.4 as well as an interpolation estimate we infer u u h 2 U + y y h 2 Ch 4 d + C y y h L. (3.20) The estimates on u u h U now follow from (2.2) and Lemma 2. respectively. Finally, in order to bound y y h H we note that a(y y h,v h ) = B(u u h )v h for all v h X h, from which one derives the desired estimates using standard finite element techniques and the bounds on u u h U. Remark 3.3. Let us note that the approximation order of the controls and states in the presence of control and state constraints is the same as in the purely state constrained case, if Bu W,s (). This assumption holds for the important example U = L 2 (), B = Id and u 0h = P h u 0, with u 0 H () and P h : L 2 () X h denoting the L 2 projection, and subsets of the form U ad = {v L 2 (),a l v a u a.e. in }, with bounds a l,a u W,s (), since u 0 H (), and p W,s (). 4 Implementation issues The numerical computation of solutions to problem (2.5) is more involved than in the purely state constrained case, i.e. U ad U. The latter is treated in [6] for the special case U = L 2 (). As already mentioned in Remark 2.3 the purely state constrained problem may be substituted by a finite dimensional one which yields the same solution u h. Therefore, common solution techniques for finite dimensional optimization problems with equality and inequality constraints can be applied for its numerical computation. In the present situation this is definitely different, since problem (2.5) really is infinite dimensional, with finitely many state constraints. It is not clear at the first instance, whether an iterative solution algorithm for this problem can be implemented on a computer without further discretization steps, while keeping the management of information overhead bounded with increasing number of iterations. We now sketch that this is indeed possible. To simplify the exposition we now assume a ij = δ ij, B = Id, U = L 2 () and u 0 = 0,y 0 = 0. We rewrite (2.5) as nonlinear, nonsmooth operator equation. To begin with, we define m m m A = Ã + M := φ i φ j + φ i φ j U := uφ i B := [b(x i )] m i=, and i,j=, i,j=, Y := [y i ] m i=, P := [p i ] m i=, M := [µ i ] m i=, i=, where y h p h µ h = m i= y i φ i p i φ i µ i δ xi 8
11 denote the corresponding finite element representations of Y, P and M, respectively, and u L 2 (). We recall that (2.8) in the present situation is equivalent to the nonsmooth equation u h = P Uad ( ) α p h, and that the complementarity system in (2.9) can be equivalently rewritten in the form M = max {0,M + c(y B)} component wise, (4.2) where c > 0 is arbitrary. Let us define the mapping G : Ê m Ê m L 2 () Ê m Ê m Ê m L 2 () Ê m by AY U AP MY M G(Y,P,u,M) := ( u P Uad α p ) h M max {0,M + c(y B)}. Then, the first order necessary optimality conditions y h = G h (u) together with (2.7) (2.9) are equivalent to the nonsmooth system G(Y,P,u h,m) = 0, where u h denotes the solution of (2.5), and Y,P, and M the nodal representations of the corresponding state y h, adjoint state p h, and multiplier µ h, respectively. Let us apply a generalized Newton step to this system. For this purpose let [Y,P,u,M] t Ê m Ê m L 2 () Ê m Ê m Ê m L 2 () Ê m be given. Let us define fmax := max (0,M + c(y B)) and dmax := diag (max (0,M + c(y B))). For the generalized Jacobian of G at [Y,P,u,M] t we obtain DG(Y,P,u,M) = A 0 [ ] m φ i i= 0 M A 0 Id M = 0 α P U ( ad α p h) IdL 2 () 0 cdmaxid Y 0 0 Id M dmaxid M so that one generalized Newton step for the computation of [Y n,p n,u n,m n ] t amounts to solving Y n 0 P n 0 DG(Y,P,u,M) = u n α P U ( ad α p ( h) ph + P Uad α p ) h M n cdmaxy dmaxm + fmax From the first equation we deduce Y n = A U n, the second yields M n = AP n MA U n, and from the third equation we obtain U n = α CP n +R, where the matrix C and the vector R are defined by C = C(p h ) := φ j φ i I( α p h) m i,j= and R = α 9 p h φ i I( α p h) m i= + P Uad ( α p h ) φ i, m i=.
12 Here, I ( α p h) denotes the subset of, where α p h is inactive. Its description clearly relies on the structure of U ad. Finally, with the substitutions from above the fourth equation leads to ( A dmaxa + ( ( M + dmax cid Y α M )) ) A C P n = = MA R + ( cid Y α M ) A R dmaxy dmaxm + fmax, (4.22) which represents a linear system for the computation of P n. We conclude, that a step of the generalized Newton method is feasible, iff the matrix A dmaxa + α ( ( M + dmax cid Y M )) A C is regular. In this case, once P n is computed, the vectors Y n,m n and the function u n are obtained by performing the corresponding re substitutions. Let us briefly comment on the cases i.) U ad U and ii.) b =, i.e. no state constraints are present. In the first case the system matrix in (4.22) takes the form A dmaxa + α ( ( M + dmax cid Y M )) A M and is positive definite if we choose c large enough. In the second case we obtain the system matrix A + α MA C which is positive definite since MA α C is positive semi definite. In both cases the Newton step is feasible. This reflects the fact that in these cases the adjoint variables p,p h and the multipliers µ,µ h are unique. Let us finally note that the constant c in (4.2) may also be replaced by an vector c Ê m containing positive components c i > 0(i =,...,m). If one now replaces the mass matrix M by its mass lumped version M, the choice c := diag( M) leads to the system matrix ( A dmaxa + ) α MA C in (4.22). We note that the adjoint states p,p h and the multipliers µ,µ h need not be unique if state and control constraints occur simultaneously. Necessary and sufficient conditions for the uniqueness of multipliers in the presence of cone constraints are provided by Shapiro in [5]. The exposition in this section shows that the infinite dimensional problem (2.5) can numerically be implemented on a computer by a generalized Newton method, say, without further discretization steps, and with keeping the management of information overhead constant in each iteration step. Further details inc. convergence analysis for generalized Newton methods applied to the numerical solution of (2.5) will be given in a forthcoming paper. Acknowledgements The authors gratefully acknowledge the support of the DFG Priority Program 253 entitled Optimization With Partial Differential Equations. 0
13 References [] Casas, E.: L 2 estimates for the finite element method for the Dirichlet problem with singular data, Numer. Math. 47, (985). [2] Casas, E.: Control of an elliptic problem with pointwise state constraints, SIAM J. Cont. Optim. 24, (986). [3] Casas, E.: Boundary control of semilinear elliptic equations with pointwise state constraints, SIAM J. Cont. Optim. 3, (993). [4] Casas, E.: Error Estimates for the Numerical Approximation of Semilinear Elliptic Control Problems with Finitely Many State Constraints, ESAIM, Control Optim. Calc. Var. 8, (2002). [5] Casas, E., Mateos, M.: Uniform convergence of the FEM. Applications to state constrained control problems. Comp. Appl. Math. 2 (2002). [6] Deckelnick, K., Hinze, M.: Convergence of a finite element approximation to a state constrained elliptic control problem, MATH-NM , Institut fr Numerischen Mathematik, TU Dresden (2006). [7] Hintermüller, M., K., Kunisch, K.: Path following methods for a class of constrained minimization methods in function spaces, Report RICAM , RICAM Linz (2004). [8] Hintermüller, M., K., Kunisch, K.: Feasible and non-interior path following in constrained minimization with low multiplier regularity, Report, Universität Graz (2005). [9] Hinze, M.: A variational discretization concept in control constrained optimization: the linear-quadratic case, J. Computational Optimization and Applications 30, (2005). [0] Hinze, M., Meyer, C.: Numerical analysis of Lavrentiev-regularized state constrained elliptic control problems, WIAS Preprint 68, Berlin (2006). [] Meyer, C.: Error estimates for the finite element approximation of an elliptic control problem with pointwise constraints on the state and the control, WIAS Preprint 59 (2006). [2] Meyer, C., Prüfert, U., Tröltzsch, F.: On two numerical methods for state-constrained elliptic control problems, Technical Report , Institut für Mathematik, TU Berlin (2005). [3] Meyer, C., Rösch, A., Tröltzsch, F.: Optimal control problems of PDEs with regularized pointwise state constraints, Preprint 4, Institut für Mathematik, TU Berlin, to appear in Computational Optimization and Applications (2004). [4] Schatz, A.H.: Pointwise error estimates and asymptotic error expansion inequalities for the finite element method on irregular grids. I: Global estimates, Math. Comput. 67, No.223, (998). [5] Shapiro, A.: On uniqueness of Lagrange multipliers in optimization problems subject to cone constraints, SIAM J. Optim. 7, (997). [6] Tröltzsch, F.: Optimale Steuerung mit partiellen Differentialgleichungen. Vieweg, 2005.
Hamburger 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 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 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 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 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 informationBang bang control of elliptic and parabolic PDEs
1/26 Bang bang control of elliptic and parabolic PDEs Michael Hinze (joint work with Nikolaus von Daniels & Klaus Deckelnick) Jackson, July 24, 2018 2/26 Model problem (P) min J(y, u) = 1 u U ad,y Y ad
More 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
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 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 informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More 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 informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
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 informationA DUALITY APPROXIMATION OF SOME NONLINEAR PDE s
ISSN 2066-6594 Ann. Acad. Rom. Sci. Ser. Math. Appl. Vol. 8, No. /206 A DUALITY APPROXIMATION OF SOME NONLINEAR PDE s Dan Tiba Abstract We discuss a discretization approach for the p - Laplacian equation
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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 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 informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationInterpolation in h-version finite element spaces
Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
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 informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
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 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 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 informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationProper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints
Proper Orthogonal Decomposition for Optimal Control Problems with Mixed Control-State Constraints Technische Universität Berlin Martin Gubisch, Stefan Volkwein University of Konstanz March, 3 Martin Gubisch,
More 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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 153 Optimization with Partial Differential Equations T. Apel, J. Pfefferer and A. Rösch Finite Element Error Estimates for Neumann Boundary Control Problems
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 informationOptimal Control of Partial Differential Equations I+II
About the lecture: Optimal Control of Partial Differential Equations I+II Prof. Dr. H. J. Pesch Winter semester 2011/12, summer semester 2012 (seminar), winter semester 2012/13 University of Bayreuth Preliminary
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationLecture Notes of the Autumn School Modelling and Optimization with Partial Differential Equations Hamburg, September 26-30, 2005
Lecture Notes of the Autumn School Modelling and Optimization with Partial Differential Equations Hamburg, September 26-30, 2005 Michael Hinze, René Pinnau, Michael Ulbrich, and Stefan Ulbrich supported
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationA Mesh-Independence Principle for Quadratic Penalties Applied to Semilinear Elliptic Boundary Control
Schedae Informaticae Vol. 21 (2012): 9 26 doi: 10.4467/20838476SI.12.001.0811 A Mesh-Independence Principle for Quadratic Penalties Applied to Semilinear Elliptic Boundary Control Christian Grossmann 1,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
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 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 informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More information1 Constrained Optimization
1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y
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 informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More 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 informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationEXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS
Electronic Journal of Differential Equations, Vol. 2008(2008), No. 98, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) EXISTENCE
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationA POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS RONALD H.W. HOPPE AND MICHAEL KIEWEG Abstract. We provide an a posteriori error
More 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 informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
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 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 informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More 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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Space-time finite element approximation of parabolic optimal control problems Wei Gong, Michael Hinze, and Zhaojie Zhou Nr. 2-22 November 2 SPACE-IME FINIE
More information