‫ﺔﻌﻣﺎﺠ ﺲﻧﻭﺗ ﺭﺎﻧﻣﻟﺍ‬

Université de Tunis El Manar

Département Génie Mécanique

Projet de fin d’année

Etude par approche CFD de fonctionnement d’un

Échangeur thermique

Réalisé par :
Attia Med Amine
Haouala Khalil
Classe :
2 AGM 3

Encadré par :
Mr Kanfoudi Hatem
Mr Bellkhal Ghazi
Année universitaire 2017/2018
 Remerciements
Nous tenons à remercier toute personne qui a contribué à l’élaboration de mon projet et qui
m’a aidé durant sa rédaction.

En premier lieu, nous tenons à remercier notre encadrant Mr « KANFOUDI Hatem »

Professeur sous le département de l’hydraulique à l’ENIT, qui ne nous a jamais privé de son
aide en nous faisant bénéficier de ses idées et ses conseils et essentiellement son bagage
impressionnant et qui nous a offert un climat de travail favorable.

Nous souhaitons également adresser un remerciement très spécial à Mr « BELKHALL Ghazi

», pour ses conseils et ses idées qui nous ont permis d’avancer énormément dans
notre projet.

Enfin, nous tenons également à remercier chaleureusement nos familles pour leur soutien et
leur support.

2
 Résumé

Mots clefs: transfert thermique, modélisation, simulation, ANSYS CFD.

3
 Table des matières

Remerciements......................................................................................................................................2
Table des matières.................................................................................................................................4
Liste des figures......................................................................................................................................6
Liste des tableaux...................................................................................................................................7
Introduction générale.............................................................................................................................7
1. Chapitre 1 : Etude bibliographique.................................................................................................9
1.1. Introduction...........................................................................................................................10
1.2. Généralités sur les échangeurs de chaleur............................................................................10
1.2.1. Définition générale.........................................................................................................10
1.2.2. Principe de fonctionnement...........................................................................................10
1.3. Critères de classification des échangeurs de chaleur.............................................................12
1.3.1. Classification selon la technologie..................................................................................12
1.3.2. Classification selon le mode de transfert de chaleur......................................................12
1.3.3. Classification suivant le procédé de transfert de chaleur...............................................13
1.3.4. Classification fonctionnelle.............................................................................................13
1.3.5. Classification suivant la nature du matériau de la paroi d’échange................................13
1.4. Principaux types d’échangeur de chaleur..............................................................................13
1.4.1. Echangeurs tubulaires....................................................................................................13
1.4.2. Echangeurs à plaques.....................................................................................................15
1.4.3. Autres types....................................................................................................................15
1.5. Hypothèse de fonctionnement..............................................................................................20
1.6. MÉTHODES DE CALCUL ANALYTIQUE....................................................................................20
1.6.1. MÉTHODES DE CALCUL ANALYTIQUE.............................................................................20
1.6.2. Détermination du coefficient d'échange global..............................................................20
1.7. Les équations fondamentales................................................................................................21
1.7.1. Equations de Navier-Stokes............................................................................................21
1.7.2. Equation de la chaleur....................................................................................................22
Conclusion............................................................................................................................................25
2. Chapitre 2 : Modélisation Numérique..........................................................................................26
2.1. Introduction...........................................................................................................................27

4
 2.2. Organisation du code CFD : Computational Fluid Dynamics..................................................27
2.2.1. Préprocesseurs...............................................................................................................29
2.2.2. Solveur............................................................................................................................30
2.2.3. Post-processeur..............................................................................................................30
2.3. Les codes de CFD...................................................................................................................30
2.3.1. Le code ANSYS CFX.........................................................................................................30
2.3.2. Le code ANSYS Fluent.....................................................................................................31
2.4. La méthode numérique.........................................................................................................31
2.4.1. Procédure numérique.....................................................................................................32
2.4.2. Maillage..........................................................................................................................32
Conclusion............................................................................................................................................35
3. Chapitre 3 : Modélisation, résultats et interprétation..................................................................36
3.1. Introduction :.........................................................................................................................37
3.2. Construction du modèle :......................................................................................................37
3.2.1. Démarche numérique par CFD :.....................................................................................37
3.3. Simulation numérique :.........................................................................................................43
3.3.1. Observations et interprétations :....................................................................................43
3.3.2. Etude de l’efficacité........................................................................................................48
4. Conclusion :..................................................................................................................................54
Conclusion générale.............................................................................................................................55
Références bibliographiques................................................................................................................58

5
 Liste des figures

Figure 1:Les différents types d'échanges..............................................................................................10

Figure 2 : différents types d’écoulements .........................................................................................11
Figure 3: divers types d'échangeurs tubulaires....................................................................................14
Figure 4 : échangeur tubulaire à tube graphite ...................................................................................15
Figure 5 : échangeur à bloc en graphite ..............................................................................................15
Figure 6 : échangeur à plaque en graphite ..........................................................................................16
Figure 7 : les avantages et les inconvénients des échangeurs en polymère.........................................17
Figure 8 : échangeur de chaleur en verre.............................................................................................18
Figure 9 : maillage structuré.................................................................................................................30
Figure 10 : maillage non structuré........................................................................................................31
Figure 11 : Adaptation de la géométrie (solides)..................................................................................36
Figure 12 : Adaptation de la géométrie (Liquides)................................................................................36
Figure 13 : Construction du maillage (premier essai)...........................................................................37
Figure 14 : Construction du maillage (essai final).................................................................................38
Figure 15 :definition du domaine d'ecoulement..................................................................................39
Figure 16 : Variation de la température pour le système étudié..........................................................42
Figure 17 : Variation de la pression pour le fluide intérieur.................................................................43
Figure 18 : Variation de la pression pour le fluide extérieur.................................................................43
Figure 19 : Profil de vitesse avant l'établissement du régime laminaire..............................................44
Figure 20 : Profil de vitesse en régime laminaire..................................................................................45

6
 Liste des tableaux

Introduction générale

7
8
1. Chapitre 1 : Etude bibliographique

9
 1.1. Introduction 
Dans notre environnement quotidien direct et dans une multitude de procédés industriels
nous
sommes amenés à rencontrer des échangeurs thermiques. Un exemple assez bien connu
d'échangeur de chaleur est le "radiateur" d'une automobile, qui permet de refroidir le moteur.
Ce sont des appareils qui servent à transférer de la chaleur entre deux fluides de températures
différentes.

1.2. Généralités sur les échangeurs de chaleur 

1.2.1. Définition générale

Comme leur nom l'indique, les échangeurs thermiques sont des dispositifs destinés à
favoriser l'échange thermique entre deux fluides, sans que les fluides ne se mélangent. 1
[CITATION And94 \l 1036 ]

Le principe le plus général consiste à faire circuler deux fluides à travers des conduits qui les
mettent en contact thermique. De manière générale, les deux fluides sont mis en contact
thermique à travers une paroi qui est le plus souvent métallique ce qui favorise les échanges
de chaleur. 2[ CITATION PTo02 \l 1036 ]

1.2.2. Principe de fonctionnement 

Il existe de nombreuses technologies d’échangeurs de chaleur, mais tous

fonctionnent suivant les mêmes processus physiques d’échanges, c’est-à-dire :

- La conduction qui représente les échanges à travers les parois (le plus souvent
métalliques),

- La convection qui représente les échanges entre les fluides et les parois,

- Le rayonnement qui représente les échanges radiatifs entre les fluides et les parois

(principalement infra rouge) bien que ce dernier soit souvent négligeable (car il est pris
en compte uniquement pour un fonctionnement à haute température).

On distingue également deux

