I.
Introduction
Un échangeur de chaleur est un dispositif permettant de transférer de l'énergie thermique d'un fluide
vers un autre avec des températures différentes sans les mélanger. Le flux thermique y traverse la
surface d'échange qui sépare les fluides.
Il existe de nombreuses technologies d’échangeurs de chaleur, mais tous fonctionnent suivant les
mêmes processus physiques d’échanges.
Les échangeurs de chaleur sont utilisés dans de nombreuses applications industrielles telles que la
climatisation, l’automobile (radiateurs de refroidissement), l’industrie chimique.
Dans ce TP, nous allons simuler des écoulements et des transferts de chaleur dans un échangeur de
chaleur coaxial à l’aide du logiciels ANSYS Fluent et de comparer les performances des échangeurs à
contre et à co-courant pour déterminer lequel est le plus efficace, ainsi de valider les expériences
qu’on a pu faire dans le cadre d’Activités Expérimentales.
A)Types d’échangeurs :
Il existe plusieurs types des échangeurs de chaleur, chacun ayant des avantages et inconvénients.
Le choix de l’échangeur dépend de domaine de l’utilisation. Voici les types les plus courants :
- Echangeur tubulaire : ces échangeurs de chaleur ont plusieurs tubes parallèles disposés dans un
faisceau. Echangeur à spirales.
- Echangeur à plaques : ce type d'échangeur de chaleur utilise des plaques minces pour
transférer la chaleur entre les deux fluides.
- Echangeurs à coque et tubes : ce type d'échangeur de chaleur est constitué d'un tube à
l'intérieur d'une coque.
- Echangeurs à condensation : ces échangeurs de chaleur sont utilisés pour condenser un fluide
de vapeur en liquide Echangeurs à récupération d'énergie.
- Échangeur à bloc.
Figure 01 : échangeur tubulaire Figure 02 : échangeur à plaque
À co-courant : les deux fluides sont disposés parallèlement et vont dans le même sens.
À contre-courant : mais les courants vont dans des sens opposés.
À tête d'épingle : un des deux fluides fait un demi-tour dans un conduit plus large, que le deuxième
fluide traverse.
À courants croisés : les deux fluides sont positionnés perpendiculairement.
À contact direct ou à mélange : les deux fluides peuvent être mis en contact comme c'est le cas
dans les tours de refroidissemen
II. Présentation du problème :
Comme on l’a mentionné précédemment, il existe plusieurs types d’échangeurs. Pour ce TP, nous
allons s’intéresser à 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. Co-courant où les deux fluides s’écoulent
dans le même sens. Contre-courant où les deux fluides s’écoulent dans des sens opposés. Ces deux
modes de fonctionnement conduisent à des performances différentes.
Les hypothèses du problème :
Ecoulement 2D axisymétrique.
Ecoulement stationnaire ( ∂∂t =0) .
Ecoulement incompressible ( ρ=cste ) .
Ecoulement visqueux newtonien ( μ=μeau=cste ) .
Description du l’étude muni dans le TP
L'objectif de cette étude consiste à modéliser et simuler le banc expérimental que nous avons étudié pour les
configurations co-courant et contre-courant. Nous comparerons les résultats avec ceux obtenus lors de
l'expérience réalisée sur le banc expérimental dans le cadre du module de planification des expériences, ainsi
qu'avec les résultats théoriques.
Pour ce faire, nous déterminerons les températures de sortie du fluide chaud et froid, ainsi que le flux échangé
entre ces deux fluides. Ces valeurs seront ensuite comparées aux résultats obtenus expérimentalement et aux
prévisions théoriques.
A) Description de l’échangeur étudier
Cette étude porte sur un échangeur coaxial composé de deux tubes concentriques,
à travers lesquels les fluides circulent. Le tube intérieur est parcouru par le fluide froid,
tandis que le fluide chaud circule dans l'espace annulaire entre les deux tubes.
Dans la configuration en co-courant, les deux fluides circulent dans la même
direction. En revanche, dans la configuration contre-courant, le sens de circulation du
fluide froid est inversé par rapport au sens du fluide chaud
6
Figure 2 : Evolution des températures des deux fluides le long de l’échangeur, (co-
courant à gauche et contre-courant à droite).
Dans notre étude, nous avons un échangeur coaxial composé principalement de deux tubes concentriques. Le
fluide froid circule à travers le tube extérieur, délimité d'un côté par une paroi en acier en contact avec le fluide
chaud, et de l'autre côté par une paroi en acrylique. Le fluide chaud, quant à lui, circule à travers le tube
intérieur, séparé du fluide froid par une paroi en acier.
Figure 3 : représentation de la géométrie étudiée
III. Etude théorique :
Avant de procéder à des mesures expérimentales sur l’échangeur coaxial à courant parallèle (Co-courant) ou à
courant opposé (Contre-courant), il est important de comprendre les équations qui interviennent dans le calcul
des performances énergétique au sein de cet échangeur, toutefois, le développement de la théorie ne fera pas
l’objet de cette étude, on va présenter dans ce qui suit, les équations nécessaire à la compréhension de la
technologie, et des principes de fonctionnement des échangeurs de chaleur, ces équations seront expliquées
lors de la soutenance de notre projet (voir les cours de Guy Lauriat sur les échangeurs de chaleur pour plus
d’informations).
Le flux Φ cédé par le fluide chaud ou gagné par le fluide froid (échangeur sans pertes) :
Φ=ṁC p ∆ T
𝜙 = − 𝑚̇𝑐𝑐𝑝𝑐∆𝑇𝑐 = 𝑚̇𝑓𝑐𝑝𝑓∆𝑇𝑓 (1)
Où, ṁ : le débit massique,
C p: la capacité thermique,
∆ T : l’écart de température.
Le débit massique pour un fluide chaud et froid :
𝑚̇𝑐 = 𝜌𝑐𝑆𝑐𝑉𝑠 𝑚̇𝑓 = 𝜌𝑓𝑆𝑓𝑉𝑓
𝜌: Masse volumique du fluide chaud (Kg/m3).
𝑆 : Surface d’écoulement du fluide chaud dans le tube.
𝑉: vitesse du fluide chaud.
Flu le flux peut s’écrire en fonction du coefficient de transfert global K :
Φ=( KS )e ∆ T ml (2)
Où, K : Coefficient de transfert global,
∆ T ml : L’écart de température moyen logarithmique.
On définit l’écart de température moyenne logarithmique à co-courant et à contre-courant :
( T c ,e −T f ,s ) −(T c , s−T f , e )
ΔT ml , contre =
ln [ ( T c, e −T f , s ) /(T c , s−T f , e )]
( T c ,e −T f ,e )−(T c, s−T f ,s )
ΔT ml , co=
ln [ ( T c, e −T f , e ) /(T c , s−T f ,s )]
a) Expressions des températures d’eau chaude et froide en sortie de
l’échangeur :
On a :
T c ,s=
[
T c ,e ( B−1 )−B {1−exp ( KS )e
( B−1 )
ṁf C p]}T f ,e
[
B exp ( KS )e
( B−1 )
]
ṁ f C p
−1
Avec le même raisonnement, l’expression de la température de sortie du circuit chaud à co-courant.
T c ,s=
[
{B+ exp ( KS )e
( B+1 )
ṁf C p ]
}T −B {exp ( KS )e
c ,e
[
( B+1 )
ṁf C p
−1 }T f , e
]
[
(B+1)exp ( KS )e
( B+1 )
ṁf C p ]
b) Expressions des températures d’eau chaude et froide en sortie de l’échangeur :
On a :
ṁf C pf T c, e −T c, s
ṁf C pf ( T f , s−T f , e ) =ṁc C pc ( T c, e −T c, s ) → = =B
ṁc C pc T f , s−T f ,e
( T c ,e −T c, s )
→ T f , s=T f , e +
B
L’expression précédente est valable pour les deux sens de circulation.
- Cas contre-courant :
( T c, e −T f , s )−( T c , s−T f , e )
ṁf C pf ( T f , s−T f , e ) =( KS )e ∆T ml= ( KS )e
ln ( ( T c ,e −T f ,s ) / ( T c , s−T f , e ) )
→ ln
( ) =
[
( T c ,e −T f , s ) ( KS )e ( T c ,e −T c ,s )
( T c ,s−T f ,e ) ṁf C pf ( T f ,s −T f , e )
−1 =
]
( KS )e
ṁ f C pf
[ B−1 ]
→
( T c ,e −T f ,s )
( T c , s−T f , e )
=exp
[
( KS )e
ṁf C pf ] ( T −T )
[ B−1 ] avec T f , s=T f ,e + c ,e c , s
B
→ T c , s=
T c , e ( B−1 )−B 1−exp
[ [ ( KS )e
ṁf C pf ]]
[ B−1 ] T f , e
B exp
[ ( KS )e
ṁf C pf ]
[ B−1 ] −1
- Cas co-courant :
( T c, e −T f , e ) −( T c, s−T f ,s )
ṁf C pf ( T f , s−T f , e ) =( KS )e ∆T ml= ( KS )e
ln ( ( T c ,e −T f ,e ) / ( T c, s−T f ,s ) )
→ ln
( ) =
[
( T c ,e −T f , e ) ( KS )e ( T c , e −T c , s )
( T c , s−T f , s ) ṁf C pf ( T f , s−T f , e )
+1 =
]
( KS )e
ṁf C pf
[ B+1 ]
→
( T c ,e −T f ,e )
( T c, s−T f ,s )
=exp
[
( KS )e
ṁf C pf ] ( T −T )
[ B+1 ] avec T f , s=T f ,e + c ,e c , s
B
→ T c , s=
[ [
B+exp
( KS )e
ṁf C pf ]] [ [
[ B+1 ] T c ,e + B exp
( KS )e
ṁf C pf ] ]
[ B+ 1 ] −1 T f , e
( B+1 ) exp
[ ( KS )e
ṁf C pf
[ B+1 ]
]
Maintenant on cherche à déterminer le coefficient global d’échange K e qui est donné par la formule
K e =1/ R S e
Avec Se la surface d’échange extérieur et R la résistance thermique.
a) Calcul de la résistance thermique équivalente R :
on peut calculer la résistance totale équivalente R par utilisation de la méthode analogie
électrique, qui est égale à la somme de trois résistances équivalentes.
Dans l’échangeur y’a deux mode de transfert de chaleur, transfert de chaleur par conduction
au niveau de la paroi en acier séparent les deux fluides et transfert par convection échange
interne et externe.
L’expression de la résistance totale :
R=
1
+
1
2 πLr i hi 2 πLk a
r
ln e +
() 1
r i 2 πL r c he
On prend : 𝑘𝑎 = 16.3 W/𝑚. 𝐾
K=1/RS
Ou
1
K e=
1 e re 1
r
+ ln +
hi k r i he
1
K i=
1 i re r i
r
+ ln +
hi k r i r e he
E) Détermination des coefficients d’échange hi et he
Si les coefficients d'échange hi et he ne sont pas connus, il n'est pas possible de calculer
le coefficient de transfert K. Cependant, on peut utiliser la corrélation de Gneilinski pour
estimer hi dans le cas d'écoulements turbulents à l'intérieur d'un cylindre. Cette
corrélation est exprimée de la manière suivante :
(C ¿¿ f /2)(ReDi −103) Pr
N uDi = ¿
1+ 12,7 ¿¿
cette corrélation est valable pour :
0,5≤ Pr ≤ 106
2300≤ ReDi ≤5. 106
Pour un écoulement turbulent laisse, le coefficient de frottement :
−1
4
C f ≈ 0,078 ReDi
Avec le nombre de Nusselt
hD
N u=
k
Calcule de nombre de Reynolds (𝑹𝒆𝑫𝒊 et 𝑹𝒆𝑫𝒉) et le nombre de
Prandtl (Pr) :
Pour vérifier la validité de la corrélation il faut calculer le nombre de Prandtl et le nombre de
Reynolds. Le domaine de validité de la corrélation
0,5 ≤ 𝑃𝑟 ≤ 106
2300 ≤ 𝑅𝑒 ≤ 5. 106
Le nombre de Reynolds dans la conduit interne ou circule le fluide chaud
u c Di
ReDi =
ν
uc : vitesse débitante du fluide chaud.
Di : diametre de la conduite interne
𝜈 : la viscosité cinématique de l’eau.
Le nombre de Reynolds dans l’espace annulaire ou circule le fluide froid
u f Di
ReDh =
ν
uf : vitesse débitante du fluide froid
4S
Dh : diametre hydraulique, D h =
P
𝜈 : la viscosite cinematique de l’eau
Le nombre de Reynolds dans l’espace annulaire
uf Deq
ReDeq =
ν
2 2
D c −Di
Le diametre equivalent est defini par : Deq = 2
Di
Le nombre de Prandtl
ν
Pr=
α
α : la diffusivité thermique
𝜈 : la viscosité cinématique de l’eau.
Etablissement de l’expression des coefficients d’échanges ℎ𝑒 et ℎi
Puisque la corrélation est valable, donc : Pour la conduite interne :
hi D i k N uDi
N uDi = hi =
k Di
((C ¿¿ f /2)(ReDi −103 ) Pr )
hi =k ¿
Di ¿ ¿
La corrélation de Gneilinski est valable tant pour l'espace annulaire que pour la conduite interne. Pour
obtenir la corrélation en fonction du diamètre équivalent dans l'espace annulaire, il suffit de l'adapter
en conséquence. À partir de cette adaptation, on peut déduire le coefficient d'échange externe he :
he D eq k N uDi
N uDi = h e=
k D eq
((C ¿¿ f /2)( ReDeq −103 ) Pr )
h e=k ¿
D eq ¿ ¿