Evaporateur à film tombant

I. But :
Le but de cette manipulation est l’étude d’un évaporateur à film tombant et de chercher les
conditions opératoires optimale de sa fonctionnement.

II. Etude théorique :

On considère un débit volumique de l’alimentation liquide Q =40 L/h.

1. Estimation de l’épaisseur du film liquide à l’intérieur des tubes de l’évaporateur :

Un bilan énergétique sur l’évaporateur donne :

Ece+Epe=Ecs+Eps (on néglige les pertes de charge le long des

tubes)

12 m V12+m g z1=12 m V22+m g z2

V12+2 g h= V22 (1)

Or V1=Qeπ4D2 et V2 =Qsπ4 (D2-(D-e)2)

Sachant que Qs=Qe=Q

(1) → (Qπ4D2)2+2 g h=(Qπ4 (D2-(D-e)2) )2 (2)

1D4+(π4)2 2 g hQ2=1(D2-(D-e)2)2

Or D2-D-e2=D-D+eD+D-e

=e 2 D-e

=2 e D 1-e2 D

e <<<D → D2-D-e2=2 e D

(2) → 1D4+(π4)2 2 g hQ2=12 e D2

4 D2(1D4+π422 g hQ2)=1e2

e=14 D2(1D4+ (π4)2 2 g h Q2)

D= 24 10-3m ; h=1.5 m ; Q= 40 L/h

AN : e=53.81 10-6 m

1. Détermination du temps de séjour moyenne:

ts=VQ Avec V=π4D2-D-e2 h=π4 2 e D h

ts=π e D h2 Q

AN : ts=π*53.8 *10-6(24 *10-3)1.52*(40 10-3)*3600=0.273 s

Evaporateur à film tombant

ts=0.273 s

2. Détermination du taux d’évaporation:

Solution sucrée 0.1 sucre

me 0.9 eau

Concentrât 0.2 sucre

Vapeur ms 0.8 eau

Un bilan massique :

➢ Bilan global : me-ms=mevaporée (1)

➢ Bilan partiel par rapport à l’eau : 0.9*me-0.8*ms=mevaporée (2)

(1)+(2) → 0.9*me-0.8*(me-mevaporée)=mevaporée

→ mevaporée=me*(0.9-0.8)(1-0.8)

AN : mevaporée=40*(0.9-0.8)(1-.8)=20 kg/h

Taux d’évaporation : ζevap = mevaporéeme=2040=50%

1. Puissance thermique nécessaire pour réaliser un taux de conversion de 50% :

On établie un bilan thermique sur l’évaporateur :

Q=mevaporée*Lv+me*Cp*(Tevap-Te)

Pour P=1 bar Tevap=99.6 °C et Lv=539.1 kcal/kg

Q= 203600*539.5*4.18*103+403600*99.6-25

Q =12529.2 W

2. Le débit nécessaire de vapeur chauffante pour fournir cette puissance :

Q=W*Lv → W=QLv

Pour P=2.5 bars ; Lv=520.7 kcal/kg

W=12529520.7*103*4.18=5.756*10-3kg/s

I. Manipulation :
Evaporateur à film tombant

1. Etalonnage du débitmètre :

On fait varier le débit de liquide pompé par la pompe en plaçant le flotteur en différentes
positions. Par la suite on mesure le débit volumique de liquide sortant du ballon tampon. On
aura le tableau des valeurs suivant :

Position Volume Temps Débit volumique Débit

du calculé (mL/s) volumique
flotteur (mL) (s) calculé (L/h)
60 285 275 355 5.34 4.9 6.5 53.3 55.2 54.3 195.54
8 3 7 2 6
50 270 330 225 5.91 7.4 5.0 45.6 44.5 44.5 161.712
1 5 8 3 5
40 230 250 249 6.4 6.7 6.8 35.9 36.9 36.2 131.052
6 4 8 9
30 210 210 200 7.5 7.1 6.7 28 29.4 29.7 104.7
2 2 9 6
20 148 180 140 7.08 9.1 9.1 20.9 19.7 19.7 72.39
2 1 4
10 130 100 90 12.1 9.2 8.6 10.7 10.7 10.3 38.22
4 9 7 1 6 8

Courbe d’étalonnage obtenue par EXCEL :

On peut donc déterminer la position correspondante à un débit volumique de 40 L/h à partir

de la courbe d’étalonnage obtenue : c’est la position 10.5.
Evaporateur à film tombant

On va travailler en circuit ouvert. On fixe le débit de l’alimentation liquide à 40L/h. La vapeur

chauffante est fournie par la chaudière avec un débit contrôlable par l’intermédiaire d’un
vanne pneumatique. On va donc mesurer les différents paramètres des fluides circulants dans
l’unité pour diverses ouvertures de cette vanne. On obtient les tableaux des valeurs suivants :