1
: P. Tobaly,≪ LES ECHANGEURS DE CHALEUR≫, France, 2002
2
 : BONTEMPS et al,≪ DESCRIPTION DES ECHANGEURS≫, Université Joseph Fourier, Institut
universitaire de Technologie, Département Génie thermique et Énergie (Grenoble), B 2 341

10
 modes d’écoulement différents :

– Celui à co-courants : écoulements parallèles des fluides et dans la même direction ;

– Celui à contre-courants : écoulements parallèles des fluides mais dans des directions
inverses ;

Le mode de circulation influe également sur la qualité et l’efficacité du transfert. Un

échangeur dit anti-méthodique (ou à co-courants) signifie que l’entrée des deux fluides (chaud
et froid) se situe du même côté de l’échangeur. La configuration dite méthodique est celle ou
les fluides circulent à contre-courants. On l’appelle méthodique car c’est cette configuration
qui permet d’obtenir les meilleures performances d’échangeurs. En effet, pour une circulation
à co-courant des fluides, la température de sortie du fluide chaud ne pourra pas être plus basse
que la température de sortie du fluide froid tandis qu’à contre-courant, il est possible
d’abaisser la température du fluide chaud à une valeur qui tendra vers la température d’entrée
du fluide froid. Des représentations graphiques illustreront plus loin ces distinctions. 3
[ CITATION WEB15 \l 1036 ]

Les deux figures ci-dessous récapitulent les différents types d’échanges thermiques ainsi que
les modes d’écoulements.

Figure 1 : Les différents types

4 d'échanges

[ CITATION DEC17 \l 1036 ]

3
 : https://www.encyclopedie-energie.org/les-echangeurs-de-chaleur/
4
 : DECLAYE Sébastien GENDEBIEN Samuel LEMORT Vincent ≪P 31 LES ÉCHANGEURS DE
CHALEUR≫., disponible à l’adresse : https://docplayer.fr/24972780-P-31-les-echangeurs-de-chaleur-declaye-
sebastien-gendebien-samuel-lemort-vincent.html

11
 5
Figure 2 : différents types d’écoulements
[ CITATION Adi08 \l 1036 ]

1.3. Critères de classification des échangeurs de chaleur

Il existe plusieurs critères de classification des échangeurs de chaleur dont les plus
importants sont cités ci-dessous 6:[ CITATION Zoé14 \l 1036 ]

1.3.1. Classification selon la technologie

Les principaux types d’échangeurs rencontrés sont les suivants :

 A tubes : monotubes, coaxiaux ou multitubulaires ;

 A plaques : à surface primaire ou à surface secondaire ;

 Autres types : contact direct, à caloducs ou à lit fluidisé, à un seul fluide...

1.3.2. Classification selon le mode de transfert de chaleur

Les trois modes de transfert de chaleur (conduction, convection, rayonnement) sont couplés
dans la plupart des applications (chambre de combustion, récupération sur les fumées, etc.) ; il
y’a souvent un mode de transfert prédominant. Pour tout échangeur avec transfert de chaleur à
travers une paroi, la conduction intervient.

5
 : Adil Ridha (université de Caen) Transferts de chaleur et de masse≫ 2008/2009
6
 : technique

12
1.3.3. Classification suivant le procédé de transfert de chaleur

Suivant qu’il y a ou non stockage de chaleur, on définit un fonctionnement en

récupérateur ou en régénérateur de chaleur :

 Transfert sans stockage : en récupérateur, avec 2 ou n passages et un écoulement en

général continu ;

 Transfert avec stockage : en régénérateur, avec un seul passage et un écoulement

intermittent, la matrice de stockage étant statique ou dynamique.

1.3.4. Classification fonctionnelle

Le passage des fluides dans l’échangeur peut s’effectuer avec ou sans changement de
phase ; suivant le cas, on dit que l’on a un écoulement monophasique ou diphasique. On
rencontre alors les différents cas suivants :

 Les deux fluides ont un écoulement monophasique ;

 Un seul fluide à un écoulement avec changement de phase, cas des évaporateurs ou

des condenseurs ;

 Les deux fluides ont un écoulement avec changement de phase, cas des évapo-
condenseurs.

1.3.5. Classification suivant la nature du matériau de la paroi d’échange

On retiendra deux types de parois :

 Les échangeurs métalliques en acier, cuivre, aluminium ou matériaux ;

 Spéciaux : superalliages, métaux ou alliages réfractaires ;

 Les échangeurs non métalliques en plastique, céramique, graphite, verre, etc.

1.4. Principaux types d’échangeur de chaleur

1.4.1. Echangeurs tubulaires

Pour des raisons historiques et économiques, les échangeurs utilisant des tubes comme
constituant principal de la paroi d’échange sont les plus répandus. On peut distinguer trois
catégories en fonction du nombre de tubes et de leur arrangement :[ CITATION Zoé14 \l
1036 ]

13
 Échangeur monotube (Figure 3a), dans lequel le tube est placé à l’intérieur d’un
réservoir et a généralement la forme d’un serpentin ou d’un drapeau ;

 Échangeur bitube ou coaxial (Figure 3b), dans lequel les deux tubes sont
concentriques et l’échangeur est souvent cintré pour limiter son encombrement ; en
général, le fluide chaud ou le fluide à haute pression s’écoule dans le tube intérieur ;

 Échangeur multitubulaire, existant sous plusieurs formes :

- Échangeur à tubes séparés (Figure 3c) : à l’intérieur d’un tube de diamètre

suffisant (de l’ordre de 100 mm) se trouvent placés plusieurs tubes de petit
diamètre (8 à 20 mm) maintenus écartés par des entretoises. L’échangeur peut
être soit rectiligne, soit enroulé,

- Échangeur à tubes rapprochés (Figure 3d) : pour maintenir les tubes et obtenir un
passage suffisant pour le fluide extérieur au tube, on place un ruban ou un fil
enroulé en spirale autour de certains d’entre eux. Les tubes s’appuient les uns sur
les autres par l’intermédiaire des rubans. Les tubes peuvent être soit lisses, soit
ailetés (figure 3e) : ces derniers permettent d’améliorer le coefficient d’échange
thermique.

- Échangeur à tubes et calandre (figure 3f) : le nombre de tubes varie alors de

quelques dizaines à plusieurs milliers ; c’est l’échangeur actuellement le plus
répandu.

14
 -

Figure 3: divers types d'échangeurs tubulaires

[ CITATION Zoé14 \l 1036 ]

Dont les avantages cités ci-dessous :

Les échangeurs de chaleur tubulaires sont des échangeurs extrêmement robustes et fiables.

C’est là un des principaux avantages de la technologie, car ils sont virtuellement utilisables
dans toutes les applications, et capables de résister à des températures et pressions extrêmes
tout en acceptant des écarts de température importants.

Or il a certains inconvénients, on cite parmi l’encombrants puisqu’Ils peuvent mesurer

jusqu’à 40 mètres de long et 5 mètres de large, et prennent donc beaucoup de place. Aussi ils
ne sont pas très performants. Du fait de la faible turbulence et de la surface d’échange
moindre (par rapport à un échangeur à plaques), la performance d’un échangeur tubulaire est
plutôt faible.

1.4.2. Echangeurs à plaques

1.4.3. Autres types

1.4.3.1. Échangeurs en graphite

Il existe différentes versions de ces échangeurs : à faisceau tubulaire, à bloc et à plaques.

Les échangeurs de chaleur en graphite sont destinés au chauffage et refroidissement de

fluides particulièrement corrosifs tels que des acides minéraux (Acide chlorhydrique,
sulfurique, ...).

Le graphite offre, outre une exceptionnelle résistance à la corrosion, de bonnes propriétés

