ARTÍCULO CIENTÍFICO / RESEARCH ARTICLE

Evaluación energética de un colector solar de placa plana de doble cubierta*

Energy assessment of a double cover flat-plate solar collector

Adrián Enrique Ávila Gómez** Jorge Mario Mendoza Fandiño*** Julio Fernando Beltrán Sarmiento****
Universidad Pontificia Bolivariana Seccional Montería (Colombia)

* Apoyos y subvenciones: Basado en el proyecto de investigación "Diseño, montaje y modelación de un secador integrado solar biomasa", en el marco de proyectos financiados internamente por la Universidad Pontificia Bolivariana de Montería al grupo de investigación OPUREB, Uso Racional de Energía y Biomasa (22 enero de 2007 - 22 de mayo de 2009).

** MSc en Ingeniería Mecánica. Director Facultad de Ingeniería Mecánica, Universidad Pontificia Bolivariana, seccional Montería (Colombia). adrian.avila@upbmonteria.edu.co

*** Ingeniero Mecánico. Docente Facultad de Ingeniería Mecánica, Universidad Pontificia Bolivariana, seccional Montería (Colombia). jorge.mendoza@upbmonteria.edu.co
Correspondencia: Universidad Pontificía Bolivariana. Km 8 vía Cereté, Córdoba (Colombia).

**** Ingeniero Mecánico. Jóven Investigador, Universidad Pontificia Bolivariana, seccional Montería. jfbeltran20@gmail.com

Fecha de recepción: 22 de julio de 2009
Fecha de aceptación: 1 de marzo de 2010

Resumen

En este artículo se establece la incidencia de la relación de aspecto (razón entre la longitud y ancho de un colector de área constante) sobre la eficiencia térmica de un colector solar de placa plana que utiliza doble cubierta transparente para el calentamiento de aire. Este estudio se realiza desde dos enfoques: el primero utiliza un modelo físico-matemático a partir de los balances de energía del colector, y el segundo utiliza un prototipo con su respectiva instrumentación para calcular su rendimiento de manera experimental. El colector está compuesto por una estructura en madera, una doble cubierta de vidrio transparente, una placa absorbedora de radiación solar y un ventilador para extraer el aire calentado. Para calcular el rendimiento se registraron de manera automática datos de temperatura, radiación solar y flujo másico de aire; durante aproximadamente 6 horas por varios días y bajo distintas condiciones climáticas de la ciudad de Montería, Córdoba (Colombia).

En general se encontró, para un área constante del colector, que la eficiencia térmica de este aumenta cuando la relación de aspecto aumenta y que el modelo teórico se ajusta de manera aceptable con los resultados experimentales.

Palabras clave: Área constante; colector solar; eficiencia; placa plana; relación de aspecto.

Abstract

This study establishes the effect of the aspect ratio (ratio between length and width of a constant collector area) on the thermal efficiency of a flat-plate solar collector that uses a double cover for air heating. The analysis has two approaches: first, using a physical-mathematical model that develop a steady state analysis from a collector energy balance and the second, using an instrumented prototype to calculate its performance in experimental way. The collector configuration consists in a wood structure, a double glass cover, an absorber plate of solar radiation and blower to extract the heated air. To calculate the performance temperature data, solar radiation and mass flow of air were automatically recorded for about 6 hours for several days and under different weather conditions in the city of Montería, Cordoba.

In general it was found, for a constant collector area, where the thermal efficiency increases when the aspect ratio increases and that the theoretical model fits in an acceptable way with the experimental results.

Keywords: Constant collector area; solar collector; efficiency; flat-plate; aspect ratio.

1. INTRODUCCIÓN

El secado es el proceso comercial más utilizado para la conservación y preservación de la calidad de los productos agrícolas [1]. Esta operación requiere de la utilización de temperaturas relativamente bajas, por lo cual son ideales para el aprovechamiento de energías renovables y disponibles en elevadas cantidades, como lo es el uso de la energía solar [2], [3].

El departamento de Córdoba (Colombia) se encuentra ubicado en una zona privilegiada en cuanto a la incidencia de radiación solar (aprox. 5Kwh/m² como promedio anual) [4], por lo cual esta región dispone de un recurso energético valioso, que potencialmente se puede convertir en eje de desarrollo tecnológico agroindustrial y beneficiar económicamente a los agricultores de esta zona que se destaca por su gran vocación agrícola.

