UNIVERSITÉ ALIOUNE DIOP DE BAMBEY

Département : Mathématiques
Filière : Mathématiques et Application
Option : Analyse Non-linéaire

Projet Contrôle des EDO

Fallou BEYE
Analyse Numérique
Année universitaire 2021/2022
Projet Contrôle des équations aux dérivées partielles

Soit Ω un ouvert non vide, borné de RN de classe C2 .

On considéré le système suivant :


−∆y = f dans Ω
(∗)
y = 0 sur ∂Ω
avec f ∈ L2 (Ω).
On prend comme ensemble des contrôles admissibles Uad = L2 (Ω).

Soit le système suivant :


−∆y(u) = f + u dans Ω
(∗∗) ,
y(u) = 0 sur ∂Ω

où yd ∈ L2 (∂Ω) est l’état de référence.

On pose
Z  2 Z
1 ∂y(u) α
J(u) = − yd dσ + u2 dx.
2 ∂Ω ∂n 2 Ω

L’objectif est de trouver ū ∈ Uad qui minimise le fonctionnel J.

Existence et unicité :
On pose l’application L définie par :

L : L2 (Ω) −→ L2 (∂Ω)
∂y
f 7−→ L(f) =
∂n
où y est solution de l’équation (∗). La fonction L ainsi définie est linéaire et continue.
En effet :
— Linéarité :
Soit λ ∈ R.
On considère y, z ∈ H10 (Ω) telles que y soit solution de

−∆y = f dans Ω
(Ey ) :
y = 0 sur ∂Ω

et z la solution de 
−∆z = g dans Ω
(Ez ) : ,
z = 0 sur ∂Ω

1
 avec f, g ∈ L2 (Ω).
On peut donc dire que λy + z est solution de l’équation

−λ∆y − ∆z = λf + g dans Ω
(E) :
λy + z = 0 sur ∂Ω

−∆(λy + z) = λf + g dans Ω
≡ , (par linéarité de lópérateur ∆).
λy + z = 0 sur ∂Ω

Ainsi on a
∂(λy + z)
L(λf + g) =
∂n
∂y ∂z
= λ + (par linéarité de lópérateur ∂)
∂n ∂n
= λL(f) + L(g)

D’où la linéarité de L.

— Continuité :
soit ϕ ∈ H10 (Ω) Z Z
− ∆yϕdx = fϕdx
Ω Ω
En appliquant la formule de Green, on obtient :

Z Z Z
∂y
∇y∇ϕdx − ϕdσ = fϕdx
∂n
ZΩ ∂Ω Ω
∂y
Or ϕdσ = 0 car ϕ|∂(Ω) = 0
∂n
Z ∂Ω Z
=⇒ ∇y∇ϕdx = fϕdx
Ω Ω

En utilisant l’inégalité de Cauchy-Schwartz, on aura

Z
fϕdx 6 kfkL2 (Ω) kϕkL2 (Ω) .
Ω

Appliquons l’inégalité de Poincaré sur ϕ ∈ H10 (Ω), on obtient

kϕkL2 (Ω) 6 C1 k∇ϕkL2 (Ω) ,

où C1 est une constante ne dépendant que de Ω.

Donc Z Z
∇y∇ϕdx = fϕdx 6 C1 kfkL2 (Ω) k∇ϕkL2 (Ω) .
Ω Ω

Puisque cela est vrai pour toute fonction test ϕ dans H10 (Ω), en choisissant ϕ = y,
on a : Z
2
k∇ykL2 (Ω) = fydx 6 C1 kfkL2 (Ω) k∇ykL2 (Ω) .
Ω

=⇒ k∇ykL2 (Ω) 6 C1 kfkL2 (Ω) .

2
 Maintenant, appliquons à ∇y l’opérateur de trace

γ : H1 (Ω) −→ L2 (∂Ω).
∂y
Donc, γ(∇y) = ∂n .
Ainsi, on peut établir l’égalité
∂y
k kL2 (∂Ω) = kγ(∇y)kL2 (∂Ω) = k∇ykL2 (Ω) .
∂n
Ceci nous permet de dire qu’il existe une constante C2 telle que :
∂y
k kL2 (∂Ω) 6 C2 k∇ykL2 (Ω) .
∂n
Donc
1 ∂y
k kL2 (∂Ω) 6 k∇ykL2 (Ω) 6 C1 kfkL2 (Ω) ,
C2 ∂n
∂y
=⇒ k kL2 (∂Ω) 6 C1 .C2 kfkL2 (Ω) .
∂n
=⇒ kL(f)kL2 (∂Ω) 6 CΩ kfkL2 (Ω) ,
avec CΩ = C1 .C2 .
D’où la continuité de L