mécaniques et une bonne conductivité thermique favorable au transfert de chaleur et à la
résistance aux chocs thermiques.

Le graphite est naturellement poreux et pour être apte à confectionner des échangeurs de
chaleur il doit être imprégné. Les résines d'imprégnation sont variées et dépendent des
fabricants. Elles déterminent la tenue en température et la résistance à la corrosion. La
température d'utilisation maximum est au moins de 180°C, mais peut être supérieure selon les
fabricants.

15
Si le graphite résiste parfaitement aux acides minéraux réducteurs, il peut par contre être
sensible aux oxydants (acide nitrique, sulfurique, chlore, ...) 7[CITATION Ech \l 1036 ]

Il existe différentes versions de ces échangeurs : à faisceau tubulaire, à bloc et à plaques.

Figure 4 : échangeur tubulaire à tube graphite 8

[ CITATION ECH \l 1036 ]

Figure 5 : échangeur à bloc en graphite 9

7
 : http://processs.free.fr/Pages/VersionMobile.php?page=4528#Voir_aussi
8
 :https://www.mersen.com/fr/produits/equipement-anticorrosion/echangeurs-de-chaleur/echangeurs-de-
chaleur-tubulaires-en-graphite
9
 : http://processs.free.fr/Pages/VersionMobile.php?page=4528#Voir_aussi

16
 [ CITATION Ech \l 1036 ]

Figure 6 : échangeur à plaque en graphite 10

[ CITATION Zoé14 \l 1036 ]

1.4.3.2. Échangeurs en céramique

Pour des gaz à haute température (supérieure à 1 000 oC), les échangeurs de chaleur
métalliques ne peuvent plus être utilisés. La figure 32 donne les plages d'utilisation des
échangeurs de chaleur en fonction de la température.

Les principaux avantages des échangeurs en céramique se résument ainsi :

 Une bonne tenue à haute température (supérieure à 1 000  C) ;

o

 Une bonne conductivité thermique ;

 Une bonne résistance à la corrosion ;

 Une masse volumique faible (3 000 kg/m ). 3

10
 : technique

17
 Les inconvénients par rapport aux échangeurs métalliques sont une certaine fragilité et une
dispersion de leurs propriétés mécaniques.

1.4.3.3. Échangeurs en polymère

Les premiers échangeurs de chaleur en matière plastique ont été fabriqués dans les années
1960 par DuPont. 11[ CITATION Jua12 \l 1036 ]

La figure ci-dessous récapitule les avantages et les inconvénients des échangeurs en

polymères.

Figure 7 : les avantages et les inconvénients des échangeurs en polymère12

[ CITATION Zoé14 \l 1036 ]

1.4.3.4. Échangeurs en verre

Les divers avantages de l'utilisation du verre comme surface d'échange sont les suivants :

Une bonne tenue à la corrosion notamment sulfurique, ce qui permet son utilisation
industrielle dans des problèmes de récupération sur des gaz à des températures inférieures à la
température de condensation des acides ; une facilité de nettoyage due à l'état de surface du
verre qui limite l'accrochage des dépôts et qui permet de retrouver l'état initial après un simple
nettoyage ; un prix peu élevé.
11
 : CEVALLOS(J.G), BERGLES (A.E.),BARCOHEN(A),RODGERS(P.),GUPTA(S.K.)-Polymer heat
exchangers-history, opportunités, and challenges.
12
 : technique

18
 Les inconvénients sont essentiellement la fragilité du matériau et la conductivité thermique
du verre qui est très inférieure à celle des métaux. De ce fait, ce matériau trouve surtout ses
applications sur des gaz, où la résistance thermique se trouve essentiellement côté fluide
gazeux.

Figure 8 : échangeur de chaleur en verre13

[ CITATION Con17 \l 1036 ]

1.5. Hypothèse de fonctionnement

Les hypothèses retenues pour modéliser un échangeur dapres [ CITATION Rou11 \l 1036 ]
sont :14

13
 :https://www.buchiglas.com/fr/produits/reacteurs-en-verre-unite-pilote-et-production/accessoires-
composants-buchiflex/condenseurs-echangeurs-de-chaleur.html
14
 : E. ROULAND (université de ROUEN) : les échangeurs thermiques 2010/2011, disponible sur l’adresse
http://gsi-energie.univ-rouen.fr/IMG/pdf/cours-iup-me-echangeurthermique-2.pdf

19
  Régime permanent

 Caractéristiques des fluides : masse volumique, viscosité, conductivité thermique et

les coefficients d'échange thermique sont des constantes

 Pas de pertes thermiques : Le transfert s’effectue au travers d’une paroi

 Le transfert thermique ne s’effectue que par convection et conduction

 Pas de changement de phase au cours du transfert

1.6. MÉTHODES DE CALCUL ANALYTIQUE

1.6.1. MÉTHODES DE CALCUL ANALYTIQUE

Plusieurs types de calcul peuvent être envisagés. Les deux principaux sont les suivants :

 La détermination de la surface d'échange connaissant la puissance échangée Φ et les

températures d'entrée et de sortie des deux fluides (fluide chaud, T 1e, T 1s ; fluide
froid, T 2e, T 2s ; ). C'est la méthode de DTML, différence de température moyenne
logarithmique. Cette méthode est particulièrement bien adaptée au mode
dimensionnement ;

 La détermination des températures de sortie des fluides, connaissant leurs

températures d'entrée et la surface d'échange. C'est la méthode de l'efficacité. Elle est
bien adaptée au mode évaluation.

Les deux méthodes peuvent être utilisées simultanément en particulier lorsqu'on désire faire
varier certaines grandeurs.

1.6.2. Détermination du coefficient d'échange global

Le flux thermique transféré du fluide chaud au fluide froid est donné par l’expression

Avec :

h c: Coefficient de convection du fluide chaud

20
 h f : Coefficient de convection du fluide froid

1
: La résistance thermique convective du fluide chaud
hc S

1
hf S
: La résistance thermique convective du fluide froid

e : l’épaisseur de la paroi

e
: La résistance thermique de conduction
λS

Soit encore :

∅=K ( T c +T f ) S

Et donc :

1.7. Les équations fondamentales

1.7.1. Equations de Navier-Stokes

En se basant sur ces hypothèses, on peut établir les équations de Navier-Stokes à l’aide des
équations de conservation : masse, la quantité de mouvement, l’énergie.

Equation de continuité :

Appelée aussi équation de conservation de la masse :

∂U i
=0
∂ xi

Avec U i est la composante spatiale de la vitesse.

Equation de conservation de la quantité de mouvement :

Appelée, également, équation de Navier-Stokes. Elle s’écrit de la manière suivante :

∂U i ∂ Ui 1 ∂p ∂ ² Ui
+Uj =- +v + Si
∂t ∂ xj ρ ∂ xi ∂² xj

Avec :

21
v: Viscosité cinématique

∂U i
 : Terme de variation temporelle
∂t

∂U i
Uj  : Terme de transport convectif
∂ xj

Si :Regroupe les forces extérieures agissant sur le domaine d’étude

1 ∂p
 : Variation de pression
ρ ∂ xi

∂ ² Ui
v  : Le frottement visqueux
∂² xj

1.7.2. Equation de la chaleur

On peut établir l’équation da la chaleur en utilisant l’équation du bilan de l’énergie interne et

les lois suivantes :

 Loi de comportement rhéologique d’un fluide Newtonien

 Loi de Fourrier

 Loi constitutive de fermeture de l’énergie interne déduite de la thermodynamique

Loi de comportement rhéologique d’un fluide Newtonien :

La loi de comportement linéaire de Newton qui s’écrit sous sa forme la plus générale :