Ouverture de la vanne :100%

Q(L/h) Te(°C) Ts(°C) P(bars)
Alimentation 40 23.9 100
liquide
Vapeur 17.78 113.3 81 1.6
chauffante
Liquide 187.7 24.6 47

De
refroidissement
Condensat 9.28 100 40
Liquide non 31.92 23.9 100
évaporé

Ouverture de la vanne :80%

Q(L/h) Te(°C) Ts(°C) P(bars)
Alimentation 40 23.9 100
liquide
Vapeur 17.54 116.92 81 1.8
chauffante
Liquide 542 24.6 37

De
refroidissement
Condensat 9.38 100 34
Liquide non 30.62 23.9 100
évaporé

Ouverture de la vanne :60%

Q(L/h) Te(°C) Ts(°C) P(bars)
Alimentation 40 23.9 100
Evaporateur à film tombant

liquide
Vapeur 17.3 116.6 80.3 1.8
chauffante
Liquide 569.45 24.6 34.6

De
refroidissement
Condensat 9.44 100 32.8
Liquide non 30.56 23.9 100
évaporé

2. Détermination de l’enthalpie du condensat :

Bilan énergétique sur l’évaporateur (calcul de l’enthalpie de condensation de vapeur

chauffante) :

Hs,;Ws hF ;m

mL; HL

Hcond ;Ws mv;Hv

-E+S=Q+W

On suppose que Q=W=0

-E+S=0

Hcond *Ws+mv*Hv+mL*HL-m*HF-WsHs=0

Hcond=1Ws*(-mv*Hv-mL*HL+m*HF+WsHs)

a. Pour une ouverture de 100%:

m=40 kg /h à T=23.9 °C HF=23.96 kcal/kg

Ws=17.78 kg/h à T=113.3°C Hs=643.7 kcal/kg

mv=9.28 kg/h à T=100°C Hv=639.1 kcal/kg

Evaporateur à film tombant

mL=31.92 kg/h à T=100°C HL=100.1 kcal/kg

AN : Hcond=117.78*40*23.96+643.7*17.78-9.28*639.1-31.92*100.1

Hcond=184.32 kcal /kg=770.49 kJ

Pour cet essai le mélange obtenu à la sortie de l’évaporateur est caractérisé par :

P=Pe=1.6 bars

T=81°C

x=Hcond-HLHg-HL=184.32-81.004629.7-81.004=18.85%

b. Pour une ouverture de 80% :

m=40 kg /h à T=23.9 °C HF=23.96 kcal/kg

Ws=17.54 kg/h à T=116.92°C Hs=644.9 kcal/kg

mv=9.38 kg/h à T=100°C Hv=639.1 kcal/kg

mL=30.62 kg/h à T=100°C HL=100.1 kcal/kg

AN :

Hcond=117.54*40*23.96+644.9*17.54-9.38*639.1-30.62*100.1

Hcond=183.017 kcal /kg=765.013 kJ

Pour cet essai le mélange obtenu à la sortie de l’évaporateur est caractérisé par :

P=Pe=1.68bars

T=81°C

x=Hcond-HLHg-HL=18.83%

c. Pour une ouverture de 60% :

m=40 kg /h à T=23.9 °C HF=23.96 kcal/kg

Ws=17.3 kg/h à T=116.92°C Hs=644.9 kcal/kg

mv=9.44 kg/h à T=100°C Hv=639.1 kcal/kg

mL=30.56 kg/h à T=100°C HL=100.1 kcal/kg

AN : Hcond=117.3*40*23.96+644.9*17.3-9.44*639.1-30.56*100.1

Hcond=174.6 kcal /kg=729.835 k

Pour cet essai le mélange obtenu à la sortie de l’évaporateur est caractérisé par :

P=Pe=1.8 bars
Evaporateur à film tombant

T=80.3°C

x=Hcond-HLHg-HL=174.6-80.3012629.42-80.3012=17.17%

1. Calcul de chaleur d’évaporation :

Qevap=Ws*hs-hcond

a. Pour une ouverture de 100% :

Qevap=17.78*643.7-184.32=8167.77 kcal /h

Qevap=9.483 kw

b. Pour une ouverture de 80% :

Qevap=17.54*644.9-183.017=8101.42 kcal /h

Qevap=9.406 kw

c. Pour une ouverture de 60% :

Qevap=17.3*644.9-174.6=8136.19 kcal /h

Qevap=9.447 kw

1. Calcul de taux d’évaporation :

ζevap=mevapm=mvm

a. Pour une ouverture de 100% :

ζevap=9.2840=23.2%

b. Pour une ouverture de 80% :

ζevap=9.3840=23.45%

c. Pour une ouverture de 60% :

ζevap=9.4440=23.6%
Evaporateur à film tombant

