Chapitre 4 : Méthode des volumes nis

Y. EL KHCHINE Mod-Syst-Energ 1 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)
• Modèle mathématique continu est discrétisé en s'appuyant sur
une(des) méthode(s) numérique(s)

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)
• Modèle mathématique continu est discrétisé en s'appuyant sur
une(des) méthode(s) numérique(s)
• Equations discrétisées sont approximées à l'aide des schémas
numériques appropriés, l'algorithme de résolution est établie

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)
• Modèle mathématique continu est discrétisé en s'appuyant sur
une(des) méthode(s) numérique(s)
• Equations discrétisées sont approximées à l'aide des schémas
numériques appropriés, l'algorithme de résolution est établie
• Algorithme est codé

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)
• Modèle mathématique continu est discrétisé en s'appuyant sur
une(des) méthode(s) numérique(s)
• Equations discrétisées sont approximées à l'aide des schémas
numériques appropriés, l'algorithme de résolution est établie
• Algorithme est codé
• Code est executé sur un ordinateur

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - méthodologie

Résolution numérique des problèmes de la Mécanique des Fluides

Grandes étapes :
• Problème physique continu est décrit par un modèle mathématique
continue (mis en équations)
• Modèle mathématique continu est discrétisé en s'appuyant sur
une(des) méthode(s) numérique(s)
• Equations discrétisées sont approximées à l'aide des schémas
numériques appropriés, l'algorithme de résolution est établie
• Algorithme est codé
• Code est executé sur un ordinateur
• Si tout va bien, la solution approchée du problème initial est obtenue.

Y. EL KHCHINE Mod-Syst-Energ 2 / 57
CFD - principales méthodes

Diérences nis
Appoximation des dérivées intervenantes dans les équations à l'aide de
devéloppement en série de Taylor

Y. EL KHCHINE Mod-Syst-Energ 3 / 57
CFD - principales méthodes

Diérences nis
Appoximation des dérivées intervenantes dans les équations à l'aide de
devéloppement en série de Taylor

Eléments nis
Détermination d'un champ local à attribue à chaque sous domaine (élément)
pour que le champ global obtenu par juxtaposition de ces champs locaux soit
proche de la solution du problème (bilan global).

Y. EL KHCHINE Mod-Syst-Energ 3 / 57
CFD - principales méthodes

Diérences nis
Appoximation des dérivées intervenantes dans les équations à l'aide de
devéloppement en série de Taylor

Eléments nis
Détermination d'un champ local à attribue à chaque sous domaine (élément)
pour que le champ global obtenu par juxtaposition de ces champs locaux soit
proche de la solution du problème (bilan global).

Volumes nis
Bilan local des ux dans un petit volume de contrôle.

Y. EL KHCHINE Mod-Syst-Energ 3 / 57
Pourquoi utiliser l'approche "Volumes Finis" ?

Diérences nis
Bien connue
Mise en ÷vre simple pour une géometrie simple
Mise en ÷vre dicile pour une géometrie complexe
Pas toujours conservative
Utilisation dans des codes de "recherche"

Y. EL KHCHINE Mod-Syst-Energ 4 / 57
Pourquoi utiliser l'approche "Volumes Finis" ?

Diérences nis
Bien connue
Mise en ÷vre simple pour une géometrie simple
Mise en ÷vre dicile pour une géometrie complexe
Pas toujours conservative
Utilisation dans des codes de "recherche"

Eléments nis
Approche très "mathématique"
S'adapte à une géometrie quelconque
Dicultées pour resoudre les termes non-lineaires
Très utilisée dans le domaine de Mécanique des Solides.

Y. EL KHCHINE Mod-Syst-Energ 4 / 57
Pourquoi utiliser l'approche "Volumes Finis" ?

Volumes nis
Approche très "physique" : bilan des ux
S'adapte à une géometrie quelconque
Plusieurs schémas pour la résolution des termes non-lineaires
hyperboliques
Conservative (par sa formulation)
La base de tout les codes généralistes en Mécanique des Fluides : Fluent
et CFX (ANSYS), StarCCM+, OpenFoam (Libre)...

Y. EL KHCHINE Mod-Syst-Energ 5 / 57
Méthode des volumes nis

Les diérentes étapes de la mise en ÷uvre sont :

Le domaine de calcul est discrétisé en un nombre ni de points (les n÷uds
du maillage), autour desquels on dénit des volumes élémentaires (appelés
volumes de contrôle) contigus, non juxtaposés et sans discontinuités aux
interfaces ;

Y. EL KHCHINE Mod-Syst-Energ 6 / 57
Méthode des volumes nis

Les diérentes étapes de la mise en ÷uvre sont :

Le domaine de calcul est discrétisé en un nombre ni de points (les n÷uds
du maillage), autour desquels on dénit des volumes élémentaires (appelés
volumes de contrôle) contigus, non juxtaposés et sans discontinuités aux
interfaces ;
Les équations modèles, sous forme conservative, sont intégrées sur chaque
volume de contrôle (VC) ;

Y. EL KHCHINE Mod-Syst-Energ 6 / 57
Méthode des volumes nis

Les diérentes étapes de la mise en ÷uvre sont :

Le domaine de calcul est discrétisé en un nombre ni de points (les n÷uds
du maillage), autour desquels on dénit des volumes élémentaires (appelés
volumes de contrôle) contigus, non juxtaposés et sans discontinuités aux
interfaces ;
Les équations modèles, sous forme conservative, sont intégrées sur chaque
volume de contrôle (VC) ;
Les intégrales sur un volume de contrôle en un n÷ud donné sont évaluées
en approchant la variation de Φ par des prols ou des lois d'interpolation
entre les n÷uds voisins du point considéré ;

Y. EL KHCHINE Mod-Syst-Energ 6 / 57
Méthode des volumes nis

Les diérentes étapes de la mise en ÷uvre sont :

Le domaine de calcul est discrétisé en un nombre ni de points (les n÷uds
du maillage), autour desquels on dénit des volumes élémentaires (appelés
volumes de contrôle) contigus, non juxtaposés et sans discontinuités aux
interfaces ;
Les équations modèles, sous forme conservative, sont intégrées sur chaque
volume de contrôle (VC) ;
Les intégrales sur un volume de contrôle en un n÷ud donné sont évaluées
en approchant la variation de Φ par des prols ou des lois d'interpolation
entre les n÷uds voisins du point considéré ;
Écriture des équations algébriques en fonction des valeurs de Φ aux n÷uds
du maillage ;

Y. EL KHCHINE Mod-Syst-Energ 6 / 57
Méthode des volumes nis

Les diérentes étapes de la mise en ÷uvre sont :

Le domaine de calcul est discrétisé en un nombre ni de points (les n÷uds
du maillage), autour desquels on dénit des volumes élémentaires (appelés
volumes de contrôle) contigus, non juxtaposés et sans discontinuités aux
interfaces ;
Les équations modèles, sous forme conservative, sont intégrées sur chaque
volume de contrôle (VC) ;
Les intégrales sur un volume de contrôle en un n÷ud donné sont évaluées
en approchant la variation de Φ par des prols ou des lois d'interpolation
entre les n÷uds voisins du point considéré ;
Écriture des équations algébriques en fonction des valeurs de Φ aux n÷uds
du maillage ;
Résolution du système algébrique linéaire obtenu.

Y. EL KHCHINE Mod-Syst-Energ 6 / 57
Méthode des volumes nis

Equations de conservation Pour un écoulement incompressible, on a

=0
 ∂vi


 ∂xi

  
∂vi ∂vi ∂p
ρ ∂t + vj ∂xj
= − ∂xi
+ ρgi + µ∆vi

   2

 ρcV ∂T

+ vi ∂T = k ∂x∂ i ∂x
T
+ϕ
∂t ∂xi i

En plus de ces équations, on ajoute les équations d'états des gaz,

ρ = ρ (p, T ) et p = p (ρ, T )

Y. EL KHCHINE Mod-Syst-Energ 7 / 57
Méthode des volumes nis

Forme diérentielle et intégrale de l'équation générale de transport

Les équations de conservation peuvent être réécrites sous la forme générale
suivante :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t
C'est l'équation de convection-diusion

Y. EL KHCHINE Mod-Syst-Energ 8 / 57
Méthode des volumes nis

Forme diérentielle et intégrale de l'équation générale de transport

Les équations de conservation peuvent être réécrites sous la forme générale
suivante :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t
C'est l'équation de convection-diusion
Φ est une variable caractéristique généralisée du système (vitesse,
température...).

Y. EL KHCHINE Mod-Syst-Energ 8 / 57
Méthode des volumes nis

Forme diérentielle et intégrale de l'équation générale de transport

