Méthodes quantitatives
pour ingénieurs
6MIG818-H22
Cours-2
Université du Québec à Chicoutimi
Département des sciences appliquées 1
3.2 Méthode des volumes finis
Soit l’équation de diffusion 1D stationnaire avec terme source
d d
( ) + S = 0
dx dx
d d
V dx ( dx ) dV + V S dV = 0
d d
xyz
( ) dxdydz + S dxdydz = 0
dx dx xyz
Application du théorème de divergence
xyz
• dxdydz = n • dS
S
2
3.2 Méthode des volumes finis
d d
xyz
( ) dxdydz + S dxdydz = 0
dx dx xyz
Application du théorème de divergence
Équation de difusion 1D stationnaire discrétisée par MVF
df df
(GA )e - (GA )w + SDxDyDz= 0
dx dx
Application de la méthode des différences finies + linéarisation terme source
G eAe ( jdEx-PEj P ) - G w Aw ( jdPx-WP
jW
) + SDxDyDz= 0
3
3.2 Méthode des volumes finis
G eAe ( jdEx-PEj P ) - G w Aw ( jdPx-WP
jW
)=0
Équation discrétisé par MVF & réarrangée dans la forme Ax=b
w e w e
− Aw W + Ae + Aw P − Ae E = 0
xWP xPE xWP xPE
aW aP aE
4
3.2 Méthode des volumes finis
Problème de conduction 1D stationnaire sans production
d dT
(k ) = 0 avec T = Ta à x = 0 & T = Tb à x = L
dx dx
Équation discrétisé par MVF
0
Noeud intérieur:équation discrétisé par MVF & réarrangée dans la forme Ax=b
kw ke kw e
− Aw TW + Ae + Aw TP − Ae TE = 0
xWP xPE xWP xPE
a
aW E
aW
aE
aP 5
3.2 Méthode des volumes finis
Problème de conduction 1D stationnaire sans production
Noeud intérieur:équation discrétisé par MVF & réarrangée dans la forme Ax=b
kw k k e
− Aw TW + e
Ae + w
Aw TP − Ae TE = 0
xWP xPE xWP xPE
aW a E
aW
aE
aP
Noeud intérieur: Forme indicielle pour les noeuds i=2,n-1
(ki −1 + ki ) (ki −1 + ki ) (ki + ki +1 ) (ki + ki +1 )
− 2 Ai Ti −1 + 2 Ai + 2
Ai Ti − 2 Ai Ti +1 = 0
( xi − xi −1 ) ( xi − xi −1 ) ( xi +1 − xi ) ( xi +1 − xi )
i−
a a i +1,i
ai ,i −1
1,i
ai ,i +1
ai ,i
6
3.2 Méthode des volumes finis
Problème de conduction 1D stationnaire sans production
Noeud frontière W (condition de Dirichlet):
kw k k e
− Aw TW + e
Ae + w
Aw TP − Ae TE = 0
xWP xPE xWP xPE
aW a E
aW
aE
aP
Forme indicielle pour les noeuds i=1
( k + k ) (k1 + k 2 )
1 2
A1 TWall − w + A1 T1 −
k1 k
− 1
A1 + 2 2 A1 T2 = 0
( x1 − xWall − w ) ( x1 − xWall − w ) ( x2 − x1 ) ( x2 − x1 )
1
a ,2
a1,1 1, 2a
7
3.2 Méthode des volumes finis
Problème de conduction 1D stationnaire sans production avec
condition de Neumann flux nul ou plan de symétrie
d dT W
(k ) = 0 avec q = 0 2 à x = 0 & T = Tb à x = L
dx dx m
Noeud intérieur: équation discrétisé par MVF & réarrangée dans la forme Ax=b
(ki −1 + ki ) (ki −1 + ki ) (ki + ki +1 ) (ki + ki +1 )
− 2 Ai Ti −1 + 2 Ai + 2 Ai Ti − 2 Ai Ti +1 = 0
( xi − xi −1 ) ( xi − xi −1 ) ( xi +1 − xi ) ( xi +1 − xi )
ai −1,i ai ,i +1
ai ,i −1
ai ,i +1
ai ,i 8
3.2 Méthode des volumes finis
Problème de conduction 1D stationnaire sans production avec
Condition de Neumann
(
kP + kE
2
) Ae (TE −TP
x PE ) −(
kW + k P
2
) Aw (
TP −TW
xW P )= 0
ke kw
( k1 + k 2 ) (k1 + k 2 )
2
A T − 2 A T = 0
( x2 − x1 ) 1 1 ( x2 − x1 ) 1 2
a
1, 2
a
a1,1 1, 2
9
3.2 Méthode des volumes finis
Soit l’équation de diffusion 1D transitoire avec terme source
d d d
= ( ) + S
dt dx dx
t + t
t
xyz
t
dxdyd zdt =
(1)
t + t t + t
t
xyz
( ) dxdydzdt + S dxdydzdt
x x t xyz
( 2)
t + t t + t
(1) :
t
xyz
t
dxdyd zdt =
t
t
dVdt
( ) ( )
t + t
E − P
e Ae xPE − w Aw xW P +t SdVV
=
P −W
=0
= ( t + t − t )dV
10
3.2 Méthode des volumes finis
Soit l’équation de diffusion 1D transitoire avec terme source
t + t
t
xyz
t
dxdyd zdt =
(1)
t + t t + t
t
xyz
(
x x
) dxdydzdt + S dxdydzdt
t xyz
( 2)
t + t t + t
(2) : (e Ae (
E − P
x PE ))dt + S Vdt
)− w Aw ( P −W
xW P
t t
( 2) : (1 − f ) A ( )− A ( ) t +
E − P P −W t
e e x PE w w xW P
f A ( )− A ( ) t + S Vt
E − P P −W t + t
e e x PE w w xW P
f = 0, 1 ,1
2
11
3.2 Méthode des volumes finis
Soit l’équation de diffusion 1D transitoire avec terme source
( P
t + t
− P )V = (1 − f ) e A
t
( Et − Pt
e x PE )− A ( )t +
w
Pt − Pt
w xW P
(
f e A e
Et +t − Pt +t
x PE )− A (
w w)t + S Vt
Pt +t −W
xW P
t + t
f = 0, 1 ,1
2
Soit l’équation de conduction 1D transitoire avec terme source
( )− k A ( )+
t + t
− TP )
t
(TP TEt −TPt TPt −TWt
c V = (1 − f ) k e Ae x PE w xW P
t
w
(
f ke A e
TEt + t −TPt + t
x PE )− k A (
w )+ S V
w
TPt + t −TWt + t
xW P T
f = 0, 1 ,1
2
12
3.2 Méthode des volumes finis
Soit l’équation de conduction 1D transitoire avec terme source
Schéma explicite (f=0); noeud intérieur
( )− k A ( )+
t + t
− TP )
t
(TP TEt −TPt TPt −TWt
c V = (1 − f ) k e Ae x PE w xW P
t
w
(
f ke A e
TEt + t −TPt + t
x PE )− k A (
w w )+ S V
TPt + t −TWt + t
xW P T
f = 0, 1 ,1
2
( )− k A ( )+ S V
t + t
− TP )
t
(TP TEt −TPt TPt −TWt
c V = ke Ae x PE w xW P
t
w T
c t + t
t V TP =
aT
kw t ke kw c
Aw TW − Ae + Aw − V TPt + e Ae TEt + ST V
xWP xPE xWP t
xPE
a
aW E aW
aT
aE
13
aP
3.2 Méthode des volumes finis
Soit l’équation de conduction 1D transitoire avec terme source
Schéma implicite (f=1); noeud intérieur
( )− k A ( )+
t + t
− TP )
t
(TP TEt −TPt TPt −TWt
c V = (1 − f ) k e Ae x PE w xW P
t
w
(
f ke A e
TEt + t −TPt + t
x PE )− k A (
w w )+ S V
TPt + t −TWt + t
xW P T
f = 0, 1 ,1
2
( )− k A ( )+ S V
t + t
− TP )
t
(TP TEt + t −TPt + t TPt + t −TWt + t
c V = ke Ae x PE xW P
t
w w T
kw t + t c kw ke t + t t + t
− Aw TW + V + Aw + Ae TP − e
Ae TE =
xWP t xWP xPE xPE
aW aT aW aE aE
aP
c
V TPt + + ST V
t
14
aT
3. Introduction aux méthodes numériques pour la simulation de phénomènes
avec diffusion
3.2 Méthode des volumes finis
Linéarisation du terme source
S V = S u + S P P
Équation de diffusion 1D discrétisée pour chaque VF au point P
e A (
E − P
e x PE )− A (
w
P −W
w xW P )+ (S u + S PP ) = 0
En réarrangeant
e w w e
Ae + Aw − S P P = Aw W + Ae E + Su
xPE xWP xWP xPE
Équation de diffusion 1D discrétisée sous la forme linéaire
a pP = aW W + aEE + Su
15
3. Introduction aux méthodes numériques pour la simulation de phénomènes
avec diffusion
3.2 Méthode des volumes finis
Équation de diffusion 1D discrétisée sous la forme linéaire
e A + w A − S = w A + e A + S
xPE e xWP w P
P xWP w W xPE e E u
aP aW aE
a pP = aW W + aEE + Su
aW aE aP S P Su
w
Aw E
AE aW + aE − SP
xWP xEP
16
3. Introduction aux méthodes numériques pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
Exemple-1
17
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
18
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
19
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
Exemple-2
20
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
21
3. Introduction aux méthodes numériques pour la simulation de phénomènes
avec diffusion
3.2 Méthode des volumes finis
22
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
23
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
24
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
25
3. Introduction aux méth. Num. pour la simulation
de phénomènes avec diffusion
3.2 Méthode des volumes finis
26
Numérotation (nœud intérieur)
système matriciel vs système global
Découpage de la géométrie(maillage en 2D)
NR=4 rangées & NC=4 colonnes
3. Introduction aux méth. Num. pour la simulation
(i=1,j=1) (i=1,j=2) (i=1,j=3) (i=1,j=4)
3.3 Représentation matricielle
#1 #2 #3 #4
de phénomènes avec diffusion
(i=2,j=1) (i=2,j=2) (i=2,j=3) (i=2,j=4)
#5 #6 #7 #8
y=1000m
(i=3,j=1) (i=3,j=2) (i=3,j=3) x=1000m
#9 #10 #11
X
28
29
Système A•T=b ou A•C=b
1 2 3 4 5 6 7 8 9 10 11 12
1 x x x T1 x
2 x x x x T2 x
3 x x x x T3 x
4 x x x T4 x
5 x x x x T5 x
x x x x x
6
•
T6
7 x
x
x x x
x x x
x T7 =
8 T8 x
9 x x x T9 x
10 x x x x T10 x
11 x x x x T11 x
12 x x x T12 x
30
Numérotation (nœud intérieur)
système matriciel vs système global
NR=3; NC=4; NR=3; NC=4;
DX=; DY= ; DX=; DY= ;
Pour i=2 jusqu’à NR-1 i=2; j=2;
Pour j=2 jusqu’à NC-1 Imat=6;
imat= j + NCx(i-1); A(6,2)=
(* NORD *) A(6,6)=
A(imat , imat - NC)= ; A(6,7)=
(* WEST *) A(7,8)=
A(imat , imat - 1)= ; A(7,10)=
(* P *)
A(imat , imat )= ; i=2; j=3;
(* EAST *) Imat=7;
A(imat , imat +1)= ; A(7,3)=
(* SUD *) A(7,7)=
A(imat , imat + NC)= ; A(7,8)=
A(7,9)=
A(7,11)=
