Cours

Méthode des différences finies

• Technique de discrétisation en 1D
• Construction du système
• Prise en compte des conditions aux limites
• Notion de convergence
• Extension au 2D

1 31/01/2020
Méthode des différences finies
Objectif : transformer une équation « continue » valable sur un domaine
continu en un système à N équations pour N inconnues associées
à un domaine discret appelé maillage

 k11 k12 k13   u1   f1 

u  2u k    
L(u, , 2 ...)  f  0
 21 k22 k23  u2    f 2 
x x
+ conditions aux limites  k31 k32 k33  u3   f3 

Méthode : écrire sous forme discrète (i-1, i, i+1 …) tous les termes de
dérivées présents dans l’équation d’équilibre appliquée en i
ainsi que dans les C.L.

2 31/01/2020
Différences finies 1D : méthode générale
d 2T  x 
Reprenons l’exemple de thermique 1D régi par : k 2
 f  0,  x   0, L 
dx
T ( x  0)  30
( L)  h T  L   Text 
dT
q  L   k
dx
1. On discrétise le domaine en « N » nœuds (maillage) :

A domaine discret, équation « discrète » !

d 2T
2. On applique alors cette équation au nœud « i » : k 2  fi  0  i  1,..., N
dx i

A ce stade, il nous faut donc discrétiser le terme de dérivée seconde !

31/01/2020 3
 Discrétisation des termes de dérivées
• Utilisation des développements limités :

notation indicielle

d T x d T x
2 2 3 3
dT
T  x  x   T  i  1  T (i )  x    x (...)
4
2 3
(1)
dx i dx i
2 dx i
6
d T x d T x
2 2 3 3
dT
T  x  x   T  i  1  T (i )  x    x (...)
4
2 3
(2)
dx i dx i
2 dx i
6

• On combine ces deux équations. Par exemple, la somme de (1) et de (2) :

d T x
2 2

T  i  1  T  i  1  2 T (i )  2 2
 x (...)
4
1   2 
dx i
2
dT
2
T  i  1  2 T (i )  T  i  1
permet d’isoler :   x (...)
2

dx
2
i
x 2
représentatif de l’ordre de tous les termes tronqués
31/01/2020 4
Principales formes discrètes à connaître
En combinant de différentes manières, on obtient ainsi les approximations
discrètes suivantes : Précision du schéma

dT T T
(1)   i 1 i x ... Décentré droit
dx i x
dT T T
(2)   i i 1 x ... Décentré gauche
dx i x
d 2T T  2Ti  Ti 1
1   2   2  i 1 x 2 ... Centré
dx i x 2
dT T T
1   2    i 1 i 1 x 2 ... Centré
dx i 2x

Nouvelle notation : T(i+1)=Ti+1 Termes Type

tronqués

5 31/01/2020
Interprétation graphique

Discrétisation centrée :
relation dans laquelle les contributions des valeurs nodales de part et d'autre du point considéré
(noeud i) sont équivalentes.

Discrétisation décentrée :
relation dans laquelle les contributions des valeurs nodales de part et d'autre du point considéré
(noeud i) ne sont pas équivalentes.

6 31/01/2020
Construction globale du système
Ti 1  2Ti  Ti 1
• La relation discrète finalement obtenue s’écrit
k
x 2 :
 fi  0

x 2
ou encore : Ti 1  2Ti  Ti 1   fi
k

0 0 0 0 ... 0 0 0   T   0 
•
1

Elle est applicablex 2 seulement  aux nœuds   …, N-1

i=2,  : 
i  2: T1  2T2  T3   f2 1 2 1 0 ... 0 0 0   T    f  x 
2

2 2
k
k     
0 1 2 1 ... 0 0 0   T    f 

2

x 2 x
i  3: T2  2T3  T4   f3     
3 3
k

k
    
    
0 0 0 0 ... 1 2 1  T   f  x 2

x k
2
i  N  1: TN  2  2TN 1  TN   f ( N 1)    
N 1 N 1


k 0 0 0 0 ... 0 0 0   T  
N
0 

Écriture sous forme matricielle