On trace la courbe représentant la variation de taux d’évaporation en fonction de débit

massique de la vapeur chauffante :

Interprétation :

On constate que la courbe a une allure parabolique. Le taux d’évaporation diminue avec
l’augmentation du débit de la vapeur chauffante.

1. Calcul de coefficient global d’échange :

Ts Te

ts te

∆Ѳ1=Te-te ; ∆Ѳ2=Ts-ts ; F=1 (car il s’agit d’un évaporateur)

U=QA*(DTLM)

Avec : DTLM=∆Ѳ1-∆Ѳ2ln(∆Ѳ1∆Ѳ2) et A=Se*nt =1.47 m2

a. Pour une ouverture de 100% :

me=187.7 kg /h Te=Ts=113.3°C

te=23.9°C et ts=100°C

DTLM=89.4-13.3ln(89.413.3)=39.9°C et U=9.4831.47*39.9=161.68 wm2°C

b. Pour une ouverture de 80% :

me=542 kg /h Te=Ts=116.92°C

te=23.9°C et ts=100°C

DTLM=93.02-16.92ln(93.0216.92)=44.65°C et
U=9.4061.47*44.65=143.306 wm2°C
Evaporateur à film tombant

c. Pour une ouverture de 60% :

me=187.7 kg /h Te=Ts=116.6°C

te=23.9°C et ts=100°C

DTLM=92.7-16.6ln(92.716.6)=44.65°C et
U=9.4471.47*44.65=143.93 wm2°C

On trace la courbe de variation de coefficient global d’échange en fonction de débit massique

de vapeur chauffante :

Interprétation : on constate que le coefficient global d’échange pour des débits de vapeur
chauffante inferieures à 17.5 kg/h est presque constante ou varie très peu alors que pour des
valeurs supérieurs de débit il augmente avec ce dernier : plus le débit de vapeur de chauffage
est important plus le coefficient global d’échange est élevé.

1. Essai condenseur :

Qcond=26.44 L/h ; Tcond=67.9°C

Tle=23.9°C ; Tls= 60°C

Qc=mliqHls-Hle=mvLv

mv=26.443600*10-3*314=3.03*10-3 kg/s ; mliq=1953600*10-3*1000=54.17

10-3kg /s

Lv=527.8 kcal/kg=2209.4 kJ/kg

Qc=3.03*10-3*2209.4=6702 w donc Qc≈6702 w

Uexp=QcA*∆Tmoy Avec A=1.47 m2

∆Tmoy=60-23.9ln(6023.9)=39.22

Uexp=67021.47*39.22=116.25 w/ m2 °C

Ucalc=1he+eλ+1hi avec he=0.945*ρ(ρ-ρv)L μf(Tv-Tl) g hfg (Kf)3(14)

Or ρv ≪ρ he=0.945*ρ2L μf(Tv-Tl) g hfg (Kf)3(14)

A tmf =(116.92+23.9)/2=70.41°C

hfg =2333 kJ/kg

ρ=977.56 kg/m3 →he=3116.978 w /m2°C

Kf=0.6683 w /mK
Evaporateur à film tombant

μf=402.94 10-6 kg/mK

hi est donné par NUSSELT : Nu=D hiλ

Re=4mlπ di μl=4*54.17*10-3π*22.8*10-3*(10-3)=3025

Re>2100 donc le régime est turbulent et Nu=4200*1.35+0.02*tmf*(vt)0.8di0.2

vt=ml/ρA

vt=54.17*10-3/977.561.47=0.0368 10-3 m/s

AN : Nu=4200*1.35+0.02*41.95*(0.0368* 10-3)0.8(22.8)0.2=1.39

Donc hi=1.39*1.20.0228=73.158 w /m2°C

1Ucalc=13116.978+53.81* 10-61.2+173.158 =0.014

Ucalc=71.43 w /m2°C

Remarque : on constate qu’il ya une différence entre la valeur calculée et celle

trouvée expérimentalement ceci est du aux divers erreurs de calcul : le calcul
des débits ,les mesures des températures…en plus des pertes.

I. Conclusion :
Dans cette TP nous avons étudié un échangeur multitubulaire à
simple effet : évaporateur à film tombant. Nous avons étalonné le débitmètre de
l’installation et réalisé des divers essais dans des conditions opératoires
différentes affin de dégager les valeurs optimales de débit de la vapeur
chauffante qui permet de maximiser le taux d’évaporation et le coefficient
d’échange : cette valeur correspond à une ouverture de 80% de la vanne
pneumatique. En plus nous avons fixé le débit de l’eau d’alimentation à une
valeur maximale pour que l’échangeur fonctionne en tant que condenseur et
nous avons calculé son coefficient global d’échange théoriquement par les
corrélations et expérimentalement par des bilans matière et enthalpique.