EPL

Avril 2023 Eléments finis

LEPL1110 Solution

So easy : Babus̆ka1 !

Ivo recherche le champ de température u(x, y) dans

un domaine carré Ω de côté de longueur deux
discrétisé au moyen de quatre éléments bilinéaires de y Température nulle
même taille unitaire. Le problème est stationnaire et
plan.

h i 6 7 8
∇ · k∇u(x, y) + f (x, y) = 0
III IV
La température des zone droite et supérieure est 3 4 5
imposée à une valeur nulle, tandis que les parties
gauche et inférieure sont parfaitement isolées. La I II
densité de production de chaleur f = 4 et la con- 0 1 2
x
ductivité thermique k = 6 sont constantes dans les 1
quatre éléments. Toutes les valeurs numériques sont
données dans des unités compatibles !

Isolant parfait
Il est possible de tirer habilement profit de la
symétrie du problème !
Plus précisément, Ivo vous demande :

1. Donner les unités de la conductivité thermique k et de la densité de production de chaleur f .

   
W W
Il suffit simplement d’écrire : Les unités de k et f sont et
mK m3

Attention, comme on écrit un bilan d’énergie,

les unités de l’équation sont celles de la dérivée temporelle de la quantité conservée : [W ] !
2. Ecrire la formulation faible du problème.

Il suffit d’écrire : Trouver u ∈ U tel que u · ∇u > = < u

< k∇b bf >, ∀b
u ∈ Ub

3. Définir l’espace de Sobolev H1 (Ω) lorsque Ω ⊂ R.

n Z  2  2 o
H1 (Ω) = v(x) : Ω → R tels que v(x) + v ′ (x) dx < ∞
Ω

1 Ivo Babus̆ka, né en 1926 nous a quitté ce 12 avril 2023. Il était un mathématicien tchéco-américain célèbre pour ses

travaux sur la méthode des éléments finis. Un de ses résultats les plus connus sur les éléments finis est la condition LBB
(Ladyzhenskaya-Babus̆ka-Brezzi) qui donne des conditions suffisantes pour assurer une formulation discrète stable.
4. Démontrer2 que la forme bilinéaire associée à la formulation faible est coercive pour cet espace.
Sans perte de généralité, on suppose qu’une condition essentielle homogène est imposée à l’origine
et on écrit :
Z x
v(x) − v(0) = v ′ (t) dt ∀x ∈ Ω
|{z} 0
=0

En utilisant l’inégalité de Cauchy : < v ′ , 1 >≤ ∥v ′ ∥ ∥1∥

? Z
h i2 x Z x h i2
v(x) ≤ dt v ′ (t) dt ∀x ∈ Ω
| 0{z } |0 {z }
≤ A ≤ ∥v ′ ∥20

En intégrant les deux expressions sur Ω

?
∥v∥20 ≤ C ∥v ′ ∥20

∥v∥20 + ∥v ′ ∥20 ≤ (C + 1) ∥v ′ ∥20

k h i
∥v∥20 + ∥v ′ ∥20 ≤ k ∥v ′ ∥20 □
1+C | {z } | {z }
2
| {z }
∥v∥1 a(v,v)
α

Et on a donc bien : ∃ α > 0 tel que a(u, u) ≥ α ∥u∥2 ∀u ∈ U

5. Donner3 les quatre fonctions de forme bilinéaires non-nulles τi (x, y) sur l’élément I.
6. Donner les dérivées de ces fonctions τi,x (x, y) et τi,y (x, y) sur l’élément I.
Les quatres fonctions de forme et leurs dérivées sont données par :

τ0 (x, y) = (1 − x)(1 − y) τ0,x = y−1 τ0,y = x−1

τ1 (x, y) = x(1 − y) τ1,x = 1−y τ1,y = −x
τ4 (x, y) = xy τ4,x = y τ4,y = x
τ3 (x, y) = (1 − x)y τ3,x = −y τ3,y = 1−x

Ces deux sous-questions étaient vraiment élémentaires et peu d’erreurs sont donc admises.
En particulier, il fallait bien veiller à obtenir les expressions correctes, en vérifiant que:
X X X
τi = 1, τi,x = τi,y = 0.

7. Calculer les matrices locales et membres de droite des éléments4 .

