1.- FUNDAMENTO TERICO 1.1.

- INTRODUCCIN La humidificacin es una operacin unitaria en la que tiene lugar una transferencia simultnea de materia y calor sin la presencia de una fuente de calor externa. De hecho siempre que existe una transferencia de materia se transfiere tambin calor. Pero para operaciones como extraccin, adsorcin, absorcin o lixiviacin, la transferencia de calor es de menor importancia como mecanismo controlante de velocidad frente a la transferencia de materia. Por otro lado, en operaciones como ebullicin, condensacin, evaporacin o cristalizacin, las transferencias simultneas de materia y calor pueden determinarse considerando nicamente la transferencia de calor procedente de una fuente externa. La transferencia simultnea de materia y calor en la operacin de humidificacin tiene lugar cuando un gas se pone en contacto con un lquido puro, en el cual es prcticamente insoluble. Este fenmeno nos conduce a diferentes aplicaciones adems de la humidificacin del gas, como son su deshumidificacin, el enfriamiento del gas (acondicionamiento de gases), el enfriamiento del lquido, adems de permitir la medicin del contenido de vapor en el gas. Generalmente la fase lquida es el agua, y la fase gas el aire. Su principal aplicacin industrial es el enfriamiento de agua de refrigeracin, que ser el objeto de estudio de la prctica que nos ocupa. A grandes rasgos, el proceso que tiene lugar en la operacin de humidificacin es el siguiente:

una corriente de agua caliente se pone en contacto con una de aire seco (o con bajo contenido en humedad), normalmente aire atmosfrico. parte del agua se evapora, enfrindose as la interfase. el seno del lquido cede entonces calor a la interfase, y por lo tanto se enfra. a su vez, el agua evaporada en la interfase se transfiere al aire, por lo que se humidifica.

En la deshumidificacin, agua fra se pone en contacto con aire hmedo. La materia transferida entre las fases es la sustancia que forma la fase lquida, que dependiendo de cmo estemos operando, o se evapora (humidificacin), o bien se condensa (deshumidificacin.) Existen diferentes equipos de humidificacin, entre los que destacamos las torres de enfriamiento por su mayor aplicabilidad. En ellas, el agua suele introducirse por la parte superior en forma de lluvia provocada, y el aire fluye en forma ascendente, de forma natural o forzada. En el interior de la torre se utilizan rellenos de diversos tipos que favorecen el contacto entre las dos fases. . 1.2.- DEFINICIONES EN LA INTERACCIN AIRE-AGUA Antes de desarrollar las ecuaciones de diseo en una torre de enfriamiento, hemos de definir una serie de variables y conceptos involucrados en la operacin de humidificacin. - Humedad absoluta: es la razn msica de vapor de agua respecto al aire seco

(1) MA: masa molecular del agua; MA= 18 g/mol

MB: masa molecular del aire; MB= 28,9 g/mol PA: presin parcial que ejerce el vapor de agua en la mezcla gaseosa P: presin total (atmosfrica.) - Humedad relativa: es la relacin molar entre la cantidad de vapor de agua presente en el aire y la cantidad mxima posible (saturacin) para esa temperatura: (2) Pv es la presin de vapor, que podramos definir como la presin que ejerce un vapor en equilibrio con su lquido. La presin de vapor es funcin de la temperatura. Para un rango de temperaturas comprendido entre 0C y 50C, podemos obtener la Pv del agua mediante la ecuacin: (3) en la que Pv est expresada en mm Hg, y la T (temperatura) en K. - Calor especfico (cP): calor necesario para aumentar 1C a una unidad de masa de una sustancia, a presin constante. Tomamos los valores medios de los calores especficos del agua y del aire entre 0C y 100C: cPaire = cP,B = 1005 J/(kgC) cPvapor de agua = cP,A = 1884 J/(kgC) cPagua lquida = cL = 4180 J/(kgC) - Calor especfico hmedo (cS): (es el cP de la mezcla gaseosa aire-agua): cS = cP,B + cP,A H = 1005 + 1884H [J/(kgC)] (4) - Entalpa especfica del gas: calor asociado a un gas a cierta temperatura y presin (referencia 0C y agua lquida): (5)

donde = 2.502.300 J/kg agua (calor latente) - Entalpa especfica de saturacin: es la entalpa del gas saturado de humedad; la entalpa de saturacin es funcin de la temperatura, y cumple la ecuacin: (6) donde la T est expresada en C. - Condiciones de saturacin adiabtica: Considrese el proceso que se lleva a cabo de forma esquemtica en la figura 1, mediante el cual se busca que el gas se sature. Para ello hacemos circular agua a lo largo de

la torre en forma de lluvia; el agua entra en contacto directo con una corriente de aire no saturado. Para que el proceso sea adiabtico ha de cumplirse que:

no se transfiera calor a la torre la corriente lquida se recircule

De esta forma la temperatura a la entrada ser igual que a la salida. La diferencia de temperaturas entre el gas y el lquido tender a ser cero tanto ms cuanto nos acerquemos a la salida del gas, y no habr transferencia de calor sensible entre las fases. El nico efecto que ejercer el lquido en la torre es que parte de l se vaporizar y pasar a la corriente gaseosa. Si suponemos que la torre es suficientemente alta para que las fases lquida y gas alcancen el equilibrio en la parte superior de la misma, la fase gaseosa estar saturada, y TLentrada = TLsalida = TGsalida = TS . As pues, la temperatura Ts ser la temperatura de saturacin adiabtica, y Hs la humedad del gas saturado a Ts , es decir, la humedad de saturacin adiabtica. A partir de balances globales de materia y entalpa entre las condiciones inicial del gas (T y H a la entrada) y la condicin de saturacin adiabtica, se llega a la ecuacin:

