1

Notes de ours de l'ISIMA, troisième année

http://www.isima.fr/∼leborgne

Comprendre la ondition inf-sup et le lo king.

Exemples : Stokes, lapla ien, bilapla ien, lo king de type
Kir hhoLove ou MindlinReissner, onditions aux limites
Gilles Leborgne
21 juin 2018

La ondition inf-sup, relative aux équations aux dérivées partielles elliptiques sous ontraintes
anes, est une manière d'exprimer le théorème de l'appli ation ouverte, voir Brézis [6℄.
La le ture peut ommen er au Ÿ 7 qui traite du problème modèle de Stokes. Et e problème
modèle justie le besoin des résultats mathématiques donnés dans les premiers Ÿ.

Table des matières

0 Introdu ution : 'est quoi la ondition inf-sup ? 3

I Théorie 4

1 Notations 4
2 Théorème de l'appli ation ouverte et ondition inf-sup 4
3 Problème mixte (problème sous ontrainte lassique) 6
3.1 Le problème : forme faible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2 Le problème : forme forte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.3 Le problème matri iel asso ié . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.4 Le théorème inf-sup . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.5 Le lagrangien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

4 Opérateurs divergen e, gradient, rotationnel 8

5 Cadre fon tionnel 9
5.1 Espa es de Hilbert utilisés dans la suite . . . . . . . . . . . . . . . . . . . . . . . . 9
5.2 Sous-espa es de Hilbert . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1
5.3 Appli ation tra e γ0 et l'Hilbert H 2 (Γ) . . . . . . . . . . . . . . . . . . . . . . . . 10
5.4 Autres appli ations tra es ontinues . . . . . . . . . . . . . . . . . . . . . . . . . . 11
5.5 Dual et Bana h . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
5.6 Intégrations par parties (IPP) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1 div
5.7 Dual de H (Ω) et H (Ω) ( ara térisations) . . . . . . . . . . . . . . . . . . . . . 12
5.8 Noyaux des appli ations tra es . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5.9 Poin aréFriedri hs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2 n
5.10 Dé omposition de L (Ω) (Helmholtz simple) . . . . . . . . . . . . . . . . . . . . . 14

6 Les surje tivités de l'opérateur divergen e 14

6.1 Surje tivité de l'opérateur divergen e div : H div (Ω) → L2 (Ω) . . . . . . . . . . . . 15
6.2 Surje tivité de l'opérateur divergen e H0div (Ω) → L20 (Ω) . . . . . . . . . . . . . . . 15
n
6.3 Surje tivité de l'opérateur divergen e L2 (Ω) → H −1 (Ω) . . . . . . . . . . . . . . . 15
n
6.4 Surje tivité de l'opérateur divergen e div : H01 (Ω) → L20 (Ω) . . . . . . . . . . . . . 15

II Appli ations 17

7 Le modèle Stokes 17
7.1 Problème sous ontrainte initial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
7.2 Problème ave ontrainte relaxée . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

1 21 juin 2018
2

8 Approximation numérique du problème de Stokes 17

8.1 Approximation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
8.2 Proje tions (méthode des éléments nis) . . . . . . . . . . . . . . . . . . . . . . . . 18
8.3 Représentation matri ielle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
8.4 Le problème du multipli ateur de Lagrange (la pression p) . . . . . . . . . . . . . . 19
8.5 Analyse à l'aide des proje tions : des pertes . . . . . . . . . . . . . . . . . . . . . . 19
8.6 Une orre tion : ompensation des pertes . . . . . . . . . . . . . . . . . . . . . . . 20
8.7 Méthode de Brezzi et Pitkäranta . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
8.8 Méthode de Hughes, Fran a et Balestra . . . . . . . . . . . . . . . . . . . . . . . . 21
8.9 Méthode de Douglas et Wang . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

9 Lapla ien 22
9.1 Problème initial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
9.2 Problème sous ontrainte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

10 Bilapla ien (problème biharmonique) 23

10.1 Problème initial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
10.2 Formulation ave la variable ∆p . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
10.3 Formulation ave la variable ~
gradp . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
10.3.1 Formulation initiale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
10.3.2 Première formulation sous forme ontrainte . . . . . . . . . . . . . . . . . . 24
10.3.3 Formulation bis sous forme ontrainte . . . . . . . . . . . . . . . . . . . . . 24
10.3.4 Formulation ter sous forme ontrainte . . . . . . . . . . . . . . . . . . . . . 25
10.3.5 Formulation quater sous forme ontrainte . . . . . . . . . . . . . . . . . . . 26

11 Lo king 27
11.1 Le problème type . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
11.2 Coer ivité nie en p . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
11.3 Le problème dis rétisé : onstante de oer ivité en ph non ontrlée . . . . . . . . . 28
11.4 Une orre tion optimale du lo king . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
11.5 Traitement lassique du lo king : problème ontinu . . . . . . . . . . . . . . . . . . 30
11.5.1 Problème ontinu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
11.5.2 Problème dis rétisé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

12 Conditions aux limites de Diri hlet faible 32

12.1 Problème initial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
12.2 Problème ontraint . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
12.3 Problème appro hé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
12.4 Eléments nisPk −C 0 : non stables . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
12.4.1 Condition inf-sup dis rète . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
12.4.2 Barbosa et Hughes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
12.4.3 Elimination du multipli ateur : méthode de Nits he [25℄ . . . . . . . . . . . 34

A Singular value de omposition (SVD) 36

B Appli ation à la ondition inf-sup 38

2 21 juin 2018
3 0. Introdu ution : 'est quoi la ondition inf-sup ?

0 Introdu ution : 'est quoi la ondition inf-sup ?

C'est une ondition qui permet d'avoir le ara tère bien posé, i.e. d'avoir une solution qui
dépend ontinûment des données, et, dans les dis rétisations, d'avoir la matri e à inverser qui est
inversible et bien onditionnée.
C'est un orollaire du théorème de l'appli ation ouverte : si T : E → G est une appli ation
−1
linéaire ontinue bije tive entre deux Bana h (espa es ve toriels normés omplets), alors T est
−1
également ontinue. Don si g ∈ G si u est solution de T u = g , alors on aura u = T g ave
||u||E ≤ C||g||G où C = ||T −1 || < ∞. Voir Ÿ 2 pour la démonstration.
Et la ondition inf-sup peut être vue omme une généralisation de la oer ivité : on rappelle
que si a(·, ·) : V × V → R est bilinéaire ontinue sur V Hilbert, elle est dite oer ive ssi :

∃α > 0, ∀v ∈ V, a(v, v) ≥ α||v||2V . (0.1)

Et alors le problème : pour ℓ ∈ V ′, trouver u∈V t.q :

∀v ∈ V, a(u, v) = hℓ, viV ′ ,V , (0.2)

est bien posé. I.e., notant A : V → V ′ l'appli ation dénie par hAu, viV ′ ,V = a(u, v) pour tout
′
u, v ∈ V , le problème : pour ℓ ∈ V , trouver u ∈ V t.q :

Au = ℓ, (0.3)

vérie u = A−1 ℓ ave :

1
∃α > 0, ∀ℓ ∈ V ′ , ||u||V ≤ ||ℓ||V ′ , (0.4)
α
1
i.e. ||A−1 || ≤α < ∞, et A
−1
est ontinue.
Généralisation de la ondition de oer ivité : soit E, F deux Bana h, et F ′ le dual de F (qui
′
est don un Bana h pour la norme usuelle), soit T : E → F ontinue. On onsidère le problème
′
(faire le parallèle ave (0.3)) : pour ℓ ∈ F , trouver u ∈ E t.q :

T u = ℓ. (0.5)

C'est le problème (faire le parallèle ave (0.2) en posant t(u, v) = hT u, viF ′ ,F ) : pour ℓ ∈ F ′ , trouver
u∈E t.q :
∀v ∈ F, hT u, viF ′ ,F = hℓ, viF ′ ,F . (0.6)

−1
Quand T est bije tive e problème admet une unique solution u=T ℓ, et est bien posé quand
(faire le parallèle ave (0.1)) :

∃k > 0, ∀u ∈ E, ∃v ∈ F, |hT u, viF ′ ,F | ≥ k||u||E ||v||F , (0.7)

soit :
hT u, viF ′ ,F
∃k > 0, ∀u ∈ E, ∃v ∈ F, sup ≥ k||u||E , (0.8)
v∈F |v||F
e qui indique que T −1 est ontinu (théorème de l'appli ation ouverte, voir Ÿ 2), et qui s'exprime
sous la forme  ondition inf-sup :

hT u, viF ′ ,F
∃k > 0, inf (sup ) ≥ k, (0.9)
u∈E v∈F ||u||E ||v||F

(On a allégé l'é riture en omettant de pré iser u 6= 0 et v 6= 0).

Et la ondition inf-sup exprime également la ondition  T est à image fermée. En eet,
on onsidère la restri tion à
( )E/KerT de l'appli ation linéaire T : E → F, sous la forme
E/KerT → ImT
Tr : , appli ation qui est linéaire ontinue (ave E/KerT muni de la
u → Tr u = T u
norme quotient et ImT muni de la norme ||.||F )
et bije tive (par dénition du noyau, de l'espa e
′
quotient, et de l'image), et l'hypothèse ImT est
fermée dans F indique que (ImT, ||.||F ′ ) est un
−1
Bana h. Don le théorème de l'appli ation ouverte s'applique à Tr , et don Tr : ImT → E/KerT
hT u, viF ′ ,F
est ontinue. D'où la ondition inf-sup : ∃k > 0, inf ( sup ) ≥ k.
u∈E/KerT v∈ImT ||u||E/KerT ||v||F

3 21 juin 2018
4 2. Théorème de l'appli ation ouverte et ondition inf-sup

Première partie

Théorie
1 Notations

Soit (E, ||.||E ) un espa e de Bana h (espa e ve toriel normé omplet). On note BE (x, ρ) la
boule de entre x et rayon ρ.
′
On note E = L(E; R) le dual topologique de E (l'ensemble des formes linéaires ontinues sur
(E, ||.||E )). Pour ℓ ∈ E ′ et x ∈ E on note ℓ(x) = ℓ.x = hℓ, xiE ′ ,E ( ro het de dualité). La norme
′ ′
(usuelle) sur E est donnée par, pour ℓ ∈ E :

|hℓ, xiE ′ ,E |
||ℓ||E ′ = sup = sup |hℓ, xiE ′ ,E |. (1.1)
x∈E ||x||E ||x||E =1

Ainsi (E ′ , ||.||E ′ )
est un Bana h.
Le bidual de
′′
est E =
déf (E ′ )′ = L(E ′ , R). Ainsi (E ′′ , ||.|| ′′ ) est un Bana h. Et on dispose
E E
′′
de l'inje tion anonique J : E → E (linéaire et isométrique) donnée par :

∀x ∈ E, ∀ℓ ∈ E ′ , J(x).ℓ = ℓ.x. (1.2)

Lorsque ette inje tion est une bije tion, on dit que E est réexif, et on identie E ′′ à E à l'aide
de J, et on note J(x) = x.
(En parti ulier les Hilbert sont réexifs, f. théorème de représentation
′′
de Riesz.) (Interprétation du bidual E : si x est un ve teur dans E alors J(x) est un opérateur
de dérivation dans la dire tion x.)
Pour E et F deux espa es ve toriels et T :E→F linéaire, on note T (x) = T.x (notation de la
distributivité).
Pour (E, ||.||E ) et (F, ||.||F ) deux espa es de Bana h, on note L(E; F ) l'ensemble des appli ations
linéaires ontinues, et pour T ∈ L(E; F ), la norme (usuelle) de T est le réel :

||T.x||F
||T ||L(E;F ) = sup = sup ||T.x||F . (1.3)
x∈E ||x||E ||x||E =1

Ainsi (L(E; F ), ||.||L(E;F ) ) est un Bana h.

Soit(V, (·, ·)V ) et (Q, (·, ·)Q ) deux espa es de Hilbert.

Soita(·, ·) : V × V → R et b(·, ·) : V × Q → R deux formes bilinéaires ontinues, i.e. il existe
des onstantes ca et cb t.q. :

∀(v, w) ∈ V × V, |a(v, w)| ≤ ca ||v||V ||w||V ,

(1.4)
∀(v, q) ∈ V × Q, |b(v, q)| ≤ cb ||v||V ||q||Q .
|a(v,w)|
On note ||a|| = supv,w∈V ||v||V ||w||V (la plus petite onstante ca ).

2 Théorème de l'appli ation ouverte et ondition inf-sup

Soit (E, ||.||E ) et (F, ||.||F ) deux espa es de Bana h.

Théorème 2.1 (de l'appli ation ouverte.) Si T ∈ L(E; F ) (linéaire et ontinue) est surje tive,
alors :
∃γ > 0 t.q. T (BE (0, 1)) ⊃ BF (0, γ). (2.1)

Et un ouvert de E est transformé en un ouvert de F (pour E et F des Bana h et T surje tif ).

Preuve.
S Voir Brézis [6℄. Etapes : 1- appli ation du lemme de Baire à l'union dénombrable de
fermés n∈N∗ T (BE (0, n)) = F T ), qui donne l'existen e d'un fermé T (BE (0, n))
(surje tivité de

qui ontient un ouvert ; d'où par linéarité T (BE (0, 1)) ontient un ouvert, d'où par linéarité il
existe γ > 0 t.q. T (BE (0, 1)) ⊃ BF (0, 2γ). 2- Puis on utilise la ontinuité qui donne T (BE (0, 1)) ⊃
BF (0, γ) (passage à la limite dans des espa es omplets). Voir Brézis [6℄.

4 21 juin 2018
5 2. Théorème de l'appli ation ouverte et ondition inf-sup

Corollaire 2.2 (Continuité de l'appli ation inverse.) Si T est linéaire, ontinue et bije tive entre
−1
deux Bana h, alors T est linéaire et ontinue ; i.e. si T ∈ L(E; F ) est bije tive (entre deux
−1
Bana h), alors T ∈ L(F ; E). Ave (utilisation de (2.1)) :
1 1
∀y ∈ F, ||T −1 .y||E ≤ ||y||F , don ||T −1 ||L(F ;E) ≤ . (2.2)
γ γ
Don :
∀x ∈ E, ||T.x||F ≥ γ||x||E , (2.3)

et ||T.x||F est minorée sur la sphère unité de E.

Preuve. BF (0, γ) ⊂ T (BE (0, 1)), f. (2.1), donne T −1 (BF (0, γ)) ⊂ BE (0, 1).
Et T
−1
est linéaire ar T est linéaire, don T
−1
(BF (0, 1)) ⊂ BE (0, γ1 ).
Don y ∈ BF (0, 1) donne ||T
−1
.y||E ≤ γ1 ||y||F , i.e. (2.2).
Puis y = T.x donne (2.3) (bije tivité).

Remarque 2.3 Don une appli ation linéaire ontinue bije tive entre deux Bana h
 se omporte

−2 0
omme en dimension nie, et par exemple T : R2 → R2 donnée par la matri e où
 1  0 3
−2 0
||T || = 3 et γ = 2 > 0, et T −1 a pour matri e 1 où ||T
−1
|| = γ1 = 21 > 0.
0 3

Corollaire 2.4 (Condition inf-sup.) Si T ∈ L(E; F ) (linéaire et ontinue) est bije tive alors :

∀x ∈ E, ∃ℓ ∈ F ′ : |hℓ, T.xiF ′ ,F | ≥ γ||ℓ||F ′ ||x||E . (2.4)

Soit :
|hℓ, T.xiF ′ ,F |
∀x ∈ E, sup ≥ γ||x||E . (2.5)
ℓ∈F ′ ||ℓ||F ′
Soit :
|hℓ, T.xiF ′ ,F |
inf ( sup ) ≥ γ. (2.6)
x∈E ℓ∈F ′ ||x||E ||ℓ||F ′
Et, quand E est réexif :
|hℓ, T.xiF ′ ,F |
∃γ ′ > 0, inf ′ (sup ) ≥ γ ′. (2.7)
ℓ∈F x∈E ||x||E ||ℓ||F ′

( )
Ry → R
Preuve. Soit x ∈ E , soit y = T.x ∈ F et soit ℓy : (forme linéaire dénie
ty → ℓy (ty) = t||y||2F
|ℓy (ty)|
sur la droite t → ty ) : on a ||ℓy || = supt∈R
||ty||F = ||y||F . À l'aide du théorème de prolongement
des formes linéaires de HahnBana h, voir Brézis [6℄, on prolonge ℓy en une forme linéaire ℓ dénie
sur tout F où ||ℓ||F ′ = ||ℓy ||F ′ = ||y||F = ||T.x||F et hℓ, T.xiF ′ ,F = hℓ, yiF ′ ,F = hℓy , 1.yiF ′ ,F =
||y||2F = ||T.x|2F = hℓy , T.xiF ′ ,F = ||ℓy ||F ′ ||T x||F = ||ℓ||F ′ ||T x||F , d'où (2.3) donne (2.4). D'où (2.5).
D'où (2.6).
Puis on onsidère le dual de T, i.e. l'appli ation linéaire :
(
′ F ′ → E′
T : (2.8)
ℓ → T ′ .ℓ, où hT ′ .ℓ, xiE ′ ,E = hℓ, T.xiF ′ ,F , ∀x ∈ E.

Et T étant ontinue, T ′ est ontinue, ar |(T ′ .ℓ).x| = |ℓ.(T.x)| ≤ ||ℓ||F ′ ||T ||L(E,F ) ||x||E donne
||T .ℓ||E ′ ≤ ||ℓ||F ′ ||T ||L(E,F ) donne ||T ′ ||L(F ′ ,E ′ ) ≤ ||T ||L(E,F ) . Don T ′ ∈ L(F ′ , E ′ ).
′
′ ′
Et T bije tif donne T bije tif : T .ℓ = 0 donne ℓ.(T.x) = 0 = ℓ.y pour tout x ∈ E don pour
′ ′ ′
tout y ∈ F ( ar T surje tif ), don ℓ = 0, don T inje tif ; et pour k ∈ E on trouve ℓ ∈ F
′ −1 −1
t.q. T .ℓ = k ssi ℓ.T.x = k.x pour tout x ∈ E , i.e. ℓ.y = k.T .y pour tout y , i.e. ℓ = k.T
−1 ′ −1
(véri ation : ℓ = k.T donne (T .ℓ).x = ℓ.(T.x) = k.T .T.x = K.x pour tout x), don T ′ est
′ ′ ′
surje tif. Et (E , ||.||E ′ ) et (F , ||.||F ′ ) sont des Bana h. Et (2.6) appliqué à T ave E réexif (don
′′
E ≃ E ) donne (2.7).

5 21 juin 2018
6 3. Problème mixte (problème sous ontrainte lassique)

Corollaire 2.5 (Condition inf-sup.) Si B ∈ L(E; F ′ ) (linéaire et ontinue) est bije tive et si F est
réexif alors :
∀x ∈ E, ∃y ∈ F : |hB.x, yiF ′ ,F | ≥ γ||y||F ||x||E . (2.9)

Soit :
|hB.x, yiF ′ ,F |
∀x ∈ E, sup ≥ γ||x||E . (2.10)
y∈F ||y||F
Soit :
|hB.x, yiF ′ ,F |
inf (sup ) ≥ γ. (2.11)
x∈E y∈F ||x||E ||y||F
Et, quand E est réexif :
|hB.x, yiF ′ ,F |
∃γ ′ > 0, inf (sup ) ≥ γ ′. (2.12)
y∈F x∈E ||x||E ||y||F

Preuve. On applique le théorème pré édent ave F′ au lieu de F.

Le noyau de T ∈ L(E; F ) est KerT = {x ∈ E : T x = 0}, et l'espa e quotient est E/KerT
(l'espa e des lasses d'équivalen e ẋ = {x + k, ∀k ∈ KerT }). On munit E/KerT de sa norme
usuelle (la norme quotient) : pour ẋ ∈ E/KerT , notant x un représentant de ẋ :

||ẋ||E/KerT = inf ||x + k||E . (2.13)

k∈R

Et, quand E est un Bana h, (E/KerT, ||.||E/KerT ) est un Bana h.

Corollaire 2.6 Soit T ∈ L(E; F ) (linéaire et ontinue). Si ImT est fermée dans F, alors on a la
ondition inf-sup :
hT x, yiF ′ ,F
∃γ > 0, inf ( sup ) ≥ γ. (2.14)
x∈E/KerT y∈ImT ||x||E/KerT ||y||F

Preuve. Immédiat ar (ImT, ||.||F ) est i i fermé, don est un Bana h.

Remarque 2.7 En dimension nie, tout sous-espa e ve toriel est fermé, et don , ImT est fermée
dans F.
C'est faux en dimension innie : un sous-espa e ve toriel n'est pas toujours fermé.
N∗
X
2 2
Exemple. Soit ℓ l'ensemble des suites de arré sommable : ℓ = {(xn ) ∈ R : x2n < ∞}.
n∈N∗
Soit E = F = ℓ2 et T : ℓ2 → ℓ2 l'opérateur de multipli ation par la suite (λn ) = ( n1 ), i.e.
T ((xn )N∗ ) = ( xnn )N∗ pour tout (xn ) ∈ ℓ2 . Cet opérateur est inje tif ar KerT = {0} (immédiat).
2
Et son image est dense dans ℓ (pour une suite (yn = T xn )N∗ dans l'image, onsidérer la suite
2
tronquée dont l'anté édent est une suite nie don dans ℓ ...). Cependant il n'est pas surje tif
1 2
(l'image n'est pas fermée) : par exemple, la suite ( n )N n'a pas d'anté édent dans ℓ .
∗

3 Problème mixte (problème sous ontrainte lassique)

3.1 Le problème : forme faible

Soit (V, (·, ·)V ) et (Q, (·, ·)Q ) deux Hilbert.
Soit a(·, ·) : V × V → R et b(·, ·) : V × Q → R deux formes bilinéaires ontinues.
′ ′
Soit f ∈ V et g ∈ Q (des formes linéaires ontinues sur V et Q).
Problème : trouver (u, p) ∈ V × Q t.q. :

(
a(u, v) + b(v, p) = hf, viV ′ ,V , ∀v ∈ V,
(3.1)
b(u, q) = hg, qiQ′ ,Q , ∀q ∈ Q.

(Exemple : voir (7.5).)

6 21 juin 2018
7 3. Problème mixte (problème sous ontrainte lassique)

3.2 Le problème : forme forte

On note A : V → V ′ l'opérateur asso ié à a(·, ·) :
(
V →V′
A: (3.2)
u → Au, t.q. hAu, viV ′ ,V = a(u, v), ∀(u, v) ∈ V × V.

On note B : V → Q′ l'opérateur asso ié à b(·, ·) :

(
V → Q′
B: (3.3)
v → Bv, t.q. hBv, qiQ′ ,Q = b(v, q), ∀(v, q) ∈ V × Q.

On note B ′ : Q → V ′ l'opérateur dual :

(
′ Q →V′
B : (3.4)
q → B ′ q, t.q. hB ′ q, viV ′ ,V = b(v, q), ∀(v, q) ∈ V × Q,

où don hB ′ q, viV ′ ,V = hBv, qiQ′ ,Q .

Et (3.1) se lit (problème fort asso ié à (3.1)) : trouver (u, p) ∈ V × Q t.q. :

(
Au + B ′ p = f dans V ′,
(3.5)
Bu =g dans Q′ .

(Exemple : voir (7.7).)

Remarque 3.1 p est appelé le multipli ateur de Lagrange de la ontrainte Bu = g . En eet,

Et
soit l'espa e ane Kg = {v ∈ V : Bv = g}, d'espa e ve toriel asso ié K0 = {v ∈ V : Bv = 0} ;
alors (3.5) se lit également : trouver u ∈ Kg t.q. hAu, vi + hp, Bvi = hf, vi pour tout v ∈ V , soit :
trouver u ∈ Kg t.q. hAu, vi = hf, vi pour tout v ∈ K0 , dit problème sous la ontrainte Bu = g .
Et p relaxe la ontrainte Bu = g à l'aide de la te hnique usuelle des multipli ateurs de Lagrange,
f. (3.5) ou (3.1).

3.3 Le problème matri iel asso ié

Le problème matri iel généralisé asso ié à (3.5) ou à (3.1) s'é rira sous une forme de type,
après é riture dans des bases : trouver U, P t.q. :

BT
    
A U F
= , (3.6)
B 0 P G

où BT est la matri e transposée de B.

3.4 Le théorème inf-sup

BT
 
A
Le problème (3.1), ou (3.5), est-il bien posé ? Autrement dit, la matri e asso iée ,
B 0
f. (3.6), est-elle inversible d'inverse bornée (matri e inversible bien onditionnée) ?

Théorème 3.2 . Hypothèses : (i) la forme bilinéaire a(·, ·) est ontinue sur V ×V et oer itive
sur KerB . On note α la onstante de oer ivité, à savoir :

a(v, v)
∃α > 0, inf ≥ α. (3.7)
v∈KerB ||v||2V

(ii) la forme bilinéaire b(·, ·) est ontinue sur V × Q, et B est surje tive sur Q′ , don ( ondition
inf-sup (2.11)) :
|b(v, q)|
∃β > 0, inf (sup ) ≥ β. (3.8)
v∈V /KerB q∈Q ||~v ||V /KerB ||q||Q
(iii) f ∈ V ′, et g ∈ ImB (la solution doit satisfaire Bu = g ).

7 21 juin 2018
8 4. Opérateurs divergen e, gradient, rotationnel

Con lusion : le problème (3.1) a une unique solution (u, p) ∈ V × Q/KerB ′ qui vérie ( ara tère
||a||
bien posé), ave Ca = (1 +
α ) :

1 Ca

 ||u||V ≤ α ||f ||V ′ + β ||g||Q′ ,


(3.9)
C  ||a|| 
 ||p||Q/KerB ′ ≤ a ||f ||V ′ + ||g||Q′ .


β β

Preuve. Soit ug ∈ V Bug = g (existe ar g ∈ ImB ). Soit u0 ∈ KerB la solution du problème :

t.q.
trouver u0 ∈ KerB a(u0 , v0 ) = hf, v0 iV ′ ,V − a(ug , v0 ) pour tout v0 ∈ KerB . Ayant a(·, ·)
t.q.
ontinue et oer itive sur KerB , et ayant F : v0 ∈ KerB → F (v0 ) = hf, v0 iV ′ ,V − a(ug , v0 ) linéaire
et ontinue sur KerB ( ar |F (v0 )| ≤ (||f ||V ′ + ||a|| ||ug ||V )||v0 ||V ), le théorème de LaxMilgram
dans (KerB, (·, ·)V ) donne l'existen e et l'uni ité.
Soit u = u0 + ug . Don a(u, v0 ) = hf, v0 iV ′ ,V . Et u est indépendant du relèvement hoisi ug ,
′ ′ ′ ′ ′
ar si ug est un autre relèvement de g , si u0 ∈ KerB est la solution orrespondante, si u = u0 + ug ,
′ ′ ′ ′ ′
alors u − u = u0 − u0 + ug − ug ∈ KerB ( ar B(ug − ug ) = g − g = 0) et a(u − u , v0 ) = 0 pour
′ ′
tout v0 ∈ KerB , ar a(·, ·) est oer itif sur KerB . Don u − u = 0, don u = u .
De plus a(u, v0 ) = hf, u0 iV ′ ,V pour tout v0 ∈ KerB , f. (3.1), don a(u0 + ug , u0 ) = hf, u0 iV ′ ,V
2
donne a(u0 , u0 ) = hf, u0 iV ′ ,V − a(ug , u0 ), d'où α||u0 ||V ≤ ||f ||V ′ ||u0 ||V + ||a|| ||ug ||V ||u0 ||V , don
1
||u0 ||V ≤ α (||f ||V ′ + ||a|| ||ug ||V ). Ave b(ug , q) = hg, qiQ′ ,Q pour tout q ∈ Q, et il existe q ′
′ ′
t.q. b(ug , q ) ≥ β||ug ||V /KerB ||q ||Q , f. (3.8) (quitte à rempla er β par β + ε où ε > 0). D'où
β||ug ||V /KerB ||q ||Q ≤ hg, q iQ′ ,Q ≤ ||g||Q′ ||q ′ ||Q . D'où ||ug ||V /KerB ≤ β1 ||g||Q′ . D'où ||u0 ||V ≤
′ ′

1 ||a|| 1 ||a|| 1 ||a|| 1

α (||f ||V + β ||g||Q ). Ave α β = β α ≤ β Cα . D'où (3.9)1 .
′ ′

◦ ′
Soit (KerB) = {m ∈ V : hm, v0 iV ′ ,V = 0, ∀v0 ∈ KerB} (l'orthogonal de KerB au sens de
la dualité). Soit L(v) =
déf
a(u, v) − hf, viV ′ ,V : on veut trouver un p ∈ Q t.q. L(v) = b(v, p) =
hB ′ .p, viV ′ ,V , i.e. t.q. L = B ′ .p. La forme linéaire L est nulle sur KerB . Don L ∈ (KerB)◦ .
◦ ′
Ave (KerB) = ImB ′ (propriété lassique entre noyau de B et image de B de démonstration
′
élémentaire, ave l'orthogonal qui est ferme, de démonstration élémentaire). Et ave ImB ′ = ImB
′ ′
(hypothèse image fermée = ondition inf-sup). Don L ∈ ImB ⊂ V , et il existe p ∈ ImB ⊂ Q t.q.
′ ′
L = B .p (unique à un élément de KerB près). D'où (3.9)2 à l'aide de (3.8).

3.5 Le lagrangien
Dans le as où a(·, ·) est symétrique, on note :

1
L(v, q) = a(v, ~v ) + b(v, q) − (f, v)L2 − (g, q)L2 , (3.10)
2
et L est appelé le lagrangien du problème sous ontrainte (3.1).
Le problème d'optimisation asso ié est : trouver (u, p) ∈ V × Q (point selle) t.q. :

L(u, p) = inf (sup L(v, q)). (3.11)

v∈V q∈Q

Véri ation qu'une solution de (3.11) est bien solution de (3.1) : on a (ave a(·, ·) symétrique) :

∂L L(u+hv, p) − L(u, , )
 ∀v ∈ V, (u, p).v = lim = a(u, v) + b(u, q) − hf, vi,


∂v h→0 h
(3.12)
∂L L(u, p+hq) − L(u, p)
 ∀q ∈ Q, (u, p).q = lim = b(u, q) − hg, qi.


∂q h→0 h
Et si (u, p) est solution de (3.11), don extremum de L, alors on retrouve (3.1).

4 Opérateurs divergen e, gradient, rotationnel

Le gradient, la divergen e et le rotationnel sont des opérateurs de dérivation usuels.

Soit (~ei ) une base eu lidienne de Rn . Soit (·, ·)Rn le produit s alaire eu lidien asso ié. On note
~ Rn =noté ~v · w
(~v , w) ~.
n
Soit Ω un ouvert orienté régulier borné de R , de bord Γ.
En x ∈ Γ, on note ~ n(x) le ve teur (·, ·)Rn -normal unitaire sortant.
Soit F (E; F ) l'ensemble des fon tions entre deux ensembles E et F .

8 21 juin 2018
9 5. Cadre fon tionnel

Dénition 4.1 Le gradient est formellement l'opérateur de dérivation donné par :

F (Ω; R) → Rn



~ : n
grad ~
X ∂f (4.1)
 f → gradf = ~ei .

i=1
∂xi

Dénition 4.2 La divergen e est formellement l'opérateur de dérivation donné par :

F (Ω; Rn ) → R



n n
div : X X ∂v i (4.2)
 ~v = v i~ei → div~v = .

i=1
∂xi i=1

Dénition 4.3 Le rotationnel est formellement l'opérateur de dérivation dans R3 :

F (Ω; R3 ) → R3



~ : n
curl X
~ v = ( ∂v3 − ∂v2 )~e1 + ( ∂v1 − ∂v3 )~e2 + ( ∂v2 − ∂v1 )~e3 . (4.3)
 ~v =
 v i~ei → curl~
i=1
∂x2 ∂x3 ∂x3 ∂x1 ∂x1 ∂x2
∂v2 ∂v1
(Dans R2 le rotationnel se réduit à la troisième omposante et est noté curl~v = ∂x1 − ∂x2 .)

5 Cadre fon tionnel

5.1 Espa es de Hilbert utilisés dans la suite

p
Dans un Hilbert (E, (·, ·)E ) la norme utilisée sera
R R ||~v ||E = (~v , ~v )E .
L'é riture générique ondensée
Ω f g dΩ signie x∈Ω f (x)g(x) dΩ.
n
Soit Ω un ouvert dans R . On note :
Z Z
2 2
L (Ω) = {f : Ω → R : f dΩ < ∞}, et (f, g)L2 = f g dΩ. (5.1)
Ω Ω
P
Soit (~ei ) une base eu lidienne. Ave ~u = i ui~ei on note :
n
Z X Z
n
L2 (Ω) = {~u : Ω → Rn : u2i dΩ < ∞}, et (~u, ~v )L2 = ~u · ~v dΩ. (5.2)
Ω i=1 Ω

2n 2 R Pn
Idem pour L2 (Ω) = {F : Ω → Rn : Ω
2
i=1 Fij dΩ < ∞}. En parti ulier ave d~u de matri e
i 2
∂u 2 ∂ f
gradient [
∂xj ] et d f de matri e hessienne [
∂xi ∂xj ] on note :
n n
∂ui ∂v i ∂2f ∂2g
Z X Z X
(d~u, d~v )L2 = dΩ, et (d2 f, d2 g)L2 = dΩ. (5.3)
Ω i,j=1 ∂xj ∂xj Ω i,j=1 ∂xi ∂xj ∂xi ∂xj

Puis on note :

~ n ~ ~
H 1 (Ω) = {f ∈ L2 (Ω) : gradf ∈ L2 (Ω) }, (f, g)H 1 = (f, g)L2 + (gradf, gradg)L2 .
2
n
H 2 (Ω) = {f ∈ H 1 (Ω) : d2 f ∈ L2 (Ω) }, (f, g)H 2 = (f, g)H 1 + (d2 f, d2 g)L2 .
n
H div (Ω) = {~v ∈ L2 (Ω) : div~v ∈ L2 (Ω)}, (~u, ~v )H div = (~u, ~v )L2 + (div~u, div~v )L2 .
curl 2 3
~ v ∈ L (Ω) }, 2 3 ~ u, curl~
~ v )L2 .
H (Ω) = {~v ∈ L (Ω) : curl~ (~u, ~v )H curl = (~u, ~v )L2 + (curl~

5.2 Sous-espa es de Hilbert

Les adhéren es de D(Ω) = Cc∞ (Ω), l'espa e des fon tions C ∞ (Ω) à support ompa t dans Ω,
sont :

H1
H01 (Ω) = {f ∈ H 1 (Ω) : f|Γ = 0} = D(Ω) , ~
(f, g)H01 = (gradf, ~
gradg)L2 .
2
H
H02 (Ω) = {f ∈ H 2 (Ω) : f|Γ = 0 et ~
gradf.~
n|Γ = 0} = D(Ω) , (f, g)H02 = (d2 f, d2 g)L2 .
H div
H0div (Ω) = {~v ∈ H div (Ω) : (~v · ~n)|Γ = 0} = D(Ω)n , (~u, ~v )H0div = (div~u, div~v )L2 .
H curl
H0curl (Ω) = {~v ∈ H curl (Ω) : (~v ∧ ~n)|Γ = 0} = D(Ω)3 , ~ u, curl~
(~u, ~v )H0curl = (curl~ ~ v )L2 .

9 21 juin 2018
10 5. Cadre fon tionnel

Pour Ω borné, dans es ensembles les semi-produits s alaires donnés sont des produits s alaires
équivalents aux produits s alaires des espa es non indi és par 0.

1
5.3 Appli ation tra e γ0 et l'Hilbert H 2 (Γ)
L'appli ation tra e γ0 est l'appli ation linéaire et ontinue (voir ours d'éléments nis) :
( )
H 1 (Ω) → L2 (Γ),
γ0 : . (5.4)
f → γ0 (f ) = f|Γ ,
n n
Et on garde la notation γ0 pour γ0 : H 1 (Ω) → L2 (Γ) , et γ0 (~v ) = ~v|Γ .
1,1
On a, supposant Ω régulier (Γ lips hitzienne C ) :

H01 (Ω) = Ker(γ0 ). (5.5)

Et pour on note :

1
Im(γ0 ) = H 2 (Γ), et ||d|| 1 = inf ||u||H 1 . (5.6)
H 2 (Γ) u∈H 1 (Ω):u|Γ =d

Proposition 5.1 On a :
||d|| 1 = ||ud ||H 1 , (5.7)
H 2 (Γ)
1
où ud ∈ H 1 (Ω) est la solution du problème de Diri hlet : pour d ∈ H 2 (Γ), trouver u ∈ H 1 (Ω) t.q. :
( )
−∆u + u = 0 dans Ω,
et (ud , v0 )H 1 = 0, ∀v0 ∈ H01 (Ω). (5.8)
u|Γ = d sur Γ,

Et ||.|| 1 dérive du produit s alaire :

H 2 (Γ)

(c, d) 1 = (uc , ud )H 1 (Ω) (5.9)

H 2 (Γ)

1
Et (H 2 (Γ), (·, ·) 1 ) est un Hilbert.
H 2 (Γ)
1
1
De plus γ0 : v ∈ (H (Ω), ||.||H 1 (Ω) ) → v|γ ∈ (H 2 (Γ), ||.|| 1 ) est (linéaire) ontinue.
H 2 (Γ)

Preuve. Soit zd ∈ H 1 (Ω) un relèvement de d ∈ H 1

(Γ) (un anté édent par l'appli ation tra e γ0 ),
2

i.e.zd ∈ H 1 (Ω) t.q. γ0 (zd ) = d. Et on note zd|Γ = d). On her he une solution u de (5.8) sous la
1 1
forme u = u0 + zd ∈ H0 (Ω) + zd (espa e ane passant par zd , d'espa e ve toriel asso ié H0 (Ω)).
1
Don u0 ∈ H0 (Ω) est solution de :

−∆u0 + u0 = ∆zd − zd dans H −1 (Ω),

( )
1 (5.10)
u0|Γ = 0 dans H 2 (Γ).

Le problème faible asso ié est : trouver u ∈ H01 (Ω) t.q. :

(u0 , v0 )H 1 (Ω) = −(zd , v0 )H 1 (Ω) , ∀v0 ∈ H01 (Ω), (5.11)

et le théorème de LaxMilgram donne l'existen e d'une solution Et la somme ud = u0 ∈ H01 (Ω).

′
u0 + zd est indépendante du relèvement zd hoisi : pour un se ond relèvement zd on trouve une
′ ′ ′ ′ ′ 1
se onde solution u0 , et la somme ud = u0 + zd vérie (ud − ud , v0 )H 1 (Ω) = 0 pour tout v0 ∈ H0 (Ω),
′ 1 ′ 1
ave ud − ud ∈ H0 (Ω), don ud − ud = 0 (LaxMilgram). Et on a, pour tout v0 ∈ H0 (Ω) :

||ud + v0 ||2H 1 (Ω) = ||ud ||2H 1 (Ω) + ||v0 ||2H 1 (Ω) + 2(ud , v0 )H 1 (Ω) = ||ud ||2H 1 (Ω) + ||v0 ||2H 1 (Ω) + 0,

1 ⊥ 1
f. (5.11). (Et (5.11) indique que ud ∈ H0 (Ω) H .)
2 2
D'où inf w∈H 1 (Ω) ||w||H 1 (Ω) = inf v0 ∈H 1 (Ω) ||ud + v0 ||H 1 (Ω)
0
= ||ud ||2H 1 (Ω) . Et on pose (5.9), et on
wΓ =d

vérie que (·, ·) 1 dénit une forme bilinéaire symétrique dénie positive (immédiat).
H2
1 1 1
On en déduit que (H 2 (Γ), ||.|| 1) est fermé : si (dn )N∗ ∈ H 2 (Γ) est de Cau hy dans H 2 (Γ),
H2
alors (udn )N∗ est de Cau hy dans (H 1 (Ω), ||.||H 1 ) omplet, don onvergente dans H 1 (Ω) vers un

10 21 juin 2018
11 5. Cadre fon tionnel

1
u ∈ H 1 (Ω). Notons d = γ0 (u) ∈ H 2 (Γ). Comme (udn , v0 )H 1 (Ω) = 0 pour tout v0 ∈ H01 (Ω), f. (5.8),
1
on a (ud , v0 )H 1 (Ω) = 0 pour tout v0 ∈ H0 (Ω) par ontinuité du produit s alaire. Ave γ0 (ud ) = d,
don ud est solution de (5.8). Don (5.8) donne ||d − dn || 1 = ||u − udn ||H 1 −→n→∞ 0.
2 H
1
Et l'appli ation linéaire γ0 : H 1 (Ω) → H 2 (Γ) vérie ||γ0 (v)|| 1 = ||ud ||H 1 (Ω) ≤ ||ud +
H 2 (Γ)
v0 ||H 1 (Ω) pour tout v0 ∈ H01 (Ω), et parti ulier ||γ0 (v)|| 1 ≤ ||v||H 1 (Ω) .
H 2 (Γ)

5.4 Autres appli ations tra es ontinues

 2 2
 H (Ω) → L (Γ),
γ1 : noté ∂f (5.12)
 ~ )|Γ · ~n =
f → γ1 (f ) = (gradf .
∂~n |Γ
− 12
H div (Ω) → H
(
(Γ),
γn : noté
(5.13)
~v → γn (~v ) = γ0 (~v ) · ~n = (~v · ~n)|Γ ,
et la divergen e permet d'avoir le ontrle de la omposante normale ~v .~n sur Γ.

 H curl(Ω) → H − 12 (Γ)3 ,
~γt : (5.14)
noté
 ~v → ~γ (~v ) = γ (~v ) ∧ ~n = (~v ∧ ~n) ,
t 0 |Γ

et le rotationnel permet d'avoir le ontrle de la omposante tangentielle ~v ∧ ~n sur Γ.

5.5 Dual et Bana h

Bana h :

(L2 (Ω))′ ≃ L2 (Ω) (identi ation usuelle) ||f ||L2 (Ω)′ = ||f ||L2 (Ω) .
Z
L20 (Ω) = {f ∈ L2 (Ω) : f dΩ = 0} ≃ L2 (Ω)/R, ||f ||L20 = inf ||f + c||L2 .
Ω c∈R

′ hf, vi (5.15)
H −1 (Ω) = H01 (Ω) , ||f ||H −1 = sup .
v∈H01 (Ω) ||v||H01
1 1 |hµ, λi|
H − 2 (Γ) = (H 2 (Γ))′ , ||µ|| 1 = sup .
H − 2 (Γ) 1 ||λ|| 12
λ∈H (Γ) 2 H (Γ)

2 ′ 2
On a identié (L (Ω)) ave L (Ω) à l'aide du théorème de représentation de Riesz, à savoir :
′
∀ℓ ∈ L (Ω) , ∃!f ∈ L (Ω), t.q. ∀g ∈ L2 (Ω), hℓ, gi = (f, g)L2 (Ω) , et on a ||ℓ||L2 (Ω)′ = ||f ||L2 (Ω) . Don
2 2
′
H01 (Ω) ⊂ H 1 (Ω) ⊂ L2 (Ω) ≃ L2 (Ω) ⊂ (H 1 (Ω))′ ⊂ H −1 (Ω).
L2 (Ω)/R est l'ensemble des fon tions L2 (Ω) dénies à une onstante près, identié à L20 (Ω). Et
2
la norme ||.||L2 donnée est la norme quotient (norme usuelle de L (Ω)/R).
0

Proposition 5.2 Pour λ ∈ H − 2 (Γ)

1
on a :

||λ|| 1 = ||wλ ||H 1 (Ω) , (5.16)

H − 2 (Γ)

où wλ ∈ H 1 (Ω) est la solution du problème de Neumann : trouver w ∈ H 1 (Ω) t.q. :

 
 − ∆w + w = 0 dans Ω, 
soit (w, v)H 1 (Ω) = hλ, vi , ∀v ∈ H 1 (Ω). (5.17)
 ∂w = λ sur Γ,
1 1
 H − 2 (Γ),H 2 (Γ)
∂n

Preuve.
1
Soit λ ∈ H 2 (Γ), et soit wλ ∈ H 1 (Ω) la solution de (5.17) (LaxMilgram : la
1
forme bilinéaire a(u, v) = (u, v)H 1 (Ω) est H (Ω)- ontinue et oer itive, et la forme linéaire

11 21 juin 2018
12 5. Cadre fon tionnel

ℓ : H 1 (Ω) → R donnée par ℓ(v) = hλ, vi

1 1 vérie |ℓ(v)| ≤ ||λ|| 1 ||γ0 (v)|| 1 ≤
H − 2 (Γ),H 2 (Γ) H − 2 (Γ) H 2 (Γ)
||λ|| − 12 ||v||H 1 (Ω) , f. prop. 5.1, don ℓ est ontinue). Et :
H (Γ)

|hλ, di 1 1 |
H − 2 (Γ),H 2 (Γ)
||λ|| 1 = sup (dénition)
H− 2 (Γ) 1 ||d|| 1
d∈H 2 (Γ) H 2 (Γ)

|hλ, ud i 1 1 |
H − 2 (Γ),H 2 (Γ)
= sup ( f. (5.9))
1 ||ud ||H 1 (Ω) (5.18)
d∈H 2 (Γ)

|(wλ , ud )H 1 (Ω) |
= sup ( f. (5.17))
1 ||ud ||H 1 (Ω)
d∈H 2 (Γ)

≤ ||wλ ||H 1 (Ω) , (Cau hyS hwarz dans H 1 (Ω)).

1
Et ave d = γ0 (wλ ) ∈ H 2 (Γ), on a

|hλ, γ0 (wλ )i 1 1
H − 2 (Γ),H 2 (Γ)
||λ|| 1 ≥ | (par dénition du sup)
H− 2 (Γ) ||γ0 (wλ )|| 1
H 2 (Γ)
|(wλ , wλ iH 1 (Ω) | (5.19)
≥ ( f. (5.17))
||γ0 (wλ )|| 12
H (Γ)

≥ ||wλ ||H 1 (Ω) ( ar ||wλ ||H 1 (Ω) ≥ ||γ0 (wλ )|| 1 f. (5.6)).
H 2 (Γ)

D'où (5.16).

5.6 Intégrations par parties (IPP)

Pour f ∈ H 1 (Ω) et ~v ∈ H div (Ω), on a :
Z Z Z
~
gradf · ~v dΩ = − f div~v dΩ + f ~v · ~n dΓ. (5.20)
Ω Ω Γ

Pour ~v ∈ H 1 (Ω) et ~ ∈ H curl (Ω), on a :

w
Z Z Z
~ v·w
curl~ ~ w
~ dΩ = + ~v · curl ~ dΩ + ~v · (w
~ ∧ ~n) dΓ. (5.21)
Ω Ω Γ

5.7 Dual de H 1 (Ω) et H div (Ω) ( ara térisations)

Théorème 5.3 Dual de H 1 (Ω).
n ~
ℓ ∈ (H 1 (Ω))′ ⇒ ∃(f, ~u) ∈ L2 (Ω) × L2 (Ω) : hℓ, ψi = (f, ψ)L2 + (~u, gradψ)L2 ∀ψ ∈ H 1 (Ω).
(5.22)
Dual de H01 (Ω).
n
ℓ ∈ H −1 (Ω) ⇒ ∃(f, ~u) ∈ L2 (Ω) × L2 (Ω) t.q. ℓ = f − div~u. (5.23)

Et on peut prendre f =0 quand Ω est borné.

Preuve. (On suit Brézis.) Cara térisation de

n
H 1 (Ω) . Soit Z = L2 (Ω) × L2 (Ω)
n
muni du
produit s alaire
( 1 ((ϕ, ~u), (ψ, ~v ))
)Z = (ϕ, ψ)L2 + (~u, ~v )L2 . Ainsi T est un Hilbert. Soit T :
H (Ω) → Z
, où don ||ψ||H 1 = ||T ψ||Z = ||(ψ, gradψ)|| ~ Z , et don T est une
~
ψ → T ψ = (ψ, gradψ)
isométrie sur son image ImT muni de la norme ||.||Z .
1 ′
Soit ℓ ∈ H (Ω) . La fon tion :
(
ImT → R
ΦImT :
~
(ψ, ~v =gradψ) → hΦImT , (ψ, ~v )iZ ′ ,Z = hℓ, T −1 (ψ, ~v )iH 1 ′ ,H 1 = hℓ, ψiH 1 ′ ,H 1
est une forme linéaire ontinue sur ImT ( ar ℓ et T −1 sont ontinues). Grâ e au théorème de Hahn
( )
Z →R
Bana h, on peut la prolonger à tout Z en une forme linéaire ontinue ΦZ :
(ψ, ~v ) → hΦZ , (ψ, ~v )i

12 21 juin 2018
13 5. Cadre fon tionnel

(ave ||ΦZ ||Z ′ = ||ℓ||). Et on peut don appliquer le R théorème de

R représentation de Riesz : il existe
(ϕ, ~u) ∈ Z t.q. hΦZ , (ψ, ~v )i = ((ϕ, ~u), (ψ, ~v ))Z = Ω ϕ ψ dΩ + Ω ~u.~v dΩ pour tout (ψ, ~v ) ∈ Z . En
parti ulier pour (ψ, ~ ~
v =gradψ) ∈ ImT , d'où (5.22).
Démonstration similaire pour (5.23).

Théorème 5.4 Dual de H div (Ω).

n
F ∈ (H div (Ω))′ ⇒ ∃(f~, ϕ) ∈ L2 (Ω) × L2 (Ω) : hF, ~v i = (f~, ~v )L2 + (ϕ, div~v )L2 , ∀~v ∈ H div (Ω).
(5.24)
Dual de H0div (Ω). En parti ulier :

′ n
F ∈ H0div (Ω) ⇒ ∃(f~, ϕ) ∈ L2 (Ω) × L2 (Ω) t.q. F = f~ − gradϕ.
~ (5.25)

Et on peut prendre f~ = 0 quand Ω est borné.

Preuve. (Similaire à la démonstration de (5.22).) Soit Z = L2 (Ω)n ×(L2(Ω) muni du produit s a- )

H div (Ω) → Z
laire ((~u, p), (~v , q))Z = (~u, ~v )L2 +(p, q)L2 qui en fait un Hilbert. Soit T : ,
~v → T ~v = (~v , div~v )
où don ||~v ||H div = ||T ~v||Z = ||(~v , div~v )||Z , et don T est une isométrie sur son image ImT muni
de la norme ||.||Z .
div ′
Soit F ∈ H (Ω) . La fon tion (~v , q=div~v ) ∈ ImT → hF, T −1 (~v , q)iH div ′ ,H div = hF, ~v iH div ′ ,H div
−1
est une forme linéaire ontinue sur ImT ( ar F et T sont ontinues). Grâ e au théorème de Hahn
Bana h, on peut la prolonger à tout Z en une forme linéaire ontinue Φ : (~ v , q) ∈ Z → hΦ, (~v , q)i
(ave ||Φ||Z ′ = ||F ||). Et le théorème de représentation de Riesz donne : il existe (~ u, p) ∈ Z t.q.
R R
hΦ, (~v , q)i = ((~u, p), (~v , q))Z = Ω ~u.~v dΩ + Ω pq dΩ pour tout (~v , q) ∈ Z . En parti ulier pour
(~v , q=div~v ) ∈ ImT , d'où (5.24).
Démonstration similaire pour (5.25).

5.8 Noyaux des appli ations tra es

Ω est supposé régulier.

1
Kerγ0 = H01 (Ω), Im(γ0 ) = H 2 (Γ) dense dans L2 (Γ).
1
Kerγ1 ∩ Kerγ0 = H02 (Ω), Im(γ1 ) = H 2 (Γ).
1 (5.26)
Kerγn = H0div (Ω), Im(γn ) = H − 2 (Γ).
1 3
Ker~γt = H0curl (Ω), Im(~γt ) = H − 2 (Γ) .

5.9 Poin aréFriedri hs

(Voir poly opié Eléments nis, ou RaviartThomas [27℄, ou Ciarlet...
Si Ω est borné (au moins dans une dire tion), on dispose de l'inégalité de Poin aré dans H01 (Ω) :
il existe une onstante cΩ telle que :

∀v ∈ H01 (Ω), ~
||v||L2 ≤ cΩ ||gradv|| L2 , (5.27)

et les normes ~
||v||H 1 (Ω) et ||gradv|| 1
L2 (Ω) sont équivalentes dans H0 (Ω) (qui est omplet par déni-
1
tion : adhéren e de D(Ω) dans H (Ω)).
Si Ω est borné alors il existe une onstante CΩ telle que :
\
∀v ∈ H01 (Ω) H 2 (Ω), ||v||H 2 (Ω) ≤ CΩ ||∆v||L2 (Ω) . (5.28)

H01 (Ω) H 2 (Ω)

T
et les normes ||v||H 2 (Ω) et ||∆v||L2 (Ω) sont équivalentes dans (qui n'est pas fermé
1
dans H (Ω)).

13 21 juin 2018
14 6. Les surje tivités de l'opérateur divergen e

n
5.10 Dé omposition de L2 (Ω) (Helmholtz simple)
Ave div : H div (Ω) → L2 (Ω), notons :
(
Ker(div) = {~v ∈ H div (Ω) : div~v = 0},
(5.29)
Ker(div)0 = Ker(div) ∩ H0div (Ω) = {~v ∈ Ker(div) : (~v .~n)|Γ = 0},

l'ensemble des fon tions in ompressibles, et son sous ensemble orrespondant à Γ imperméable.

Théorème 5.5 Soit Ω ⊂ Rn un ouvert régulier borné. On a ave div : H div (Ω) → L2 (Ω) :


⊥L 2
 L2 (Ω)n = grad(H
 ~ 1
(Ω)) ⊕ Ker(div),
0
(5.30)
n L ⊥ 2
2 ~ 1

L (Ω) = grad(H (Ω)) ⊕ Ker(div)0 ,


n
i.e, pourf~ ∈ L2 (Ω) il existe (ϕ, w)
~ respe tivement dans ~
grad(H 1
0 (Ω)) × Ker(div) et dans
~ 1
grad(H (Ω)) × Ker(div)0 t.q. :

f~ = gradϕ
~ + w,
~ ave ~
(gradϕ, w)
~ L2 = 0. (5.31)

De plus ave :
∃C > 0, ~
||gradϕ|| ~ L2 ≤ C||f~||L2 .
L2 (Ω) + ||w|| (5.32)

Preuve. Soit
n
f~ ∈ L2 (Ω) .
ϕ la solution du
Soit problème de Diri hlet homogène dans H01 (Ω) : trouver ϕ ∈ H01 (Ω) t.q.
∆ϕ = divf~. Formulation ~ ~ ~ ~
gradϕ.gradψ dΩ = Ω f .gradψ dΩ pour tout ψ ∈ H01 (Ω). On a
R R
faible :
Ω
n
1
existen e et uni ité de la solution dans H0 (Ω) (LaxMilgram). On note w ~ = f~ − gradϕ
~ ∈ L2 (Ω) ,
et on a w ~
~ ⊥L2 grad(H 1
~ = divf~ − ∆ϕ, ψ = 0, d'où (5.30)1 .
0 (Ω)). Et divw
~
Et ||gradϕ|| L2 (Ω) ≤ ||f ||L2 (prendre ψ = ϕ), don ||w||
~
~ L2 ≤ ||f ||L2 + ||gradϕ|| L2 ≤ 2||f ||L2 ,
d'où (5.32) ave C = 3.
1 1
Soit ϕ une solution du problème de type Neumann homogène dans H (Ω) : trouver ϕ ∈ H (Ω)
~ ~ ~ ~ 1
R R
t.q.
Ω gradϕ.gradψ dΩ = Ω f .gradψ dΩ pour tout ψ ∈ H (Ω). On a existen e et uni ité de la
1
solution dans H (Ω)/R (LaxMilgram), où don ϕ est dénie à une onstante près, et gradϕ ~ est
n
unique. On note w ~ = f~ − gradϕ
~ ∈ L2 (Ω) , et on a w ~
~ ⊥L2 grad(H 1
(Ω)). Et ave ψ ∈ D(Ω) on
obtient (au sens des distributions) divw ~ = divf~ − ∆ϕ = 0 dans H −1 (Ω), don divw ~ = 0 dans L2 (Ω)
2 div ~
~ gradψ dΩ = 0 pour tout ψ ∈ H 1 (Ω),
R
( ar 0 ∈ L (Ω)), don w ~ ∈ H (Ω) et w ~ ∈ Ker(div), ave Ω w.
1 1
~ n) ψ dΓ = 0 pour tout ψ ∈ H 1 (Ω), don w.~~ n = 0 dans H − 2 (Γ) ( ar H 2 (Γ) est
R
et don (IPP)
Ω
(w.~
2
dense dans L (Γ)). Don w ~ ∈ Ker(div)0 et (5.30)2 .
~
Et ||gradϕ||L2 (Ω) ≤ ||f ||L2 (prendre ψ = ϕ), don ||w|| ~
~ L2 ≤ ||f ||L2 + ||gradϕ|| L2 ≤ 2||f ||L2 ,
d'où (5.32) ave C = 3.

Remarque 5.6 On a :
n ⊥H 1
H01 (Ω) = Ker(div) + Ker(div) 0 , (5.33)

⊥H 1 ~ ~ n
et Ker(div) 0 = ∆−1 (grad(L 2
(Ω))) = {~v = ∆−1 (gradq), q ∈ L2 (Ω)}, où ∆−1 : H −1 (Ω) →
n
H01 (Ω) est l'opérateur de résolution du problème de Diri hlet ∆~u = ℓ ave ~u|Γ = ~0.

6 Les surje tivités de l'opérateur divergen e

Même si la ontrainte est de type div~u = 0, on a besoin de regarder le as div~u = g pour

g 6= 0 : ela permettra de voir si le problème est bien posé (matri e asso iée bien onditionnée).
Plus pré isément on a besoin de savoir si Im(div) est fermé dans un espa e donné.

14 21 juin 2018
15 6. Les surje tivités de l'opérateur divergen e

6.1 Surje tivité de l'opérateur divergen e div : H div (Ω) → L2 (Ω)

( )
H div (Ω) → L2 (Ω)
Théorème 6.1 Soit Ω borné. L'opérateur div : est ontinu et surje tif.
~v → div~v ,
Et don :
|(div~v , p)L2 |
∃β > 0, inf sup ≥ β. (6.1)
v ∈H div (Ω) ||~v ||H div ||p||L2
p∈L 2(Ω) ~

Preuve. Comme ||div~ v ||L2 ≤ ||~v ||H div , la ontinuité est immédiate.
Soitf ∈ L2 (Ω). Soit p ∈ H01 (Ω) la solution de ∆p = f (Diri hlet homogène), f. LaxMilgram.
On a ~
div(gradp) = f ∈ L2 (Ω). Don ~v = gradp ~ ∈ H div (Ω), et div~v = f . D'où la surje tivité.
Et (6.1) est donné par (3.8).

6.2 Surje tivité de l'opérateur divergen e H0div (Ω) → L20 (Ω)

( )
H0div (Ω) → L20 (Ω)
Théorème 6.2 Soit Ω borné. L'opérateur : div : est ontinu et surje tif.
~v → div~v ,
Et don :
|(div~v , p)L2 |
∃β > 0, inf sup ≥ β. (6.2)
p∈L20 (Ω) ~
v ∈H div (Ω) 0
||~v ||H0div ||p||L20

Preuve. ~v ∈ H0div (Ω), don Ω div~v dΩ = Γ ~v .~n dΓ = 0, don div~v ∈ L20 (Ω).
R R
Soit
Comme ||div~ v ||L2 = ||~v ||H0div ≤ ||~v ||H0div , la ontinuité est immédiate.
2 1
Soit f ∈ L0 (Ω). Soit p ∈ H (Ω)/R la solution de (gradp, ~ ~
gradq)L2 = (f, q)L2 pour tout q ∈
1
H (Ω), f. LaxMilgram (Neumann homogène, i i p est dénie à une onstante près, et pour q = 1
2
on a (f, 1)L2 = 0 = la ondition de ompatibilité f ∈ L0 (Ω)). Et div(gradp) ~ = f ∈ L2 (Ω) ave
~ ~
(gradp.~n)|Γ = 0. Don gradp = noté div
~v ∈ H (Ω), ave (~v .~n)|Γ = 0. D'où la surje tivité.
Et (6.1) est donné par (3.8).

n
6.3 Surje tivité de l'opérateur divergen e L2 (Ω) → H −1 (Ω)
Théorème 6.3 Soit Ω borné. div : L2 (Ω) → H −1 (Ω)
n
est ontinu et surje tif.
Et ~
hdiv~v , piH −1 ,H01 = (~v , gradp)L2 donne :

~
|(~v , gradp)L2 |
∃β > 0, inf1 sup ≥ β. (6.3)
p∈H0 (Ω) ~
v ∈L2 (Ω)n ||v||L2 (Ω) ||p||H01

~
Preuve. Continuité : ||div~u||H −1 = supϕ∈H01 (Ω) |hdiv~u,ϕ|i
||ϕ||H 1 = supϕ∈H01 (Ω)
|(~
u,gradϕ)
||ϕ||H 1
L2
≤ ||~u||L2
0 0
2
(Cau hyS hwarz dans L (Ω)).
−1 n
Soit ℓ ∈ H (Ω) : de la forme ℓ = f + div~u où f ∈ L2 (Ω) et ~u ∈ L2 (Ω) , f. (5.23).
n
On pose w~ ∈ H div (Ω) t.q. divw~ = f , f. thm. 6.1. D'où ℓ = div(w ~ ∈ L2 (Ω) .
~ + ~u) ave ~u + w
D'où la surje tivité.
Et (6.3) est donné par (3.8).

n
6.4 Surje tivité de l'opérateur divergen e div : H01 (Ω) → L20 (Ω)
Théorème 6.4 Soit Ω borné. L'opérateur divergen e déni sur H01 (Ω)
n
est ontinu et surje tif
sur L20 (Ω) :
n
( )
H01 (Ω) → L20 (Ω)
div : est ontinu et surje tif. (6.4)
~v → div~v ,
n
~v ∈ H01 (Ω) Im(div) ⊂ L20 (Ω).)
R R
(Noter que
Ω
div~v dΩ = Γ
~v .~n dΓ = 0 pour , ar ~v|Γ = 0, don
n
∃β > 0, ∀p ∈ L20 (Ω), ∃~v ∈ H01 (Ω) , (div~v , p)L2 ≥ β ||p||L20 ||~v ||H01 . (6.5)

Soit ( ondition inf-sup) :

div~v p
∃β > 0, inf ( sup ( , )L2 ) ≥ β. (6.6)
2
p∈L0 (Ω) ~
v ∈H01 (Ω)n ||~v ||H01 ||p||L20

15 21 juin 2018
16 6. Les surje tivités de l'opérateur divergen e

Preuve. C'est la surje tivité di ile à établir. Voir par exemple GiraultRaviart [17℄, ou poly-
opié Problèmes spe traux....

( )
H −1 (Ω) → H01 (Ω)
Remarque 6.5 Soit ∆ −1
: l'appli ation linéaire ontinue (Lax
f → ∆−1 f = u
Milgram), où don u ∈ H01 (Ω) et ∆u = f . On a ave (6.4) :

~ )
Ker(div)⊥ = {~v ∈ H01 (Ω) : ~v = ∆−1 (gradf où f ∈ L2 (Ω)}, (6.7)

i.e. ~v ∈ Ker(div)⊥ ssi il existe f ∈ L2 (Ω) t.q. ~ ,

∆~v = gradf i.e. ssi ~ v , grad
(grad~ ~ w)~ L2 (Ω) =
(f, divw) ~ ∈ H01 (Ω).
~ L2 (Ω) pour tout w

Remarque 6.6 Soit Ω borné. Le dual de div donné en (6.4) est l'opérateur gradient :

 2 −1 n 
 L (Ω) → H (Ω)
 

~
grad : n
∂f et ~ = {les
Ker(grad) onstantes}, (6.8)
~
X
 f → gradf = ~ei , 
∂xi
 
i=1

et l'opérateur ~ : L2 (Ω) → H −1 (Ω)n

grad est inje tif.
0

16 21 juin 2018
17 8. Approximation numérique du problème de Stokes

Deuxième partie

Appli ations
7 Le modèle Stokes

7.1 Problème sous ontrainte initial

Soit :
n
V = {~v ∈ H01 (Ω) : div~v = 0} (7.1)
n
muni du produit s alaire (·, ·)H01 qui en fait un Hilbert (i i on onsidère div : H01 (Ω) → L20 (Ω)).
n
Problème initial de type Diri hlet : pour f~ ∈ H −1 (Ω) , trouver ~u ∈ V t.q. :

−∆~u = f~, où don u|Γ = 0. (7.2)

Problème variationnel asso ié : trouver ~u ∈ V t.q. :

(grad~u, grad~v )L2 = (f~, ~v )L2 , ∀~v ∈ V. (7.3)

Le théorème de LaxMilgram indique que e problème est bien posé dans (V, (·, ·)H01 )
Problème de minimum asso ié : ave J(~v ) =déf 12 ||grad~
~ v ||2 2 − (f~, ~v )L2 ,
L trouver u∈V t.q. :

J(~u) = min J(~v ). (7.4)

~
v ∈V

7.2 Problème ave ontrainte relaxée

On impose la ontrainte div~u = 0 à l'aide d'un multipli ateur de Lagrange, noté p : le pro-
blème (7.2) est modifée en : trouver (~u, p) ∈ H01 (Ω) × L20 (Ω) t.q. :
n
(grad~u, grad~v )L2 − (p, div~v )L2 = (f~, ~v )L2 ,
(
∀~v ∈ H01 (Ω) ,
(7.5)
− (div~u, q)L2 = 0, ∀q ∈ L20 (Ω).

(On her he p ∈ L20 (Ω) ar Im(div) = L20 (Ω).) Notons :

( n n
a(~u, ~v ) = (grad~u, grad~v )L2 sur H01 (Ω) × H01 (Ω) ,
n (7.6)
b(~v , q) = −(div~v , q)L2 sur H01 (Ω) × L20 (Ω).
n n
Ainsi (7.5) prend la forme (3.1) ave V = H01 (Ω) , Q = L20 (Ω), g = 0 et B = div : H01 (Ω) → L20 (Ω).
1 n 1 2
I i a(·, ·) est ontinue oer ive sur tout H0 (Ω) , don sur KerB . Et B = div : H0 (Ω) → L0 (Ω)
est ontinue surje tive, f. (6.4). On peut appliquer le théorème 3.2 : on a existen e et uni ité de u
et p ave un ontrle ||; ||H 1 de ~
u et un ontrle ||; ||L20 de p, f. (3.9).
0
Le problème fort asso ié est :
(
~
−∆~u + gradp = f~ ∈ H −1 (Ω),
(7.7)
div~u = 0.
n
Le lagrangien asso ié est, sur H01 (Ω) × L20 (Ω) :

1
L(~v , q) = ||grad~v ||2L2 − (q, div~v )L2 − (f~, ~v )L2 . (7.8)
2

8 Approximation numérique du problème de Stokes

Un des buts de e Ÿ est de omprendre le problème ausé par ertains hoix d'approximations.
On se limite au adre onforme : Vh et Qh sont deux espa es de dimension nie t.q. Vh ⊂ H01 (Ω)
2
et Qh ⊂ L0 (Ω).

17 21 juin 2018
18 8. Approximation numérique du problème de Stokes

8.1 Approximation
Le problème (7.5) dis rétisé est : trouver ~uh ∈ (Vh )n et ph ∈ Q h t.q. :

(
(grad~uh , grad~vh )L2 − (ph , div~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ (Vh )n ,
(8.1)
(div~uh , qh )L2 = 0, ∀qh ∈ Qh ,

8.2 Proje tions (méthode des éléments nis)

( )
L2 (Ω) → Xh
Si Xh est un sous espa e ve toriel de L2 (Ω), on note ΠXh : l'opérateur de
f → ΠXh f
proje tion (·, ·)L2 -orthogonal sur Xh :

∀f ∈ L2 (Ω), ΠXh f ∈ Xh et (ΠXh f, xh )L2 = (f, xh )L2 , ∀xh ∈ Xh , (8.2)

i.e. la fon tion ΠP1 f ∈ Xh est la meilleure approximation Xh de f au sens L2 (Ω). Et on généralise
2 n
ette notation à Xh un sous espa e ve toriel de L (Ω) .
On note :
n
(
H01 (Ω) → Qh
divh = ΠQh ◦ div : (8.3)
~v → divh~v = ΠVh (div~v ),
où don (divh~vh , qh )L2 = (div~v , qh )L2 pour tout qh ∈ Qh . Don (8.1) se lit aussi :
(
(grad~uh , grad~vh )L2 − (ph , divh~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,
(8.4)
(divh ~uh , qh )L2 = 0, ∀qh ∈ Qh ,

On note :
L2 (Ω) → Vh
(
~ déf ~
gradh = ΠVh ◦ grad : (8.5)
~ p = ΠV (gradp).
p → grad ~
h h

où don ~ p, ~vh )L2 = (gradp,

(grad ~ ~vh )L2 pour tout ~vh ∈ Vh . Et ayant Vh ⊂ H01 (Ω), (8.1) se lit aussi :
h
(
~ ph , ~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,
(grad~uh , grad~vh )L2 + (grad h
(8.6)
(divh ~uh , qh )L2 = 0, ∀qh ∈ Qh ,

8.3 Représentation matri ielle

Soit dim Vh = nV et (ϕi )i=1,...,nV une base de Vh (base d'éléments nis).
Soit dim Qh = nQ et (ψi )i=1,...,nQ une base de Qh (base d'éléments nis).P
nQ
Pour ph ∈ Qh on note pi ses omposantes sur la base (ψi ), don ph = i=1 pi ψi , soit don
PnQ
ph (x) = i=1 pi ψi (x) pour tout x ∈ Ω.
(n )
ei Q ) la base anonique de RnQ . Soit l'appli ation linéaire de hangement de base LQ :
Soit (~
(n )
Qh → RnQ dénie par LQ (ψi ) = ~ei Q : permet de passer de Qh (espa e de fon tions) à RnQ(espa e 
p1
n  . 
ve toriel usuel), et de faire les al uls dans R Q . Pour ph ∈ Qh on note LQ (ph ) = p̃h =  .. ,

pnQ
PnQ (nQ ) nQ
où don p̃h = i=1 pi~
ei ∈R .  
u1h
uh = k=1 ukh =noté
P2
Idem ave ~ ∈ R2 dans l'espa e géométrique R2 (dans R3 on ajoute
u2h
PnV k PnV k (nV )
u3h ), où k k k
 uhk ∈Vh (fon tion) ; pour uh = i=1 ui ψi ∈ Vh on note ũh = i=1 ui ~ei ∈ RnV , noté
u1  1
k  ..  ũh
ũh =  . . Et on pose ũh = ∈ R2nV .
ũ2h
uknV
Le problème (8.1) se met sous la forme matri ielle :

BhT f˜h
     
Ah ũh
. = . (8.7)
Bh 0 p̃h 0

18 21 juin 2018
19 8. Approximation numérique du problème de Stokes

8.4 Le problème du multipli ateur de Lagrange (la pression p)

Le al ul de ~uh ne pose pas de problème : ayant supposé Vh ⊂ V , la forme bilinéaire donnée
par a(~uh , ~vh ) est ontinue oer itive sur KerBh , et la matri e A asso iée est inversible et bien
onditionnée.
Le terme qui peut poser problème est le terme (ph , div~vh )L2 : la pression appro hée ph ne peut
être al ulée qu'à l'aide de e terme (le seul où apparaît ph dans (8.1)).
Notant ~ uh , grad~
(grad~ ~ vh )L2 − (f~, ~vh )L2 = (g, ~vh )L2 , le problème sous-ja ent à résoudre est :
trouver ph ∈ Qh t.q. :
(ph , div~vh )L2 = (g, ~vh )L2 , ∀~vh ∈ Vh , (8.8)

soit, f. (8.3) et/ou (8.4) :

(ph , divh~vh )L2 = (g, ~vh )L2 , ∀~vh ∈ Vh , (8.9)

Don :
1- Cas favorable : divh = ΠQh div : Vh → Qh est surje tif, ave de plus la onstante inf-sup
indépendante de h :
(div~vh , ph )L2 (Ω)
∃k > 0, inf sup ≥ k. (8.10)
~
vh ∈Vh ph ∈Qh ||~vh ||H01 ||ph ||L2
(Cas des éléments nis  P2 − P1  de TaylorHood en vitesse-pression.) Alors le problème (8.4),
i.e. (8.7), est inversible et bien onditionné, f. théorème 3.2.
~ uh , grad~
~ vh )L2 assure que la matri e ( Ah B t ) a ses lignes
Pour le problème (8.1), le terme
 (grad~ h
Ah BhT
 
Ah
indépendantes, et que la matri e a ses olonnes indépendantes, et la matri e
Bh Bh 0
est inversible et bien onditionnée i i.

2- En revan he, si divh ; Vh → Ph est surje tif :

(div~vh , ph )L2 (Ω)

∃kh > 0, inf sup ≥ kh , mais kh −→ 0, (8.11)
ph ∈Qh ~
vh ∈Vh ||~vh ||H01 ||ph ||L2 h→0

alors bien que le problème (8.4) soit mathématiquement inversible, il ne l'est plus numériquement
(mal onditionné, f. (3.9)2 , omme souvent pour les éléments nis  P1 − P1  ontinus en vitesse-
pression).

3- Ou pire divh : Vh → Qh n'est pas surje tif (et il n'existe pas de kh > 0). Exemple lassique
de l'instabilité en  her kerboard (é hiquier) donnée en 2-D par les éléments nis sur un maillage
de arrés ave Vh = Q1 - ontinus en vitesse, Qh = Q0 en pression.

8.5 Analyse à l'aide des proje tions : des pertes

Reprenons (8.6) : notant g(~vh ) = (f~, ~vh )L2 − (grad~uh , grad~vh )L2 on a :

~ h , ~vh )L2 = g(~vb ),

(ΠVhgradp ∀~vh ∈ Vh (8.12)

Quand un hoix parti ulier d'un Vh donne ~

ΠVhgrad(Q ~
h ) ( grad(Qh ), on a dégradé l'opérateur
gradient ~
grad : il est deveni l'opérateur ~ ~ =noté ΠV (grad)
gradient dis rétisé gradh = ΠVh ◦ grad ~ . Et
h
alors, omme :
~ h = ΠV gradp
~ h + gradp
~ h − ΠV gradp
~ h ,


gradp h h
(8.13)

~ h,
gradp ~
grad ~
relativement à l'information souhaité l'opérateur h = ΠVh(grad) fait perdre :

~ h − ΠV gradp
(gradp ~ h ). (8.14)
h

Et suivant les as, on se retrouve dans les situations du Ÿ pré édent : perte admissible ( f. as 1)
ou perte non admissible ( f. as 2- et 3-).

19 21 juin 2018
20 8. Approximation numérique du problème de Stokes

8.6 Une orre tion : ompensation des pertes

Pour ré upérer l'information (8.14) perdue (à ause de l'utilisation de la méthode des éléments
nis), on la réintroduit dans (8.1) : nouveau problème : trouver ~uh ∈ Vh et ph ∈ Qh t.q. :

~ uh , grad~
 (grad~
 ~ vh )L2 − (ph , div~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,
~ h −ΠV gradp ~ h , gradq~ h −ΠV gradq~ h )L2 (K) = 0, ∀qh ∈ Qh ,
X
 − (div~uh , qh )L2 −
 h2K (gradp h h
K
(8.15)
2
P
la somme K étant une somme sur les éléments K du maillage.
√ On a mis le oe ient hK (où
hK est la taille moyenne de la maille K , par exemple hh = aire de K de la maille en 2-D) pour
des problèmes de dimensionnement et de onvergen e optimale ( 'est ph qui nous intéresse et qui
est ontrlé par le terme ajouté, grâ e aux estimations de type  ||ph || ~ h || dans le as
≃ h ||gradp
de maillages réguliers). Voir Leborgne [22℄. Ainsi pour des éléments nis P1 ontinus en vitesse
et pression, on obtient des résultats de onvergen e à l'ordre 1 ( omme lassiquement pour des
éléments nis P1 ).

Remarque 8.1 Le problème initial était de trouver un point selle du lagrangien :

1
L(~vh , ph ) = ||grad~vh ||2L2 − (ph , div~vh )L2 − (f, vh )L2 . (8.16)
2
Et maintenant un point selle du lagrangien :

1 1X 2 ~ h −ΠV gradp
~ h ||2 2 .
Lmod (~vh , ph ) = ||grad~vh ||2L2 − (ph , div~vh )L2 − (f, vh )L2 − hK ||gradp h L (K)
2 2
K
(8.17)

Noter que l'équation (8.15) s'é rit plus simplement : trouver ~uh ∈ Vh et ph ∈ Q h t.q. :

~ uh , grad~
~ vh )L2 − (ph , div~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,
(
(grad~
(8.18)
~ h −ΠV gradp
− (div~uh , qh )L2 − h2 (gradp ~ h , gradq~ h )L2 = 0, ∀qh ∈ Qh ,
h

puisque ~ h −ΠV gradp

(gradp ~ h, w~h) = 0 pour tout w
~ h ∈ Vh par dénition de la proje tion ΠVh, en
h

parti ulier pour w ~ h.

~ h = ΠVhgradq
Et pour résoudre (8.15) il faut al uler une nouvelle in onnue ~ h ∈ Vh ,
~zh = ΠVhgradp i.e.
résoudre : trouver ~uh , ~zh ∈ Vh et ph ∈ Q h t.q. :


 ~ uh , grad~
(grad~ ~ vh )L2 − (ph , div~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,

~ h , gradq
− (div~uh , qh )L2 − h2 (gradp ~ h )L2 + h2 (~zh , gradq
~ h )L2 = 0, ∀qh ∈ Qh , (8.19)

 2 ~
h (gradph , ~zh′ )L2 − h2 (~zh , ~zh′ )L2 = 0, ∀~zh′ ∈ Vh .


Mais le al ul de ~zh
est presque gratuit ave des éléments nis P1 en vitesse et pression ar la
~zh′ peut être rendue diagonale à l'aide de la méthode
matri e (de masse) relative à l'équation en
de mass lumping (matri e de masse approximée à l'ordre 1 par une matri e diagonale). Et alors
la dernière équation donne expli itement ~zh en fon tion de ~ h
gradp (ave une pré ision d'ordre 1).

8.7 Méthode de Brezzi et Pitkäranta

Une méthode antérieure proposée par Brezzi et Pitkäranta [9℄ onsiste à pénaliser ( ontrle des
os illations de p en ontrolant le lapla ien de p) : trouver ~uh ∈ Vh et ph ∈ Q h t.q. :

~ uh , grad~
~ vh )L2 − (ph , div~vh )L2 = (f~, ~vh )L2 , ∀~vh ∈ Vh ,
(
(grad~
(8.20)
~ h , gradq
− (div~uh , qh )L2 − εh2 (gradp ~ h )L2 = 0, ∀qh ∈ Qh ,

où ε > 0 est petit. Cependant l'intégration par partie (interprétation du problème) donne la
∂p
ondition aux limites naturelle
∂n = 0, ondition qui n'est pas dans le problème à résoudre (7.5)
ou (8.1). Et la programmation donne ee tivement des résultats t.q. ~ h · ~n ≃ 0
gradp ( 'est la

20 21 juin 2018
21 8. Approximation numérique du problème de Stokes

ondition aux limites attendue pour le terme ~

(gradp, ~
gradq)L2 qui intervient dans la dis rétisation
du problème du lapla ien  (−∆p, q)L2 .)
Cette ondition aux limites de Neumann homogène non souhaitée (étrangère au problème ini-
tial) est atténuée dans le as de la programmation de (8.19).

Remarque 8.2 L'intérêt d'une orre tion de type (8.20) (ou (8.15)) peut se omprendre par
exemple en observant les os illations numériques non ontrlées de la pression (phénomène de type
he kerboard) : pour ontrler les variations de pressions, on ontrle son gradient par l'intermé-
diare de la norme ~
||gradq|| 2
L2 . Au lieu de onsidérer le lagragien :
1
L(~v , q) = ||grad~u||2 − (div~v , q) − (f~, ~v )L2 ,
2
dont la re her he d'un extremum donne (8.20), on onsidère le lagrangien pénalisé :

1 1 ~
Lε (~v , q) = ||grad~u||2 − (div~v , q) − f − ε||gradq|| 2
.
2 2
L'extremum de Lε est alors donné par :


 ~u ∈ V et p ∈ Q t.q. :
trouver

~ u, grad~
(grad~ ~ v )L2 − (p, div~v )L2 = (f~, ~v )L2 , ∀~v ∈ V, (8.21)

~ ~

 − (div~u, q) 2 − ε(gradp,
L gradq) L2 = 0, ∀q ∈ Q,
d'où la dis rétisation (8.20) ave ε rempla é par εh2 pour avoir un résultat de onvergen e optimale
en norme ||ph ||L2 (Ω) .

8.8 Méthode de Hughes, Fran a et Balestra

Cette méthode dite de Galerkin Least-squares [21℄ onsiste à stabiliser à l'aide de l'équation
satisfaite par la solution (méthode naturellement onsistante) : au lieu de (8.16), on onsidère :

1 εX 2 ~
L(~v , p) = ||grad~v ||2L2 (Ω) − (p, div~v )L2 (Ω) − (f, v)L2 (Ω) − h || − ∆u+gradp−f ||2L2 (K) , (8.22)
2 2
K
P
la somme K étant une somme sur les éléments K du maillage (même démar he que pré édem-
ment). Cette méthode né essite une valeur de ε susamment petite pour ne pas détruire l'ellipti ité
en u (ellipti ité du terme (gradu, ~ ~
gradv) L2 (Ω) − ε
P 2
k h (∆u, ∆v)L2 (K) ) : à l'aide de l'inégalité in-
verse
~ h ||L2 ,
||∆uh ||L2 ≤ Ch||gradu ∀uh ∈ Vh ,
1
on supposera don que 0<ε< √ . Et la méthode est programmée à l'aide de :
C

~
(gradu, ~
gradv) 2 ~ ~
L2 − (p, divv) − (q, divu) − εh (−∆u + gradp − f, −∆v + gradq)L2 = (f, v)L2 , (8.23)

pour tout v et tout q (don deux équations).

Remarque 8.3 Dans le as Vh = P1 , ette méthode s'apparente à elle de BrezziPitakäranta (il

n'y a plus de terme en ∆u et l'analyse de stabilité est la même).

8.9 Méthode de Douglas et Wang

Cette méthode [14℄ est in onditionnellement stable, mais non symétrique. Pour ne pas détruire
l'ellipti ité en u, l'équation (8.23) est modiée en, pour tout v et q :

~
(gradu, ~
gradv) 2 ~ ~
L2 − (p, divv) + (q, divu) + εh (−∆u + gradp − f, −∆v + gradq)L2 = (f, v)L2 . (8.24)

En parti ulier, en prenant v=u et q=p, la forme bilinéaire asso iée (membre de gau he) vérie :

~
(gradu, ~
gradu) 2 ~ ~
L2 − (p, divu) + (p, divu) + εh (−∆u + gradp, −∆u + gradp)L2
~
= ||gradu|| 2 2 ~
2 + εh || − ∆u + gradp||L2 ,
L

et la oer ivité en u ne dépend plus de ε>0. Cette méthode est adaptée aux modi ations des
équations de Stokes, lorsque le problème est de nature non symétrique.

Remarque 8.4 Même remarque que remarque 8.3.

21 21 juin 2018
22 9. Lapla ien

9 Lapla ien

9.1 Problème initial

Problème type : pour f ∈ H −1 (Ω), trouver p ∈ H01 (Ω) t.q. :

−∆p = f. (9.1)

Forme faible : trouver p ∈ H01 (Ω) t.q. :

~
(gradp, ~
gradq)L2 = hf, qi, ∀q ∈ H01 (Ω). (9.2)

Problème de minimum asso ié : ave J(q) =déf 12 ||gradq||

~ 2 2
L2 − hf, qi, trouver p ∈ H0 (Ω) t.q. J(p) =
minq∈H02 (Ω) J(q).

9.2 Problème sous ontrainte

Dans (9.1) on introduit (par exemple si on a besoin de al uler ~ )
gradp :

~
~u = gradp, et don − div~u = f. (9.3)

Et on réé rit (9.1) sous la forme : trouver (~u, p) ∈ L2 (Ω) × H01 (Ω) t.q. :

~ ~v )L2 = 0, ∀~v ∈ L2 (Ω),

(
(~u, ~v )L2 − (gradp,
(9.4)
~
− (~u, gradq) L2 = −hf, qiH −1 ,H 1 , ∀q ∈ H 1 (Ω).
0 0

Et p est devenu le multipli ateur de Lagrange de la ontrainte div~u + f = 0, f. IPP.

Interprétation de (9.4). u, p) ∈ L2 (Ω)

Si (~ ~
× H01 (Ω) est solution, alors ~u = gradp ∈ L2 (Ω)
n

−1
dans Ω, et div~u = −f dans H (Ω). Don ∆p = f dans H −1 (Ω) pour p ∈ H01 (Ω).
Notons :
n n
a(~u, ~v ) = (~u, ~v )L2 sur L2 (Ω) × L2 (Ω) ,
(
n (9.5)
~
b(~v , q) = −(~v , gradq) L2
2 1
sur L (Ω) × H (Ω). 0
2 n 1
Et (9.4) prend la forme (3.1) ave V = L (Ω) , Q = H0 (Ω).
2 n −1
I i B : ~
v ∈ L (Ω) → B~v = div~v ∈ H (Ω) est surje tif, f. (6.3).

Remarque 9.1 Quand f ∈ L2 (Ω), modions (9.4) en : trouver (~u, p) ∈ H div (Ω) × L2 (Ω) t.q. :

(
(~u, ~v )L2 + (p, div~v )L2 = 0, ∀~v ∈ H div (Ω),
(9.6)
(div~u, q)L2 = (f, q)L2 , ∀q ∈ L2 (Ω).

Interprétation de (9.6). Si(~u, p) ∈ H div (Ω) × L2 (Ω) est solution, alors div~u = f ∈ L2 (Ω),
~
~u = −gradp ∈ H −1 (Ω), don −∆p = f ∈ L2 (Ω), ave et ∂n ∂p 1
= 0 ( ar Imγn = H − 2 (Γ)). Don p
|Γ
∂p
est solution du problème de Neumann −∆p = f dans Ω et
∂n =0 sur Γ.
Notons : (
a(~u, ~v ) = (~u, ~v )L2 sur H div (Ω) × H div (Ω),
(9.7)
b(~v , q) = (div~v , q)L2 sur H div (Ω) × L2 (Ω).
div
Et (9.6) prend la forme (3.1) ave V = H (Ω), Q = L2 (Ω).
On a B : ~ v ∈ H (Ω) → B~v = div~v ∈ L2 (Ω) est surje tif, et a(·, ·) est H div - oer itif sur
div

KerB = {~v ∈ H div (Ω) : div~v = 0} = H0div (Ω). Don la formulation mixte (9.6) rentre dans le adre
lassique, f. théorème (3.2). ~
u est la variable prin ipale et p est le multipli ateur de Lagrange de
la ontrainte div~u = f .

22 21 juin 2018
23 10. Bilapla ien (problème biharmonique)

10 Bilapla ien (problème biharmonique)

10.1 Problème initial

Problème type : pour f ∈ H −2 (Ω) = (H02 (Ω))′ , trouver p ∈ H02 (Ω) t.q. :

∂p
∆(∆p) = f, où don p|Γ = 0 et = 0. (10.1)
∂n |Γ
Forme faible : trouver p ∈ H02 (Ω) t.q. :

(∆p, ∆q)L2 = hf, qiH −2 ,H02 , ∀q ∈ H02 (Ω). (10.2)

Interprétation : si p ∈ H02 (Ω) est solution alors deux IPP su essives donnent (10.1).
Le problème (10.2) est bien posé (LaxMilgram)
Problème de minimisation asso ié : ave J(q) =déf 12 ||∆q||2L2 − hf, qiH −2 ,H02 , trouver p ∈ H02 (Ω)
t.q. :
J(p) = min J(q). (10.3)
q∈H02 (Ω)

10.2 Formulation ave la variable ∆p

Une fon tion p t.q. ∆p ∈ L2 (Ω) est peu agréable à appro her pour la méthode des éléments
1
nis (éléments nis C d'Argyris par exemple). On se ramène à des gradients. Notons :

ϕ = ∆p (10.4)
(
ϕ = ∆p,
Le problème (10.1) est réé rit : trouver (ϕ, p) ∈ L2 (Ω) × H02 (Ω) t.q. :
∆ϕ = f.
Forme variationnelle proposée dans le as f ∈ H −1 (Ω) : trouver (ϕ, p) ∈ H 1 (Ω) × H01 (Ω) t.q. :

~ ~ ∀ψ ∈ H 1 (Ω),
(
(ϕ, ψ)L2 + (gradp, gradψ)L2 = 0,
(10.5)
~
(gradϕ, ~
gradq)L2 = −hf, qiH −1 ,H01 , ∀q ∈ H01 (Ω).
1 1
Interprétation : soit (ϕ, p) ∈ H (Ω) × H0 (Ω) une solution. (10.5)1 est en parti ulier vraie
1 −1
pour tout ψ ∈ H0 (Ω), don ϕ − ∆p = 0 dans H
R ∂p (Ω), don ∆p = ϕ ave ϕ ∈ H 1 (Ω), et
∂p
p ∈ H (Ω) au moins, et Γ ∂n ψ dΓ = 0 pour tout ψ ∈ H 1 (Ω) donne ∂n
2
|Γ
= 0. Et (10.5)2 donne
−1 2 1
∆ϕ = f ∈ H (Ω), d'où ∆ p = f , ave p ∈ H0 (Ω), d'où (10.1).
Notons :
a(ϕ, ψ) = (ϕ, ψ)L2 sur H 1 (Ω) × H 1 (Ω),
(
(10.6)
~
b(ϕ, v) = (gradϕ, ~
gradv)L2
1 1
sur H (Ω) × H (Ω). 0
Et (10.5) s'é rit sous la forme (3.1).
1
I i B : ϕ ∈ H (Ω) → Bϕ = −∆ϕ ∈ H −1 (Ω) est surje tif (LaxMilgram pour g ∈ H −1 (Ω)
~ ~ ~ 1 1
et (gradϕ, gradψ)L2 = (g, gradψ)L2 et ϕ ∈ H (Ω)/R, i.e. ϕ ∈ H (Ω) à une onstante près). Mais
1
KerB = {ψ ∈ H (Ω) : ∆ψ = 0} (l'ensemble des fon tions harmoniques), et a(·, ·) n'est pas (·, ·)H 1 -
~
oer itif dans KerB ( ar ||gradψ|| L2 n'est pas ontrlé : prendre par exemple ψ(x, y) = cos(nx)e
ny

sur Ω la boule unité, qui vérie bien ∆ψ = 0). Don (10.5) est mal posé au sens lassique du
théorème 3.2. Et une perte de pré ision, pour une approximation ave la méthode des éléments
nis, est attendue.

10.3 Formulation ave la variable ~

gradp
10.3.1 Formulation initiale
Dans (10.2) on introduit :

~
~u = gradp, don ∆p = div~u. (10.7)

Ayant p∈ H01 (Ω) dans (10.1), on notera pour les ~v 1

dérivant d'un potentiel H0 (Ω) :

~ 1 ~ noté
quand ~v ∈ grad(H 0 (Ω)) et ~v = gradq ave q ∈ H01 (Ω), on note ~v = ~vq . (10.8)

~
~u = ~vp = gradp ~
∈ grad(H 1
Le problème (10.2) devient : trouver 0 (Ω)) t.q. :

(div~u, div~vq )L2 = hf, qiH −2 ,H02 , ∀q ∈ H02 (Ω). (10.9)

23 21 juin 2018
24 10. Bilapla ien (problème biharmonique)

10.3.2 Première formulation sous forme ontrainte

~
grad(H 1
H01 (Ω)
Pour éviter de travailler ave le petit espa e 0 (Ω)), on se pla e dans tout en
ajoutant la ontrainte ~
~u − gradp = 0, f. (10.7), ontrainte mise en oeuvre ave un multipli ateur
de Lagrange ~λ.
Et pour ~
p ∈ H02 (Ω) solution de (10.2), posant ~u = gradp 1
u = ∆p ∈ L2 (Ω)
(dans H (Ω)) on a div~
∂p div
et ~u.~n|Γ = ∂n |Γ
= 0 , don ~
u ∈ H 0 (Ω). On note :

X = H0div (Ω) × H01 (Ω), (10.10)

muni du produit s alaire :

~
((~u, p), (~v , q))X = (div~u, div~v )L2 + (gradp, ~
gradq)L2 (Ω) , (10.11)

qui en fait un Hilbert.

n
Quand f ∈ H
−1
(Ω), on transforme alors (10.9) en : trouver ((~u, p), ~λ) ∈ X × L2 (Ω) t.q. :

(div~u, div~v )L2 + (~λ, ~v − gradq)

~
(
L2 = hf, qiH −1 ,H01 , ∀(~v , q) ∈ X,
(10.12)
~ n
(~u − gradp, µ)L2 = 0,
~ µ ∈ L2 (Ω) ,
∀~

soit :
(div~u, div~v )L2 + (~λ, ~v )L2 = 0, ∀~v ∈ H0div (Ω),




− (~λ, gradq)
~ L2 = hf, qiH −1 ,H01 , ∀q ∈ H01 (Ω), (10.13)

n
~ µ ∈ L2 (Ω) .

(~u, ~
µ)L2 − (gradp, µ)L2 = 0, ∀~
~


n
Interprétation : soit ((~u, p), ~λ) ∈ X × L2 (Ω) une solution. Dans Ω : ~λ = grad(div~
~ u), div~λ = f ,
~ ~
~u = gradp, don div~u = ∆p, et div(grad(∆p)) = f , i.e. ∆ p = f . Et sur Γ : ~u ∈ H0div (Ω) donne
2

~ n = 0, qui, ave p ∈ H01 (Ω), donne p ∈ H02 (Ω).

~u.~n = 0, don gradp.~
I i (10.12) orrespond à (3.1) ave :

(
a((~u, p), (~v , q)) = (div~u, div~v )L2 sur X × X,
n (10.14)
b((~v , q), ~ ~
µ) = (~v − gradq, µ)L2 sur X × L2 (Ω) .
~
n
H0div (Ω) × H01 (Ω) → L2 (Ω)
( )
Et on a B: . Et KerB = {(~v , q) ∈ H0div (Ω) × H01 (Ω) :
~
(~v , q) → B(~v , q) = ~v − gradq
~
~v = gradq} , don pour (~ v , q) ∈ KerB on a ∆p ∈ L2 (Ω) et a((~u, p), (~v , q)) = 12 (div~u, div~v )L2 +
1 2 1
2 (∆p, ∆q)L . Et pour p ∈ H (Ω) ∩ H0 (Ω) on dispose de ||∆p||L ≥ C ||p||H ≥ C ||p||H , f. (5.28).
2 2 2 1

Don a(·, ·) est oer itif sur (KerB, (·, ·)X ).

2 n′ n
Mais B n'est pas surje tif a priori (pour ℓ ∈ L (Ω) ≃ L2 (Ω) on doit trouver (~v , q) ∈
~
H0div (Ω) × H01 (Ω) t.q. ~v − gradq = ℓ, voir (5.30)). Et ~λ n'est pas ontrlé de manière lassique,
f. (3.9)2 . Pour la dis rétisation, on peut don s'attendre à une perte de pré ision.

10.3.3 Formulation bis sous forme ontrainte

Au lieu de X donné en (10.10), on onsidère :

n
X+ = H01 (Ω) × H01 (Ω), (10.15)

muni du produit s alaire produit :

~
((~u, p), (~v , q))X+ = (grad~u, grad~v )L2 + (gradp, ~
gradq)L2 (Ω) (10.16)

qui en fait un Hilbert.

24 21 juin 2018
25 10. Bilapla ien (problème biharmonique)

P R ∂2p ∂2q
On remarque que (dans une base eu lidienne) on a (∆p, ∆q)L2 = ij Ω ∂x2i ∂x2j dΩ, et, pour

p, q ∈ H02 (Ω) :

∂ 2 p ∂ 2q ∂ 3 p ∂q ∂2p ∂2q
Z Z Z
dΩ = − dΩ = dΩ. (10.17)
Ω ∂x2i ∂x2j Ω ∂x2i ∂xj ∂xj Ω ∂xi ∂xj ∂xi ∂xj

~
grad(H 1 ~ 1
Don , dans 0 (Ω)), pour tout ~vp , ~vq ∈ grad(H 0 (Ω)) f. (10.8), on a :

~
(div~vp , div~vq )L2 = (∆p, ∆q)L2 = (grad(gradp), ~
grad(gradq))L2 = (grad~
vp , grad~vq )L2 . (10.18)

~
~u = ~vp ∈ grad(H 1
Don (10.9) s'é rit aussi : trouver 0 (Ω)) t.q., ave (10.8) :

(grad~u, grad~vq )L2 = (f, q)L2 , ∀q ∈ H01 (Ω). (10.19)

n
Et quand f ∈ H −1 (Ω), on transforme (10.9) en : trouver ((~u, p), ~λ) ∈ X+ × L2 (Ω) t.q. :

(grad~u, grad~v )L2 + (~λ, ~v − gradq)

~
(
L2 = hf, qiH −1 ,H01 , ∀(~v , q) ∈ X+ ,
(10.20)
~ n
(~u − gradp, µ)L2 = 0,
~ µ ∈ L2 (Ω) ,
∀~

où ~λ est en ore le multipli ateur de Lagrange de la ontrainte (10.7). Soit :

n
(grad~u, grad~v )L2 + (~λ, ~v )L2 = 0, ∀~v ∈ H01 (Ω) ,




− (~λ, gradq)
~ L2 = hf, qiH −1 ,H01 , ∀q ∈ H01 (Ω), (10.21)

n
~ µ ∈ L2 (Ω) .

(~u, ~
µ)L2 − (gradp, µ)L2 = 0, ∀~
~


I i (10.20) orrespond à (3.1) ave :

(
a((~u, p), (~v , q)) = (grad~u, grad~v )L2 sur X+ × X+ ,
n (10.22)
b((~v , q), ~
µ) = (~v , ~ ~
µ)L2 − (gradp, µ)L2 sur X+ × L2 (Ω) ,
~
( 1 n n
H0 (Ω) × H01 (Ω) → L2 (Ω)
)
~ n
Et B : . Et (~
v , q) ∈ KerB ssi gradq = ~v ∈ H01 (Ω) ,
~
(~v , q) → B(~v , q) = ~v − gradq
don a(·, ·) est oer itif sur (KerB, (·, ·)X+ ). (La diéren e ave (10.14) est qu'i i on a un
ontrle ||.||H 1 pour ~
u, et pas uniquement un ontrle ||.||H div .)
Mais B n'est pas surje tif a priori. Et ~ λ n'est pas ontrlé de manière lassique.

10.3.4 Formulation ter sous forme ontrainte

Soit :
Y = H0div (Ω) × H 1 (Ω), (10.23)

muni du produit s alaire :

((~u, p), (~v , q))Y = (div~u, div~v )L2 + (p, q)H 1 (10.24)

qui en fait un Hilbert.

2
Quand f ∈ L (Ω) (ou dans (H 1 (Ω))′ ), on transforme (10.12) en : trouver ((~u, p), ~λ) ∈ Y ×
div
H (Ω) t.q. :

(div~u, div~v )L2 + (~λ, ~v )L2 + (q, div~λ)L2 = (f, q)L2 ,

(
∀(~v , q) ∈ Y,
(10.25)
µ)L2 + (~λ, div~
(~u, ~ µ) = 0, ∀~ µ ∈ H div (Ω),

i.e. t.q. : 
~ div
 (div~u, div~v )L + (λ, ~v )L = 0, ∀~v ∈ H0 (Ω),

 2 2

(div~λ, q)L2 = (f, q)L2 , ∀q ∈ H 1 (Ω), (10.26)


µ ∈ H div (Ω).

(~u, ~
µ)L2 + (p, div~
µ)L2 = 0, ∀~


Interprétation : dans Ω : ~λ = grad(div~

~ u), div~λ = f , ~u = gradp
~ , don ~
div(grad(div ~
gradp)) = f,
2
µ ∈ H div (Ω), don
R
i.e. ∆ p = f . Et sur Γ : p ~
µ .~
n dΓ = 0 = hp, µ
~ .~
n i 1 −1 pour tout ~
Γ H (Γ),H
2 2 (Γ)

25 21 juin 2018
26 10. Bilapla ien (problème biharmonique)

1 1
p|Γ = 0 dans H 2 (Γ) ( ar l'opérateur de tra e ~v ∈ H div (Ω) → ~v .~n ∈ H − 2 (Γ) est surje tif ). Don
~
p ∈ H01 (Ω). Ave ~u ∈ H0div (Ω), don gradp.~
n = ~u.~n = 0, don p ∈ H02 (Ω).
I i (10.25) orrespond à (3.1) ave :

(
a((~u, p), (~v , q)) = (div~u, div~v )L2 sur Y × Y,
(10.27)
b((~v , q), ~
µ) = (~v , ~
µ)L2 + (q, div~
µ)L2 sur Y × H div (Ω),
′
( )
H0div (Ω) × H 1 (Ω) → H div (Ω)
Et B : où hB(~v , q), µ
~ i = (~v , ~
µ)L2 + (q, div~
µ)L2 . Et B est sur-
(~v , q) → B(~v , q),
je tif, f. (5.24).
Et b((~v , q), ~
µ) = (~v , ~ ~
µ)L2 (Ω) − (gradq, µ)L2 (Ω) + (q, µ
~ ~ .~n)L2 (Γ) donne (~v , q) ∈ KerB ssi (~v , q) ∈
H0div (Ω) 1 ~
× H0 (Ω) ave ~v = gradq . Don a(·, ·) est oer itif sur (KerB, (·, ·)Y ).

Remarque 10.1 Ave Z = H div (Ω) × H01 (Ω), onsidérons (10.25) modié en : ((~u, p), ~λ) ∈ Z ×
H0div (Ω) t.q. : 
~ div
 (div~u, div~v )L + (λ, ~v )L = 0, ∀~v ∈ H (Ω),

 2 2

(div~λ, q)L2 = (f, q)L2 , ∀q ∈ H01 (Ω), (10.28)


µ ∈ H0div (Ω).

(~u, ~
µ)L2 + (p, div~
µ)L2 = 0, ∀~


Interprétation : dans Ω : ~λ = grad(div~

~ u), div~λ = f , ~u = gradp
~ , don div(grad(div
~ ~
gradp)) = f,
2
i.e. ∆ p = f . Et sur ~ ~
Γ : grad(div~u).~n = 0, don grad(∆p).~n = 0 ( ondition de Neumann), ave
p ∈ H01 (Ω).

10.3.5 Formulation quater sous forme ontrainte

Idem ave ( omme au Ÿ 10.3.3) :

Y+ = H01 (Ω) × H 1 (Ω), (10.29)

muni du produit s alaire :

~ u, grad~
((~u, p), (~v , q))Y+ = (grad~ ~ v )L2 + (p, q)H 1 . (10.30)

On rempla e (10.26) par :


~ div
 (grad~u, grad~v )L + (λ, ~v )L = 0, ∀~v ∈ H0 (Ω),

 2 2

(div~λ, q)L2 = (f, q)L2 , ∀q ∈ H 1 (Ω), (10.31)


µ ∈ H div (Ω).

(~u, ~
µ)L2 + (p, div~µ)L2 = 0, ∀~


Et on rempla e (10.27) par :

(
a((~u, p), (~v , q)) = (grad~u, grad~v )L2 sur Y+ × Y+ ,
(10.32)
b((~v , q), µ
~ ) = (~v , ~
µ)L2 + (q, div~
µ)L2 sur Y+ × H div (Ω)

26 21 juin 2018
27 11. Lo king

11 Lo king

Le lo king (blo age) intervient quand l'ellipti ité numérique explose (alors que l'ellipti ité du
problème initiale est nie) : dans e as, la solution numérique trouvée est pro he de la solution
identiquement nulle : la solution numérique est absurde.

11.1 Le problème type

Soit Ω un ouvert borné de Rn .
Soit λ ∈ R très grand
(mais xé).
On her he un hamp de ve teurs ~ u : Ω → Rn et un hamp de s alaires p : Ω → R qui
1 n 1
minimisent sur X = H0 (Ω) × H0 (Ω) la fon tionnelle :

1 λ ~ ~ v )L2 (Ω) − (g, q)L2 (Ω) .

M (~v , q) = ||grad~v ||2L2 (Ω) + ||~v − gradq|| L2 (Ω) − (f , ~ (11.1)
2 2
Remarque 11.1 ||grad~v ||2L2 (Ω) par a(~v , ~v ) où a(·, ·) est une forme bilinéaire
On peut rempla er
1
ontinue oer itive sur H0 (Ω), omme 'est le as pour les équations de MindlinReissner.

On a supposé Ω borné, et on dispose don de l'inégalité de Poin aré :

∃β > 0, ∀z ∈ H01 (Ω), ~

||gradz|| L2 (Ω) ≥ β||z||L2 (Ω) . (11.2)

Et H01 (Ω) est ~

alors muni de sa norme ||z||H 1 (Ω) = ||gradz||
1
L2 (pour tout z ∈ H0 (Ω)). De même
0
1 n 1 n 1
pour H0 (Ω) . On munit X = H0 (Ω) × H0 (Ω) du produit s alaire asso ié à la norme :

1
~
||(~v , q)||X = (||grad~v ||2L2 + ||gradq|| 2
L2 )
2

qui fait de X un Hilbert. La solution (~u, p) ∈ X réalisant le min de M vérie :

~ n
~v )L2 = (f~, ~v ),
(
(grad~u, grad~v )L2 + λ(~u − gradp, ∀~v ∈ H01 (Ω) ,
(11.3)
~
λ(~u − gradp, ~
gradq) L2 = (g, q), ∀q ∈ H 1 (Ω).
On va regarder e problème sous la forme LaxMigram en posant :

~
Φ((~u, p), (~v , q)) = (grad~u, grad~v )L2 + λ(~u − gradp, ~
~v − gradq)L2 , (11.4)

et en her hant (~u, p) ∈ X tel que pour tout (~v , q) ∈ X on ait (réé riture de (11.3)) :

Φ((~u, p), (~v , q)) = (f~, ~v )L2 + (g, q)L2 . (11.5)

Proposition 11.2 La forme bilinéaire Φ : X ×X → R est oer itive et ontinue sur X :

ave (11.2) :

1

 ∃αΦ > 0,
 ∀(~v , q) ∈ X, Φ((~v , q), (~v , q)) ≥ αΦ ||(~v , q)||2X , et αΦ ∼ ,
λ→∞ β
 ∃C > 0,
 ∀(~u, p), (~v , q) ∈ X, Φ((~u, p), (~v , q)) ≤ C||(~u, p)||X ||(~v , q)||X , et C = O(λ).
λ→∞
(11.6)
Et le problème (11.5) est bien posé.

Preuve. Bilinéarité. On a Φ symétrique : Φ((~u, p), (~v , q)) = Φ((~v , q), (~u, p)). Il sut don de mon-
trer que Φ((~u1 , p1 )+α(~u2 , p2 ), (~v , q)) = Φ((~u1 , p1 ), (~v , q))+αΦ((~u2 , p2 ), (~v , q)), e qui est immédiat.
Coer ivité. Pour tout κ > 0, on a ave (11.2) :
~
Φ((~v , q), (~v , q)) = ||grad~v ||2L2 + λ||~v − gradq|| 2
L2
~
≥ [(1−κ)||grad~v ||L2 + βκ||~v ||2L2 ] + λ[||~v ||2L2 + ||gradq|| 2
L2 − 2||~
v ||L2 ||q||L2 ].

On pose x = ||~v ||L2 et ~

y = ||gradq|| L2 , et on a :

λβκ 2
βκx2 + λ(x − y)2 ≥ y ,
λ + βκ
λβκ
la onstante c= λ+βκ étant la plus grande des onstantes vériant ette inégalité pour tout x, y > 0

27 21 juin 2018
28 11. Lo king

( al ul du c max vériant βκx2 + λ(x − y)2 ≥ cy 2 immédiat). On a don :

λβκ ~
Φ((~v , q), (~v , q)) ≥ (1−κ)||grad~v ||L2 + ||gradq|| L2 .
λ + βκ
λβκ 2
Le meilleur αΦ est obtenu en hoisissant κ 1−κ =
λ+βκ , i.e. pour κ solution de κ +
tel que

bκ − βλ = 0 où on a posé b = λ β+1 2 λ
β − 1. Le dis riminant vaut b + 4 β , et la ra ine positive vaut
q q
κ = 2 (−1 + 1 + 4 b2 β ). Pour λ >> 1, on a b ≃ λ, d'où 4 b2 β ≃ 4 λβ , d'où −1 + 1 + 4 b2λβ ≃ λβ
b λ λ 1 2
,
1
d'où κ ≃ β (quand λ >> 1, i.e. au voisinage de λ = ∞).
Continuité. Immédiat.

Remarque 11.3 Pour le problème numérique par éléments nis, quand

est grand on s'attend λ
C
à quelques di ultés. En eet, le onditionnement des matri es asso iées est d'ordre
αΦ = 0(λ),
et e onditionnement explose ave λ. Dans la pratique, même pour des valeurs de λ relativement
petites, le phénomène de lo king apparaît (voir plus loin la onstante de oer ivité numérique dans
la proposition 11.5), problème non lié au onditionnement, dû à la méthode des éléments nis.

11.2 Coer ivité nie en p

Pour l'analyse du lo king, regardons la oer ivité en p (terme où λ est présent) :

Proposition 11.4 Pour tout (~v , q) ∈ X on a, ave (11.2) :

λ ~ 2
Φ((~v , q), (~v , q)) ≥ β ||gradq|| L2 . (11.7)
λ+β
En parti ulier, quand λ est grand, on a :

λ
β ≃ β. (11.8)
λ+β

Preuve. On modie la démonstration pré édente : on a :

~
Φ((~v , q), (~v , q)) ≥ β||~v ||2L2 + λ[||~v ||2L2 + ||gradq|| 2 ~ 2
L2 − 2||~
v ||L2 ||q||L2 ] ≥ αp ||gradq|| L2 ,

2
où on her heαp le plus grand possible, i.e. le αp le plus grand possible dans l'inéquation  βx +
λ
λ(x − y) ≥ αp y 2 . On obtient αp = β λ+β
2
.

11.3 Le problème dis rétisé : onstante de oer ivité en ph non ontrlée

n n
Soit Vh ⊂ H01 (Ω) de dimension nie. On note ΠVh : L2 (Ω) → Vh la proje tion sur Vh , don :

n
∀~v ∈ H01 (Ω) , ∀w
~ h ∈ Vh , (ΠVh ~v , w
~ h )L2 (Ω) = (~v , w
~ h )L2 (Ω) .

Et soit Qh ⊂ H01 (Ω) de dimension nie, et on note Xh = Vh × Qh (don de dimension nie).

Le problème (11.5) dis rétisé s'é rit : trouver (~uh , ph ) ∈ Xh tel que :

Φ((~uh , ph ), (~vh , qh )) = (f~, ~vh ) + (g, qh ), ∀(~vh , qh ) ∈ Xh . (11.9)

Proposition 11.5 Pour tout (~vh , qh ) ∈ Xh on a :

λ ~ h ||2 2 + λ λ ||gradq
~ h − ΠV gradq
~ h ||2 2 ,
Φ((~vh , qh ), (~vh , qh )) ≥ β ||gradq L h L (Ω) (11.10)
λ+β λ+β

à omparer ave (11.7).

Illustration : siVh est susamment petit relativement à Qh  pour que dans ertaines dire tions
qh on ait à la fois ~ h 6= 0 et ΠV gradq
gradq ~ h = 0, alors la onstante de oer ivité en ph est de l'ordre
h
de λ, à omparer ave (11.8). Et don la solution dans Xh est majorée par une onstante de l'ordre
1
de ( λ ), i.e. s'appro he de la solution nulle quand λ est grand. C'est le phénomène de lo king.

28 21 juin 2018
29 11. Lo king

Preuve. Pour (~uh , ph ) et (~vh , qh ) ∈ Xh , on a :

~ h − ~uh , gradq
Φ((~uh , ph ), (~vh , qh )) = (grad~uh , grad~vh )L2 + λ(gradp ~ h − ~vh )L2
~ h − ΠV gradp
= (grad~uh , grad~vh )L2 + λ(gradp ~ h , gradq
~ h − ΠV gradq
~ h )L2
h h

+ λ(ΠV gradp~ h − ~uh , ΠV gradq

~ h − ~vh )L2 .
h h

D'où :

λ ~ h ||2 2 + λ||gradq
~ h − ΠV gradq
~ h ||2 2 ,
Φ((~vh , qh ), (~vh , qh )) ≥ β ||ΠVh gradq L h L
λ+β
λ λ
al ul fait au Ÿ pré édent. Et λ = β λ+β + λ λ+β et Pythagore donnent (11.5).

Remarque 11.6 Pour mémoire, on rappelle que le terme ~ h −ΠV gradq

gradq h
~ h est le terme qui pose
problème dans la stabilisation de la ondition inf-sup dans les équations de Stokes, voir Ÿ 8.6.

11.4 Une orre tion optimale du lo king

Comme on a éventuellement trop de oer ivité, due à l'approximation numérique hoisie (élé-
ments nis), f. (11.5), on l'enlève : on onsidère le problème (11.1) modié : trouver le minimum
(~u, p) ∈ Vh × Qh dans Vh × Qh de :

1 λ ~ h ||2 2 − λ λ ||gradq
~ h − ΠV gradq
~ h ||2 2 )
Mh (~v , q) = ||grad~vh ||2L2 + (||~vh − gradq L h L
2 2 λ+β (11.11)

− (f~, ~vh )L2 − (g, qh )L2 .

On a don modié (11.4) en :

~ h , ~v − gradq
Φh ((~uh , ph ), (~vh , qh )) = (grad~uh , grad~vh )L2 + λ(~uh − gradp ~ h )L2
λ2 ~ h − ΠV gradp
~ h , gradq
~ h − ΠV gradq
~ h )L2 ,
− (gradp h h
λ+β

et on résout : trouver (~uh , ph ∈ Vh × Qh t.q. pour tout (~vh , qh ∈ Vh × Qh :

Φh ((~uh , ph ), (~vh , qh )) = (f~, ~vh )L2 + (g, qh )L2 .

Il va falloir al uler ~ h,
ΠVh gradp mais si on hoisit Vh = P1 ( ontinus), e al ul se réduit à
inverser une matri e diagonale (mass-lumping), et n'est don pas oûteux.
Programmation : il s'agit de trouver (~uh , ph ) ∈ Vh × Qh t.q., pour tout (~vh , qh ) ∈ Vh × Qh :

~ h , ~vh )L2 = (f~, ~vh )L2 ,


 (grad~uh , grad~vh )L2 + λ(~uh − gradp
 −λ(~uh − gradp ~ h )L2 − λ λ (gradp
~ h , gradq ~ h −ΠV gradp ~ h , gradq
~ h )L2 = (g, qh )L2 .
h
λ+β

Et don , ave la variable auxiliaire w ~h ~ h

= ΠVh gradp : trouver (~uh , ph , w
~ h ) ∈ Vh × Qh × Vh t.q.,
pour tout (~ ~ h′ ) ∈ Vh × Qh × Vh :
vh , qh , w

~ h , ~vh )L2 = (f~, ~vh )L2 ,



 (grad~uh , grad~vh )L2 + λ(~uh , ~vh )L2 − λ(gradp

~ h )L2 + βλ (gradp ~ h )L2 + λ λ (w

−λ(~uh , gradq ~ h , gradq ~ h )L2 = (g, qh )L2 ,
~ h , gradq (11.12)

 λ+β λ+β

 ~ h, w
(gradp ~ h′ )L2 − (w ~ h′ )L2 = 0.
~h, w

(Ave des éléments nis Vh = P1 , la dernière équation s'élimine à l'aide de mass-lumping.)

Cette méthode donne des résultats optimaux de onverge. En parti ulier, pour des éléments
nis P1 en ~u et q, la onvergen e est en O(h).

29 21 juin 2018
30 11. Lo king

11.5 Traitement lassique du lo king : problème ontinu

11.5.1 Problème ontinu
Voir par exemple Chapelle [10℄, Brezzi et Fortin [8℄.
On introduit la variable :
~
~γ = λ(~u − gradp). (11.13)
n
Le problème (11.3) devient : trouver (~u, p, ~γ ) ∈ H01 (Ω) × H01 (Ω) × Y t.q. :

 a((~u, p), (~v , q)) + b((~v , q), ~γ ) = (f~, ~v )L2 + (g, q)L2 ,

∀(~v , q) ∈ X,
(11.14)
 b((~u, p), ~δ) − 1 (~γ , ~δ)L2 = 0, ∀~δ ∈ Y,
λ
où on a posé :

~ −1 n ~ −1

 Y = {δ ∈ (H (Ω)) : divδ ∈ H (Ω)},

a(·, ·) : X × X → R : a((~u, p), (~v , q)) = (grad~u, grad~v )L2 , (11.15)

b((~u, p), ~δ) = h~u − gradp,
~ ~δi = h~u, ~δi + hp, div~δi,

b(·, ·) : X × Y → R :

(apparen e d'un problème pénalisé) et où on montre que Y est un Bana h lorsqu'il est muni de la
norme :
déf
||~δ||Y = ||~δ||H −1 (Ω)n + ||div~δ||H −1 (Ω) .

Remarque 11.7 Vu l'espa e Y i-dessus, l'expression (~γ , ~δ)L2 n'a pas de sens dans (11.14)2 : il
aurait fallu é rire, au lieu de (11.14) :

 a((~u, p), (~v , q)) + b((~v , q), ~γ ) = (f~, ~v )L2 + (g, q)L2 ,

∀(~v , q) ∈ X,
~
 ~u − gradp 1
− ~γ = 0.
λ
1
L'espa e Y est l'espa e qui orrespond au as
λ =0 (i.e. λ inniment grand), f. le modèle de
plaque de Kir hhoLove :

a((~u, p), (~v , q)) + b((~v , q), ~γ ) = (f~, ~v )L2 + (g, q)L2 ,
(
∀(~v , q) ∈ X,
(11.16)
b((~u, p), ~δ) = 0, ∀~δ ∈ Y,

De plus, pour le problème dis rétrisé par éléments nis, on hoisit en général Yh ⊂ L2 (Ω), et (11.14)
a un sens dans Vh × Qh × Yh , d'où l'abus de notation dans (11.14).

Si B :X → Y′ est l'opérateur asso ié à b(·, ·), i.e. hB(~v , q), ~δi =déf b((~v , q), ~δ), i.e. :

~
B(~v , q) = ~v − gradq

alors il est simple de montrer que a(·, ·) est oer itif sur ~
KerB = {(~v , q) ∈ X : ~v = gradq} (à l'aide
de l'inégalité de Poin aré).
Puis on montre que b(·, ·) satisfait la ondition inf-sup :

b((~v , q), ~δ)

∃k > 0, ∀~δ ∈ Y, sup ≥ k.
(~
v ,q)∈X ||(~v , q)||X ||~δ||Y

Voir par exemple Chapelle [10℄, Brezzi et Fortin [8℄.

11.5.2 Problème dis rétisé

Pour le problème dis rétisé, la ondition inf-sup dis rète doit être satisfaite. C'est le traitement
de ette ondition dis rète qui a donné naissan e à de nombreux arti les.
Noter qu'i i il y a deux di ultés :
1- la stabilisation de ~γh , ou le hoix adequat d'éléments nis ompatibles pour satisfaire la
ondition inf-sup,

30 21 juin 2018
31 11. Lo king

2- la oer ivité sur le noyau KerBh où Bh est l'opérateur dis rétisé. Ce problème peut être
régler simplement en modiant (11.1) en :

1 1 ~ λ−1 ~ ~ v )L2 (Ω) − (g, q)L2 (Ω) ,

M̃ (~v , q) = ||grad~v ||2L2 (Ω) + ||~v − gradq|| L2 (Ω) + ||~v − gradq|| L2 (Ω) − (f , ~
2 2 2
(11.17)
et en remplaçant a(·, ·), f.-(11.15), en :

~
ã((~u, p), (~v , q)) = (grad~u, grad~v )L2 + (~u − gradp, ~
~v − gradq)L2 .

n
Et on onsidère don le problème : trouver (~u, p, ~γ ) ∈ H01 (Ω) × H01 (Ω) × Y t.q. :

 ã((~u, p), (~v , q)) + b((~v , q), ~γ ) = (f~, ~v )L2 + (g, q)L2 ,

∀(~v , q) ∈ X,
(11.18)
 b((~u, p), ~δ) − 1 (~γ , ~δ)L2 = 0, ∀~δ ∈ Y.
λ−1

Maintenant ã(·, ·) oer itif sur X (grâ e à (11.2)), et on a ~

KerB = {(~v , q) : ~v = gradq} : on s'est
ramené au as onnu des équations de Stokes : hoix des éléments nis, ou hoix des stabilisations.

31 21 juin 2018
32 12. Conditions aux limites de Diri hlet faible

12 Conditions aux limites de Diri hlet faible

12.1 Problème initial

1
On onsidère le problème, pour f ∈ H −1 (Ω) et d ∈ H 2 (Γ) donnés :
 1
 trouver u ∈ H (Ω) t.q. :

− ∆u + u = f dans Ω, (12.1)

u|Γ = d sur Γ.


Cher hant u ∈ d + H01 (Ω)(espa e ane d'espa e ve toriel asso ié H01 (Ω)), la démar he lassique
1
onsiste à é rire la formulation variationnelle dans H0 (Ω) :
(
trouver u ∈ d + H01 (Ω) t.q. :
(12.2)
(u, v)H 1 (Ω) = hf, viH −1 (Ω),H 1 (Ω) , ∀v ∈ H01 (Ω).

(Voir poly opié éléments nis pour l'appli ation de LaxMilgram : on prend un relèvement
1
zd ∈ H 1 (Ω) 1 1
de d ∈ H 2 (Γ), on pose u0 = u − zd ∈ H0 (Ω), et LaxMilgram dans H0 (Ω) donne un
problème bien posé en u0 = u0 (zd ). Et u = u0 + zd = u0 (zd ) + zd ne dépend pas du relèvement zd
hoisi.)

12.2 Problème ontraint

(Voir Babu²ka [2℄ pour l'arti le initial.)
I i on souhaite imposer la ontrainte  ondition aux limites de Diri hlet u|Γ − d = 0 à l'aide
1
d'un multipli ateur de Lagrange λ : trouver (u, λ) ∈ H 1 (Ω) × H − 2 (Γ) t.q. :

∀v ∈ H 1 (Ω),

 (u, v)H 1 (Ω) + hλ, viH − 12 (Γ),H 12 (Γ) = (f, v)L2 (Ω) ,
1 (12.3)
 hu, µi 1 1 = hd, µi 1 1 , ∀µ ∈ H − 2 (Γ),
H 2 (Γ),H − 2 (Γ) H 2 (Γ),H − 2 (Γ)

orrespondant au problème de re her he du point selle du Lagrangien :

1
L(u, λ) = ||u||2H 1 (Ω) + hu − d, λ) 12 1 − hf, viH −1 (Ω),H 1 (Ω) , (12.4)
2 H (Γ),H − 2 (Γ)

Quand ela a un sens, par exemple quand f ∈ L2 (Ω) pour avoir u ∈ H 2 (Ω), l'interprétation de e
problème donne (par intégration par parties) :

− ∆u + u = f dans L2 (Ω),




 1
u = d dans H 2 (Γ), (12.5)

 λ = − ∂u dans H 12 (Γ).



∂n
L'interprétation du multipli ateur de Lagrange apparaît : 'est la for e (au signe près) qu'il faut
appliquer sur le bord Γ pour que la position u vaille bien d.
On pose :
(
a(u, v) = (u, v)H 1 (Ω) ,
(12.6)
b(v, λ) = hv, λi 21 1 ,
H (Γ),H − 2 (Γ)
( 1
)
H 1 (Ω) → H 2 (Γ)
et (12.3) est de la forme (3.1), ave B : surje tif et KerB = H01 (Ω). On
v → Bv = γ0 (v)
dispose don de la ondition inf-sup :

1
∃k > 0, ∀λ ∈ H − 2 (Γ), ∃v ∈ H 1 (Ω) : |b(v, λ)| ≥ k||v||H 1 (Ω) ||λ|| 1 , (12.7)
H − 2 (Γ)

soit :
v λ
∃k > 0, inf1 sup |b( , )| ≥ k, (12.8)
λ∈H − 2 (Γ) v∈H 1 (Ω) ||v||H 1 (Ω) || λ|| 1
H − 2 (Γ)

Et la forme bilinéaire a(·, ·) est ontinue et oer itive sur H 1 (Ω) ( 'est le produit s alaire), don
1
sur KerB (= H0 (Ω)).

32 21 juin 2018
33 12. Conditions aux limites de Diri hlet faible

Remarque 12.1 Le al ul de λ en lui-même peut être important, mais le ontrle de λ est i i

− 12 2
uniquement un ontrle H (Γ) (même pas un ontrle L (Γ)). On peut don s'attendre à des
problèmes numériques.

Remarque 12.2 Cette appro he onduit aux problèmes de transmission aux parois traités par
la méthode des éléments nis ave joints (mortar nite elements) de Bernardi, Maday et Patera.

12.3 Problème appro hé

1
On se donne des espa es éléments nis Vh ⊂ H 1 (Ω) et Λh ⊂ H − 2 (Γ). Et le problème appro-
hant (12.3) à résoudre est :


 trouver (uh , λh ) ∈ Vh × Λh t.q. :

(uh , vh )H 1 (Ω) + (vh , λh )L2 (Γ) = (f, vh )L2 (Ω) , ∀vh ∈ Vh , (12.9)

(uh , µh )L2 (Γ) = (d, µh )L2 (Γ) , ∀µh ∈ Λh .


On obtient un problème numérique à résoudre sous la forme (8.7) :

Bt f~
     
A ~u
. ~ = . (12.10)
B 0 λ ~g

12.4 Eléments nis Pk −C 0 : non stables

12.4.1 Condition inf-sup dis rète
L'appro he élémentaire, pour la méthode des éléments nis, onsiste à prendre pour l'espa e
d'éléments nis de Lagrange Vh = Pk −C 0 dans
( Ω (don né essairement k ≥ 1), et) Λh = γ0 (Vh )
(Vh , ||.||H 1 (Ω) ) → (Λh , ||.|| 1 )
H − 2 (Γ)
Pk −C 0 sur Γ. Don l'opérateur Tr = γ0|Γ : est surje tif
vh → γ0 (vh ) = vh|Γ
(immédiat pour des éléments nis de Lagrange), et don :

|(λh , vh )L2 (Γ) |

∃kh > 0, ∀λh ∈ Λh , sup ≥ kh ||λh || − 12 . (12.11)
vh ∈Vh ||vh ||H 1 (Ω) H (Γ)

Soit :

∃kh > 0, ∀λh ∈ Λh , ∃vh ∈ Vh : |(λh , vh )L2 (Γ) | ≥ kh ||λh || 1 ||vh ||H 1 (Ω) . (12.12)
H − 2 (Γ)

1
Mais le ontrle de λ est i i uniquement un ontrle H − 2 (Γ), et kh dépend de h a priori.

12.4.2 Barbosa et Hughes

(Voir Barbosa et Hughes [3℄, Pitkäranta [26℄, Stenberg [30℄).
On veut utiliser omme maillage sur Γ la tra e du maillage Th déni dans Ω.
L'idée de Barbosa et Hughes est de stabiliser le multipli ateur de Lagrange λ à l'aide de sa valeur
∂u
souhaitée λ = −
∂n sur Γ. Le lagrangien formel initial (12.4) est transformé en ( f. Stenberg [30℄),
dans Vh × Λh :
∂u 2
Lh (u, λ) = L(u, λ) − α h ||λ+ || 2 , (12.13)
∂n L (Γ)
pour α une onstante à déterminer. I.e. le problème appro hé (12.9) est transformé en le problème :

 trouver (uh , λh ) ∈ Vh × Λh t.q. pour tout (vh , µh ) ∈ Vh × Λh :


∂uh ∂vh

 Z Z
(uh , vh )H 1 (Ω) + vh λh dΓ − αh (λh + ) dΓ = (f, vh )L2 (Ω) ,

∂n ∂n (12.14)
 Γ Γ
∂uh

 Z Z Z
uh µh dΓ − αh (λh + )µh dΓ = dµh dΓ.


∂n

Γ Γ Γ

Et le problème est ainsi transformé en le problème matri iel sous forme pénalisée :

Bt f~
     
A ~u
. = . (12.15)
B −αC p ~g

33 21 juin 2018
34 12. Conditions aux limites de Diri hlet faible

Théorème 12.3 (Barbosa et Hughes [3℄, Pitkäranta [26℄.) On suppose que le maillage Th est
quasi-uniforme, ou plus pré isement qu'il est tel que l'inégalité inverse suivante est satisfaite :

1 ∂vh ~ h ||L2 (Ω) .

∃Ci > 0, ∀vh ∈ Vh , h 2 || ||L2 (Γ) ≤ Ci ||gradv (12.16)
∂n
On suppose dans (12.13) que α est assez petit (pour ne pas détruire la oer ivité en u), à savoir :

1
0<α< . (12.17)
Ci
Alors le problème appro hé stabilisé (12.14) est bien posé, et de plus on a l'estimation d'erreur
optimale, dès que la solution exa teu ∈ H k+1 (Ω) :

||u − uh ||H 1 (Ω) ≤ Chk ||u||H k+1 (Ω) ,

où C est une onstante indépendante de h.

Preuve. Voir Barbosa et Hughes [3℄ et Stenberg [30℄.

Remarque 12.4 L'inégalité (12.16) se réé rit :

Z Z
h ~ n)2 dΓ ≤ Ci2
(gradv.~ ~
||gradv|| 2
Rn dΩ,
Γ Ω
R R
et le fa teur h est attendu : (h Γ ) a même dimension que ( Ω ), dimension d'un volume.

12.4.3 Elimination du multipli ateur : méthode de Nits he [25℄

Stenberg [30℄ a montré que la Méthode de Barbosa et Hughes [3℄ était équivalente à la méthode
plus an ienne de Nits he [25℄, dans le as où Λh = P0 (Γh ) et Vh = P1 (Ωh ) (quand Γh est la tra e
de Ωh ) :



 trouver uh ∈ Vh ∀vh ∈ Vh :
t.q.

∂uh ∂vh X 1
 (uh , vh )H 1 (Ω) − h ∂n , vh iΓ − h ∂n , uh − diΓ + γ

hE
huh − d, vh iE = (f, vh )L2 (Ω) ,
E∈Eh

i.e.
trouver uh ∈ Vh t.q. ∀vh ∈ Vh :



∂u ∂v X 1

 (uh , vh )H 1 (Ω) − h h , vh iΓ − h h , uh iΓ + γ

huh , vh iE


∂n ∂n hE (12.18)
E∈Eh


 ∂vh X 1
= (f, vh )L2 (Ω) − h , diΓ + γ hd, vh iE .


∂n hE


E∈Eh

On retrouve alors uΓ = d. Cette méthode est plus simple à programmer ar sans multipli ateur de
Lagrange.

Proposition 12.5 Sous l'hypothèse (12.16), quand γ > Ci , le problème (12.18) est bien posé :

||u − uh ||H 1 (Ω) ≤ Chk ||u||H k+1 (Ω) ,

pour une approximation Pk quand u ∈ H k+1 (Ω).

Preuve. Voir Stenberg [30℄.

34 21 juin 2018
35 12. Conditions aux limites de Diri hlet faible

Comparaison de la méthode de Nits he et de la méthode de Barbosa et Hughes : on

her he à éliminer λh dans (12.14). A partir de (12.14)2 on obtient :

∂uh 1
λh = −ΠΛh ( )+ (uh −d).
∂n αh
Reporté dans (12.14)1 , on obtient :

∂uh ∂vh 1
Z Z Z
(uh , vh )H 1 (Ω) − ΠΛh ( )vh dΓ − ΠΛh ( )uh dΓ + uh vh dΓ
Γ ∂n Γ ∂n αh Γ
∂uh ∂vh ∂uh ∂vh
Z
− αh( − ΠΛh ΠΛh dΓ)
Γ ∂n ∂n ∂n ∂n
∂vh 1
Z Z
= (f, vh )L2 (Ω) − g ΠΛh dΓ + gvh dΓ,
Γ ∂n αh Γ

On se pla e alors dans le as Λh = P0 et Vh = P1 : on a alors ΠΛh ( ∂u

∂n ) =
h ∂uh
∂n . Et on re-
trouve (12.18).

35 21 juin 2018
36 A. Singular value de omposition (SVD)

A Singular value de omposition (SVD)

Il s'agit d'estimer la onstante k (plus exa tement kh ) de la ondition inf-sup (numérique) : elle
donne une idée du onditionnement.

Rappel. Dans le as d'une matri e A arrée symétrique (réelle) n ∗ n, elle est diagonalisable,
ave D = P −1 .A.P = diag(λ1 , ..., λn ). En eet, é rire A.P = P.D revient à é rire que les olonnes
sont égales, i.e. que A.~vi = λi~vi où les ~vi sont les ve teurs (propres) olonnes de P . Et on peut
hoisir n ve teurs propres ~vi onstituant une base orthonormée, et ainsi P −1 = P T , i.e. P est
matri e orthogonale n ∗ n.
On s'intéresse i i à des matri es B re tangulaires m ∗ n et à leurs valeurs propres, appelés
valeurs singulières. On veut trouver une matri e Σ diagonale m∗n, au sens que tous les éléments
sont nuls sauf (éventuellement) les éléments Σii pour i = 1, ..., p = min(m, n), et deux matri es
U m∗m et V n∗n telles que :
Σ = U T .B.V.
On notera Σ = diagm,n (σ1 , ..., σp ) ette matri e m∗n diagonale.
Exemple, si m < n, on a p = m et

σ1 0 ... 0 0 ... 0
 
. .
. .
 0 σ1 0 0 . .
 

Σ = diagm,n (σ1 , ...σp ) =  .. .. ..
 
 . . .


 0 
0 ... σp 0 ... 0

et si m > n, on a p=n et diagm,n (σ1 , ...σp ) est la transposée de la matri e i-dessus.

On va se servir du :

Lemme A.1 Si B m ∗ n, et si λi est valeur propre de la matri e n ∗ n B T .B ,

est une matri e réelle
T
alors λi est valeur propre de la matri e m ∗ m B.B .
T
Et si λi est valeur propre de la matri e m ∗ m B.B , alors λi est valeur propre de la matri e
T
n ∗ n B .B .

Preuve. T
La matri e B .B (taille n ∗ n) est symétrique réelle, don diagonalisable. De plus elle
x .(B .B).~x = (B.~x)T .(B.~x) = ||B.~x||2 ≥ 0. On note λ1 , ..., λn ses valeurs propres
est positive ar ~
T T

(qui sont don ≥ 0),

et quitte à les renuméroter, on suppose λ1 ≥ ... ≥ λn (≥ 0). Et on note ~ vi
des ve teurs propres asso iés formant une b.o.n. de Rn , et V la matri e n ∗ n orthogonale dont les
olonnes sont les ve teurs propres ~vi . Comme Rang(B) ≤ p = min(m, n), on a Rang(B T .B) ≤ p,
et on a λp+1 = ... = λn = 0, d'où :

diagn,n (λ1 , ..., λp , 0, ..., 0) = V T .B T .B.V matri e n ∗ n.

Même démar he pour la matri e B.B T de valeurs propres µ1 ≥ .. ≥ µm ≥ 0 asso iée aux
ve teurs propres orthonormés ~u1 , ..., ~um formant la matri e orthogonale U : on a

diagm,m (µ1 , ..., µp , 0, ..., 0) = U T .B.B T .U matri e m ∗ m.

T
En parti ulier B.B .~ui = µi ~ui , d'où B T .B.B T .~ui = µi B T .~ui , d'où B T .~ui est valeur propre
T T
de B .B asso ié à la valeur propre µi . Et omme B .B.~ vj = λj ~vj on en déduit que µi est l'un
des λj .

Corollaire A.2 Si B est une matri e réelle m ∗ n, toujours ave p = min(m, n), on a la dé ompo-
sition en valeurs singulières de B :

Σ = U T .B.V, i.e. B = U.ΣT .V T , (A.1)

où Σ = diagm,n (σ1 , ..., σp ) ave les σi ≥ 0, où les olonnes de U sont les ve teurs propres de B.B T
T 2 T T
et elles de V de B .B , et les λi = σi sont les valeurs propres de B .B ou de B.B .

36 21 juin 2018
37 A. Singular value de omposition (SVD)

Don , si Rang(B) = r, et quitte à renuméroter supposant σ1 ≥ ... ≥ σr > 0 (et don σi = 0

pour i > r) :
r
X
B= σi ~ui .~viT . (A.2)
i=1

 De plus, si les olonnes de

 U et de V sont les ve teurs notés ~uj ∈ Rm et ~vj ∈ Rn , alors
~uj 0 B
∈ Rm+n pour j = 1, ..., p est ve teur propre de la matri e pour la valeur
~vj BT 0
propre σj .
Et si B est une matri e arré n ∗ n symétrique positive (réelle), alors les σi sont ses valeurs
propres, et U = V .

Preuve. Remarque : notons que si la dé omposition Σ = diagm,n (σ1 , ..., σp ) = U T .B.V existe, alors
on a Σ.ΣT = diagm,m (σ12 , ..., σp2 ) = (U T .B.V ).(U T .B.V )T = U T .(B.B T ).U et don les λi = σi2 sont
T
bien valeurs propres de B.B asso iées aux ve teurs propres les olonnes de U (par identivation
T 2 2 T
des olonnes dans l'égalité (B.B ).U = U.diagm,m (σ1 , ..., σp ) qui donne (B.B ).~ ui = σi2 ~ui ).
T
Idem pour B .B . Et les matri es U et V ontiennent bien les ve teurs propres annon és.
Montrons don que la dé omposition existe. On note λi pour i = 1, ..., n les valeurs propres de
T
la matri e B .B qui est symétrique positive. Quitte à renuméroter, on suppose λ1 ≥ ... ≥ λr > 0,
λr+1 = ... = λn = 0. On pose V = la matri e dont les olonnes sont les ve teurs propres ~ vj ∈ Rn
n
p
asso iés aux λj formant une b.o.n de R . On pose σj = λj .
B.~vj
Puis on pose ~ uj = σj ∈ R pour 1 ≤ j ≤ r. Les ~uj sont des ve teurs propres de B.B T 2 à
m

B.(B T .B).~
vj λj B.~vj viT .(B T .B~
~ vj )
2 orthonormaux. En eet, (B.B T ).~uj = σj = σj = λj ~uj . Et ~uTi .~uj = σi σj =
~v T .~
v
λj σii σjj = δij (et don les ve teurs ~uj sont les ve teurs B.~vj normalisés).
Puis on omplète la famille ~ uj ainsi obtenue pour obtenir une b.o.n. (~uj ) de Rm (par exemple
à l'aide de pro édé de GramS hmidt). On note U la matri e dont les olonnes sont les (~ uj ).
T
uTi .B.~vj ] = [σj ~uTi .~uj ] = [σj δij ] pour j ≤ r, et est nul si j
Soit Σ = U .B.V : on a [Σij ] = [~ >r
( as où B.~vj = 0). D'où le résultat (A.1).
On en déduit immédiatement (A.2).
Par ailleurs par onstru tion ~ uj est donnée par B~vj = σj ~uj pour j ≤ r et don B T .B.~vj =
σj B T .~uj = λj ~vj pour j ≤ r, d'où B T .~uj = σj ~vj 
pour j ≤ r. Et pour j >
    r on a B T .~uj = 0 ar un
0 B ~uj ~uj
tel ~uj ∈ (ImB)⊥ = Ker(B T ). Et don . = σj omme annon é.
BT 0 ~vj ~vj
T T 2 T
Enn, si B est une matri e symétrique (réelle) positive, alors B .B = B = B.B
et B .B
T
et B.B ont la même dé omposition. Et omme B est positive, ses valeurs propres sont positives,
√
σi = + λi , et sont ses valeurs singulières.

Remarque A.3 Si m > n, les ve teurs ~uj pour j ≥ m + 1 ne servent à rien, et on al ule
uniquement U omme matri e m∗ n (ave don U T matri e n∗ m). Et Σ est alors une matri e n∗ n.
La méthode de al ul ne al ulant sto kant que les ~
ui pour j = 1, ..., p (don pas d'algorithme de
Grams hmidt pour ompléter la base) est appelée la méthode Thin SVD. Cette méthode Thin
n'a de sens que si p = min(m, n) = n < m.

Corollaire A.4 On a Rang(B) = r, KerB = Vect{~vr+1 , ...~vn } et ImB = Vect{~u1 , ...~ur }.

Preuve. À l'aide de (A.2).

Corollaire A.5 Soit k ≤ r−1. On a :

min ||B − Z|| = σk+1 = ||B − Bk ||

Z:RangZ=k

Pk
où Bk = ui .~viT et ||Z|| = sup~x6=~0 ||Z.~
i=1 σi ~
x||Rm
x||Rn est la norme usuelle.
||~
Cela donne une mesure numérique du rang de B : quand σk+1 est de l'ordre de grandeur de la
pré ision ma hine, on en déduit que le rang numérique de B est k.

37 21 juin 2018
38 Référen es bibliographiques

Preuve. On a immédiatement U T .Bk .V = diagm,n (σ1 , ..., σk , 0, ...). D'où U T .(B − Bk ).V =
diagm,n (0, ..., 0, σk+1 , ..., σr , 0, ...) d'où ||B − Bk || = σk+1 , et en parti ulier min ||B − Z|| ≤
Z:RangZ=k
σk+1 . T
Soit Z une matri e m∗n de rang k . Don dim KerZ = n−k . Soit E = KerZ Vect{~v1 , ...~vk+1 }.
Comme E est l'interse tion d'un espa e de dimension n−k et d'un espa e de dimension k+1
n
dans R , et espa e n'est pas réduit à {0}.
x ∈ E tel que ||~x||Rn = 1 : on a ||(B−Z).~x||2Rm = ||B.~x||2Rm = || k+1 σi (~viT .~x)~ui ||2 =
P
Soit ~ i=1 P
Pk+1 2 T 2 k+1
i=1 σi (~ vi .~x) ar les ~ui sont orthonormaux, et don ||(B−Z).~x||2Rm ≥ σk+1 i=1 (~viT .~x)2 =
σk+1 ||~x||2 = σk+1 . D'où min ||B − Z|| ≥ σk+1 .
Z:RangZ=k

B Appli ation à la ondition inf-sup

Cas d'une ontrainte omme div~u = 0 qui s'é rit B.~u = 0 ave B matri e re tangle m ∗ n
asso iée à la forme bilinéaire b(vh , qh ) = (Bvh , qh ) lorsqu'on travaille ave des sous-espa es de
dimension nie Qh ⊂ Q et Vh ⊂ V .
T T
I i on a B au lieu de B , et on al ule les valeurs propres singulières de B (ou de B ), i.e. les
T
 
0 B
valeurs propres de la matri e .
B 0

Proposition B.1 Valeur de la onstante de la ondition inf-sup. On a :

b(vh , qh )
inf sup = σr ,
qh ∈Qh vh ∈Vh ||vh ||V ||qh ||Q

où σr est la plus petite valeur singulière non nulle de B.

Preuve. On aPB =
Pr Pr Pr
ui .~viT , d'où B.~x = P
i=1 σi ~ viT .~x) ~ui , d'où ~y T .B.~x = i=1 σi (~viT .~x) (~y T .~ui ).
i=1 σi (~
n j n
Soit~x = j=1 x ~vj ave ||~x||2 = 1 et ~y = k=1 y k ~uk ave ||~x||2 = 1 : on obtient ~y T .B.~x =
Pr i i
i=1 σi x y .
i
Et à ~ y xé le sup en ~x est immédiatement donné par xi = P σi2y 2 1 et donne ~y T .B.~x =
( i σi yi ) 2
1
Pr 2 2
P 2 2 1
P 1 σ y
i=1 i i = ( σ y
i i i ) 2 . D'où l'inf en ~
y est donné pour y
~ = ~ur et donne ~y T .B.~x = σr .
( i σi2 yi2 ) 2

Référen es bibliographiques

1 Arnold D., Brezzi F., Fortin M. : A stable nite element for the Stokes equations. Cal olo, 21,
pp. 337-344 (1984).

2 Babu²ka I. : The Finite Element Method with Lagrangian Multipliers. Numer. Math. 20, 179-
192 (1973).

3 Barbosa H., Hughes T. : The Finite Element Method with Lagrange Multipliers on the Boun-
dary : ir umventing the Babu²kaBrezzi Condition. Comp. Meths. in Appl. Me h. and Eng.,
85, 109128, 1991.

4 Ber ovier M., Pironneau O. : Error estimates for nite element method solution of the Stokes
problem in the primitive variables. Numer. Math., 33, pp. 211-224 (1977).

5 Brenner S.C., S ott L.R. : The Mathemati al Theory of Finite Element Methods. Text in
Applied Mathemati s, Springer-Verlag.

6 Brézis H. : Analyse fon tionnelle, théorie et appli ations. Colle tion mathématiques appliquées
pour la maîtrise, Masson (1983).

7 Brezzi F. : On the uniqueness, existen e and approximation of saddle point problem arising
from lagrangian multipliers. RAIRO Anal. Numer, 8, 129151 (1974).

8 Brezzi F., Fortin M. : Mixed and Hybrid Finite Element Methods. SpringerVerlag, 1991.

9 Brezzi F., Pitkäranta J. : On the stabilization of nite element approximations of the Stokes
equations. E ient Solutions of Ellipti Systems, Notes on Numeri al Fluid Me hani s, 10,
Ha kbush ed. (1984).

38 21 juin 2018
39 Référen es bibliographiques

10 Chapelle D. : Une formulation mixte de plaques où l'eort tran hant est appro hé dans son
espa e naturel. INRIA, rapport de re her he 2248, novembre 1993.

11 Ciarlet P.G. : Introdu tion à l'analyse numérique matri ielle et à l'optimisation. Masson, 1990.

12 Ciarlet P.G. : The Finite Element Method for Ellipti Problems. North-Holland (1977).

13 Crouzeix M., Raviart P.A. : Conforming and non- onforming nite element methods for solving
the stationary Stokes equations. R.A.I.R.O. Anal. Numer., 7, pp. 33-76 (1973).

14 Douglas J. Jr, Wang J. : An absolutely stabilized nite element method for the Stokes problem.
Math. Comp. 52, pp. 495-508, 1989.

15 Ern A., Guermond J.L. : Eléments nis : théorie, appli ations, mise en oeuvre. Mathématiques
& Appli ations 36, Springer 2002.

16 Fortin M. : Utilisation de la méthode des éléments nis en mé anique des uides. Cal alo, 12,
pp. 405-441 (1975).

17 Girault, V., Raviart, P. A. : Finite Element Methods for Navier-Stokes Equations. Springer-
Verlag (1986).

18 Golub H.G., Van Loan C.F. : Matrix Computations. John Hopkins, 3rd edition, 1996.

19 Goulaoui C. : Analyse fon tionnelle et al ul diérentiel. Ed. de l'E ole Polyte hnique, 1982.

20 Goulaoui C., Meyer Y. : Analyse fon tionnelle et al ul diérentiel. Ed. de l'E ole Polyte h-
nique, 1984.

21 Hughes T.J.R., Fran a L.P., Balestra M. : A new nite element formulation for omputational
uid dynami s : V. Cir umventing the Babu²kaBrezzi ondition : A stable PetrovGalerkin
formulation of the Stokes problem a omodating equal-order interpolations. C.M.A.M.E. 59,
pp. 85-99, 1986.

22 Leborgne G. : An Optimally Consistent Stabilization of the Inf-Sup Condition. Numer. Math.,

91, pp. 35-56 (2002).

23 Lions J.L., Magenes E. : Problèmes aux limites non homogènes et appli ations, Vol 1. Dunod
(1968).
◦
24 Nédéle J.C. : Méthode des éléments nis. Mathématiques et appli ations n 7. Ellipses (1991).

25 Nits he J. : Über ein Variationsprinzip zur Lösung von Diri hlet-Problemen bei Verwendung
von Teilräumen, die keinen Randbedingungen unterworfen sind. Abh. Math. Sem. Univ. Ham-
burg, 36 :915, 1970/1971.

26 Pitkäranta J. : Boundary Subspa es for the Finite Element Method with Lagrange Multipliers.
Numer. Math. 33, pp. 273189, 1979.

27 Raviart P.A., Thomas J.M. : Introdu tion à l'analyse numérique des équations aux dérivées
partielles. Colle tion mathématiques appliquées pour la maîtrise, Masson (1983).

28 S hwartz L. : Méthodes mathématiques pour les s ien es physiques. Colle tion enseignement
des s ien es, Hermann (1993).

29 S hweizer B. : www.mathematik.uni-dortmund.de/lsi/s hweizer/Preprints/perfe t ondrings-

preprint.pdf

30 Stenberg R. : On some te hniques for approximating boundary onditions in the nite ele-
ment method. Journal of Computational and Applied Mathemati s 63 (1995) 139-148. Et
http ://math.tkk./ rstenber/publi ations.htm.

31 Strang G. : Linear Algebra and its Appli ations. Har ourt Bra e & Company, 3rd edition
(1988).

32 Strang G., Fix J.F. : An Analysis of the Finite Element Method. Wellesley-Cambridge Press
(1973).

39 21 juin 2018