UNIVERSIDAD NACIONAL DE CRDOBA

FACULTAD DE CIENCIAS EXACTAS,

FSICAS Y NATURALES

TRANSFERENCIA DE
CALOR Y MASA
Torres de enfriamiento
Trabajo Prctico

PROFESORES
Ing. Sirena Alberto
Ing. Galeasso Angel

ALUMNOS
AMBROSINO, Gastn
32683182
FERNANDEZ HERBAS, Pablo Alejandro
32408702
GAY, Gabriel
32795937

MATRICULA
MATRICULA
MATRICULA

FCEFyN

Transferencia de calor y Masa

CARRERA: Ingeniera Mecnica

AO 2009

INTRODUCCIN
Las torres para enfriar agua son aparatos que tienen la funcin de
disminuir la temperatura del agua que circula por su interior y en las cuales el
proceso de transferencia de calor, entre el medio lquido caliente (agua) y el
fluido refrigerante (aire hmedo), se realiza por contacto directo entre ambos
fluidos. En las torres el proceso de transferencia de calor se realiza en parte
por un proceso de transferencia de calor sensible entre la fase liquida y la
gaseosa y por otra parte debido a una transferencia de calor latente
(formacin de vapor de agua en la interfase liquido/gas) como consecuencia
de un fenmeno de difusin de masa, es decir por transferencia de masa, el
primero se debe a un diferencia de temperaturas y el segundo a la diferencia
de concentracin de vapor de agua en el aire.
La utilizacin de torres de enfriamiento produce un impacto directo en
la eficiencia de los procesos en los cuales participa, como es le caso de la
generacin de energa elctrica. Se utilizan estas torres cuando es necesario
enfriar agua a un costo razonablemente reducido y el salto trmico requerido
no sea elevado (5 a 10 C).
Cuando el caudal de agua que se necesita enfriar es muy grande se
pueden disponer de torres de enfriamiento en paralelo o bien cuando el salto
trmico es elevado se suelen instalar torres de enfriamiento en serie.
En general las torres de enfriamiento de agua estn conformadas por
una carcasa exterior o gabinete en el cual se posicionan los distintos
componentes que la integran, a saber: sistema de control y distribucin de
agua que incluye las toberas o picos rociadores, sistema impulsor del aire
(motor y ventilador), relleno, separagotas y la cuba o pileta de recoleccin del
agua fra.

CONCEPTOS GENERALES
En trabajo prctico que se realiza a continuacin sobre una torre de
enfriamiento de agua consta de dos partes. En la primera se realiza en
ensayo sobre la torre desde el punto de vista del fabricante donde se busca
obtener la curva de funcionamiento de la torre (nmero de unidades de
difusin vs. relacin de caudales agua-aire) que da una ideal de las
prestaciones termodinmicas de la misma. La segunda parte es el tipo de
ensayo que realiza un comprador de una torre de enfriamiento, donde busca
verificar si la misma posee realmente la capacidad de enfriamiento que el
fabricante especifica.

Pag 2

FCEFyN

Transferencia de calor y Masa

ESQUEMA DE LA TORRE DE ENFRIAMIENTO

Para el ensayo en estudio utilizaremos la torre mecnica de tiro
inducido de flujo en contracorriente. El flujo a contracorriente significa que el
aire se mueve verticalmente a travs del relleno, de manera que los flujos de
agua y de aire se realizan en la misma direccin pero en sentido opuesto y los
ventiladores estn ubicados en la zona de salida del aire. El ingreso de aire se
realiza por aberturas laterales ubicadas en la parte inferior de las torres. El
rea de ingreso del aire debe ser lo suficientemente grande para lograr un
flujo de baja velocidad, evitando de esta manera la aspiracin de cuerpos
extraos y una perdida de carga elevada.
La ventaja que tiene este tipo de torres es que el agua ms fra se pone
en contacto con el aire ms seco, logrndose un mximo rendimiento. El
relleno que posee la torre es de tipo de pelcula, ste empaque tiene como
misin mejorar el proceso de transferencia de calor incrementando la
superficie de intercambio y aumentando el tiempo de contacto entre el agua
y el aire.
En el interior de la torre se tienen zonas bien diferenciadas bajo el
punto de vista de las condiciones de funcionamiento: el relleno y dos zonas
de flujo libre de agua, una desde los picos rociadores hasta el relleno y la otra
desde ste ltimo hasta la pileta colectora, a esta ltima zona se la conoce
como zona de lluvia.
El esquema de la instalacin de la torre es el siguiente, donde se puede
observar la distribucin de los sensores de presin caudal y temperatura:

Pag 3

FCEFyN

Transferencia de calor y Masa

Tanto la temperatura del bulbo hmedo como la temperatura del bulbo

seco, son medidas a un lado de la instalacin torre, ya que ambas dependen
de las condiciones atmosfricas del lugar donde se realizan los ensayos.

Pag 4

FCEFyN

Transferencia de calor y Masa

PRIMERA PARTE:
(fabricante)

Obtencin

de

la

curva

de

funcionamiento

Las mediciones realizadas durante el ensayo fueron las siguientes:

N
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43

Thw
[C]
32,176
33,647
34,199
34,468
34,648
34,777
34,857
34,93
34,981
35,063
35,248
35,532
35,496
35,474
35,452
35,452
35,518
35,539
35,554
35,583
35,627
35,656
35,656
35,656
35,656
35,656
35,656
35,656
35,612
35,612
35,598
35,627
35,598
35,598
35,576
35,569
35,569
35,569
35,561
35,539
35,554
35,532
35,292

Tcw
[C]
21,801
22,688
22,922
22,987
23,109
23,064
23,109
23,039
22,842
23,169
23,169
23,522
24,688
25,31
25,656
25,432
25,773
25,365
25,308
25,173
25,343
25,395
25,417
26,134
26,383
26,555
26,573
26,286
26,286
26,458
26,222
26,478
26,508
27,071
27,098
26,896
27,566
27,357
27,32
26,58
22,82
22,476
22,658

Twb
[C]
10,345
10,291
10,276
10,276
10,261
10,253
10,238
10,268
10,322
10,337
10,337
10,337
10,337
10,337
10,276
10,276
10,276
10,276
10,246
10,246
10,261
10,261
10,246
10,246
10,246
10,246
10,268
10,246
10,246
10,246
10,238
10,215
10,215
10,192
10,144
10,174
10,184
10,184
10,184
10,144
10,121
10,113
10,106

Tdb
[C]
19,269
19,231
19,216
19,193
19,147
19,109
19,078
19,063
19,033
19,017
18,997
18,946
18,908
18,885
18,847
18,77
18,618
18,47
18,279
18,142
18,111
18,134
18,172
18,233
18,264
18,218
18,124
17,994
17,918
17,887
17,91
17,88
17,686
17,549
17,488
17,442
17,412
17,373
17,305
17,205
17,175
17,198
17,236

L
[l/h]
5583,9
5670
5662,2
5687
5665,9
5626
5684,8
5674,4
5637
5665,9
5637,3
7226
8039,3
8067,9
8017,6
7983,6
7904,7
7948,1
7937,7
7792,1
7976,1
7933,6
9158,1
9743,1
9866
9828
9799,4
9807,2
9717,3
9782,4
9700,9
9742,1
10741
11482
11519
11588
11512
11532
11489
5446,7
1926,6
1835,7
728,05

Dp
Patm
[Pa]
[Pa]
-138,6 96995
-135,6 96993
-137,9 96991
-136,4 96992
-135,9 96993
-142,1 96995
-149,1 97006
-144,7 97004
-140,7 96997
-140,5 97000
-148,8 96997
-150,4 96989
-149,8 96989
-153,2 96982
-124,6 96987
-145,6 96988
-148,8 96996
-145,9 96992
-140,9 96993
-140,9 96993
-146,6 96995
-143,2 96995
-146,2 96994
-149,7 96991
-149,4 96986
-155,5 96985
-148,9 96986
-149,8 96985
-151,6 96976
-150,2 96973
-152 96977
-150,5 96976
-152,4 96987
-150,1 96977
-155,8 96979
-149,4 96985
-148,5 96987
-150,6 96994
-151,8 96994
-137,3 97003
-124,3 97010
-126,9 97013
-122,8 97013