Pour tout u ∈ Uad , on a

∂y(u)
= L(f + u) = L(u) + L(f)
∂n
∂y(0)
=⇒ = L(f + 0) = L(f)
∂n
∂y(u) ∂y(0)
=⇒ = + L(u)
∂n ∂n
∂y(u) ∂y(0)
On a donc L(u) = −
∂n ∂n

Z  2 Z
1 ∂y(u) α
J(u) = − yd u2 dx
dσ +
2 ∂Ω ∂n 2 Ω
Z  2 Z
1 ∂y(u) ∂y(0) ∂y(0) α
= − + − yd dσ + u2 dx
2 ∂Ω ∂n ∂n ∂n 2 Ω
Z  2 Z
1 ∂y(0) α
= L(u) + − yd dσ + u2 dx
2 ∂Ω ∂n 2 Ω
Z Z   Z  2 Z
1 2 ∂y(0) 1 ∂y(0) α
= [L(u)] dσ + L(u) − yd dσ + − yd dσ + u2 dx
2 ∂Ω ∂Ω ∂n 2 ∂Ω ∂n 2 Ω
Z Z Z   Z  2
1 2 α 2 ∂y(0) 1 ∂y(0)
J(u) = [L(u)] dσ + u dx + L(u) − yd dσ + − yd dσ
2 ∂Ω 2 Ω ∂Ω ∂n 2 ∂Ω ∂n

On peut donc écrire

1
J(u) = a(u, u) + l(u) + C,
2
3
 Z Z Z  
2 2 ∂y(0)
avec a(u, u) = [L(u)] dσ + α u dx, l(u) = L(u) − yd dσ et
∂n
Z  ∂Ω
2 Ω ∂Ω
1 ∂y(0)
C= − yd dσ
2 ∂Ω ∂n

∀u, v ∈ L2 (Ω), on a
Z Z
a(u, v) = L(u)L(v)dσ + α uvdx
Z∂Ω ZΩ
= L(v)L(u)dσ + α vudx
∂Ω Ω
= a(v, u)
Donc a(., .) est symétrique.

∀u ∈ L2 (Ω), on a si de plus u 6= 0,
Z Z Z
2 2
a(u, u) = [L(u)] dσ + α u dx > α u2 dx > 0,
∂Ω Ω Ω

=⇒ a(u, u) > 0.
a(., .) est donc définie positive.
Comme a(., .) est symétrique, alors pour conclure de sa bilinéarité, il suffit d’étudier la
linéarité par rapport à une de ses deux variables.
Soit λ ∈ R et soient u, v et w des éléments de L2 (Ω).
Z Z
a(λu + v, w) = L(λu + v)L(w)dσ + α (λu + v)wdx
Z∂Ω Ω Z

= [λL(u) + L(v)]L(w)dσ + α (λu + v)wdx

Z
∂Ω Z Ω Z Z
= λ L(u)L(w)dσ + λα uwdx + L(v)L(w)dσ + α vwdx
∂Ω Ω ∂Ω Ω
= λa(u, w) + a(v, w).
Donc a(., .) est bilinéaire.

Pour démontrer la continuité de a(., .), on choisit u, v ∈ L2 (Ω) et on montre qu’il existe
une constante M > 0 telle que
|a(u, v)| 6 MkukL2 (Ω) kukL2 (Ω)
. On a
Z Z
|a(u, v)| = L(u)L(v)dσ + α uvdx
Z ∂Ω Ω
Z
6 L(u)L(v)dσ + α uvdx
Z ∂Ω
Z Ω
6 |L(u)L(v)|dσ + α |uv|dx
∂Ω Ω
6 kL(u)kL2 (∂Ω) kL(v)kL2 (∂Ω) + αkukL2 (Ω) kvkL2 (Ω)
6 CkukL2 (Ω) C 0 kvkL2 (Ω) + αkukL2 (Ω) kvkL2 (Ω) (par continuité de L),
=⇒ |a(u, v)| 6 (C.C 0 + α)kukL2 (Ω) kvkL2 (Ω)