En los últimos años se han estudiado y desarrollado diversas técnicas y equipos con el propósito de utilizarlos en el secado de especies hortícolas, granos y demás productos agrícolas; dentro de estos estudios se destacan sistemas como los de Bistoni et al. [5], que consideran a los colectores de aire como óptimos para el uso de sistemas de secado de productos agrícolas que utilicen energía solar, y describen un modelo numérico que incluye el análisis de la transferencia de calor entre las distintas partes constitutivas del colector, las temperaturas y la eficiencia del mismo. También se destaca Costa et al. [6], quienes combinaron el secado mecánico con el natural al utilizar la energía colectada como radiación solar y convección forzada del aire en el procedimiento de secado de granos, para lo cual implementaron un colector de placa plana. Otros autores, como Belen et al. [7], quienes utilizan sistemas híbridos solar-eléctricos, diseñaron un prototipo de colector solar para establecer el efecto del secado en las propiedades bromatológicas de la palma coroba. En el desarrollo de los colectores solares para el calentamiento de aire se han estudiado diversos tipos como son los de placa corrugada en "V", de doble paso y los de placa plana [8]. Por lo general se utilizan en aplicaciones de temperatura moderada (hasta 100 °C por encima de la temperatura ambiente), emplean tanto la radiación solar directa como la difusa, no requieren seguimiento solar y necesitan poco mantenimiento. Estas características los hacen ideales para el diseño de sistemas solares de energía al mínimo costo. Bajo esta perspectiva se justifica el diseño de un colector de placa plana con una eficiencia baja, si los costos se reducen sustancialmente.

Los colectores de placa plana se han estudiado desde hace varios años, principalmente en las aplicaciones de calentamiento de agua. En cuanto al calentamiento de aire se destacan los estudios realizados por Ho et al. [9], en el cual se estudia el efecto de la utilización de bafles y doble paso de aire en el desempeño de un colector solar de placa plana. Gupta et al. [10] realizaron una evaluación energética de un colector solar de placa plana, en el cual se estableció la relación de aspecto y la profundidad del ducto óptimas del colector. Karwa et al. [11] evaluaron analíticamente el rendimiento termo-hidráulico de un calentador solar de aire, con base en la igualdad de poder de bombeo de colectores o subcolectores unidos en serie o en paralelo, variando la profundidad del ducto y evaluando la eficiencia energética de los módulos de subcolectores en serie o en paralelo para deducir la mejor configuración para el conjunto de parámetros investigados.

Los colectores de placa tienen una gran aplicabilidad en sectores de la sociedad donde el factor costo es más importante que el de eficiencia; sin embargo, se requiere que la implementación de estos sistemas de baja eficiencia aprovechen la mayor cantidad de energía posible. Con base en esto, el propósito de la investigación fue realizar un estudio teórico - experimental de la influencia que tiene la relación de aspecto (razón entre la longitud y ancho de un colector de área constante) sobre la eficiencia de un colector solar de placa plana que utiliza doble cubierta transparente para el calentamiento de aire.

2. DESCRIPCIÓN TEÓRICA

Colector solar para el calentamiento de aire

Su forma más simple consiste de una lámina de vidrio o material transparente situado en la parte superior, una placa plana que actúa como un cuerpo negro para absorber calor y una estructura aislada térmicamente que dé soporte a los elementos anteriores. Para su funcionamiento los rayos del sol pasan a través del vidrio y son atrapados en el espacio entre la cubierta y la placa o absorbedor. El calor es posteriormente usado para la calefacción de las casas, para el calentamiento de piscinas o para calentar el aire.

El uso de aire como fluido de trabajo hace que estos colectores ocupen un importante puesto entre los sistemas de calentamiento solar, puesto que reduce el número de componentes requeridos en estos.

Calor útil

De acuerdo con Duffle et. al [12], el calor útil para el aire se puede calcular mediante las siguientes ecuaciones:

Donde A_c, es el área del colector, S es la radiación solar absorbida por unidad de área, U es el coeficiente global de transferencia de calor, T_pm, es la temperatura media de la placa, T_a es la temperatura ambiente, T_g es la transmitancia del vidrio, a_p es la absorbancia de la placa, U_l es el coeficiente global de perdidas del colector, m es el flujo másico, C_p es el calor especifico, T_o es la temperatura de salida y T_i es la temperatura de entrada.