(7) La ecuacin 7 relaciona la temperatura y humedad de un gas para cualquier condicin de entrada con las condiciones correspondientes para el mismo gas con su temperatura de saturacin adiabtica (slo se puede aplicar para esos dos puntos, no para describir la trayectoria seguida por el gas a medida que se satura). Si en nuestro sistema aire-agua consideramos que la Ts es igual a la temperatura de bulbo hmedo TH, y mediante las ecuaciones de transferencia correspondientes, se obtiene que (8.a)

Para el sistema aire-agua, se cumple que , por lo que la ecuacin 8.a es: (8.b) HG y TG son la humedad y temperatura del gas - Diagrama de humedad: tambin llamado carta psicomtrica, permite la obtencin mediante lectura directa de la mayora de las propiedades de las mezclas aire-vapor de agua que son necesarias en los clculos a realizar en la operacin de humidificacin, para una presin determinada. En la figura 2 se representa el diagrama de humedad para la presin de 1 atm.

Figura 2: Diagrama de humedad La limitacin del diagrama de humedad es la curva de saturacin. Los puntos representados a la derecha representan mezclas definidas de aire no saturado ( a la izquierda se representan mezclas de aire saturado y agua lquida.) Las curvas situadas entre el eje de abcisas y la curva de saturacin son las lneas correspondientes a diferentes humedades relativas. Las rectas inclinadas con pendiente negativa corresponden a las lneas de temperatura de saturacin adiabtica, que prcticamente coincide con la temperatura de bulbo hmedo para el sistema aire-agua, luego estas lneas siguen la ecuacin 8.b. 1.3.- OPERACIN DE HUMIDIFICACIN. MODELO MATEMTICO 1.3.1.- Proceso general de humidificacin Como hemos explicado en el apartado 1.1, en un proceso de humidificacin el agua caliente se pone en contacto con aire seco o poco hmedo; la materia (agua) y el calor se transfieren hacia la fase gaseosa a travs de la interfase, como se muestra en la figura 3, la cual se enfra.

Figura 3: Condiciones de la interfase en una altura determinada de una torre de humidificacin Como resultado, el lquido cede calor a la interfase, enfrindose tambin. Y el gas se humidifica. Podemos esquematizar el proceso de humidificacin por el dispositivo mostrado en la figura 4. L: caudal msico de lquido G: caudal msico de gas H: humedad del gas TL: temperatura del lquido TG: temperatura del gas q: calor transferido a la columna desde el exterior z: altura de la columna A: superficie interfacial a: rea interfacial del volumen de la columna (m2/m3) S: seccin transversal de la torre Figura 4. Dispositivo para el proceso general de humidificacin hL y hG son los coeficientes de transferencia de calor por conveccin de la fase lquida y gaseosa respectivamente. kG es el coeficiente de transferencia de materia de la fase gaseosa.

(es evidente que no hay coeficiente de transferencia de materia de la fase lquida, ya que en ella no hay transferencia de materia, al tratarse de un componente puro.) La columna opera adiabticamente (q=0) en rgimen estacionario. La aproximacin a la operacin adiabtica ser mayor a medida que el dimetro de la columna aumente. Escribiremos pues los balances para la altura diferencial dz. Balance de materia para el componente que se transfiere (agua): dL= GdH (9) Balance de entalpa: GdiG = d (L L) (10.a) donde L es la entalpa de la fase lquida. Como la proporcin de transferencia de agua es muy pequea comparada con la corriente del flujo total, utilizamos el valor promedio L, y el cambio de entalpa de la fase lquida L puede expresarse como la resultante del cambio de temperatura con calor especfico constante: d(L L ) = LCLdTL (10.b) por lo que el balance de entalpa queda como: LCLdTL = GdiG (10) Como la entalpa del gas sigue la ecuacin 5: GdiG = Gd(cSTG + H G 0) = GcSdTG + 0dH (11)

Consideremos ahora las ecuaciones de velocidad de transferencia. En la fase lquida, la transferencia de calor se produce desde el total de la fase lquida a la interfase, como resultado de una diferencia de temperaturas. Pero desde la interfase hasta la fase gaseosa hay dos tipos de transferencia de calor (como se indica en la figura 3):

transferencia debida a un gradiente de temperaturas (calor sensible) calor latente asociado con la masa que se transfiere como resultado de un gradiente de concentraciones (diferencia entre la humedad en la interfase y en el seno del gas.)

El calor transferido mediante estos dos mecanismos se separan en los dos trminos de la derecha de la ecuacin 11. Teniendo en cuenta estas ideas escribimos las ecuaciones de velocidad de transferencia para nuestro dz:

Transferencia de calor en la fase lquida: LcLdTL = hLa(TL-Ti)Sdz (12)

Transferencia de calor sensible en la fase gaseosa:

GcSdTG = hGa(Ti-TG)Sdz (13)

o
G

Transferencia de calor latente en la fase gaseosa: 0dH = 0kGa(Hi-H)Sdz (14)

Este modelo general es vlido para cualquier sistema gas-lquido. 1.3.2.- Ecuaciones de diseo En el caso del sistema aire-agua, podemos suponer que se cumple que

