29.03.2013 Views

Optimal Control of Partial Differential Equations

Optimal Control of Partial Differential Equations

Optimal Control of Partial Differential Equations

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

8.4. THE HEAT EQUATION 89<br />

For every z0 ∈ Z, f ∈ L 2 (Q), and u ∈ L 2 (0, T ; L 2 (Γ)) equation (8.4.16) admits a unique<br />

solution in H1 (0, T ; D(A∗ ) ′ ) which is<br />

t<br />

z(t) = e t(A∗ ) ∗<br />

z0 +<br />

0<br />

(I − (A ∗ ) ∗ )e (t−s)(A∗ ) ∗<br />

Nu(s) ds +<br />

t<br />

0<br />

e (t−s)A f(s) ds. (8.4.17)<br />

When u ∈ C 1 ([0, T ]; H 1<br />

2 (Γ)), the solutions given by Theorem 5.3.4 and by the formula (8.4.17)<br />

coincide. Henceforth, by density arguments it follows that the solution defined by (8.4.17) and<br />

the one <strong>of</strong> Theorem 5.3.4 also coincide when u ∈ L 2 (0, T ; L 2 (Γ)). In this case to simplify the<br />

writing, we <strong>of</strong>ten write<br />

z ′ = Az + f + (I − A)Nu, z(0) = z0, (8.4.18)<br />

in place <strong>of</strong> (8.4.16). Since N ∈ L(L2 (Γ); D((I −A) α )) for all α ∈]0, 3[,<br />

the operator (I −A)N =<br />

4<br />

(I − A) 1−α (I − A) αN can be decomposed in the form (I − A)N = (I − A) 1−αB1, where<br />

B1 = (I − A) αN belongs to L(L2 (Γ); L2 (Ω)). This decomposition will be very useful to study<br />

the Riccati equation corresponding to problem (P3) (see [26]).<br />

8.4.2 Dirichlet boundary control<br />

We set Z = L 2 (Ω) and we define the unbounded operator A in Z by<br />

D(A) = H 2 (Ω) ∩ H 1 0(Ω), Az = ∆z.<br />

We know that A is the infinitesimal generator <strong>of</strong> an analytic semigroup on Z and that (see<br />

[2])<br />

D((−A) α ) = H 2α (Ω) if α ∈]0, 1<br />

4 [<br />

and<br />

D((−A) α ) = {z ∈ H 2α (Ω) | z = 0 on Γ} if α ∈] 1<br />

, 1[.<br />

4<br />

We define the Dirichlet operator G from L2 (Γ) into L2 (Ω) by G :<br />

solution to<br />

u ↦→ w , where w is the<br />

∆w = 0 in Ω, w = u on Γ.<br />

From [12] we deduce that G ∈ L(L 2 (Γ); H 1<br />

2 (Ω)). This implies that G ∈ L(L 2 (Γ); D((−A) α ))<br />

for all α ∈]0, 1<br />

4<br />

[. We also have G ∈ L(H 3<br />

2 (Γ); H 2 (Ω)).<br />

Suppose that u ∈ C 1 ([0, T ]; H 3<br />

2 (Γ)) and denote by z the solution to equation (5.4.20). Set<br />

y(x, t) = z(x, t)−(Gu(t))(x). Since u ∈ C 1 ([0, T ]; H 3<br />

2 (Γ)), Gu(·) belongs to C 1 ([0, T ]; H 2 (Ω)).<br />

Thus we have<br />

∂y<br />

∂t<br />

− ∆y = f − ∂Gu<br />

∂t<br />

in Q, y = 0 on Σ, y(x, 0) = (z0 − Gu(0))(x) in Ω.<br />

As for Neumann controls we can check that, when u is regular enough, equation (5.4.20) may<br />

be written in the form<br />

z ′ = Az + f + (−A)Gu, z(0) = z0.<br />

And in the case when u ∈ L 2 (0, T ; L 2 (Γ)) we still continue to use the above formulation even<br />

if for a correct writing A should be replaced by its extension (A ∗ ) ∗ .

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

Saved successfully!

Ooh no, something went wrong!