31/01/2020 7
Condition à la limite de type DIRICHLET
On a la condition suivante : T ( x  0)  T  30
1

Méthode : on ajoute :
1. un terme unité « 1 » sur la diagonale du nœud concerné
2. la valeur connue dans le 2nd membre

1 0 0 0 ... 0 0  T
0 1   30 
    
1 2 1 0 ... 0 0 0  T  f  x 2

2 2
k
    
0 0   T x 
0 1 2    f
2
1 ... 0
 
3
k

3


    
    
0 0 0 2 1  T   f  x 2

k
0 ... 1
 
N 1 N 1

  
0 0 0 0 ... 0 0 0   T N
  0 

8 31/01/2020
 Condition à la limite de type CAUCHY (1/2)
avec noeud fictif !
( L)  h T  L   Text 
dT
On a la condition suivante : q  L   k
dx

Méthode : on discrétise le terme de dérivée présent dans la condition à la

limite (aussi appelée condition de type « flux »).

• Avec noeud fictif : plus long

T T On applique la relation d’équilibre
k précis !h T  T 
mais N 1 N 1
discrète en N car le nœud N+1
2x
N ext

existe :
hx x 2
 TN 1  TN 1  2 TN  Text  i  N : TN 1  2TN  TN 1   fN
k k
 hx  x 2
x
 2TN 1   2  2  N
T   f N  2 h Text
 k  k k
31/01/2020 9
 Condition à la limite de type CAUCHY (2/2)
sans noeud fictif !

• Sans noeud fictif : rapide mais perte en précision !

TN  TN 1
k  h TN  Text  (précis ordre 1)
On a recours à une formule décentrée pour la CL : x

conduisant ainsi à :  x  x
TN 1  1  h  TN   h Text
 k  k

+ : rapide à mettre en oeuvre

- : on diminue la précision globale du schéma

31/01/2020 10
Système final à résoudre
1 0 0 0 ... 0 0 0   30 
  T 1   
1 2 1     f x
2
0 ... 0 0 0 
  T   k
2
2 
0 1 2     f x 
2
1 ... 0 0 0
   T   k 
3
3

   
    
0 0 0 0 ... 1 2 1     f x
2

  T N 1   N 1
k 
  hx      x 2
x 
0 0 0 0 ... 0 2  2  2    T   f  2h T 
  k  
N
 k
N
k
ext


Rem : ce système est basé sur le traitement de la CL avec nœud fictif

11 31/01/2020
Affichage et post-traitement de la solution
Pour des systèmes de tailles supérieures à 3-4, on a généralement recours à des outils informatiques dédiés
à la résolution et l’affichage.

Apprentissage de l’outil Matlab lors des séances TP de NF04

12 31/01/2020
Pour résumer …
• Mailler le domaine

• Discrétiser l’équation d’équilibre et les conditions aux limites :

• En remplaçant toutes les dérivées par leur forme discrète

• Construire le système global

• En appliquant les équations discrètes sur les nœuds concernés

• Résoudre le système (voir TP et TD encadrés sous Matlab)

• Post-traiter :
• Tracer la solution
• Calculer les variables dérivées : flux (thermique), contrainte (méca) …

13 31/01/2020
Fiabilité du modèle : notion de convergence
Modèle Modèle
mathématique numérique
(continu) (algébrique)
d 2T  x  Ti 1  2Ti  Ti 1
k  f 0 k  fi  0
dx 2 x 2

Erreur introduite en négligeant les termes des

développements limités à partir d’un certain ordre

Question : comment s’assurer que l’équation discrète est représentative,

en termes de phénomènes physiques, de l’équation de départ ?

Méca. Flu., thermique : transport, diffusion …

MMC : traction, flexion, dynamique …

Idée : le comportement du modèle numérique doit converger vers le

comportement du modèle mathématique (censé être proche du réel …).

14 31/01/2020
Notion de convergence
Méthode : s’assurer de la propriété de CONVERGENCE de l’équation discrète.

Théorème de LAX :

Convergence = consistance + stabilité

Comportement numérique
proche du « réel » Absence d’oscillations
parasites

