Sujet TP Convection
Sujet TP Convection
Sujet TP Convection
2 Présentation du problème
Un échangeur de chaleur permet de transférer de l’énergie sous forme de chaleur d’un fluide
vers un autre. Ce procédé est régulièrement utilisé dans l’industrie pour le refroidissement ou le
réchauffement de liquides ou de gaz. On peut citer comme exemple le refroidissement de l’eau cir-
culant dans les radiateurs automobiles ou dans les générateurs de vapeur des centrales nucléaires
entre les circuits primaires et secondaires.
On distingue plusieurs types de géométrie pour les échangeurs classés en fonction de leur construc-
tion, de la nature et du sens de l’écoulement des fluides. On s’intéresse ici à l’un des modèles
le plus simple d’échangeur thermique, l’échangeur coaxial ou tubulaire, constitué de deux tubes
cylindriques coaxiaux.
Pour ce type d’échangeur, deux configurations existent. Soit les deux fluides s’écoulent dans le
même sens et l’échangeur thermique est dit à co-courant, soit ils s’écoulent dans des sens opposés
et l’échangeur est dit à contre-courant. Ces deux modes de fonctionnement conduisent à des per-
formances différentes. Le mode d’utilisation le plus fréquent est celui dit à contre-courant, car il a
un rendement plus élevé.
Dans le cadre de l’UE Activité expérimentale en énergétique, vous avez effectué des travaux pra-
tiques sur un banc expérimental permettant d’étudier un échangeur de chaleur coaxial dans lequel
circulent eau chaude et eau froide pour des configurations à co et contre-courant. Vous avez eu à
déterminer l’efficacité de l’échangeur thermique dans ses deux modes de fonctionnement en cal-
culant les flux de chaleur transférés des côtés chaud et froid et le coefficient de transfert global
rapporté à la surface extérieure du tube. Avec les températures relevées, vous avez comparé le profil
de température obtenu expérimentalement pour les fluides chaud et froid le long de l’échangeur au
profil théorique.
L’objectif de ce travail sera de modéliser et de simuler le banc expérimental que vous avez étudié
pour des configurations à co et contre-courant à l’aide du logiciel Ansys/Fluent afin de confronter
vos résultats numériques aux résultats expérimentaux et théoriques. Lors de ces deux séances de
travaux pratiques, vous allez dans un premier temps réaliser une géométrie représentant l’échan-
geur de chaleur que vous avez étudié expérimentalement. Vous construirez le maillage associé à
cette géométrie puis vous simulerez l’écoulement des fluides chaud et froid dans cet échangeur avec
le solveur Fluent. Vous pourrez ainsi déterminer à partir de la simulation numérique les flux de
chaleur transférés des côtés chaud et froid et les températures en sortie des fluides chaud et froid
1
Université Gustave Eiffel, Master 1 MFT 2021-2022
dans les configuration à co et contre-courant. Ces températures seront comparées aux températures
obtenues expérimentalement et théoriquement.
3 Etude théorique
On rappelle que le flux Φ cédé par le fluide chaud ou gagné par le fluide froid peut s’écrire sous
la forme suivante :
Φ = ṁCp ∆T (1)
où ṁ est le débit massique, Cp la capacité thermique et ∆T un écart de température. On montre
que le flux de chaleur (voir [1]) peut s’écrire aussi en fonction du coefficient de transfert global K
entre l’entrée et la sortie du tube considéré rapporté à la surface extérieure :
et à co-courant :
(Tc,e − Tf,e ) − (Tc,s − Tf,s )
∆Tml = (4)
ln[(Tc,e − Tf,e ) − (Tc,s − Tf,s )]
— En utilisant les relations (1), (2), (3) et (4), vous établirez l’expression des températures
d’eau chaude et froide (Tc,s et Tf,s ) en sortie de l’échangeur dans les configurations à co et
contre-courant.
— En utilisant la méthode de l’analogie électrique, calculez la résistance totale thermique
équivalente R au niveau de la paroi d’acier inoxydable de l’échangeur thermique conductivité
k, de coefficient d’échange interne à la paroi, hi , et de coefficient externe à la paroi, he . Vous
en déduirez l’expression du coefficient de transfert global Ke rapporté à la surface externe
Se sachant que K = 1/RS.
En l’absence d’informations relatives aux deux coefficients d’échange hi et he , on ne peut évaluer le
coefficient de transfert K. Vous pourrez utiliser la corrélation de Gnielinski [2] pour des écoulements
turbulents dans un cylindre pour déterminer hi . Elle s’écrit :
Avec
hi Di
N uDi = (6)
k
Cf représente le coefficient de frottement. Pour écoulement turbulent "lisse", il s’écrit :
2
Université Gustave Eiffel, Master 1 MFT 2021-2022
0, 5 ≤ P r ≤ 106
2300 ≤ Re ≤ 5.106
ū D
Dans l’espace annulaire le nombre de Reynolds ReDh = fν h s’exprime en fonction du diamètre
hydraulique Dh = 4S
P
où P est le périmètre de la section mouillée S et ūf est la vitesse débitante
du fluide froid.
La corrélation (5) dans un cylindre reste applicable dans un espace annulaire à condition d’uti-
D2 −D2
liser un diamètre équivalent tel que : Deq = eDi i . Il faut ainsi définir un nombre de Reynolds
ū D
associé à ce diamètre équivalent ReDeq = f ν eq . Le domaine de validité de la corrélation pour
l’espace annulaire est le même que celui de la conduite interne.
Eventuellement, vous pourrez comparer les coefficients d’échange obtenus avec la corrélation (5)
et avec les corrélations de Dittus Boelter ou de Sieder et Tate (voir [1]).
4 Géométrie du problème
Le problème est représenté sur les figures 1a et 1b. L’étude porte sur un échangeur thermique
formé de tubes coaxiaux en forme de U pour des raisons de compacité. On rappelle que l’eau
chaude circule dans le tube interne en acier inoxydable alors que l’eau froide circule dans l’espace
annulaire entre le tube métallique et la paroi en acrylique afin de limiter les pertes de chaleur du
fluide externe vers le milieu ambiant. On rappelle les caractéristiques géométriques de l’échangeur :
3
Université Gustave Eiffel, Master 1 MFT 2021-2022
Vous allez considérer le domaine simplifié 2D présenté sur la figure 1b pour simuler le problème
sous Ansys/FLuent. Vous prendrez obligatoirement l’axe (O~x) comme axe de symétrie.
Tracer l’esquisse
Pour ce problème, choisissez le plan XY (Dans Arborescence, sélectionnez plan XY), et le re-
père apparait. Dans la barre d’outils, cliquez sur l’outil voir face, pour une vue perpendiculaire
au plan.
Pour construire la géométrie, cliquez sur Esquisse et sélectionnez Dessiner. A l’aide des outils
proposés, réalisez la géométrie présentée sur la figure 1b en respectant les dimensions données.
4
Université Gustave Eiffel, Master 1 MFT 2021-2022
5
Université Gustave Eiffel, Master 1 MFT 2021-2022
pour la résolution. Il est donc nécessaire de comprendre l’écoulement que l’on résout ainsi que la
signification des options que l’on coche sur Fluent.
Les écoulements dans les tubes étant turbulents, incompressibles et stationnaires en moyenne
temporelle, vous résoudrez les équations de Navier-Stokes incompressibles pour calculer les écou-
lement moyens et vous prendrez un modèle de turbulence de type RANS (Reynolds Average
Navier-Stokes), appelé k-ω-SST, pour modéliser l’influence des fluctuations turbulentes sur ces
écoulements. L’équation de conservation de la masse sera résolue ainsi que l’équation de convec-
tion/diffusion de l’énergie dans les deux tubes. Dans les deux parois solides, en acier et en acrylique,
l’équation de conduction de la chaleur sera résolue. La résolution de ces équations nécessitera d’im-
poser des conditions aux limites en adéquation avec les conditions expérimentales. Ces conditions
aux limites sont données dans la figure 4.
Remarque : Au niveau de la paroi extérieure en contact avec l’air ambiant, un échange de chaleur
par convection et rayonnement peut être imposé si la paroi extérieure du tube externe n’est pas
complètement isolé. Dans ce cas, il faut connaître le coefficient d’échange convectif par convection
naturelle hconv , agissant sur la paroi cylindrique externe mais également le coefficient radiatif hr
sur cette surface.
6
Université Gustave Eiffel, Master 1 MFT 2021-2022
Setup :
— General :
• Solver → Cochez Pressure-Based, Absolute, Steady, axisymmetric
— Models :
• Energy → Activez la résolution de l’équation d’énergie
• Viscous → Edit, choisissez k-ω model et sélectionnez SST
— Materials → Create/Edit, cherchez water-liquid dans l’étiquette Fluent Database,
cliquez sur copy pour qu’un nouveau fluide apparaissent et modifiez les caractéristiques
de l’eau. Vous penserez à différencier eau chaude et froide en fonction des caractéristiques
de l’eau ρ(T ), µ(T ) et Cp (T ) en rentrant par exemple un profil polynomial de ces deux
variables plutôt que constant dans Fluent. Vous ajouterez également les solides acier
inoxydable et acrylique en conservant les paramètres par défaut. Si le matériau n’existe
pas, vous avez la possibilité de le créer dans Fluent en y rentrant les paramètres du matériau
souhaité.
— Boundary conditions
• inlet eau froide changez le type à mass flow inlet → Edit, entrez dans mass flow
rate le débit massique d’eau froide, sélectionnez Direction Vector pour Direction
specification method et entrez 1 pour Axial component of flow direction. Dans
la partie turbulence, laissez le taux de turublence pr défaut et rentrez les diamètre
hydraulique de la conduite froide. Dans l’onlget Thermal rentrez la température de
l’eau froide en entrée.
• inlet eau chaude effectuez la même opération que pour l’eau froide en adaptant les
propriétés
• Outlet eau froide changez le type à Pressure Outlet
• Outlet eau chaude changez le type à Pressure Outlet
• Wall acier gauche changer le type à Wall
• Wall acier droite changer le type à Wall
• Wall acry gauche changer le type à Wall
• Wall acry droite changer le type à Wall
• Wall acry inf changer le type à Wall, sélectionnez Stationary wall dans Wall mo-
7
Université Gustave Eiffel, Master 1 MFT 2021-2022
Solution :
— Solution Methods
• Pour Pressure-Velocity coupling → SIMPLE
• Pour Pressure → PRESTO ! et laissez les autres paramètres par défaut
— Monitors Décocher Check convergence
— Calculation activities → Run calculation et 100 itérations dans Number of Time
Steps pour commencer.
Cliquez ensuite sur Calculate et observez l’évolution des résidus au cours du temps. Pensez à
sauvegarder toute la configuration du modèle avec File → Write → Case dans un fichier *.cas
(ces fichiers contiennent l’ensemble des paramètres physiques et numériques de la simulation).
Effectuez ces calculs pour des débits d’eau froide Qvf différents (1L/mn, 2L/mn, 3L/mn).
8
Université Gustave Eiffel, Master 1 MFT 2021-2022
A la fin de vos calculs, pensez à sauvegarder toute la configuration du modèle ainsi que les
données de calculs avec File → Write → Case and Data dans des fichiers *.cas et *.dat (ces
fichiers contiennent l’ensemble des paramètres physiques et numériques de la simulation ainsi que
les résultats).
Le tableau 5 présente les valeurs de températures relevées expérimentalement ainsi que les
incertitudes sur ces températures pour un débit d’eau chaude Qvc fixé à 2L/mn et un débit d’eau
froide Qvf allant de 1 à 3L/mn dans les configurations à co et contre-courant.
La figure 6 présente l’emplacement des thermocouples où sont relevées les température en entrée,
9
Bibliographie
[1] G. Lauriat, Echangeurs de chaleur : Licence de génie des procédés licence de physique et
application
[2] A. Bejan, Convection heat transfer, Wiley, 1995
[3] S. M’Zali et H. Iabbadene, Rapport de projet long : simulation numérique d’un échangeur
thermique coaxial
[4] M. Tavares, Rapport de stage : étude expérimentale et numérique d’un échangeur de chaleur à
co et contre courant
10