Les équations de conservation peuvent être réécrites sous la forme générale
suivante :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t
C'est l'équation de convection-diusion
Φ est une variable caractéristique généralisée du système (vitesse,
température...). Γ est le coecient de diusion généralisé (propre au
phénomène traité). Il exprime la proporiété de transport comme la viscosité ou
la diusivité, conjointement avec le gradient de la variable Φ.

Y. EL KHCHINE Mod-Syst-Energ 8 / 57
Méthode des volumes nis

Forme diérentielle et intégrale de l'équation générale de transport

Les équations de conservation peuvent être réécrites sous la forme générale
suivante :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t
C'est l'équation de convection-diusion
Φ est une variable caractéristique généralisée du système (vitesse,
température...). Γ est le coecient de diusion généralisé (propre au
phénomène traité). Il exprime la proporiété de transport comme la viscosité ou
la diusivité, conjointement avec le gradient de la variable Φ. S est le terme
source spécique à chaque variable. Il exprime le processus de production
interne comme la génération de la chaleur dans le uide.

Y. EL KHCHINE Mod-Syst-Energ 8 / 57
Méthode des volumes nis
Équation diérentielle de convection diusion :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t

Variables
Équations Φ Γ S
Continuité 1 0 0
Conser. Quant. Mvt. x u µ −∂p
∂x + Sx
Conser. Quant. Mvt. y v µ −∂p
∂y + Sy
Conser. Quant. Mvt. z w µ −∂p
∂z + Sz
Conservation de l'énergie T λ
Cp Si

Y. EL KHCHINE Mod-Syst-Energ 9 / 57
Méthode des volumes nis
Équation diérentielle de convection diusion :
∂(ρΦ)
+ div (ρΦ⃗
u ) = div (ΓgradΦ) + SΦ
∂t

Variables
Équations Φ Γ S
Continuité 1 0 0
Conser. Quant. Mvt. x u µ −∂p
∂x + Sx
Conser. Quant. Mvt. y v µ −∂p
∂y + Sy
Conser. Quant. Mvt. z w µ −∂p
∂z + Sz
Conservation de l'énergie T λ
Cp Si

En intégrant sur le volume de contrôle 3D :

Z Z Z Z
∂
ρΦdV + div (ρΦ⃗
u )dV = div (ΓgradΦ)dV + SΦ dV
∂t VC VC VC VC
Y. EL KHCHINE Mod-Syst-Energ 9 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
L'équation de la diusion stationnaire dans le cas monodimensionnel s'écrit
div (ΓgradΦ) + SΦ = 0

Y. EL KHCHINE Mod-Syst-Energ 10 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
L'équation de la diusion stationnaire dans le cas monodimensionnel s'écrit
div (ΓgradΦ) + SΦ = 0
Cas 1D :  
d dΦ
Γ + SΦ = 0
dx dx

Y. EL KHCHINE Mod-Syst-Energ 10 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
L'équation de la diusion stationnaire dans le cas monodimensionnel s'écrit
div (ΓgradΦ) + SΦ = 0
Cas 1D :  
d dΦ
Γ + SΦ = 0
dx dx
Cas typique de la conduction de chaleur dans une barre métallique d'extrémité A et
B.

Y. EL KHCHINE Mod-Syst-Energ 10 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
L'équation de la diusion stationnaire dans le cas monodimensionnel s'écrit
div (ΓgradΦ) + SΦ = 0
Cas 1D :  
d dΦ
Γ + SΦ = 0
dx dx
Cas typique de la conduction de chaleur dans une barre métallique d'extrémité A et
B.
Domaine de calcul 1D : Le domaine d'étude est subdivisé en un nombre ni de
n÷uds. Chaque volume ni entoure un n÷ud principal "P". Les n÷uds voisins sont
"E" côté Est et "W" côté West.

Y. EL KHCHINE Mod-Syst-Energ 10 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Les conditions aux limites sont supposées connues aux extrémités de la tige.

Y. EL KHCHINE Mod-Syst-Energ 11 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Les conditions aux limites sont supposées connues aux extrémités de la tige.
□ Étape 1 : Génération du maillage
Notations utilisées pour le maillage :

Y. EL KHCHINE Mod-Syst-Energ 11 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

□ Étape 2 : Intégration de l'équation de diusion sur le VC

L'équation précédente est intégrée sur le volume de controle limité par les faces (w,e).
Pour un problème monodimensionnel, le volume d'un volume ni est ∆x × 1 × 1.

Y. EL KHCHINE Mod-Syst-Energ 12 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

□ Étape 2 : Intégration de l'équation de diusion sur le VC

L'équation précédente est intégrée sur le volume de controle limité par les faces (w,e).
Pour un problème monodimensionnel, le volume d'un volume ni est ∆x × 1 × 1.
Z e   Z e
d dΦ
⇔ Γ dx + SΦ dx = 0
w dx dx w

Y. EL KHCHINE Mod-Syst-Energ 12 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

□ Étape 2 : Intégration de l'équation de diusion sur le VC

L'équation précédente est intégrée sur le volume de controle limité par les faces (w,e).
Pour un problème monodimensionnel, le volume d'un volume ni est ∆x × 1 × 1.
Z e   Z e
d dΦ
⇔ Γ dx + SΦ dx = 0
w dx dx w

 e Z e
dΦ
Γ + SΦ dx = 0
dx w w

Y. EL KHCHINE Mod-Syst-Energ 12 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

□ Étape 2 : Intégration de l'équation de diusion sur le VC

L'équation précédente est intégrée sur le volume de controle limité par les faces (w,e).
Pour un problème monodimensionnel, le volume d'un volume ni est ∆x × 1 × 1.
Z e   Z e
d dΦ
⇔ Γ dx + SΦ dx = 0
w dx dx w

 e Z e
dΦ
Γ + SΦ dx = 0
dx w w
   
dΦ dΦ
⇒ Γ − Γ + S̄Φ ∆x = 0
dx e dx w

S̄Φ est la valeur moyenne du terme source sur le VC.

C'est la forme discrétisée de l'équation de diusion. Elle peut être écrite sur des n÷uds
P du domaine de calcul.

Y. EL KHCHINE Mod-Syst-Energ 12 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Signication physique :
Le bilan de ϕ sur le VC ⇒ le ux diusif de ϕ quittant la face EST du VC - le ux
diusif de ϕ entrant de la face WEST du VC = génération de ϕ dans le VC.
   
dΦ dΦ
Γ − Γ + S̄Φ ∆x = 0
dx e dx w

Y. EL KHCHINE Mod-Syst-Energ 13 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Signication physique :
Le bilan de ϕ sur le VC ⇒ le ux diusif de ϕ quittant la face EST du VC - le ux
diusif de ϕ entrant de la face WEST du VC = génération de ϕ dans le VC.
   
dΦ dΦ
Γ − Γ + S̄Φ ∆x = 0
dx e dx w

Calcul des gradients (ux diusif) sur les faces du VC :

Le prol linéaire de température, entre les n÷uds du maillage, dénie les dérivées aux
interfaces par :

Y. EL KHCHINE Mod-Syst-Energ 13 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Signication physique :
Le bilan de ϕ sur le VC ⇒ le ux diusif de ϕ quittant la face EST du VC - le ux
diusif de ϕ entrant de la face WEST du VC = génération de ϕ dans le VC.
   
dΦ dΦ
Γ − Γ + S̄Φ ∆x = 0
dx e dx w

Calcul des gradients (ux diusif) sur les faces du VC :

