Calcul
Calcul
Calcul
Q
EV = ————
1000
E
D= ————
C-1
-La purge se fait uniquement quand le taux de concentration « C » est supérieur à 3.
Salinité du circuit
C = —————————
Salinité d’appoint
Habituellement,ce taux de concentration est déterminé par la mesure des chlorures, facilement
dosable qui du fait de leur solubilité sont les plus stables dans le circuit.
A = D + E + EV + fuites
a=S/V
Simulation de la circulation d’eau de refroidissement avec un appoint d’eau fraîche :
1- Introduction :
1- La vitesse maximale du vent ne dépassant pas 5 m/s,elle est mesurée à une hauteur
de 1,5 à 2 m du sol.
2- La température du bulbe humide de l’air à l’entrée ne dépassant pas ± 5 °C par
rapport au design.
La quantité de l’eau circulante ne doit pas contenir des matières dissoutes plus que 5000
ppm, et pas plus de 10 ppm de traces d’hydrocarbures.
- le niveau d’eau refroidie dans le bassin ne dépasse pas la hauteur normale.
- Le débit d’appoint est suffisant pour récupérer toutes les pertes qui sont rencontrées
durant le fonctionnement.
- La quantité d’eau circulante ne doit pas être supérieure ou inférieure à 10% par
rapport au design.
3- Calcul massique et thermique :
Données théoriques :
Te = 25 °C : température au bulbe d’air sec à l’entrée.
Ts = 30 °C : température de sortie de l’air.
D’après le diagramme psychométrique on aura :
à Te = 25 °C, on a l’humidité relative ε1 = 70%.
à Ts = 30 °C, on a l’humidité à la saturation ε2 = 100%.
Sachant que :
Pi εi 18
Xi = —————— ——
P atm – Pi εi 29
Pe = 23,756 mmHg
Ps = 31,824 mmHg
Hs = 23,706 kcal/kg
Donc : ΔH = HS - He
ΔH = 9,333 kcal/kg
Calcul de la quantité d’eau évaporée au niveau de la tour 33 E101 :
Pour la détermination de la quantité d’eau évaporée il y’a lieu d’effectuer un bilan de matière
global au niveau de la tour, tels que représenté par le schéma ci-dessous.
G ,HS ,TS,E
La tour de
Q, h1, ta Refroidissement Q, h2, ta
33 E101
G ,He ,Te
Le débit massique d’air sec est constant, tandis que le débit massique d’air humide augmente
par suite de l’évaporation d’une partie de l’eau.
L’air entre avec un peu d’humidité et il devient saturé à la sortie après qu’il ait absorbé une
quantité d’eau, pour cela on a effectué un bilan de matière pour déduire la quantité d’eau
évaporée par la relation suivante :
G . X1 + E = G . X2 E = G (X2- X1)
Où :
H entrée = H sortie
Q - Δt
G = —————————
(Hs - He) - td ( Xs – Xe)
G = 8,694.10 kg/h
G = 7,286.10 m³/h
E = G ( X2 – X1)
E = 115,64 m³/h
Calcul de la quantité de la chaleur évacuée au niveau de la tour 33 E101 :
Qc = m C p . ΔT
Qc = 5200 .10³.1.15
Remarque :
1 thermies = 1000 Kcal
Q
EV = ————
1000
EV = 5,2 m³/h
A = D + E + EV + fuites
Où :
C = déconcentration .
A = débit d’appoint.
E
D= ————
Et C-1
Salinité du circuit
C = —————————
Salinité d’appoint
( CT ) du circuit
C = ————————— « C » varie entre : 1≤ C ≤ 3
( CT ) d’appoint
C = 2,078
115,64
D= ————
D = 107,27 m³/h
2,078 - 1
Nous supposons que les fuites sont négligeables par rapport au débit circulé.
A = E + D + EV
A = 228,11 m³/h
Remarque :
L’alimentation de la tour par l’eau brute a nécessité une déconcentration d’une seule fois en
trois mois.
Donc : A = E + EV
A = 120,84 m³/h
Ecart de température :
Où 1000 représente le nombre de calories à enlever à l’eau par Kg d’air et par degré.
W = 7,8 10 Kcal /h
640 - tm
Ce = ( Xs – Xe) ————
1000
L’eau à refroidir est finement pulvérisée pour offrir le maximum de surface de contact avec
l’air .
Le refroidissement de l’eau se fait par contact avec l’air atmosphérique fourni par l’action
directe d’un ventilateur,les deux actions qui rentrent en jeu sont :
a-l’air se réchauffe au contact de l’eau à laquelle il enlève ainsi une certaine quantité de
chaleur.soit « te », la température de l’air au contact de l’eau ; sa température devient « ts »
et il a enlevé une quantité de chaleur Cc = 0,2375( ts – te), la claleur spécifique de l’air étant Cp
= 0,2375 cal /gr,à pression atmosphérique.
b-l’air atmosphérique est rarement saturé, il absorbe donc de la vapeur d’eau, d’autant plus
que sa température s’est élevée par le contact avec l’eau à refroidir.
Donc en tout 1kg d’air a élevé a l’eau (on aura aussi des petites calories) :
640 - tm
Cc + Ce = 0,2375 ( Ts – Te) ÷ ( Xs – Xe) ————
1000
On a :
42 + 27
640 - tm Où tm = ———— = 34,5°C
2
Ce = ( Xs – Xe) ———— 2
1000
640 – 34,5
Ce = ( 0,0271 – 0,0138 ) —————
1000
Ce = 8,05315 cal / Kg
Cc = 1,1875 cal / Kg
Donc en tout 1 Kg d’air a enlevé à l’eau :
640 - tm
Cc + Ce = 0,2375 ( T2 – T1) + ( X2 – X1) ————
1000
Cc + Ce = 9,24065 cal / Kg
Cc + Ce = 69,6.103 Kcal / h
Le rendement du réfrigérant
ta - td
ŋ = ——————
ta - Th
ŋ = 76,92 %
Remarque :
Nous pouvons conclure que le rendement est proportionnel à l’augmentation
De la température entre l’entrée et la sortie de l’eau de refroidissement.
ŋ = f ( ∆t )
1-Introduction :
Le choix d’une méthode de calcul d’un réfrigérant atmosphérique est basé sur plusieurs
facteurs tels que l’humidité, le débit d’air aspiré, la surface d’échange et le débit d’eau à
refroidir.
Le choix s’est porté sur cette méthode parce qu’elle est la plus utilisée. Elle permettra de
donner des valeurs appropriées des différentes caractéristiques thermiques du réfrigérant
atmosphérique de type contre courant à tirage forcé.
Etant donné que la quantité d’eau transportée par l’écoulement de l’air à partir des gouttelettes
d’eau est supposée négligeable et que le Cp de l’eau est égal à l’unité.
b) Transfert de chaleur :
Il a été démontré que dans toute portion de volume dv de la tour ayant une surface « a » par
unité de volume, la chaleur dégagée par kilogramme d’eau est mise en relation avec la chaleur
transférée à l’air par convection et à la chaleur perdue par vaporisation selon la relation
suivantes :
Où :
∆Hv : L’enthalpie de vaporisation de l’eau ( Kcal/ Kg).
a : La surface d’échange (m2 /m3).
dv : Variation de volume.
Xs : La quantité d’eau dans l’air saturé à la sortie ( Kg d’eau / Kg d’air).
Xe : La quantité d’eau dans contenue dans l’air à l’entrée de la tour ( Kg d’eau / Kg d’air).
K : Coefficient de transfert de chaleur .
Cp : La chaleur spécifique de l’air ( Kcal/ Kg .°C).
h : Coefficient de transfert de chaleur par convection .
Remarque :
Les résultats des tests pratiques [4] indiquent que la Coefficient de transfert de chaleur perdue par
vaporisation à partir de l’eau vers l’écoulement de l’air est approximativement égal au coefficient
conventionnel de transfert de chaleur par convection (h) divisé par le coefficient de chaleur
spécifique de l’air (Cp),soit :
K = h / Cp
Lw dtw = K ( Hs - Ha ) . a . dv - - - - - - - - - - - -(3)
∫ dtw ∫ K . a . dv K.a.V
————— = ——————— = ——————— - - - - - - - - - - - - - (4)
Hs – Ha Lw Lw
Nous pouvons constater d’après l’équation (4) la variation d’enthalpie ne varie pas
linéairement avec la température de l’eau ; le graphe V.1 ci-après montre cette variation.
Entre autre sur ce même graphe il est tracé la variation linéaire de l’enthalpie de l’air humide
en fonction toujours de la température de l’eau , nous obtenons de ce graphe pour chaque
valeur de température de l’eau, la différence d’enthalpie (Hs – Ha) telle que l’indique
l’intégration de l’équation (4) .
Fig V.1 : Effets de la température de l’eau sur l’enthalpie du débit d’air dans une
tour de refroidissement à contre courant.
Une seconde courbe reprenant [1/ (Hs – Ha)] est tracée en fonction de la température de l’eau tw
et nous obtenons la courbe suivantes :
Remarque :
La position de cette droite pourrait être définie par l’introduction d’un facteur de correction
δh donné par l’équatiion (5).
Remarque :
La valeur du facteur de correction δh est donnée par la Fig V.4 qui représente la variation de
δh en fonction de la température de l’eau.
Où A et n sont des constantes pour les différents types de garnissages qui sont donnés dans le
tableau V.1 et N est le nombre de garnissages contenus dans la tour de refroidissement.
Notons que la caractéristique de la tour (Ka V/ Lw) est proportionnelle à la hauteur de
l’ensemble de garnissages dans la tour.
Type géométrie de A n
garnissage
A 0,06 0,99
B 0,07 0,99
C 0,092 0,98
D 0,019 0,92
E 0,110 0,86
F 0,100 0,79
G 0,104 0,82
H 0,127 0,84
I 0,135 0,87
J 0,103 0,81
Remarque :
La géométrie du garnissage tel que le représente la figure en annexe V, est une
variable dans la détermination de la caractéristique de la performance (K.a .V/ Lw).
Pi εi 18
Xi = —————— ——
P atm – Pi εi 29
Hs = Cp . Ts + ∆Hv . Xs
Ha = Cp . Ta + ∆Hv . Xe
Nous calculons [1/ (Hs – Ha)] dont les résultats sont donnés dans le tableau V.3 ci-après :
MATIN SOIR
Cellule1 Cellule2 Cellule3 Cellule4 Cellule1 Cellule2 Cellule3 Cellule4
Ta 25,67 25,66 25,68 25,73 28,87 28,98 28,98 29,00
%s 58,06 58,47 58,20 58,27 51,21 51,22 51,11 51,12
P1 24,720 24,705 24,735 24,808 29,818 30,008 30,008 30,043
Xe 0,01194 0,01202 0,01198 0,01203 0,01272 0,01281 0,01278 0,01280
Ts 29,10 28,45 28,82 29,02 32,25 31,79 31,18 32,17
%ε 100,0 100,0 100,0 100,0 100,0 100,0 100,0 100,0
P2 30,217 29,099 29,731 30,077 36,170 35,241 36,027 36,007
Xs 0,02569 0,02471 0,02526 0,02557 0,03101 0,03018 0,03088 0,03086
Ha 13,4091 13,4552 13,4358 13,4785 14,6697 14,7516 14,7333 14,7504
Hs 22,6209 21,6808 22,2947 22,5277 26,6609 26,0376 26,5638 26,5481
∆H 9,2118 8,4056 8,8589 9,0492 11,9912 11,2860 11,8305 11,7987
1 0,1085 0,1189 0,1128 0,1105 0,08339 0,08860 0,08452 0,08475
———
(Hs- Ha )
Par les équations suivantes, nous calculons ∆Hm ( moyenne logarithmique ), dont les résultats
sont donnés dans le tableau ci-après.
Avec :
∆H1 = Hs1 – Ha1
∆H2 = Hs2 – Ha2
Hsm = ( 0,24 + 0,47 . Xsm ) Tsm + 595 . Xsm
Où :
Ha1 : L’enthalpie de l’air à l’entrée au matin ( Kcal / Kg).
Ha2 : L’enthalpie de l’air à l’entrée au soir ( Kcal / Kg).
Hs1 : L’enthalpie de l’air saturé à la sortie au matin ( Kcal / Kg).
Hs2 : L’enthalpie de l’air saturé à la sortie au soir ( Kcal / Kg).
∆Hm : La moyenne logarithmique de L’enthalpie de l’air ( Kcal / Kg).
δh : Facteur de correction de L’enthalpie ( Kcal / Kg).
Hsm : L’enthalpie moyenne de saturation entre le matin et le soir ( Kcal / Kg).
Tsm : Température moyenne de saturation entre le matin et le soir ( °C).
Xsm : L’humidité absolue moyenne de saturation entre le matin et le soir ( Kg d’eau / Kg d’air).
Ka .V / Lw ( Lw /G)moy ηth
Cellule1 1,251 0,8376 66,95
Cellule2 1,286 0,8056 62,64
Cellule3 1,272 0,8219 64,66
Cellule4 1,269 0,8235 64,89
Tableau V.6 : valeurs du rendement théorique.
Les résultats des calculs de la caractéristique ( Ka .V / Lw) sont donnés dans le tableau V.6 ci-
après :
MATIN SOIR
3
Lw (10 ) G (106 ) Ka .V / Lw 3
Lw (10 ) G (106 ) Ka .V / Lw
Cellule1 1300 1,786 1,215 1300 1,372 1,553
Cellule2 1300 1,898 1,149 1300 1,405 1,519
Cellule3 1300 1,861 1,170 1300 1,375 1,550
Cellule4 1300 1,821 1,193 1300 1,393 1,531
Il est démontré que le rendement de la tour de refroidissement est en relation liée avec la
caractéristique de performance selon :
η = ( Lw /G) / ( Ka .V / Lw)
MATIN SOIR
Lw / G K aV / Lw η Lw / G KaV / Lw η
Cellule1 0,7278 1,143 63,67 0,9475 1,463 64,76
Cellule2 0,6849 1,080 63,41 0,9252 1,431 64,65
Cellule3 0,6985 1,100 63,50 0,9454 1,460 64,75
Cellule4 0,7138 1,122 63,61 0,9332 1,442 64,71
Ka .V / Lw = 0,07 + A . N( Lw /G)ⁿ
MATIN SOIR
Lw / G K . a . V / Lw η Lw / G K . a . V / Lw η
Cellule1 0,7278 1,143 63,67 0,9475 1,463 64,76
Cellule2 0,6849 1,080 63,41 0,9252 1,431 64,65
Cellule3 0,6985 1,100 63,50 0,9454 1,460 64,75
Cellule4 0,7138 1,122 63,61 0,9332 1,442 64,71
III – Conclusion :
IV – Tableau récapitulatif :
Le tableau ci-dessous ( V.10 ), donne les valeurs exprimées dans la méthode précédente :