(15) Si sustituimos la relacin 15 en la ecuacin 13, y combinamos sta con la 11 y la 14, obtenemos que: Gd(cSdTG + (cSTG + 0dH) = kGa[(cSTi + 0 HG)]Sdz (16) 0 Hi) -

Vemos que los trminos de dentro de los parntesis son entalpas, definidas en la ecuacin 5, por lo que expresamos la ecuacin 16 como: GdiG = kGa (ii - iG )Sdz (17) Reordenando trminos e integrando desde la entrada de la torre hasta la salida, obtenemos que:

(18) El (NUT)G , nmero de unidades de transferencia, representa el nmero de veces que la fuerza impulsora promedio (iG2 - iG1) se divide entre el cambio de la entalpa. Es una medida de la dificultad de la transferencia de entalpa. (HUT)G es la altura de una unidad de transferencia. Podemos expresar la ecuacin 18, combinndola con la ecuacin 10, de esta otra forma. (19) Las ecuaciones 18 y 19 son las ecuaciones de diseo de una torre de humidificacin. En general, para el clculo de (HUT)G se deben los caudales y los parmetros de la columna que permiten determinar el coeficiente de transferencia. Para el clculo del nmero de unidades de transferencia se debe resolver la integral conociendo los valores de entalpa interfacial para que correspondan a los valores de entalpa del gas. Como los valores de entalpa interfacial dependen de la temperatura del lquido y de la entalpa del gas se deben tener una relacin entre estas variables a lo largo de la columna. Esta relacin se denomina recta de operacin, resultante del balance de entalpa (ec.10):

(20)

Por otro lado, igualando las ecuaciones 12 y 17, obtenemos la recta de unin:

(21) La recta de unin permite determinar la entalpa interfacial a partir de la temperatura del lquido y la entalpa del gas, asumiendo que la interfase se encuentra en equilibrio, es decir, que la entalpa interfacial es la entalpa de saturacin a la temperatura interfacial: ii = is (Ti) " ec.6. 1.3.3.- Resolucin de la integral (NUT)G. Mtodo Mickley El mtodo Mickley es un mtodo grfico para la obtencin de las condiciones de la interfase, mediante la ecuacin 21. Se basa en una grfica de entalpas de la fase gaseosa frente a las temperaturas de la fase lquida, como se representa en la figura 5. La curva de equilibrio es la curva de saturacin, que sigue la ecuacin 6, y que es la misma que la representada en el diagrama de humedad (figura 1). Figura 5. Representacin grfica de una operacin de contacto adiabtico gas-lquido. En esta curva estn localizadas todas las condiciones posibles de Ti, iGi a travs de la columna. La lnea ABC es la lnea de operacin, ecuacin 20, que contiene todos los valores de la iG correspondientes a la temperatura del lquido en toda la columna. Se obtiene a partir de las condiciones de iG, TL a la entrada y la salida de la torre. Sobre la lnea de operacin, el punto B representa un punto arbitrario dentro de la columna, en el cual la temperatura del lquido y la entalpa del gas tienen los valores TLn e iGn. Las condiciones de la interfase en este punto pueden encontrarse con la recta de unin, ec.21, por lo que una recta de unin que empiece en el punto B y que tenga una pendiente igual a -(hL/kG ) intersecar a la curva de equilibrio en las condiciones de la interfase que corresponde al punto B, en el punto I. De esta forma las condiciones de la interfase pueden encontrarse correspondiendo a cualquier punto entre A y C sobre la recta de operacin. Por lo tanto, si dividimos la recta de operacin en distintos intervalos (iguales) como se muestra en la figura 6, y los unimos con la curva de saturacin mediante lar rectas de unin, obtenemos los valores necesarios para la integracin numrica del NUT. Figura 6. Rectas de unin a lo largo de la columna. Paralelamente podremos determinar tambin la evolucin de la temperatura del gas a travs de la columna, como se indica en la figura 7. El punto F (TG1, iG1) representa la entrada del gas en la columna. Unimos este punto F con el punto D, que corresponde a las condiciones de la interfase a las condiciones de entrada de la columna. La recta FD permite determinar la evolucin de la temperatura del gas en un intervalo en el que se suponga que la pendiente de FD permanece constante. Esta TG es la correspondiente al punto G (TGG,iGG=iGH). En resumen, el punto G queda determinado por la recta FD, sobre la cual se encuentra, y el valor de la entalpa iG del punto H (lnea discontinua que une H y G), ya que los puntos H y G corresponden a una misma altura en la columna, por lo que tienen la misma iG. Procederemos de esta forma en todos los intervalos, hasta determinar TG2, temperatura del gas a la salida. As, la temperatura del gas a cualquier altura de la columna se determina a partir de los datos del intervalo anterior, siguiendo la ecuacin:

(22) La ec. 22 se alejar de la realidad conforme ms grandes sean los intervalos. Si no conocemos el valor de la relacin hL/ kG es decir, la pendiente de la recta de unin, tendremos que aplicar el mtodo de forma iterativa suponiendo su valor. El mtodo termina cuando la temperatura del gas a la salida de la torre calculada coincida con la obtenida experimentalmente. 2.- OBJETIVOS El objetivo de esta prctica es determinar los valores de:

coeficiente de transferencia de calor por conveccin de la fase lquida: hL coeficiente de transferencia de materia la fase gaseosa: kG nmero de unidades de transferencia en la fase gaseosa