F_R es el factor de remoción de calor y es dado por:

Para la ejecución de ecuación anterior (4) es conveniente definir el factor de flujo del colector F" como la relación entre F_R y F'.

Este factor de flujo del colector está en función de la capacitancia adimen-sional, la cual está dada como:

La temperatura de placa media de la ecuación (1) está dada por

La temperatura del fluido media puede ser encontrada integrando la siguiente ecuación [12] para valores de 0 a L

Al resolver esta integral y sustituyendo FR y Qu de las ecuaciones (4) y (3) respectivamente, la temperatura de fluido media demostrada por Klein et al. [13] es:

Cálculo de coeficiente de pérdidas de calor total

A lo largo de la dirección de flujo existen puntos en los cuales la energía solar absorbida calienta la placa a una temperatura TP. Además, se transfiere energía de la placa al aire ambiente por medio del coeficiente de pérdida trasero Ub, hacia el fluido Tf a través del coeficiente de transferencia de calor por convección h2, y a la parte inferior de la cubierta de vidrio debido al coeficiente de transferencia de calor por radiación linealizado hr. La energía es transferida a la cubierta de vidrio desde el fluido por medio del coeficiente de transferencia de calor por convección hl. Se pierde energía hacia el ambiente a través del efecto combinado de los coeficientes de convección y radiación Ut. Es importante anotar que el Ut puede incluir los efectos de la utilización de múltiples cubiertas.

A partir del balance de energía en la cubierta, la placa y el del fluido se obtienen las siguientes ecuaciones:

Estas tres ecuaciones se resuelven de manera que el calor útil esté expresado como una función de U_t, h₁, h_r, T_f y T_a, Es decir, se eliminan, T_p y T_c; luego, al sustituir esto dentro de la ecuación q_u, se obtiene que:

Por consiguiente, el coeficiente de pérdidas de calor total y el factor de eficiencia para colectores solares de placa plana se muestra en la figura. 2

Con h =h₁ =h₂ están dados por las siguientes ecuaciones [14]:

El coeficiente de pérdidas de calor inferior está dada por la ecuación:

Y el coeficiente de pérdidas de calor superior por la ecuación

Donde e = 0,430(1-100/T_pm)

f = (1 + 0,089h_W - 0,1166h_wε_p)(1 + 0,07866β²) C = 520(1 - 0,000051β²) para 0^o < β < 70 ° l = 0,00591Nh_w

Y el coeficiente de transferencia de calor por convección del aire que fluye por encima de la superficie exterior de la cubierta de vidrio, el cual depende de la velocidad del viento [15] está dado por la ecuación

El coeficiente de transferencia de calor por radiación h_r entre la placa ab-sorbedora y la cubierta, está dada por la siguiente ecuación:

Muchas aplicaciones de ingeniería implican transporte por convección en ductos no circulares. Sin embargo, muchos de los resultados del ducto circular se vuelven aplicables mediante el uso de un diámetro efectivo como longitud característica [16]. A este se le denomina diámetro hidráulico y se define como

Este es el diámetro que se debe utilizar para calcular parámetros como R_e y N_u.

Donde el flujo másico es calculado con la ecuación

Cálculo del coeficiente convectivo

De acuerdo con Kays [17], para régimen de flujo turbulento completamente desarrollado (D_h > 30) para la condición en la cual un lado es calentado y el otro permanece aislado, se tiene la siguiente correlación:

Cuando se tienen ductos cortos, el desarrollo térmico e hidrodinámico resaltará en un aumento significativo en el coeficiente de transmisión de calor cerca de la entrada. Para este caso, Rohsenow et al. [18] presenta números de Nusselt locales para el caso de temperatura de pared constante. La correlación correspondiente en este caso es:

Donde las constantes a, b, m y n, se encuentran en la tabla 1.

De este modo se tiene que los coeficientes de transferencia de calor por convección dentro del ducto, h₁ y h₂ equivalen a:

Cálculo de la eficiencia del colector

El rendimiento del colector está determinado por su eficiencia, la cual se define como la relación entre el calor útil en períodos de tiempo específicos sobre la energía solar incidente en un área determinada en el mismo período [12]. A partir de esta definición, la representación conveniente de la eficiencia para colector solar de lámina y tubo es:

