Optimal Control of Partial Differential Equations

Optimal Control of Partial Differential Equations Optimal Control of Partial Differential Equations

math.univ.toulouse.fr
from math.univ.toulouse.fr More from this publisher
29.03.2013 Views

68 CHAPTER 6. CONTROL OF THE WAVE EQUATION Let us establish the optimality conditions for (P2). As usual we set F2(u) = J2(z(u), u), where z(u) is the solution to (6.3.6). We have F ′ 2(ū2)u = Ω (∇¯z2(T ) − ∇zd(T ))∇wu(T ) + β Q χωū2u, where wu is the solution to Since ∂2w ∂t2 − ∆w = χωu in Q, w = 0 on Σ, w(x, 0) = 0 and ∂w (x, 0) = 0 in Ω. ∂t (∇¯z2(T ) − ∇zd(T ))∇wu(T ) = 〈−∆(¯z2(T ) − zd(T )), wu(T )〉 H−1 (Ω)×H1 0 (Ω), Ω applying formula (6.2.3) to p2 and wu, we obtain F ′ 2(ū2)u = (χω(βū2 + p2)u = 0 Q for every u ∈ L 2 (0, T ; L 2 (ω)). Thus the optimality condition for (P2) is proved. The proof of the other results is left to the reader. Comments. As for the heat equation with distributed controls, equation (6.3.6) is of the form y ′ = Ay + F + Bu, y(0) = y0, with Ay = A y1 y2 = y2 ∆y1 0 , F = f z0 , y0 = z1 0 , and Bu = χωu Thus problem (P1) is a particular case of control problems studied in Chapter 7. 6.4 Neumann boundary control We first study the equation ∂2z − ∆z = f in Q, ∂t2 ∂z ∂n = 0 on Σ, z(x, 0) = z0 and ∂z ∂t (x, 0) = z1 in Ω. (6.4.10) We set D(A) = {y1 ∈ H 2 (Ω) | ∂y1 ∂n = 0} × H1 (Ω), Y = H 1 (Ω) × L 2 (Ω), and Ay = A y1 y2 y2 = ∆y1 − y1 , Ly = Equation (6.4.10) may be written in the form dy dt 0 y1 , F = = (A + L)y + F, y(0) = y0. 0 f , and y0 = Theorem 6.4.1 The operator (A, D(A)) is the infinitesimal generator of a strongly continuous semigroup of contractions on Y . Proof. We leave the reader adapt the proof of Theorem 6.2.1. z0 z1 . .

6.4. NEUMANN BOUNDARY CONTROL 69 Theorem 6.4.2 For every f ∈ L 2 (Q), every z0 ∈ H 1 (Ω), every z1 ∈ L 2 (Ω), equation (6.4.10) admits a unique weak solution z(f, z0, z1), moreover the operator (f, z0, z1) ↦→ z(f, z0, z1) is linear and continuous from L 2 (Q)×H 1 (Ω)×L 2 (Ω) into C([0, T ]; H 1 (Ω))∩C 1 ([0, T ]; L 2 (Ω)). To study the wave equation with nonhomogeneous boundary conditions, we set D( A) = H1 (Ω) × L2 (Ω), Y = L2 (Ω) × (H1 (Ω)) ′ , and Ay = y1 y2 A = , y2 Ãy1 − y1 where Ãy1, ζ (H1 (Ω)) ′ = − ∇y1 · ∇(−∆ + I) Ω −1 ζ. Theorem 6.4.3 The operator ( A, D( A)) is the infinitesimal generator of a semigroup of contractions on Y . Proof. The proof is similar to the one of Theorem 6.2.3. Now, we consider the wave equation with a control in a Neumann boundary condition: ∂2z ∂z − ∆z = f in Q, ∂t2 ∂n = u on Σ, z(x, 0) = z0 and ∂z ∂t (x, 0) = z1 in Ω. (6.4.11) For any u ∈ L2 (Γ), the mapping ζ ↦→ Γ uζ is a continuous linear on H1 (Ω). Thus it can be identified with an element of (H1 (Ω)) ′ . Thus for u ∈ L2 (Σ), the mapping ζ ↦→ u(·)ζ is an element of L2 (0, T ; (H1 (Ω)) ′ ). Let us denote this mapping by û. We set y1 0 V = 0 û , F = 0 f , L Equation (6.4.11) may be written in the form y2 = y1 dy dt = ( A + L)y + F + V, y(0) = y0, , and y0 = z0 with F and V belong to L 2 (0, T ; L 2 (Ω)) × L 2 (0, T ; (H 1 (Ω)) ′ ), y0 ∈ L 2 (Ω) × (H 1 (Ω)) ′ . Theorem 6.4.4 For every (f, u, z0, z1) ∈ L 2 (Q)×L 2 (Σ)×L 2 (Ω)×(H 1 (Ω)) ′ , equation (6.4.11) admits a unique weak solution z(f, u, z0, z1) in C([0, T ]; L 2 (Ω))∩C 1 ([0, T ]; (H 1 (Ω)) ′ ). Moreover the mapping (f, u, z0, z1) ↦→ z(f, u, z0, z1) is continuous from L 2 (Q)×L 2 (Σ)×L 2 (Ω)×(H 1 (Ω)) ′ into C([0, T ]; L 2 (Ω)) ∩ C 1 ([0, T ]; (H 1 (Ω)) ′ ). Proof. The result is a direct consequence of Theorem 6.4.3. We consider the control problem (P4) inf{J4(z, u) | (z, u) ∈ C([0, T ]; L 2 (Ω)) × L 2 (0, T ; L 2 (Ω)), (z, u) satisfies (6.4.11)}, with J4(z, u) = 1 (z − zd) 2 Q 2 + 1 (z(T ) − zd(T )) 2 Ω 2 + β u 2 Σ 2 . z1 Γ .

6.4. NEUMANN BOUNDARY CONTROL 69<br />

Theorem 6.4.2 For every f ∈ L 2 (Q), every z0 ∈ H 1 (Ω), every z1 ∈ L 2 (Ω), equation (6.4.10)<br />

admits a unique weak solution z(f, z0, z1), moreover the operator<br />

(f, z0, z1) ↦→ z(f, z0, z1)<br />

is linear and continuous from L 2 (Q)×H 1 (Ω)×L 2 (Ω) into C([0, T ]; H 1 (Ω))∩C 1 ([0, T ]; L 2 (Ω)).<br />

To study the wave equation with nonhomogeneous boundary conditions, we set D( A) =<br />

H1 (Ω) × L2 (Ω), Y = L2 (Ω) × (H1 (Ω)) ′ , and<br />

Ay = <br />

y1<br />

y2<br />

A =<br />

,<br />

y2 Ãy1 − y1<br />

where <br />

Ãy1, ζ<br />

(H1 (Ω)) ′<br />

<br />

= − ∇y1 · ∇(−∆ + I)<br />

Ω<br />

−1 ζ.<br />

Theorem 6.4.3 The operator ( A, D( A)) is the infinitesimal generator <strong>of</strong> a semigroup <strong>of</strong> contractions<br />

on Y .<br />

Pro<strong>of</strong>. The pro<strong>of</strong> is similar to the one <strong>of</strong> Theorem 6.2.3.<br />

Now, we consider the wave equation with a control in a Neumann boundary condition:<br />

∂2z ∂z<br />

− ∆z = f in Q,<br />

∂t2 ∂n = u on Σ, z(x, 0) = z0 and ∂z<br />

∂t (x, 0) = z1 in Ω. (6.4.11)<br />

For any u ∈ L2 (Γ), the mapping ζ ↦→ <br />

Γ uζ is a continuous linear on H1 (Ω). Thus it can be<br />

identified with an element <strong>of</strong> (H1 (Ω)) ′ . Thus for u ∈ L2 (Σ), the mapping ζ ↦→ <br />

u(·)ζ is an<br />

element <strong>of</strong> L2 (0, T ; (H1 (Ω)) ′ ). Let us denote this mapping by û. We set<br />

<br />

y1<br />

<br />

0<br />

V =<br />

0<br />

û<br />

<br />

, F =<br />

0<br />

f<br />

<br />

, L<br />

Equation (6.4.11) may be written in the form<br />

y2<br />

=<br />

y1<br />

dy<br />

dt = ( A + L)y + F + V, y(0) = y0,<br />

<br />

, and y0 =<br />

z0<br />

with F and V belong to L 2 (0, T ; L 2 (Ω)) × L 2 (0, T ; (H 1 (Ω)) ′ ), y0 ∈ L 2 (Ω) × (H 1 (Ω)) ′ .<br />

Theorem 6.4.4 For every (f, u, z0, z1) ∈ L 2 (Q)×L 2 (Σ)×L 2 (Ω)×(H 1 (Ω)) ′ , equation (6.4.11)<br />

admits a unique weak solution z(f, u, z0, z1) in C([0, T ]; L 2 (Ω))∩C 1 ([0, T ]; (H 1 (Ω)) ′ ). Moreover<br />

the mapping (f, u, z0, z1) ↦→ z(f, u, z0, z1) is continuous from L 2 (Q)×L 2 (Σ)×L 2 (Ω)×(H 1 (Ω)) ′<br />

into C([0, T ]; L 2 (Ω)) ∩ C 1 ([0, T ]; (H 1 (Ω)) ′ ).<br />

Pro<strong>of</strong>. The result is a direct consequence <strong>of</strong> Theorem 6.4.3.<br />

We consider the control problem<br />

(P4) inf{J4(z, u) | (z, u) ∈ C([0, T ]; L 2 (Ω)) × L 2 (0, T ; L 2 (Ω)), (z, u) satisfies (6.4.11)},<br />

with<br />

J4(z, u) = 1<br />

<br />

(z − zd)<br />

2 Q<br />

2 + 1<br />

<br />

(z(T ) − zd(T ))<br />

2 Ω<br />

2 + β<br />

<br />

u<br />

2 Σ<br />

2 .<br />

z1<br />

Γ<br />

<br />

.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!