2
τ =2 μD+( μ ¿ ¿ θ− μ)( ∇ . v ) I ¿
3

Avec :
τ : Le tenseur des taux de cisaillement

D : le tenseur des taux de déformation

μ: la viscosité dynamique

μθ : La viscosité volumique ou deuxième viscosité rencontrée dans le cas des fluides

compressibles.

22
 On utilise l’équation de conservation de la masse pour réduire cette loi dans le cas d’un
fluide incompressible à l’expression suivante :

τ =2 μD

Loi de Fourier :

Le vecteur courant de chaleur q pour un fluide newtonien est écrit sous la forme :

q = -λ ∇ T

Avec :

λ : La conductivité thermique (Wm−1 k−1)

Loi constitutive de fermeture de l’énergie interne déduite de la thermodynamique :

À l’aide de la relation de Gibbs :

de ds
ρ =ρT −ρ(∇ . v)
dt dt

Comme le fluide est incompressible l’équation se réduit à :

de= T ( dSdt ) dT=c ( T ) dT

v

Où Cv est la chaleur spécifique à volume constant.

On remplace cette relation dans l’équation de l’énergie interne on trouve :

ρ cv ( ∂∂Tt +v . ∇ T )=−∇ q+ ρr +τ :( ∇ v )

Et à l’aide de la loi de Fourier on trouve :

ρ cv ( ∂∂Tt +v . ∇ T )=−λ ∆ T + ρr+ μ [ ∆ v +( ∆ v ) ] :( ∇ v )

t

Avec :

r: Densité massique du rayonnement (Wm−2)

𝜆 : Conductivité thermique (Wm−1 k−1)

23
𝐶v : Chaleur spécifique à volume constant (𝐽kg−1 K −1)
∂T
: variation temporelle
∂t

v . ∇ T : le transport convectif

μ [ ∆ v+ ( ∆ v )t ] :(∇ v): le terme source

Conclusion

Dans ce chapitre, on a présenté le contexte général de notre projet en réalisant une synthèse
bibliographique concernant le phénomène de transfert de chaleur qui se produit au sein d’un
échangeur thermique ainsi que les types, les critères de classifications des échangeurs et les
équations fondamentales pour les études des échangeurs.

24
25
2. Chapitre 2 : Modélisation
Numérique

26
 2.1. Introduction
De nos jours, la simulation numérique est devenue nécessaire dans la résolution des
problèmes mathématiques insolubles analytiquement ou les problèmes non linéaires.
Dans notre cas, on veut résoudre un problème en mécanique de fluide en cherchant à
gagner le temps et le coût de travail, pour cela on a utilisé le logiciel de dynamique de
fluide ANSYS CFX dans une première partie pour simplifier notre problème, et dans
une deuxième partie, on a travaillé avec le code ANSYS Fluent puisqu’on a trouvé cette
méthode plus simple à exécuter.

D’abord, ANSYS est un logiciel de modélisation par éléments finis permettant de

résoudre numériquement une grande variété de problèmes mécaniques. Il peut effectuer
une analyse structurelle (linéaire et non linéaire), une analyse thermique, des problèmes
électromagnétiques, une analyse de fluide, etc. 15[CITATION lea \l 1036 ]

Nous allons présenter la structure du code ANSYS CFX en décrivant le

fonctionnement des outils du logiciel de simulation et nous allons définir le code
ANSYS Fluent.

2.2. Organisation du code CFD : Computational Fluid

Dynamics

Le terme « CFD » est une abréviation de « Computational Fluid Dynamics », qui

définit une approche de mécanique de fluide numérique et présente des codes
numériques permettant d’analyser les équations mathématiques pour un problème de
mécanique de fluide et les résoudre dans un espace bien déterminé, défini par un
maillage structuré ou non structuré.

Ce logiciel nous permet de suivre l’évolution des paramètres instantanés définissant

notre problème tels que la vitesse, la pression et la température dans chaque point du
maillage où se fait le calcul.

En fait, la qualité du maillage intervient dans la limitation du temps de calcul et la

précision du résultat présenté, c’est-à-dire pour un bon maillage le temps de calcul est
optimisé et le résultat est précis.

15
 : http://fulllibrary.com/Mechanical/ANSYS/Learn-Ansys.html
 Le CFD présente aujourd’hui un logiciel de simulation numérique assez efficace pour
la résolution des problèmes numériques pour plusieurs domaines, tels que la recherche
scientifique et l’industrie.

En général, la résolution d’un problème de mécanique de fluide par ce logiciel passe

pour les étapes décrites ci-dessous :

C’est la première La deuxième C’est la dernière

Définition du problème

Exploitation des résultats

Résolution numérique du problème
étape qui permet étape permet de étape qui permet
de définir la résoudre le de visualiser les
géométrie étudiée problèmeen résultats de calcul
et de générer le évaluant la correspondants
maillage adéquat convergence du aux différents
ainsi que les problème défini paramètres définis
conditions aux en exploitant les
limites et initiaux différents outils de
afin de procéder CFD
au calcul.

Le schéma suivant présente la structure de CFD :

Préprocesseur

CFD Solveur

Post-
processeur

2.2.1. Préprocesseurs
 . Dans le préprocesseur, le modèle est configuré. Il comprend un certain nombre
d'étapes et généralement dans l'ordre suivant :

 Construire une géométrie. Selon que la géométrie du problème est un, deux ou
En trois dimensions, la géométrie consiste à créer des lignes, des zones ou des
volumes. Ces géométries peuvent ensuite, si nécessaire, être utilisées pour créer
d’autres géométries en l'utilisation d'opérations booléennes. L'idée clé lors de la
construction de la géométrie comme c'est pour simplifier la génération du
maillage de l'élément. Par conséquent, cette étape est facultative mais le plus
souvent utilisé. Des nœuds et des éléments peuvent cependant être créés à partir
de coordonnées seulement.

 Définir les matériaux. Un matériau est défini par ses constantes de matériau.
Chaque l'élément doit être affecté à un matériau particulier.

 Générer un maillage d'élément. Le problème est discrétisé avec des points

nodaux. les nœuds sont connectés pour former des éléments finis, qui forment
ensemble le matériau le volume. En fonction du problème et des hypothèses
retenues, le type d'élément doit être déterminé. Les types d’éléments communs
sont les suivants : treillis, poutre, plaque, coquille et éléments solides. Chaque
type d'élément peut contenir plusieurs sous-types, par exemple Éléments solides
2D à 4 nœuds, éléments solides 3D à 20 nœuds. Par conséquent, soin doit être
pris lorsque le type d'élément est choisi.

Dans ANSYS, le maillage des éléments peut être créé de plusieurs manières. Le plus
moyen habituel est qu'il soit créé automatiquement, mais plus ou moins contrôlé. Par
exemple, vous pouvez spécifier un certain nombre d’éléments dans une zone spécifique,
ou vous pouvez forcer le générateur de maillage à conserver une taille spécifique de
l'élément dans une zone. Certaines formes ou tailles d’éléments ne sont pas
recommandé et si ces limites sont violées, un avertissement sera généré dans

ANSYS. Il appartient à l'utilisateur de créer un maillage capable de générer des

résultats avec un degré suffisant de précision.16[ CITATION INT \l 1036 ]

2.2.2. Solveur

16
 : http://innomet.ttu.ee/martin/mer0070/wb/ws1/ansys_classic_tutorial_lab1_.pdf
 Le rôle principal du solveur consiste à exécuter le calcul pour un problème étudié. Les
résultats sont générés par la méthode des volumes finis en se basant sur les équations de
Navier Stockes dans chaque maille. Pour notre cas, on s’est basé également sur
l’équation de la chaleur en adoptant le bilan énergétique pour étudier le transfert
thermique dans notre système.

