Cycle d’ingénieur

Filière Génie mécatronique d’automobile

Transfert thermique et mecanique des

fluides
--------------------------
Pr. A. SAAD

Année universitaire: 2021/2022

CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent.
- Géométrie prismatique (le mur plan).
a-Le mur où la conductivité est constante.
b-Le mur où la conductivité varie avec la température.
c- Le mur avec dissipation d ’énergie.
- Géométrie cylindrique.
a-Cylindre plein.
b-Cylindre creux à surface latérale isotherme.
c- Cylindre avec dissipation d ’énergie.
2. Conduction bidimensionnelle en régime permanent.
- Concept de facteur de forme entre deux surfaces isothermes
(surfaces équivalentes).
a-Facteur de forme pour une paroi prismatique (Mur 3D).
b-Facteur de forme pour une paroi cylindrique.
- Méthode numérique (Différences finies).
3. Conduction en régime transitoire.
- Conduction avec températures limites imposées.
- Méthode numérique (Différences finies) en régime transitoire.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent

Ou encore
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
- Géométrie prismatique (le mur plan).
a- Le mur où la conductivité est constante.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
- Géométrie prismatique (le mur plan).
a- Le mur où la conductivité est constante.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
b- Le mur où la conductivité varie avec la température.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
b- Le mur où la conductivité varie avec la température.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
b- Le mur où la conductivité varie avec la température.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c-
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c-
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c-
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
- Géométrie cylindrique.
a-Cylindre plein.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
- Géométrie cylindrique.
a. Cylindre plein.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c-Cylindre plein avec dissipation d’ énergie.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c- Cylindre plein avec dissipation d’ énérgie.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
c-Cylindre plein avec dissipation d’ énérgie.
CHAPITRE 3:
Résolution de l’équation de la chaleur
1. Conduction unidimensionnelle en régine permanent
CHAPITRE 3:
Résolution de l’équation de la chaleur
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.

- Concept de facteur de forme entre deux surfaces isothermes (surfaces

équivalentes)
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
a-Facteur de forme pour une paroi prismatique (Mur 3D).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
a-Facteur de forme pour une paroi prismatique (Mur 3D).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
a- Facteur de forme pour une paroi prismatique (Mur 3D).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
a- Facteur de forme pour une paroi prismatique (Mur 3D).

b- Facteur de forme pour une paroi cylindrique.

CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
b-Facteur de forme pour une paroi cylindrique.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).

Principe : Cas 1D

On partage l’intervalles [0,1] en N intervalles identiques.

xm-1 xm xm+1
0 1

Soit xm=mh telle que h=xm+1-xm=1/N et 0<m<N-1, h représente le pas d’espace.

CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Formule de Taylor (Rappel):
Une fonction plusieurs fois dérivable au voisinage d'un point peut être approximée
par une fonction polynôme dont les coefficients dépendent uniquement des
dérivées de la fonction en ce point. Cette fonction polynôme est parfois
appelée polynôme de Taylor.
Taylor
I un intervalle réel, a un élément de I, n un entier naturel, f une fonction de I dérivable
en a jusqu'à l'ordre n.
Alors pour tout nombre réel x appartenant à I, on a :

où le reste Rn(x) est une fonction négligeable par rapport à (x – a)n au voisinage
de a.
Par simple changement de variable, la formule de Taylor peut aussi s'exprimer
sous la forme :
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Cas 1D
a) Formule de Taylor :
u∈C4(0,1)
i- u(xm-1)= u(xm)-h u’(xm)+h2/2!* u(2)(xm)- h3/3!* u(3) (xm)+h4/4!* u4(xm)………
ii- u(xm+1)= u(xm)+h u’(xm)+h2/2!* u(2) (xm)+ h3/3!* u(3) (xm)+ h4/4!* u4 (xm )………
iii- u(xm+2)= u(xm)+2h u’(xm)+(2h)2/2* u(2) (xm)+(2h)3/3!* u(3) (xm)+(2h)4/4!* u4(xm)…
Pour les équations du premier ordre :
u (x m +1 ) − u (x m )
u '(x m ) = lim
h →0 h

u (x m +1 ) − u (x m )
A la limite la dérivée devient u '(x m ) =
: h
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
a) Technique pour trouver une approximation des dérivées u(k)(Xm) :

Faire disparaître toutes les dérivées de u(l)(xm) pour l<K, en faisant les
combinaisons de i, ii et iii.

Garder la donnée u(K)(xm).

Faire disparaître un certain nombre de dérivées u(l)(xm ) avec l >K,
pour obtenir une certaine précision.

b) Applications:
∂u (x m ) (u m +1 − u m ) Schéma décentré à droite (décentré aval)
= + θ (h )
∂x h

∂u (x m ) (u m − u m −1 ) Schéma décentré à gauche (décentré amont)

= + θ (h )
∂x h
∂u (x m ) (u m +1 − u m −1 ) Schéma centré
= + θ (h 2 )
∂x 2.h

(u − 2u m + u m −1 ) + θ (h 2 ) h2 (4)
(i)+(ii) donne : u (2) (x ) = m +1 L’erreur est de l’ordre de h2 = u (xm )
m
h 2 12
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
 Les conditions aux limites imposant sur un bord une température de surface
s’expriment simplement en fixant la valeur de la température T(m,n) à la valeur
imposée pour tout couple (m, n) représentant un point de ce bord.
 Les conditions aux limites avec transfert convectif ou flux imposé s’expriment de
la manière suivante :

• Convection dans un coin intérieur