Il faut évidemment choisir les intégrales les plus simples à calculer et tirer profit de symétries !
En utilisant une numérotation locale circulaire, on calcule que :
2 On se limite bien au cas unidimensionnel, comme cela a été fait au cours et en séances d’exercices.
3 Dans ce problème, il est vraiment plus simple de ne pas utiliser l’isomorphisme vers l’élément parent !
4 Bien noter que les matrices locales et les membres de droites locaux des quatre éléments sont identiques :-)
 Z 1 Z 1
Ae11 = 2 x2 dxdy = 4
6
0 0

Z 1 Z 1 Z 1 Z 1
Ae12 = Ae14 = − y 2 dxdy + x(1 − x) dxdy = − 61
0 0 0 0

Z 1Z 1
Ae13 = 2 x(x − 1) dxdy = − 62
0 0

Z 1 Z 1
1
Bie = xy dxdy = 4
0 0

En tirant ensuite profit de toutes les symétries,

on en déduit directement que les matrices et membres de droites locaux s’écrivent :

   
4 −1 −2 −1 1
1  −1 4 −1 −2 
  1 
e
Aij = 6   Bie = 41  
−2 −1 4 −1   1 
−1 −2 −1 4 1

Observer que la somme de tous les éléments d’une ligne de la matrice vaut bien zéro !
8. Quelles sont les composantes de la matrice de raideur requises pour obtenir les quatre valeurs
nodales inconnues5 , en tenant compte des conditions essentielles homogènes du problème ?
En imposant les conditions essentielles sur cinq valeurs nodales,
il ne reste plus que quatre inconnues !
    
■ ■ ■ ■ U0 ■
 ■
 ■ ■ ■ 
 U1  
  ■ 


 1 
 U2  
  0 

 ■
 ■ ■ ■ 
 U3  
  ■ 

 ■
 ■ ■ ■ 
 U4 =
  ■ 


 1 
 U5  
  0 


 1 
 U6  
  0 

 1  U7   0 
1 U8 0

9. Assembler ces composantes (et uniquement celles-là !) du système linéaire global.

On assemble les matrices locales des 4 éléments.
Attention, de ne pas oublier les éléments extra-diagonaux !
    
4 −1 −1 −2 U0 1
k  −1 4+4 −2 −1 − 1
   U1  f
 1+1 
6  −1
  =  
−2 4+4 −1 − 1   U3  4  1+1 
−2 −1 − 1 −1 − 1 4 + 4 + 4 + 4 U4 1+1+1+1

Comme k = 6 et f = 4,

    
4 −1 −1 −2 U0 ? 1
 −1 8 −2 −2   U1   2 

 −1
  =  
−2 8 −2   U3   2 
−2 −2 −2 16 U4 4
5 Il est possible de répondre en représentant graphiquement la matrice et le vecteur et en cochant les termes adéquats !
*** 10. Obtenir6 les trois fractions correspondant aux quatre valeurs nodales inconnues U0 , U1 , U3 et U4 .

Par symétrie, on a que U1 = U3 , et on peut donc réduire le système :

    
4 −2 −2 U0 1
 −1 6 −2   U1  =  2 
−2 −4 16 U4 4

En remplacant L1 par L1 + 4L2 ,

En remplacant L3 par L3 − 2L2 ,

   ?  
0 22 −10 U0 9
 −1 6 −2   U1  =  2 
0 −16 20 U4 0

En remplacant L1 par L1 + 11L3 /8 ,

En remplacant L3 par L3 /4 ,

   ?  
0 0 17.5 U0 9
 −1 6 −2   U1  =  2 
0 −4 5 U4 0

?
18
U4 =
35

A partir de la dernière équation 4U1 = 5U4 , on déduit :

5 18 9
U1 = × =
4 35 14

Ensuite, on utilise la seconde équation U0 = −2 + 6U1 − 2U4 pour obtenir :

54 36 −70 + 135 − 36 −70 + 135 − 36 29

U0 = −2 + − = = =
14 35 35 35 35

29 9 18
Et, on peut donc conclure : U0 = U1 = U4 =
35 14 35

Bravo aux 12 étudiant(e)s sur 332 qui ont trouvé les 3 fractions : c’était encore trop facile !

6 Question réservée pour les étudiants qui trouvent l’interrogation trop simple et voudraient obtenir plus que 20/20 :-)