Le prol linéaire de température, entre les n÷uds du maillage, dénie les dérivées aux
interfaces par :
−ΦP E −ΦP
Γ dΦ = Γe ΦxEE −x = Γe Φδx
(

dx e P PE

Γ dΦ = Γw ΦxPP −Φ ΦP −ΦW


−xW = Γw δxPW
W
dx w

Y. EL KHCHINE Mod-Syst-Energ 13 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

Signication physique :
Le bilan de ϕ sur le VC ⇒ le ux diusif de ϕ quittant la face EST du VC - le ux
diusif de ϕ entrant de la face WEST du VC = génération de ϕ dans le VC.
   
dΦ dΦ
Γ − Γ + S̄Φ ∆x = 0
dx e dx w

Calcul des gradients (ux diusif) sur les faces du VC :

Le prol linéaire de température, entre les n÷uds du maillage, dénie les dérivées aux
interfaces par :
−ΦP E −ΦP
Γ dΦ = Γe ΦxEE −x = Γe Φδx
(

dx e P PE

Γ dΦ = Γw ΦxPP −Φ ΦP −ΦW


−xW = Γw δxPW
W
dx w
et on utilise une approximation linéaire pour calculer le terme Γ sur les faces du VC :

ΓP +ΓW
(
Γw = 2
ΓP +ΓE
Γe = 2
Y. EL KHCHINE Mod-Syst-Energ 13 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
Linéarisation du terme source :
Le terme source peut être constant ce qui ne pose aucun problème particulier. Mais souvent
il dépend explicitement de la variable Φp . Cela peut causer des problèmes au niveau du
traitement numérique (instabilités, divergence de la méthode numérique). Pour éviter ce
genre de problème il faut linéariser ce terme en le décomposant en deux parties :

S̄Φ = Su + Sp Φp
Où Su est la partie constante du terme source et Sp et le coecient de Φp .

Y. EL KHCHINE Mod-Syst-Energ 14 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
Linéarisation du terme source :
Le terme source peut être constant ce qui ne pose aucun problème particulier. Mais souvent
il dépend explicitement de la variable Φp . Cela peut causer des problèmes au niveau du
traitement numérique (instabilités, divergence de la méthode numérique). Pour éviter ce
genre de problème il faut linéariser ce terme en le décomposant en deux parties :

S̄Φ = Su + Sp Φp
Où Su est la partie constante du terme source et Sp et le coecient de Φp .
   
ΦE − ΦP ΦP − ΦW
Γe − Γw + Su + Sp ΦP = 0
δxPE δxPW

Y. EL KHCHINE Mod-Syst-Energ 14 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
Linéarisation du terme source :
Le terme source peut être constant ce qui ne pose aucun problème particulier. Mais souvent
il dépend explicitement de la variable Φp . Cela peut causer des problèmes au niveau du
traitement numérique (instabilités, divergence de la méthode numérique). Pour éviter ce
genre de problème il faut linéariser ce terme en le décomposant en deux parties :

S̄Φ = Su + Sp Φp
Où Su est la partie constante du terme source et Sp et le coecient de Φp .
   
ΦE − ΦP ΦP − ΦW
Γe − Γw + Su + Sp ΦP = 0
δxPE δxPW
     
Γe Γw Γw Γe
+ − SP ΦP = ΦW + Φ E + Su
δxPE δxPW δxPW δxPE
| {z } | {z } | {z }
aP aW aE

Y. EL KHCHINE Mod-Syst-Energ 14 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
Linéarisation du terme source :
Le terme source peut être constant ce qui ne pose aucun problème particulier. Mais souvent
il dépend explicitement de la variable Φp . Cela peut causer des problèmes au niveau du
traitement numérique (instabilités, divergence de la méthode numérique). Pour éviter ce
genre de problème il faut linéariser ce terme en le décomposant en deux parties :

S̄Φ = Su + Sp Φp
Où Su est la partie constante du terme source et Sp et le coecient de Φp .
   
ΦE − ΦP ΦP − ΦW
Γe − Γw + Su + Sp ΦP = 0
δxPE δxPW
     
Γe Γw Γw Γe
+ − SP ΦP = ΦW + Φ E + Su
δxPE δxPW δxPW δxPE
| {z } | {z } | {z }
aP aW aE

Ainsi l'équation algébrique nale s'écrit sous la forme :

=⇒ aP ΦP = aW ΦW + aE ΦE + Su

Y. EL KHCHINE Mod-Syst-Energ 14 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.

Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.
aP ΦP = aW ΦW + aE ΦE + Su

Cette équation s'écrit autant de fois qu'il y a de n÷uds où la variable Φ est inconnue.

Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.
aP ΦP = aW ΦW + aE ΦE + Su

Cette équation s'écrit autant de fois qu'il y a de n÷uds où la variable Φ est inconnue.
Dans le cas d'un maillage uniforme, ∆x = (δx)PW = (δx)PE .

Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.
aP ΦP = aW ΦW + aE ΦE + Su

Cette équation s'écrit autant de fois qu'il y a de n÷uds où la variable Φ est inconnue.
Dans le cas d'un maillage uniforme, ∆x = (δx)PW = (δx)PE .

alors on aura aW = Γw
∆x
, et aE = Γe
∆x

Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.
aP ΦP = aW ΦW + aE ΦE + Su

Cette équation s'écrit autant de fois qu'il y a de n÷uds où la variable Φ est inconnue.
Dans le cas d'un maillage uniforme, ∆x = (δx)PW = (δx)PE .

alors on aura aW = Γw
∆x
, et aE = Γe
∆x

L'écriture de l'équation pour tous les n÷uds du domaine conduit à un système algébrique
qui s'écrit sous la forme matricielle suivante :

[M]{ϕ} = {b}

Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D
□ Étape 3 : Résolution du système algébrique
La valeur de Φp au centre du VC est exprimée en fonction des valeurs de la grandeur aux
n÷uds voisins à l'aide de l'équation discrète.
aP ΦP = aW ΦW + aE ΦE + Su

Cette équation s'écrit autant de fois qu'il y a de n÷uds où la variable Φ est inconnue.
Dans le cas d'un maillage uniforme, ∆x = (δx)PW = (δx)PE .

alors on aura aW = Γw
∆x
, et aE = Γe
∆x

L'écriture de l'équation pour tous les n÷uds du domaine conduit à un système algébrique
qui s'écrit sous la forme matricielle suivante :

[M]{ϕ} = {b}

Méthodes de résolution du système algébrique :

Gauss-Seidel,
Tri-Diagonal Matrix Algorithm (TDMA)
Y. EL KHCHINE Mod-Syst-Energ 15 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Méthode de Gauss-Seidel :
Une méthode itérative de résolution d'un système linéaire (de dimension nie)
de la forme [M]{ϕ} = {b}. On calcule les composantes Φk+ i de Φk+ de
1 1

manière séquentielle pour i = 1, . . . , n par :

Y. EL KHCHINE Mod-Syst-Energ 16 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Méthode de Gauss-Seidel :
Une méthode itérative de résolution d'un système linéaire (de dimension nie)
de la forme [M]{ϕ} = {b}. On calcule les composantes Φk+ i de Φk+ de
1 1

manière séquentielle pour i = 1, . . . , n par :

 
1  i−1
X n
X
Φk+
i
1
= bi − mij Φk+
j
1
− mij Φkj 
mii
j=1 j=i+1

Y. EL KHCHINE Mod-Syst-Energ 16 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Méthode de Gauss-Seidel :
Une méthode itérative de résolution d'un système linéaire (de dimension nie)
de la forme [M]{ϕ} = {b}. On calcule les composantes Φk+ i de Φk+ de 1 1

manière séquentielle pour i = 1, . . . , n par :

 
1  i−1
X n
X
Φk+
i
1
= bi − mij Φk+
j
1
− mij Φkj 
mii
j=1 j=i+1

En donnent aux inconnues ϕki des valeurs arbitraires initiales ϕi (pour k = 0)

0

Le processus sera arrêté si : |ϕik+ − ϕki |< ε, avec ε est la précision absolue
1

imposée par l'utilisateur.

Y. EL KHCHINE Mod-Syst-Energ 16 / 57
 Volumes nis pour les problèmes de diusion stationnaire 1D
• Tri-diagonal matrix algorithm (TDMA)
Cet algorithme permet de calculer la solution d'un système linéaire lorsque la matrice est tri-diagonale.

C'est notre cas, car les équations discrétisées dans le système linéaire s'écrivent :

Y. EL KHCHINE Mod-Syst-Energ 17 / 57
 Volumes nis pour les problèmes de diusion stationnaire 1D
• Tri-diagonal matrix algorithm (TDMA)
Cet algorithme permet de calculer la solution d'un système linéaire lorsque la matrice est tri-diagonale.

C'est notre cas, car les équations discrétisées dans le système linéaire s'écrivent :

ai Φi = bi Φi+1 + ci Φi−1 + di 2 ≤i ≤N

Y. EL KHCHINE Mod-Syst-Energ 17 / 57
 Volumes nis pour les problèmes de diusion stationnaire 1D
• Tri-diagonal matrix algorithm (TDMA)
Cet algorithme permet de calculer la solution d'un système linéaire lorsque la matrice est tri-diagonale.

C'est notre cas, car les équations discrétisées dans le système linéaire s'écrivent :

ai Φi = bi Φi+1 + ci Φi−1 + di 2 ≤i ≤N
où le maillage est présenté comme suit :

Y. EL KHCHINE Mod-Syst-Energ 17 / 57
 Volumes nis pour les problèmes de diusion stationnaire 1D
• Tri-diagonal matrix algorithm (TDMA)
Cet algorithme permet de calculer la solution d'un système linéaire lorsque la matrice est tri-diagonale.

C'est notre cas, car les équations discrétisées dans le système linéaire s'écrivent :

ai Φi = bi Φi+1 + ci Φi−1 + di 2 ≤i ≤N
où le maillage est présenté comme suit :

Soit un système d'équation ayant la forme tri-diagonale :

Φ1 = C1 (1)

−β2 Φ1 + D2 Φ2 − α2 Φ3 = C2 (2)

−β3 Φ2 + D3 Φ3 − α3 Φ4 = C3 (3)

−β4 Φ3 + D4 Φ4 − α3 Φ4 = C4 (4)

··· ··· ··· ··· =· (5)

−βN ΦN−1 + DN ΦN − αN ΦN+1 = CN (6)

ΦN+1 = CN+1 (7)

Φ1 et ΦN+1 sont donnés par les conditions aux limites.

Y. EL KHCHINE Mod-Syst-Energ 17 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en amont (forward elimination)

Ce processus commence par remplacer Φ2 obtenue par l'équation (2) dans l'équation (3) :

Y. EL KHCHINE Mod-Syst-Energ 18 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en amont (forward elimination)

Ce processus commence par remplacer Φ2 obtenue par l'équation (2) dans l'équation (3) :
  
β2

! C2
α3 β3 Φ
D2 1
+D2
+ C3
Φ3 = 2
Φ4 +  2

D 3 − β3 α
D2
D3 − β3 α
D2

Y. EL KHCHINE Mod-Syst-Energ 18 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en amont (forward elimination)

Ce processus commence par remplacer Φ2 obtenue par l'équation (2) dans l'équation (3) :
  
β2

! C2
α3 β3 Φ
D2 1
+D2
+ C3
Φ3 = 2
Φ4 +  2

D 3 − β3 α
D2
D3 − β3 α
D2

Supposons que :
α2 β C
A2 = et C2′ = 2 Φ1 + 2
D2 D2 D2

Y. EL KHCHINE Mod-Syst-Energ 18 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en amont (forward elimination)

Ce processus commence par remplacer Φ2 obtenue par l'équation (2) dans l'équation (3) :
  
β2

! C2
α3 β3 Φ
D2 1
+D2
+ C3
Φ3 = 2
Φ4 +  2

D 3 − β3 α
D2
D3 − β3 α
D2

Supposons que :
α2 β C
A2 = et C2′ = 2 Φ1 + 2
D2 D2 D2
On aura :
β3 C2′ + C3
   
α3
Φ3 = Φ4 +
D3 − β3 A2 D3 − β3 A2

Φ3 = A3 Φ4 + C3′

Y. EL KHCHINE Mod-Syst-Energ 18 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en amont (forward elimination)

Ce processus commence par remplacer Φ2 obtenue par l'équation (2) dans l'équation (3) :
  
β2

! C2
α3 β3 Φ
D2 1
+D2
+ C3
Φ3 = 2
Φ4 +  2

D 3 − β3 α
D2
D3 − β3 α
D2

Supposons que :
α2 β C
A2 = et C2′ = 2 Φ1 + 2
D2 D2 D2
On aura :
β3 C2′ + C3
   
α3
Φ3 = Φ4 +
D3 − β3 A2 D3 − β3 A2

Φ3 = A3 Φ4 + C3′
Cette équation est utilisée à son tour dans l'équation (4). Cette procédure est répété
jusqu'à la dernière équation.

Y. EL KHCHINE Mod-Syst-Energ 18 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

•Tri-diagonal matrix algorithm (TDMA)

□ Élimination en arrière (backward elimination)
On utilise la forme généralisé de l'équation précédente :

Φj = Aj Φj+1 + Cj′

avec : !
  ′
βj Cj− + Cj
αj
Aj = et Cj′ =
1

Dj − βj Aj−1 Dj − βj Aj−1
sur les frontières :
j =1 on a A1 = 0 et C ′ = Φ
1 1

j = N + 1 on a AN+1 = 0 et CN+ ′
= ΦN+
1 1

Y. EL KHCHINE Mod-Syst-Energ 19 / 57
Volumes nis pour les problèmes de diusion stationnaire 1D

• Tri-diagonal matrix algorithm (TDMA)

□ Élimination en arrière (backward elimination)

Y. EL KHCHINE Mod-Syst-Energ 20 / 57
Application 1 :

On considère une barre cylindrique, sans source de chaleur, ayant l'aire

transversale A et la longueur L = 0.5 m. Les extrémités, A et B de la barre
sont maintenues aux températures constantes de 100◦ C et de 500◦ C
respectivement.
Calculer la distribution de la température le long de la barre. On connait la
conductivité thermique λ = 1000 W/m.K.

Utiliser 6 n÷uds le long de la barre en plus des conditions aux limites.

Y. EL KHCHINE Mod-Syst-Energ 21 / 57
Application 1 :
La distribution de la température est gouvernée par l'équation : d
dx = 0
λ dT


dx
On considère six points le long de la barre avec ∆x = 0.1m.

Si l'on considère que λe = λw = λ pour les noeuds intérieurs (3 et 4),

l'équation discritisée est :
aP TP = aW TW + aE TE où aP = aE + aW et aW = aE = ∆x λ
.
Pour les noeuds 2 et 5 on utilise la même équation que pour un noeud
intérieur mais on tient compte que pour les noeuds voisins 1 et 6 les
températures sont connues, TA et TB respectivement.
Y. EL KHCHINE Mod-Syst-Energ 22 / 57
Application 1 :
Pour le noeud 2, l'équation discrétisée est la suivante : aP TP = aE TE + aW TA
Pour le noeud 5, l'équation discrétisée est la suivante :
aP TP = aW TW + aE TB
Les équations à résoudre sont les suivantes (le nombre d'équations est égal à
4) :
ai Ti = ai− Ti− + ai+ Ti+ ; i = 2 → 5
1 1 1 1

En tenant compte que

∆x = , = 10000; aP = aW + aE = 10000 + 10000 = 20000
λ 1000
0 1

Le système d'équations à résoudre est :

2T = T + TA
2 3

2T = T + T
3 2 4

2T = T + T
4 3 5

2T = T + TB
5 4

Y. EL KHCHINE Mod-Syst-Energ 23 / 57
Application 1 :

Y. EL KHCHINE Mod-Syst-Energ 24 / 57
Application 1 :

Y. EL KHCHINE Mod-Syst-Energ 25 / 57
Application 2 :

On considère une plaque très longue d'épaisseur L = 2 cm, ayant la conductivité

thermique constante λ = 0.5 W/m.K et une source de chaleur uniforme, S =
1000 kW/m . Les faces de la plaque se trouvent à la température constante de
3

100 ◦ C et 200 ◦ C respectivement.

En supposant que les dimensions de la plaque dans les directions "y" et "z"
soient très grandes et donc le gradient de la température est signicatif dans la
direction "x" seulement, calculer la distribution de la température et comparer
les résultats numériques avec la solution analytique :
 
TB − TA S
T (x) = + (L − x) x + TA
L 2λ

Y. EL KHCHINE Mod-Syst-Energ 26 / 57
Application 2 :

Y. EL KHCHINE Mod-Syst-Energ 27 / 57
Application 2 :

Y. EL KHCHINE Mod-Syst-Energ 28 / 57
Application 2 :

Y. EL KHCHINE Mod-Syst-Energ 29 / 57
Application 2 :

Y. EL KHCHINE Mod-Syst-Energ 30 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D

L'équation de la diusion stationnaire dans le cas 2D est :

   
∂ ∂ϕ ∂ ∂ϕ
Γ + Γ + Sϕ = 0
∂x ∂x ∂y ∂y

Le discétisation de domaine 2D

Y. EL KHCHINE Mod-Syst-Energ 31 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D

L'intégrale de l'équation de la diusion stationnaire sur le volume de contrôle

donne

Y. EL KHCHINE Mod-Syst-Energ 32 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D

L'intégrale de l'équation de la diusion stationnaire sur le volume de contrôle

donne
Z   Z   Z
∂ ∂Φ ∂ ∂Φ
Γ dx · dy + Γ dx · dy + SΦ dV = 0
VC ∂x ∂x VC ∂y ∂y VC

Y. EL KHCHINE Mod-Syst-Energ 32 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D

L'intégrale de l'équation de la diusion stationnaire sur le volume de contrôle

donne
Z   Z   Z
∂ ∂Φ ∂ ∂Φ
Γ dx · dy + Γ dx · dy + SΦ dV = 0
VC ∂x ∂x VC ∂y ∂y VC

Cela conduit à
     
∂Φ ∂Φ
Γe ∆y − Γw ∆y
∂x e ∂x w
    
∂Φ ∂Φ
+ Γn ∆x − Γs ∆x + S̄Φ ∆V = 0
∂y n ∂y s

Y. EL KHCHINE Mod-Syst-Energ 32 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

∂Φ (ΦE − ΦP )
Γe = Γe
∂x e δxPE

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

∂Φ (ΦE − ΦP )
Γe = Γe
∂x e δxPE

Flux à travers la face Sud

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

∂Φ (ΦE − ΦP )
Γe = Γe
∂x e δxPE

Flux à travers la face Sud

∂Φ (ΦP − ΦS )
Γs = Γs
∂y s δySP

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

∂Φ (ΦE − ΦP )
Γe = Γe
∂x e δxPE

Flux à travers la face Sud

∂Φ (ΦP − ΦS )
Γs = Γs
∂y s δySP

Flux à travers la face Nord

Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Calcul des gradients (ux diusif) sur les faces du VC
Flux à travers la face West
∂ϕ (ΦP − ΦW )
Γw = Γw
∂x w δxWP

Flux à travers la face Est

∂Φ (ΦE − ΦP )
Γe = Γe
∂x e δxPE

Flux à travers la face Sud

∂Φ (ΦP − ΦS )
Γs = Γs
∂y s δySP

Flux à travers la face Nord

∂Φ (ΦN − ΦP )
Γn = Γn
∂y n δyPN
Y. EL KHCHINE Mod-Syst-Energ 33 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Système algébrique
En remplaçons les ux dans l'équation :
     
∂Φ ∂Φ
Γe ∆y − Γw ∆y
∂x e ∂x w
    
∂Φ ∂Φ
+ Γn ∆x − Γs ∆x + S̄Φ ∆V = 0
∂y n ∂y s

Y. EL KHCHINE Mod-Syst-Energ 34 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Système algébrique
En remplaçons les ux dans l'équation :
     
∂Φ ∂Φ
Γe ∆y − Γw ∆y
∂x e ∂x w
    
∂Φ ∂Φ
+ Γn ∆x − Γs ∆x + S̄Φ ∆V = 0
∂y n ∂y s
Et en linéarisant le terme source, on obtient le système algébrique suivant :

Y. EL KHCHINE Mod-Syst-Energ 34 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Système algébrique
En remplaçons les ux dans l'équation :
     
∂Φ ∂Φ
Γe ∆y − Γw ∆y
∂x e ∂x w
    
∂Φ ∂Φ
+ Γn ∆x − Γs ∆x + S̄Φ ∆V = 0
∂y n ∂y s
Et en linéarisant le terme source, on obtient le système algébrique suivant :

aP Φp = aW ΦW + aE ΦE + aN ΦN + aS ΦS + Su ∆x∆y

Y. EL KHCHINE Mod-Syst-Energ 34 / 57
Volumes nis pour les problèmes de diusion stationnaire 2D
- Système algébrique
En remplaçons les ux dans l'équation :
     
∂Φ ∂Φ
Γe ∆y − Γw ∆y
∂x e ∂x w
    
∂Φ ∂Φ
+ Γn ∆x − Γs ∆x + S̄Φ ∆V = 0
∂y n ∂y s
Et en linéarisant le terme source, on obtient le système algébrique suivant :

aP Φp = aW ΦW + aE ΦE + aN ΦN + aS ΦS + Su ∆x∆y
avec :
Γw ∆y Γe ∆y Γs ∆x Γn ∆x
aW = δxWP ; aE = δxPE ; aS = δySP ; aN = δyPN
ap = aW + aE + aS + aN − Sp ∆x∆y
Y. EL KHCHINE Mod-Syst-Energ 34 / 57
Volumes nis pour les problèmes de convection-diusion

Dans les problèmes où le mouvement du uide est mis en évidence et joue un

rôle signicatif, il faut tenir en compte l'eet de convection. Ceci est pris en
considération dans l'équation de convection-diusion (stationnaire) :

Y. EL KHCHINE Mod-Syst-Energ 35 / 57
Volumes nis pour les problèmes de convection-diusion

Dans les problèmes où le mouvement du uide est mis en évidence et joue un

rôle signicatif, il faut tenir en compte l'eet de convection. Ceci est pris en
considération dans l'équation de convection-diusion (stationnaire) :

div(ρΦ→
−
u ) = div(Γ grad Φ) + SΦ

Y. EL KHCHINE Mod-Syst-Energ 35 / 57
Volumes nis pour les problèmes de convection-diusion

Dans les problèmes où le mouvement du uide est mis en évidence et joue un

rôle signicatif, il faut tenir en compte l'eet de convection. Ceci est pris en
considération dans l'équation de convection-diusion (stationnaire) :

div(ρΦ→
−
u ) = div(Γ grad Φ) + SΦ

En intégrant sur le VC de volume ∆V nous obtenons le bilan du ux de la

propriété Φ à l'intérieure du VC :

Y. EL KHCHINE Mod-Syst-Energ 35 / 57
Volumes nis pour les problèmes de convection-diusion
Pour la diusion, la propagation de la propriété Φ s'eectue
équitablement dans toutes les directions,

Y. EL KHCHINE Mod-Syst-Energ 36 / 57
Volumes nis pour les problèmes de convection-diusion
Pour la diusion, la propagation de la propriété Φ s'eectue
équitablement dans toutes les directions,
Pour la convection, la propagation de la propriété Φ s'eectue seulement
suivant le sens de l'écoulement du uide ; elle dépend fortement du sens
et vitesse de l'écoulement.

Y. EL KHCHINE Mod-Syst-Energ 36 / 57
Volumes nis pour les problèmes de convection-diusion
Pour la diusion, la propagation de la propriété Φ s'eectue
équitablement dans toutes les directions,
Pour la convection, la propagation de la propriété Φ s'eectue seulement
suivant le sens de l'écoulement du uide ; elle dépend fortement du sens
et vitesse de l'écoulement.
□ Discrétisation de l'équation de convection-diusion

Y. EL KHCHINE Mod-Syst-Energ 36 / 57
Volumes nis pour les problèmes de convection-diusion
Pour la diusion, la propagation de la propriété Φ s'eectue
équitablement dans toutes les directions,
Pour la convection, la propagation de la propriété Φ s'eectue seulement
suivant le sens de l'écoulement du uide ; elle dépend fortement du sens
et vitesse de l'écoulement.
□ Discrétisation de l'équation de convection-diusion
Le terme diusif : déjà discrétisé et montre une grande stabilité.
Le terme convectif : des dicultés apparaissent lors de l'approximation
de la propriété Φ sur les faces de chaque VC ainsi que son ux convectif.

Y. EL KHCHINE Mod-Syst-Energ 36 / 57
Volumes nis pour les problèmes de convection-diusion
Pour la diusion, la propagation de la propriété Φ s'eectue
équitablement dans toutes les directions,
Pour la convection, la propagation de la propriété Φ s'eectue seulement
suivant le sens de l'écoulement du uide ; elle dépend fortement du sens
et vitesse de l'écoulement.
□ Discrétisation de l'équation de convection-diusion
Le terme diusif : déjà discrétisé et montre une grande stabilité.
Le terme convectif : des dicultés apparaissent lors de l'approximation
de la propriété Φ sur les faces de chaque VC ainsi que son ux convectif.
C'est pourquoi, un grand nombre de schémas, se diérenciant par leur précision
sont inventées pour traiter ce problème. Parmi ces schémas, on peut citer :
Schéma Centré
Schéma décentré (Upwind)
Schéma exponentiel
Schéma hybride
Y. EL KHCHINE Mod-Syst-Energ 36 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques
Application sur l'équation de convection-diusion 1-D :
 
d d dΦ
(ρΦu) = Γ
dx dx dx

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques
Application sur l'équation de convection-diusion 1-D :
 
d d dΦ
(ρΦu) = Γ
dx dx dx
Le uide doit aussi satisfaire l'équation de continuité :
d
(ρu) = 0
dx

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques
Application sur l'équation de convection-diusion 1-D :
 
d d dΦ
(ρΦu) = Γ
dx dx dx
Le uide doit aussi satisfaire l'équation de continuité :
d
(ρu) = 0
dx
En intégrant sur le VC on obtient :
(ρuΦ)e − (ρuΦ)w = Γ dΦ − Γ dΦ



dx e dx w
(ρu)e − (ρu)w = 0

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques
Application sur l'équation de convection-diusion 1-D :
 
d d dΦ
(ρΦu) = Γ
dx dx dx
Le uide doit aussi satisfaire l'équation de continuité :
d
(ρu) = 0
dx
En intégrant sur le VC on obtient :
(ρuΦ)e − (ρuΦ)w = Γ dΦ − Γ dΦ



dx e dx w
(ρu)e − (ρu)w = 0
soit : F = ρu et D = Γ
δx

donc sur les faces du VC on a :

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
□ Schémas numériques
Application sur l'équation de convection-diusion 1-D :
 
d d dΦ
(ρΦu) = Γ
dx dx dx
Le uide doit aussi satisfaire l'équation de continuité :
d
(ρu) = 0
dx
En intégrant sur le VC on obtient :
(ρuΦ)e − (ρuΦ)w = Γ dΦ − Γ dΦ



dx e dx w
(ρu)e − (ρu)w = 0
soit : F = ρu et D = Γ
δx

donc sur les faces du VC on a :

(
Fw = (ρu)w ; Fe = (ρu)e
Γw Γe
Dw = δxWP ; De = δxPE

Y. EL KHCHINE Mod-Syst-Energ 37 / 57
Volumes nis pour les problèmes de convection-diusion
En utilisant le schéma centré pour représenter le terme de diusion on
obtient :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (Φp − ΦW )
Fe − Fw = 0

Y. EL KHCHINE Mod-Syst-Energ 38 / 57
Volumes nis pour les problèmes de convection-diusion
En utilisant le schéma centré pour représenter le terme de diusion on
obtient :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (Φp − ΦW )
Fe − Fw = 0

1. Schéma Centré

Y. EL KHCHINE Mod-Syst-Energ 38 / 57
Volumes nis pour les problèmes de convection-diusion
En utilisant le schéma centré pour représenter le terme de diusion on
obtient :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (Φp − ΦW )
Fe − Fw = 0

1. Schéma Centré
Ce schéma a été utilisé pour représenter le terme diusif.

Y. EL KHCHINE Mod-Syst-Energ 38 / 57
Volumes nis pour les problèmes de convection-diusion
En utilisant le schéma centré pour représenter le terme de diusion on
obtient :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (Φp − ΦW )
Fe − Fw = 0

1. Schéma Centré
Ce schéma a été utilisé pour représenter le terme diusif.
Dans ce schéma, on exprime la variable Φ sur la face du volume de contrôle
par une approximation linéaire comme suit :
ΦE + ΦP ΦW + Φ P
Φe = et Φw =
2 2

Y. EL KHCHINE Mod-Syst-Energ 38 / 57
Volumes nis pour les problèmes de convection-diusion
En utilisant le schéma centré pour représenter le terme de diusion on
obtient :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (Φp − ΦW )
Fe − Fw = 0

1. Schéma Centré
Ce schéma a été utilisé pour représenter le terme diusif.
Dans ce schéma, on exprime la variable Φ sur la face du volume de contrôle
par une approximation linéaire comme suit :
ΦE + ΦP ΦW + Φ P
Φe = et Φw =
2 2
En remplaçant dans l'équation :
Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

Y. EL KHCHINE Mod-Syst-Energ 38 / 57
Volumes nis pour les problèmes de convection-diusion
On obtient
Fw Fe Fw Fe





Dw + 2 + De − 2 + (Fe − Fw ) ΦP = Dw + 2 Φw + D e − 2 ΦE

Y. EL KHCHINE Mod-Syst-Energ 39 / 57
Volumes nis pour les problèmes de convection-diusion
On obtient
Fw Fe Fw Fe





Dw + 2 + De − 2 + (Fe − Fw ) ΦP = Dw + 2 Φw + D e − 2 ΦE

=⇒ aP ΦP = aW ΦW + aE ΦE
□ Quelques propriétés des schémas numériques

Y. EL KHCHINE Mod-Syst-Energ 39 / 57
Volumes nis pour les problèmes de convection-diusion
On obtient
Fw Fe Fw Fe





Dw + 2 + De − 2 + (Fe − Fw ) ΦP = Dw + 2 Φw + D e − 2 ΦE

=⇒ aP ΦP = aW ΦW + aE ΦE
□ Quelques propriétés des schémas numériques
Théoriquement, la solution numérique est indépendante du schéma de
discrétisation utilisé quand le nombre de n÷uds considérés tend vers l'inni,
alors qu'en pratique ce nombre est toujours ni. Ceci nous conduit à une
étude plus détaillée des schémas de convection (puisque nous avons montré
précédemment que pour la diusion le schéma centré marche bien).
Pratiquement, un bon schéma de convection doit vérier les trois propriétés
suivantes :

Y. EL KHCHINE Mod-Syst-Energ 39 / 57
Volumes nis pour les problèmes de convection-diusion
On obtient
Fw Fe Fw Fe





Dw + 2 + De − 2 + (Fe − Fw ) ΦP = Dw + 2 Φw + D e − 2 ΦE

=⇒ aP ΦP = aW ΦW + aE ΦE
□ Quelques propriétés des schémas numériques
Théoriquement, la solution numérique est indépendante du schéma de
discrétisation utilisé quand le nombre de n÷uds considérés tend vers l'inni,
alors qu'en pratique ce nombre est toujours ni. Ceci nous conduit à une
étude plus détaillée des schémas de convection (puisque nous avons montré
précédemment que pour la diusion le schéma centré marche bien).
Pratiquement, un bon schéma de convection doit vérier les trois propriétés
suivantes :
Propriété 1 : Être conservative
Cette propriété vérie la conservation de Φ sur les faces du VC :
Y. EL KHCHINE Mod-Syst-Energ 39 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC

Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC
Cette propriété est généralement satisfaite par la formulation de la
méthode des volumes nis elle-même.

Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC
Cette propriété est généralement satisfaite par la formulation de la
méthode des volumes nis elle-même.
Le ux Φ doit être exprimer par une et même expression sur l'ensemble
du domaine du calcul.
Propriété 2 : Être bornée

Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC
Cette propriété est généralement satisfaite par la formulation de la
méthode des volumes nis elle-même.
Le ux Φ doit être exprimer par une et même expression sur l'ensemble
du domaine du calcul.
Propriété 2 : Être bornée
Critère de Scarborough :
≤ 1 sur tous les noeuds
P 
|avs |
|ap | < 1 au moins sur noeuds

Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC
Cette propriété est généralement satisfaite par la formulation de la
méthode des volumes nis elle-même.
Le ux Φ doit être exprimer par une et même expression sur l'ensemble
du domaine du calcul.
Propriété 2 : Être bornée
Critère de Scarborough :
≤ 1 sur tous les noeuds
P 
|avs |
|ap | < 1 au moins sur noeuds

Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Le ux de Φ Le ux de Φ
entrant à travers = quittant la face
une face du VC du même VC
Cette propriété est généralement satisfaite par la formulation de la
méthode des volumes nis elle-même.
Le ux Φ doit être exprimer par une et même expression sur l'ensemble
du domaine du calcul.
Propriété 2 : Être bornée
Critère de Scarborough :
≤ 1 sur tous les noeuds
P 
|avs |
|ap | < 1 au moins sur noeuds

Il s'agit d'une condition susante pour que le processus itérative utilisé

pour résoudre le système algébrique converge vers la solution et que ce
critère soit vérié.
Y. EL KHCHINE Mod-Syst-Energ 40 / 57
Volumes nis pour les problèmes de convection-diusion

Par exemple la valeur calculé de Φ sur une face du volume de contrôle doit être
bornée entre les valeurs de Φ aux limites. En général cette propriété est vériée
si tous les coecients son positifs.

Y. EL KHCHINE Mod-Syst-Energ 41 / 57
Volumes nis pour les problèmes de convection-diusion

Par exemple la valeur calculé de Φ sur une face du volume de contrôle doit être
bornée entre les valeurs de Φ aux limites. En général cette propriété est vériée
si tous les coecients son positifs.
Un schéma non borné, sa solution numérique ne converge pas en tous n÷uds
du maillage, ceci se manifeste par l'apparition des oscillations autour de la
solution de l'équation.

Y. EL KHCHINE Mod-Syst-Energ 41 / 57
Volumes nis pour les problèmes de convection-diusion

Par exemple la valeur calculé de Φ sur une face du volume de contrôle doit être
bornée entre les valeurs de Φ aux limites. En général cette propriété est vériée
si tous les coecients son positifs.
Un schéma non borné, sa solution numérique ne converge pas en tous n÷uds
du maillage, ceci se manifeste par l'apparition des oscillations autour de la
solution de l'équation.
Propriété 3 : Être transportée
Cette propriété est liée à la nature du transport qui est généralement traduit par le
nombre de Peclet :
Pe = ρu
Γ/δx

Le nombre de Peclet répresente le rapport du transfert par convection forcée et du

transport par diusion.

Y. EL KHCHINE Mod-Syst-Energ 41 / 57
Volumes nis pour les problèmes de convection-diusion

Par exemple la valeur calculé de Φ sur une face du volume de contrôle doit être
bornée entre les valeurs de Φ aux limites. En général cette propriété est vériée
si tous les coecients son positifs.
Un schéma non borné, sa solution numérique ne converge pas en tous n÷uds
du maillage, ceci se manifeste par l'apparition des oscillations autour de la
solution de l'équation.
Propriété 3 : Être transportée
Cette propriété est liée à la nature du transport qui est généralement traduit par le
nombre de Peclet :
Pe = ρu
Γ/δx

Le nombre de Peclet répresente le rapport du transfert par convection forcée et du

transport par diusion.
Les faibles valeurs du nombre de Peclet reètent un transport dominé par la
diusion alors que l'inverse caractérise un transport dominé par la convection.

Y. EL KHCHINE Mod-Syst-Energ 41 / 57
Volumes nis pour les problèmes de convection-diusion

Par exemple la valeur calculé de Φ sur une face du volume de contrôle doit être
bornée entre les valeurs de Φ aux limites. En général cette propriété est vériée
si tous les coecients son positifs.
Un schéma non borné, sa solution numérique ne converge pas en tous n÷uds
du maillage, ceci se manifeste par l'apparition des oscillations autour de la
solution de l'équation.
Propriété 3 : Être transportée
Cette propriété est liée à la nature du transport qui est généralement traduit par le
nombre de Peclet :
Pe = ρu
Γ/δx

Le nombre de Peclet répresente le rapport du transfert par convection forcée et du

transport par diusion.
Les faibles valeurs du nombre de Peclet reètent un transport dominé par la
diusion alors que l'inverse caractérise un transport dominé par la convection.
Un bon schéma de convection doit reéter cette propriété.
Y. EL KHCHINE Mod-Syst-Energ 41 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :

Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :
Φ = 1 à x = 0 ; ΦL = 0 à x = L
0

Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :
Φ = 1 à x = 0 ; ΦL = 0 à x = L
0

En utilisant un maillage à 6 n÷uds et le schéma centré (diérences centrales),

calculer la distribution de Φ en fonction de x pour les cas suivants :

Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :
Φ = 1 à x = 0 ; ΦL = 0 à x = L
0

En utilisant un maillage à 6 n÷uds et le schéma centré (diérences centrales),

calculer la distribution de Φ en fonction de x pour les cas suivants :
1) Cas 1 : u = 0.1 m/s et comparer avec la solution analytique ;

Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :
Φ = 1 à x = 0 ; ΦL = 0 à x = L
0

En utilisant un maillage à 6 n÷uds et le schéma centré (diérences centrales),

calculer la distribution de Φ en fonction de x pour les cas suivants :
1) Cas 1 : u = 0.1 m/s et comparer avec la solution analytique ;
2) Cas 2 : u = 2.5 m/s et comparer avec la solution analytique ;

Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :
Une propriété Φ (par exemple la température) est transportée par convection et
diusion à travers le domaine unidimensionnel schématisé à la gure ci-dessous.
Les conditions aux limites sont les suivantes :
Φ = 1 à x = 0 ; ΦL = 0 à x = L
0

En utilisant un maillage à 6 n÷uds et le schéma centré (diérences centrales),

calculer la distribution de Φ en fonction de x pour les cas suivants :
1) Cas 1 : u = 0.1 m/s et comparer avec la solution analytique ;
2) Cas 2 : u = 2.5 m/s et comparer avec la solution analytique ;