4
En posant M = C.C 0 + α > 0, on obtient la continuité de a(., .).
On peut donc conclure que a(., .) est strictement convexe et continue.
R  
∀u ∈ L2 (Ω), l(u) = ∂Ω L(u) ∂y(0) ∂n
− y d dσ est linéaire (grâce à la linéarité de l’inté-
grale).
Pour la continuité, soit u ∈ L2 (Ω)
Z  
∂y(0)
|l(u)| = L(u) − yd dσ
∂Ω ∂n
∂y(0)
6 kL(u)kL2 (∂Ω) k − yd kL2 (∂Ω)
∂n
∂y(0)
6 CkukL2 (Ω) k − yd kL2 (∂Ω) .
∂n

Alors on obtient la continuité de l(.) en posant M 0 = Ck ∂y(0)

∂n
− yd kL2 (∂Ω) > 0.

On peut donc conclure que J strictement convexe et continue et donc elle est semi-
continue inférieurement (s.c.i).
De plus J est coercive.
En effet,
α
∀u ∈ L2 (Ω), J(u) > kuk2L2 (Ω) ,
2
=⇒ lim J(u) = +∞
kuk→+∞

Alors
∃! ū ∈ L2 (Ω), J(ū) = min J(u).
u∈L2 (Ω)

Système d’optimalité
J est différentiable et pour tout v ∈ L2 (Ω), on a

dJ(ū).(v) = 0

J(ū + tv) − J(ū)

dJ(ū).(v) = lim
t→0 t
1
[ a(ū + tv, ū + tv) + l(ū + tv) + C] − J(ū)
= lim 2
t→0 t
2
J(ū) + 2 a(ū, v) + 2 a(v, ū) + t2 a(v, v) + tl(v) − J(ū)
t t
= lim
t→0 t
1 1 t
= lim a(ū, v) + a(v, ū) + a(v, v) + l(v)
t→0 2 2 2
1
= [a(ū, v) + a(v, ū)] + l(v)
2
1
= .2a(ū, v) + l(v) (par symétrie de a(., .)),
2
=⇒ ∀v ∈ L2 (Ω), a(ū, v) + l(v) = 0 (Equation d’Euler-Lagrange).

5
 Z Z Z  
∂y(0)
=⇒ L(ū)L(v)dσ + α ūvdx + L(v) − yd dσ = 0
∂n
Z 
∂Ω Ω
 ∂Ω
Z Z  
∂y(ū) ∂y(0) ∂y(0)
=⇒ − L(v)dσ + α ūvdx + L(v) − yd dσ = 0
∂Ω ∂n ∂n Ω ∂Ω ∂n
∂z
L(v) ∈ L2 (∂Ω) est tel qu’il existe z ∈ H10 (Ω) qui vérifie L(v) = ∂n
et qui est solution de

−∆z = v dans Ω
z = 0 sur ∂Ω

Soit p̄ ∈ H1 (Ω) l’unique solution de


−∆p̄ = 0 dans Ω 
∂y(ū) ∂y(0)
p̄ = ∂n
− ∂n
sur ∂Ω

On peut donc écrire

Z Z
∂z
p̄ dσ + α ūvdx = 0
∂Ω ∂n Ω

En appliquant la formule de Green, on obtient

Z Z Z
∇p̄∇zdx + p̄∆zdx + α ūvdx = 0
ZΩ Z Ω Z Ω Z
∂p̄
=⇒ z dσ − z∆p̄dx + p̄∆zdx + α ūvdx = 0
∂Ω ∂n Ω Ω Ω

Or z|∂Ω = 0, −∆p̄ = 0 dans Ω et aussi −∆z = v dans Ω. Alors on a

Z Z
− p̄vdx + α ūvdx = 0
ZΩ Ω

=⇒ (αū − p̄)vdx = 0
Ω
=⇒ αū − p̄ = 0
1
=⇒ ū = p̄, avec α > 0.
α
Comme p̄ ∈ H1 (Ω) alors ū ∈ H1 (Ω).
D’où le système optimal suivant :


 1 1 1
ū ∈ H (Ω), y(ū) ∈ H0 (Ω), p̄ ∈ H (Ω)



 −∆y(ū) = f + ū dans Ω



y(ū) = 0 sur ∂(Ω)
(SO) :

 −∆p̄= 0 dans Ω



p̄ = ∂y( ū)
− yd sur ∂(Ω)

 ∂n

ū = 1 p̄ ∈ H1 (Ω)
α