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

70 CHAPTER 6. CONTROL OF THE WAVE EQUATION Theorem 6.4.5 Assume that f ∈ L 2 (0, T ; L 2 (Ω)), z0 ∈ L 2 (Ω), z1 ∈ (H 1 (Ω)) ′ , and zd ∈ C([0, T ]; L 2 (Ω)). Problem (P2) admits a unique solution (¯z, ū). Moreover the optimal control ū is defined by ū = − 1 β p|Σ, where p is the solution to the equation ∂ 2 p ∂t 2 − ∆p = ¯z − zd in Q, p(T ) = 0, ∂p ∂n = 0 on Σ = Γ×]0, T [, ∂p ∂t (T ) = ¯z(T ) − zd(T ) in Ω. Proof. We leave the reader adapt the proof of Theorem 6.3.1. 6.5 Trace regularity (6.4.12) To study the wave equation with a control in a Dirichlet boundary condition, we have to establish a sharp regularity result stated below. Theorem 6.5.1 Let y be the solution to the equation ∂2y ∂t2 − ∆y = f in Q, y = 0 on Σ, y(x, 0) = y0 and ∂y ∂t (x, 0) = y1 in Ω. (6.5.13) We have ∂y ∂n L2 (Σ) ≤ C fL2 (Q) + y0H 1 0 (Ω) + y1L2 (Ω) . (6.5.14) The proof can be found in [14, Theorem 2.2]. 6.6 Dirichlet boundary control We consider the wave equation with a control in a Dirichlet boundary condition ∂2z ∂t2 − ∆z = f in Q, z = u on Σ, z(x, 0) = z0 and ∂z ∂t (x, 0) = z1 in Ω. (6.6.15) As for the heat equation with a Dirichlet boundary control, the solution to equation (6.6.15) is defined by the transposition method. Definition 6.6.1 A function z ∈ C([0, T ]; L2 (Ω)) ∩ C1 ([0, T ]; H−1 (Ω)) is called a weak solution to equation (6.6.15) if, and only if, fy dxdt = −〈 ∂z ∂t (0), y(0)〉 H −1 (Ω)×H 1 0 (Ω) − Q Q Ω zϕ dxdt + 〈 ∂z ∂t (T ), yT 〉 H −1 (Ω)×H 1 0 (Ω) z(T )νT dx + Ω z(0) ∂y ∂y (0) dx + u dsdt (6.6.16) ∂t Σ ∂n for all (ϕ, yT , νT ) ∈ L 1 (0, T ; L 2 (Ω)) × H 1 0(Ω) × L 2 (Ω), where y is the solution to ∂ 2 y ∂t 2 − ∆y = ϕ in Q, y = 0 on Σ, y(x, T ) = yT and ∂y ∂t (x, T ) = νT in Ω. (6.6.17)

6.6. DIRICHLET BOUNDARY CONTROL 71 Theorem 6.6.1 For every (f, u, z0, z1) ∈ L 1 (0, T ; H −1 (Ω)) × L 2 (Σ) × L 2 (Ω) × H −1 (Ω), equation (6.6.15) admits a unique weak solution z(f, u, z0, z1) in C([0, T ]; L 2 (Ω))∩C 1 ([0, T ]; H −1 (Ω)). The mapping (f, u, y0, y1) ↦→ z(f, u, z0, z1) is linear and continuous from L 1 (0, T ; H −1 (Ω)) × L 2 (Σ) × L 2 (Ω) × H −1 (Ω) into C([0, T ]; L 2 (Ω)) ∩ C 1 ([0, T ]; H −1 (Ω)). Proof. This existence and regularity result can be proved by the transposition method with Theorem 6.5.1. (i) Due to Theorem 6.2.3, the mapping (f, z0, z1) ↦−→ z(f, 0, z0, z1) is linear and continuous from L 1 (0, T ; H −1 (Ω)) × L 2 (Ω) × H −1 (Ω) into C([0, T ]; L 2 (Ω)) ∩ C 1 ([0, T ]; H −1 (Ω)). Thus we have only to consider the case where (f, z0, z1) = (0, 0, 0). (ii) Denote by Λ the mapping ϕ ↦→ ∂y , where y is the solution to equation (6.6.17) cor- ∂n responding to (yT , νT ) = (0, 0). Due to Theorem 6.5.1, Λ is a linear operator from L2 (Q) into L2 (Σ). If we set z = Λ∗u, with u ∈ L2 (Σ), we observe that z ∈ L2 (Q), and z is a solution to equation (6.6.15) with (f, z0, z1) = (0, 0, 0), in the sense of definition 6.6.1. This solution is unique in L2 (Q). Indeed if z1 and z2 are two solutions to equation (6.6.15) with (f, z0, z1) = (0, 0, 0), in the sense of definition 6.6.1, we have (z1 − z2)ϕ = 0 for all ϕ ∈ L Ω 2 (Q). To prove that z belongs to C([0, T ]; L 2 (Ω))∩C 1 ([0, T ]; H −1 (Ω)), we proceed by approximation. Let (un)n be a sequence of regular functions such that zn = z(0, un, 0, 0) be regular. For τ ∈]0, T ] and (yτ, ντ) ∈ H 1 0(Ω) × L 2 (Ω), we denote by y(yτ, ντ) the solution to the equation ∂2y ∂t2 − ∆y = 0 in Q, y = 0 on Σ, y(x, τ) = yτ and ∂y ∂t (x, τ) = ντ in Ω. With Theorem 6.5.1, we have ∂y(yτ, ντ) L2 (Σ) ≤ C yτH 1 0 (Ω) + ντL2 (Ω) , (6.6.18) ∂n where the constant C depends on T , but is independent of τ. equation (6.6.15), according to definition 6.6.1, we have Since zn is the solution to 〈 ∂zn ∂t (τ), yτ〉 H−1 (Ω)×H1 0 (Ω) = ∂zn ∂t (τ)yτ dx = − ∂y(yτ, 0) u dsdt ∂n and From which we deduce: Ω zn(τ)νT dx = Ω zn − zm C([0,T ];L 2 (Ω)) = sup τ∈]0,T ]sup ντ L 2 (Ω) =1 Σ ≤ Cun − um L 2 (Σ) Σ ∂y(0, ντ) un dsdt. ∂n Σ ∂y(0, ντ) ∂n (un − um) dsdt