La solution analytique est donnée par :

Φ − Φ0 exp(ρux/Γ) − 1
=
ΦL − Φ0 exp(ρuL/Γ) − 1
Y. EL KHCHINE Mod-Syst-Energ 42 / 57
Application 3 :

3) Cas 3 : Recalculer la solution numérique pour u = 2.5 m/s en utilisant 21

n÷uds et comparer avec la solution analytique.
On donne : L = 1 m, ρ = 1 kg/m , Γ = 0.1 kg/m.s.
3

Y. EL KHCHINE Mod-Syst-Energ 43 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

L'équation qui gouverne le transport de la propriété ϕ est :

Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

L'équation qui gouverne le transport de la propriété ϕ est :

 
d d dϕ
(ρuϕ) = Γ
dx dx dx

Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

L'équation qui gouverne le transport de la propriété ϕ est :

 
d d dϕ
(ρuϕ) = Γ
dx dx dx
Pour un noeud intérieur (les noeuds 3 et 4), l'équation discrétisée peut être
écrite sous la forme suivante :

Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

L'équation qui gouverne le transport de la propriété ϕ est :

 
d d dϕ
(ρuϕ) = Γ
dx dx dx
Pour un noeud intérieur (les noeuds 3 et 4), l'équation discrétisée peut être
écrite sous la forme suivante :
aP ϕP − aW ϕW − aE ϕE = 0

Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :
On utilise le maillage présenté à la gure qui donne δx = 0,2 m. On rappel
que :
Γ
F = ρu, D = δx , Fe = Fw = F , De = Dw = D