Pest
[Pa]
-268,1
-273,2
-275,8
-277,3
-282,0
-281,9
-285,9
-285,4
-289,7
-283,2
-290,4
-291,5
-291,9
-291,2
-220,2
-281,5
-285,5
-281,0
-278,3
-283,9
-284,4
-279,9
-285,9
-284,4
-291,7
-287,7
-285,2
-286,9
-287,5
-289,7
-287,4
-285,3
-290,6
-289,8
-288,6
-289,1
-288,2
-288,7
-286,3
-273,0
-264,2
-260,5
-260,5

Frecuencia del motor: 45 [Hz]

Potencia consumida: 1860 [W]
Pag 5

FCEFyN

Transferencia de calor y Masa

En donde:
Thw: temperatura del agua caliente.
Tcw: temperatura del agua fra.
Twb: temperatura del bulbo hmedo.
bulbo seco
Dp: diferencia de presin antes y despus
salida.

L: caudal de agua.
Patm: presin atmosfrica.
Tdb: temperatura del
del cambio de seccin a la

Si se grafica el caudal de agua segn las diferentes mediciones, puede

observarse que existen 4 zonas de la curva en donde est bien definido un
valor de caudal:

Promediado los tramos rectos y despreciando los tramos ascendentes, se

obtienen 4 valores de caudal a los cuales se les asigna un valor promedio de
todas las mediciones hechas:
Caudal 1: L1 = 5654 [litros/hora]
Thw1
[C]

Tcw1
Twb1
Tdb1
L1
Dp1
Patm1
[C]
[C]
[C]
[l/h]
[Pa]
[Pa]
22,89990 10,29127
5654,036
96996,63
34,454
9
3
19,123
4
140,9364
6

Caudal 2: L2 = 7951 [litros/hora]

Thw2
Tcw2
Twb2
Tdb2
L2
[C]
[C]
[C]
[C]
[l/h]
35,5351 25,3443 10,2792 18,5164 7960,07

Dp2
[Pa]
-143,95

Patm2
[Pa]
96991

Dp3

Patm3

Caudal 3: L3 = 9808 [litros/hora]

Thw3

Tcw3

Twb3

Tdb3

L3

Pag 6

FCEFyN

Transferencia de calor y Masa

[C]
[C]
[C]
[C]
[l/h]
[Pa]
[Pa]
35,63655
10,24411 18,04755 9776,266
96981,66
6
26,375
1
6
7
32181,17
7

Caudal 4: L4 = 11520 [litros/hora]

Thw4
Tcw4
Twb4
Tdb4
[C]
[C]
[C]
[C]
35,57714 27,11657 10,18242
3
1
9
17,465

L4
[l/h]
11409

Dp4
Patm4
[Pa]
[Pa]
96986,14
151,2286
3

Otra grfica importante de ver, es como evolucionan las temperaturas del

agua caliente y fra del agua a lo largo del ensayo:
Evolucin de las temperaturas del agua
40
35
30
T [C]

25
Thw

20

Tcw

15
10
5
0
0

10

15

20

25

30

35

40

45

50

N de mediciones

Es necesario corregir el valor a cero de los sensores, porque al

cambiar su temperatura cambia su cero. Para ello se calibra el sensor,
teniendo en cuenta los valores de los errores a cero Dp previo y posterior al
ensayo. La forma de hacerlo es tomar un promedio de las mediciones previas
al ensayo y otro de las mediciones posteriores al ensayo, y luego hacer un
nuevo promedio con estos subpromedios anteriores de forma tal de obtener
el valor de correccin Dp0.

Prom

