Formule de Green, formulation variationnelle (séance 2)
Laurent Monasse∗
∗ Université Côte d’Azur, Inria, CNRS, LJAD, Nice
16 décembre 2024
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 1/8
Ouvert régulier
On va donner un sens aux formules d’intégration par partie :
Z Z
∂f
= fni .
Ω ∂xi ∂Ω
Première question : qu’est-ce que la normale ?
Hypothèse : Ω est un domaine borné de Rd de frontière ∂Ω suffisamment
régulière.
Définition : ∂Ω est suffisamment régulière si
Ω est localement d’un seul côté de ∂Ω ;
On peut localement cartographier Ω par une application lipschitzienne : pour
tout point x ∈ ∂Ω, il existe un voisinage W de x dans Rd , un voisinage V de
0 dans Rd , et une bijection lipschitzienne Ψ : V → W telle que
Ψ−1 (Ω ∩ W ) = {u ∈ V , u1 > 0}.
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 2/8
Ouvert régulier
Définition : ∂Ω est suffisamment régulière si
Ω est localement d’un seul côté de ∂Ω ;
On peut localement cartographier Ω par une application lipschitzienne : pour
tout point x ∈ ∂Ω, il existe un voisinage W de x dans Rd , un voisinage V de
0 dans Rd , et une bijection lipschitzienne Ψ : V → W telle que
Ψ−1 (Ω ∩ W ) = {u ∈ V , u1 > 0}.
Exemples :
→
−
n
Ψ
x
Ω
n̂
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 2/8
Ouvert régulier
Définition : ∂Ω est suffisamment régulière si
Ω est localement d’un seul côté de ∂Ω ;
On peut localement cartographier Ω par une application lipschitzienne : pour
tout point x ∈ ∂Ω, il existe un voisinage W de x dans Rd , un voisinage V de
0 dans Rd , et une bijection lipschitzienne Ψ : V → W telle que
Ψ−1 (Ω ∩ W ) = {u ∈ V , u1 > 0}.
Exemples :
Ψ
x
Ω
n̂
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 2/8
Ouvert régulier
Définition : ∂Ω est suffisamment régulière si
Ω est localement d’un seul côté de ∂Ω ;
On peut localement cartographier Ω par une application lipschitzienne : pour
tout point x ∈ ∂Ω, il existe un voisinage W de x dans Rd , un voisinage V de
0 dans Rd , et une bijection lipschitzienne Ψ : V → W telle que
Ψ−1 (Ω ∩ W ) = {u ∈ V , u1 > 0}.
Contre-exemples :
Ψ
x
Ω
n̂
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 2/8
Ouvert régulier
Définition : ∂Ω est suffisamment régulière si
Ω est localement d’un seul côté de ∂Ω ;
On peut localement cartographier Ω par une application lipschitzienne : pour
tout point x ∈ ∂Ω, il existe un voisinage W de x dans Rd , un voisinage V de
0 dans Rd , et une bijection lipschitzienne Ψ : V → W telle que
Ψ−1 (Ω ∩ W ) = {u ∈ V , u1 > 0}.
Contre-exemples :
Ψ
x
Ω
n̂
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 2/8
Formules de Green (sur C 1 )
→
−
n
Soit f ∈ C 1 (Ω̄). Alors Z Z
∂f
dx = fni ds.
Ω ∂xi ∂Ω
Corollaire 1 : pour u, v ∈ C 1 (Ω̄),
Z Z
∂u ∂v
v +u dx = uvni ds
Ω ∂xi ∂xi ∂Ω
Corollaire 2 : pour →
−
u ∈ C 1 (Ω̄)d , v ∈ C 1 (Ω̄),
→
−
Z Z
(div →
−
uv +→ −u · ∇v )dx = v→
−u ·→−n ds.
Ω ∂Ω
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 3/8
Formules de Green (sur C 1 )
→
−
n
Soit f ∈ C 1 (Ω̄). Alors Z Z
∂f
dx = fni ds.
Ω ∂xi ∂Ω
Corollaire 3 : pour u ∈ C 2 (Ω̄), v ∈ C 1 (Ω̄),
Z 2 Z
∂ u ∂u ∂v ∂u
2
v+ dx = vni ds
Ω ∂x i ∂x i ∂x i ∂Ω ∂x i
Corollaire 4 : pour u ∈ C 2 (Ω̄), v ∈ C 1 (Ω̄),
→
− →
− →
−
Z Z
(∆uv + ∇u · ∇v )dx = v ∇u · →
−
n ds.
Ω ∂Ω
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 3/8
Formulation variationnelle
Formulation forte : trouver u ∈ C 2 (Ω̄) tel que
−∆u = f dans Ω
u = u0 sur ΓD
→
−
∇u · n = g sur ΓN .
Multiplication par une fonction v ∈ C 1 (Ω̄) et intégration :
Z Z
−∆uv = fv .
Ω Ω
Formule de Green :
Z Z Z
∇u · ∇v − v ∇u · →
−
n = fv .
Ω ∂Ω Ω
Conditions aux limites : on impose v = 0 sur ΓD , et on utilise la condition de
Neumann sur u sur ΓN :
Z Z Z
∇u · ∇v = fv + gv
Ω Ω ΓN
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 4/8
Formulation variationnelle
Formellement :
Trouver u ∈ C 2 (Ω̄) tel que
−∆u = f dans Ω
u = u0 sur ΓD
∇u · →
−
n = g sur ΓN .
⇓
Trouver u ∈ C (Ω̄) tel que u = u0 sur ΓD et tel que pour tout v ∈ C 1 (Ω̄)
1
avec v = 0 sur ΓD , Z Z Z
∇u · ∇v = fv + gv
Ω Ω ΓN
Quelle régularité pour la solution ?
Les deux problèmes sont-ils équivalents ?
Existence et unicité ?
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 5/8
Existence et unicité
Trouver u ∈ C 1 (Ω̄) tel que u = u0 sur ΓD et tel que pour tout v ∈ C 1 (Ω̄)
avec v = 0 sur ΓD , Z Z Z
∇u · ∇v = fv + gv
Ω Ω ΓN
Unicité : soient u1 et u2 deux solutions, alors
Z
1
∀v ∈ C (Ω̄), v|ΓD = 0, (∇u1 − ∇u2 ) · ∇v = 0.
Ω
Si on prend v = u1 − u2 , v = 0 sur ΓD car u1 = u2 = u0 sur ΓD :
Z
|∇(u1 − u2 )|2 = k∇(u1 − u2 )k2L2 = 0.
Ω
Donc u1 − u2 = C avec C une constante sur Ω (connexe). Comme
u1 = u2 = u0 sur ΓD , C = 0 et la solution est unique.
Existence : on verra dans la suite le théorème de Lax-Milgram qui donne
l’existence. Mais l’espace sur lequel la formulation variationnelle est posée
doit être complet (espace de Hilbert).
Mais C 1 (Ω̄) n’est pas complet pour la norme k · kL2 . . .
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 6/8
Espaces fonctionnels
En général, on n’a pas existence d’une solution dans C 1 (Ω̄).
En revanche, on va avoir existence de la solution dans H 1 (Ω) (espace de
Hilbert pour k · kH 1 )
Équivalence entre formulation forte et formulation faible dans H 1 (Ω).
Ingrédients :
Densité de D(Ω̄) dans H 1 (Ω) ou de D(Ω) dans H01 (Ω) (selon les conditions
aux limites)
Continuité dans H 1 (Ω) de la forme bilinéaire (u, Rv ) 7→ ΩR∇u · ∇v
R
Continuité dans H 1 (Ω) de la forme linéaire v 7→ Ω fv + Γ gv
N
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 7/8
Cas simple : Dirichlet homogène
Trouver u ∈ H 1 (Ω) tel que
−∆u = f dans D0 (Ω)
u = 0 sur ∂Ω.
~
Trouver u ∈ H01 (Ω) tel que pour tout v ∈ H01 (Ω),
Z Z
∇u · ∇v = fv .
Ω Ω
D(Ω) est dense dans H01 (Ω) : on peut prendre ϕn ∈ D(Ω) avec ϕn →n→∞ u
dans H 1 (Ω),
R et on va tester contre ϕn .
(u, v ) 7→ Ω ∇u · ∇v est continu dans H 1 (Ω) :
Z
∇u · ∇v ≤ k∇ukL2 k∇v kL2 ≤ kukH 1 kv kH 1
Ω C.S.
fv est continu dans H 1 (Ω) :
R
v 7→ Ω
Z
fv ≤ kf kL2 kv kL2 ≤ kf kL2 kv kH 1
Ω C.S.
L. Monasse (Inria & LJAD, Nice) Méthodes numériques pour les EDPs 16/12/2024 8/8