L'équation qui gouverne le transport de la propriété ϕ est :

 
d d dϕ
(ρuϕ) = Γ
dx dx dx
Pour un noeud intérieur (les noeuds 3 et 4), l'équation discrétisée peut être
écrite sous la forme suivante :
aP ϕP − aW ϕW − aE ϕE = 0
Fw Fe
aW = D w + 2
, aE = De + 2
, aP = aW + aE + (Fe − Fw )
Y. EL KHCHINE Mod-Syst-Energ 44 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ − 0, 45ϕ = 0, 55ϕA
2 3

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ − 0, 45ϕ = 0, 55ϕA
2 3

−Noeud 3 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

2 3 4

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ − 0, 45ϕ = 0, 55ϕA
2 3

−Noeud 3 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

2 3 4

−Noeud 4 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

3 4 5

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ − 0, 45ϕ = 0, 55ϕA
2 3

−Noeud 3 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

2 3 4

−Noeud 4 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

3 4 5

−Noeud 5 : −0, 55ϕ + 1ϕ = 0, 45ϕB

4 5

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Ainsi, pour le noeud 2, l'équation discrétisée est :aP ϕP − aE ϕE = aW ϕA

Pour le noeud 5, l'équation discrétisée devient :−aW ϕW + aP ϕP = aE ϕB
Cas 1 : F = ρu = 1 × 0, 1 = 0, 1, D = δx Γ
= 00,,12 = 0, 5
aW = Dw + F2w = 0, 55, aE = De + F2e = 0, 45, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ − 0, 45ϕ = 0, 55ϕA
2 3

