This paper investigates the sensitivity of optimal control problems governed by nonlinear evolution inclusions, establishing well-posedness and continuity properties, with an application to nonlinear parabolic systems.
Contribution
It provides a general framework for analyzing the solution set and sensitivity of nonlinear evolution inclusion-based control problems, including continuity of solutions and value functions.
Findings
01
Solution set nonemptiness and continuous selection established
02
Hadamard well-posedness of the control problem proved
03
Continuity of the optimal solution multifunction demonstrated
Abstract
We consider a nonlinear optimal control problem governed by a nonlinear evolution inclusion and depending on a parameter λ. First we examine the dynamics of the problem and establish the nonemptiness of the solution set and produce continuous selections of the solution multifunction ξ↦S(ξ) (ξ being the initial condition). These results are proved in a very general framework and are of independent interest as results about evolution inclusions. Then we use them to study the sensitivity properties of the optimal control problem. We show that we have Hadamard well-posedness (continuity of the value function) and we establish the continuity properties of the optimal multifunction. Finally we present an application on a nonlinear parabolic distributed parameter system.
Equations659
\displaystyle\displaystyle\left\{\begin{array}[]{l}J(x,u,\lambda)=\int^{b}_{0}L(t,x(t),\lambda)dt+\int^{b}_{0}H(t,u(t),\lambda)dt+\hat{\psi}(\xi,x(b),\lambda)\rightarrow\inf=\\
\hskip 284.52756ptm(\xi,\lambda),\\
-x^{\prime}(t)\in A_{\lambda}(t,x(t))+F(t,x(t),\lambda)+G(t,u(t),\lambda)\ \mbox{for almost all}\ t\in T,\\
x(0)=\xi,u(t)\in U(t,\lambda)\ \mbox{for almost all}\ t\in T,\lambda\in E.\end{array}\right\}
\displaystyle\displaystyle\left\{\begin{array}[]{l}J(x,u,\lambda)=\int^{b}_{0}L(t,x(t),\lambda)dt+\int^{b}_{0}H(t,u(t),\lambda)dt+\hat{\psi}(\xi,x(b),\lambda)\rightarrow\inf=\\
\hskip 284.52756ptm(\xi,\lambda),\\
-x^{\prime}(t)\in A_{\lambda}(t,x(t))+F(t,x(t),\lambda)+G(t,u(t),\lambda)\ \mbox{for almost all}\ t\in T,\\
x(0)=\xi,u(t)\in U(t,\lambda)\ \mbox{for almost all}\ t\in T,\lambda\in E.\end{array}\right\}
\left\{\begin{array}[]{l}-x^{\prime}(t)\in A(t,x(t))+E(t,x(t))\ \mbox{for almost all}\ t\in T,\\
x(0)=\xi.\end{array}\right\}
\left\{\begin{array}[]{l}-x^{\prime}(t)\in A(t,x(t))+E(t,x(t))\ \mbox{for almost all}\ t\in T,\\
x(0)=\xi.\end{array}\right\}
∣∣h∣∣∗⩽a1(t)+c1∣∣x∣∣p−1
∣∣h∣∣∗⩽a1(t)+c1∣∣x∣∣p−1
⟨h∗,x⟩⩾c2∣∣x∣∣p−a2(t),
⟨h∗,x⟩⩾c2∣∣x∣∣p−a2(t),
A(x)=−div∂φ(Dx)−divξ(Dx)
A(x)=−div∂φ(Dx)−divξ(Dx)
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Sensitivity analysis for optimal control problems governed by nonlinear evolution inclusions
Nikolaos S. Papageorgiou
Department of Mathematics, National Technical University,
Zografou Campus, 15780 Athens, Greece
Institute of Mathematics “Simion Stoilow” of the Romanian Academy, P.O. Box 1-764,
014700 Bucharest, Romania &
Department of Mathematics, University of Craiova, Street A. I. Cuza 13,
200585 Craiova, Romania
We consider a nonlinear optimal control problem governed by a nonlinear evolution inclusion and depending on a parameter λ. First we examine the dynamics of the problem and establish the nonemptiness of the solution set and produce continuous selections of the solution multifunction ξ↦S(ξ) (ξ being the initial condition). These results are proved in a very general framework and are of independent interest as results about evolution inclusions. Then we use them to study the sensitivity properties of the optimal control problem. We show that we have Hadamard well-posedness (continuity of the value function) and we establish the continuity properties of the optimal multifunction. Finally we present an application on a nonlinear parabolic distributed parameter system.
Key words and phrases:
Evolution triple, evolution inclusion, PG and G convergence, compact embedding, Hadamard well-posedness.
aa 2010 AMS Subject Classification: Primary 49J24, 49K20, Secondary 34G20
1. Introduction
One of the important problems in optimal control theory, is the study of the variations of the set of optimal state-control pairs and of the value of the problem, when we perturb the dynamics, the cost functional and the initial condition of the problem. Such a sensitivity analysis (also known in the literature as “variational analysis”) is important because it gives information about the tolerances which are permitted in the specification of the mathematical models, it suggests ways to solve parametric problems and also can be useful in the computational analysis of the problem. For infinite dimensional systems (distributed parameter systems), such investigations were conducted by Buttazzo and Dal Maso [8], Denkowski and Migorski [13], Ito and Kunisch [24], Papageorgiou [31] (linear systems), Papageorgiou [30], Sokolowski [38] (semilinear systems) and Hu and Papageorgiou [23], Papageorgiou [32, 33] (nonlinear systems). We also mention
the books of Buttazzo [7], Dontchev and Zolezzi [17], Ito and Kunisch [25], Sokolowski and Zolezio [39] (the latter for shape optimization problems). In this paper we conduct such an analysis for a very general class of systems driven by nonmonotone evolution inclusions.
So, let T=[0,b] be the time interval and (X,H,X∗) an evolution triple of spaces (see Section 2). We assume that X↪H compactly. The space of controls is modelled by a separable reflexive Banach space Y and E is a compact metric space and corresponds to the parameter space. As we have already mentioned, we consider systems monitored by evolution inclusions. These inclusions represent a way to model systems with deterministic uncertainties, see the books of Aubin and Frankowska [2], Fattorini [18], and Roubicek [37].
The problem under consideration is the following:
[TABLE]
In this problem
[TABLE]
and the precise conditions on them will be given in Section 4. For every initial state ξ∈H and every parameter λ∈E, we denote the set of admissible state-control pairs (that is, pairs (x,u) which satisfy the dynamics and the constraints of problem (5)) by Q(ξ,λ). We investigate the dependence of Q(ξ,λ) on the two variables (ξ,λ)∈H×E. Also, Σ(ξ,λ) denotes the set of optimal state-control pairs (that is, (x∗,u∗)∈Σ(ξ,λ) such that J(x∗,u∗,ξ,λ)=m(ξ,λ)). So, Σ(ξ,λ)⊆Q(ξ,λ). We establish the nonemptiness of the set Σ(ξ,λ) and examine the continuity properties of the value function (ξ,λ)↦m(ξ,λ) and of the multifunction (ξ,λ)↦Σ(ξ,λ).
The nonemptiness and other continuity and structural properties of the set Q(ξ,λ) are consequences of general results about evolution inclusions which we prove in Section 3 and which are of independent interest. The class of evolution inclusions considered in Section 3 is more general than the classes studied by Chen, Wang and Zu [11], Denkowski, Migorski and Papageorgiou [16], Liu [28], Papageorgiou and Kyritsi [34].
In the next section, for the convenience of the reader, we review the main mathematical tools which we will need in this paper.
2. Mathematical Background
Suppose that V and Z are Banach spaces and assume that V is embedded continuously and densely into Z (denoted by V↪Z). Then it is easy to check that
•
Z∗ is embedded continuously into V∗;
•
if V is reflexive, then Z∗↪V∗.
Having this observation in mind, we can introduce the notion of evolution triple of spaces, which is central in the class of evolution equations considered here.
Definition 1**.**
A triple (X,H,X∗) of spaces is said to be an “evolution triple” (or “Gelfand triple” or “spaces in normal position”), if the following hold:
(a)
X* is a separable reflexive Banach space and X∗ is its topological dual;*
(b)
H* is a separable Hilbert space identified with its dual H∗=H (pivot space);*
(c)
X↪H.
According to the remark made in the beginning of this section, we also have H∗=H↪X∗. In this paper we also assume that the embedding of X into H is compact. Hence by Schauder’s theorem (see, for example, Gasinski and Papageorgiou [20, Theorem 3.1.22, p. 275]), so is the embedding of H∗=H into X∗. In what follows, by ∣∣⋅∣∣ (resp. ∣⋅∣,∣∣⋅∣∣∗) we denote the norm of the space X (resp. H,X∗). By ⟨⋅,⋅⟩ we denote the duality brackets for the pair (X∗,X) and by (⋅,⋅) the inner product of the Hilbert space H. We know that
[TABLE]
Also, let β>0 be such that
[TABLE]
We introduce the following space which has a central role in the study of the evolution inclusions. So, let 1<p<∞ and set
[TABLE]
In this definition the derivative of x is understood in the sense of vectorial distributions (weak derivative). In fact, if we view x as an X∗-valued function, then x(⋅) is absolutely continuous, hence strongly differentiable almost everywhere. Therefore
[TABLE]
The space Wp(0,b), equipped with the norm
[TABLE]
becomes a separable reflexive Banach space. We know that
[TABLE]
The following integration by parts formula is very helpful:
Proposition 2**.**
If x,y∈Wp(0,b), then t↦(x(t),y(t)) is absolutely continuous and
[TABLE]
We know that for all 1⩽p<∞,
[TABLE]
with p′=+∞ if p=1 (see Gasinski and Papageorgiou [20, Theorem 2.2.9, p. 129]).
Now, let (Ω,Σ) be a measurable space and V a separable Banach space. We introduce the following hyperspaces:
[TABLE]
Given a multifunction F:Ω→2V\{∅}, the “graph” of F is the set
[TABLE]
We say that F(⋅) is “graph measurable” if GrF∈Σ×B(V) with B(V) being the Borel σ-field of V. If μ(⋅) is a σ-finite measure on Σ and F:Ω→2V\{∅} is graph measurable, then the Yankov-von Neumann-Aumann selection theorem (see Hu and Papageorgiou [22, Theorem 2.14, p. 158]) implies that F(⋅) admits a measurable selection, that is, there exists a Σ-measurable function f:Ω→V such that f(ω)∈F(u)μ-almost everywhere. In fact, there is a whole sequence {fn}n⩾1 of such measurable selections such that F(ω)⊆{fn(ω)}μ-almost everywhere (see Hu and Papageorgiou [22, Proposition 2.17, p. 159]). Moreover, the above results are valid if V is only a Souslin space. Recall that a Souslin space need not be metrizable (see Gasinski and Papageorgiou [21, p. 232]). A multifunction F:ΩrightarrowPf(V) is said to be “measurable” if for all y∈V, the function
[TABLE]
is Σ-measurable. A multifunction F:Ω→Pf(V) which is measurable is also graph measurable. The converse is true if (Ω,Σ) admits a complete σ-finite measure μ. If (Ω,Σ,μ) is a σ-finite measure space and F:Ω→2V\{∅} is a multifunction, then for 1⩽p⩽∞ we introduce the set
[TABLE]
Evidently, SFP=∅ if and only if ω↦inf[∣∣v∣∣V:v∈F(ω)] belongs to Lp(Ω). Moreover, the set SFP is “decomposable”, that is, if (A,f1,f2)∈Σ×SFP×SFP, then
[TABLE]
Here, for C∈Σ, by χC we denote the characteristic function of the set C∈Σ.
For every D⊆Σ,D=∅, we define
[TABLE]
Here, ⟨⋅,⋅⟩V denotes the duality brackets of the pair (V∗,V). The function σ(⋅,D):V∗→Rˉ=R∪{+∞} is known as the “support function” of D.
Let Z,W be Hausdorff topological spaces. We say that a multifunction G:Z→2W\{∅} is “upper semicontinuous” (usc for short), respectively “lower semicontinuous” (lsc for short), if for all U⊆W open, the set
[TABLE]
is open in Z. If G(⋅) is both usc and lsc, then we say that G(⋅) is continuous. On a Hausdorff topological space (W,τ) (τ being the Hausdorff topology), we can define a new topology τseq whose closed sets are the sequentially τ-closed sets. Then topological properties with respect to this topology have the prefix “sequential”. Note that τ⊆τseq and the two are equal, if τ is first countable (see Buttazzo [7, p. 9] and Gasinski and Papageorgiou [21, p. 808]). We say that G:Z→2W\{∅} is “closed”, if the graph GrG⊆Z×W is closed.
For any Banach space V, on Pf(V) we can define a generalized metric, known as the “Hausdorff metric”, by setting
[TABLE]
Recall that (Pf(V),h) is a complete metric space (see Hu and Papageorgiou [22, p. 6]). If Z is a Hausdorff topological space, a multifunction G:Z→Pf(V) is said to be “h-continuous”, if it is continuous from Z into (Pf(V),h).
Also, if E,M⊆V are nonempty, bounded, closed and convex subsets, then
[TABLE]
(Hörmander’s formula).
Let (W,τ) be a Hausdorff topological space with topology τ and let {En}n⩾1⊆2W\{∅}. We define
[TABLE]
Sometimes we drop the Kseq-symbol and simply write τ−n→∞limsupEn and τ−n→∞liminfEn.
Returning to the setting of an evolution triple, we consider a sequence of multivalued maps an,a:Lp(T,X)→2Lp′(T,X∗)\{∅} (n∈N) such that for every h∗∈Lp′(T,X∗) the inclusions
[TABLE]
have unique solutions en(h∗),e(h∗)∈Wp(0,b).
We say that dtd+an “PG-converges” to dtd+a (denoted by dtd+anPGdtd+a as n→∞), if for every h∗∈Lp′(T,X∗) we have
[TABLE]
In what follows, by Xw (respectively Hw,Xw∗) we denote the space X (respectively H,X∗) furnished with the weak topology. Also, by ∣⋅∣1 we denote the Lebesgue measure on R and by ((⋅,⋅)) we denote the duality brackets for the pair (Lp′(T,X∗),Lp(T,X)). So, we have
[TABLE]
Next, let us recall some useful facts from the theory of nonlinear operators of monotone type.
So, let V be a reflexive Banach space, L:D(L)⊆V→V∗ a linear maximal monotone operator and a:V→2V∗. We say that a(⋅) is “L-pseudomonotone” if the following conditions hold:
(a)
For every v∈V,a(v)∈Pwkc(V∗).
(b)
a(⋅) is bounded (that is, maps bounded sets to bounded sets).
(c)
If {vn}n⩾1⊆D(L),vn→wv∈D(L) in V, L(vn)→wL(v) in V∗, vn∗∈a(vn) for all n∈N, vn∗→wv∗ in X∗ and n→∞limsup⟨vn∗,vn−v⟩V⩽0, then v∗∈a(v) and ⟨vn∗,vn⟩V→⟨v∗,v⟩V
Such maps have nice surjectivity properties. The next result is due to Papageorgiou, Papalini and Renzacci [35] and it extends an earlier single-valued result of Lions [27, Theorem 1.2, p. 319].
Proposition 3**.**
Assume that V is a reflexive Banach space which is strictly convex, L:D(L)⊆V→V∗ is a linear maximal monotone operator and A:V→2V∗ is L-pseudomonotone and strongly coercive, that is,
[TABLE]
Then R(L+V)=V∗ (that is, L+V is surjective).
In the next section we obtain some results about a general class of evolution inclusions, which will help us study problem (5) (see Section 4).
3. Nonlinear Evolution Inclusions
Let T=[0,b] and let (X,H,X∗) be an evolution triple with X↪H compactly (see Definition 1). In this section we deal with the following evolution inclusion:
[TABLE]
The hypotheses on the data of (9) are the following:
H(A)1:A:T×X→2X∗ is a map such that
(i)
for all x∈X,t↦A(t,x) is graph measurable;
(ii)
for almost all t∈T,GrA(t,⋅) is sequentially closed in Xw×Xw∗ and x↦A(t,x) is pseudomonotone;
(iii)
for almost all t∈T, all x∈X and all h∗∈A(t,x), we have
[TABLE]
with 2⩽p,a1∈Lp′(T) and c1>0;
(iv)
for almost all t∈T, all x∈X and all h∗∈A(t,x), we have
[TABLE]
with c2>0,a2∈L1(T).
Remark 1**.**
If A(⋅,⋅) is single-valued, then in hypothesis H(A)(ii) we can drop the condition on the graph of GrA(t,⋅) and only assume that for almost all t∈T,x↦A(t,x) is pseudomonotone. Similarly, if for almost all t∈T,A(t,⋅) is maximal monotone. An example of where the condition on the graph of A(t,⋅) is satisfied is the following. For simplicity we drop the t-dependence
[TABLE]
where φ:Lp(Ω,RN)→R is continuous, convex and ξ:Lp(Ω,RN)→R is continuous and ∣ξ(y)∣⩽c^(1+∣y∣τ−1) for all y∈RN, some c^>0 and with 1⩽τ<p. Then recalling that W1,p(Ω)↪W1,τ(Ω) compactly (see Zeidler [41, p. 1026]), we easily see that GrA is sequentially closed in W1,p(Ω)w×W1,p(Ω)w∗.
H(F)1:F:T×H→Pfc(H) is a multifunction such that
(i)
for all x∈H, t↦F(t,x) is graph measurable;
(ii)
for almost all t∈T,GrF(t,⋅) is sequentially closed in H×Hw;
(iii)
for almost all t∈T, all x∈H and all h∈F(t,x)
[TABLE]
with a3∈L2(T),c3>0 and if p=2, then β2c3<c2 (see (6)).
By a solution of problem (9) we understand a function x∈Wp(0,b) such that
[TABLE]
with h∗∈Lp′(T,X∗),f∈L2(T,H) such that
[TABLE]
By S(ξ) we denote the set of solutions of problem (9). In the sequel we investigate the structure of S(ξ).
Consider the multivalued map a:Lp(T,X)→2Lp′(T,X∗) defined by
[TABLE]
Note that a(⋅) has values in Pwkc(Lp′(T,X∗)) (see hypotheses H(A)(i),(iii) and use the Yankov-von Neumann-Aumann selection theorem (see Hu and Papageorgiou [22, Theorem 2.15, p. 518])).
Lemma 4**.**
If hypotheses H(A)1 hold, xn→wx in Wp(0,b), xn(t)→wx(t) in X for almost all t∈T, hn∗→wh∗ in Lp′(T,X∗) and hn∗∈a(xn) for all n∈N, then h∗∈a(x).
Proof.
Let v∈X and consider the function x↦σ(v,A(t,x)) (see Section 2). We will show that it is sequentially upper semicontinuous. To this end we need to show that given λ∈R, the superlevel set
[TABLE]
is sequentially closed in Xw. So, we consider a sequence {xn}n⩾1⊆Eλ and assume that
[TABLE]
Let hn∗∈A(t,xn) (n∈N) be such that
[TABLE]
(recall A(t,xn)∈Pwkc(X∗)). Evidently, {hn∗}n⩾1⊆X∗ is bounded (see hypothesis
H(A)1(iii)) and so by passing to a subsequence if necessary, we may assume that
[TABLE]
Then we have
[TABLE]
This proves the upper semicontinuity of the map x↦σ(v,A(t,x)).
Now let v∈Lp(T,X). We have
[TABLE]
∎
Lemma 5**.**
If hypotheses H(A)1 hold, then the multivalued map a:Lp(T,X)→2Lp′(T,X∗) defined by (10) is L-pseudomonotone.
Proof.
Suppose xn→wx in Wp(0,b),hn∗∈a(xn) for all n∈N, hn∗→wh∗ in Lp′(T,X∗) and
We set ϑn(t)=⟨hn∗(t),xn(t)−x(t)⟩. Let N be the Lebesgue-null set in T=[0,b] outside of which hypotheses H(A)1(ii),(iii),(iv) hold. Using hypotheses H(A)1(iii),(iv) we have
[TABLE]
We introduce the Lebesgue measurable set D⊆T defined by
[TABLE]
Suppose that ∣D∣1>0. If t∈D∩(T\N), then from (15) we see that
From the above proof it is clear why in the case of a single-valued map A(t,x), in hypothesis H(A)1(ii) we can drop the condition on the graph of A(t,⋅) and only assume that for almost all t∈Tx↦A(t,x) is pseudomonotone. Indeed, in this case, from (21) we have (at least for a subsequence) that
[TABLE]
In the multivalued case, there is no canonical way to identify the pointwise limit of the sequence {hn∗(t)}n⩾1⊆X∗. If for almost all t∈T,A(t,⋅) is maximal monotone, then again, we do not need the graph hypothesis on A(t,⋅). In this case a(⋅) is also maximal monotone and then the lemma is a consequence of (13) and Lemma 1.3, p. 42 of Barbu [4]. It is worth mentioning that a similar strengthening of the topology in the range space was used by Defranceschi [12], while studying G-convergence of multivalued operators.
Without loss of generality, invoking the Troyanski renorming theorem (see Gasinski and Papageorgiou [21, Remark 2.115, p. 241]), we may assume that both X and X∗ are locally uniformly convex, hence Lp(T,X) and Lp′(T,X∗) are strictly convex.
We are now ready for the first result concerning the solution set S(ξ).
Theorem 6**.**
If hypotheses H(A)1,H(F)1 hold and ξ∈H, then the solution set S(ξ) is nonempty, weakly compact in Wp(0,b) and compact in C(T,H).
Proof.
First suppose that ξ∈X. We define
[TABLE]
Evidently, A1(t,x) and F1(t,x) have the same measurability, continuity and growth properties as the multivalued maps A(t,x) and F(t,x). So, we may equivalently consider the following Cauchy problem
[TABLE]
Note that if x∈Wp(0,b) is a solution of (22), then x^=x−ξ is a solution of (9) (when ξ∈X, that is, the initial condition is regular). Consider the linear densely defined operator L:D(L)⊆Lp(T,X)→Lp′(T,X∗) defined by
[TABLE]
(the evaluation y(0)=0, makes sense by virtue of (7)).
Consider the multivalued maps a1,G1:Lp(T,X)→2Lp′(T,X∗)\{∅} defined by
[TABLE]
We set K(x)=a1(x)+G1(x) for all x∈Lp(T,X). Then
[TABLE]
Claim 1**.**
K* is L-pseudomonotone.*
Clearly, K has values in Pwkc(Lp′(T,X∗)) and it is bounded (see hypotheses H(A)(iii), H(F)(iii)).
Next, we consider a sequence {xn}n⩾1⊆D(L) such that
[TABLE]
Then we have
[TABLE]
Hypotheses H(A)(iii) and H(F)(iii) imply that
[TABLE]
So, we may assume (at least for a subsequence) that
It follows from (30) that K is coercive. This proves Claim 2.
Now Claims 1 and 2 permit the use of Proposition 3 to find x∈Wp(0,b) solving problem (9) when ξ∈X.
Next, we remove the restriction ξ∈X. So, suppose ξ∈H. We can find {ξn}n⩾1⊆X such that ξn→ξ in H (recall that X is dense in H). From the first part of the proof, we know that we can find xn∈S(ξn)⊆Wp(0,b) for all n∈N. We have
[TABLE]
It follows that
[TABLE]
We have
[TABLE]
From (31), (3) and hypotheses H(A)1(iii),H(F)1(iii) it follows that
Hypothesis H(F)1(iii) implies that {fn}n⩾1⊆L2(T,H) is bounded. Hence
[TABLE]
By hypothesis H(A)1(iii) we see that
[TABLE]
So, we may assume that
[TABLE]
From (33), (36), (37), we see that we can use Lemma 5 and infer that
[TABLE]
As we have already mentioned {fn}n⩾1⊆L2(T,H) is bounded and so we may assume that
[TABLE]
Using Proposition 6.6.33 on p. 521 of Papageorgiou and Kyritsi [34], we have
[TABLE]
In (31) we pass to the limit as n→∞ and use (33), (37), (39) to obtain
[TABLE]
So, we have proved that when ξ∈H, the solution set S(ξ) is a nonempty subset of Wp(0,b).
Next, we will prove the compactness of S(ξ) in Wp(0,b)w and in C(T,H). Let x∈S(ξ). For every t∈T we have
[TABLE]
Then let rM:H→H be the M-radial retraction defined by
[TABLE]
Because of the a priori bound (41), we can replace F(t,x) by
[TABLE]
Note that for all x∈H,t↦F^(t,x) is graph measurable (hence also measurable, see Section 2) and for almost all t∈T, x↦F^(t,x) has a graph which is sequentially closed in H×Hw. Moreover, we see that
[TABLE]
We introduce the set
[TABLE]
We consider the following Cauchy problem
[TABLE]
Let H:C→2C(T,H) be the map (in general, multivalued) that assigns to each f∈C the set of solutions of problem (42). It is a consequence of Proposition 3 and Lemma 5, that H(⋅) has nonempty values.
Returning to (47), passing to the limit as n→∞ and using (48), (49), (50) we obtain
[TABLE]
From the proof of Lemma 5 (see (21)), we know that
[TABLE]
In a similar fashion we also have
[TABLE]
Also, by (8), (3), (46) and Vitali’s theorem, we have
[TABLE]
For every t∈T and every n∈N, using Proposition 2, we have
[TABLE]
However, from the previous parts of the proof it is clear that S(ξ)⊆H(C) is weakly closed in Wp(0,b) and closed in C(T,H). Therefore we conclude that S(ξ) is weakly compact in Wp(0,b) and compact in C(T,H).
∎
Next, we want to produce a continuous selection of the multifunction ξ↦S(ξ) (we refer to Repovš and Semenov [36] for more details about continuous selections of multivalued mappings). Note that S(⋅) is in general not convex-valued and so the Michael selection theorem (see Hu and Papageorgiou [22, Theorem 4.6, p. 92]) cannot be used. To produce a continuous selection of the solution multifunction ξ↦S(ξ), we need to strengthen the conditions on the multimap A(t,⋅), in order to guarantee that certain Cauchy problems admit a unique solution.
The new hypotheses on the map A(t,x) are the following:
H(A)2:A:T×X→2X∗\{∅} is a multivalued map such that
(i)
for every x∈X,t↦A(t,x) is graph measurable;
(ii)
for almost all t∈T,x↦A(t,x) is maximal monotone;
(iii)
for almost all t∈T, all x∈X and all h∗∈A(t,x), we have
[TABLE]
with a1∈Lp′(T),c1>0,2⩽p;
(iv)
for almost all t∈T, all x∈X and h∗∈A(t,x), we have
[TABLE]
with c2>0,a2∈L1(T)+.
Remark 3**.**
As we have already mentioned in an earlier remark, since now for almost all t∈T, A(t,⋅) is maximal monotone, we do not need the condition on the graph of A(t,⋅) (see hypothesis H(A)1(ii) and Barbu [4, Lemma 1.3, p. 42]).
Also, we strengthen the condition on the multifunction F(t,⋅).
H(F)2:F:T×H→Pfc(H) is a multifunction such that
(i)
for every x∈H, t↦F(t,x) is graph measurable;
(ii)
for almost all t∈T and all x,y∈H we have
[TABLE]
with k∈L1(T)+;
(iii)
for almost all t∈T, all x∈H and all h∈F(t,x), we have
[TABLE]
with a3∈L2(T)+,c3>0 and if p=2, then β2c3<c2 (see (6)).
Remark 4**.**
Hypothesis H(F)2(ii) is stronger than condition H(F)1(ii). Indeed, suppose that H(F)2(ii) holds and we have
[TABLE]
By the definition of the Hausdorff metric (see Section 2), we have
[TABLE]
The function y↦d(y,F(t,x)) is continuous and convex, hence weakly lower semicontinuous. Therefore by (54) we have
[TABLE]
This proves that condition H(F)1(ii) holds.
So, we can use Theorem 6 and establish that given any ξ∈H, the solution set S(ξ) is nonempty, weakly compact in Wp(0,b) and compact in C(T,H). The next result extends an earlier result of Cellina and Ornelas [9] for differential inclusions in RN with A≡0.
Proposition 7**.**
If hypotheses H(A)2,H(F)2 hold, then there exists a continuous map ϑ:H→C(T,H) such that
[TABLE]
Proof.
Consider the following auxiliary Cauchy problem
[TABLE]
This problem has a unique solution x0(ξ)∈Wp(0,b) (see Proposition 3 and use the monotonicity of A(t,⋅) and Proposition 2 to check the uniqueness of this solution).
We consider the multifunction Γ0:H→Pwkc(L1(T,H)) defined by
[TABLE]
We have
[TABLE]
Also, Γ0(⋅) has decomposable values. So, we can apply the selection theorem of Bressan and Colombo [5] (see also Hu and Papageorgiou [22, Theorem 8.7, p. 245]) and find a continuous map γ0:H→L1(T,H) such that γ0(ξ)∈Γ0(ξ) for all ξ∈H. Evidently γ0(ξ)∈L2(T,H) for all ξ∈H.
We consider the following auxiliary Cauchy problem:
[TABLE]
This problem has a unique solution x1(ξ)∈Wp(0,b).
By induction we will produce two sequences
[TABLE]
which satisfy:
(a)
xn(ξ)∈Wp(0,b) is the unique solution of the Cauchy problem
[TABLE]
(b)
ξ↦γn(ξ) is continuous from H into C(T,H);
(c)
γn(ξ)(t)∈F(t,xn(ξ)(t)) for almost all t∈T, all ξ∈H;
(d)
∣γn(ξ)(t)−γn−1(ξ)(t)∣⩽k(t)βn(ξ)(t) for almost all t∈T, all ξ∈H, where
[TABLE]
with ϵ>0,η(ξ)(t)=a2(t)+c2∣x0(ξ)(t)∣,τ(t)=∫0tk(s)ds.
Note that the maps ξ↦η(ξ) and ξ↦βn(ξ) are continuous from H into L1(T).
So, suppose we have produced {xk(ξ)}k=1n and {γk(ξ)}k=1n (induction hypothesis). Let xn+1(ξ)∈Wp(0,b) be the unique solution of the Cauchy problem
By (3) and Lemma A.5, p. 157 of Brezis [6], we infer that
[TABLE]
Using the induction hypothesis (see (c)) and hypothesis H(F)2(ii), we have
[TABLE]
Consider the multifunction Γn+1:H→2L1(T,H) defined by
[TABLE]
By (62) and Lemma 8.3 on p. 239 of Hu and Papageorgiou [22], we have that
[TABLE]
We can apply the selection theorem of Bressan and Colombo [5, Theorem 3] to find a continuous map γn+1:H→L1(T,H) such that γn+1(ξ)∈Γn+1(ξ) for all ξ∈H.
This completes the induction process and we have produced two sequences {xn(ξ)}n⩾1, {γn(ξ)}n⩾1 which satisfy properties (a)→(d) stated earlier.
Evidently, ξ↦x(ξ) is continuous from H into C(T,H), while because of hypothesis H(F)2(iii), we have γn(ξ)→γ(ξ) in L2(T,H). Let x^(ξ)∈Wp(0,b) be the unique solution of
[TABLE]
As before, exploiting the monotonicity of A(t,⋅) (see hypothesis H(A)2(ii)), we have
[TABLE]
So, x(ξ)∈S(ξ) and the map ϑ:H→C(T,H) defined by ϑ(ξ)=x(ξ) is a continuous selection of the solution multifunction ξ↦S(ξ).
∎
An easy, but useful consequence of Proposition 7 and of its proof, is a parametric version of the Filippov-Gronwall inequality (see Aubin and Cellina [1, Theorem 1, pp. 120-121] and Frankowska [19]) for differential inclusions.
So, we consider the following parametric version of problem (9):
[TABLE]
The parameter space D is a complete metric space. The hypotheses on the parametric vector field F(t,x,λ) and the initial condition ξ(λ) are the following:
H(F)2′:F:T×H×D→Pfc(H) is a multifunction such that
(i)
for all (x,λ)∈H×D,t↦F(t,x,λ) is graph measurable;
(ii)
for almost all t∈T, all x,y∈H, all λ∈D, we have
[TABLE]
with k∈L1(T)+;
(iii)
for almost all t∈T, all x∈H, all λ∈D and all h∈F(t,x), we have
[TABLE]
with a3∈L2(T)+,c3>0 and if p=2, then β2c3<c2 (see (6));
(iv)
for almost all t∈T, all x∈H, the multifunction λ↦F(t,x,λ) is lsc.
H0: the mapping λ↦ξ(λ) is continuous from D into H.
Assume that λ↦(u(λ),h(λ)) is a continuous map from D into C(T,H)×L2(T,H). We can find a continuous map p:D→L2(T) such that
[TABLE]
(see hypothesis H(F)2′(iii)).
In what follows, by e(h,λ)∈Wp(0,b) we denote the unique solution of the Cauchy problem
[TABLE]
with h∈L2(T,H).
We have the following approximation result.
Proposition 8**.**
Assume that hypotheses H(A)2,H(F)2′,H0 hold, λ↦(u(λ),h(λ)) is a continuous map from D into C(T,H)×L2(T,H) with u(λ)=e(h(λ),λ), ϵ>0 and p:D→L2(T)+ is a continuous map such that
[TABLE]
Then there exists a continuous map λ↦(x(λ),f(λ)) from D into C(T,H)×L2(T,H) such that
[TABLE]
and ∣x(λ)(t)−u(λ)(t)∣⩽bϵeτ(t)+∫0tp(λ)(s)eτ(t)−τ(s)ds for all t∈T with τ(t)=∫0tk(s)ds.
Proof.
Consider the multifunction Rϵ:D→2L1(T,H) defined by
[TABLE]
This multifunction has nonempty, decomposable values and it is lsc (see Hu and Papageorgiou [22, Lemma 8.3, p. 239]). Hence λ↦Rϵ(λ) has the same properties. So, we can find a continuous map γ0:D→L1(T,H) such that
[TABLE]
Let x1(λ)∈Wp(0,b) be the unique solution of the following Cauchy problem
[TABLE]
Then as in the proof of Proposition 7, we can generate by induction two sequences
[TABLE]
satisfying properties (a)→(d) listed in the proof of Proposition 7.
As before (see the proof of Proposition 7), we have
[TABLE]
From this inequality and property (d) of the sequences (see the proof of Proposition 7), we infer that
[TABLE]
are both Cauchy uniformly in λ∈K⊆D compact (recall that λ↦p(λ) is continuous, hence locally bounded). So, we have
[TABLE]
and both maps D∋λ↦x^(λ)∈C(T,H) and D∋λ↦γ^(λ)∈L1(T,H) are continuous. Moreover, we have γ^(λ)∈SF(⋅,x^(λ)(⋅),λ)2 (see the proof of Theorem 6) and that λ↦γ^(λ) is continuous from D into L2(T,H). If x(λ)=e(γ^(λ),λ), then
[TABLE]
From the triangle inequality, we have
[TABLE]
Using property (d) (see the proof of Proposition 7), we have
[TABLE]
So, finally we can write that
[TABLE]
∎
We want to strengthen Proposition 7, and require that the selection ϑ(⋅) passes through a preassigned solution. We mention that an analogous result for differential inclusions in RN with A≡0, was proved by Cellina and Staicu [10].
We start with a simple technical lemma.
Lemma 9**.**
If {uk}k=0N⊆L1(T,H) and {Tk(ξ)}k=0N is a partition of T=[0,b] with endpoints which depend continuously on ξ∈H, then there exists d^∈L1(T)+ for which the following holds:
We set d^(t)=k=0∑N∣uk(t)∣∈L1(T)+. From the hypothesis concerning the partition {Tk(ξ)}k=0N of T, we see that given ϵ>0, we can find δ>0 such that
[TABLE]
with C⊆T measurable, ∣C∣1⩽ϵ. Then
[TABLE]
The proof is now complete.
∎
With this lemma, we can produce a continuous selection of the solution multifunction ξ↦S(ξ), which passes through a preassigned point.
Proposition 10**.**
If hypotheses H(A)2,H(F)2 hold, K⊆H is compact, ξ0∈K and v∈S(ξ0), then there exists a continuous map ψ:K→C(T,H) such that
[TABLE]
Proof.
Since v∈S(ξ0), we have
[TABLE]
with f∈SF(⋅,v(⋅))2. Given g∈L2(T,H), we consider the unique solution of the Cauchy problem
[TABLE]
In what follows, by e(g,ξ)∈Wp(0,b) we denote the unique solution of problem (68) and we set μ0(ξ)=e(f,ξ). An easy application of the Yankov-von Neumann-Aumann selection theorem (see Hu and Papageorgiou [22, Theorem 2.14, p. 158]), gives γ0(ξ)∈L2(T,H) such that
[TABLE]
Let ϑ>0 we define
[TABLE]
The family {Bδ(ξ)(ξ)}ξ∈K is an open cover of the compact set K. So, we can find {ξk}k=0N⊆K such that {Bδ(ξk)(ξk)}k=0N is a finite subcover of K. Let {ηk}k=0N be a locally Lipschitz partition of unity subordinated to the finite subcover. We define
[TABLE]
The endpoints in these intervals are continuous functions of ξ. We consider the following Cauchy problem
[TABLE]
Problem (69) has a unique solution μ1(ξ)∈Wp(0,b). Let
[TABLE]
Using Lemma 9, we can find d^∈L1(T)+ such that, for any given ϵ>0, we can find δ>0 for which we have
[TABLE]
with C⊆T measurable, ∣C∣1⩽ϵ. We have μ1(ξ′)=e(λ0(ξ′),ξ′). As before, exploiting the monotonicity of A(t,⋅) (see hypothesis H(A)2(ii)) and using Lemma A.5, p. 157, of Brezis [6], we have
[TABLE]
Let ϵ>0 be given. By the absolute continuity of the Lebesgue integral, we can find δ1>0 such that
with C1⊆T measurable, ∣C1∣1⩽δ1. So, returning to (71) and using (72) and (73), we see that
[TABLE]
Therefore ξ↦μ1(ξ) is continuous from H into C(T,H). Again, with an application of the Yankov-von Neumann-Aumann selection theorem, we obtain γ1(ξ)∈L2(T,H) such that
[TABLE]
As in the proof of Proposition 7, we produce inductively two sequences
[TABLE]
which satisfy the following properties:
(a)
μn(ξ)=e(λn−1(ξ),ξ) with λn−1(ξ)=k=0∑NχTk(ξ)γn−1(ξk)(t),γ−1(ξ)=f;
(b)
ξ↦μn(ξ) is continuous from K into C(T,H);
(c)
∣μn(ξ)(t)−μn−1(ξ)(t)∣⩽2n+2n!ϑ(∫0tk(s)ds)n for all ξ∈K;
(d)
γn(ξ)(t)∈F(t,μn(ξ)(t)) for almost all t∈T and
[TABLE]
So, as induction hypothesis, suppose that we have produced
[TABLE]
which satisfy properties (a)→(d) stated above. We set
[TABLE]
As above (see in the first part of the proof the argument concerning the map ξ↦μ1(ξ)), we can show that ξ↦μn+1(ξ) is continuous from K into C(T,H). Also, by the monotonicity of A(t,⋅) (see hypothesis H(A)2(ii) and Lemma A.5, p. 157, of Brezis [6]), we have
[TABLE]
Moreover, a standard measurable selection argument, produces a measurable map
γn+1(ξ):T→H,ξ∈K, such that
[TABLE]
This completes the induction process.
Note that
[TABLE]
Therefore, we can say that
[TABLE]
It follows that ξ↦ψ(ξ) is continuous from K into C(T,H).
Note that T0(ξ0)=T=[0,b] and so μ0(ξ0)=e(f,ξ0)=v (see (67)). Hence ψ(ξ0)=v. It remains to show that ψ is a selection of the solution multifunction ξ↦S(ξ). By property (d) and hypothesis H(F)2(ii), we have
[TABLE]
Let μ^(ξ)=e(k=0∑N∫0tχTk(ξ)(s)γ^(ξk)(s)ds,ξ). Because
[TABLE]
we have
[TABLE]
∎
4. Optimal Control Problems
In this section we deal with the sensitivity analysis of the optimal control problem (5).
Let Q(ξ,λ)⊆Wp(0,b)×L2(T,Y) be the admissible “state-control” pairs. First we investigate the dependence of this set on the initial condition ξ∈H and the parameter λ∈E. Recall that the control space Y is a separable reflexive Banach space and the parameter space E is a compact metric space. To have a useful result on the dependence of Q(ξ,λ) on (ξ,λ)∈H×E, we introduce the following conditions on the data of the evolution inclusion in problem (5) (the dynamical constraint of the problem).
H(A)3:A:T×X×E→2X∗\{∅} is a multifunction such that
(i)
for every (x,λ)∈X×E,t↦Aλ(t,x) is graph measurable;
(ii)
for almost all t∈T, all λ∈E, x↦Aλ(t,x) is maximal monotone;
(iii)
for almost all t∈T, all x∈X, all λ∈E and all h∗∈Aλ(t,x), we have
[TABLE]
with {aλ}λ∈E⊆Lp′(T) bounded {cλ}λ∈E⊆(0,+∞) bounded and 2⩽p<∞;
(iv)
for almost all t∈T, all x∈X, all λ∈E and all h∗∈Aλ(t,x), we have
[TABLE]
with c^>0,a^∈L1(T)+;
(v)
if λn→λ in E, then dtd+aλnPGdtd+aλ as n→∞.
Hypotheses H(A)3(i)→(iv) are the same as hypotheses H(A)2(i)→(iv) for every map Aλ, λ∈E. The new condition is hypothesis H(A)3(v), which requires elaboration. In the examples that follow, we present characteristic situations where this hypothesis is satisfied.
Example 11**.**
(a)* First, we present a situation which will be used in Section 5.*
So, let Ω⊆RN be a bounded domain with Lipschitz boundary ∂Ω. Let X=W01,p(Ω)(2⩽p<∞), H=L2(Ω), X∗=W−1,p′(Ω). Evidently, (X,H,X∗) is an evolution triple (see Definition 1), with compact embeddings. We consider a map a(t,z,ξ) satisfying the following conditions:
H(a):* a:T×Ω×RN→RN is a map such that*
(i)
∣a(t,z,0)∣⩽c0* for almost all (t,z)∈T×Ω;*
(ii)
for every ξ∈RN,(t,z)↦a(t,z,ξ) is measurable;
(iii)
∣a(t,z,ξ1)−a(t,z,ξ2)∣⩽c^(1+∣ξ1∣+∣ξ2∣)p−1−α∣ξ1−ξ2∣α* for almost all (t,z)∈T×Ω, all ξ1,ξ2∈RN, with c^1>0, α∈(0,1];*
(iv)
(a(t,z,ξ1)−a(t,z,ξ2),ξ1−ξ2)RN⩾c^2∣ξ1−ξ2∣p* for almost all (t,z)∈T×Ω, all ξ1,ξ2∈RN, ξ1=ξ2, with c^2>0.*
We consider the operator A:T×X→X∗ defined by
[TABLE]
Using the nonlinear Green’s identity (see Gasinski and Papageorgiou [20, p. 210]), we have
[TABLE]
with B(t,x)(⋅)=a(t,⋅,Dx(⋅))∈Lp′(Ω,RN) for all (t,x)∈T×X.
Now consider a sequence {an(t,z,ξ)}n⩾1 of such maps satisfying
[TABLE]
with ϑ:R+→R+ being an increasing function which is continuous at r=0 and ϑ(0)=0. We assume that for almost all t∈T, an(t,⋅,⋅)→Ga(t,⋅,⋅) in the sense of Defranceschi [12]. By Svanstedt [40] we have
[TABLE]
(b)* We can allow multivalued maps, provided that we drop the t-dependence. So, we consider multivalued maps a(z,ξ) which satisfy the following conditions:*
H(a)′:* a:Ω×RN→2RN\{∅} is a measurable map such that*
(i)
a(⋅,⋅)* is measurable;*
(ii)
for almost all z∈Ω,ξ↦a(z,ξ) is maximal monotone;
(iii)
for almost all z∈Ω, all ξ∈RN and all y∈a(z,ξ), we have
[TABLE]
We again consider the evolution triple
[TABLE]
and consider the multivalued map A:X→2X∗\{∅} defined by
[TABLE]
We consider a sequence {an(z,ξ)}n⩾1 of such maps and assume that an→Ga in the sense of Defranceschi [12]. Then by Denkowski, Migorski and Papageorgiou [14] we have
[TABLE]
(c)* A third situation leading to hypothesis H(A)3(v) is the following one. We consider maps Aλ(t,x) satisfying the following conditions:*
H(A)3′:* A:T×X×E→X∗ is a map such that*
(i)
∣∣Aλ(t+τ,x)−Aλ(t,x)∣∣⩽O(τ)(1+∣∣x∣∣p−1)* for all t,t+τ∈T, all x∈X, all λ∈E;*
(ii)
for all (t,λ)∈T×E,x↦Aλ(t,x) is hemicontinous;
(iii)
for all t∈T, all x,u∈X, all λ∈E, we have
[TABLE]
(iv)
if λn→λ in E, then for all t∈T, Aλn(t,⋅)→GAλ(t,⋅) (this means that for all x∗∈X∗, Aλn−1(t,x∗)→wAλ−1(t,x), see Denkowski, Migorski and Papageorgiou **[16, Definition 3.8.20, p. 478]**).
Next, we introduce the conditions on the multifunctions F and G involved in the dynamics of (5).
H(F)3:F:T×H×E→Pfc(H) is a multifunction such that
(i)
for all (x,λ)∈H×E, t↦F(t,x,λ) is graph measurable;
(ii)
for almost all t∈T, all x,y∈H and all λ∈E, we have
[TABLE]
(iii)
for almost all t∈T, all x∈H and all λ∈E, we have
[TABLE]
with {aλ}λ∈E⊆L2(T) and {cλ}λ∈E⊆(0,+∞) bounded;
(iv)
for almost all t∈T, all x∈H and all λ,λ′∈E, we have
[TABLE]
with β(r)→0+ as r→0+ and w(t,⋅) bounded on bounded sets.
H(G):G:T×Y×E→Pfc(H) is a multifunction such that
(i)
for all (u,λ)∈Y×E,t↦G(t,u,λ) is graph measurable;
(ii)
for almost all t∈T, all λ∈E, u↦G(t,u,λ) is concave (that is, GrG(t,⋅,λ)⊆Y×H is concave, see Hu and Papageorgiou [22], Definition 1.1 and Remark 1.2, p. 585) and (u,λ)↦G(t,u,λ) is h-continuous;
(iii)
for almost all t∈T, all u∈U(t,λ), all λ∈E
[TABLE]
with {a^λ}λ∈E⊆L2(T) bounded.
Remark 5**.**
A typical situation resulting to a concave multifunction u↦G(t,u,λ), is when
[TABLE]
with Bλ(t)∈L(Y,H) and C(t,λ)∈Pfc(H) for all (t,λ)∈T×E.
Another situation, leading to the concavity of G(t,⋅,λ), is when H is an ordered Hilbert space and gλ,g~λ:T×Y→H are two Carathéodory maps such that for almost all t∈T
[TABLE]
We set G(t,u,λ)={h∈H:gλ(t,u)⩽h⩽g~λ(t,u)}. Then G(t,⋅,λ) is concave.
Finally we impose conditions on the control constraint U(t,λ).
H(U):U:T×E→Pfc(Y) is a multifunction such that
(i)
for all λ∈E,t↦U(t,λ) is graph measurable;
(ii)
for almost all t∈T,λ↦U(t,λ) is h-continuous;
(iii)
∣U(t,λ)∣⩽a~λ(t) for almost all t∈T, all λ∈E, with {a~λ}λ∈E⊆L2(T) bounded.
Proposition 12**.**
If hypotheses H(A)3,H(F)3,H(G),H(U) hold and (ξn,λn)→(ξ,λ) in H×E, then
[TABLE]
Proof.
Let (x,u)∈Kseq(s×w)−n→∞limsupQ(ξn,λn). By definition (see Section 2), we can find a subsequence {m} of {n} and (xm,um)∈Q(ξm,λm), m∈N such that
[TABLE]
For every m∈N, we have
[TABLE]
with fm,gm∈L2(T,H) such that
[TABLE]
We deduce by hypotheses H(F)3(iii),H(G)(iii) and Theorem 6 and its proof that
By (78) and hypotheses H(F)3(iii),H(G)(iii) it is clear that
[TABLE]
Hence, we may assume (at least for a subsequence), that
[TABLE]
Proposition 6.6.33 on p. 521 of Papageorgiou and Kyritsi [34], implies that
[TABLE]
Fix t∈T\N and let y∈w−m→∞limsupF(t,xm(t),λm). By definition, we know that there exists a subsequence {k} of {m} and yk∈F(t,xk(t),λk) for all k∈N such that yk→wy in H as k→∞. The function v↦d(v,F(t,x(t),λ)) is continuous and convex, hence weakly lower semicontinuous. Therefore
Let h∈L2(T,H) and let (⋅,⋅)L2(T,H) denote the inner product of L2(T,H) (recall that L2(T,H)∗=L2(T,H)). Then
[TABLE]
The concavity of G(t,⋅,λ) (see hypothesis H(G)(ii)), implies that the function u↦σ(h(t),G(t,u,λ)) is concave. Since E is a complete metric space, it can be isometrically embedded, by the Arens-Eells theorem (see Gasinski and Papageorgiou [21, Theorem 4.143, p. 655]), as a closed subset of a separable Banach space (recall that E is compact). So, by Balder [3], we have
[TABLE]
Since h∈L2(T,H) is arbitrary, it follows that
[TABLE]
Let ym∈Wp(0,b) be the unique solution of the Cauchy problem
[TABLE]
Hypothesis H(A)3(v) implies that
[TABLE]
with y∈Wp(0,b) being the unique solution of the Cauchy problem
[TABLE]
(see Section 2). From (77) and (85) and the monotonicity of Aλm(t,⋅) (see hypothesis H(A)3(ii)), we have
Next, we will prove the second convergence of the proposition.
So, let (x,u)∈Q(ξ,λ). By definition we have
[TABLE]
with g∈L2(T,H) satisfying
[TABLE]
For every v∈L2(T,Y), we have
[TABLE]
Hypothesis H(U)(ii) and the dominated convergence theorem imply that
[TABLE]
Hence Proposition 6.6.22 on p. 518 of Papageorgiou and Kyritsi [34] implies that we can find un∈SU(⋅,λn)2 (n∈N) such that
[TABLE]
Then hypothesis H(G)(ii) guarantees that we can find
[TABLE]
such that
[TABLE]
Given ξ′∈H, let S(ξ′)⊆Wp(0,b) be the set of solutions of the Cauchy problem
[TABLE]
Let K={ξn,ξ}n⩾1⊆H. This is a compact set in H. Invoking Proposition 10 (with ξ0=ξ), we produce a continuous map ψ:K→C(T,H) such that
[TABLE]
Let yn=ψ(ξn) (n∈N) and use Proposition 8 to find xn∈Wp(0,b) solution of the Cauchy problem
[TABLE]
for which we have
[TABLE]
with ϵ>0, τ(t)=∫0tk(s)ds,ηn∈L1(T),ηn→0 in L1(T).
So, we obtain
[TABLE]
Since ϵ>0 is arbitrary, it follows that
[TABLE]
Finally, we have
[TABLE]
Since (xn,un)∈Q(ξn,λn) (n∈N) and un→u in L2(T,Y), we conclude that
[TABLE]
∎
An immediate consequence of the above proposition is the following corollary concerning the multifunction (ξ,λ)↦Q(ξ,λ) of admissible state-control pairs.
Corollary 13**.**
If hypotheses H(A)3,H(F)3,H(G),H(U) hold, then the multifunction Q:H×E→2C(T,H)×L2(T,Y)\{∅} is lsc and sequentially closed in C(T,H)×L2(T,Y)w (that is, GrQ⊆H×E×C(T,H)×L2(T,Y)w is sequentially closed).
Now we bring the cost functional into the picture. The hypotheses on the integrands L(t,x,λ) and H(t,u,λ) are the following.
H(L):L:T×H×E→R is an integrand such that
(i)
for every (x,λ)∈H×E, t↦L(t,x,λ) is measurable;
(ii)
if λn→λ in E, then for all x∈H we have L(⋅,x,λn)→wL(⋅,x,λ) in L1(T);
(iii)
for almost all t∈T, all x,y∈H and all λ∈E, we have
[TABLE]
where ∣x∣∨∣y∣=max{∣x∣,∣y∣} and ρ(t,r) is a Carathéodory function on T×R+ with values in (0,+∞) such that
[TABLE]
with βϑ∈L1(T)+,ϑ>0.
H(H):H:T×Y×E→R is an integrand such that
(i)
for all (u,λ)∈Y×E,t↦H(t,u,λ) is measurable;
(ii)
for almost all t∈T, all v∈E, u↦H(t,u,λ) is convex and for almost all t∈T, all u∈Y, λ↦H(t,u,λ) is continuous;
(iii)
for almost all t∈T and all (u,λ)∈Y×E, we have
[TABLE]
with a∈L∞(T).
H(ψ^):ψ^:H×E→R is a continuous function.
Using the direct method of the calculus of variations, we can produce optimal admissible state-control pairs for problem (5).
Proposition 14**.**
If hypotheses H(A)3,H(F)3,H(G),H(U),H(L),H(H) and H(ψ^) hold, then for every (ξ,λ)∈H×E we can find (x∗,u∗)∈Q(ξ,λ) such that
[TABLE]
Proof.
Let {(xn,un)}n⩾1⊆Q(ξ,λ) be a minimizing sequence for problem (5). So, we have
is relatively w×w-compact (respectively, s×w-compact). So, by the Eberlein-Smulian theorem and by passing to a suitable subsequence if necessary, we can say that
Next, we estimate the second term on the right-hand side of (4).
Let ϑ>2∣∣x∣∣C(T,H) and let βϑ∈L1(T)+ as postulated by hypothesis H(L)(iii). Given ϵ>0, we can find δ>0 such that
“if C⊆T is measurable with ∣C∣1⩽δ,
[TABLE]
Here, we use the absolute continuity of the Lebesgue integral. Invoking the Scorza-Dragoni theorem (see Papageorgiou and Kyritsi [34, Theorem 6.2.9, p. 471]), we can find T1⊆T closed with ∣T\T1∣⩽2δ and ρ∣T1×R+ is continuous. Since ρ(t,0)=0, we can find δ1>0 such that
[TABLE]
Recall that simple functions are dense in Lp(T,H). Using this fact, the property that Lp(T,H)-convergence implies pointwise convergence for almost all t∈T for at least a subsequence and invoking Egorov’s theorem, we can find T2⊆T closed and s:T→H a simple function such that
[TABLE]
We set T3=T1∩T2. This is a closed subset of T with ∣T\T3∣1⩽δ. We have
[TABLE]
Similarly, we show that
[TABLE]
Let s(t)=k=1∑NvkχCk(t) with vk∈H,Ck⊆T measurable. Using hypothesis H(L)(ii), we can find n0∈N such that
From Proposition 14 we know that for every n∈N, we can find (xn,un)∈Q(ξn,λn) such that
[TABLE]
As in the proof of Theorem 6, we can show that {xn}n⩾1⊆Wp(0,b) is bounded. In addition, hypothesis H(U) implies that {un}n⩾1⊆L2(T,Y) is bounded. So, by passing to a suitable subsequence if necessary, we may assume that
For every (ξ,λ)∈H×E, we introduce the set Σ(ξ,λ) of optimal state-control pairs, that is,
[TABLE]
By Proposition 14, we know that for every (ξ,λ)∈H×E, Σ(ξ,λ)=∅. For this multifunction we can prove the following useful continuity property.
Theorem 16**.**
If hypotheses H(A)3,H(F)3,H(U),H(L),H(H) and H(ψ^) hold, then the multifunction Σ:H×E→2C(T,H)×L2(T,Y)\{∅} is sequentially usc into C(T,H)×L2(T,Y)w.
Proof.
Let C⊆C(T,H)×L2(T,Y)w be sequentially closed. We need to show that
[TABLE]
is closed in H×E (see Section 2). To this end, let {(ξn,λn)}n⩾1⊆Σ−(C) and assume that
[TABLE]
Let (xn,un)∈Σ(ξn,λn)∩C, n∈N. We know from the proof of Theorem 15 that at least for a subsequence, we have
Moreover, from (115) and since C⊆C(T,H)×L2(T,Y)w is sequentially closed, we deduce that (x,u)∈Σ(ξ,λ)∩C. Therefore Σ−(C)⊆H×E is closed and this proves the desired sequential upper semicontinuity of the multifunction (ξ,λ)↦Σ(ξ,λ).
∎
5. Application to Distributed Parameter Systems
In this section we present an application to a class of multivalued parabolic optimal control problems.
So, let T=[0,b] and let Ω⊆RN be a bounded domain with a Lipschitz boundary ∂Ω. We examine the following nonlinear, multivalued parabolic optimal control problem:
[TABLE]
Here, aλ:Ω×RN→2RN (λ∈E) is a family of multifunctions as in Example 11(b). For the other data of problem (123), we introduce the following conditions:
H(F1):F1:T×Ω×R×E→Pfc(R) is a multifunction such that
(i)
for all (x,λ)∈R×E,(t,z)↦F1(t,z,x,λ) is measurable;
(ii)
for almost all (t,z)∈T×Ω, all x,y∈R, all λ∈E, we have
[TABLE]
with k1∈L1(T,L∞(Ω));
(iii)
for almost all (t,z)∈T×Ω, all x∈R, all λ∈E, we have
[TABLE]
with a^1∈L2(T×Ω),c^1>0;
(iv)
for almost all (t,z)∈T×Ω, all x∈R, all λ,λ′∈E, we have
[TABLE]
with β(r)→0 as r→0+ and w∈Lloc∞(Ω×R+).
Remark 6**.**
Consider the multifunction F(t,z,x,λ) defined by
[TABLE]
with f,f^:T×Ω×R×E→R two functions such that
•
for all (x,λ)∈R×E, (t,z)↦f(t,z,x,λ),f^(t,z,x,λ) are both measurable;
•
for almost all (t,z)∈T×Ω, all x,x′∈R, all λ,λ′∈E, we have
[TABLE]
with k,k^∈L1(T,L∞(Ω)).
Then this multifunction satisfies hypotheses H(F1).
H(g):g:T×Ω×E→R is a Carathéodory function (that is, for all λ∈E, (t,z)→g(t,z,λ) is measurable and for almost all (t,z)∈T×Ω, λ→g(t,z,λ) is continuous) and for almost all (t,z)∈T×Ω and all λ∈E, we have ∣g(t,z,λ)∣⩽M with M>0.
H(r):r:T×E→R+ is a Carathéodory function (that is, for all λ∈E, t↦r(t,λ) is measurable and for almost all t∈T, λ→r(t,λ) is continuous) and for almost all t∈T, all λ∈E, we have
[TABLE]
with a∈L2(T).
Now, we introduce the conditions on the two integrands involved in the cost functional problem (123).
H(L1):L:T×Ω×R×E→R is an integrand such that
(i)
for all (x,λ)∈R×E,(t,z)↦L1(t,z,x,λ) is measurable;
(ii)
if λn→λ in E, then for all x∈L2(Ω) we have L1(⋅,⋅,x(⋅),λn)→wL1(⋅,⋅,x(⋅),λ) in L1(T×Ω);
(iii)
for almost all (t,z)∈T×Ω, all x,y∈R, all λ∈E
[TABLE]
with ρ(t,z,r) Carathéodory, ρ(t,z,0)=0 for almost all (t,z)∈T×Ω and for almost all (t,z), all r∈[0,ϑ] we have
[TABLE]
with βϑ∈L1(T×Ω).
H(H)1:H1:T×Ω×R×E→R is an integrand such that
(i)
for all (x,λ)∈R×E, (t,z)↦H1(t,z,x,λ) is measurable;
(ii)
for almost all (t,z)∈T×Ω, u↦H1(t,z,u,λ) is convex for all λ∈E, while λ↦H1(t,z,u,λ) is continuous for all u∈R;
(iii)
for almost all (t,z)∈T×Ω, all ∣u∣⩽rλ(t,z), all λ∈E, we have
[TABLE]
with {a^λ}λ∈E⊆L2(T×Ω) bounded.
We consider the following evolution triple:
[TABLE]
Since 2⩽p<∞, the Sobolev embedding theorem implies that in this triple the embeddings are compact.
For every λ∈E, let Aλ:X→2X∗\{∅} be the multivalued map defined by
[TABLE]
This map is maximal monotone and if λn→λ in E, then
[TABLE]
(see Example 11(b)). So, hypotheses H(A)3 hold. In fact, we can have t-dependence at the expense of assuming that aλ is single-valued. So, we assume that aλ(t,z,ξ) satisfies the conditions of Example 11 (a). Then the map Aλ:T×X→X∗ is defined by
[TABLE]
In fact, by the nonlinear Green’s identity (see Gasinski and Papageorgiou [20, p. 210]), we have
[TABLE]
As we have already mentioned in Example 11(a), we know from Svanstedt [40] that if λn→λ in E, then
[TABLE]
and so hypotheses H(A)3 hold.
As a special case of interest, we consider the situation where the elliptic differential operator is a weighted p-Laplacian, that is,
[TABLE]
Here, for every λ∈E,aλ:T×Ω→R is a measurable function such that
•
0<c^1⩽aλ(t,z)⩽c^2 for almost all (t,z)∈T×Ω, all λ∈E;
•
if λn→λ in E, then for almost all t∈T,
[TABLE]
For this case we consider the following parametric (with parameter λ∈E) family of convex (in ξ∈RN) integrands:
[TABLE]
Then the convex conjugate of φλ(t,z,⋅) is given by
[TABLE]
By hypothesis we have that
[TABLE]
We introduce the integral functional Φλ defined by
[TABLE]
By Marcellini and Sbordone [29], we know that (124) implies
[TABLE]
with Γseq(w) denoting the sequential Γ-convergence of Φλn(t,⋅) on W01,p(Ω)w (see Buttazzo [7]). Then it follows from Defranceschi [12, Theorem 3.3] that
Then hypotheses H(F1),H(g),H(r) imply that conditions H(F)3,H(G),H(U) hold. So, the dynamics of (123) are described by an evolution inclusion similar to the one in problem (5).
Finally let
[TABLE]
Hypotheses H(L1),H(H1) imply that conditions H(L),H(H) respectively hold.
So, we can apply Theorems 15 and 16 and obtain the following result concerning the variational stability of problem (123).
Proposition 17**.**
If the maps aλ are as above and hypotheses H(F1),H(g),H(r),H(L1), H(H1) hold, then for every (ξ,λ)∈L2(Ω)×E, problem (123) admits optimal pairs (that is, Σ(ξ,λ)=∅) and
[TABLE]
Acknowledgements. This research was supported in part by the SRA grants P1-0292-0101, J1-6721-0101 and J1-7025-0101.
Bibliography41
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] J.-P. Aubin, A. Cellina, Differential Inclusions , Springer-Verlag, Berlin (1984).
2[2] J.-P. Aubin, H. Frankowska, Set Valued Analysis , Birkhäuser, Boston (1990).
3[3] E. Balder, Necessary and sufficient conditions for L 1 subscript 𝐿 1 \displaystyle L_{1} -strong-weak lower semicontinuity of integral functionals, Nonlinear Anal. 11 (1987), 1399-1404.
4[4] V. Barbu, Nonlinear Semigroups and Differential Equations in Banach Spaces , Noordhoff International, Leyden, The Netherlands (1976).
5[5] A. Bressan, G. Colombo, Extensions and selections of maps with decomposable values, Studia Math. 90 (1988), 69-85.
6[6] H. Brezis, Opérateurs Maximaux Monotones et Semi-Groupes de Contractions dans les Espaces de Hilbert , North Holland, Amsterdam (1973).
7[7] G. Buttazzo, Semicontinuity, Relaxation and Integral Representation in the Calculus of Variations , Pitman Research Notes in Mathematics Series, Vol. 207, Longman Scientific & Technical, Harlow; copublished in the United States with John Wiley & Sons, Inc., New York (1989).
8[8] G. Buttazzo, G. Dal Maso, Γ Γ \displaystyle\Gamma -convergence and optimal control problems, J. Optim. Theory. Appl. 38 (1982), 385-407.