Le solveur adopte une solution approchée pour la résolution du problème vu la non

linéarité de système d’équations, et pour aboutir au résultat final, le système doit
converger, c’est-à-dire la différence entre deux itérations successives sera inférieure à
un critère défini par l’utilisateur, si ceci est achevé, le calcul s’arrête automatiquement,
sinon le calcul s’arrête lorsqu’il atteint le nombre d’itérations maximal désigné.

Une fois le calcul est terminé, le programme enregistre les informations nécessaires
pour les transmettre au Post-processeur.

2.2.3. Post-processeur 

Après avoir terminé le calcul, CFD-Post permet d’analyser et de traiter les résultats
observés.

Cet outil nous permet également de visualiser les profils de vitesse, d’étudier la
variation de pression ou de température selon le problème étudié, ainsi que les lignes de
courant, les contours et les iso-surfaces pour chaque variable, on peut aussi observer une
animation pour mieux comprendre le phénomène décrit.

2.3. Les codes de CFD

2.3.1. Le code ANSYS CFX

ANSYS CFX est un code industriel de la CFD, ce code est l’un des plus utilisés
actuellement pour la modélisation numérique des problèmes de mécanique de fluide, il
permet de traiter ces problèmes en exploitant les différents schémas numériques
présentés.

Le traitement des problèmes numériques en CFX est effectué dans le but d’optimiser le
résultat qu’on veut trouver, c’est-à-dire un résultat précis pour une durée de réalisation
limitée.
 Il est aussi compatible avec les systèmes CAO, c’est-à-dire on peut extraire les
données des logiciels comme CATIA par exemple pour l’exploiter pour un problème de
dynamique de fluide exécuté par ANSYS Fluent.

La modélisation des écoulements à l’aide de CFX se fait en trois étapes principales :

 La réalisation du maillage en définissant les conditions initiales et aux limites,

ceci est réalisé à l’aide de (CFX-Pré)

 La résolution des équations à l’aide de (CFX-Solver)

 Le traitement des résultats à l’aide de (CFX-Post)

2.3.2. Le code ANSYS Fluent

Ce code est l’un des codes les plus performants comme ANSYS CFX, et il offre toutes
les capacités physiques qui permettent de résoudre n’importe quel problème de
dynamique de fluide, les écoulements des fluides, les transferts thermiques et même les
réactions chimiques.

Il est utilisé pour des applications industrielles, comme la modélisation de l’action du

vent sur un structure ou l’action de l’eau sur une hélice et même la simulation d’un
transfert de chaleur pour un échangeur thermique, qui est le cas étudié ici.

De même que ANSYS CFX, il est aussi compatible avec les logiciels CAO ce qui
justifie son efficacité et sa performance pour l’étude des projets de grande importance.

2.4. La méthode numérique

Plusieurs phénomènes de la mécanique du fluide peuvent être modélisés par des
équations différentielles et non linéaires.

Ces équations sont difficiles à résoudre donc on utilise la méthode numérique pour
transformer une équation différentielle en une équation algébrique facile à résoudre
(équation discrétisée).

Les méthodes de discrétisation les plus couramment utilisées sont :

 La méthode des différences finies.

  La méthode des éléments finis.

 La méthode des volumes finis.

2.4.1. Procédure numérique

Dans ce projet la méthode utilisée par le code ANSYS Fluent est la méthode des
volumes finis.

Pour résoudre les équations de Navier-Stokes par la méthode des volumes finis, on
divise le domaine en volumes de contrôle discrète en utilisant une grille de calcul
(maillage) et pour chacun de ces volumes on intègre le flux de masse, de quantité de
mouvement et d’énergie sur toutes ses faces afin de construire les équations algébriques
pour les variables discrètes dépendantes, les inconnues, telle que la vitesse, la pression
et la température. Il est alors possible de résoudre un système matriciel pour obtenir les
variables physiques de chaque volume.

En régime permanent, la forme linéaire des équations discrètes est la suivante :

Avec :

a p et a np sont des coefficients, ∅p est la valeur d’une variable scalaire au centre de la

cellule et 𝑛𝑏 représente les indices des cellules voisines.

2.4.2. Maillage

Le maillage est la première étape de la méthode des éléments finis. Cette étape consiste
à discrétiser l’espace en un ensemble de petits volumes.

3 types de maillage existent :

2.4.2.1. Maillage structuré

Un maillage structuré est un maillage régulier, c’est un maillage utilisant des éléments
de type hexaédriques en volume, quadrilatère (carrée) en surface.

Ce type de maillage est caractérisé par l’économie en nombre d'éléments, il présente

un nombre inférieur de mailles par rapport à un maillage non structuré équivalent.
 Ce type de maillage est difficile à générer dans le cas d'une géométrie complexe. De
plus, ce type de maillage pose des problèmes au niveau du raffinement, car le
raffinement d’une zone du maillage pourrait influencer l’entièreté du maillage, ce qui
peut induire à d’autres problématiques.

Figure 9 : maillage structuré

2.4.2.2. Maillage non structuré

C’est un maillage utilisant des éléments de type tétraèdres, prismes, d'hexaèdres et de

pyramides en volume, type triangle ou encore quadrilatère en surface…

Ce maillage a l’avantage de créer des maillages triangulaires ou tétraédriques dans des

géométries quelconques et complexes.

Ce type de maillage pose des difficultés pour contrôler la densité des points dans une
zone précise.
 Figure 10 : maillage non structuré
2.4.2.3. Maillage hybride

Maillage généré par un mélange d'éléments de différents types, triangulaires ou

quadrilatéraux en 2D, tétraédriques, prismatiques, ou pyramidaux en 3D. Il combine les
avantages du maillage structuré et ceux du maillage non structuré.
 Conclusion 

Dans ce chapitre, on a introduit la méthode de simulation numérique en expliquant les

différentes étapes établies pour la réalisation d’un projet par le logiciel de modélisation
CFD et les différents codes utilisés.

On a parlé des différents outils utilisés pour l’obtention d’une géométrie, la réalisation
d’un maillage et l’obtention des résultats.

De plus, on a donné les différents types de maillage ainsi que les méthodes utilisées
pour sa réalisation.

Ce chapitre explique de façon générale les étapes de réalisation d’un projet numérique,
maintenant on va procéder pour expliquer notre travail
3. Chapitre 3 : Modélisation,
résultats et interprétation
 3.1. Introduction :
Dans ce dernier chapitre, nous allons décrire les différentes étapes établies pour la
réalisation de notre projet, par suite nous allons exploiter les résultats trouvés par CFD-
Post pour interpréter le phénomène étudié.

Nous allons conclure le comportement de notre système (échangeur thermique) en

visualisant les variations des différents paramètres (température et pression) pour un
écoulement de deux fluides dans deux tubes coaxiaux.

Dans cette partie nous allons voir deux cas différents, les résultats sont établis en
régime laminaire pour deux fluides visqueux, dans le premier cas nous allons fixer un
débit constant pour le fluide froid et nous allons varier le débit du fluide chaud, dans le
deuxième cas nous allons varier le débit du fluide chaud et faire varier celui du fluide
froid.

3.2. Construction du modèle :

3.2.1. Démarche numérique par CFD :

Maintenant, nous allons décrire les étapes établies pour la résolution du problème
numérique.

Dans la première partie, nous allons décrire le problème physiquement, par suite nous
allons expliquer la méthode de génération du maillage, et finalement nous allons
introduire les conditions aux limites pour que le problème soit bien exposé.

On a utilisé le code ANSYS Fluent pour la réalisation de notre projet.