−Noeud 3 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

2 3 4

−Noeud 4 : −0, 55ϕ + 1ϕ − 0, 45ϕ = 0

3 4 5

−Noeud 5 : −0, 55ϕ + 1ϕ = 0, 45ϕB

4 5

La
 forme matricielle du système à 
résoudre
  est:
1 −0, 45 0 0 0, 55 0, 871
    
ϕ2 ϕ2
 −0, 55 1 − 0, 45 0  ϕ3   0  ϕ3  0, 714
 = ⇒ =
0 −0, 55 1 −0, 45  ϕ4   0  ϕ4  0, 522
 

0 0 −0, 55 1 ϕ5 0 ϕ5 0, 287

Y. EL KHCHINE Mod-Syst-Energ 45 / 57
Application 3 :

Y. EL KHCHINE Mod-Syst-Energ 46 / 57
Application 3 :

Cas 2 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,12 = 0, 5
aW = Dw + Fw
2 = 1, 75, aE = De − F2e = −0, 75, aP = aW + aE + (Fe − Fw ) = 1

Y. EL KHCHINE Mod-Syst-Energ 46 / 57
Application 3 :

Cas 2 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,12 = 0, 5
aW = Dw + Fw
2 = 1, 75, aE = De − F2e = −0, 75, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :

Y. EL KHCHINE Mod-Syst-Energ 46 / 57
Application 3 :

Cas 2 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,12 = 0, 5
aW = Dw + Fw
2 = 1, 75, aE = De − F2e = −0, 75, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ2 + 0, 75ϕ3 = 1, 75ϕA

Y. EL KHCHINE Mod-Syst-Energ 46 / 57
Application 3 :

Cas 2 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,12 = 0, 5
aW = Dw + Fw
2 = 1, 75, aE = De − F2e = −0, 75, aP = aW + aE + (Fe − Fw ) = 1
Les équations à résoudre sont les suivantes :
−Noeud 2 : 1ϕ2 + 0, 75ϕ3 = 1, 75ϕA
−Noeud 3 : −1, 75ϕ2 + 1ϕ3 + 0, 75ϕ4 = 0
Y. EL KHCHINE Mod-Syst-Energ 46 / 57
Application 3 :
−Noeud 4 : −1, 75ϕ3 + 1ϕ4 + 0, 75ϕ5 = 0

Y. EL KHCHINE Mod-Syst-Energ 47 / 57
Application 3 :
−Noeud 4 : −1, 75ϕ3 + 1ϕ4 + 0, 75ϕ5 = 0
−Noeud 5 : −1, 75ϕ4 + 1ϕ5 = −0, 75ϕB

