Error analysis for global minima of semilinear optimal control problems
Dedicated to Eduardo Casas on the occasion of his 60th birthday.
Abstract: In  we consider an optimal control problem subject to a semilinear elliptic PDE together with its variational discretization, where we provide a condition which allows to decide whether a solution of the necessary first order conditions is a global minimum. This condition can be explicitly evaluated at the discrete level. Furthermore, we prove that if the above condition holds uniformly with respect to the discretization parameter the sequence of discrete solutions converges to a global solution of the corresponding limit problem. With the present work we complement our investigations of  in that we prove an error estimate for those discrete global solutions. Numerical experiments confirm our analytical findings.
Mathematics Subject Classification (2000): 49J20, 35K20, 49M05, 49M25, 49M29, 65M12, 65M60
Keywords: Optimal control, semilinear PDE, uniqueness of global solutions, error estimates
In this work we are concerned with the error analysis of a variational discretization of the control problem
and the pointwise constraints
where the precise assumptions on the data of the problem will be given in Section 2.1. In  the authors considered the same class of problems and established a sufficient condition for the global minima of assuming particular types of growth conditions for the nonlinearity . The same result was established for the variational discrete counterpart of , and it was shown that a sequence of the computed discrete global minima converges to a global minimum of the continuous control problem but without discussing the corresponding rate of convergence. Hence, our aim in this study is to investigate this convergence rate.
The organization of the paper is as follows: in § 2.1 we formulate the control problem and give the exact assumptions on the data. In § 2.2 and in § 2.3 we review the results concerning the state equation and the control problem , respectively. The variational discretization of is considered in § 2.4 while § 3 is devoted to the error analysis. Finally, in § 4 we verify our theoretical findings by a numerical example.
Before starting, we give a short list of literature considering the problem . For a broad overview, we refer the reader to the references of the respective citations. In  the problem is studied when the controls are of boundary type, and the necessary first order conditions are established. Compare  where the function is linear, and  where the pointwise constraints are imposed on the gradient of the state.
The regularity of the optimal controls of and their associated multipliers are investigated in  and , where also the sufficient second order conditions are discussed. Compare [9, 6, 7] for second order conditions when the set contains finitely/infinitely many points, and  for the role of those conditions in PDE constrained control problems.
Finite element discretization of problem under more general setting is studied in , and in  where a wider class of perturbations are considered. The convergence of the discrete solutions to the continuous solutions is verified there but without rates. However, when the set contains finitely many points, convergence rates are established in  for finite dimensional controls, and in  for control functions. Only in  error analysis is studied for general pointwise state constraints in . There, Pfefferer at al. prove an error estimate for discrete solutions in the vicinity of a local solution which satisfies a quadratic growth condition. Error analysis for linear-quadratic control problems can be found in e.g. ,  and . A detailed discussion of discretization concepts and error analysis in PDE-constrained control problems can be found in [20, 21] and[17, Chapter 3].
2 Problem Setting and discretization
is a bounded, convex and polygonal domain.
is a (possibly empty) compact subset of .
and with .
are given functions that satisfy , .
and are given.
is of class and monotonically increasing.
There exist and such that
where and denote the first and second derivative of , respectively.
2.2 The State Equation
The existence and uniqueness of the solution in follows from the monotone operator theorem. Using the method of Stampacchia one can show, in addition, that . Utilizing the boundedness of and the properties of the nonlinearity , one can show and the estimate (5) using the regularity results from [16, Chapter 4]. For a detailed proof compare for instance . ∎
In the light of Theorem 2.1, we introduce the control–to–state operator
such that is the solution to (4) for a given .
Let be the mapping introduced in (6). Then there exists depending only on such that
Given let and . Using Poincaré’s inequality, the monotonicity of and (4) we have
which implies the result. ∎
Let be the mapping introduced in (6). Then for any there exists such that
for all , with , .
Defining again we infer from Theorem 2.1 and the continuous embedding that , for some depending on . Clearly, belongs to and satisfies
Using a standard a–priori estimate, the Lipschitz continuity of on bounded sets and Lemma 2.2 we infer that
where is a constant depending on . This completes the proof. ∎
2.3 The Optimal Control Problem
Using the control-to-state operator defined in (6), the reduced form of our optimal control problem reads
It is well–known that admits at least one solution provided that a feasible point exists (compare ). Moreover, if a solution of satisfies some constraint qualification, then one can guarantee the existence of a multiplier associated with the pointwise state constraints and the necessary first order conditions can be established. A typical constraint qualification for a local solution of problem is the linearized Slater condition which reads: there exist and such that
The next result is a consequence of [4, Theorem 5.2].
Let be a local solution of problem satisfying (7). Then there exist for and a regular Borel measure such that with there holds
Note that in view of (10) is the –projection of onto so that
Since for it follows from [22, Corollary A.6] that for as well. Furthermore, it is well known that the multiplier associated with the pointwise state constraints is concentrated at the points in where the state constraints are active. We state this more precisely in the next proposition whose proof can be found in . Compare also the proof in  when the bounds , are constant functions.
Let and satisfy (11). Then there holds
where with is the Jordan decomposition of .
We note that the problem is in general nonconvex since the state equation is not linear. In other words, the problem can have several solutions. A decision of which of these solution is a global minimum proves difficult in general. However, it is shown in  that if the nonlinearity of the state equation enjoys certain growth conditions, namely (3), then one can establish a condition that helps to decide if a given point satisfying the first order conditions is a global minimum. We state this condition of global optimality in the next result, but before that we first need to introduce the following constant:
Here, , while and appear in (3). Furthermore, is an upper bound on the optimal constant in the Gagliardo-Nirenberg inequality
For sharp upper bounds for the constant , see for instance [1, Theorem 7.3].
2.4 Variational Discretization
Let be an admissible triangulation of the polygonal domain with
Here is the maximum mesh size, while stands for the diameter of the triangle . We introduce the following spaces of linear finite elements:
The Lagrange interpolation operator is defined by
where denote the nodes in the triangulation and are the basis functions of the space which satisfy .
The finite element discretization of (4) reads: for a given , find such that
such that is the solution of (14).
The variational discretization (see ) of Problem reads:
where we define
with the set of nodes
We remark that provided that is small enough. This follows from the fact that and are continuous functions with .
In an analogous way to that of problem , one can show that admits at least one solution, denoted by , provided that a feasible point exists. In practice one calculates candidates for solutions of by solving the system of necessary first order conditions which reads: find and such that
where is the Dirac measure at . We can easily deduce from (19) the following result about the support of the measure .
3 Error Analysis
Let be a sequence of admissible triangulations of . We assume that the sequence is quasi-uniform in the sense that each is contained in a ball of radius and contains a ball of radius for some independent of . In addition we make the following assumption concerning the set :
For every there exists a set of triangles such that
for some that is independent of . We immediately infer from Theorem 2.8 that is the unique global minimum of and we are interested in the convergence properties of these solutions as . It is shown in  (see Theorem 4.2 and its proof) that there exist , and such that
Theorem 2.6 implies that is the unique global optimum of . The aim in the remaining part of this paper is to prove error estimates for and the corresponding optimal states . Our main results read:
Suppose that (23) holds and let be the unique global minima of and respectively. Then we have for any that
In  Pfefferer at al. for problems in two and three dimensions present a similar error estimate for discrete (local) solutions in the vicinity of a local solution which satisfies a quadratic growth condition. Assuming (23) we here use different techniques to prove an error estimate for the unique global discrete solutions which converge to the unique global solution of our optimization problem.
Before we start presenting the proof of this result we collect some results concerning the uniform boundedness of the discrete optimal control , its state and the associated multipliers and .
Next, let us introduce the auxiliary functions , , as the solutions of
Let and be as above and an open set such that and . Then we have
from which we obtain