E.S.I.

M
Ecole Supérieure des Ingénieurs de Medjez El Beb

Compte rendu :

Transfer thermique

Elaboré par : Wissem Ben Mohamed

Enseignant : Mr. Chtourou Wajdi
Classe : 1ére année HM

A.U :2023/2024
Sommaire :

 Exercice …………………………………………………………………………… :

 Objective ………………………………………………………………………… :

 Hypothèse ……………………………………………………………………… :

 Discrétisation de l’équation de conduction de chaleur en 2D… :

 Conditions aux limites…………………………………………………………… :

 Résolution

numérique…………………………………………………………… :

 Code python …………………………………………………………………… :

 Les étapes de python ……………………………………………………… :

 Code Matlab………………………………………………………………………… :

 Les étapes de Matlab …………………………………………………… :

 Premier Principe de la Thermodynamique …………………… :

 Deuxième Principe de la Thermodynamique………………… :

 Conclusion …………………………………………………………………… :
 Travail demandé

1) Exercice :
Vous avez une cheminée industrielle de forme carrée
creuse. Les températures aux faces de la cheminée sont
spécifiées comme suit :

- Les faces extérieures parallèles à (a) et (a’) sont à une

température de 200°C.
- Les faces intérieures (n) et (n’) sont à une
température de 360°C.
- Les faces extérieures parallèles à (b) et (b’) sont à une
température de 235°C.
- Les faces intérieures (m) et (m’) sont à une
température de 380°C.
2)Objective :
Déterminer la distribution de température (T
(x, y)) dans la cheminée en régime
stationnaire. Le matériau de la cheminée est
homogène avec une conductivité thermique
uniforme (lambda).

3)Hypothèse :
-Le régime est stationnaire (pas de variation
temporelle de la température).
- La différence de température entre les
points adjacents de la grille (ΔT) est nulle.
- Le problème est formulé dans un repère
cartésien, et nous utiliserons une approche
discrétisée.
-pas de production de chaleur à l’intérieur
=0

4)Discrétisation de l'équation de conduction de

chaleur en 2D :
Nous allons utiliser la méthode des différences finies
pour discrétiser l'équation de conduction de chaleur en
2D. Supposons que nous avons une grille N×M où N est
le nombre de points dans la direction x et M est le
nombre de points dans la direction y.
La discrétisation de l'équation de conduction de chaleur
donne :

5)Conditions aux limites :

  Sur les faces extérieures parallèles à a et a’ :
T(x,0) = T (x, M−1) =200°C.
 Sur les faces intérieures n et n’ : T (0, y) =T
(N−1, y) =360°C
 Sur les faces extérieures parallèles à bet b’:
T (0, y) =T (N−1, y) =235°C
 Sur les faces intérieures m et m′ : T(x,0) =T (x, M−1)
=380°C

6)Résolution numérique :
Nous allons résoudre le système d'équations
discrètes de manière itérative jusqu'à ce qu'une
solution convergente soit obtenue. Pour cela,
nous pouvons utiliser des méthodes itératives
telles que la méthode de Gauss-Seidel ou la
méthode de relaxation.

-Partie calcule :
Nous définissons les températures sur les bords
de la cheminée en utilisant les conditions aux
limites spécifiées :

o Sur les bords parallèles à a et a′ :

T(i,0) =T (i, M−1) =200°C pour i=1, 2, ..., N−2
 o Sur les bords parallèles à b et b′ :
T (0, j) =T (N−1, j) =235°C pour j=1, 2 ,..., M−2

o Sur les bords intérieurs n et n′

T(0,j)=T(N−1,j)=360°Cpour j=1,2,...,M−2

o Sur les bords intérieurs m et m′ :

T(i,0) =T (i, M−1) =380°Cpour i=1, 2,..., N−2

Itérations :

À chaque itération, nous mettons à jour la température

T(xi,yi) à l'intérieur de la cheminée en utilisant la
moyenne des températures de ses voisins, selon
l'équation suivante :

: est la température à l'itération k au point (xi

,yj),

: est la température mise à jour à l'itération

k+1 au point (xi, yi).
7)Code python :
8)Les étapes de python :

1)Initialisation des températures :

 La fonction commence par initialiser les
températures sur les bords de la cheminée en
fonction des conditions aux limites spécifiées dans
l'énoncé du problème. Ces conditions aux limites
sont définies en fonction des températures sur les
bords parallèles à a, ′a′, b, ′b′, n, et m.

2)Boucle d'itération pour la résolution :

 Ensuite, la fonction entre dans une boucle

d'itération pour résoudre le problème de
distribution de température dans la cheminée.
 À chaque itération de la boucle, les
températures à l'intérieur de la cheminée sont
mises à jour en utilisant la méthode de
relaxation de Jacobi. Pour chaque point de la
grille à l'exception des bords, la nouvelle
température est calculée comme la moyenne
pondérée des températures de ses voisins.
 Cette mise à jour itérative se poursuit jusqu'à
ce que la différence entre les températures de
deux itérations successives soit inférieure à
une tolérance prédéfinie ou que le nombre
maximal d'itérations soit atteint.

3)Retour de la grille de température résolue :

 Une fois que la solution convergée est obtenue, la

fonction retourne la grille de température résolue,
qui représente la distribution de température dans
la cheminée en régime stationnaire.
 9) Code Matlab :
