Ministère de l'Enseignement Supérieur et de la Recherche Scientifique
ECOLE NATIONALE D’INGENIEURS DE MONASTIR
DEPARTEMENT GENIE ENERGETIQUE
COMPTES RENUS DE
TRAVEAUX PRATIQUES
METHODES NUMERIQUES
Filière :Génie énergétique /Niveau : 2 ère année /TP :2
Réalisés par : Mejdi Tesnim
Enseignante : [Link] Ben Ltaifa
Année Universitaire : 2025-2026
Table des matières
Manipulation 1: ................................................................................................................ 2
I. Introduction ................................................................................................................. 2
II. But de la manipulation: ............................................................................................... 2
III. Problème étudié: ....................................................................................................... 3
1-Calcul de solution analytique :..................................................................................... 3
IV. Algorithme et programme : ............................................................................................ 6
V. Conclusion.................................................................................................................10
Manipulation 2: ...............................................................................................................11
I. Introduction ................................................................................................................11
II. But de la manipulation: ..............................................................................................11
III. Problème étudié: ......................................................................................................11
1. Calculer la solution analytique du profil de la température :.......................................12
IV. Algorithme et programme : ...........................................................................................15
VI. Graphe et Interprétation : ..........................................................................................21
VII. Conclusion ..............................................................................................................22
Manipulation 3: ...............................................................................................................24
I. Objectif : ........................................................................................................................24
II. Problème étudié : ..........................................................................................................24
[Link] de la solution analytique : ................................................................................25
III. Algorithme et programme : ...........................................................................................28
IV. Conclusion................................................................................................................32
1
Manipulation 1:
I. Introduction
Dans cette manipulation, nous étudions le transfert de chaleur stationnaire le long
d’une barre cylindrique soumise à une conduction unidimensionnelle. L’objectif est de
déterminer le profil de température en résolvant l’équation de la chaleur à l’aide d’une
solution analytique puis d’une méthode numérique par différences finies. Cette étude
permet de comparer les deux approches et de comprendre l’influence des conditions
aux limites et du maillage sur la précision des résultats.
II. But de la manipulation:
Cette manipulation vise à analyser le transfert de chaleur stationnaire en régime
unidimensionnel dans une barre cylindrique dépourvue de source interne. L’objectif est
de déterminer la distribution de température le long de la barre en résolvant l’équation
de la diffusion thermique, d’une part par une approche analytique et d’autre part par
une méthode numérique, puis de confronter les résultats obtenus.
2
III. Problème étudié:
On considère une barre cylindrique de longueur L = 0,5 m, dont les extrémités A et B
sont maintenues à des températures constantes de 100°C et 500°C respectivement. Le
transfert thermique est supposé stationnaire et monodimensionnel. La conductivité
thermique du matériau est λ = 1000 W/m.K. Pour la résolution numérique, on adopte un
maillage 1D à six nœuds, uniformément répartis le long de la barre, avec un pas Δx = 0,1
m.
1-Calcul de solution analytique :
On a l'équation de diffusion stationnaire sans source est :
Après tout intégration faite :
Avec les conditions aux limites :
Équations nodales
Pour chaque nœud intérieur, l’équation discrète du bilan thermique s’écrit :
• Nœud 2 :
• Nœud 3 :
3
• Nœud 4
• Nœud 5 :
Système d’équations linéaires
Le système peut être écrit sous forme matricielle :
le profil de Température complet :
Noeud Position (m) Température (°C) Température (K)
1 0 100 273
2 0.1 180 453
3 0.2 260 533
4 0.3 340 613
5 0.4 420 693
6 0.5 500 773
4
Équation discrétisée :
L’équation de conduction 1D stationnaire :
𝑑2𝑇
=0
𝑑𝑥 2
se discrétise en différences finies centrales :
𝑇𝑖−1 + 2𝑇𝑖 + 𝑇𝑖+1
=0
(Δ𝑥)2
Donc :
𝑇𝑖−1 − 2𝑇𝑖 + 𝑇𝑖+1 = 0
On réarrange pour obtenir la forme souhaitée :
𝑎𝑖 𝑇𝑖 = 𝑎𝑖1 𝑇𝑖−1 + 𝑎𝑖2 𝑇𝑖+1 + 𝑏𝑖
Ici :
• 𝑎𝑖 = 2
• 𝑎𝑖1 = 1
• 𝑎𝑖2 = 1
• 𝑏𝑖 = 0
Donc pour les nœuds intérieurs :
2𝑇𝑖 = 𝑇𝑖−1 + 𝑇𝑖+1
ou sous forme standard :
5
𝑇𝑖−1 + 2𝑇𝑖 + 𝑇𝑖+1 = 0
On obtient :
Noeud Pi Qi
2 0.5 50
3 0.66 33
4 0.75 24.75
5 0.8 10.625
IV. Algorithme et programme :
Le code utilisé :
6
7
Suite au lancement du programme Fortran, le fichier de résultats a été ouvert à l'aide de
l'éditeur de texte Bloc-notes. Ce fichier a permis d'obtenir et de comparer directement
les valeurs de températures calculées par la solution numérique et celles fournies par la
solution analytique. L'affichage côte à côte de ces deux jeux de données, qui présentent
un accord très satisfaisant, permet de valider immédiatement l'exactitude et l'efficacité
de notre implémentation de la méthode des volumes finis. La parfaite concordance des
résultats numériques avec la référence analytique confirme la robustesse de la
simulation et la correcte résolution du système d'équations.
8
On obtient Le résultat analytique suivant :
On obtient Le résultat numériques suivant :
Afin de comparer visuellement les performances des solutions analytique et
numérique, les valeurs de température obtenues ont été entrées dans un tableur Excel.
L'abscisse a été définie par la position normalisée des nœuds (delta X), représentant la
discrétisation du domaine. Les ordonnées correspondent aux températures analytique
et numérique.
9
Interprétation :
Les résultats numériques et théoriques évoluent de la même manière linéairement, ce
qui est logique. Cependant, la courbe de la simulation est toujours en dessous de la
courbe de référence : pour chaque mesure, sa valeur est trop faible. Cet écart n'est pas
dû au hasard. Il prouve qu'il y a une erreur systématique (un biais) dans le code de
calcul, qui fait qu'il donne systématiquement une température inférieure à la réalité, et
ce, partout.
V. Conclusion
La solution numérique reproduit correctement la tendance linéaire attendue, mais elle présente
systématiquement des températures inférieures à celles de la solution théorique. Cet écart
constant suggère une erreur dans un paramètre ou une formule du programme Fortran, par
exemple au niveau des températures aux limites ou de la conductivité thermique. Il est donc
nécessaire de vérifier ces valeurs afin de corriger le modèle et d’obtenir des résultats plus
précis.
10
Manipulation 2:
I. Introduction
Cette manipulation porte sur l’étude du transfert thermique stationnaire unidirectionnel
par conduction dans une plaque soumise à un terme source volumique constant. Le
système considéré présente une géométrie simple : une plaque d’épaisseur 𝐿, dont les
deux faces sont maintenues à des températures imposées.
L’objectif du travail est de résoudre l’équation de la chaleur associée, à la fois
analytiquement et numériquement, à l’aide d’un maillage unidimensionnel régulier.
Cette étude permet de comparer les résultats obtenus par les deux approches et
d’évaluer l’influence du terme source sur le profil de température à travers la plaque.
II. But de la manipulation:
Le but de cette manipulation est d’étudier le transfert thermique stationnaire
unidirectionnel dans une plaque soumise à un terme source volumique constant. Elle
vise à :
• résoudre analytiquement le profil de température à partir de l’équation de la
chaleur,
• établir une solution numérique à l’aide d’un maillage unidimensionnel régulier,
• comparer les deux approches afin d’évaluer l’effet du terme source sur la
distribution de température.
III. Problème étudié:
On considère une plaque d’épaisseur L = 20 mm, de conductivité thermique constante
λ = 0,5 W/m·K, soumise à une source volumique uniforme S = 1000 kW/m³.
Les deux faces de la plaque sont maintenues à des températures imposées de 100 °C et
200 °C.
Le transfert thermique est supposé stationnaire et strictement monodimensionnel
selon l’axe x.
La plaque est discrétisée en six nœuds selon un maillage 1D uniforme de pas
Δx = 0,004 m.
11
1. Calculer la solution analytique du profil de la température :
L'équation de la chaleur stationnaire avec terme source est :
Équations nodales
Pour chaque nœud intérieur, l’équation discrète du bilan thermique s’écrit :
• Nœud 2
2𝑇2 = 𝑇3 + 𝑇1 + 𝑏
En remplaçant 𝑇1 = 100∘ 𝐶et 𝑏 = 32:
2𝑇2 = 𝑇3 + 100 + 32
2𝑇2 = 𝑇3 + 132
• Nœud 3
2𝑇3 = 𝑇4 + 𝑇2 + 32
• Nœud 4
2𝑇4 = 𝑇5 + 𝑇3 + 32
• Nœud 5
2𝑇5 = 200 + 𝑇4 + 32
2𝑇5 = 𝑇4 + 232
12
Système d’équations linéaires
Le système peut être écrit sous forme matricielle :
2 −1 0 0 𝑇2 132
−1 2 −1 0 𝑇3 32
[ ][ ] = [ ]
0 −1 2 −1 𝑇4 32
0 0 −1 2 𝑇5 232
Résolution du système
Ligne 1 :
2𝑇2 − 𝑇3 = 132
𝑇2 = 0.5 𝑇3 + 66
Ligne 2 :
−𝑇2 + 2𝑇3 − 𝑇4 = 32
Substitution :
−(0.5𝑇3 + 66) + 2𝑇3 − 𝑇4 = 32
1.5𝑇3 − 𝑇4 = 98
2 196
𝑇3 = 𝑇4 +
3 3
Ligne 3 :
−𝑇3 + 2𝑇4 − 𝑇5 = 32
Substitution :
2 196
− (3 𝑇4+ 3 ) + 2𝑇4 − 𝑇5 = 32
4 292
𝑇4 − 𝑇5 =
3 3
𝑇4 = 0.75 𝑇5 + 73
Ligne 4 :
−𝑇4 + 2𝑇5 = 232
Substitution :
13
−(0.75𝑇5 + 73) + 2𝑇5 = 232
1.25𝑇5 = 305
𝑇5 = 244
le profil de Température complet :
Noeud Position (m) Température (°C) Température (K)
1 0 100 273
2 0.004 184 457
3 0.008 236 509
4 0.012 256 529
5 0.016 244 517
6 0.02 200 473
Équations discrétisées pour les nœuds intérieurs :
En appliquant le schéma aux différences finies au modèle stationnaire 1D, on obtient
les équations nodales suivantes :
Réduction du système tridiagonal
Le système obtenu est tridiagonal. Il peut être transformé en un système à deux diagonales en
utilisant une méthode de type Thomas, menant à une relation de récurrence de la forme :
𝑇𝑖 = 𝑃𝑖 𝑇𝑖+1 + 𝑄𝑖
où les coefficients 𝑃𝑖 et 𝑄𝑖 sont définis par :
𝑎𝑖+1
𝑃𝑖 =
𝑎𝑖 − 𝑎𝑖−1 𝑃𝑖−1
𝑏𝑖 + 𝑎𝑖−1 𝑄𝑖−1
𝑄𝑖 =
𝑎𝑖 − 𝑎𝑖−1 𝑃𝑖−1
14
On obtient :
Noeud Pi Qi
2 0.5 66
3 0.66 65.33
4 0.74 72.63
5 0.79 83.04
IV. Algorithme et programme :
Pour résoudre numériquement ce problème de transfert thermique, nous avons
développé un programme en Fortran utilisant la méthode des différences finies.
L’algorithme mis en place repose sur une démarche progressive comprenant les étapes
suivantes :
• Initialisation des paramètres physiques (conductivité, terme source,
températures imposées) et définition du maillage unidimensionnel.
• Construction du système tridiagonal issu de la discrétisation de l’équation de
la chaleur.
• Application des conditions aux limites, qui modifient les premiers et derniers
termes du second membre.
• Résolution du système linéaire à l’aide d’une méthode de type Thomas, c’est-à-
dire une version optimisée de l’élimination de Gauss adaptée spécifiquement
aux matrices tridiagonales.
15
Le code :
16
17
On obtient Le résultat analytique suivant :
Et le Résultat Numérique suivant :
V. Graphe et interprétation :
À l’aide d’Excel, nous traçons les courbes de température obtenues par les solutions
analytique et numérique.
18
Interprétation :
Convergence des solutions :
On constate que la solution numérique se rapproche progressivement de la solution
analytique lorsque la valeur de Δx diminue, ce qui confirme la convergence de la
méthode des différences finies.
Précision numérique :
Pour des pas Δx très faibles, les deux courbes deviennent quasiment confondues. Cela
montre la grande précision de la méthode numérique. Les écarts résiduels observés
proviennent des erreurs de troncature liées à la discrétisation.
Cas de 𝑺 = −𝟏𝟎𝟎𝟎 kW/m³ :
Afin de compléter l’étude, nous avons appliqué la même démarche numérique en
considérant un terme source négatif 𝑆 = −1000 kW/m³. Ce cas correspond
physiquement à un puits de chaleur uniforme dans le matériau.
La résolution du système tridiagonal modifié met en évidence un profil de température
totalement différent, marqué par l’apparition d’un minimum thermique significatif au
centre du domaine.
Un léger ajustement a également été apporté à la solution analytique (T_anal) pour tenir
compte du changement de signe du terme source.
19
On obtient comme solution Analytique :
20
Et pour la solution Numérique :
VI. Graphe et Interprétation :
Interprétation :
L’introduction d’un terme source négatif modifie profondément le comportement
thermique du système. Au lieu du profil linéaire attendu, la température adopte une
courbure présentant un minimum marqué au centre du domaine.
Ce phénomène met en évidence l’équilibre entre l’absorption d’énergie par le puits de
chaleur et les conditions aux limites imposées. Cette configuration illustre de manière
claire les mécanismes de conduction en présence d’un terme d’absorption thermique .
21
Interprétation :
Lorsque l’on superpose les deux solutions numériques correspondant aux cas S = +1000 W/m³
et S = –1000 W/m³, on observe deux comportements thermiques opposés :
• Pour S = +1000 W/m³ (source de chaleur) :
Le profil de température présente un maximum au centre du domaine. La chaleur
générée à l’intérieur de la plaque provoque une accumulation énergétique, ce qui élève
la température et crée un bombement net de la courbe.
• Pour S = –1000 W/m³ (puits de chaleur) :
La courbe montre au contraire un minimum au centre. Le terme source négatif agit
comme un puits absorbant l’énergie thermique, entraînant une chute de température
vers la zone centrale.
Cette comparaison illustre clairement l’influence du signe du terme source sur l’allure du profil
thermique :
• S > 0 augmente la température interne,
• S < 0 l’abaisse,
tout en respectant les conditions aux limites imposées.
VII. Conclusion
Cette manipulation a permis d’étudier l’effet du terme source sur la distribution de
température au sein d’un milieu solide soumis à un régime permanent. Les résultats
numériques montrent que :
• un terme source positif (S > 0) provoque une élévation interne de la température,
conduisant à un maximum au centre du domaine ;
22
• un terme source négatif (S < 0) engendre au contraire une diminution interne,
caractérisée par un minimum central.
La comparaison des deux cas illustre clairement la symétrie entre génération et
absorption de chaleur et confirme la cohérence de la méthode numérique utilisée.
Cette manipulation démontre ainsi l’importance du terme source dans la modélisation
des transferts thermiques et la fiabilité de la résolution par différences finies.
23
Manipulation 3:
I. Objectif :
L’objectif de cette manipulation est de résoudre numériquement l’équation de transfert
de chaleur stationnaire unidimensionnelle dans une barre cylindrique, avec une
extrémité maintenue à température constante et l’autre isolée (flux nul). L’étude inclut
également l’effet de la convection à la surface de la barre, en comparant les résultats
analytiques et numériques pour différents pas de maillage.
II. Problème étudié :
On considère une barre cylindrique de longueur 𝐿 = 1 met de section 𝐴. L’extrémité 𝑥 =
0est maintenue à 𝑇 = 100 °Ctandis que l’extrémité 𝑥 = 𝐿est isolée (𝑞 = 0). La barre
échange de la chaleur par convection avec l’air ambiant à 𝑇∞ = 20 °C, caractérisé par
un coefficient de convection ℎ.
Le maillage utilisé pour la résolution numérique est unidimensionnel à 6 nœuds, avec
un pas initial Δ𝑥 = 0.2 m.
24
[Link] de la solution analytique :
En conduction 1D avec convection latérale, l’équation stationnaire est :
𝑑 2 𝑇 ℎ𝑃
− (𝑇 − 𝑇∞ ) = 0
𝑑𝑥 2 𝑘𝐴
On pose :
ℎ𝑃 ℎ𝑃
𝑛2 = ⇒ 𝑛 = √ = √25 = 5 m−1
𝑘𝐴 𝑘𝐴
En posant 𝜃(𝑥) = 𝑇(𝑥) − 𝑇∞ , l’équation devient :
𝑑2𝜃
− 𝑛2 𝜃 = 0
𝑑𝑥 2
La solution générale est :
𝜃(𝑥) = 𝐶1 cosh(𝑛𝑥) + 𝐶2 sinh(𝑛𝑥)
Donc :
𝑇(𝑥) = 𝐶1 cosh(𝑛𝑥) + 𝐶2 sinh(𝑛𝑥) + 𝑇∞
Application des conditions limites
À 𝒙 = 𝟎: 𝑻(𝟎) = 𝑻𝑩
𝑇𝐵 = 𝐶1 + 𝑇∞ ⇒ 𝐶1 = 𝑇𝐵 − 𝑇∞
À 𝒙 = 𝑳: extrémité isolée
𝑑𝑇
∣ =0
𝑑𝑥 𝑥=𝐿
25
Or :
𝑑𝑇
= 𝜆[𝐶1 sinh(𝜆𝑥) + 𝐶2 cosh(𝜆𝑥)]
𝑑𝑥
À 𝑥 = 𝐿:
𝐶1 sinh(𝜆𝐿) + 𝐶2 cosh(𝜆𝐿) = 0 ⇒ 𝐶2 = −𝐶1 tanh(𝜆𝐿)
Après avoir effectué les calculs (voir document joint), le profil de température analytique
le long de la barre est donné par :
On trouve alors les températures analytiques aux nœuds :
Noeud Position (m) Température (°C) Température (K)
1 0 100 273.15
2 0.2 49.43 322.58
3 0.4 30.85 304
4 0.6 24.05 297.20
5 0.8 21.66 294.81
6 1 21.07 294.22
Équation discrétisée :
Conditions aux limites et coefficients :
26
27
On obtient :
Noeud Pi Qi
2 0.33 40
3 0.375 22.5
4 0.38 16.19
5 0.381 13.81
III. Algorithme et programme :
Le code :
28
29
Résultats analytiques et numériques :
Interprétation :
Les courbes numérique et analytique coïncident pratiquement en tout point, traduisant
une erreur négligeable et validant la fiabilité du modèle. Le profil thermique présente
une chute rapide près de la source, puis tend vers une stabilisation graduelle.
Même travail avec maillage=0.1m :
30
Même travail avec maillage=0.05m :
31
IV. Conclusion
L’étude a montré que la précision de la solution numérique dépend
fortement du pas de maillage utilisé. En raffinant le maillage, de Δ𝑥 =
0.2 mà Δ𝑥 = 0.05 m, on observe une amélioration significative de
l’approximation numérique par rapport à la solution analytique. Cette
convergence progressive confirme la stabilité et la robustesse de la
méthode numérique employée, en l’occurrence la méthode des différences
finies.
Le raffinement du maillage permet de mieux capturer les variations de
température le long de la barre, notamment près des extrémités où les
gradients sont les plus prononcés. Le pas le plus fin (Δ𝑥 = 0.05 m) offre
ainsi la meilleure précision, avec une erreur pratiquement négligeable, ce
qui souligne l’importance d’un maillage adapté pour obtenir des résultats
fiables.
Enfin, cette étude met en évidence l’efficacité des méthodes numériques
dans la résolution des problèmes de transfert de chaleur, en permettant de
32
traiter des configurations plus complexes, comme l’effet de la convection,
tout en fournissant des résultats comparables aux solutions analytiques
lorsque les paramètres de discrétisation sont correctement choisis.
33