République Tunisienne
Ministère de l’Agriculture, des E.S.I.M Ministère de l’Enseignement Supérieur
Ressources Hydrauliques et de la Pêche et de la Recherche Scientifique
--*-- --*--
Institution de la Recherche et de Université de Jendouba l'Enseignement Supérieur Agricoles
Ecole Supérieur d’Ingénieurs de Medjez El
Beb
Département du Génie Mécanique et Agro-Industriel
TP : Transfert Thermique
Elaboré par :
Alaimi Meriam
Encadré par :
[Link] Wajdi
1ére année GMAI
Année Universitaire : 2023 /2024
Enoncé de l’Exercice : on désir refroidir 30 tonne/heure
d’aniline
De 17°c à 32°c dans un changeur turbulaire ,le fluide refreinant est de
l’eau de mer à 18°c disponible à raison de 40 tonne/heure
On donne les propriétés de deux fluides :
Aniline:
Cp=2100j/kg*k
ρ = 980 g/L
λ=0.17 w/mk
ɳ=0.7 centipoise
Eau :
Cp=4180j/kg*k
ρ = 1000 g/L
λ=0.6 w/mk
ɳ=0.8 centipoise
1) Calculer le flux de chaleur Q ?
2) calculer le moyenne logarithmique △TLM ?
3) calculer la surface d’échange AO ?
4) calculer le nombre de tube par passage de( NTPP) ?
5) faire une correction de vitesse ?
6) )le nombre de passage dans le tube n
7) le nombre de tube totale Ntt
8)faire une correction sur le surface d’échange et le coefficient
d’échange
9)calculer le NU à l’intérieur pour une viscosité apparente
initiale =1
10)calculer de diamètre dc de la calendre
11)calculer la vitesse longitudinale et transversale moyen de
la calendre
12)calculer le NU dans la calandre pour une viscosité
apparente=1
13)calculer le coefficient de transfert à l’extérieur et à
l’intérieur HI ET HE
14)calculer le temp (Te moy )et (Ttmoy),Tfmoy
15)calculer la viscosité à la temp moyenne de tube
16)faire une correction sur le nombre NU à l’intérieur et le
NU extérieur
17)faire une correction sur hi he
18)calculer le coefficient de transfert propre Ip
19)calculer Rs
20) vérifier si l’appareil et sous dimensionne ou non
21) calculer le perte de charge et cote disque ,cote
calendre à travers le fuseau tubulaire
Travail demandé : un code Python qui donne les résultats de cet exercice.
Le code :
import math as mt
# Données de cette problèmes
m_aniline=30 # Tonnes/heures
Tfe=57 # °C
Tfs=32 # °C
Tce=18 # °C
m_eau=40 # Tonnes/heurs
cp_aniline=2100 #J/Kg.K
P_aniline=980 # g/l
lambda_aniline=0.17 # W/m.K
u_aniline=0.7 # centipoise
cp_eau=4180 #J/Kg.K
P_eau=1000 # g/l
lambda_eau=0.6 # W/m.K
u_eau=0.8 # centipoise
L=6 #m
d=0.05 #m
# Question 1
Q=(m_aniline*1000*cp_aniline)*(Tfe-Tfs)/3600
print('la reponse de la question 1 (le flux de chaleur) est Q=',Q,'Watt' )
# Question 2
Tcs=Tce-(Q/(m_eau*1000*cp_eau/3600))
delta_T1=Tfs-Tce
delta_T2=Tfe-Tcs
moy_log=(delta_T2-delta_T1)/[Link]((delta_T2/delta_T1))
print('la reponse de la question 2 (moyenne logarithmique est
delta(LTM)=',moy_log)
# Question 3
S=Q/((1000*u_aniline*moy_log))
print('la reponse de la question 3 (surface d','échange) est S=',S,'m²')
# Question 4
NTTP=S/([Link]*L*d)
print('la reponse de la question 4 (nbrs des tubes par passage) est
NTTP=',NTTP)
#Question5
Facteur_correction_V=1.2
vitesse_corrigee_aniline = Facteur_correction_V * (m_aniline / ([Link] *
(d/2)**2 * P_aniline))
vitesse_corrigee_eau = Facteur_correction_V * (m_eau / ([Link] * (d/2)**2 *
P_eau))
print('la reponse de la question 5 (correction de la vitesse) est: correction de
la vitesse aniline=',vitesse_corrigee_aniline,'m/s correction de la vitesse
eau=',vitesse_corrigee_eau)
# Question 6
n=S/(NTTP*[Link]*(d**2))
print('la réponse de la question 6 (nbrs de passages par tube) est n=',n)
# Question 7
NTT=n*NTTP
print('la réponse de la question 7 (nbrs des tubes total) est: NTT=',NTT)
# Question 8
facteur_correction_surface = 1.1
facteur_correction_coefficient = 0.9
S_corrigée=S*facteur_correction_surface
print('la réponse de la question 8 (correction de la surface) est: S=',S)
# Question 9
T_moy_aniline=(Tfe+Tfs)/2
Re=P_aniline*vitesse_corrigee_aniline*d/(0.7*1000)
Pr=cp_aniline*u_aniline/lambda_aniline
print('la reponse de la question 9 est Re=',Re,'Pr=',Pr)
# Question 10
if Re > 2300:
Nu = 0.023 * Re**0.8 * Pr**0.4
else:
Nu = 3.66
print ('la reponse de la question 10(nbr de NUsselt) est Nu=',Nu)
# Question 11
Dc=0.1#m
print('la reponse de question 11 (diametre de la calandre) est Dc=',Dc)
# Question 12
aire_section_transversale = 0.1
vitesse_longitudinale_moyenne = m_eau / aire_section_transversale
print('la reponse de la question 12 (vitesse longitudinale)est
Vmlong=',vitesse_longitudinale_moyenne,'m/s')
# Question 13
L = 5 # Longueur du tube en mètre
Nu_intérieur = 150 # Nombre de Nusselt à l'intérieur (supposé)
Nu_extérieur = 80 # Nombre de Nusselt à l'extérieur (supposé)
lambda_tube = 0.6 # Conductivité thermique du matériau du tube en W/mK
HI = Nu_intérieur * lambda_tube /d
HE= Nu_extérieur*lambda_tube/d
print('la reponse de la question 13 (coeff de transfert inter) HI=',HI,'W/m²K et
HE=',HE,'W/m²K')
# Question 14
Te_moyenne = (Tfe+Tfs) / 2
Tt_moyenne = (Tce+Tcs) / 2
Tf_moyenne = (Te_moyenne+Tt_moyenne) / 2
print('la reponse de la question 14 Temperature moyenne de surface
dechange=',Te_moyenne,'°C temperature moyenne du
fluide=',Tt_moyenne,'°C la temperature moyenne de film=',Tf_moyenne,'°C')
# Question 15
Tm_aniline = (Tfe+Tfs) / 2
viscosite_aniline = 0.001458*((Tm_aniline**1.5)/(Tm_aniline+110.4))
Tm_eau =(Tce+Tcs)/2
viscosite_eau=0.001458*(Tm_eau**1.5)/(Tm_eau+110.4)
print('la reponse de la question 15 (viscosité)
V_aniline=',viscosite_aniline,'Pa.s')
# Question 16
Tm_eau =(Tce+Tcs)/2
viscosite_eau=0.001458*(Tm_eau**1.5)/(Tm_eau+110.4)
print('la reponse de la question 16 (viscosité) V_eau=',viscosite_eau,'Pa.s')
# Question 18
Ip = 1 / ((1 / NTT) + (1 / (NTT * Nu * lambda_eau * Dc)))
print("la reponse de question 18 (Coefficient de transfert propre)
(Ip)=",Ip,"W/m^2*K")
# Question 19
Rs = 1 / (Ip * [Link] * Dc * L)
print("la reponse de la question 19 Facteur de résistance (Rs) :", Rs)
# Question 20
seuil_sous_dimensionnement = 0.1
if Rs > seuil_sous_dimensionnement:
print("la reponse de question 20: L'échangeur est sous-dimensionné.")
else:
print("la reponse de question 20 L'échangeur n'est pas sous-
dimensionné.")
Exécution du code :
la reponse de la question 1 (le flux de chaleur) est Q= 437500.0 Watt
la reponse de la question 2 (moyenne logarithmique est delta(LTM)=
27.738874760103034
la reponse de la question 3 (surface d échange) est S=
22.531555638260453 m²
la reponse de la question 4 (nbrs des tubes par passage) est NTTP=
23.90672303586144
la reponse de la question 5 (correction de la vitesse) est: correction de
la vitesse aniline= 18.70882596345545 m/s correction de la vitesse
eau= 24.44619925891512
la réponse de la question 6 (nbrs de passages par tube) est n=
119.99999999999999
la réponse de la question 7 (nbrs des tubes total) est: NTT=
2868.8067643033723
la réponse de la question 8 (correction de la surface) est: S=
22.531555638260453
la reponse de la question 9 est Re= 1.3096178174418815 Pr=
8647.05882352941
la reponse de la question 10(nbr de NUsselt) est Nu= 3.66
la reponse de question 11 (diametre de la calandre) est Dc= 0.1
la reponse de la question 12 (vitesse longitudinale)est Vmlong= 400.0
m/s
la reponse de la question 13 (coeff de transfert inter) HI= 1800.0
W/m²K et HE= 960.0 W/m²K
la reponse de la question 14 Temperature moyenne de surface
dechange= 44.5 °C temperature moyenne du fluide=
13.29007177033493 °C la temperature moyenne de film=
28.895035885167466 °C
la reponse de la question 15 (viscosité) V_aniline=
0.002794126875870177 Pa.s
la reponse de la question 16 (viscosité) V_eau=
0.0005711019879012076 Pa.s
la reponse de question 18 (Coefficient de transfert propre) (Ip)=
516.5545797318962 W/m^2*K
la reponse de la question 19 Facteur de résistance (Rs) :
0.0012324346687585302
la reponse de question 20 L'échangeur n'est pas sous-dimensionné.