10)Les étapes de Matlab :

1. Initialisation des températures sur les bords de la

cheminée :
 La fonction `solve_heat_equation` commence par
initialiser les températures sur les bords de la
cheminée en fonction des conditions aux limites
spécifiées dans l'énoncé du problème. Ces
conditions aux limites sont définies en fonction
des températures sur les bords parallèles à a, a', b,
b', n, et m.

2. Boucle d'itération pour la résolution :

 Ensuite, la fonction entre dans une boucle
d'itération pour résoudre le problème de
distribution de température dans la cheminée.
 À chaque itération de la boucle, les températures à
l'intérieur de la cheminée sont mises à jour en
utilisant la méthode de relaxation de Jacobi. Pour
chaque point de la grille à l'exception des bords, la
nouvelle température est calculée comme la
moyenne pondérée des températures de ses
voisins.
 Cette mise à jour itérative se poursuit jusqu'à ce
que la différence entre les températures de deux
itérations successives soit inférieure à une
tolérance prédéfinie ou que le nombre maximal
d'itérations soit atteint.
3. Retour de la grille de température résolue :
 Une fois que la solution convergée est obtenue, la
fonction retourne la grille de température résolue,
qui représente la distribution de température dans
la cheminée en régime stationnaire.

4. Affichage de la distribution de température :

 Après avoir résolu le problème, vous pouvez
afficher la distribution de température en utilisant
`surf` pour créer un tracé surfacique en 3D. Cela
vous permet de visualiser la variation de
température dans la cheminée en fonction de sa
position.

11)Premier Principe de la Thermodynamique :

Le premier principe de la thermodynamique, également connu sous
le nom de principe de conservation de l'énergie, énonce que l'énergie
totale d'un système isolé est constante. Autrement dit, l'énergie peut
être transformée d'une forme à une autre (par exemple, de l'énergie
thermique à l'énergie mécanique), mais elle ne peut ni être créée ni
détruite. Matériellement, cela s'exprime par la formule suivante :
ΔU=Q−W
où :
 ΔU est la variation de l'énergie interne du système,
 Q est la chaleur ajoutée au système,
 W est le travail effectué par le système.
12)Deuxième Principe de la Thermodynamique :
Le deuxième principe de la thermodynamique stipule que la chaleur
ne peut pas spontanément passer d'un corps froid à un corps chaud
sans travail externe. Ce principe introduit la notion d'entropie, une
mesure de l'ordre ou du désordre d'un système. L'entropie totale
d'un système isolé ne peut jamais diminuer au fil du temps. Ce
principe est souvent résumé par les énoncés suivants :

Énoncé de Clausius : Il est impossible de construire une machine

fonctionnant en cycle qui transfère de la chaleur d'un réservoir froid
à un réservoir chaud sans l'apport de travail externe.

Énoncé de Kelvin-Planck : Il est impossible de construire une

machine fonctionnant en cycle qui convertit entièrement la chaleur
en travail sans aucune perte d'énergie sous forme de chaleur
résiduelle.

Machine de Carnot :
La machine de Carnot est une machine thermique idéale qui
fonctionne selon un cycle réversible appelé cycle de Carnot. Elle est
composée de deux réservoirs de température (un chaud et un froid)
et réalise quatre transformations thermodynamiques réversibles :

Expansion isotherme : La machine absorbe de la chaleur Q1

d'un réservoir chaud à température constante T1 et le gaz se dilate.
 Expansion adiabatique : Le gaz continue de se dilater sans
échange de chaleur avec l'extérieur, ce qui diminue sa température
de T1 à T2.

Compression isotherme : Le gaz est comprimé à température

constante T2et rejette de la chaleur Q2 dans un réservoir froid.

Compression adiabatique : Le gaz continue d'être comprimé

sans échange de chaleur, augmentant sa température de T2 à T1.

L'efficacité (ou rendement) d'une machine de Carnot est donnée par

la formule :
Conclusion :

Dans cet exercice, nous avons déterminé la

distribution de température dans une cheminée
industrielle carrée creuse en utilisant une approche
numérique. Voici les principaux points :

1. Méthode Utilisée :
- Nous avons discrétisé l'équation de conduction de
la chaleur en 2D avec la méthode des différences
finies.
- Les conditions aux limites spécifiques ont été
appliquées aux différentes faces de la cheminée.

2. Résolution:
- La méthode de Gauss-Seidel a été utilisée pour
obtenir une solution itérative, convergeant vers la
distribution de température en régime stationnaire.

3. Visualisation :
- La distribution de température a été visualisée,
montrant comment la chaleur se répartit à l'intérieur
de la cheminée.
Principes Thermodynamiques

En complément, nous avons rappelé deux principes

essentiels de la thermodynamique et la machine de
Carnot :

- Premier Principe : Conservation de l'énergie dans un

système.
- Deuxième Principe : Augmentation de l'entropie et
impossibilité de transférer spontanément de la
chaleur d'un corps froid à un corps chaud.
- Machine de Carnot : Machine thermique idéale avec
un rendement maximum théorique.

Résultat Pratique

Cet exercice démontre l'utilisation de techniques

numériques pour résoudre des problèmes de
conduction thermique, cruciales pour la conception et
l'optimisation de systèmes industriels.