NF04 Cours2
NF04 Cours2
NF04 Cours2
• Technique de discrétisation en 1D
• Construction du système
• Prise en compte des conditions aux limites
• Notion de convergence
• Extension au 2D
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.
d 2T
2. On applique alors cette équation au nœud « i » : k 2 fi 0 i 1,..., N
dx i
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
d T x
2 2
T i 1 T i 1 2 T (i ) 2 2
x (...)
4
1 2
dx i
2
permet d’isoler : dT
2
T i 1 2 T (i ) T i 1
x (...)
2
dx
2
i
x 2
représentatif de l’ordre de tous les termes tronqués
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 2x
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.
x 2
ou encore : Ti 1 2Ti Ti 1 fi
k
Elle est applicable seulement aux nœuds i=2, …, N-1 :
0 0 0 0 ... 0 0 T
0 1 0
x 2
i 2: T1 2T2 T3 f2 1 2 1 0 ... 0 0 0 T f x 2
2 2
k
k
0 0 T x
0 1 2 f
2
x 2 1 ... 0
k
i 3: T2 2T3 T4 f3
3
3
k
x 2 0 0 0 0 ... 1 2 1 T f x 2
i N 1: TN 2 2TN 1 TN
N 1 N 1
k
f ( N 1)
k 0 0 0 0 ... 0 0 0 T N
0
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 0 T 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
TN TN 1
On a recours à une formule décentrée pour la CL : k h TN Text (précis ordre 1)
x
conduisant ainsi à :
x x
TN 1 1 h TN h Text
k k
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
hx x 2
x
0 0 0 0 ... 0 2 2 2 T f 2h T
k
N
k
N
k
ext
Mailler le domaine
Post-traiter :
Tracer la solution
Calculer les variables dérivées : flux (thermique), contrainte (méca) …
Théorème de LAX :
Comportement numérique
proche du « réel » Absence d’oscillations
parasites
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
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
T T0 K (Dirichlet)
x y