3.2.1.1. Définition du problème :

Le but principal de notre projet est d’étudier le comportement d’un système

comportant de deux fluides circulants dans un même sens dans deux tubes coaxiaux,
c’est-à-dire le premier fluide ayant une température élevée circule dans une conduite
entourée par un autre fluide froid qui circule dans une autre conduite de diamètre
supérieur au premier. Les deux fluides sont séparés par la conduite dans laquelle circule
le premier.

L’objectif est d’étudier l’efficacité de notre système pour un transfert de l’énergie

thermique du premier fluide vers le deuxième.
 Dans cette partie, nous allons présenter la simulation numérique de la convection pour
les deux fluides et la conduction de chaleur pour les deux tubes, pour cela nous allons
suivre les étapes suivantes afin de définir notre projet :

 Etude géométrique

 Construction du maillage

 Introduction des conditions initiaux et aux limites

3.2.1.2. Etude géométrique :

La première partie pour la résolution du problème consiste à créer une géométrie

adéquate pour définir les frontières qui limitent les domaines d’écoulement pour les
deux fluides, tout en adoptant une géométrie réduite afin de simplifier le calcul, on a
pris seulement une longueur égale à L=2.5 m pour optimiser le temps de calcul tout en
adoptant un maillage bien précis, puisque la réalisation d’une géométrie assez simple
permet d’obtenir un maillage convenable pour aboutir à des résultats efficaces.

Pour réaliser cette géométrie, on a créé 4 cercles coaxiaux qui définissent les 2 tubes.
Chaque tube est défini par 2 cercles et la différence de deux rayons de cercles successifs
désigne l’épaisseur de chaque tube.

Le premier tube a comme dimensions d1=16 mm pour le diamètre intérieur et d2=20

mm pour le diamètre extérieur, le deuxième est de diamètre intérieur d3=30 mm et de
diamètre extérieur d4=34 mm.

Après avoir dessiné les 4 cercles, on a extrudé la géométrie selon l’axe z pour obtenir
la forme désirée.

Finalement, nous avons définis les 4 corps qui constituent cette géométrie, 2 solides et
2 liquides
 Figure 11 : Adaptation de la géométrie (solides)
Cette figure présente les deux tubes.

Figure 12 : Adaptation de la géométrie (Liquides)

Cette figure présente le liquide intérieur chaud et le liquide extérieur froid colorés en
vert.

3.2.1.3. Construction du maillage :

Après avoir généré la géométrie désirée, on procède maintenant à la construction du

maillage pour la résolution du problème par la méthode des volumes finis.
 En effet, l’efficacité du calcul est dépendante de la qualité du maillage puisque le
maillage réalisé assure la convergence du problème en détectant les variations des
différents paramètres pour le domaine spécifié.

Pour notre problème on a adopté plusieurs maillages qui ont été testés pour choisir
finalement un maillage qui nous a paru le plus adéquat en termes de temps de calcul et
de précision.

Le maillage généré est non structuré pour le volume intérieur puisqu’on a conçu deux
tubes coaxiaux qui sont présentés par un ensemble de cercles successifs qui ont une
géométrie convenable à ce type de maillage.

Pour les surfaces latérales des cylindres coaxiaux, le maillage est structuré puisque la
géométrie est simple et peut être présentée par ce type de maillage.

On a défini ce cercle comme un polygone ayant 30 arêtes, ce nombre a été choisi après
plusieurs essais puisque les premiers essais avec un nombre d’arêtes inférieur à celui-ci
nous a donné des résultats mois efficaces qui ne permettent pas de bien évaluer le
comportement de notre système, alors que pour un maillage plus précis, le simulateur
exige plus de temps de calcul pour aboutir au résultat final, pour cela on a choisi ce
nombre d’arêtes.31

Cependant, dans le premier essai indiqué dans cette figure, on remarque que le
maillage réalisé au niveau du domaine correspondant au fluide froid n’est pas adéquat
pour obtenir des résultats précis.

Figure 13 : Construction du maillage (premier essai)

 Pour cela on a effectué un raffinement dans cette zone pour plus de précision.
Finalement, on a établi le maillage suivant.

Figure 14 : Construction du maillage (essai final)

On a obtenu un maillage comportant 1610304 nœuds et 449424 éléments pour
procéder au calcul numérique.

3.2.1.4. Conditions aux limites :

Une fois le maillage réalisé, on passe maintenant à l’introduction des conditions aux
limites pour définir le domaine d’écoulement et décrire l’état réel de notre projet.

Pour définir les conditions aux limites, on choisit d’abord les différentes sections qui
définissent l’entrée et la sortie de l’écoulement, et les différents obstacles qui désignent
le «wall» et sont représentés par des tubes dans lesquels circulent les fluides. Dans notre
cas on a nommé les sections attribuées à l’entrée et la sortie de l’écoulement pour le
fluide chaud simultanément « inner inlet » et « inner outlet », et on a défini celles du
liquide froid comme « inlet » et « outlet ».

Pour les tubes, on a nommé le tube intérieur « inner wall » et le tube extérieur « outer
pipe ».
 Figure 15 :definition du domaine d'ecoulement

Par suite, on définit les matériaux utilisés, on a choisi l’acier pour les deux conduites et
on a modifié les caractéristiques de deux fluides différents pour obtenir deux fluides
virtuels visqueux afin de créer un écoulement en régime laminaire. Le fluide intérieur
nommé « fluide virtuel1 » est plus visqueux que le deuxième nommé « fluide virtuel2 ».

Après avoir réalisé cette étape, on passe maintenant à l’introduction des différents
paramètres qui caractérisent l’écoulement.

On a défini la pression relative à la sortie qui est égale à zéro pour les deux fluides, la
vitesse à la paroi est aussi nulle.

Dans une première partie, on a fixé la vitesse au niveau de « inlet » à une valeur
constante v0=1 m/s et on a varié la vitesse de « inner inlet », et dans une deuxième
partie, on a fixé la vitesse au niveau de « inner inlet » à une valeur constante v1=0.6 m/s
et on a varié celle de « inlet ».

Pour la température, on a injecté un fluide chaud ayant une température égale à 340°K
à travers « inner inlet » et un fluide froid ayant une température égale à 290°K passant à
travers la section « inlet ».

3.2.1.5. Critère de convergence du problème :

 On peut déterminer la convergence d’un calcul lorsque l’erreur tend à décroître après
chaque itération.

La figure suivante donne l’évolution du résidu qui représente la somme d’erreurs sur
toutes les mailles pour les équations de continuité et d’énergie et les différentes
composantes de vitesse.

Le calcul est dit convergent lorsque le résidu correspond à une valeur indiquée égale à
10-6 pour l’équation énergétique et 10-3 pour l’équation de continuité et celles des
vitesses.

Si les résidus calculés atteignent cette valeur prédéfinie, le calcul s’arrête

automatiquement puisqu’on on a atteint le critère de convergence désiré, sinon le calcul
est déclenché jusqu’à aboutir à une valeur prédéfinie qui limite le nombre d’itérations à
effectuer.

Dans notre projet le calcul s’est arrêté avant d’atteindre le nombre d’itérations
maximal désigné pour tous les cas étudiés.

****

3.3. Simulation numérique :

Dans cette partie nous allons étudier et interpréter les résultats obtenus, on va voir deux
cas différents, on va fixer une vitesse constante pour un fluide et la varier pour l’autre
fluide dans chaque cas en notant les variations de la température et de la pression dans
chaque cas afin de déterminer le débit pour lequel le système atteint une efficacité
maximale.

3.3.1. Observations et interprétations :

Dans un premier cas, nous avons étudié un transfert thermique pour un écoulement
laminaire afin d’interpréter son effet sur l’efficacité du système.