Dp
Dp
posterior
anterior
[Pa]
[Pa]
-9,635
-8,354
-9,767
-8,821
-9,731
-8,749
-9,695
-8,929
-9,647
-8,57
-9,516
-8,666
-9,444
-8,845
-9,504
-9,013
-9,3
-8,977
-9,1757

Dpo= -9,1757
Finalmente, se corrige el valor de presin medido de la forma:
Pag 7

FCEFyN

Transferencia de calor y Masa

Dp corregido = Dp medido Dp 0

Finalmente, los valores corregidos son los siguientes:

Dp1
Dp2
Dp3
Dp4
[Pa]
[Pa]
[Pa]
[Pa]
-131,761 -134,774 -141,669 -142,053

Para obtener la curva de funcionamiento, es necesario calcular el

caudal de aire Gi para cada caudal de agua Li y as poder obtener las
relaciones de caudales L/G.
Operamos de la siguiente manera:
La presin dinmica antes del ventilador vale:
q =1,5 Dp [Pa ]

y a su vez tambin es igual a

q=

1
V 2
2

Donde es la densidad del aire [kg/m3] y V es la velocidad [m/s].

Utilizando la ecuacin de estado P = R T, donde P es la presin atmosfrica
[Pa],
R = 286,98 [Pa m3/K kg] es la constante del aire y T la temperatura
del bulbo seco [K], se calcula la densidad del aire .
As, la velocidad del aire vale: V =

2q

[m s ]

kg
Luego, el caudal de aire G se calcula como: G = 3600 A V h

Donde A es el rea de la seccin de salida de dimetro d = 0,45m, por lo que:

A=

[ ]

d2
= 0,15904 m 2
4

Finalmente, se calcula la relacin de caudales L/G, donde L [Kg/h].

Para el primer caudal:

q = 1,5 (131,76 ) = 197,64 [ Pa ]
V =

L = 5654,03 kg
h

[ ]

2 (197,64 )
= 18,48 m
s
1,157

96996,6
= 1,157 kg 3
m
286,98 (19,12 + 273)

G = 3600 1,157 0,159 18,65 = 12238,3 kg

h

L
5654,06
=
= 0,462
G 1 12238,3

Para el segundo caudal:

Pag 8

FCEFyN

Transferencia de calor y Masa

[ ]

G = 3600 1,159 0,159 18,67 = 12390,8 kg

h

V =

2 ( 202,2)
= 18,67 m
s
1,159

L = 7960,07 kg
h

L
7960,07
=
= 0,642
G 2 12390,8

Para el tercer caudal:

96981,6
= 1,161 kg 3
m
286,98 (18,04 + 273)

q = 1,5 (141,67 ) = 212,5 [ Pa ]

[ ]

G = 3600 1,161 0,159 19,14 = 12719,6 kg

h

V =

2 ( 212,5)
= 19,14 m
s
1,161

L = 9776,2 kg
h

96991
= 1,159 kg 3
m
286,98 (18,51 + 273)

q = 1,5 (134,77 ) = 202,2 [ Pa ]

L
9776,2
=
= 0,768
G 3 12719,6

Para el cuarto caudal:

q = 1,5 (142,05) = 213,77 [ Pa ]
V =

[ ]

2 ( 213,77 )
= 18,9 m
s
1,163

L = 11409 kg
h

97008,2
= 1,163 kg 3
m
286,98 (17,46 + 273)

G = 3600 * 1,163 * 0,159 * 19,3 = 12758,7 kg

h

L
11409
=
= 0,891
G 4 12758,7

El nmero de unidades de difusin K para cada caudal, se calcula

utilizando el software KAVLC, ingresando los parmetros en la siguiente
secuencia: Twb, Thw, Tcw, L/G y Patm, con la salvedad de que la presin
atmosfrica de ser ingresada en mm de Hg. Se obtienen as los valores para
cada caudal:
Caudal
Caudal
Caudal
Caudal

1
2
3
4

K
K
K
K

=
=
=
=

0,9058
0,7002
0,6085
0,5417

Con las relaciones de caudal L/G y sus respectivas unidades de difusin K se