Y. EL KHCHINE Mod-Syst-Energ 47 / 57
Application 3 :
−Noeud 4 : −1, 75ϕ3 + 1ϕ4 + 0, 75ϕ5 = 0
−Noeud 5 : −1, 75ϕ4 + 1ϕ5 = −0, 75ϕB
La
 forme matricielle du système 
à résoudre
  est  :
1 0, 75 0 0 1, 75 0, 952
    
ϕ2 ϕ2
 −1, 75 1 0, 75 0   ϕ3  =  0  ⇒ ϕ3  = 1, 063
       
0 −1, 75 1 0, 75 0 0, 805

   ϕ4    ϕ4  
0 0 −1, 75 1 ϕ5 0 ϕ5 1, 408

Y. EL KHCHINE Mod-Syst-Energ 47 / 57
Application 3 :
−Noeud 4 : −1, 75ϕ3 + 1ϕ4 + 0, 75ϕ5 = 0
−Noeud 5 : −1, 75ϕ4 + 1ϕ5 = −0, 75ϕB
La
 forme matricielle du système 
à résoudre
  est  :
1 0, 75 0 0 1, 75 0, 952
    
ϕ2 ϕ2
 −1, 75 1 0, 75 0   ϕ3  =  0  ⇒ ϕ3  = 1, 063
       
0 −1, 75 1 0, 75 0 0, 805

   ϕ4    ϕ4  
0 0 −1, 75 1 ϕ5 0 ϕ5 1, 408

Y. EL KHCHINE Mod-Syst-Energ 47 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Les équations à résoudre sont les suivantes :

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Les équations à résoudre sont les suivantes :

−Noeud 2 : 4ϕ − 0, 75ϕ = 3, 25ϕA
2 3

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Les équations à résoudre sont les suivantes :

−Noeud 2 : 4ϕ − 0, 75ϕ = 3, 25ϕA
2 3

−Noeud 3 − 19 : −3, 25ϕi− + 4ϕi − 0, 75ϕi+ = 0

1 1

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Les équations à résoudre sont les suivantes :

−Noeud 2 : 4ϕ − 0, 75ϕ = 3, 25ϕA
2 3

−Noeud 3 − 19 : −3, 25ϕi− + 4ϕi − 0, 75ϕi+ = 0

1 1

−Noeud 20 : −3, 25ϕ + 4ϕ = 0, 75ϕB

19 20

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :
Cas 3 : u = 2, 5m/s, F = ρu = 1 × 2, 5 = 2, 5, D = δx
Γ
= 00,,05
1
=2
aW = Dw + F2w = 3, 25, aE = De − 2 = 0, 75, aP = aW + aE + (Fe − Fw ) = 4
Fe

Les équations à résoudre sont les suivantes :

−Noeud 2 : 4ϕ − 0, 75ϕ = 3, 25ϕA
2 3

−Noeud 3 − 19 : −3, 25ϕi− + 4ϕi − 0, 75ϕi+ = 0

1 1

−Noeud 20 : −3, 25ϕ + 4ϕ = 0, 75ϕB

19 20

Y. EL KHCHINE Mod-Syst-Energ 48 / 57
Application 3 :

Conclusion :
On constate une bonne concordance entre la solution numérique et la solution
analytique cette [Link] comparaison des résultats des cas 2 et 3 montre qu'un
maillage plus n (cas 3) a déterminé la réduction du nombre de Peclet de 5 à
1,5. Le schéma aux diérences centrales fournis de bons résultats quad le
nombre de Peclet est petit.

Y. EL KHCHINE Mod-Syst-Energ 49 / 57
Volumes nis pour les problèmes de convection-diusion
2. Schéma décentré (Upwind Dierencing Scheme)
Le schéma dit Upwind prend en considération la direction du ux lors de la
détermination de la propriété Φ sur la face du VC.

Y. EL KHCHINE Mod-Syst-Energ 50 / 57
Volumes nis pour les problèmes de convection-diusion
2. Schéma décentré (Upwind Dierencing Scheme)
Le schéma dit Upwind prend en considération la direction du ux lors de la
détermination de la propriété Φ sur la face du VC.
Lorsque l'écoulement est dans le sens positive, on a :
uw > 0, ue > 0 (Fw > 0, Fe > 0)

Y. EL KHCHINE Mod-Syst-Energ 50 / 57
Volumes nis pour les problèmes de convection-diusion
2. Schéma décentré (Upwind Dierencing Scheme)
Le schéma dit Upwind prend en considération la direction du ux lors de la
détermination de la propriété Φ sur la face du VC.
Lorsque l'écoulement est dans le sens positive, on a :
uw > 0, ue > 0 (Fw > 0, Fe > 0)

Y. EL KHCHINE Mod-Syst-Energ 50 / 57
Volumes nis pour les problèmes de convection-diusion
2. Schéma décentré (Upwind Dierencing Scheme)
Le schéma dit Upwind prend en considération la direction du ux lors de la
détermination de la propriété Φ sur la face du VC.
Lorsque l'écoulement est dans le sens positive, on a :
uw > 0, ue > 0 (Fw > 0, Fe > 0)

Alors l'équation précédente :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

Y. EL KHCHINE Mod-Syst-Energ 50 / 57
Volumes nis pour les problèmes de convection-diusion
2. Schéma décentré (Upwind Dierencing Scheme)
Le schéma dit Upwind prend en considération la direction du ux lors de la
détermination de la propriété Φ sur la face du VC.
Lorsque l'écoulement est dans le sens positive, on a :
uw > 0, ue > 0 (Fw > 0, Fe > 0)

Alors l'équation précédente :

Fe Φe − Fw Φw = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

Y. EL KHCHINE Mod-Syst-Energ 50 / 57
Volumes nis pour les problèmes de convection-diusion
devient :
Fe ΦP − Fw ΦW = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

Y. EL KHCHINE Mod-Syst-Energ 51 / 57
Volumes nis pour les problèmes de convection-diusion
devient :
Fe ΦP − Fw ΦW = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

(Fe + De + Dw )ΦP = (Fw + Dw )ΦW + De ΦE

Y. EL KHCHINE Mod-Syst-Energ 51 / 57
Volumes nis pour les problèmes de convection-diusion
devient :
Fe ΦP − Fw ΦW = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

(Fe + De + Dw )ΦP = (Fw + Dw )ΦW + De ΦE

aP ΦP = aW ΦW + aE ΦE
Résumé du schéma Upwind :
aW aE
Fw > 0, Fe > 0 D w + Fw De
Fw < 0, Fe < 0 Dw De − F e

Y. EL KHCHINE Mod-Syst-Energ 51 / 57
Volumes nis pour les problèmes de convection-diusion
devient :
Fe ΦP − Fw ΦW = De (ΦE − ΦP ) − Dw (ΦP − ΦW )

(Fe + De + Dw )ΦP = (Fw + Dw )ΦW + De ΦE

aP ΦP = aW ΦW + aE ΦE
Résumé du schéma Upwind :
aW aE
Fw > 0, Fe > 0 D w + Fw De
Fw < 0, Fe < 0 Dw De − F e

Ou bien sous la forme compacte suivante en tenant compte les deux directions de
l'écoulement,

aW aE aP
Dw + max(Fw , 0) De + max(0, −Fe ) Dw + max(−Fw , 0) + De + max(Fe , 0) = aW + aE + Fe − Fw

Y. EL KHCHINE Mod-Syst-Energ 51 / 57
Volumes nis pour les problèmes de convection-diusion

Conclusions :
Le schéma UDS vérie la propriété de conservation.
Les coecients sont toujours positifs, d'où la propriété d'être bornée.
La propriété d'être transportée est aussi assurée par ce schéma.

Y. EL KHCHINE Mod-Syst-Energ 52 / 57
Application 4

Résoudre le problème considéré à l'exemple 1 en utilisant le schéma de discrétisatio

"upwind" pour les cas :

1 u = 0.1 m/s ;
2 u = 2.5 m/s ,

et pour le même maillage (à six n÷uds). Comparer les résultats avec la solution
analytique.

Y. EL KHCHINE Mod-Syst-Energ 53 / 57
Application 4

Y. EL KHCHINE Mod-Syst-Energ 54 / 57
Application 4

Y. EL KHCHINE Mod-Syst-Energ 55 / 57
Application 4

Y. EL KHCHINE Mod-Syst-Energ 56 / 57
Application 4

Y. EL KHCHINE Mod-Syst-Energ 57 / 57