en una torre de enfriamiento. Para ello se realizan experimentos de enfriamiento de agua en una columna a escala de laboratorio para diferentes caudales de agua y de aire. 3.- DISPOSITIVO EXPERIMENTAL El dispositivo experimental que se utiliza para la realizacin de la prctica de humidificacin se muestra esquemticamente en la figura 8. LG TL2 TG2, 2 TL1 TG1, 1 Figura 8. Dispositivo experimental El dispositivo est formado por una columna y por dos circuitos, de agua y de aire. - La corriente de agua proviene de la red. Pasa por un calentador, despus por un rotmetro, y se alimenta en la columna por la parte superior. - La corriente de aire proviene de la lnea general de la facultad (aire atmosfrico que pasa por un compresor). El aire pasa por un rotmetro, y se alimenta en la columna por la parte inferior, pasando a su travs de forma ascendente. Regulamos ambos caudales mediante las vlvulas que preceden a los respectivos rotmetros, gracias a los cuales medimos dichos caudales.

El circuito de agua tiene un medidor de temperatura a la entrada y otro a la salida de la columna; el de aire tiene un medidor de temperatura y de humedad relativa a la entrada de la columna y otro a la salida. - La columna est formada por dos tubos cilndricos concntricos. El agua se desliza por la superficie externa del tubo interior en forma de pelcula descendente. El aire circula de forma ascendente por el espacio comprendido entre los dos tubos. Las dimensiones de la columna son las siguientes: Dimetro del tubo interior: di = 6 mm Dimetro del tubo exterior: de = 3.5 cm Altura de la columna: z = 36 cm 4.- PROCEDIMIENTO EXPERIMENTAL Haremos dos experimentos diferentes: uno en el que mantenemos el caudal de agua constante y variamos el de gas; y otro en el que mantenemos el caudal de gas constante y variamos el de agua. 1.- Abrimos las vlvulas V1 y V2 para poner en marcha nuestro dispositivo. 2.- Conectamos el calentador, ponemos el termostato en la posicin 3, y cerramos la vlvula V2. Dejamos que se caliente el agua alrededor de los 40C (lo comprobamos abriendo de vez en cuando la vlvula V2 y mirando la temperatura en el medidor de la entrada de agua en la columna.) El volumen de agua que cabe en el calentador no es elevado, estamos alimentndolo continuamente con agua fra de la red, y sale de manera continua agua caliente. Al proceder de esta forma hacemos que el agua alcance rpidamente la temperatura deseada, y que sta no vare demasiado a lo largo de la experimentacin. 3.- Seleccionamos un caudal de aire (entre 500 y 1500 L/h) y uno de agua (entre 40 y 100 L/h) mediante las vlvulas V1 y V2. 4.- Controlamos el experimento de forma que los caudales se mantengan constantes, y la pelcula de agua que se forma no se rompa. Anotamos entonces las lecturas de los respectivos rotmetros. 5.- Anotamos las temperaturas de entrada y salida de agua y aire, as como las humedades de entrada y salida de aire cuando observemos que estos datos permanezcan constantes con el tiempo , durante unos 10-15 segundos. Suponemos as que durante este periodo de tiempo trabajamos en estado estacionario. 6.- Manteniendo el caudal de agua constante, y variando el de aire desde 500 a 1500 L/h, realizamos 5 experimentos, repitiendo los puntos 4 y 5. 7.- Manteniendo el caudal de aire constante, y variando el de agua desde 40 a 75 L/h, realizamos 5 experimentos, repitiendo los puntos 4 y 5. A lo largo de la experimentacin, hemos de vigilar que la temperatura del agua no baje de 35C. Si es as, repetimos el punto 2. Tambin hay que repetir el punto 2 al pasar de un tipo de experimento a otro, es decir, al pasar del punto 6 al 7. 8.- Una vez finalizada la prctica, cerramos las vlvulas V1 y V2. 5.- RESULTADOS Y DISCUSIN En nuestra prctica se han realizado dos tipos de experimentos:

Experimento A: caudal de agua constante e igual a 55 L/h, variando el de aire desde 725 a 1450 L/h, con una temperatura del agua a la entrada de la columna de 38 - 40C.

Experimento B: caudal de aire constante e igual a 550 L/h, variando el de agua desde 50 a 72,5 L/h, con una temperatura del agua a la entrada de la columna de 40 - 47C.

