have demonstrated that an augmented Lagrangian techniques combined with a SQP approach lead to first order conditions and provide an efficient numeric
|
|
- Richard Cunningham
- 5 years ago
- Views:
Transcription
1 Optimization Techniques for Solving Elliptic Control Problems with Control and State Constraints. Part : Boundary Control HELMUT MAURER Westfalische Wilhelms-Universitat Munster, Institut fur Numerische und instrumentelle Mathematik, Einsteinstrasse 6, 4849 Munster, Germany HANS D. MITTELMANN Arizona State University, Department of Mathematics, Tempe, A , USA Abstract: We study optimal control problems for semilinear elliptic equations subject to control and state inequality constraints. In a rst part we consider boundary control problems with either Dirichlet or Neumann conditions. By introducing suitable discretization schemes, the control problem is transcribed into a nonlinear programming problem. It is shown that a recently developed interior point method is able to solve theses problems even for high discretizations. Several numerical examples with Dirichlet and Neumann boundary conditions are provided that illustrate the performance of the algorithm for dierent types of controls including bang{bang and singular controls. The necessary conditions of optimality are checked numerically in the presence of active control and state constraints. Key words: Elliptic control problems, boundary control, control and state constraints, discretization techniques, interior point optimization methods Introduction This is the first part of two papers in which we develop nonlinear programming techniques for solving elliptic control problems under general control and state inequality constraints. In the rst part, we study boundary control problems with boundary conditions of either Dirichlet or Neumann type. The second part is devoted to elliptic problems with distributed control. There are several recent papers on both the theoretical treatment and numerical solution methods for elliptic control problems. First order necessary and second order sufficient optimality conditions for Neumann boundary conditions have been given in Bonnans [6], Casas [9] and Casas et al. [, ]. First order necessary conditions for linear operators and Dirichlet or Neumann boundary conditions are obtained in Bergounioux, Kunisch [3], Ito, Kunisch [6], Kunisch, Volkwein [8]. These authors Supported by the Deutsche Forschungsgemeinschaft (DFG) under MA 69/8-
2 have demonstrated that an augmented Lagrangian techniques combined with a SQP approach lead to first order conditions and provide an efficient numerical algorithm. Despite this work on elliptic problems, we feel it worthwile to consider these problems from a more systematic numerical point of view. We treat semilinear elliptic operators and concentrate on handling possibly nonlinear control and state constraints jointly. Our numerical approach will be able to capture also controls of bang{bang or singular type for which the Legendre condition is not satisfied. This type of control is well studied for ODE control problems, but we are not aware of any numerical example for elliptic problems although Hettich et al. [4, 5] present some theoretical work on the subject. Moreover, we analyze adjoint variables corresponding to equality and inequality constraints in the discretized problem. This enables us to check first order necessary conditions explicitly in the presence of active control and state constraints. As a byproduct, we give a informal form of first order necessary conditions for problems with Dirichlet boundary control. Such conditions have not been given in the literature to full extent so far. In the application of NLP{techniques to optimal control, there are two components that have been extensively worked out for ODE control problems; cf., e.g., Barclay et al. [], Betts [4], Betts, Humann [5], Buskens [7], Buskens, Maurer [8], Grachev, Evtushenko [3], Teo et al. []. The first aspect concerns the suitable choice of a discretization scheme while the other is the selection of the NLP{method. One has two options for the discretization scheme. The first one is to discretize both the control and state variables and to incorporate the integration method as an explicit equality constraint at each gridpoint. This approach leads to a high{ dimensional NLP{problem with a sparse structure of the Jacobian; cf. Barclay et al. []. The other discretization approach consists in treating the discretized control variables as the only optimization variables while the state variable is expressed and computed as a function of the control variable. This leads to a lower-dimensional NLP with a rather dense Jacobian. However, in this approach derivatives usually can not be calculated explicitly but only through a numerical differentiation scheme. In this paper, we formulate NLP{problems using a full discretization scheme where the optimization variables comprise both the control and state variables. The resulting NLP{problems may contain up to 4; variables. To solve such a high{ dimensional and sparse NLP{problem, the interior point method developed by Vanderbei, Shanno [3] has turned out to be particularly efficient and reliable. The organization of the rst part is the following. In Section we discuss necessary optimality conditions for elliptic problems with Neumann boundary control. Necessary optimality conditions for problems with Dirichlet boundary conditions have not yet been developed in the literature for the general problem considered here. In Section 3 we state an informal form of such necessary conditions. Section 4 formulates NLP{problems associated with discretized versions of the elliptic problems. Necessary conditions of Kuhn{Tucker type are discussed both for Dirichlet and Neumann boundary conditions. Finally, in Section 5 we present several numerical
3 examples for both types of boundary conditions. Example 5. exhibits a singular control while Examples 5.4, 5.6 and 5.8 present bang{bang controls. Elliptic control problems with Neumann boundary conditions The following elliptic control problem with control and state constraints constitutes a generalization of elliptic problems considered in Bonnans [6], Casas [9], Casas et al.[, ], Ito, Kunisch [6], Kunisch, Volkwein [8]. The problem is to minimize the functional F (y; u) = f(x; y(x)) dx + subject to the state equation? g(x; y(x); u(x)) dx (.)?y(x) + d(x; y(x)) = ; for x y(x) = b(x; y(x); u(x)) ; for x? ; (.) and the inequality constraints on control and state C(x; y(x); u(x)) ; for x? ; (.3) S(x; y(x)) ; for x : (.4) In this setting, IR is a bounded domain with piecewise smooth boundary? The derivative in the direction of the outward unit normal of? is denoted in (.). Note that the state inequality constraints (.4) are supposed to hold on the closure of. The Laplacian? in (.) can be replaced by an elliptic operator Ay(x) =? xk (a kj xj y)(x) ; where the coefficients a kj C () satisfy the following coercivity condition with some c > : k;j= a kj (x)v k v j c(v + v ) 8x ; v IR : However, in the sequel we restrict the discussion to the operator A =? which simplifies the form of the necessary conditions. The functions f : IR! IR; g :?IR! IR; d : IR! IR; b :?IR! IR; C :?IR! IR and S : IR! IR are assumed to be C -functions. It is straightforward to include more than one inequality constraint into (.3) or (.4). However, since both the state and control 3
4 variable are scalar variables, the active sets for different inequality constraints are disjoint and hence can be treated separately. Then under appropriate assumptions on the function d it can be shown that the state equation (.) admits for each u L (?) a weak solution y Y = C() \ H () (cf. Casas et al. []), i.e., it holds [ y(x) v(x) + d(x; y(x))v(x) ] dx =? b(x; y(x); u(x)) v(x) dx for all v H() : An optimal solution of problem (.){(.4) will be denoted by u and y. From [] we infer the further assumption that the function b in the Neumann condition (.) is sufficiently smooth and satisfies the following inequality with suitable >, b y (x; y; u) for all x?; j y? y(x) j < ; j u? u(x) j < : (.5) Questions of existence of optimal solutions will not be discussed here. The active sets for the inequality constraints (.3) and (.4) are defined by A(C) := f x? j C(x; y(x); u(x)) = g ; A(S) := f x j S(x; y(x)) = g : It is required that the following regularity conditions hold: C u (x; y(x); u(x)) 6= 8 x A(C) ; S y (x; y(x)) 6= 8 x A(S) : (.6) Here and in the following, partial derivatives are denoted by subscripts. First order optimality conditions for a local optimal solution u and y can be derived by generalizing the line of proof in Casas [9], Casas et al. [, ]. Problem (.){(.4) is considered as a mathematical programming problem in Banach spaces to which the first order Kuhn{Tucker conditions are applicable. In particular, this approach requires that the regularity condition given in owe, Kurcyusz [4] is satisfied; cf. Casas et al. []. We do not discuss this regularity condition in detail although condition (.6) forms part of it. The first order necessary conditions imply that there exist an adjoint state q W ; (), a multiplier L (?), and a function of bounded variation BV ( ) such that the following three conditions hold,. adjoint equation:?q(x) + q(x) d y (x; y(x)) + f y (x; y(x)) + d(x)s y (x; y(x)) = on ; q(x)? q(x) b y (x; y(x); u(x)) + g y (x; y(x); u(x)) + +d(x)s y (x; y(x)) + (x)c y (x; y(x); u(x)) = on? ; (.8). minimum condition on? : g u (x; y(x); u(x))? q(x) b u (x; y(x); u(x)) + (x) C u (x; y(x); u(x)) = ; (.9) 4
5 3. complementarity condition: (x) on A(C) ; (x) = on? n A(C) d(x) on ; (x) = constant on n A(S) : (.) The adjoint equations (.7), (.8) are understood in the weak sense, cf. Casas et al. []. In later applications, we shall mostly deal with cost functionals of tracking type (cf. Ito, Kunisch [6]), F (y; u) = (y(x)? y d (x)) dx +? (u(x)? u d (x)) dx (.) with given functions y d C(); u d L (?), and a nonnegative weight : The control and state constraints (.3) and (.4) are supposed to be box constraints of the simple type y(x) (x) on ; u (x) u(x) u (x) on? ; (.) with functions C( ) and u ; u L (?). For these data the adjoint equation (.7), (.8) become?q(x) + q(x) d y (x; y(x)) + y(x)? y d (x) + d(x) = on q(x)? q(x)b y (x; y(x); u(x)) + (x)c y (x; y(x); u(x)) = on? ; (.3) If the function b in (.) has the form b(x; y; u) = u + b (x; y), i.e., if b u holds, then the minimum condition (.9) reduces to [(u(x)? u d (x))? q(x)] (u? u(x)) 8 x?; u [u (x); u (x)] : (.4) More explicitly, for > this condition yields the control law: u(x) = 8 >< >: u d (x) + q(x)= ; if u d (x) + q(x)= (u (x); u (x)) ; u (x) ; if u d (x) + q(x)= u (x) ; u (x) ; if u d (x) + q(x)= u (x) : 9 >= >; (.5) In the case = we obtain an optimal control of bang{bang or singular type: u(x) = 8 >< >: u (x) ; if q(x) < ; u (x) ; if q(x) > ; singular ; if q(x) = on? s? ; meas(? s ) > : 9 >= >; (.6) Thus, for = the adjoint function q(x) on the boundary plays the role of a switching function. The isolated zeros of q(x)j? are the switching points of a bang{ bang control. 5
6 3 Elliptic control problems with Dirichlet boundary conditions The following elliptic control problem with control and state constraints generalizes the elliptic problems considered in Bergounioux, Kunisch [3], Ito, Kunisch [6]. When admitting Dirichlet boundary conditions, the functions g in the cost functional (.), b in the boundary conditions (.), and C in the constraints (.3) do not depend on the state variable y. This leads to the problem of minimizing the functional F (y; u) = f(x; y(x)) dx + subject to the state equation? g(x; u(x)) dx (3.)?y(x) + d(x; y(x)) = ; for x ; y(x) = b(x; u(x)) ; for x? ; (3.) and the inequality constraints on control and state C(x; u(x)) ; for x? ; (3.3) S(x; y(x)) ; for x : (3.4) When treating Dirichlet boundary conditions, an intrinsic difficulty arises from the fact that the solution operator u()j?! y() for (3.) is not suciently smooth as to give an appropriate form of first order necessary conditions: cf., Lions [9], Lions, Magenes [], Chapter. A weak formulation of rst order necessary conditions for linear elliptic equations may be found in Bergounioux, Kunisch [3]. Instead of trying to prove rst order necessary conditions under strong assumptions we content ourselves with deriving first order necessary conditions in a purely formal way. This form of the necessary conditions will be justified by its analogy in the rst order necessary conditions for the discretized version of the elliptic problem; cf. Section 4.. Denote an optimal solution of problem (3.){(3.4) by u and y. The active sets of inequality constraints (3.3) and (3.4) are A(C) := f x? j C(x; y(x); u(x)) = g ; A(S) := f x j S(x; y(x)) = g : We require the following regularity conditions: C u (x; u(x)) 6= on A(C) ; S y (x; y(x)) 6= on A(S) : (3.5) Let q and p be multipliers associated with the elliptic equation and the Dirichlet boundary condition in (3.), and let and be multipliers associated with the 6
7 control and state inequality constraints (3.3) and (3.4). Then the Lagrangian for problem (3.){(3.4) becomes L(y; u; q; p; ; ) : = + +? f(x; y(x)) dx +? g(x; u(x)) dx [?y(x) + d(x; y(x) ] q(x) dx + [ y(x)? b(x; u(x)) ] p(x)dx +? d(x)s(x; y(x)) (x)c(x; u(x)) dx: The rst order necessary condition with respect to the state variable gives L y (y; u; q; p; ; )y = for all y : Function evaluation along a stationary solution (y; u; q; p; ; ) will be denoted henceforth by a bar, e.g., f = f(x; y); etc. Using partial integration (Green's theorem), we rewrite L y y = as L y y = (?y q + y q) dx = [?q + q d y + f y + d S y ] y dx +?? y + q) dx ; q? p ] y dx?? y dx = which holds for all functions y. From this we obtain the adjoint equations :?q(x) + q(x) d y (x; y(x)) + f y (x; y(x)) + d(x)s y (x; y(x)) = on ; (3.6) q(x) = on? : (3.7) Moreover, we find the multiplier p =?@ q on?. Then the optimality condition for the control variable is evaluated as which yields the minimum condition: L u u = [ g u + (@ q) b u + C u ]u = for all u ; g u (x; u(x)) q(x) b u (x; u(x)) + (x) C u (x; u(x)) = on? : (3.8) The complementarity condition is similar to (.): (x) on A(C) ; (x) = on? n A(C) d(x) on ; (x) = constant on n A(S) : (3.9) Note that the minimum condition (3.8) agrees with the minimum (.9) if we replace q formally by?@ q. Hence, assuming functionals of tracking type (.), box constraints (.) for control and state and the property b u, we obtain the 7
8 counterparts of the control laws (.5), (.6): for > the control is determined by u(x) = 8 >< >: u d q(x)= ; if u d q(x)= (u (x); u (x)) ; u (x) ; if u d q(x)= u (x) ; u (x) ; if u d q(x)= u (x) ; while for = the optimal control is of bang{bang or singular type u(x) = 8 >< >: u (x) ; q(x) > ; u (x) ; q(x) < ; singular ; q(x) = on? s? ; meas(? s ) > : 9 >= >; 9 >= >; (3.) (3.) Thus, for = the outward normal q(x) on the boundary plays the role of a switching function. The isolated zeros q(x)j? are the switching points of a bang{bang control. 4 Discretization and optimization techniques The discussion of discretization schemes is restricted to the standard situation where the elliptic operator is the Laplacian A =? and the domain is the unit square = (; )(; ) : The generalization to a general elliptic operator is straightforward. However, the modifications for an arbitrary domain depend essentially on the geometry of the boundary?. The purpose of this section is to develop discretization techniques by which the problem (.){(.4) with Neumann boundary conditions resp. problem (3.){(3.4) with Dirichlet boundary conditions is transformed into a nonlinear programming problem (NLP-problem) of the form Minimize F h (z) subject to G h (z) = ; H(z) : (4.) The functions F h ; G h and H are suciently smooth and of appropriate dimension. The upper subscript h denotes the dependence on the stepsize. The optimization variable z will comprise both the state and the control variables. The form (4.) will be achieved by solving the elliptic equation (.) with the standard ve-point-star discretization scheme. Choose a number N IN + and the stepsize h := =(N + ) : Consider the mesh points x ij = (ih; jh) ; i; j N + ; and consider the following sets of indices (i; j) residing either in the domain or on the four edges of the boundary? : I() := f (i; j) j i; j N g ; I(?) := f (i; j) j i = ; :::; N; j = or j = N; j = ; :::; N; i = or i = N g ; I() := I() [ I(?) : 8 (4.)
9 We have #I() = N and #I(?) = 4 N. Denote approximations for the values y(x ij ) of the state variables by y ij for (i; j) I() and denote approximations of the values u(x ij ) of the control variables by u ij for (i; j) I(?). Now we shall specify the functions F h ; G h ; H for the optimization problem (4.) both for Neumann and Dirichlet boundary conditions. 4. Neumann boundary conditions The optimization variable z in (4.) is taken as the vector z := ( (y ij ) (i;j)i( ) ; (u ij ) (i;j)i(?) ) IR N +4N +4N : Equality constraints are obtained by applying the ve{point{star to the elliptic equation?y(x) + d(x; y(x)) = in (.) in all points x ij with (i; j) I(): G h ij(z) := 4y ij? y i+;j? y i?;j? y i;j+? y i;j? + h d(x ij ; y ij ) = : (4.3) The y(x ij ) in the direction of the outward normal is approximated by the expression yij=h where y ij := 8 >< >: y i? y i ; for j = ; i = ; :::; N y j? y j ; for i = ; j = ; :::; N y Nj? y N?;j ; for i = N; j = ; :::; N y in? y i;n? ; for j = N; i = ; :::; N 9 >= >; : (4.4) Then the discrete form of the Neumann boundary condition in (.) leads to the equality constraints B h (z) := y ij? h b(x ij ; y ij ; u ij ) = for (i; j) I(?) : (4.5) The control and state inequality constraints (.3) and (.4) yield the inequality constraints S(x ij ; y ij ) ; (i; j) I() ; (4.6) C(x ij ; y ij ; u ij ) ; (i; j) I(?) : (4.7) Observe that these inequality constraints do not depend on the meshsize h. Later on, this fact will require a scaling of the Lagrange multipliers. The discretized form of the cost function (.) is F h (z) := h (i;j)i() f(x ij ; y ij ) + h (i;j)i(?) g(x ij ; y ij ; u ij ) : (4.8) 9
10 Then the relations (4.3){(4.8) define an NLP{problem of the form (4.). The Lagrangian function for this NLP{problem becomes L(z; q; ; ) := h (i;j)i() + (i;j)i() + (i;j)i(?) f(x ij ; y ij ) + h q ij G h ij(z) + (i;j)i(?) (i;j)i( ) g(x ij ; y ij ; u ij ) (4.9) ij S(x ij ; y ij ) [ q ij B h (z) + ij C(x ij ; y ij ; u ij ) ] : The Lagrange multipliers q = (q ij ) (i;j)i( ), = ( ij ) (i;j)i( ) resp. = ( ij ) (i;j)i(?) are associated with the equality constraints (4.3) and (4.5), the inequality constraints (4.6), resp. the inequality constraints (4.7). The multipliers and satisfy complementarity conditions corresponding to (.): ij and ij C(x ij ; y ij ; u ij ) = for all (i; j) I(?) ; ij and ij S(x ij ; y ij ) = for all (i; j) I() : The discussion of the necessary conditions of optimality = L z = ( (L yij ) (i;j)i( ) ; (L uij ) (i;j)i(?) ) will be performed for different combinations of indices (i; j). For indices (i; j) I() we obtain the relations = L yij = 4q ij? q i+;j? q i?;j? q i;j+? q i;j? +h q ij d y (x ij ; y ij ) + h f y (x ij ; y ij ) + ij S y (x ij ; y ij ) : Hence, the Lagrange multipliers q = (q ij ) satisfy the ve{point{star dierence equations for the adjoint equation?q + q d y + f y + d S y = in (.7) if we make the identification d(x ij ) ij =h for (i; j) I() (4.) for the function of bounded variation. For indices (i; j) I(?) on the boundary we get, e.g. for j = ; i = ; :::; N : = L yi =?q i + q i? q i hb y (x i ; y i ; u i ) + hg y (x i ; y i ; u i ) + + i S y (x i ; y i ) + i C y (x i ; y i ; u i ) : This is just the discrete version of the Neumann boundary condition (.8) if we identify (x i ) i =h; d(x i ) i =h : (4.) Similar relations hold for other indices (i; j) I(?). Finally, necessary conditions with respect to the control variables u ij ; (i; j) I(?) are determined, e.g., for indices j = ; i = ; :::; N ; by = L ui = hg u (x i ; y i ; u i )? q i hb u (x i ; y i ; u i ) + i C u (x i ; y i ; u i ) : This is the discrete version of the optimality condition (.9) for the control, if we use again the identification (x i ) i =h.
11 4. Dirichlet boundary conditions The Dirichlet boundary conditions in (3.) are incorporated by the relations y ij = b(x ij ; u ij ) for all (i; j) I(?) : (4.) Hence it suffices to work with a reduced number of optimization variables in (4.), namely we can take z := ( (y ij ) (i;j)i() ; (u ij ) (i;j)i(?) ) IR N +4N : The equality constraints agree with those from (4.3), G h ij(z) := 4y ij? y i+;j? y i?;j? y i;j+? y i;j? + h d(x ij ; y ij ) = (4.3) for all indices (i; j) I() where the values y ij on the boundary are substituted from the relations (4.). The control and state inequality constraints (3.3) resp. (3.4) give rise to the inequality constraints: S(x ij ; y ij ) ; (i; j) I() ; (4.4) C(x ij ; u ij ) ; (i; j) I(?) : (4.5) Finally, the cost function is derived from (3.) as F h (z) := h (i;j)i() f(x ij ; y ij ) + h (i;j)i(?) g(x ij ; u ij ) : (4.6) By means of (4.3){(4.6) we have obtained an NLP{problem of the form (4.). Let us evaluate the first order optimality conditions of Kuhn{Tucker type for problem (4.3){(4.6). The Lagrangian function becomes L(z; q; ; ) := h (i;j)i() + (i;j)i() f(x ij ; y ij ) + h (i;j)i(?) [ q ij G h ij(z) + ij S(x ij ; y ij ) ] + g(x ij ; u ij ) (4.7) (i;j)i(?) ij C(x ij ; u ij ) ; where the Lagrange multipliers q = (q ij ) (i;j)i(), = ( ij ) (i;j)i() and = ( ij ) (i;j)i(?) which are associated with the equality constraints (4.3), the inequality constraints (4.4), and the inequality constraints (4.5). The multipliers and satisfy the complementarity conditions corresponding to (3.9): ij and ij C(x ij ; u ij ) = for all (i; j) I(?) ; ij and ij S(x ij ; y ij ) = for all (i; j) I() : In the next step we discuss the necessary conditions of optimality = L z = ( (L yij ) (i;j)i() ; (L uij ) (i;j)i(?) )
12 for different combinations of indices (i; j). For indices i; j N? we obtain the relations = L yij = 4q ij? q i+;j? q i?;j? q i;j+? q i;j? +h q ij d y (x ij ; y ij ) + h f y (x ij ; y ij ) + ij S y (x ij ; y ij ) : Hence, for this set of indices we see that the Lagrange multipliers q = (q ij ) satisfy the ve{point{star dierence equations for the adjoint equation?q+q d y +f y +d S y = in (3.6) if we make the further identification d(x ij ) ij =h (4.8) involving the function of bounded variation. The situation is different for indices with either i = ; j = ; i = N or j = N. E.g., for indices j = ; i = ; :::; N? ; we get = L yi = 4q i? q i+;? q i?;? q i; +h q i d y (x i ; y i ) + h f y (x i ; y i ) + i S y (x i ; y i ) : Note that this equation does not involve the variable q i on the boundary. Defining q i = we arrive again at the discretized form of the adjoint equation (3.6). Similar relations hold for index combinations i = ; j = ; ::; N?; or i = N; j = ; ::; N?; or j = N; i = ; ::; N?. For the remaining indices we find, e.g. for i = j = : = L y = 4q? q? q + h q d y (x ; y ) + h f y (x ; y ) + S y (x ; y ) : These relations do not involve the boundary variables q and q. Again, defining q = and q = we recover the adjoint equations. In summary, by requiring the Dirichlet boundary condition for the adjoint variables, q ij = for (i; j) I(?); we see that the adjoint equation holds for all indices (i; j) I() : The necessary conditions with respect to the control variables u i for the indices j = ; i = ; :::; N; by are determined = L ui = hg u (x i ; u i )? q i b u (x i ; u ii ) + i C u (x i ; u i ) = h [ g u (x i ; u i )? q i h b u(x i ; u ii ) + i h C u(x i ; u i ) ] Observing q i = and the approximation of the normal derivative in (4.4), the minimum condition (3.8) holds with the q(x i )?q i =h ; (xi ) i =h (4.9) Similar identifications hold for the other indices (i; j) I(?).
13 4.3 Optimization codes and modeling environment For the numerical solution of all problems considered in the following section a combination of the AMPL [] algebraic modeling language and the interior point solver LOQO [3] proved to be both convenient and powerful. In order to make the formulation of mathematical optimization problems generic and independent of both the actual solver used and the programming language it is written in, modeling languages were developed. AMPL provides interfaces to a large number of solvers, both commercial and free-for-research codes. One of the latter ones is LOQO which grew out of an interior point LP optimizer to a convex QP and very recently to a general NLP solver implementing an interior point approach. Although the code is currently still being perfected it proved to be very ecient for the solution of large-scale nonlinear problems in the benchmarks of []. It was thus chosen for the following computations; see also the comparison in Example 5. below. It should be remarked that LOQO implements an infeasible primal{dual path{following method. The KKT necessary conditions are essentially solved as a system of nonlinear equations with a Newton{like method. Therefore, it causes no problem if iterates are not feasible because it only means that residuals or right{hand sides corresponding to the equality constraints are not zero. At least asymptotically feasibility will be attained. Another feature that makes AMPL attractive and that was exploited is its automatic dierentiation capability. Only functions for objective and constraints need to be provided. 5 Numerical examples We consider elliptic problems with the following specifications: the cost functional is of tracking type (.), the elliptic operator in (.) is the Laplacian A =? on the unit square = (; ) (; ) ; and the control and state constraints are box constraints given in (.3). The choice of symmetric functions y d (x) and u d (x) in the tracking functional, implies that the optimal control is the same on every edge of? : However, we have treated the discretized controls on every edge of? as independent optimization variables. The symmetry of the optimal control will then be a result of the optimization procedure. Tests for the following examples were run with different stepsizes and starting values. For convenience, in the sequel we shall report on results obtained for fixed stepsize and starting values: N = 99 ; h = =(N + ) = = ; u ij = ; y ij = : First, we report on 4 examples with Dirichlet boundary conditions for which partial numerical results are available in the literature. 3
14 .65 u x Figure : Optimal control for Example Dirichlet boundary conditions Example 5.: To enable a comparison, we consider first the example in Bergounioux, Kunisch [3], Section 5., with the following data: on :?y(x) = ; y(x) 3:5; y d (x) = 3 + 5x (x? )x (x? ); on? : y(x) = u(x); u(x) ; u d (x) ; = :: The following results were obtained, they are explained below: Cost functional : F (y; u) = :9655, CPU seconds : 96 The optimal control is shown in Figure. The optimal state and adjoint variable are not shown here, because they are very similar to those in Figure 3 where = : is replaced by = :. It should be noted that the output of Lagrange multipliers is provided by AMPL. It is no diculty for primal{dual methods to produce multipliers since, in fact they are variables that are computed simultaneously with the primal problem variables. The control constraints are not active while the state variable attains its upper bound only in the center x := x ij = (:5; :5) of the unit square with dual variable ij = :46. This is in agreement with the results in [3]. It can be checked that the adjoint equation (3.6) holds in the discretized version with h denoting the discretized Laplacian,? h q(x ij ) + y(x ij )? y d (x ij ) = for x ij 6= x ;? h q(x) + y(x)? y d (x) + d = where d(x) := ij =h in view of the scaling (4.8). These equations can be shown to be satised in all gridpoints x := x ij. We verify the second equation in the active 4
15 gridpoints of which there is only one, namely the center of the domain. To be more specific, the following data were obtained at the center point x = x ij = (:5; :5) : q ij =?:3; q i?;j =?:56; y ij = 3:5; y ij? y d (x ij ) = :875; d = ij =h = 46: : Note that the adjoint function is symmetric around x. Hence the adjoint equations holds with? h q(x) = 4 (q ij? q i?;j )=h =?46:4 and y ij? y d (x ij ) + ij =h = 46:39. This example will also be used to illustrate the numerical process in more detail. n it AMPL LOQO Acc F (y; u) y(:4; :5) u(; :5) Table : Detailed information on solution of Example 5. In Table we summarize the data: size of grid, number of iterations of LOQO, times in seconds for the AMPL compilation and the solution by LOQO, the accuracy as the number of correct signicant digits in the objective function value. The primal{ dual approach yields upper and lower bounds for this value and thus permits such a statement. In all the following examples this accuracy measure was at least 8 and is therefore not listed. Finally, we list the objective function value and one representative value of both the state and the control variable. While the state variable not far from the active center point has a small error, the control at the midpoints of the edges still varies in the fourth digit. Thus, n = was chosen for the subsequent calculations. Since the AMPL times were in the range of a few seconds or a few percent of the solution times they will not be given below. n LANCELOT LOQO MINOS SNOPT BPMPD (unbd) (inf) (unbd) Table : Comparison of dierent solvers on Example 5. We compare LOQO and other solvers with AMPL interface available to us in Table. For n = not all the codes considered (LANCELOT-A, MINOS-5.5, SNOPT-5.3, BPMPD-.; links to all codes in []) were successful. Therefore, results are listed for several coarser grids as well as for the nest grid used with LOQO. Given are CPU times on the same platform used above, but scaled to LOQO's time. MINOS reports one problem as infeasible while SNOPT reports two problems 5
16 q.34 u x x.7.59 Y x Figure : Optimal control and adjoint variable for Example 5. as unbounded. The reason for these messages may also be a "bad starting guess". All solvers were given the same AMPL le and AMPL initializes variables with zero. For n = all solvers except LOQO and BPMPD exceeded the available memory (56MB). The convex QP code BPMPD solves with increasing eciency compared to LOQO but this is only possible, because, in fact, Example 5. is a convex QP problem. BPMPD is not applicable to most of the other examples solved below. For another comparison see Example 5.7. Example 5.: All data are the same as in Example 5. except that we choose = : instead of = :. According to the optimal control law (3.) we can expect either a bang{bang or a singular control. We find the following results: Cost functional : F (y; u) = :96695, CPU seconds : 78 The optimal control, adjoint variable and optimal state are shown in Figures and 3. Both the control and state constraints do not become active. Hence, the optimal control is totally singular on?. This is in accordance with the control law (3.) since the normal derivative of the adjoint variable qj? which follows from the numerical results q i in view of (4.4). In Figure 3 as in Figures 5 and 7 below the vertical axis is labeled y=u since for these examples both the state and the control are plotted in these gures. Example 5.3: The data are the same as in Example 5. except that we choose more 6
17 y/u x x.8.56 Y Figure 3: Optimal state for Example 5. restrictive state and control constraints: on :?y(x) = ; y(x) 3:; y d (x) = 3 + 5x (x? )x (x? ); on? : y(x) = u(x); :6 u(x) :3; u d (x) ; = :: We obtain the following results: Cost functional : F (y; u) = :3, CPU seconds : 3 The optimal control, state and adjoint variable are shown in Figures 4 and 5. The optimal control is continuous and has two boundary arcs with u(x) = :3 and one boundary arc with u(x) = :6. The junction points with the boundary are the points x = (:; ); x = (:8; ); x 3 = (:3; ); x 4 = (:77; ); x 5 = (:8; ); x 6 = (:98; ) on the bottom edge of?. By inspecting the switching function in Figure 4, we can verify the control law (3.). The assumption b u underlying (3.) obviously holds. Let us check the control condition on the bottom edge of? at the points x i. Observe q(x i )?q i =h according to (4.9) since we have q i =. Then in view of = h = : the control law (3.) takes the form u(x i ) u i = 8 >< >: q i 4 ; if q i 4 (:6; :3) ; :6 ; if q i 4 :6 ; :3 ; if q i 4 :3 : The active set for the state constraint y(x) 3: is the center point x = x ij = (:5; :5). The dual variable for this active inequality constraint is ij = :647. Again, with these data the validity of the adjoint equation (3.6) can be verified in 9 >= >; 7
18 u.9 q_i*^ x x Figure 4: Optimal control and switching function for Example y/u.37 q x x x x Y Y Figure 5: Optimal state and adjoint variable for Example 5.3 8
19 u.9 q_i*^ x x Figure 6: Optimal control and switching function for Example 5.4 its discretized form. Example 5.4: The data are those from Example 5.3 but now we choose = : expecting to obtain a bang{bang control in contrast to the totally singular control in Example 5.. This expectation is met with the following results: Cost functional : F (y; u) = :4978, CPU seconds : 6 The optimal control shown in Figure 6 is indeed bang{bang. The switching function on the bottom edge of? q(x i )?q i =h according to (4.9). Then Figure 6 clearly illustrates the control law (3.): ( ) :6 ; if qi < ; u(x i ) u i = :3 ; if q i > : The switching points on the bottom edge are x = (:; ); x = (:8; ). Again, the optimal state displayed in Figure 7 is active at the center point x = x ij = (:5; :5). The dual variable for this active inequality constraint is ij = : Neumann boundary conditions Example 5.5: The rst problem has a linear partial dierential equation and nonlinear Neumann boundary conditions with data: on :?y(x) = ; y(x) :7; y d (x) =? (x (x? ) + x (x? )); on? y(x) = u(x)? y(x) ; 3:7 u(x) 4:5; u d (x) ; = :: 9
20 y/u q x x x x Y Y Figure 7: Optimal state and adjoint variable for Example 5.4 The following results were obtained: Cost functional : F (y; u) = :55334, CPU seconds : 494 The optimal control and state are shown in Figures 8 and 9. Note that the sign condition (.5) holds since b y (x; y(x); u(x)) =?y(x). The optimal control is continuous and has two boundary arcs with u(x) = 3:7 and one boundary arc with u(x) = 4:5. The junction points with the boundary are the points x = (:7; ); x = (:34; ); x 3 = (:66; ); x 4 = (:83; ) on the bottom edge of?. The control law (.5) was derived under the assumption b u which obviously holds in this example. By looking at the adjoint variable in Figure 8, the reader may verify that the control law (.5) is satisfied: u(x) = 8 >< >: q(x) ; if q(x) (3:7; 4:5) ; 3:7 ; if q(x) 3:7 ; 4:5 ; if q(x) 4:5 : The active set for the state constraint y(x) :7 are the midpoints of the edges. The dual variable for this active inequality constraint is ij = :4574. It can be checked that the adjoint equations (.7) and (.8) hold observing the scaling (4.) and (4.). As an example, let us test the Neumann boundary condition (.8) at the active point x = x i = (:5; :). Hence, we have to verify the relation q i?q i +q i hy i + i = which corresponds to the q(x)+ q(x)y(x)+d(x) =. We find q i =?:4656; q i =?:48885; y i = :7; i =?:4574 und can 9 >= >;
21 u 4. q x x Figure 8: Optimal control and switching function for Example 5.5 check indeed the Neumann condition with q i? q i + q i h y i = :458 which agrees with? i = :4574. Example 5.6: The data for the second Neumann problem are: on :?y(x) = ; y(x) :835; y d (x) =? (x (x? ) + x (x? )); on? y(x) = u(x)? y(x) ; 6 u(x) 9; u d (x) ; = : According to the optimal control law (.6) we can expect either a bang{bang or a singular control. This is confirmed by the following results: Cost functional : F (y; u) = :578, CPU seconds : 864 The optimal control in Figure is indeed bang{bang. The switching function q on? as shown in Figure obeys the optimal control law (.6): ( ) 6 ; if q(x) < ; u(x) = 9 ; if q(x) > : The switching points are approximately x = (:33; ); x = (:67; ). Again, the optimal state displayed in Figure is active at the midpoints of the edges. The dual variable for this active inequality constraint is ij = :99. Example 5.7: In this example we choose a nonlinear partial dierential equation and linear Neumann boundary conditions: on :?y(x)? y(x) + y(x) 3 = ; y(x) :7; y d (x) =? (x (x? ) + x (x? )); on? y(x) = u(x); :8 u(x) :5; u d (x) ; = ::
22 y. q x x x x.37.3 Y Y Figure 9: Optimal state and adjoint variable for Example u q x x Figure : Optimal control and switching function for Example 5.6
23 y q x x x x Y Y Figure : Optimal state and adjoint variable for Example 5.6 These equations represent a simplified Ginzburg{Landau model for super{conductivity in the absence of internal magnetic fields with y the wave function; cf. Ito, Kunisch [6] and Kunisch, Volkwein [8] where tracking functions and control or state constraints have been considered that are different from those used here. LOQO and AMPL provide the results: Cost functional : F (y; u) = :6463, CPU seconds : 37 The optimal control and state are shown in Figures and 3. The optimal control is continuous and has two boundary arcs with u(x) = :8 and one boundary arc with u(x) = :5. The junction points with the boundary are the points x = (:5; ); x = (:8; ); x 3 = (:7; ); x 4 = (:85; ) on the bottom edge of?. The adjoint variable in Figure shows that the control law (.5) is satisfied: u(x) = 8 >< >: q(x) ; if q(x) (:8; :5) ; :8 ; if q(x) :8 ; :5 ; if q(x) :5 : The active set for the state constraint y(x) :7 comprises the points adjacent to the corners of the domain. The dual variable for this active inequality constraint is ij = : For this problem the QP-code BPMPD is not applicaple. Both MINOS and SNOPT return with error messages as in Example 5., while LANCELOT solves the example on grids of n = 5; 7; 9 and needs respectively 8:4 respectively 7: times as long as LOQO. Example 5.8: The data are those from Example 5.7 but now we choose = : 3 9 >= >;
24 u q x x Figure : Optimal control and switching function for Example 5.7 y q x x x x Y Y Figure 3: Optimal state and adjoint variable for Example 5.7 4
25 expecting to obtain a bang{bang control. The numerical results are: Cost functional : F (y; u) = :6553, CPU seconds : 57 The optimal control shown in Figure 4 is indeed bang{bang. The switching function q on? as shown in Figure 4 yields the optimal control law (.6): ( ) :8 ; if q(x) < ; u(x) = :5 ; if q(x) > : The switching points are approximately x = (:; ); x = (:79; ). Again, the optimal state displayed in Figure 5 is active at the points adjacent to the corners of the domain. The dual variable for this active inequality constraint is ij = :38. 6 Conclusions Numerical techniques have been developed for solving semilinear elliptic control problems subject to control and state constraints. In the first part of two papers on this subject, boundary control problems with either Dirichlet or Neumann boundary conditions have been treated. We have proposed discretization schemes by which the elliptic control problem is transcribed into a nonlinear programming problem. Since both control and state variables are treated as optimization variables, the resulting NLP{problem is a high{dimensional problem with a sparse structure. Testing dierent optimization codes, we found that the interior point method LOQO of Vanderbei, Shanno [3] is particularly efficient and reliable. The algorithm provided solution for all types of controls: regular controls, bang{bang controls and singular controls. In each case, we have checked the necessary conditions of optimality for the adjoint variables. After nishing this paper, we learned about the paper of Bergounioux et al. [] where interior point methods are compared with methods developed in [3]. We are currently extending our approach to include a test of second order sufficient conditions (SSC). The test of (SSC) can be performed by showing that the projected Hessian of the Lagrangian is positive definite. The importance of (SSC) lies in the fact that it provides a basis for sensitivity analysis and real{time control of perturbed systems; cf., e.g., [7, 8]. Presently, we are carrying over these techniques to elliptic control problems. Acknowledgement: We would like to thank Professor Fredi Troltzsch for discussions on the necessary conditions for Dirichlet boundary conditions. References [] A. Barclay, P.E. Gill, and J. B. Rosen, \SQP methods and their applications to numerical optimal control", to appear. 5
26 u q x x Figure 4: Optimal control and switching variable for Example y.3.3 q x x x x Y Y Figure 5: Optimal state and adjoint variable for Example 5.8 6
27 [] M. Bergounioux, M. Haddou, M. Hintermuller and K. Kunisch, \A comparison of interior point methods and a Moreau{Yosida based active set strategy for constrained optimal control problems", Preprint, 998. [3] M. Bergounioux and K. Kunisch, \Augmented Lagrangian techniques for elliptic state constrained optimal control problems", SIAM J. Control Optim., vol. 35, pp , 997. [4] J.T. Betts, \Issues in the direct transcription of optimal control problems to sparse nonlinear programs", in Control Applications of Optimization, R. Bulirsch and D. Kraft, eds, vol. 5, Int. Series Numer. Mathematics, Basel, Birkhauser, 994, pp [5] J.T. Betts and W.P. Humann, \The application of sparse nonlinear programming to trajectory optimization", J. of Guidance, Control and Dynamics, vol. 4, pp , 99. [6] F. Bonnans, \Second order analysis for control constrained optimal control problems of semilinear elliptic systems", to appear in Applied Mathematics and Optimization. [7] C. Buskens, \Optimierungsmethoden und Sensitivitatsanalyse fur optimale Steuerprozesse mit Steuer- und ustandsbeschrankungen", Dissertation, Universitat Munster, Institut fur Numerische Mathematik, Munster, Germany, 998. [8] C. Buskens and H. Maurer, \Real{time control of an industrial robot using nonlinear programming methods", to appear in Automatica. [9] E. Casas, \Boundary control with pointwise state constraints", SIAM J. Control Optim., vol. 3, pp , 993 [] E. Casas, F. Troltzsch and A. Unger, \Second order sufficient optimality conditions for a nonlinear elliptic control problem", J. for Analysis and its Applications, vol 5, pp , 996. [] E. Casas, F. Troltzsch and A. Unger, \Second order sufficient optimality conditions for some state constrained control problems of semilinear elliptic equations", Fakultat fur Mathematik, Technische Universitat Chemnitz, Preprint 97-9, to appear in SIAM J. Control Optim. [] R. Fourer, D. M. Gay, and B. W. Kernighan, \AMPL: A modeling Language for Mathematical Programming", Duxbury Press, Brooks-Cole Publishing Company, 993 [3] I.I. Grachev and Yu.G. Evtushenko, \A library of programs for solving optimal control problems", U.S.S.R. Computational Maths. Math. Physics, vol. 9, pp. 99-9, 979. [4] R. Hettich, A. Kaplan and R. Tischatschke, \Regularized penalty methods for illposed optimal control problems with elliptic equations. Part I: Distributed control with bounded control set and state constraints", Control and Cybernetics, vol. 6, pp. 5-7,
28 [5] R. Hettich, A. Kaplan and R. Tischatschke, \Regularized penalty methods for illposed optimal control problems with elliptic equations. Part II: Distributed and boundary control with unbounded control sets and state constraints", Control and Cybernetics, vol. 6, pp. 9-43, 997. [6] K. Ito and K. Kunisch, \Augmented Lagrangian{SQP methods for nonlinear optimal control problems of tracking type", SIAM J. Optim., vol 6, pp. 96-5, 996. [7] D. Kraft, \On converting optimal control problems into nonlinear programming problems", in Computational Mathematical Programming, K. Schittkowski, ed., NATO ASI Series F: Computer and Systems Science, vol. 5, Springer Verlag, Berlin und Heidelberg, 985, pp [8] K. Kunisch and S. Volkwein, \Augmented Lagrangian-SQP techniques and their approximations", Contemporary Mathematics, vol. 9, pp , 997. [9] J.L. Lions, \Optimal control of systems governed by partial dierential equations", Grundlehren der mathematischen Wissenschaften, Vol. 7, Springer{Verlag, Berlin, New York, 97. [] J.L. Lions and E. Magenes, \Non{Homogeneous Boundary Value Problems and Applications, Volume I", Grundlehren der mathematischen Wissenschaften, Vol. 8, Springer{Verlag, Berlin, New York, 97. [] H.D. Mittelmann and P. Spellucci, \Decision Tree for Optimization Software", World Wide Web, (998). [] K.L. Teo, C.J. Goh and K.H. Wong, \A Unied Computational Approach to Optimal Control Problems", Longman Scientic and Technical, New York, 98. [3] R. S. Vanderbei and D. F. Shanno, \An interior point algorithm for nonconvex nonlinear programming", to appear. [4] J. owe and S. Kurcyusz, \Regularity and stability for the mathematical programming problem in Banach spaces", Appl. Math. Optimization, vol. 5, pp. 49-6,
Optimization Techniques for Solving Elliptic Control Problems with Control and State Constraints. Part 2: Distributed Control
Computational Optimization and Applications, 18, 141 160, 2001 c 2001 Kluwer Academic Publishers. Manufactured in The Netherlands. Optimization Techniques for Solving Elliptic Control Problems with Control
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 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 informationTheory and Applications of Constrained Optimal Control Proble
Theory and Applications of Constrained Optimal Control Problems with Delays PART 1 : Mixed Control State Constraints Helmut Maurer 1, Laurenz Göllmann 2 1 Institut für Numerische und Angewandte Mathematik,
More informationZentrum für Technomathematik
Zentrum für Technomathematik Fachbereich 3 Mathematik und Informatik Computational Parametric Sensitivity Analysis of Perturbed PDE Optimal Control Problems with State and Control Constraints Christof
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 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 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 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 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 informationWhat s New in Active-Set Methods for Nonlinear Optimization?
What s New in Active-Set Methods for Nonlinear Optimization? Philip E. Gill Advances in Numerical Computation, Manchester University, July 5, 2011 A Workshop in Honor of Sven Hammarling UCSD Center for
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 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 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 informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization 1 Origins Instructor: Michael Saunders Spring 2015 Notes 9: Augmented Lagrangian Methods
More informationWe consider the application of an SQP method to an optimal control. programming problems have to be treated. The solution of the constrained
Optimal Control: Theory, Algorithms, and Applications, pp. 1-2 W. W. Hager and P. M. Pardalos, Editors c1997 Kluwer Academic Publishers B.V. On a SQP-Multigrid Technique for Nonlinear Parabolic Boundary
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 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 informationINTERIOR-POINT METHODS FOR NONCONVEX NONLINEAR PROGRAMMING: COMPLEMENTARITY CONSTRAINTS
INTERIOR-POINT METHODS FOR NONCONVEX NONLINEAR PROGRAMMING: COMPLEMENTARITY CONSTRAINTS HANDE Y. BENSON, DAVID F. SHANNO, AND ROBERT J. VANDERBEI Operations Research and Financial Engineering Princeton
More informationBang-Bang and Singular Controls in Optimal Control Problems with Partial Differential Equations
Bang-Bang and Singular Controls in Optimal Control Problems with Partial Differential Equations Hans Josef Pesch, Simon Bechmann, and Jan-Eric Wurst Abstract This paper focuses on bang-bang and singular
More informationTutorial on Control and State Constrained Optimal Control Problems
Tutorial on Control and State Constrained Optimal Control Problems To cite this version:. blems. SADCO Summer School 211 - Optimal Control, Sep 211, London, United Kingdom. HAL Id: inria-629518
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationOptimization Problems with Constraints - introduction to theory, numerical Methods and applications
Optimization Problems with Constraints - introduction to theory, numerical Methods and applications Dr. Abebe Geletu Ilmenau University of Technology Department of Simulation and Optimal Processes (SOP)
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationOptimality, Duality, Complementarity for Constrained Optimization
Optimality, Duality, Complementarity for Constrained Optimization Stephen Wright University of Wisconsin-Madison May 2014 Wright (UW-Madison) Optimality, Duality, Complementarity May 2014 1 / 41 Linear
More informationNonlinear Optimization Solvers
ICE05 Argonne, July 19, 2005 Nonlinear Optimization Solvers Sven Leyffer leyffer@mcs.anl.gov Mathematics & Computer Science Division, Argonne National Laboratory Nonlinear Optimization Methods Optimization
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationSolving Dual Problems
Lecture 20 Solving Dual Problems We consider a constrained problem where, in addition to the constraint set X, there are also inequality and linear equality constraints. Specifically the minimization problem
More informationAn Inexact Sequential Quadratic Optimization Method for Nonlinear Optimization
An Inexact Sequential Quadratic Optimization Method for Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with Travis Johnson, Northwestern University Daniel P. Robinson, Johns
More informationA NEW COMPUTATIONAL METHOD FOR OPTIMAL CONTROL OF A CLASS OF CONSTRAINED SYSTEMS GOVERNED BY PARTIAL DIFFERENTIAL EQUATIONS
A NEW COMPUTATIONAL METHOD FOR OPTIMAL CONTROL OF A CLASS OF CONSTRAINED SYSTEMS GOVERNED BY PARTIAL DIFFERENTIAL EQUATIONS Nicolas Petit Mark Milam Richard Murray Centre Automatique et Systèmes École
More information2 A. BARCLAY, P. E. GILL AND J. B. ROSEN where c is a vector of nonlinear functions, A is a constant matrix that denes the linear constraints, and b l
SQP METHODS AND THEIR APPLICATION TO NUMERICAL OPTIMAL CONTROL ALEX BARCLAY y, PHILIP E. GILL y, AND J. BEN ROSEN z Abstract. In recent years, general-purpose sequential quadratic programming (SQP) methods
More informationGeneralization to inequality constrained problem. Maximize
Lecture 11. 26 September 2006 Review of Lecture #10: Second order optimality conditions necessary condition, sufficient condition. If the necessary condition is violated the point cannot be a local minimum
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 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 informationAlgorithms for Constrained Optimization
1 / 42 Algorithms for Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University April 19, 2015 2 / 42 Outline 1. Convergence 2. Sequential quadratic
More informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
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 informationSolving generalized semi-infinite programs by reduction to simpler problems.
Solving generalized semi-infinite programs by reduction to simpler problems. G. Still, University of Twente January 20, 2004 Abstract. The paper intends to give a unifying treatment of different approaches
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 informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More information4. Algebra and Duality
4-1 Algebra and Duality P. Parrilo and S. Lall, CDC 2003 2003.12.07.01 4. Algebra and Duality Example: non-convex polynomial optimization Weak duality and duality gap The dual is not intrinsic The cone
More informationPacific Journal of Optimization (Vol. 2, No. 3, September 2006) ABSTRACT
Pacific Journal of Optimization Vol., No. 3, September 006) PRIMAL ERROR BOUNDS BASED ON THE AUGMENTED LAGRANGIAN AND LAGRANGIAN RELAXATION ALGORITHMS A. F. Izmailov and M. V. Solodov ABSTRACT For a given
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 informationInterior-Point Methods as Inexact Newton Methods. Silvia Bonettini Università di Modena e Reggio Emilia Italy
InteriorPoint Methods as Inexact Newton Methods Silvia Bonettini Università di Modena e Reggio Emilia Italy Valeria Ruggiero Università di Ferrara Emanuele Galligani Università di Modena e Reggio Emilia
More informationConstrained Optimization
1 / 22 Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University March 30, 2015 2 / 22 1. Equality constraints only 1.1 Reduced gradient 1.2 Lagrange
More informationThe use of second-order information in structural topology optimization. Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher
The use of second-order information in structural topology optimization Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher What is Topology Optimization? Optimize the design of a structure
More informationSupport Vector Regression with Automatic Accuracy Control B. Scholkopf y, P. Bartlett, A. Smola y,r.williamson FEIT/RSISE, Australian National University, Canberra, Australia y GMD FIRST, Rudower Chaussee
More informationPart 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL)
Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization Nick Gould (RAL) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
More information460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses
Statistica Sinica 5(1995), 459-473 OPTIMAL DESIGNS FOR POLYNOMIAL REGRESSION WHEN THE DEGREE IS NOT KNOWN Holger Dette and William J Studden Technische Universitat Dresden and Purdue University Abstract:
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 informationLimit Analysis with the. Department of Mathematics and Computer Science. Odense University. Campusvej 55, DK{5230 Odense M, Denmark.
Limit Analysis with the Dual Ane Scaling Algorithm Knud D. Andersen Edmund Christiansen Department of Mathematics and Computer Science Odense University Campusvej 55, DK{5230 Odense M, Denmark e-mail:
More informationNumerical Optimal Control Part 2: Discretization techniques, structure exploitation, calculation of gradients
Numerical Optimal Control Part 2: Discretization techniques, structure exploitation, calculation of gradients SADCO Summer School and Workshop on Optimal and Model Predictive Control OMPC 2013, Bayreuth
More informationLagrange duality. The Lagrangian. We consider an optimization program of the form
Lagrange duality Another way to arrive at the KKT conditions, and one which gives us some insight on solving constrained optimization problems, is through the Lagrange dual. The dual is a maximization
More informationAlgorithms for constrained local optimization
Algorithms for constrained local optimization Fabio Schoen 2008 http://gol.dsi.unifi.it/users/schoen Algorithms for constrained local optimization p. Feasible direction methods Algorithms for constrained
More informationLIPSCHITZ STABILITY OF SOLUTIONS OF LINEAR-QUADRATIC PARABOLIC CONTROL PROBLEMS WITH RESPECT TO PERTURBATIONS. FREDI TR OLTZSCH y
LIPSCHIT STABILITY OF SOLUTIONS OF LINEAR-QUADRATIC PARABOLIC CONTROL PROBLEMS WITH RESPECT TO PERTURBATIONS FREDI TR OLTSCH y Abstract. We consider a class of control problems governed by a linear parabolic
More informationNew Identities for Weak KAM Theory
New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities
More informationSQP METHODS AND THEIR APPLICATION TO NUMERICAL OPTIMAL CONTROL
c Birkhäuser Verlag Basel 207 SQP METHODS AND THEIR APPLICATION TO NUMERICAL OPTIMAL CONTROL ALEX BARCLAY PHILIP E. GILL J. BEN ROSEN Abstract. In recent years, general-purpose sequential quadratic programming
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
More informationNUMERICAL OPTIMIZATION. J. Ch. Gilbert
NUMERICAL OPTIMIZATION J. Ch. Gilbert Numerical optimization (past) The discipline deals with the classical smooth (nonconvex) problem min {f(x) : c E (x) = 0, c I (x) 0}. Applications: variable added
More informationNumerisches Rechnen. (für Informatiker) M. Grepl P. Esser & G. Welper & L. Zhang. Institut für Geometrie und Praktische Mathematik RWTH Aachen
Numerisches Rechnen (für Informatiker) M. Grepl P. Esser & G. Welper & L. Zhang Institut für Geometrie und Praktische Mathematik RWTH Aachen Wintersemester 2011/12 IGPM, RWTH Aachen Numerisches Rechnen
More informationLecture 3 Interior Point Methods and Nonlinear Optimization
Lecture 3 Interior Point Methods and Nonlinear Optimization Robert J. Vanderbei April 16, 2012 Machine Learning Summer School La Palma http://www.princeton.edu/ rvdb Example: Basis Pursuit Denoising L
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More informationSecond order sucient optimality condition for a. nonlinear elliptic control problem. Eduardo Casas and Fredi Troltzsch and Andreas Unger
Second order sucient optimality condition for a nonlinear elliptic control problem by Eduardo Casas and Fredi Troltzsch and Andreas Unger Technical University of Chemnitz wickau Faculty of Mathematics
More informationLagrange Relaxation and Duality
Lagrange Relaxation and Duality As we have already known, constrained optimization problems are harder to solve than unconstrained problems. By relaxation we can solve a more difficult problem by a simpler
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationReal-time Constrained Nonlinear Optimization for Maximum Power Take-off of a Wave Energy Converter
Real-time Constrained Nonlinear Optimization for Maximum Power Take-off of a Wave Energy Converter Thomas Bewley 23 May 2014 Southern California Optimization Day Summary 1 Introduction 2 Nonlinear Model
More information10 Numerical methods for constrained problems
10 Numerical methods for constrained problems min s.t. f(x) h(x) = 0 (l), g(x) 0 (m), x X The algorithms can be roughly divided the following way: ˆ primal methods: find descent direction keeping inside
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 informationPrimal-dual Subgradient Method for Convex Problems with Functional Constraints
Primal-dual Subgradient Method for Convex Problems with Functional Constraints Yurii Nesterov, CORE/INMA (UCL) Workshop on embedded optimization EMBOPT2014 September 9, 2014 (Lucca) Yu. Nesterov Primal-dual
More informationPreprint Stephan Dempe and Patrick Mehlitz Lipschitz continuity of the optimal value function in parametric optimization ISSN
Fakultät für Mathematik und Informatik Preprint 2013-04 Stephan Dempe and Patrick Mehlitz Lipschitz continuity of the optimal value function in parametric optimization ISSN 1433-9307 Stephan Dempe and
More informationA VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION 1. INTRODUCTION
A VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION GERARD AWANOU AND LEOPOLD MATAMBA MESSI ABSTRACT. We give a proof of existence of a solution to the discrete problem
More informationImplicitly and Explicitly Constrained Optimization Problems for Training of Recurrent Neural Networks
Implicitly and Explicitly Constrained Optimization Problems for Training of Recurrent Neural Networks Carl-Johan Thore Linköping University - Division of Mechanics 581 83 Linköping - Sweden Abstract. Training
More informationTMA947/MAN280 APPLIED OPTIMIZATION
Chalmers/GU Mathematics EXAM TMA947/MAN280 APPLIED OPTIMIZATION Date: 06 08 31 Time: House V, morning Aids: Text memory-less calculator Number of questions: 7; passed on one question requires 2 points
More informationSparsity Matters. Robert J. Vanderbei September 20. IDA: Center for Communications Research Princeton NJ.
Sparsity Matters Robert J. Vanderbei 2017 September 20 http://www.princeton.edu/ rvdb IDA: Center for Communications Research Princeton NJ The simplex method is 200 times faster... The simplex method is
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationCHAPTER 2: QUADRATIC PROGRAMMING
CHAPTER 2: QUADRATIC PROGRAMMING Overview Quadratic programming (QP) problems are characterized by objective functions that are quadratic in the design variables, and linear constraints. In this sense,
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 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 informationIE 5531 Midterm #2 Solutions
IE 5531 Midterm #2 s Prof. John Gunnar Carlsson November 9, 2011 Before you begin: This exam has 9 pages and a total of 5 problems. Make sure that all pages are present. To obtain credit for a problem,
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 informationA CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION
Trends in Mathematics Information Center for Mathematical Sciences Volume 4, Number 2, December 21, Pages 39 44 A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationFirst-Order Optimally Conditions in Generalized Semi-Infinite Programming 1
JOURNAL OF OPTIMIZATION THEORY AND APPLICATIONS: Vol. 101, No. 3, pp. 677-691, JUNE 1999 First-Order Optimally Conditions in Generalized Semi-Infinite Programming 1 J. J. RUCKMANN 2 AND A. SHAPIRO 3 Communicated
More informationPoisson Solvers. William McLean. April 21, Return to Math3301/Math5315 Common Material.
Poisson Solvers William McLean April 21, 2004 Return to Math3301/Math5315 Common Material 1 Introduction Many problems in applied mathematics lead to a partial differential equation of the form a 2 u +
More informationThe Lifted Newton Method and Its Use in Optimization
The Lifted Newton Method and Its Use in Optimization Moritz Diehl Optimization in Engineering Center (OPTEC), K.U. Leuven, Belgium joint work with Jan Albersmeyer (U. Heidelberg) ENSIACET, Toulouse, February
More information10-725/ Optimization Midterm Exam
10-725/36-725 Optimization Midterm Exam November 6, 2012 NAME: ANDREW ID: Instructions: This exam is 1hr 20mins long Except for a single two-sided sheet of notes, no other material or discussion is permitted
More informationKunisch [4] proved that the optimality system is slantly dierentiable [7] and the primal-dual active set strategy is a specic semismooth Newton method
Finite Dierence Smoothing Solution of Nonsmooth Constrained Optimal Control Problems Xiaojun Chen 22 April 24 Abstract The nite dierence method and smoothing approximation for a nonsmooth constrained optimal
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationApplications of Linear Programming
Applications of Linear Programming lecturer: András London University of Szeged Institute of Informatics Department of Computational Optimization Lecture 9 Non-linear programming In case of LP, the goal
More informationCONVEXIFICATION SCHEMES FOR SQP METHODS
CONVEXIFICAION SCHEMES FOR SQP MEHODS Philip E. Gill Elizabeth Wong UCSD Center for Computational Mathematics echnical Report CCoM-14-06 July 18, 2014 Abstract Sequential quadratic programming (SQP) methods
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 informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationINTERIOR-POINT ALGORITHMS, PENALTY METHODS AND EQUILIBRIUM PROBLEMS
INTERIOR-POINT ALGORITHMS, PENALTY METHODS AND EQUILIBRIUM PROBLEMS HANDE Y. BENSON, ARUN SEN, DAVID F. SHANNO, AND ROBERT J. VANDERBEI Abstract. In this paper we consider the question of solving equilibrium
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 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 information