15 31/01/2020
Notion de consistance
• Définition : on appelle erreur de troncature t, l’ensemble des termes
négligés dans les développements limités lors de
l’obtention d’une équation (ou schéma) discrète

Il est en effet possible d’écrire :

Équation continue = Équation discrète + t

 Définition : un schéma est dit consistant si son erreur de troncature tend
vers 0 lorsque le pas x tend vers 0
lim t  0
x 0

31/01/2020 16
Exemple de calcul de l’erreur de troncature
Considérons les développements limités suivants :

d T x d T x d T x
2 2 3 3 4 4
dT
Ti 1  Ti  x     x (...)
5

x 2
2 3 4
dx dx 2 dx 6 dx 24
Ti 1  2Ti  Ti 1  
i i i i
fi
d T x d T x d T x
2 2 3 3 4 4
dT k
Ti 1  Ti  x  2
 3
 4
 x 5 (...)
dx i dx i
2 dx i
6 dx i
24

que l’on injecte dans l’équation discrète.

Ce qui conduit à : dT
2
x 2
d T x
4 4

2
x  4  x (...)   fi 2 6

dx i dx i 12k k
soit : dT
2
x 2  d 4T 1 
k 2  fi    x ...   0
2

dx i k  dx 4 i 12 
Equation continue en i Erreur de troncature

Conclusion : le schéma est bien consistant avec l’équation de départ

Remarque : la solution par différences finies sera mathématiquement
exacte dans ce cas précis. La solution math. est quadratique d’où t = 0 !
31/01/2020 17
Effets « visibles » de l’erreur de troncature
Le comportement graphique de la solution est un indicateur des
effets de l’erreur de troncature

Le schéma est dit DISPERSIF si Le schéma est dit DIFFUSIF si

des dérivées impaires apparaissent. des dérivées paires apparaissent.
Effets néfastes pouvant entraîner Effets bénéfiques mais pouvant
l’instabilité des résultats diminuer la précision des résultats

18 31/01/2020
« Notion » sur la stabilité d’un schéma
Définition : la stabilité est la propriété de contrôler toute
perturbation (numérique dans notre cas)
introduite de manière accidentelle.
 Un schéma est dit STABLE si la perturbation diminue ou mieux, disparaît.
 Un schéma est dit INSTABLE si la perturbation augmente.

Concrètement, apparition d’oscillations parasites (changement du signe

de la pente d’un nœud à l’autre).
(L’étude de la stabilité sera développée ultérieurement.)
19 31/01/2020
Extension à 2 dimensions (2D)
Thermique : exemple d’une plaque rectangulaire soumises à différentes
conditions aux limites.

Définition du contour du domaine et génération d’un maillage quadrillé :

q.n  h T  Text  W / m 2  (Cauchy)

T  T0  K  (Dirichlet)

x  y

q.n   W / m 2  (Neumann)

Rem : [Link]
n le flux normal à la paroi (normale vers l’extérieur)
31/01/2020 20
Différences finies 2D : T  x, y  Ti, j
L’équation de la chaleur 2D est la suivante
.q : f  0,   x, y   S (1)
Divergence

La loi de comportement est : q  k grad T  k T  x, y  (2)

Flux thermique

Insertion de éq.(2) dans2éq.(1) : 2

  T  x, y   T  x , y  
k   f   0,  x   0, L 
 x 2
y 2


 2T Ti 1, j  2Ti , j  Ti 1, j  2T Ti , j 1  2Ti , j  Ti , j 1

  x 2
...   y 2 ...
x 2 i, j
x 2
y 2 i, j
y 2

 T  2Ti , j  Ti 1, j Ti , j 1  2Ti , j  Ti , j 1 

k  i 1, j    fi , j  0,
 x 2
y 2

 i  2,..., N  1,  j  2,..., M  1
31/01/2020 21
Construction du système
• Balayer les lignes les unes
après les autres et appliquer    T1,1   
l’équation discrète si possible     
   T1,2   
    
   T1,3   
 K    F 
    
• Appliquer les conditions aux 

   
 TN , M 1   
limites discrètes     
   TN , M   

• Résoudre et post-traiter les

solutions
31/01/2020 22