realiza la curva de funcionamiento de la torre:

Curva de Funcionamiento de la Torre

Pag 9

FCEFyN

Transferencia de calor y Masa

Pag 10

FCEFyN

Transferencia de calor y Masa

SEGUNDA PARTE: Ensayo de recepcin (comprador)

Con los parmetros de diseo de la torre dados por el fabricante, se
construye la curva de demanda. Para ella es necesario introducir diferentes
valores de L/G al programa KAVLC y obtener los respectivos valores del
nmero de unidades de difusin K.
Las condiciones de diseo segn el fabricante son:
L/G = 0,83
L = 8500 [l/h]
Thw = 34 [C] = 93,2 [K]
Tcw = 26 [C] = 78,8 [K]
Twb = 17 [C] = 62,6 [K]
Tdb = 25 [C] = 77 [K]
Patm = 760 [mmHg]
N = 1860 [W]
Utilizando el software KAVLC, se construy la tabla de los K y con sta la
curva de demanda dando valores cualesquiera a L/G y en las condiciones de
diseo, como se ve a continuacin:
L/G
0,458
0,631
0,759
0,889

k
0,905
0,713
0,620
0,546

Curva de Funcionamiento de la Torre y Curva de la Demanda

Pag 11

FCEFyN

Transferencia de calor y Masa

Dado que el ensayo se inici a una frecuencia de motor de 40 [Hz] y a

partir de la medicin N 25 se continu con 50 [Hz], se debern analizar los
dos ensayos por separado y luego juntar los resultados para hacer una
conclusin global.

REDUCCIN A CONDICIONES DE DISEO.

1. Ensayo con una frecuencia de motor de 40 Hz
Segn las normas, la relacin de caudales del ensayo cumple con la siguiente
expresin:

( )

( )

3
L
W
3 vt L
L
= t d t
G t Ld Wt d v d G

(1)
d

Donde el subndice t significa que los parmetros son del ensayo y el

subndice d significa que son de diseo. L es el caudal de agua, W es la
potencia del motor, es la densidad del aire hmedo y v es el volumen
especfico del aire hmedo.
Pag 12

FCEFyN

Transferencia de calor y Masa

Por otro lado, utilizamos la siguiente ecuacin de balance entlpico:

( G) T

h2 = L

agua

+ h1

(2)

Donde:
h1=entalpa del aire a la entrada.
h2=entalpa del aire a la salida.
Con los valores de diseo se calculan d y vd:
Con Twb = 17 [C] = 62,6 [F] se entra en la tabla F17 y se encuentra h 1 =
28,28 [BTU/lb].
Utilizando la ecuacin (2), se obtiene h2:

( G)

h2 = L

( Thw d Tcw d ) + h 1 = 0,83 ( 93,2 78,8) + 28,28 = 40,23 [ BTU lb]

Con h2 se vuelve a entrar en la tabla F17 y se obtiene una nueva temperatura

del bulbo hmedo de salida Twb = 77,7 [F], con la que entrando en la tabla
H19 y J21 se obtienen la densidad y el volumen especfico del aire
respectivamente:

d = 0,07293 lb 3
ft

3
v d = 13,993 ft
lb

Los valores promedios obtenidos durante el ensayo son los siguientes:

Thw
Twb
Patm
[C] Tcw [C] [C] Tdb [C] L [l/h]
[Pa]
34,0369
6
25,257 9,54256 17,43 8043,6 96923,2

WTotal
[watts]
1140

Con los dos valores calculados junto con las condiciones de diseo y los
valores obtenidos en el ensayo, se puede expresar la ecuacin (1) de la
forma:

( L G ) = 0,15816
t

1
3

vt

(3)

La ecuacin (3) se resuelve en forma iterativa hasta encontrar en

valor de L/G de ensayo que converja. La secuencia que debe seguirse es la
siguiente, se supone un valor de Twb de entrada, con ese valor de entrada se
calculan t y vt con las tablas H19 y J21 respectivamente. Luego se calcula L/G
con la ecuacin (3) y con este valor h 2 con la ecuacin (2) donde h1 = 19,85
que sale utilizando la tabla F17 con Twb = 9,542 [C] = 49,18 [F] y es un
valor que no cambia.