6.6. DIRICHLET BOUNDARY CONTROL 71<br />

Theorem 6.6.1 For every (f, u, z0, z1) ∈ L 1 (0, T ; H −1 (Ω)) × L 2 (Σ) × L 2 (Ω) × H −1 (Ω), equation<br />

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

The mapping (f, u, y0, y1) ↦→ z(f, u, z0, z1) is linear and continuous from L 1 (0, T ; H −1 (Ω)) ×<br />

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

Pro<strong>of</strong>. This existence and regularity result can be proved by the transposition method with<br />

Theorem 6.5.1.<br />

(i) Due to Theorem 6.2.3, the mapping<br />

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

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

C 1 ([0, T ]; H −1 (Ω)). Thus we have only to consider the case where (f, z0, z1) = (0, 0, 0).<br />

(ii) Denote by Λ the mapping ϕ ↦→ ∂y<br />

, where y is the solution to equation (6.6.17) cor-<br />

∂n<br />

responding to (yT , νT ) = (0, 0). Due to Theorem 6.5.1, Λ is a linear operator from L2 (Q)<br />

into L2 (Σ). If we set z = Λ∗u, with u ∈ L2 (Σ), we observe that z ∈ L2 (Q), and z is a<br />

solution to equation (6.6.15) with (f, z0, z1) = (0, 0, 0), in the sense <strong>of</strong> definition 6.6.1. This<br />

solution is unique in L2 (Q). Indeed if z1 and z2 are two solutions to equation (6.6.15) with<br />

(f, z0, z1) = (0, 0, 0), in the sense <strong>of</strong> definition 6.6.1, we have<br />

<br />

(z1 − z2)ϕ = 0 for all ϕ ∈ L<br />

Ω<br />

2 (Q).<br />

To prove that z belongs to C([0, T ]; L 2 (Ω))∩C 1 ([0, T ]; H −1 (Ω)), we proceed by approximation.<br />

Let (un)n be a sequence <strong>of</strong> regular functions such that zn = z(0, un, 0, 0) be regular. For<br />

τ ∈]0, T ] and (yτ, ντ) ∈ H 1 0(Ω) × L 2 (Ω), we denote by y(yτ, ντ) the solution to the equation<br />

∂2y ∂t2 − ∆y = 0 in Q, y = 0 on Σ, y(x, τ) = yτ and ∂y<br />

∂t (x, τ) = ντ in Ω.<br />

With Theorem 6.5.1, we have<br />

∂y(yτ, ντ)<br />

<br />

<br />

L2 (Σ) ≤ C yτH 1<br />

0 (Ω) + ντL2 (Ω) , (6.6.18)<br />

∂n<br />

where the constant C depends on T , but is independent <strong>of</strong> τ.<br />

equation (6.6.15), according to definition 6.6.1, we have<br />

Since zn is the solution to<br />

〈 ∂zn<br />

∂t (τ), yτ〉<br />

<br />

H−1 (Ω)×H1 0 (Ω) =<br />

∂zn<br />

∂t (τ)yτ<br />

<br />

dx = −<br />

∂y(yτ, 0)<br />

u dsdt<br />

∂n<br />

and <br />

From which we deduce:<br />

Ω<br />

zn(τ)νT dx =<br />

Ω<br />

zn − zm C([0,T ];L 2 (Ω)) = sup τ∈]0,T ]sup ντ L 2 (Ω) =1<br />

<br />

Σ<br />

≤ Cun − um L 2 (Σ)<br />

Σ<br />

∂y(0, ντ)<br />

un dsdt.<br />

∂n<br />

<br />

<br />

<br />

<br />

Σ<br />

∂y(0, ντ)<br />

∂n<br />

<br />

<br />

(un − um) dsdt

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

Saved successfully!

Ooh no, something went wrong!