3. MATERIALES Y MÉTODOS Evaluación teórica del colector

El procedimiento para el cálculo de los valores teóricos de la eficiencia r se describe a continuación. Primero, con geometrías del colector conocidas (Longitud, L=1,3; 1,8385; 2,2517 y 2,6m; Ancho, B=0,65; 0,4596; 0,3753 y 0,325m; Altura, H=0,09m y propiedades del sistema; τ_g, α, C_p, ρ(densidad), m (viscosidad absoluta), k(conductividad térmica), κ_b, δ_b (espesor), Є_p, Є_c, Є_g, (emitancia de la placa, cubierta y vidrio), así como las condiciones de operación dadas; Io (radiación solar incidente), Ta, V_w (Velocidad del viento), m, T_i, se estiman los valores de las temperaturas T_pm, T_fm y se calcula un valor temporal de r por medio de la ecuación (26). Los valores de T_pm y T_fm son entonces comprobados con las ecuaciones (8) y (9), respectivamente; en caso de que esto valores no coincidan se itera nuevamente hasta que los valores asumidos sean iguales a los calculados. Con base en estos resultados, el valor correspondiente de r es finalmente obtenido.

Evaluación experimental del colector

El colector solar de 0,845 m² de área está constituido por cuatro elementos fundamentales: una placa absorbedora, una doble cubierta transparente, un ventilador y una carcasa aislada que encierra y soporta cada uno de los elementos ya mencionados. Para la construcción del sistema se ha utilizado lámina de acero galvanizado pintada de negro mate, placas de vidrio de 4mm de espesor, láminas de icopor como material aislante y madera. Todos estos materiales son de bajo costo y de fácil consecución.

Los equipos utilizados en la toma de datos fueron: termómetros de mercurio, termocuplas tipo K, un piranómetro marca LI-COR 200, un anemómetro y una tarjeta de adquisición de datos. Las mediciones fueron tomadas en la ciudad de Montería - Córdoba, localizada al noroeste de Colombia a 08° 45' 27" de latitud Norte y 75° 53' 24" de longitud Oeste. Para calcular el rendimiento se registraron de manera automática datos de temperaturas, radiación solar global y flujo másico de aire; durante aproximadamente 6 horas por varios días y bajo distintas condiciones climáticas. Para esto se distribuyeron uniformemente seis termocuplas a lo largo y ancho de la placa. A la salida del colector se ubicó un anemómetro que permitía el registro de la temperatura y velocidad de salida del aire. El piranómetro se ubicó en la parte central del colector para establecer la radiación global incidente. Todos estos elementos estaban conectados a una tarjeta de adquisición de datos marca DaqPRO de ocho canales, la cual guardaba los datos para ser procesados en un computador posteriormente.

Las condiciones experimentales y las propiedades físicas empleadas en este estudio son: A_c=0,845m2; L=1,3m; B=0,65m; H=0,09m; t_g=0,77; a=0,96; Є_c=0,90; Є_p=0,84; Є_g=0,0.4; Io (Ver Figuras. 4, 5, 6 y 7); Ta=303K; V_w=0,74m/s; m=0,012Kg/s. 210000. Los valores de emisividad, transmitancia y absorbancia se tomaron con base en la configuración del colector y los valores suministrados por [19] y [20]. Se asume que la emisividad espectral de las superficies y la absorbancia no dependen de la longitud de onda, es decir que se comportan como un cuerpo gris; así mismo, que la absorbancia de la cubierta es despreciable con respecto a las pérdidas que por esta se pueda generar [12].

4. RESULTADOS Y DISCUSIONES Evaluación teórica

Se estudió el efecto que tiene la relación de aspecto sobre la eficiencia de un colector solar de placa plana. Para esto se estableció el desempeño del colector bajo regímenes de flujo turbulento y laminar, utilizando distintas relaciones de aspecto (2:1, 4:1, 6:1 y 8:1). En la figura 2 se muestra el efecto que tuvo la variación de la relación de aspecto cuando se evaluó el colector bajo régimen turbulento. En esta figura se ve claramente que la eficiencia del colector se ve fuertemente afectada por la relación de aspecto y el flujo másico de aire, de modo que a medida que aumenta la relación de aspecto la eficiencia aumenta, incluso para un flujo másico fijo. Por otro lado, a medida que el flujo másico incrementa la eficiencia del colector también se aumenta.

También se estudió el efecto que tiene la relación de aspecto sobre la eficiencia de un colector de placa plana cuando se encuentra bajo régimen laminar utilizando las relaciones de aspecto establecidas anteriormente. En la figura 3 se muestra el efecto que tuvo la variación de la relación de aspecto cuando se evaluó el colector bajo régimen laminar. En esta figura se ve que la eficiencia del colector no es afectada por la relación de aspecto en la misma proporción que para el flujo turbulento, excepto para el punto que corresponde a la mayor relación de aspecto y el máximo flujo másico. La razón para este comportamiento atípico de este punto se puede explicar en el hecho de que las condiciones de flujo y la geometría del colector en este punto hacen que se acerque mucho a la transición entre el régimen laminar y el turbulento. Por otro lado, al igual que en el análisis anterior, el flujo másico de aire afecta de modo que a medida que este aumenta, la eficiencia también aumenta.

Los resultados encontrados en el estudio teórico se pueden explicar en el hecho de que a medida que la relación de aspecto aumenta, el área transversal del ducto por donde fluye el aire se reduce, lo cual incrementa la velocidad de flujo del aire y también aumenta la tasa de transferencia de calor convectiva desde la superficie absorbedora de radiación hacia el aire.

Este fenómeno es más evidente para regímenes de flujo turbulento debido a que esta condición, a su vez, implica mayores tasas de transferencia de calor. El incremento de la eficiencia con el aumento de flujo másico se explica de manera análoga, ya que el régimen de flujo depende directamente del flujo másico, y tiende a turbulento cuando este se incrementa

Evaluación experimental

Una vez realizado el estudio teórico del colector, se procedió a realizar una evaluación experimental utilizando un prototipo, cuya geometría es equivalente la relación de aspecto 2:1 usada en la evaluación teórica. En este colector se estableció la variación de la eficiencia instantánea a distintas horas del día y para diferentes flujos másicos.

Estos resultados se compararon con los cálculos de eficiencia utilizando el modelo teórico antes descrito con los parámetros de radiación promedio registrados por el piranómetro a las distintas horas.

En las figuras 4, 5, 6 y 7 se muestra que la eficiencia experimental al inicio de la puesta en marcha del colector tiende a ser más baja que la teórica; sin embargo, a lo largo de la operación esta tiende a aumentar hasta el punto que al final de la toma de datos (2:00 PM) es superior a la teórica. Esta situación se da en los dos primeros flujos másicos evaluados, pero difiere del más alto en el cual la eficiencia siempre es más alta que en la teórica.

Las inconsistencias de estos resultados radican en la limitación inherente que tiene el hecho de comparar un análisis en estado estable con uno completamente dinámico, como es la realización de pruebas al aire libre, donde las condiciones meteorológicas son tan fluctuantes, especialmente si se tiene en cuenta que la época del año en la cual se hicieron las mediciones en la ciudad de Montería se caracteriza por los súbitos cambios de clima. Para sobrellevar estos inconvenientes se requiere realizar las mediciones en épocas del año más estables y en mayor cantidad, o realizar un análisis del colector en estado transitorio. En la segunda parte de este estudio se tendrán en cuenta estas consideraciones con el fin de desarrollar un análisis más adecuado para su comparación con el modelo teórico.

5. CONCLUSIONES

Al estudiar el desempeño del colector bajo régimen de flujo turbulento y laminar y utilizando distintas relaciones de aspecto se encontró que: La eficiencia del colector se ve más afectada cuando se encuentra bajo régimen de flujo turbulento que bajo régimen de flujo laminar, debido a que la primera condición implica mayores tasas de transferencia de calor. El incremento de la eficiencia con el aumento de flujo másico de aire también es más evidente para el flujo turbulento debido a que al aumentar la relación de aspecto, el área transversal del ducto por donde fluye el aire se reduce, lo cual incrementa la velocidad de flujo del aire y con ello el flujo másico, tendiendo a turbulento cuando este se incrementa.

Por otro lado, se tiene que al realizar la evaluación experimental preliminar utilizando un prototipo de colector solar equivalente a una de las relaciones de aspecto evaluadas anteriormente y al comprarlas con los cálculos teóricos, se encontró que: Las inconsistencias de los resultados radica en el hecho de comparar un análisis en estado estable con uno completamente dinámico, como es la realización de pruebas al aire libre donde las condiciones meteorológicas son bastantes fluctuantes. Para sobrellevar estos inconvenientes se requiere realizar mediciones en épocas del año más estables o realizar un análisis del colector en estado transitorio.

Referencias

[1] M. Karim, M. Hawlader, "Development of solar air collectors for drying applications", Energy Conversion and Management, vol. 45, n.° 3, 2004, p. 329.        [ Links ]

[2] S. Nandwani. "Energía solar. Conceptos básicos y su utilización". Universidad Nacional, Heredia (Costa Rica). Jun. 2005 [Online]. Disponible: http://www.una.ac.cr/fisica/index.php?option=com_content&task=view&id=15&Itemid=28 [Consultado: Mar. 20, 2009]        [ Links ]

[3] G.N. Tiwari. Solar energy. Fundamentals, design, modelling and applications. New Delhi: Nova Science Publishers, 2004, p. 94.        [ Links ]

[4] H. Rodríguez y F. González, Manual de radiación solar en Colombia, Bogotá: Universidad Nacional de Colombia, 1992.        [ Links ]

[5] S. Bistoni, A. Iriarte y L. Saravia. "Análisis teórico de un colector solar de aire", Revista Avances en energía renovables y medio ambiente, vol. 7, pp. 201-206, 2003.        [ Links ]

[6] A. Costa y S. Ferreira. "Sistema de secado solar para frutos tropicales", Información Tecnológica, vol. 18, n.° 4, pp, 49-58, 2007.        [ Links ]

[7] D. Belen et. al, "Efecto del secado solar en los contenidos de humedad, carbohidratos, carotenoi-des totales e índice de peróxidos del mesocarpio de la palma coroba", Interciencia, vol. 32, n.° 4, pp. 257-261, 2007.        [ Links ]

[8] A. Karim, and M. Hawlader, "Development of solar air collectors for drying applications", Energy Conversion and Management, vol. 45, 2004.        [ Links ]

[9] C. Ho, H. Yeh, T. Cheng, T.C. Chen and R. Wang, "The influences of recycle on performance of baffled double-pass flat-plate solar air heaters with internal fins attached", Applied Energy, vol. 86, n.° 9, pp. 1470-1478, 2009,        [ Links ]

[10] M. Gupta, S. Kaushik, "Performance evaluation of solar air heater having expanded metal mesh as artificial roughness on absorber plate", International Journal of Thermal Sciences, vol. 48, n.° 5, pp. 1007-1016, 2009.        [ Links ]

[11] R. Karwa R, S. Garg, A. Arya, "Thermo-hydraulic performance of a solar air heater with n-subcollectors in series and parallel configuration", Energy, vol. 27, n.° 9, pp. 807-812, 2002.        [ Links ]

[12] J. A. Duffie, and W. A. Beckman, Solar engineering of thermal processes, 3^rd ed. New York: Wiley, 1980, pp. 251-299.        [ Links ]

[13] S. Klein, J. A. Duffie, and W. A. Beckman,"Transient considerations of flat-plate solar collectors", J Eng Power Trans ASME, vol. 96, pp. 109-113 9, 1974.         [ Links ]

[14] J. A. Duffle, and W. A. Beckman, Solar engineering of thermal processes, 3^rd ed. New York: Wiley, 1980, pp. 282-285.        [ Links ]

[15] W. H. McAdams, Heat transmission, 3rd ed. New York: McGraw-Hill,1954.        [ Links ]

[16] B. Petukhov, T. F. Irvine, and J. P. Hartnett, Advances in Heat Transfer, vol. 6, New York: Academic Press, 1970.        [ Links ]

[17] W. Kays W. Convective heat and mass transfer. New York: McGraw-Hill, 1966.        [ Links ]

[18] W. M. Rohsenow, and H. Choi, Heat mass and momentum transfer. New Jersey: Prentice-Hall, 1961.        [ Links ]

[19] J. Kreider, C. Hoogendoorn, F. Kreith, Solar design: Components, systems, economic. New York: Hemisphere Publishing Corporation, 1989, pp. 44-55        [ Links ]

[20] ASHRAE, Handbook of HVAC Aplications, vol 4, Ch 32. Atlanta: ASHRAE, 1997.        [ Links ]