Pag 13

FCEFyN

Transferencia de calor y Masa

Por ltimo, se vuelve a entrar al la tabla F17 con h 2 y se obtiene un nuevo

Twb que se toma como nuevo valor de entrada. El proceso converge con L/G.
Se toma como valor de entrada:
Twb =

Thw + Tcw 34,03 + 25,25

=
= 29.64 [ C] = 85.35 [ F]
2
2

Los pasos de iteracin se resuelven en la siguiente tabla:

Twb
85,35
69
68,5

0,0716
0,07438
0,07446

v
14,32
13,65
13,631

L/G
0,94
0,908
0,907

h2
33,38
32,92

Twb
69
68,5

As se ve que el valor de L/G de ensayo vale 0,907. Con este valor junto al los
restantes parmetros de ensayo se entra en el software KAVLC y se obtiene el
K de ensayo:

( L G) = 0,907
t

k t = 0,605

Lo que sigue consiste en entrar con estos parmetros en el grfico

presentado en al comienzo de la segunda parte y obtener en l un punto
llamado punto de ensayo o de funcionamiento. Luego se taza una paralela al
la curva de funcionamiento que pase por ste punto y se lee el valor de L/G
ajustado en el eje de abcisas correspondiente al punto de corte de la curva
paralela y la curva de demanda:

Curva de Funcionamiento de la Torre

Finalmente, para calcular la capacidad de la torre se emplea la siguiente

expresin:

Pag 14

FCEFyN

( L G)
Cap. torre =
( L G)

Transferencia de calor y Masa

ajustado

diseo

0,617
= 0,745 = 74,5%
0,83

Al tener una capacidad de 0,745, la torre no alcanza por un 25,5% las

prestaciones con las que fue solicitada.

2. Ensayo con una frecuencia de motor de 50 Hz

Los procedimientos a seguir son los mismos que en el caso de 40 Hz.
Los valores promedios obtenidos durante el ensayo son los siguientes:
Patm
Thw [C]Tcw [C]Twb [C] Tdb [C] L [l/h]
[Pa]
34,2322 23,6264 9,43224 16,0622 7971,06 96922,9
76
83
14
07
21
66

WTotal
[watts]
2590

La ecuacin (1) los nico parmetro que cambia es W t que ahora vale 2590
[watts], y Lt por lo que la ecuacin (3) ahora tiene la forma:

(L G)

= 0,11922 t 3 v t

(4)

Tomando nuevamente Twb de entrada como un promedio entre Thw y Tcw

del ensayo y h1 = 19,75 que sale utilizando la tabla F con Twb = 9,432[C] =
48,98 [F] se vuelve a iteratuar para obtener el valor de L/G de ensayo:

v
0,071 14,26
84,06
91
3
0,0756
61,1
6
13,368
0,0756
61,6
1
13,38
Twb

L/G

h2

Twb

0,707

27,25

61,1

0,674

26,9

61,6

0,674

Con L/G junto a los restantes parmetros de ensayo se obtiene K utilizando

KAVLC:

( L G) = 0,674
t

k t = 0,824

Continuando con la misma secuencia, se obtiene el parmetro grficamente

para luego calcular la capacidad de la torre:

Pag 15

FCEFyN

Transferencia de calor y Masa

La capacidad de la torre ser:RECALCULAR

( L G)
Cap. torre =
( L G)

ajustado
diseo

0,64
( RECALCUAAR ) = 0,818 = 81,8%
0,83

Curva de Funcionamiento de la Torre

La torre no alcanza por un 18,2% las prestaciones.

Para concluir decimos que, dado que la capacidad de la torre fue menor que
1 ensayndola con ambas frecuencia de motorizacin (40 [Hz] y 50 [Hz]), la
torre adquirida no cumple la capacidad requerida.

Pag 16