Pour établir ce régime, on a d’abord défini les conditions aux limites et on a créé par
suite deux fluides imaginaires qu’on a nommé « fluide virtuel 1 » et « fluide virtuel 2 »
afin de simplifier les calculs, tout en vérifiant les profils de vitesse et de pression pour
s’assurer que les conditions données sont convenables pour ce régime.

En effet, un régime d’écoulement dépend essentiellement du nombre de Reynolds Re,

qui s’exprime selon la relation suivante :
 v∗D
ℜ=
ϑ

v : vitesse d’écoulement en (m/s)

D : diamètre de la conduite dans laquelle s’établit l’écoulement en (m)

ϑ : viscosité cinématique du fluide en (m²/s)

Pour la résolution de notre problème on a choisi de modifier la valeur de la viscosité

Cinématique puisqu’une valeur très faible de vitesse est défavorable pour le
fonctionnement de notre système.

Pour le diamètre, on a choisi une valeur de l’ordre de 10 mm, donc on ne peut pas
diminuer cette valeur, pour cela, on a choisi de modifier la valeur de la viscosité
cinématique.

Pour un régime laminaire, le nombre de Reynolds est dans la plupart des cas inférieur à
une valeur Re=2000, on a choisi une valeur de l’ordre de 100 pour éviter d’atteindre un
régime transitoire, ceci nécessite une diminution de la valeur de la viscosité cinématique
pour obtenir le résultat voulu.

Donc on a modifié la valeur de la masse volumique et la viscosité dynamique du

fluide.

On a choisi ρ=1 kg/m3 et μ=0.01 Pa.s pour le premier fluide et ρ=10 kg/m3 et μ=0.05
Pa.s pour le deuxième (les valeurs ont étaient choisis par hasard).

Par suite on a fixé une vitesse v0=1 m/s pour le fluide extérieur « fluide virtuel 2 » et
on a pris une valeur v1= 0.1 m/s pour le fluide intérieur comme première valeur, une
température Tin=340 °K et Text=290 °K.

Après avoir saisi les différents paramètres et initialiser les conditions aux limites, on a
exécut le calcul.

************

Pour le calcul du premier cas, la solution a convergé après 60 itérations.

Les figures suivantes présentent les différents résultats observés : la variation de la

température pour tout le système, la variation de la pression au niveau des fluides
intérieurs et extérieurs et les profils de vitesse.
 Figure 16 : Variation de la température pour le système étudié

Cette première figure montre la variation de la température dans notre système, on

remarque que la température au niveau du fluide chaud tend à diminuer alors que celle
au niveau du fluide froid tend à augmenter, ceci est simplement expliqué par un flux de
chaleur dirigé du premier fluide vers le deuxième, donc le premier fluide réchauffe le
deuxième alors que le deuxième se refroidit par perte de chaleur.

Variation de la pression :
 Figure 17 : Variation de la pression pour le fluide intérieur

Figure 18 : Variation de la pression pour le fluide extérieur

Ces deux courbes montrent l’évolution de la pression pour les fluides intérieur et
extérieur, dans notre cas l’écoulement est unidimensionnel (suivant z), donc la variation
de la pression doit vérifier l’équation de conservation de quantité de mouvement qui est
simplifiée par la relation suivante :
 −1 δP δ2 w
+ϑ =0
ρ δz δ x2

Avec :

p : La pression en (Pa)

w : La composante de vitesse selon z en (m/s) Pour cela cette variation est linéaire.

Dans ce cas, l’écoulement varie dans le sens inverse de la composante z, donc

l’abscisse z=2.5 m représente l’entrée de l’écoulement et z=0 la sortie, pour cela la
courbe est croissante.

Profil de vitesse :

La courbe suivante traduit le profil de vitesse avant l’établissement du régime

laminaire.

En effet la longueur d’établissement d’un régime laminaire est exprimée par la relation
suivante

¿ ≈ 0.06 ℜ
D

Pour un écoulement en régime laminaire ou écoulement de Poiseuille, la variation de

vitesse d’un fluide visqueux dans une conduite est exprimée selon cette relation :
max

Figure 19 : Profil de vitesse avant l'établissement du régime laminaire

 Puisqu’on n’a pas encore établi ce régime, on remarque que le profil de vitesse n’est
pas
parabolique, qui est le cas d’un écoulement de Poiseuille.

Figure 20 : Profil de vitesse en régime laminaire

Après l’établissement du régime laminaire, le profil parabolique est vérifié, les
mesures sont prises pour une longueur L=1m.

3.3.2. Etude de l’efficacité

Maintenant, on va s’intéresser à la partie la plus importante de notre projet, on va

étudier l’efficacité de notre système, un échangeur thermique comprenant deux fluides
circulants dans un même sens puis dans le sens contraire

On a défini l’efficacité de notre système par la relation suivante (D’après M Ghazi

Bellakhal) :

***********

Cette relation est le rapport entre la puissance calorifique produite le fluide extérieur et
la puissance dissipée produite par la variation de pression.

Notre objectif est d’optimiser l’efficacité de notre système, c’est-à-dire maximiser ce

rapport, donc obtenir une variation importante de température en réduisant les pertes de
charges.
 Etude de l’efficacité du système (premier scénario) :

Dans ce cas, on a varié la vitesse du fluide intérieur tout en gardant la vitesse de l’autre
fluide constante sachant que ces deux fluides s’écoulent dans un même sens (co-
courant) afin de conclure à partir d’une courbe de variation de l’efficacité en fonction de
la vitesse imposée, à quel vitesse correspond la valeur maximale de l’efficacité du
système.

Les tableaux suivants indiquent les différentes valeurs calculées pour l’obtention des
résultats :

Tableau 1 : Variation de la temperature et la pression et calcul de l’efficacité (scénario

1)
V froid (m/s) V chaud (m/s) T fe (K) T ce (K) Pfe (Pa) Pce (Pa)

1 0,1 290 340 9228.2 1550.1

1 0,2 290 340 9228.4 1546.5

1 0,4 290 340 9229.8 6207.5

1 0,6 290 340 9229.8 9289.5

1 0,8 290 340 9229.8 12411

1 1 290 340 9229.8 15582

1 1,2 290 340 9229.8 16425

 Qfroid(kg/s) Qchaud(kg/s) T fs (K) T cs (K) Efficacité

0,000392493 0,000199772 308.92 306.98 0.96737

0,000392507 0,000399545 308.91 306.81 1.7939

0,000392507 0,00079909 317.8 329.03 2.6049

0,000392507 0,00119863 320.22 331.34 2.4706

0,000392507 0,00159818 321.95 332.75 2.1908

0,000392507 0,00199772 323.23 333.7 1.9227

0,000392507 0,00239727 324.24 334.45 1.7118

La courbe ci-dessus montre la variation de l’efficacité en fonction de la vitesse du

fluide intérieur.

ɛ=f(Q)
3

2.5

2
efficacité

1.5

0.5

0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
vitesse (𝑚 𝑠^−1)
 On remarque d’après cette courbe que la solution optimale correspond à une valeur de
vitesse v1=0.4 m/s. Ceci explique que le système étudié atteint une efficacité maximale
pour cette valeur.

La variation de la température dans ce cas est ∆T=27.690 °K. Cette valeur est
importante mais n’est pas la plus élevée. Pour la variation de la pression, on a obtenu
∆P= 6164.310 Pa pour le fluide intérieur qui n’est pas la valeur minimale, mais on a
obtenu le rapport le plus élevé.

