Exercice
Utilisation de la méthode de volumes finis pour résoudre le problème de diffusion de
l’équation de la chaleur dans le cas stationnaire de mur fini.
Intégrer et discrétiser l’équation de conservation sur chaque volume pour obtenir le système
matriciel ?
Solution
L’équation de la chaleur :
+ ∇.⃗ U⃗ = ∇.⃗ ∇⃗ + (01)
T : Température
ρ : La masse volumique
Cp : Capacité calorifique
t : Le temps
U : Vitesse
k : Conductivité thermique
S : source de chaleur
Régime stationnaire et pas de transport alors, il reste que le terme de diffusion et le terme
source.
L’équation (01) de la chaleur stationnaire en diffusion pur devient :
∇.⃗ ∇⃗ + =0 (02)
Le transfert de la chaleur a lieur uniquement dans la direction perpendiculaire au mur
Les plans de mur sont des isothermes
= =0
Mur 1D
Le transfert de la chaleur en diffusion pur 1D
+ =0 (03)
1
Discrétisation en volume finis
Le mur est découpé en volumes qui sont des tranches
Les nœuds sont placés au centre des volumes.
Des nœuds sont placés sur les deux surfaces du mur pour faciliter la formulation des
conditions aux limites. Dans ce cas les nœuds ne sont pas au centre des volumes.
Sur chaque volume on intègre l’équation (02) de conservation de la chaleur
Par convention on notera en MAGISCULE les grandeurs relatives aux volumes et aux
nœuds et on notera Minuscule les grandeurs relatives aux frontières des volumes. Les faces
est et ouest() sont désigné par w et e, les distances séparant les nœuds et ses voisins sont
notées δxw, δxe
w e
W P E
Δx
⃗ ⃗
δxw δxe
Intégrer sur un volume l’équation de la chaleur (02)
∇.⃗ ∇⃗ + =0 (04)
2
∇.⃗ ∇⃗ = n.⃗ ∇⃗
∇.⃗ ∇⃗ = ⃗. ∇⃗ = ⃗. ∇⃗
∇.⃗ ∇⃗ = − − −
L’équation (04) devient pour un volume:
− − − + =0 (05)
Discrétisation de l’intégrale :
Les flux diffusifs sur les faces sont approximés en discrétisant les expressions des gradients
sur le flux.
( − )
=
( − )
=
L’équation (05) reformulée sous la forme suivante :
( − ) ( − )
− + =0 (06)
Le terme source est discrétisé par sa valeur moyenne :
3
= ̅∆
L’équation (06) est devient alors
( − ) ( − )
− + ̅∆ = 0 (07)
Les coefficients.
Formalisation de l’équation sur le volume
= + + (08)
= = = ̅∆
= +
, , et sont appelés les coefficients de l’équation de la chaleur discrétisé sur le volume
Pour chaque volume i de mur on a une équation de type :
= + +
() ()= ( ) ( + 1) + ( ) ( − 1) + ( )
4
Les conditions aux limites
La formulation précédente non valable pour les volumes aux limites
Cas particulier des conditions aux limites côté ouest (West)
Le flux sur la face West ne s’exprime plus comme un flux diffusif et Δx est remplacé par Δxb
comme il est indiqué dans la figure suivante :
( − ) ( − )
− + ̅∆ = 0
( − )
− ⃗. Φ⃗ + ̅∆ = 0 (09)
5
Le flux sur la face West est le flux sur la face extérieur du mur
( − )
− Φ + ̅∆ = 0 (10)
Si
La condition à la limite de convection + flux imposé sur la face West alors :
=ℎ ( − ) +
Avec :
hw est le coefficient d’échange convectif
TP est la température de la surface du mur
Text est la température du fluide à l’infini
φw est le flux imposé sur la face West
SMur est la surface de mur qui est égale 1 m2 (Mur 1D)
L’équation (10) devient :
( − )
−ℎ ( − )+ + ̅∆ =0 (11)
Et les constantes :
6
= =0 = +
= ℎ + + ∆
= + +
Les conditions aux limites coté est.
C’est la même méthode qu’on a utilisé pour le côté West :
( − ) ( − )
− + ̅∆ = 0
Construction du système linéaire de N équations et N inconnues
Pour chaque volume i du mur on a une équation de type :
() ()= ( ) ( + 1) + ( ) ( − 1) + ( )
= 1 1 1 = 1 2 + 1 0 + 1 ( 1
= 0)
= 2 2 2 = 2 3 + 2 1 + 2
.
.
.
= +1 + −1 +
= −1 = + +
= = + + ( = 0)