Consideramos que nuestro sistema es adiabtico, y que trabaja en rgimen estacionario. Los datos experimentales obtenidos en ambos experimentos se recogen en las tablas I y II respectivamente: Tabla I: Evolucin de las condiciones de salida de una torre de humidificacin al variar el caudal de aire (experimento A) Exp. 1 2 3 4 5 qG (L/h) 725 875 960 1390 1450 TG1(C) 24,6 24,4 24,6 24,6 24,5 TG2(C) 24,7 24,5 23,8 22,6 21,9 ( 7,3 7,2 7,3 7,3 7,2 ( 49,1 50,9 53,4 58,6 60,8 qL(L/h) 55 55 55 55 55 TL2(C) 39,9 39,1 39,1 38,4 37,9 TL1(C) 39,3 38,6 38,1 37,7 37,2

Tabla II: Evolucin de las condiciones de salida de una torre de humidificacin al variar el caudal de agua (experimento B) Exp. 1 2 3 4 5 qG (L/h) 550 550 550 550 550 TG1(C) 24,7 24,7 24,7 24,6 24,7 TG2(C) 28,0 26,5 25,9 23,7 26,0 ( 7,3 7,1 7,3 7,2 7,2 ( 13,8 22,1 26,9 38,1 26,6 qL(L/h) 50,0 57,0 62,5 70,0 72,5 TL2(C) 47,3 47,5 41,5 41,3 40,7 TL1(C) 46,6 46,6 40,9 40,6 40,2

Observamos que la humedad del aire aumenta a la salida de la torre (columna), y que la temperatura del agua disminuye, para cada uno de los experimentos. Ntese que los subndices 1 y 2 indican en el pie de la columna y en la cabeza de la columna respectivamente, tal y como se indica en la figura 8 que representa el dispositivo experimental. La figura 9 no es ms que la representacin de la columna aislada de nuestro dispositivo experimental; la finalidad de esta representacin es clarificar el significado de los subndices, que pueden resultar confusos. Figura 9. Corrientes de entrada y salida en la torre de humidificacin experimental. qG y qL son los caudales volumtricos de aire y de agua. Para poder obtener los coeficientes de transferencia y el NUT, objetivo de la prctica, hemos de transformar nuestros datos experimentales: - Transformamos qG y qL en caudales msicos considerando que las densidades del aire y del agua son constantes para el rango de temperaturas de trabajo e iguales a: aire= 1,293 kg/m3 agua= 1000 kg/m3

- Transformamos la humedad relativa en absoluta: mediante la ecuacin 3 obtenemos la Pv a la temperatura del gas correspondiente, introducimos este valor en la ecuacin 2, y obtenemos la PA y con ella obtenemos la H humedad absoluta aplicando la ecuacin 1, donde P=760 mmHg.

- Calculamos la entalpa del gas a la entrada y la salida de la columna mediante las ecuaciones 4 y 5. Temperatura terica del agua a la salida Con las condiciones de entrada y salida de la corriente de aire, y las de entrada de agua, calculamos mediante el balance de entalpa, ecuacin 10, la temperatura terica del agua a la salida TL1'. Las TL1' calculadas se recogen en la tabla III: Tabla III: Aplicacin del balance de entalpa para la obtencin de TL1' en los experimentos A y B. Experimento A Exp. G (kg/h) 1 2 3 4 5 0,937 1,131 1,241 1,797 1,875 iG1 [J/kg] 28257 27965 28257 28257 28086 iG2 [J/kg] 49046 49440 48873 48168 47304 L (kg/h) 55 55 55 55 55 TL2(C) 39,9 39,1 39,1 38,4 37,9 T'L1(C) 39,82 38,99 38,99 38,24 37,74 TL1(C) 39,3 38,6 38,1 37,7 37,2

Experimento B Exp. G (kg/h) 1 2 3 4 5 0,711 0,711 0,711 0,711 0,711 iG1 [J/kg] 28379 28281 28379 28208 28330 iG2 [J/kg] 36353 38691 40207 41434 40233 L (kg/h) 50,0 57,0 62,5 70,0 72,5 TL2(C) 47,3 47,5 41,5 41,3 40,7 T'L1(C) 47,27 47,47 41,47 41,27 40,67 TL1(C) 46,6 46,6 40,9 40,6 40,2

Si comparamos los valores de la temperatura terica del agua a la salida con la experimental vemos que su diferencia oscila entre 0,4 y 0,9C. Dicha diferencia, que, considerando la temperatura a la entrada, es relativamente grande, es debida a que, por un lado estamos suponiendo estado estacionario, pero ste se alcanza en un corto perodo de tiempo, por lo que es posible que en el tiempo que se tarda en anotar todos los datos, alguno de ellos cambie, luego las condiciones del sistema estn cambiando y ya no estamos en estado estacionario. Hay que tener en cuenta que la temperatura de entrada del agua est cambiando de temperatura constantemente, a causa de lo comentado en el apartado 2 del procedimiento experimental. Por otro lado, estamos considerando que nuestro sistema es completamente adiabtico. El diseo de la columna hace que esta suposicin sea correcta en ella, dado que el agua circula por la superficie externa del tubo interior, estando sta en contacto slo con la corriente de aire, por lo que slo intercambiar calor con dicha corriente. A su vez, el aire podra intercambiar calor con el ambiente a travs de la pared de la columna, pero no lo hace, ya que el aire se encuentra a temperatura ambiente. Pero fuera de la columna, el agua caliente est circulando por una tubera que est en contacto con el ambiente. Los medidores de temperatura no estn colocados justo a la entrada y la salida de la columna, sino que estn a unos centmetros de ellas. As que en esta pequea distancia de tubera se producen pequeas prdidas de calor debidas a la diferencia de temperaturas entre la corriente de agua y el ambiente. Clculo del NUT: resolucin numrica de la integral

Para determinar los coeficientes de transferencia hemos de determinar previamente el NUT, definido en la ecuacin

de diseo (ec. 18). Por lo tanto hemos de resolver la integral , en la que desconocemos las condiciones de la interfase. Lo hacemos mediante el mtodo Mickley, descrito en el apartado 1.3.3, pero no de forma grfica, sino numrica, aplicando las ecuaciones. Puesto que desconocemos el valor de hL/kG , tendremos que suponerlo, hasta que la prediccin de temperatura de salida del gas coincida con el dato experimental. Para cada uno de los experimentos, tendremos que seguir el siguiente esquema de clculo iterativo: Tabla IV: Mtodo numrico para la resolucin del (NUT)G TL(C) TL2 TL2 + 1(TL1-TL2) /10 ... ... ... TL2 + 5(TL1-TL2) /10 ... ... ... TL2 + 9(TL1-TL2) /10 TL1 ` iG [J/kg] iG2 iGa,1 ... ... ... iGa,5 ... ... ... iGa,9 iG1 Ti(C) Ti2 Tia1 ... ... ... Tia5 ... ... ... Tia9 Ti1 iGi (J/kg) ii2 iia1 ... ... ... iia5 ... ... ... iia9 ii1 TGa5 ... ... ... TGa9 TG1 comprobacin TG (C) TG2 . TGa1 1/(iiiG) f2 fa1 ... ... ... fa5 fa9 f1

*Suponemos un valor de *Por otro lado suponemos Ti ii Tengo Ti, ii para esa TL Hacemos esto para cada TL. (TG2)calculada = (TG2)experimental Tenemos a) Dividimos la diferencia entre la temperatura de entrada y de salida del agua en 10 intervalos iguales. La temperatura de salida del agua que tomamos es la calculada mediante el balance de entalpa: TL1'. b) Para cada intervalo, calculamos la entalpa del gas mediante la recta de operacin, ecuacin 20. c) Suponemos un valor de , que es la pendiente de la recta de unin (ec.21). d) Por otro lado, para la primera TL (primera fila de la tabla), suponemos un valor de la temperatura en la interfase Ti. Con l, calculamos la entalpa de interfase ii aplicando la ecuacin 6.

e) Comprobamos que se cumpla la recta de unin (en la tabla, ecuacin de comprobacin). Si es as, ya tenemos las condiciones de interfase para nuestro intervalo (nuestra TL). Si no se cumple, hemos de volver a suponer otro Ti . Hemos de hacer esto para cada TL. f) Cuando tengamos las condiciones de interfase para todos los intervalos, calculamos las temperaturas del gas para toda la columna con la ecuacin 22, empezando por TG1 hasta TG2, es decir, en la tabla, de abajo a arriba. Entonces comparamos el valor de la temperatura de salida TG2 calculado con el experimental: si coinciden, ya tenemos nuestra pendiente de las rectas de unin ; si no coinciden, tendremos que suponer otro valor de este cociente, repitiendo los pasos de c) a f). g) Una vez obtenido , calculamos el valor de la funcin de dentro de la integral del NUT para todos los intervalos, e integramos la funcin por Simpson. La aplicacin de este mtodo implica un nmero muy elevado de iteraciones. Para agilizar los clculos, en la hoja de clculo, utilizamos una Macro que realiza rpidamente las iteraciones de los puntos d) y e) mediante Solver, para distintos valores de que nosotros introducimos, a la vez que nos da las TG2's para cada pendiente, con las que determinamos qu valor de la pendiente de la recta de unin es el correcto. La Macro ha sido programada por David Curc, profesor del departamento de Ingeniera Qumica. La finalidad de este mtodo es la resolucin de una integral, por lo tanto, se supone que a mayor nmero de intervalos, elementos diferenciales, mejor se describir la funcin, luego ms exacto ser el valor de la integral obtenido. Pero consideramos que con 10 intervalos es suficiente, porque al realizar el mismo mtodo con 30 intervalos las diferencias que se obtienen no son significativas. Adems, al ejecutar la Macro, con un nmero elevado de intervalos, el ordenador se colapsa. Tras la aplicacin del mtodo Mickley para cada experimento, obtenemos los siguientes resultados, representados en la tabla V: Tabla V: Valores de la pendiente de la recta de unin y del nmero de unidades de transferencia para cada experimento. Exp.A: Caudal de agua constante e igual a 55 kg/h Exp. 1 2 3 4 5 G (kg/h) 0,937 1,131 1,241 1,797 1,875 hL/kG [J/(kgC)] 2443 2474 1838 1333 1152 (NUT)G 0,566 0,601 0,699 0,869 0,948 Exp.B: Caudal de aire constante e igual a 0.711 kg/h L (kg/h) 50,0 57,0 62,5 70,0 72,5 hL/kG [J/(kgC)] ----13907 1573 51470 (NUT)G ----0,127 0,439 0,100

En los experimentos B.1 y B.2 no ha sido posible encontrar un valor de hL/kG para el que la temperatura de salida del gas calculada coincida con la experimental (al ir probando con valores de hL/kG nunca se alcanza la TG2exp , a mayor hL/kG mayor TG2calc , pero para determinado valor del cociente de coeficientes, la TG2calc alcanza un mximo, y a partir de aqu disminuye. Este mximo est por debajo de la temperatura experimental). Este hecho es debido a que al tomar las lecturas experimentales probablemente no se haba alcanzado el estado estacionario, por lo que el modelo matemtico descrito a lo largo de este informe no se cumple para estos puntos. Coeficientes de transferencia hL y kG Una vez hallados el NUT y el cociente , podemos determinar los valores de los coeficientes mediante la ecuacin de diseo: (18)

en la que ahora nuestra nica incgnita ser kG. S y a son la seccin transversal de la columna y el rea especfica del volumen de la columna respectivamente, iguales a: As, despejamos kG de la ecuacin 18, y con el cociente obtenemos hL. Recogemos los resultados obtenidos en la tabla VI: Tabla VI: Evolucin de los coeficientes de transferencia de materia de la fase gas y de calor por conveccin de la fase lquida al aumentar los caudales de gas y agua respectivamente. Exp.A: Caudal de agua constante e igual a 55 kg/h Exp. 1 2 3 4 5 G (kg/h) 0,937 1,131 1,241 1,797 1,875 kG (kg/h.m2) 78 100 128 230 262 hL (kJ/C.h.m2.) 191 248 235 307 302 Exp.B: Caudal de aire constante e igual a 0.711 kg/h L (kg/h) 50,0 57,0 62,5 70,0 72,5 kG (kg/h.m2) ----13,3 46,0 10,5 hL (kJ/C.h.m2.) ----185 72 539

Para observar mejor la evolucin de los coeficientes, hacemos las siguientes representaciones grficas:

Figura 10. Evolucin de kG en el experimento A.

Figura 11. Evolucin de hL en el experimento B Miremos cada experimento por separado: - Experimento A: como vemos en la figura 10, el coeficiente de transferencia de materia de la fase gas aumenta a medida que lo hace el caudal de aire. Es decir, que mientras mayor sea el caudal de aire, ms se humidificar. El coeficiente de transferencia de materia es inversamente proporcional al espesor de la pelcula laminar que se forma al poner en contacto dos fases que circulan advectivamente. Al aumentar la densidad de flujo del gas, estamos disminuyendo este espesor, por lo que el coeficiente de transferencia de materia aumenta. Respecto al coeficiente de transferencia de calor de la fase lquida hL debera permanecer constante, ya que depende de la geometra del sistema, las propiedades fsicas de la sustancia, y su caudal, que no varan a lo largo del experimento. Respecto al hG , tambin ha de aumentar con el caudal de gas. Si consideramos que para el sistema aire agua se

cumple que , su determinacin es inmediata; seguir la misma tendencia que el coeficiente de transferencia de materia. - Experimento B: como vemos en la figura 11, slo tenemos tres puntos experimentales que no siguen una tendencia clara. En teora, puesto que el caudal de aire permanece constante a lo largo del experimento:

hL debe aumentar a medida que lo hace el caudal de lquido (cosa que favorece el intercambio de calor, luego favorece al descenso de temperatura del lquido) hG debe permanecer constante, kG tambin debe hacerlo.

De los tres puntos obtenidos, consideramos que el del experimento B.4 no es vlido porque el valor de hL no sigue una tendencia ascendente, adems de que el valor de kG (tabla VI) se aleja mucho de los valores de las otras kG's. As que nuestra experimentacin queda reducida a la obtencin de nicamente dos puntos, que indican el aumento de hL a medida que lo hace el caudal de agua, como tendra que ser en teora. Pero dos puntos experimentales son insuficientes para extraer conclusiones. En este caso lo aconsejable sera repetir la experimentacin. El hecho de que en el experimento A se hayan obtenido buenos resultados y en el experimento B no puede ser debido a que en el B estamos variando el caudal de agua, que desciende por la columna en forma de pelcula, y

debido a esto el espesor de la masa del lquido es muy pequeo. Al variar su caudal, estamos desestabilizando ms el sistema que cuando varibamos slo el de aire, cuya masa tiene un espesor en comparacin mucho mayor. As, la variacin del caudal de agua produce un cambio brusco en las condiciones de la interfase. Por lo tanto, se necesita ms tiempo para alcanzar el estado estacionario. Al realizar la prctica no tuvimos en cuenta este hecho, y por ello slo nos conducen a un resultado lgico, en principio, 2 de 5 experimentos realizados. Por ltimo, el dispositivo experimental utilizado en esta prctica es adecuado para los objetivos que se desean alcanzar; no obstante el calentador no mantiene la temperatura de entrada a la columna del lquido constante, haciendo variar continuamente el estado del sistema. Aconsejamos la instalacin de otro calentador con, por ejemplo, mayor capacidad, para que no se produzca este problema; alternativamente, si se instalara otro calentador de las mismas caractersticas en serie con el que ya existe, se reducira mucho el cambio de temperatura a lo largo de la experimentacin. Tambin aconsejamos aislar las conducciones para evitar posibles prdidas de calor en la parte que va desde el medidor de temperatura hasta la columna, y viceversa. 6. CONCLUSIONES Los resultados obtenidos en este estudio nos conducen a las siguientes conclusiones:

en una torre de humidificacin o de enfriamiento de agua, al aumentar el caudal de gas aumentamos el coeficiente de transferencia de materia, por lo que estamos favoreciendo la humidificacin del gas. al aumentar el caudal de lquido, el coeficiente de transferencia de calor por conveccin de la fase lquida aumenta. Entonces el lquido pierde calor con mayor facilidad, luego estamos favoreciendo el enfriamiento del lquido.

7.- NOMENCLATURA a rea interfacial del volumen de la columna A superficie interfacial cP calor especfico cS calor especfico hmedo Deq dimetro equivalente de la columna G caudal msico del aire H humedad absoluta del aire ( hG coeficiente de transporte de calor por conveccin de la fase gaseosa hL coeficiente de transporte de calor por conveccin de la fase lquida HUT altura de la unidad de transferencia iG entalpa especifica del gas kG coeficiente de transporte de materia de la fase gaseosa L caudal msico del agua M masa molecular

NUT nmero de unidades de transferencia P presin total PA presin parcial PV presin de vapor q calor transferido a la columna S seccin transversal de la columna T temperatura z altura de la columna humedad relativa del aire densidad 0 calor latente de vaporizacin del agua

L entalpa de la fase lquida Subndices: A: agua B: aire G: gas i: interfase L: lquido s: condiciones de saturacin adiabtica 8.- BIBLIOGRAFA

qL TL1 qL

Treybal, R.E. (1980) Operaciones de transferencia de masa Ed. Mc. Graw-Hill. Costa, Cervera, Cunill, Espulgas, Mans, Mata (1999) Curso de ingeniera qumica Ed. Revert. Foust, Wenzel, Clump, Maus, Andersen (1979) Principios de operaciones unitarias Ed. C.E.C.S.A-

TL2 qG TG1 qG TG2

L1, TL1 L2, TL2 Figura 1: Operacin de saturacin adiabtica TS HS,TS HE,TE G, H2, TG2 G, H1, TG1 G H TG G H + dH TG + dTG L + dL TL + dTL

L TL hL [J/(m2sC)] hG [J/(m2sC)] kG [kg/(m2s)] Figura 7. Construccin grfica para determinar el punto geomtrico de la fase gaseosa

Introduccin En las operaciones bsicas controladas por la transferencia de materia y de calor se da simultneamente estas dos; hay, por tanto, dos fuerzas impulsoras separadas que pueden ir en el mismo sentido o en sentidos opuestos. Las operaciones bsicas ms importantes son: el acondicionamiento de gases o enfriamiento de lquidos (fases presentes lquido/gas); y el secado (fases presentes slidos o lquidos/gas). El acondicionamiento de gases es una operacin que tiene por objeto modificar las condiciones de humedad de una corriente de aire por interaccin con una corriente de agua; se conoce como humidificacin o deshumidificacin, se puede aplicar a cualquier sistema gas/lquido aunque la aplicacin ms extendida es al sistema aire/agua. Industrialmente tiene gran importancia el enfriamiento de una corriente de agua por interaccin con aire. La humidificacin como una operacin unitaria en la cual se da una transferencia simultnea de materia y calor sin la presencia de una fuente de calor externa, tiene lugar cuando un gas se pone en contacto con un lquido puro, en el cual es prcticamente insoluble. Este fenmeno nos conduce a diferentes aplicaciones adems de la humidificacin del gas, como son su deshumidificacin, el enfriamiento del gas (acondicionamiento de gases), el enfriamiento del lquido, adems de permitir la medicin del contenido de vapor en el gas. Generalmente la fase lquida es el agua, y la fase gas el aire. Su principal aplicacin industrial es el enfriamiento de agua de refrigeracin. En la deshumidificacin, agua fra se pone en contacto con aire hmedo. La materia transferida entre las fases es la sustancia que forma la fase lquida, que dependiendo de cmo estemos operando, o se evapora (humidificacin), o bien se condensa (deshumidificacin.) Existen diferentes equipos de humidificacin, entre los que destacamos las torres de enfriamiento por su mayor aplicabilidad. En ellas, el agua suele introducirse por la parte superior en forma de lluvia provocada, y el aire fluye en forma ascendente, de forma natural o forzada. En el interior de la torre se utilizan rellenos de diversos tipos que favorecen el contacto entre las dos fases.

Definicin de humidificacin y deshumidificacin

Humidificacin

La humidificacin es una operacin que consiste en aumentar la cantidad de vapor presente en una corriente gaseosa; el vapor puede aumentar pasando el gas a travs de un lquido que se evapora en el gas. Esta transferencia hacia el interior de la corriente gaseosa tiene lugar por difusin y en la interfase hay, simultneamente, transferencia de calor y de materia. A grandes rasgos, el proceso que tiene lugar en la operacin de humidificacin es el siguiente:

Una corriente de agua caliente se pone en contacto con una de aire seco (o con bajo contenido en humedad), normalmente aire atmosfrico. Parte del agua se evapora, enfrindose as la interfase. El seno del lquido cede entonces calor a la interfase, y por lo tanto se enfra. A su vez, el agua evaporada en la interfase se transfiere al aire, por lo que se humidifica. Deshumidificacin

La deshumidificacin es una operacin que consiste en reducir la cantidad de vapor presente en una corriente gaseosa, mediante una condensacin parcial del vapor, que se separa.

Aplicacin de la humidificacin/deshumidificacin Aplicacin general

La principal aplicacin de la humidificacin y de la deshumidificacin es en el acondicionamiento de aire y en el secado de gases. Un aspecto de inters relacionado con esta operacin bsica es el enfriamiento de aguas despus de un proceso industrial, con el fin de poder ser utilizada nuevamente; el equipo utilizado tiene forma de torre en la que el agua caliente se introduce por la parte superior y fluye sobre un relleno en contracorriente con aire que entra por la parte inferior de la torre de enfriamiento.

Aplicacin en la industria de alimentos

El conocimiento de los procesos de humidificacin y deshumidificacin, as como sus clculos implicados en ella, sern tiles en el diseo y anlisis de diferentes sistemas de almacenamiento y procesado de alimentos. As mismo, resulta imprescindible conocer las propiedades de las mezclas aire vapor de agua en el diseo de sistemas tales como equipos de aire acondicionado para conservar alimentos frescos, secaderos de granos de cereal y torres de enfriamiento en plantas de procesado de alimentos.