Donc pour ce cas, pour une vitesse d’écoulement du fluide extérieur Vext=1 m/s, on
peut choisir une vitesse Vint=0.4 m/s pour le fluide intérieur comme vitesse
d’écoulement pour le fonctionnement le plus efficace du système en vérifiant les
paramètres caractéristiques de chaque fluide (la viscosité, la masse volumique et la
capacité calorifique).

Etude de l’efficacité du système (deuxième scénario) :

Pour le deuxième cas, on a aussi fixé une vitesse au niveau du fluide extérieur et on a
varié la vitesse du fluide intérieur (le fluide chaud) mais tout en inversant son sens
d’écoulement (contre-courant) Les tableaux suivants indiquent les résultats de calcul
pour le deuxième cas :

Tableau 2 : Tableau 1 : Variation de la température et la pression et calcul de l’efficacité

(scénario 2)
V froid V chaud T fe (K) T ce (K) Pfe (Pa) Pce (Pa)
(m/s) (m/s)
1 0,1 290 340

1 0,2 290 340 9229,75 3108,4

1 0,4 290 340 9229,75 6207,19

1 0,6 290 340 9229,73 9093,23

 1 0,8 290 340 9229,75 12458,4

1 1 290 340 9229,75 15499,8

1 1,2 290 340 9229,75 18631

Qfroid Qchaud T sortie froid T sortie chaud Efficacité

(kg/s) (kg/s) (K) (K)

0,00039251 0,00039955 318,837 2,38362

0,00039251 0,00079909 323,989 3,18482
0,00039251 0,00119863 332,964 334 3,569

0,00039251 0,00159818 328,149 2,60741

0,00039251 0,00199772 329,202 2,27882

0,00039251 0,00239727 329,949 1,99612

ꜫ=f(v)
4

3.5

2.5
efficacité

1.5

0.5

0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
vitesse(ms^-1
 Cette courbe montre la variation de l’efficacité en fonction de la vitesse du fluide
intérieur (chaud). De même dans ce cas, on a obtenu le même résultat que le précèdent,
une valeur de vitesse v2=0.4 m/s pour le système le plus efficace.

De même dans ce cas on a trouvé que la vitesse optimale ne correspond ni à la

différence de température la plus importante ni à la différence de pression la plus faible
comme le 1 er scénario mais pour les mêmes vitesses de fluide le cas du contre-courant
nous donne une meilleure efficacité.
 4. Conclusion :
Dans ce dernier chapitre, on a étudié le comportement d’un échangeur thermique en
décrivant les différentes étapes de réalisation de ce système et en visualisant les résultats
établis par CFDPost.

Dans une première partie, nous avons donné les conditions nécessaires pour obtenir un
résultat en régime laminaire, et on a affirmé ces conditions par observation du profil de
vitesse dans chaque cas.

Dans une deuxième partie, nous avons déterminé l’efficacité de notre système en
calculant les variations de la température et de la pression, en variant la vitesse de
l’écoulement pour établir une courbe de variation de l’efficacité en fonction de vitesse
pour chaque atype d’écoulements.

On a conclu finalement que la variation de la température dépend essentiellement de la

vitesse d’écoulement, plus la vitesse est faible plus la variation est importante pour le
fluide considéré, mais ceci ne traduit pas généralement que l’efficacité du système est
maximale pour la variation la plus importante de la température, il faut aussi prendre en
considération les dimensions utilisés pour réaliser la géométrie, qui affectent aussi
l’efficacité.

On doit donc vérifier tous ces conditions, tout en tenant compte des autres paramètres
qui caractérisent chaque fluide (la masse volumique, la capacité calorifique, la
viscosité…) qui sont aussi des facteurs importants pour l’évolution du comportement de
l’échangeur thermique
 Conclusion générale
La réalisation de ce projet nous a permis d’approfondir nos connaissances dans le
domaine de la mécanique de fluide et de comprendre les différents phénomènes étudiés
en maîtrisant le travail sur le logiciel de simulation CFD ainsi que les notions théoriques
sur lesquelles est basée notre application.

D’abord, nous avons appris à modéliser un écoulement en régime laminaire en

choisissant les conditions nécessaires pour aboutir à ce résultat.

Par suite, nous avons exploité ces résultats pour étudier le transfert thermique dans
notre système en variant la vitesse de l’écoulement dans chaque cas pour déterminer son
influence sur la variation de la température et de pression.

Les résultats nous ont montré que plus la vitesse est faible plus cette variation est
importante.

Cependant, pour l’efficacité du système, on a conclu que les caractéristiques

intrinsèques de chaque fluide affectent ce comportement de façon plus importante, ceci
a été traduit en deuxième cas lorsqu’on a obtenu l’efficacité la plus faible pour le débit
le plus faible. Finalement, après avoir testé les résultats pour deux fluides circulants
dans un même sens, on peut aussi tester un autre cas permettant d’évaluer le
comportement de deux fluides circulants dans deux sens opposés et comparé son
efficacité par rapport au mécanisme actuel.
 Références bibliographiques

ADIL, R. (2008). transferts de chaleur et de masse. université de Caen.

André, B., Alain, G., Charles, G., Jacques, H., Christophe, M., Pierre, M., & Roland, V. (1994).
Description des échangeurs. Université Joseph Fourier, Institut Universitaire de
Technologie, Département Génie thermique et Énergie (Grenoble).

Condenseurs / Echangeurs de chaleur . (2017). Récupéré sur buchiglas:

https://www.buchiglas.com/fr/produits/reacteurs-en-verre-unite-pilote-et-
production/accessoires-composants-buchiflex/condenseurs-echangeurs-de-
chaleur.html

DECLAYE , S., Gendebien, S., & LEMORT, V. (2017). les echangeurs de chaleur. Récupéré sur doc
player: https://docplayer.fr/24972780-P-31-les-echangeurs-de-chaleur-declaye-
sebastien-gendebien-samuel-lemort-vincent.html

ECHANGEURS DE CHALEUR TUBULAIRES EN GRAPHITE. (s.d.). Récupéré sur mersen:

https://www.mersen.com/fr/produits/equipement-anticorrosion/echangeurs-de-
chaleur/echangeurs-de-chaleur-tubulaires-en-graphite

Echangeurs en graphite. (2013). Récupéré sur Process guide des technologies de procédes:
http://processs.free.fr/Pages/VersionMobile.php?page=4528#Voir_aussi

INTRODUCTION TO A FINITE ELEMENT. (s.d.). 5.

Juan Gabriel, C., Arthur, E. B., Avram, B.-C., Peter, R., & Satyandra, K. G. (2012, mars 09). Heat
Transfer Engineering. Polymer Heat Exchangers—History, Opportunities, and
Challenges.

learn ANSYS. (s.d.). Récupéré sur FULLLIBRARY:

http://fulllibrary.com/Mechanical/ANSYS/Learn-Ansys.html

P.Tobaly. (2002). les echangeurs de chaleur.

Rouland, E. (2010/2011). les echangeurs thermiques. Uniersite de Rouen, http://gsi-

energie.univ-rouen.fr/IMG/pdf/cours-iup-me-echangeurthermique-2.pdf. Récupéré
sur http://gsi-energie.univ-rouen.fr/IMG/pdf/cours-iup-me-echangeurthermique-2.pdf
WEBER, C., & MANIFICAT, A. (2015, octobre 29). LES ÉCHANGEURS DE CHALEUR. Récupéré sur
encyclopedie de l'énergie: https://www.encyclopedie-energie.org/les-echangeurs-de-
chaleur/

Zoé, M., Nadia, C., Patrice, C., Philippe, B., Philippe, M., Jean Antoine, G., & Claude, R. (2014,
janvier 10). Échangeurs de chaleur- Description. teechnnique de l'ingénieur.