Bilan Nœud T(m,n)

CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
• Convection sur la surface

Bilan Nœud T(m,n)

CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
• Convection sur un coin extérieur

Bilan Nœud T(m,n)

Dans le cas d'une surface isolée, il suffit de prendre les relations avec
convection et de poser h=0.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
• Convection sur un coin extérieur
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Conditions aux limites:
• Convection sur un coin extérieur
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
Résolution: Le système d’équations linéaires étant obtenu, on peut le résoudre par:
• Une méthode directe (inversion, pivot, Gauss, Gauss - Jordan).
• Une méthode itérative (Gauss - Seidel).
- Cette dernière est facile à programmer.

- On effectue des itérations successives consistant à remplacer la valeur de la température

en chaque nœud du maillage par la valeur calculée par l’équation aux différences finies qui
lui est associée.

- Une itération consiste à effectuer un balayage complet de tous les nœuds, ligne après ligne
et de gauche à droite pour chaque ligne par exemple.

- Les valeurs recalculées sont immédiatement prises en compte pour le calcul de la valeur de
la température T aux points d’ordre supérieur .

Critère de convergence: On peut par exemple arrêter le calcul dès que la variation la plus
grande de T (m,n) au cours d’une itération reste inférieure à une valeur ε donnée.
CHAPITRE 3:
Résolution de l’équation de la chaleur

2. Conduction bidimensionnelle en régime permanent.

- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies).
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.
- Méthode numérique (Différences finies) en régime transitoire.
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.
- Méthode numérique (Différences finies) en régime transitoire.
CHAPITRE 3:
Résolution de l’équation de la chaleur

3. Conduction en régime transitoire.

- Méthode numérique (Différences finies) en régime transitoire.

Merci pour
votre attention
CHAPITRE 3:
Transfert de chaleur par conduction

En physique, un milieu semi-infini est un milieu hypothétique

occupant tout un demi-espace, c'est-à-dire s'étendant à
l'infini dans les trois directions mais d'un seul côté d'un plan.
Dans un repère orthonormé, il est par exemple défini par x
quelconque, y quelconque et z>0.
CHAPITRE 3:
Transfert de chaleur par conduction
Résolution de l’équation de conduction
En physique, un milieu semi-
infini est un milieu hypothétique
occupant tout un demi-espace,
c'est-à-dire s'étendant à l'infini
dans les trois directions mais
d'un seul côté d'un plan. Dans
un repère orthonormé, il est par
exemple défini par x
quelconque, y quelconque et
z>0.
CHAPITRE 3:
Transfert de chaleur par conduction
Résolution de l’équation de conduction
CHAPITRE 3:
Transfert de chaleur par conduction
Résolution de l’équation de conduction
CHAPITRE 3:
Transfert de chaleur par conduction
Résolution de l’équation de conduction
CHAPITRE 3:
Transfert de chaleur par conduction
Résolution de l’équation de conduction
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies). Cas 2D

 −∆ u = f x∈Ω ⊂ 2
y
(2 D ) 
 u Γ = 0 Γ = ∂Ω

 ∂ 2u ∂ 2u 2
− 2 − 2 = f (x , y ) x ∈Ω ⊂
(2D )  ∂x ∂y
u (x , y ) = 0 Γ =∂Ω
 Γ x

Nous avons pour le cas où :

y
 −u m +1,n − u m −1,n + 4u m ,n − u m ,n +1 − u m ,n −1
−∆ u m ,n = 2
= f (x i , y j )
 h
u m ,0 = u1,n = u m ,1 = u0,n = 0 conditions aux limites


Avec hx = x m +1 − x m
hy = y n +1 − y n et hx = hy
x
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode numérique (Différences finies). Cas 2D

Le système (2D) peut s’écrire donc sous la forme: Aui = b

 4 -1 0 -1 0 0 0 0 0  u11   f11 
  
 -1 4 -1 0 -1 0 0 0 0  


u 21   f 21 
 0 -1 4 0 0 -1 0 0 0  u 31   f 31 
     
-1 0 0 4 -1 0 -1 0 0 .   . 
 
1
A= 2 0 -1 0 -1 4 -1 0 -1 0 u = .  b= . 
h      
 0 0 -1 0 -1 4 0 0 -1   .   . 
  . 
0 0 0 -1 0 0 4 -1 0  .   
   
 0 0 0 0 -1 0 -1 4 -1   .   . 
  f 33 
 0 0 0 0 0 -1 0 -1 4  
 u 33  

x
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.

- Concept de facteur de forme entre deux surfaces isothermes (surfaces

équivalentes)

a-Facteur de forme pour une paroi prismatique (Mur 3D).

CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
a- Facteur de forme pour une paroi prismatique (Mur 3D).
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.

- Conduction avec températures limites imposées.

CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.
- Conduction avec températures limites imposées.
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.
- Conduction avec températures limites imposées.
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.
- Conduction avec températures limites imposées.
CHAPITRE 3:
Résolution de l’équation de la chaleur
3. Conduction en régime transitoire.

- Conduction avec températures limites imposées.

CHAPITRE 3:
Résolution de l’équation de la chaleur

3. Conduction en régime transitoire.

- Conduction avec températures limites imposées.
CHAPITRE 3:
Résolution de l’équation de la chaleur

3. Conduction en régime transitoire.

- Conduction avec températures limites imposées.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.
CHAPITRE 3:
Résolution de l’équation de la chaleur
2. Conduction bidimensionnelle en régime permanent.
- Méthode de séparation des variables.