HERRAMIENTA COMPUTACIONAL PARA OPTIMIZAR EL DESEMPEÑO Y

USO DE UN EVAPORADOR CONCRETO EN AIRE ACONDICIONADO

YELMIS BARROS GALVIS

DIANA CAROLINA JAIMES BURGOS

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERÍAS FÍSICOMECÁNICAS
ESCUELA DE INGENIERÍA MECÁNICA
BUCARAMANGA
2010
HERRAMIENTA COMPUTACIONAL PARA OPTIMIZAR EL DESEMPEÑO Y
USO DE UN EVAPORADOR CONCRETO EN AIRE ACONDICIONADO

YELMIS BARROS GALVIS

DIANA CAROLINA JAIMES BURGOS

Trabajo de grado para optar al título de Ingeniero Mecánico

Director:
OMAR ARMANDO GELVEZ AROCHA
Ingeniero Mecánico

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERÍAS FÍSICOMECÁNICAS
ESCUELA DE INGENIERÍA MECÁNICA
BUCARAMANGA
2010
3
4
5
6
 DEDICATORIA

A mi niña Angie Vanessa Barros por esos momentos alegres y aquel entusiasmo
que me ha transmitido

A mi vieja Angela Galvis que lucha, lucha y lucha por sacar a sus hijos adelante;
gracias por ese gran apoyo y comprensión que me ha brindado.

A mi esposa Leisy De La Hoz por su compresión, cariño y respeto que me ha

tenido.

A mi hermanos Nelson Barros y Sauddy Barros por esa fuerza y apoyo que me
han dado.

A mi abuelita linda Dominga Benavides por esos consejos, compresión y amor que
siempre me ha tenido.

A todos aquellos que con sus conocimientos me ayudaron a realizar este sueño.

Yelmis Barros Galvis

7
 DEDICATORIA

A Papá Dios, infinitas gracias por todo lo que me ha dado, me da y me dará.

A mis padres Ana y Jesús, por darme cada día un pedacito de ellos y por
demostrarme cada día el amor incondicional que de Dios emana.
A mis hermanos Jesús y Jebel, por ser mis amigos y compañeros de vida.
A la memoria de los que ya no están con nosotros pero que su recuerdo
permanece intacto y viven en nosotros, nono Hernán y Hernando.
A mis amigos, profesores y compañeros de la U… caray! Muchas, muchas gracias
por su apoyo y lecciones de vida. Creo que hoy soy mucho más madura y mi
visión de mundo se amplió gracias a ustedes.
A mi compañero de proyecto Yelmis, su Esposa Leisy y su hija Angie por su
paciencia y apoyo; mi gratitud y respeto.
Al director de proyecto Omar Gelvéz, mi respeto y admiración por su
profesionalismo y su esfuerzo por dar lo mejor y exigir lo mejor.
A los que de algún modo nos brindaron asesoría o simplemente nos resolvieron
alguna inquietud para lograr este proyecto, gracias.

Diana Carolina Jaimes Burgos

8
 AGRADECIMIENTOS

Los autores expresan sus agradecimientos a:

El ingeniero Omar Armando Gélvez Arocha, director de proyecto, por su

dedicación, orientación y conocimientos aportados durante el desarrollo de este
trabajo.

A Don Nelson, Andrés Salinas, Héctor Navarro, Fabio Sangregorio, Luis Bernal y
Libardo por su valiosa colaboración en el desarrollo de este proyecto.

A todos los profesores por sus enseñanzas en el transcurso de nuestra carrera.

9
 TABLA DE CONTENIDO

PAG.
INTRODUCCIÓN 23

1. FUNDAMENTOS DEL AIRE ACONDICIONADO 26

1.1 AIRE HÚMEDO COMO GAS IDEAL 26
1.2 PROPIEDADES DEL AIRE ATMOSFÉRICO 27
1.2.1 Temperatura de punto de rocío 27
1.2.2 Temperatura de bulbo húmedo 27
1.2.3 Temperatura de bulbo seco 27
1.2.4 Humedad específica y relativa del aire 28
1.3 LA CARTA PSICROMÉTRICA 30
1.3.1 Ciclo de evolución del aire 32
1.3.2 Confort térmico y acondicionamiento de aire 33
1.3.3 Condiciones del local 34
1.3.4 Procesos de acondicionamiento de aire 35
1.3.4.1 Enfriamiento del aire húmedo. 36
1.3.4.2 Enfriamiento con deshumidificación del aire húmedo. 37
1.3.3.3 Mezcla adiabática de corrientes de aire o ventilación 39
1.4 INTERCAMBIADORES DE CALOR 40
1.4.1 Tipos de intercambiadores para aire acondicionado. 42
1.4.1.1 Evaporador enfriado por agua 42
1.4.1.2 Evaporador de expansión directa DX o enfriado por refrigerante. 43
1.4.2 Análisis de los intercambiadores de calor 45
1.4.2.1 Método de la LMTD. 45
1.4.2.2 El coeficiente global de la transferencia de calor. 47
1.4.3 Transferencia de calor de intercambiadores con aletas 47
1.4.4 Parámetros característicos de una aleta 48

10
1.4.4.1 Eficiencia de una aleta 48
1.4.4.2 Eficiencia superficial 49
1.4.4.3 Efectividad de una aleta 49
1.4.5 Coeficiente de transferencia de calor dentro de los tubos 51
1.4.5.1 Flujo turbulento del líquido dentro de los tubos 51
1.4.5.2 Coeficiente de transferencia de calor de evaporación del refrigerante. 52
1.4.6 Coeficiente de transferencia de calor externo en bancos de tubos 54
1.5 TRANSFERENCIA SIMULTÁNEA DE CALOR Y MASA 56
1.6 CAPACIDAD DE UN SISTEMA DE REFRIGERACIÓN 57
1.6.1 Efecto refrigerante 57
1.6.2 El subenfriamiento del refrigerante 58

2. MODELAMIENTO MATEMÁTICO PARA UNIDADES MANEJADORAS DE

AIRE 60
2.1 MODELAMIENTO MATEMATICO PARA UNIDAD MANEJADORA
ENFRIADO POR AGUA 61
2.1.1 Parámetros conocidos 61
2.1.2 Análisis del fenómeno de transferencia de masa y calor 62
2.1.2.1 Algoritmo global del modelamiento de la unidad manejadora enfriada
por agua 65
2.1.2.2 Diagrama de flujo del modelamiento para evaporador enfriado por
agua 66
2.2 MODELAMIENTO MATEMATICO PARA EVAPORADOR ENFRIADO
POR REFRIGERANTE 70
2.2.1 Parámetros conocidos. 70
2.2.2 Análisis del fenómeno de transferencia de masa y calor 71
2.2.2.1 Algoritmo global del modelamiento de la unidad manejadora enfriada
por refrigerante 73
2.2.2.2 Algoritmo del modelamiento para evaporador enfriado por refrigerante. 74

11
3. HERRAMIENTA COMPUTACIONAL PARA ESTUDIAR LAS UNIDADES
MANEJADORAS EN AIRE ACONDICIONADO 80
3.1 ANÁLISIS DE UNIDADES MANEJADORAS 82
3.2 SELECCIÓN DE UNIDADES. 83
3.3 COMPROBACION DE LA UNIDAD MANEJADORA 91
3.4 CALCULO PERSONALIZADO DE LA UNIDAD MANEJADORA 102

4. ANÁLISIS DE LOS RESULTADOS DE LA HERRAMIENTA

COMPUTACIONAL 105
4.1 CORROBORACIÓN DE LOS RESULTADOS DE LA HERRAMIENTA
CON OTRAS FUENTES 105
4.1.1 Comparación de resultados del programa Vs. Catálogo Modelo YGH
de York para el Agua 105
4.1.1.1 Comparación de la carga sensible entre el catálogo de la UMA para
agua de enfriamiento vs. la herramienta computacional 106
4.1.1.2 Comparación de la carga latente entre el catálogo de la UMA para
agua de enfriamiento vs. la herramienta computacional. 109
4.1.1.3 Comparación de temperatura de bulbo seco del aire de salida entre
el catálogo de la UMA para agua de enfriamiento vs. la herramienta
computacional. 112
4.1.1.4 Comparación de la humedad relativa del aire salida entre el catálogo
de la UMA para agua de enfriamiento vs. la herramienta computacional. 114
4.1.2 Mediciones tomadas de la UMA en el Edificio de la Facultad de Ciencias
Humanas- Primer piso Sala de Computo. 116
4.1.3 Comparación de resultados del programa Vs. Catálogo Modelo RMX
de York para el Refrigerante R-22 126
4.1.3.1 Comparación de la humedad relativa del aire salida entre el catálogo
de la UMA para vs. la herramienta computacional para refrigerante R-22 126

12
4.2 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA
HERRAMIENTA TENIENDO EN CUENTA VARIACIONES EN LOS
PARÁMETROS CONOCIDOS 131
4.2.1 Análisis del Efecto de las condiciones de entrada de los fluidos de
trabajo 131
4.3 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA
HERRAMIENTA TENIENDO EN CUENTA LA EVOLUCIÓN DEL AIRE AL
PASO POR LAS FILAS DEL SERPENTÍN 143
4.3.1 Análisis del Efecto de las condiciones de entrada de los fluidos de
trabajo 143
4.4 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA
HERRAMIENTA TENIENDO EN CUENTA LA VARIACION DEL NUMERO
DE FILAS 146

CONCLUSIONES 149

RECOMENDACIONES 151

REFERENCIAS BIBLIOGRAFICAS 152

ANEXOS 155

13
 LISTA DE FIGURAS

PAG.
Figura 1. Carta psicrométrica 31
Figura 2. Varios procesos de acondicionamiento de aire 36
Figura 3. Durante un proceso de enfriamiento simple, la humedad absoluta
es constante y la humedad relativa aumenta 37
Figura 4. Proceso de enfriamiento y deshumectación 38
Figura 5. Cuando dos corrientes de aire en los estados 1 y 2 se mezclan
adiabáticamente 40
Figura 6: Influencia sobre las condiciones del aire de salida de un serpentín
evaporador al incrementarse alguno de los parámetros geométricos o de
condiciones del fluido de trabajo 41
Figura 7. Clases de evaporadores (el de arriba, enfriado por agua y el de
abajo, enfriado por refrigerante) 42
Figura 8. Evaporador enfriado por agua 43
Figura 9. Evaporador enfriado por refrigerante 44
Figura 10. Configuraciones y evolución de temperaturas en un intercambiador
de calor 46
Figura 11. Aleta anular de espesor constante 48
Figura 12. Aleta anular de placas continuas 50
Figura13: Arreglo de tubos rectangular y hexagonal 51
Figura 14. Constantes de C para el coeficiente de transferencia de calor
evaporativo del refrigerante 53
Figura 15. Valores de para cada refrigerante 53
Figura 16: Remoción del calor latente y sensible de un serpentín 56
Figura 17: Subenfriamiento en el refrigerante durante la condensación 58
Figura 18: Condición curva del serpentín 61
Figura 19: Geometría de unidad manejadora enfriada por agua 62

14
Figura 20: Evolución del aire al pasar por cada fila del serpentín de UMA
enfriada por agua 63
Figura 21: Balance de energía en una aleta circular 64
Figura 22: Geometría de unidad manejadora enfriada por refrigerante 71
Figura 23: Evolución del aire al pasar por cada fila del serpentín de UMA
enfriada por agua 72
Figura 24. Interfaz principal de la herramienta computacional 80
Figura 25 Selección del análisis de unidades manejadoras 82
Figura 26 Selección del medio refrigerante 83
Figura 27 Pasos a seguir para el análisis de las unidades manejadoras por
catálogos 84
Figura 28 Selección de unidad enfriadora 85
Figura 29 selección del numero de filas 86
Figura 30 seleccionar unidad 87
Figura 31 Ilustración para ver la geometría 88
Figura 32 Geometria de la unidad 88
Figura 33 Ilustración para continuar en el proceso 89
Figura 34 Selección de la unidad enfriadora 90
Figura 35. Pasos a seguir para el análisis de las unidades manejadoras por
catálogos 2 91
Figura 36. Ilustración para la comprobación de la unidad manejadora 92
Figura 37. Fase de comprobación para la unidad seleccionada 93
Figura 38 Fase de comprobación para la unidad seleccionada 2 94
Figura 39 Fase de espera para la comprobación 95
Figura 40. Fase de comprobación para la unidad selecciona 3 96
Figura 41. Fase de comprobación para la unidad selecciona 4 97
Figura 42 Ilustración de la tabla de evolución del aire 98
Figura 43 Ilustración para guardar la tabla de evolución del aire 99
Figura 44 Fase de comprobación para la unidad selecciona 5 100
Figura 45. Fase de proceso de análisis de unidades para calcular 101

15
Figura 46 Fase de proceso de análisis de unidades para calcular 2 102
Figura 47 Cálculo de la unidad manejadora 103
Figura 48. Ilustración del diagrama psicométrico en la unidad a calcular 104
Figura 49. Geometría Modelo 2 YG del catálogo de York 107
Figura 50. Geometría Modelo 3 YG del catálogo de York 108
Figura 51. Geometría Modelo 9 YG del catálogo de York 109
Figura 52. Geometría Modelo 7 YG del catálogo de York 112
Figura 53. Geometría Modelo 7 YG del catálogo de York 115
Figura 54: Unidad Manejadora de Aire marca Paramo 116
Figura 55: Medidor de T y 119
Figura 56: Pistola Infrarrojo 119
Figura 57: Anemómetro 120
Figura 58. Geometría Modelo 4 YG del catalogo de York con NF=3 132
Figura 59. Flujo de aire vs Calor sensible 135
Figura 60. Flujo de aire vs Calor total 136
Figura 61. Geometría del modelo 2 YG del catalogo de York con NF=3 144

16
 LISTA DE TABLAS

PAG.
Tabla 1. Comparación de carga sensible del catálogo modelo 02 Vs Carga
sensible del programa 106
Tabla 2. Comparación de carga sensible del catalogo modelo 03 Vs Carga
sensible del programa 107
Tabla 3. Comparación de carga sensible del catálogo modelo 09 Vs Carga
sensible del programa 108
Tabla 4. Comparación de carga latente del catalogo modelo 02 Vs Carga
latente del programa 110
Tabla 5. Comparación de carga latente del catalogo modelo 03 Vs Carga
latente del programa 111
Tabla 6. Comparación de carga latente del catalogo modelo 07 Vs carga
latente del programa 111
Tabla 7. Comparación de temperatura de bulbo seco del aire de salida del
catalogo modelo 02 Vs temperatura de bulbo seco del aire de salida del
programa 113
Tabla 8. Comparación de temperatura de bulbo seco del aire de salida del
catalogo modelo 03 Vs temperatura de bulbo seco del aire de salida del
programa 113
Tabla 9. Comparación de temperatura de bulbo seco del aire de salida del
catalogo modelo 07 Vs temperatura de bulbo seco del aire de salida del
programa 114
Tabla 10. Comparación de humedad relativa del aire de salida del catálogo
modelo 12 Vs humedad relativa del aire de salida del programa 115
Tabla 11. Comparación de humedad relativa del aire de salida del catalogo
modelo 07 Vs humedad relativa del aire de salida del programa 116
Tabla12. Primera tabla de datos experimentales 121

17
Tabla 13. Cálculos de los datos experimentales 121
Tabla 14. Primera tabla de datos experimentales 122
Tabla 15. Cálculos de los datos experimentales 122
Tabla 16. Tablas de comparación de carga latente de la unidad Paramo
Vs Carga latente del programa 123
Tabla 17. Tablas de comparación de carga total de la unidad Paramo
Vs Carga total del programa 123
Tabla 18. Tablas de comparación para la temperatura de la unidad Paramo
Vs temperatura de salidad del programa 124
Tabla 19. Tablas de comparación para la humded relativa de la unidad
Paramo Vs temperatura de salidad del programa 124
Tabla 20. Tablas de comparación de carga latente de la unidad Paramo
Vs Carga latente del programa 125
Tabla 21. Tablas de comparación de carga total de la unidad Paramo
Vs Carga total del programa 125
Tabla 22. Tablas de comparación para la temperatura de la unidad Paramo
Vs temperatura de salidad del programa 126
Tabla 23. Tablas de comparación para la humedad relativa de la unidad
Paramo Vs humedad relativa de salida del programa 126
Tabla 24. Comparación de carga latente del catalogo Vs Carga latente del
programa para Tev=4.44 ºC 127
Tabla 25. Comparación de carga latente del catalogo Vs Carga latente del
programa para Tev=7.22 ºC 127
Tabla 26. Comparación de Temperatura de salida del catalogo Vs
Temperatura de salida del programa Tev=4,44 ºC 128
Tabla 27. Comparación de Temperatura de salida del catalogo Vs
Temperatura de salida del programa Tev=7,22 ºC 128
Tabla 28. Comparación de carga sensible del catalogo Vs Carga sensible del
programa con Tev=4,44 ºC 129

18
Tabla 29. Comparación de carga sensible del catalogo Vs Carga sensible
del programa con Tev=7,22 ºC 129
Tabla 30. Comparación de humedad relativa del catalogo Vs humedad
relativa del programa con Tev=4,44 ºC 130
Tabla 31. Comparación de humedad relativa del catalogo Vs humedad
relativa del programa con Tev=7,22 ºC 130
Tabla 32. Efecto de la variación del flujo de aire sobre las condiciones de
salida del aire y de los calores extraídos por la unidad 134
Tabla 33. Efecto de la variación del la temperatura del bulbo seco de aire
sobre las condiciones de salida del aire y de los calores extraídos por la
unidad 137
Tabla 34. Efecto de la variación del bulbo húmedo de aire sobre las
condiciones de salida del aire y de los calores extraídos por la unidad 139
Tabla 35. Efecto de la variación del flujo del agua sobre las condiciones de
salida del aire y de los calores extraídos por la unidad 141
Tabla 36. Variación de las condiciones de salida del aire y de los calores
extraídos al paso por cada fila 145
Tabla 37. Variación de la geometría para la unidad YG York 02 147
Tabla 38. Variación de la geometría para R1MX 018 148

19
 LISTA DE ANEXOS

PAG.
ANEXO A. MANUAL TECNICO DE COMFORT PACKAGE FAN & COIL
UNIT- YG CAPACIDADES DE 1360 A 10200 m3/h 156
ANEXO B. CATALOGO DE YORK UNIDADES FANCOILS DE EXPANSIÓN
DIRECTA CAPACIDADES DE 12000 A 36000 BTU/H 191
ANEXO C. CORRELACIÓN GENERAL PARA FLUJO DE DOS FASES
SATURADO PARA TRANSFERENCIA DE CALOR DENTRO DE TUBOS
VERTICALES U HORIZONTALES 195

20
 RESUMEN

TÍTULO:
HERRAMIENTA COMPUTACIONAL PARA EL DESEMPEÑO Y USO DE UN EVAPORADOR
CONCRETO EN AIRE ACONDICIONADO

AUTORES:
YELMIS BARROS GALVIS
DIANA CAROLINA JAIMES BURGOS

PALABRAS CLAVES:
Unidad Manejadora de Aire Acondicionado, Evolución del Aire, Superficie Seca-húmeda,
Temperatura de Bulbo Seco, Temperatura de Bulbo Húmedo, Calor Sensible, Calor Latente.

DESCRIPCIÓN:

Este documento detalla una herramienta computacional desarrollada en Matlab y basada en el

modelo matemático que permite la evolución del estado del aire (temperatura y humedad) y
determina el calor sensible y latente a su paso por el serpentín de la unidad cuando el medio de
enfriamiento es agua (Chiller) o refrigerante (DX de expansión directa) bajo condiciones de
superficie de la unidad manejadora completamente seca o seca-húmeda. El análisis del
modelamiento se realiza por filas del serpentín de la unidad manejadora (dividiendo el
intercambiador en módulos con aletas circulares) donde las salidas de los fluidos de trabajo (aire y
medio de enfriamiento) de una fila es la entrada de la fila consecutiva.

También, se comprueban los resultados de la herramienta computacional con los catálogos de

unidades manejadoras enfriadas con agua y con las unidades manejadoras enfriadas por
refrigerante. Por otra parte, se analizan los resultados teniendo en cuenta el efecto de cambiar las
condiciones de entrada del aire (flujo volumétrico- CFM, temperatura de bulbo seco-TBS,
temperatura de bulbo húmedo-TBH) y las condiciones de entrada de agua (flujo volumétrico-GPM).

Además, para una mejor comprensión de la herramienta computacional se presenta un manual de

funcionamiento ilustrado con las interfaces del programa en donde se indica el procedimiento para
el cálculo de las condiciones de salida de los fluidos y los calores sensible, latente y total de la
unidad manejadora.

Proyecto de grado
Facultad de Ingenierías Físico-mecánicas, Escuela de Ingeniería Mecánica, Director Omar
Armando Gélvez Arocha

21
 SUMMARY

TITLE:
COMPUTATIONAL TOOL FOR THE PERFORMANCE AND USE OF A AIR HANDLE UNIT IN
AIR-CONDITIONING

AUTHORS:
YELMIS BARROS GALVIS
DIANA CAROLINA JAIMES BURGOS

KEY WORDS:
Air-conditioning Handle Unit, Air Evolution, Dry-Wet Surface, Temperature Dry Bulb, Temperature
Wet Bulb, Sensible Heat, Latent Heat.

DESCRIPTION:

This document details a computational tool developed in Matlab based on the mathematical model
that allows for the evolution of the state of air (temperature and humidity) and determines the
sensible heat and latent heat removed from it passes through the coil of the air handle unit when
the cooling medium is water (Chiller) or refrigerant (DX Direct Expansion) under surface conditions
of the unit completely dry or dry-wet. The modeling analysis performed by row coil of air handler unit
(dividing the heat exchanger into many circular fins module) where the output of the working fluids
(air and cooling medium) of a row is the entrance to the next row.

Also, compare the results of the computational tool with catalogs air handler units with water cooled
and air handler units with refrigerant. On the other hand, the results are analyzed taking into
account the effect of changing the air inlet conditions (flow rates, CFM, dry-bulb temperature TBS,
wet-bulb temperature TBH) and water inlet conditions (volumetric flow-GPM).

Additionally, for better comprehension of the computational tool is presented a manual operation of
the computational tool which interfaces the program illustrates the procedure for calculating the
output conditions of the fluids and sensible heat, latent and total of the air handler unit.

Graduation Project
Faculty of Physics and Mechanical Engineering, School of Mechanical Engineering,
Director Omar Armando Gélvez Arocha

22
 INTRODUCCIÓN

Que dentro de un recinto cerrado se encuentre el adecuado confort para que las
personas se sientan a gusto se requiere que todos los elementos del sistema
funcionen adecuadamente, sin embargo es usual encontrar que en lugares donde
pareciera que durante una parte del día el equipo funcionara bien después de
cierto tiempo no se siente igual.

Surge así la pregunta ¿será que no es capaz el equipo de sacar la carga

instalada?.

Y qué pasaría si se trasladara esa misma unidad a otra ciudad, cuál sería el
comportamiento, que calor latente o que calor sensible puede extraer dicha
unidad.

Es a partir de estas preguntas que suscita la necesidad de plantear una solución

para calcular el desempeño del corazón del sistema de refrigeración de diferentes
unidades manejadoras, teniendo en cuenta la transferencia de calor y el fenómeno
de deshumidificación que usualmente tienen lugar en ella. Razón por la cual este
proyecto cobija el cálculo de capacidades tanto para la unidad manejadora de aire
enfriada por agua así como la enfriada por refrigerante, usando una herramienta
computacional desarrollada en Matlab.

Este libro de proyecto de grado está compuesto así: en el capítulo 1, contiene todo
el marco teórico que soporta los modelos matemáticos propuestos para

23
Ofrecer confort a las personas que se encuentran en un recinto cerrado requiere
que todos los elementos del sistema de aire acondicionado funcionen
adecuadamente. Sin embargo, es usual encontrarse en lugares en donde
pareciere que el equipo que durante una parte del día funciona bien luego de un
tiempo funciona mal ya que no es capaz de sacar las cargas con la unidad
instalada. O con esa misma unidad trasladada a otra ciudad, qué calor sensible,
qué calor latente puede extraer dicha unidad. Estas junto con otras preguntas
suscitaron la necesidad de plantear una solución para el calcular el desempeño
del corazón del sistema de refrigeración, la unidad manejadora teniendo en cuenta
la transferencia de calor y el fenómeno de deshumidificación que usualmente tiene
lugar en ella. Razón por la cual este proyecto cobija el cálculo de capacidades
tanto para la unidad manejadora de aire enfriada por agua así como la enfriada
por refrigerante, usando una herramienta computacional desarrollada en Matlab.

Este libro de proyecto de grado está compuesto así: en el capítulo 1, contiene todo
el marco teórico que soporta los modelos matemáticos propuestos para
acondicionar el aire a su paso por la unidad manejadora de acuerdo al los fluidos
de trabajo: agua o refrigerante. En el capítulo 2, se encuentra desarrollado el
modelo (los parámetros conocidos, los sistemas de ecuaciones así como el
algoritmo y diagrama de flujo). Mientras que en el capítulo 3, se detalla el manual
de usuario que consigna los pasos a seguir para usar adecuadamente la
herramienta computacional. Y en el último capítulo, se analizan los resultados de
extraídos de la herramienta con los datos obtenidos de catálogos y de las pruebas
realizadas en la Unidad Manejadora de la sala de computación Edificio de
Ciencias Humanas. Posteriormente, se encuentran las conclusiones y
recomendaciones del proyecto de grado.

Ya desarrollados los objetivos del proyecto, se espera que los resultados fruto de
la herramienta puedan servir para compañeros como fuente de consulta para
trabajos académicos posteriores.

24
25
 1. FUNDAMENTOS DEL AIRE ACONDICIONADO

1.1 AIRE HÚMEDO COMO GAS IDEAL1

El aire en la atmósfera contiene una cierta cantidad de vapor de agua (o humedad)

y se conoce como aire atmosférico ó aire húmedo, el cual se trata
convenientemente como un gas ideal en el intervalo donde el nitrógeno, el oxígeno
(aire seco) y el vapor de agua se comportan como un gas ideal. Por lo cual
obedecen las siguientes relaciones para la mezcla ideal:

Donde:
Masa molar promedio de la mezcla
Masa molar de cada gas
Fracción molar de cada gas
Fracción másica de cada gas
P: Presión parcial de cada gas
Presión total de la mezcla
Calor específico de la mezcla
Calor especifico de cada gas

1
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 12 Mezclas de gases.

26
Es conveniente tratar al aire como una mezcla de vapor de agua y de aire seco,
porque la composición del aire seco permanece relativamente constante, pero la
cantidad de vapor de agua cambia por la condensación y evaporación de océanos,
ríos, etc. A pesar de que el contenido del vapor del agua en el aire es pequeña,
desempeña un importante papel para la comodidad humana.

1.2 PROPIEDADES DEL AIRE ATMOSFÉRICO2

1.2.1 Temperatura de punto de rocío

El exceso de humedad en el aire se condensa en las superficies frías y forma el

rocío. La temperatura a la cual se inicia la condensación si el aire se enfría a
presión constante, se le llama punto de rocío (Tdp).

1.2.2 Temperatura de bulbo húmedo

Temperatura que se obtiene de un termómetro cuyo bulbo está cubierto con una
mecha de algodón saturada con agua y sopla aire sobre ella. Se representa con la
abreviatura Tbh.

1.2.3 Temperatura de bulbo seco

Temperatura que se obtiene de un termómetro cuyo bulbo está seco. Se

representa con la abreviatura Tbs.

2
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 13 Mezcla de gas-vapor y
Acondicionamiento de aire, y de la Termodinámica Aplicada de Maradey, capítulo 2 Equilibrio de Fases

27
1.2.4 Humedad específica y relativa del aire

La cantidad de vapor de agua en el aire puede determinarse a partir de la

humedad absoluta y se le representa por medio de :

Donde:
Masa del vapor presente en el aire atmosférico
Masa del aire seco
Presión parcial de vapor

La cantidad de humedad en el aire tiene un efecto definitivo en las condiciones de

comodidad que ofrece un ambiente. Sin embargo, el nivel de comodidad depende
más de la cantidad humedad que el aire contiene ( ) en relación con la cantidad
máxima que puede obtener a la misma temperatura ( ). La relación entre estas
dos cantidades se conoce con el nombre de humedad relativa ( ).

Presión de saturación medida a la temperatura de Tbs

Si se combinan (5) y (6) se obtiene:

28
Otro modo del calcular la humedad del aire húmedo [kgv/kgas], es usando la
siguiente ecuación:

Conociendo:
La temperatura de bulbo húmedo [K]
La temperatura de bulbo seco [K]
: La humedad absoluta medida a la temperatura de bulbo húmedo y con
humedad relativa de 100% [kgv/kgas]

Una forma aproximada para calcular la presión de vapor de agua [kPa],

conociendo la temperatura de bulbo seco [K], es la ecuación de Antoine:

También, la presión de vapor se puede calcular según la ASHRAE3, de acuerdo a

la siguiente fórmula:

Donde:

C1= -5.8002206e+03
C2= 1.3914993e00
C3= -4.8640239e-02
C4= 4.1764768e-05
C5= -1.4452093e-08

3
Tomada de ASHRAE 1997 Fundamentals Handbook HVAC. Chapter 6 Psychrometrics. Pág. 93

29
 C6= 6.5459673e00

La entalpía total del aire atmosférico es la suma de las entalpias del aire seco y del
vapor del agua:

Sin embargo, como el contenido de vapor de agua cambia es mejor dejarlo

expresado por unidad de masa de aire seco. Otro modo de calcular
aproximadamente la entalpía del aire húmedo [kJ/kgas], conociendo la humedad
absoluta [kgv/kgas] y la temperatura de bulbo seco [K], es:

1.3 LA CARTA PSICROMÉTRICA4

El estado del aire atmosférico a una presión específica se establece por completo
mediante dos propiedades intensivas independientes. El resto de las propiedades
se calcula fácilmente a partir de las relaciones anteriores. La presentación de los
cálculos luego de efectuados, es con gráfica de la carta psicométrica (Ver figura
1).

Las temperaturas de bulbo seco se muestran sobre el eje horizontal y la humedad

específica sobre el eje vertical. En el extremo izquierdo de la carta hay una curva
(llamada línea de saturación) en lugar de una línea recta. Todos los estados de
aire saturado se localizan en esta curva. Por tanto, es también la curva del 100%

4
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 13 Mezcla de gas-vapor y
Acondicionamiento de aire.

30
de humedad relativa. Otras curvas de humedad relativa constantes tienen la
misma forma en general. Las líneas de temperatura de bulbo húmedo constante
tienen una apariencia descendente hacia la derecha. Las líneas de volumen
específico constante (en m3/kg de aire seco) parecen similares, salvo que son más
inclinadas. Las líneas de entalpía constante (en kJ/kg de aire seco) están casi
paralelas a las líneas de temperatura de bulbo húmedo constante. Para aire
saturado las temperaturas de bulbo seco, de bulbo húmedo y de punto de rocío
son idénticas.

Figura 1. Carta psicrométrica

Fuente: Manual de aire acondicionado. Capítulo 8. Empleo del Diagrama Psicrometrico. Carrier Air
Conditioning Company. Boixareu Editores.1980. Págs. 109-125.

31
1.3.1 Ciclo de evolución del aire5

El ciclo de evolución del aire climatizado puede representarse sobre el diagrama

psicométrico, tal como se indica en la figura 3. El aire en el estado (3), mezcla de
aire exterior (2) y de retorno (1), pasa a través del aparato acondicionador, y su
evolución se representa por la línea (3-4). Abandona el aparato en (4) y es
impulsado hacia el local donde absorbe calor y humedad, según la transformación
(4-1). En general, gran parte del aire impulsado vuelve a recogerse para su mezcla
con el aire exterior. La mezcla pasa a través del aparato donde abandona la
humedad y calor recibidos, al objeto de mantener las condiciones deseadas.

Para calcular las cargas requeridas tanto para el recinto a acondicionar y para el
aire atmosférico, suele hacerse a través de las siguientes fórmulas aproximadas:

Razón del calor sensible:

Calor sensible del aire exterior:

Razón del calor latente:

5
Manual de aire acondicionado. Capítulo 8. Empleo del Diagrama Psicrometrico. Carrier Air Conditioning
Company. Boixareu Editores.1980. Págs. 109-125.

32
 Calor latente del aire exterior

Total calor sensible

Total calor latente

Calor total

1.3.2 Confort térmico y acondicionamiento de aire6

Los seres humanos tienen una debilidad inherente: necesitan sentirse cómodos.
Quieren vivir en un ambiente que no sea ni caliente ni frío, ni muy húmedo, ni muy
seco. Sin embargo, la comodidad no se obtiene con facilidad ya que las
necesidades del cuerpo humano y el clima suelen no ser muy compatibles. Para
alcanzar la comodidad es necesario luchar constantemente contra los factores que
causan la incomodidad. Se puede definir el confort ambiental como un estado de
satisfacción del ser humano, respecto al medio en que vive, denominándose
ergonomía a la ciencia que estudia los problemas de adaptación del hombre al
ambiente. Las variables de las cuales depende esa sensación de bienestar
pueden ser:

Individuales

6
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 13 Mezcla de gas-vapor y
Acondicionamiento de aire y Sistemas de Aire Acondicionado, Primera edición, Néstor Quadri, capítulo 2
Confort Térmico.

33
 Ambientales

Las variables individuales del confort dependen de la característica particular del

ser humano y de muchos factores como ser el nivel de actividad, tipo de ropa,
sexo, edad, estado de salud y para su análisis se debe considerar el intercambio
de calor del cuerpo humano con el ambiente. La forma de que se vale el cuerpo
humano para eliminar el calor remanente es la siguiente:

Conducción: a través de la piel y los vestidos del individuo.

Convección: desde la periferia de la piel al aire que la circunda.
Radiación: mediante la emisión de calor del cuerpo a las superficies frías del
entorno del local.
Evaporación: por la exudación de la piel y una pequeña parte contenida en el
aire de respiración.

En cambio, las variables ambientales están relacionadas con las modificaciones a

producir al clima del local a acondicionar y los parámetros básicos que debe
controlar un sistema de climatización, a fin de lograr el bienestar son la
temperatura del aire y superficiales, la humedad relativa, el movimiento y la
calidad del aire interior.

1.3.3 Condiciones del local7

Las características del ambiente para lograr las condiciones de confort son muy
importantes, si bien el cuerpo humano a causa de su mecanismo de regulación
mantiene el equilibrio térmico dentro de un gran margen de variaciones. El confort
térmico está íntimamente relacionado con las condiciones del clima del local a

7
Tomado de Sistemas de Aire Acondicionado, Primera edición, Néstor Quadri, capítulo 2 Confort Térmico.

34
acondicionar y los parámetros básicos que debe controlar un sistema de
climatización, a fin de lograrlo son:

Temperatura del aire y superficiales

Humedad relativa
Movimiento del aire

1.3.4 Procesos de acondicionamiento de aire8

Para mantener una vivienda o una construcción industrial a la temperatura y

humedad deseadas son necesarios algunos procesos llamados acondicionamiento
de aire. Estos procesos incluyen el calentamiento simple (elevar la temperatura),
el enfriamiento simple (reducir la temperatura), la humidificación (agregar
humedad) y la deshumidificación (eliminar humedad). Algunas veces dos o más de
estos procesos son necesarios para llevar el aire al nivel de temperatura y
humedad que se desea (Ver figura 2).

8
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 13 Mezcla de gas-vapor y
Acondicionamiento de aire

35
 Figura 2. Varios procesos de acondicionamiento de aire

Fuente: Termodinámica de Cengel, Quinta edición en inglés, capitulo 14 Mezcla de gas-vapor y

aire acondicionado (pág. 729)

1.3.4.1 Enfriamiento del aire húmedo. Durante un proceso de enfriamiento a

humedad específica constante del aire, la temperatura de bulbo seco disminuye y
la humedad relativa aumenta, sin incluir deshumidificación del aire (Ver figura 4).
Esto se produce al pasar el aire sobre algunos serpentines por los cuales fluye un
refrigerante o agua fría. Durante el proceso la humedad se mantiene constante:

36
 Figura 3. Durante un proceso de enfriamiento simple, la humedad absoluta
es constante y la humedad relativa aumenta

Fuente: Termodinámica de Cengel, Quinta edición en inglés, capitulo 14 Mezcla de gas-vapor y

aire acondicionado (pág. 730)

1.3.4.2 Enfriamiento con deshumidificación del aire húmedo. El proceso de

enfriamiento con deshumidificación se muestra de manera esquemática en la carta
psicrométrica. El aire caliente y húmedo entra en la sección de enfriamiento en el
estado 1(Tbs=30ºC y ф=80%). Cuando pasa por los serpentines de enfriamiento,
su temperatura disminuye y su humedad relativa aumenta a humedad especifica
constante. Si la sección de enfriamiento tiene la longitud suficiente, el aire
alcanzará su punto de rocío (estado x, aire saturado, wx=w1 y ф=100%). El
enfriamiento adicional del aire originará la condensación de parte de la humedad
en el aire. Éste permanece saturado durante todo el proceso de condensación,
que sigue una línea de 100% de humedad relativa hasta que alcanza el estado
final (estado 2). El vapor de agua, que se condensa fuera del aire durante el
proceso, se elimina de la sección por medio de un canal independiente. Suele
considerarse que el condensado sale de la sección de enfriamiento a T 2.

El aire frío y saturado en el estado 2 casi siempre se envía directamente a la

habitación, donde se mezcla con el aire que hay ahí. Puede ocurrir que el estado
del aire a la salida de los serpentines tenga una temperatura demasiado baja por

37
lo cual, se le hace pasar por una sección de calentamiento, antes de ser enviado
al recinto.

Figura 4. Proceso de enfriamiento y deshumectación

Fuente: Fuente: Termodinámica de Cengel, Quinta edición en inglés, capitulo 14 Mezcla de gas-
vapor y aire acondicionado (pág. 733)

Aplicado la ecuación de continuidad al flujo de aire seco y al vapor de agua se

tiene:

Tomando el volumen de control de la figura anterior, y realizando el balance de

energía (subíndice i es entrada y subíndice e es salida) se tiene:

Para lograr el confort en los locales, es necesario enfriar y además, deshumidificar

el aire, porque el porcentaje de humedad relativa aumentaría en forma
considerable, provocando una sensación de molestia y pesadez. Esta función se
efectúa en un proceso único y simultáneo por medio de un serpentín o batería de

38
enfriamiento en las cuales se absorbe el calor sensible y también el calor latente
del aire, por efecto de la condensación sobre su superficie del vapor de agua que
contiene, debido a que se lo enfría por debajo de su punto de rocío.

1.3.3.3 Mezcla adiabática de corrientes de aire o ventilación9. En muchas

aplicaciones del acondicionamiento de aire es necesaria la mezcla de dos
corrientes de aire. Esto es particularmente cierto en los grandes edificios, la mayor
parte de las plantas de producción y proceso, y los hospitales, donde se requiere
que el aire acondicionado se mezcle con una fracción del aire fresco exterior antes
de enviarse a espacios habitados. En un sistema de aire acondicionado, la función
de ventilación, consiste en la entrada de aire exterior, para renovar
permanentemente el aire de recirculación del sistema en las proporciones
necesarias a fin de lograr un adecuado nivel de pureza, dado que como resultado
del proceso respiratorio, se consume oxigeno y se exhala anhídrido carbónico, por
lo que debe suministrarse siempre aire nuevo a los locales para evitar que se
produzcan malos olores.

El aire nuevo del exterior del edificio o aire de ventilación penetra a través de una
reja de toma de aire en un recinto llamado pleno de mezcla, donde se mezcla el
aire nuevo con el aire de retorno de los locales, regulándose mediante persianas
de accionamiento manual o automático.

9
Tomado de la Termodinámica de Cengel, Cuarta Edición, capítulo 13 Mezcla de gas-vapor y
Acondicionamiento de aire y Sistemas de Aire Acondicionado, Primera edición, Néstor Quadri, capítulo 2
Confort Térmico.

39
 Figura 5. Cuando dos corrientes de aire en los estados 1 y 2 se mezclan
adiabáticamente

Fuente: Fuente: Termodinámica de Cengel, Quinta edición en inglés, capitulo 14 Mezcla de gas-
vapor y aire acondicionado (pág. 733)

Al realizar balances de masa y energía para el volumen de control tomado en la

figura se puede obtener la siguiente proporción para las corrientes de aire:

1.4 INTERCAMBIADORES DE CALOR10

Los intercambiadores de calor son aparatos que facilitan el intercambio de calor

entre dos fluidos que se encuentran a diferentes temperaturas y evitan al mismo
tiempo que se mezclen entre sí. Para la aplicación del aire acondicionado se
requiere un intercambiador de tipo compacto, en el cual se tiene una gran área
superficial de transferencia de calor por unidad de volumen. En ellos, se logran
velocidades de transferencia de calor entre dos fluidos en un volumen pequeño y
son de uso común en aplicaciones con limitaciones estrictas con respecto al peso

10
Tomado de la Transferencia de Calor de Cengel, segunda edición, capítulo 13 Intercambiadores de Calor y
del Handbook de Refrigeración de ASHRAE 81678, capítulo 15 Sistemas de aire: componentes-ventiladores,
serpentines, filtros y deshumidificadores.

40
y el volumen de esos aparatos. Por ejemplo, los evaporadores refrigeradores y los
evaporadores bombas de calor.

Pese a ser intercambiadores símiles, entre ellos existen dos marcadas diferencias
como son la densidad de aletas y la dirección del flujo de aire. En refrigeradores,
las densidades de aletas son generalmente menos que los evaporadores bombas
de calor. A diferencia de las bombas de calor donde el aire fluye por el evaporador
en flujo cruzado, el flujo de aire en los refrigeradores principalmente es a
contraflujo donde el aire fluye en dirección vertical en el evaporador. También
sobre estos equipos térmicos, el diseñador y operador del sistema tiene algún
control sobre las condiciones de salida del aire a acondicionar tal como lo muestra
la siguiente tabla:

Figura 6: Influencia sobre las condiciones del aire de salida de un serpentín

evaporador al incrementarse alguno de los parámetros geométricos o de
condiciones del fluido de trabajo

Fuente: Handbook de Refrigeración Industrial de Stoecker Wilbert. Pág.203

41
1.4.1 Tipos de intercambiadores para aire acondicionado.

De acuerdo al tipo de medio de enfriamiento (Ver figura 7):

Figura 7. Clases de evaporadores (el de arriba, enfriado por agua y el de

abajo, enfriado por refrigerante)

Fuente: Handbook de Refrigeración de ASHRAE 81678, capítulo 15 Sistemas de aire:

componentes-ventiladores, serpentines, filtros y deshumidificadores (Pág. 34)

1.4.1.1 Evaporador enfriado por agua. Para mantener una alta rata de
transferencia de calor, el aire y el agua normalmente sigue un arreglo a cruzado.
La diferencia entre un serpentín enfriado con agua y el enfriado con refrigerante
está principalmente en la carencia de una región supercalentada en el serpentín
enfriado con agua y la variación de la temperatura del agua enfriada, lo cual es
muy menudo más grande que la variación en la temperatura de evaporación en la
región de dos fases en un serpentín DX.

42
 Figura 8. Evaporador enfriado por agua

Fuente: Handbook de Refrigeración de ASHRAE 81678, capítulo 15 Sistemas de aire:

componentes-ventiladores, serpentines, filtros y deshumidificadores (Pág. 35)

1.4.1.2 Evaporador de expansión directa DX o enfriado por refrigerante.El

refrigerante es evaporado y expandido directamente dentro de los tubos para
enfriar y deshumidificar el aire que fluye por él, y la condensación ocurre sobre la
superficie externa del evaporador. Para igual distribución del refrigerante y reducir
la caída de presión, los senderos de flujo de refrigerante en un serpentín DX están
siempre divididos dentro de un número de circuitos, de acuerdo al ancho aleteado
del serpentín. El refrigerante es usualmente suplido por la válvula de expansión

43
termostática. Entonces fluyen por el distribuidor, el cual distribuye refrigerante de
modo uniforme. Los distribuidores tipo Venturi son usualmente usados.

Figura 9. Evaporador enfriado por refrigerante

Fuente: Handbook de Refrigeración de ASHRAE 81678, capítulo 15 Sistemas de aire:

componentes-ventiladores, serpentines, filtros y deshumidificadores (Pág. 36)

44
1.4.2 Análisis de los intercambiadores de calor11

En la práctica de los intercambiadores de calor son de uso común y un ingeniero

se encuentra a menudo en la posición de seleccionar un intercambiador de calor
que logre un cambio de temperatura específica de una corriente de fluido de gasto
de masa conocido, o bien, de predecir las temperaturas de salida de las corrientes
de de fluido caliente y del frío en un intercambiador de calor especifico.

Para el cálculo de la transferencia de calor, teniendo en cuenta el potencial de

temperatura se deriva la siguiente fórmula:

Donde:
= Coeficiente global de la transferencia de calor
= Área de transferencia de calor
= Diferencia media entre la temperatura de los fluidos

Sin embargo, no es práctico el uso de esta ecuación para un evaporador pues

tanto como son variables a su paso por el serpentín. Para usarla, se
requiere asumir un valor promedio de y la determinación de una diferencia
promedio conveniente.

1.4.2.1 Método de la LMTD. En este método se considera constante y se

desprecian las pérdidas de calor, ganancias externas desde y hacia el evaporador.
Y se obtiene la misma expresión para LMTD, tanto en flujo paralelo, como en
contracorriente (ver figura 10) y es:

11
Tomado de Transferencia de Calor de Incropera, cuarta edición, capítulo 11 Intercambiadores de calor.

45
Según la configuración de flujo:

En caso de flujo paralelo: DT1 = Th,ent - Tc,ent ; DT2 = Th,sal - Tc,sal

En caso de contraflujo: DT1 = Th,ent - Tc,sal ; DT2 = Th,sal - Tc,ent

Donde:
Th,ent: temperatura de entrada del fluido caliente
Th,sal: temperatura de salida del fluido caliente
Tc,ent: temperatura de entrada del fluido frio
Tc,sal: temperatura de salida del fluido frío

Figura 10. Configuraciones y evolución de temperaturas en un

intercambiador de calor

Fuente:http://www.telecable.es/personales/albatros1/calor/transferencia_de_calor_07_intercambiad
ores.htm

46
1.4.2.2 El coeficiente global de la transferencia de calor. Por lo común un
intercambiador de calor está relacionado con dos fluidos que fluyen separados por
una pared sólida. En primer lugar, el calor se transfiere del fluido caliente hacia la
pared por convección, después a través de la pared por conducción y, por último
de la pared hacia el fluido frío de nuevo por convección. Cualesquiera efectos de
la radiación suelen incluirse en los coeficientes de transferencia de calor por
convección.

Si el intercambiador sin aletas es una pared bañada por dos fluidos a diferente
temperatura, y cuyo coeficiente global de transferencia tiene la siguiente
expresión:

Siendo Ai el área de la superficie interior, Am es el área afectada por la

conducción el y Ao el área de la superficie exterior de transferencia de calor.
-Siendo hi el coeficiente de película interior y ho el coeficiente de película
exterior.
Siendo k el coeficiente de conducción del material de la pared

1.4.3 Transferencia de calor de intercambiadores con aletas12

La velocidad de la transferencia de calor desde una superficie que está a una

temperatura Ts hacia el medio circundante que está a T∞ se expresa por la ley de
Newton del enfriamiento como:

12
Tomado de la Transferencia de Calor de Cengel, capítulo 3 Conducción de calor en estado estable y
Transferencia de Calor, capítulo X Superficies ampliadas de sección transversal variable.

47
Donde es el área superficial de transferencia de calor y h es el coeficiente de
transferencia de calor por convección. Cuando las temperaturas y se fijan
por diseño, como por consideraciones de diseño, como con frecuencia es el calor,
existen dos maneras de incrementar la velocidad de la transferencia de calor:
aumentar el coeficiente de transferencia de calor por convección, h, o aumentar el
área superficial . La alternativa es aumentar el área superficial al agregar unas
superficiales extendidas llamadas aletas, hechas de materiales intensamente
conductores como el aluminio. En general, el coeficiente h varía a lo largo de la
aleta así como de su circunferencia y que su valor en un punto es una fuerte
función del movimiento del fluido en ese punto.

Figura 11. Aleta anular de espesor constante

Fuente: Transferencia de calor de Pedro Fernández Diez, capítulo X Superficies

ampliadas de sección transversal variable (pág. 144)

1.4.4 Parámetros característicos de una aleta

1.4.4.1 Eficiencia de una aleta. Es la relación entre la potencia térmica (Q-punto)

que se disipa en la misma y la potencia térmica que se disiparía si toda la aleta
estuviese a una temperatura igual a la de la base (la temperatura de la aleta será
inferior a la de la base):

48
1.4.4.2 Eficiencia superficial. Es la relación del calor total transmitido y del que
se transmitirá si la superficie total de la aleta se mantuviese a la misma
temperatura de la base:

1.4.4.3 Efectividad de una aleta. Es la relación entre la potencia térmica (Q-

punto) que se disipa en la misma y la potencia térmica que se disipa sin aleta
desde el área de la base que ocupa ésta en la superficie primaria:

La eficiencia de la aleta circular está determinada por:

En donde:

49
Donde:
k: conductividad térmica de la aleta
p: perímetro de la sección transversal de la aleta
h: coeficiente de película
At: área de la sección transversal de la aleta.
Las placas de aletas continuas son muy usadas en los intercambiadores de calor.
En este caso cada aleta se extiende de tubo a tubo en todas las direcciones.

Figura 12. Aleta anular de placas continuas

Fuente: Handbook de Refrigeración de ASHRAE 81678, capítulo 15 Sistemas de aire:

componentes-ventiladores, serpentines, filtros y deshumidificadores (Pág. 37)

No hay solución analítica precisa para este tipo de aleta. Sin embargo, se utilizan
los siguientes tipos de aproximación: arreglo rectangular y hexagonal.

50
 Figura13: Arreglo de tubos rectangular y hexagonal

Fuente: Heating, Ventilating, and Air Conditioning Analysis and Design, Fourth Edition. McQuiston
Faye C, chapter 14 Extended surface heat exchangers (Pág. 579)

1.4.5 Coeficiente de transferencia de calor dentro de los tubos

Para aplicaciones de HVAC los intercambiadores de calor suelen usar como fluido
de trabajo agua, refrigerantes, entre otros.

1.4.5.1 Flujo turbulento del líquido dentro de los tubos. La correlación de

mayor exactitud para hallar el coeficiente de convección interna es la ecuación de
Gnieliski:

(1.38)

Y el factor de fricción es:

(1.39)

Donde:
: Reynolds interno

51
 : Número de Prandtl
: Conductividad del fluido de trabajo
: Diámetro interno

1.4.5.2 Coeficiente de transferencia de calor de evaporación del

refrigerante13. Una sencilla relación para describir el desempeño de la
transferencia de calor en un evaporador sobre el completo rango de calidad del
refrigerante. Es dada por Kandlikar 1990, como sigue:

Donde:
Coeficiente de transferencia de calor convectivo solo líquido

Las constantes con los subíndices de C, se extraen de la siguiente tabla:

13
Boilers, evaporators, and Condenser. Chapter 12. Evaporators and condensers for refrigeration and air-
conditioning systems.Editorial John Wiley & Sons, Inc. 1991. Págs. 647-678.

52
 Figura 14. Constantes de C para el coeficiente de transferencia de calor
evaporativo del refrigerante

Fuente: Boilers, evaporators, and Condenser. Chapter 12. Evaporators and condensers for
refrigeration and air-conditioning systems.Editorial John Wiley & Sons, Inc. 1991. Págs. 647-678.

Y de acuerdo a la clase de refrigerante se obtiene de la tabla como se muestra:

Figura 15. Valores de para cada refrigerante

Fuente: Boilers, evaporators, and Condenser. Chapter 12. Evaporators and condensers for
refrigeration and air-conditioning systems.Editorial John Wiley & Sons, Inc. 1991. Págs. 647-678.

El número de Boiling se obtiene así:

53
El número de convección es:

La calidad del refrigerante suele tomarse , de acuerdo a la Figura 12.9.

(Libro de Evaporator coils, reboilers and condensers. Capítulo 12 Evaporator coils.
Págs.676-678).

El número de Froude es:

El número de Reynolds para todo el fluido de refrigerante:

(1.45)

1.4.6 Coeficiente de transferencia de calor externo en bancos de tubos14

El coeficiente de transferencia de calor de la primera fila de tubos es mayor que el

que corresponde a un solo tubo en un flujo cruzado ya que el fluido debe
acelerarse para pasar por los espacios intermedios.

En un haz de tubos el coeficiente de transferencia de calor va en aumento desde

la primera fila hasta alrededor de la quinta. El número de Nusselt promedio
en un haz de tubos con 10 o más filas se puede calcular a partir de la

relación:

14
Tomado de Transferencia de Calor de Mills, capítulo 4 Fundamentos y correlaciones de la convección.

54
Donde es el número de Nusselt de la primera fila y es un factor de arreglo.
Definimos el paso transversal adimensional como , el paso longitudinal
adimensional como y un factor como:

Si

Si

Entonces, los factores de arreglo están dados por:

Si el haz tiene menos de 10 filas podemos aplicar una sencilla fórmula de

interpolación, (donde N es el número de filas):

La correlación de transferencia de calor propuesta por Churchill y Bernstein:

(Para ReD<10000)

(Para ReD>10000)

55
1.5 TRANSFERENCIA SIMULTÁNEA DE CALOR Y MASA15

Cuando la superficie de un intercambiador de calor que entra en contacto con aire

húmedo se encuentra a una temperatura inferior a la temperatura de rocía del aire,
se presenta condensación de agua. Usualmente la temperatura de bulbo seco y la
humedad relativa disminuyen a medida que el aire fluye a través del
intercambiador. Así pues la transferencia de calor sensible y latente ocurre
simultáneamente. Los fenómenos de transporte tienen lugar en aquellos
procesos, conocidos como procesos de transferencia, en los que se establece el
movimiento de una propiedad (masa, momentum o energía) en una o varias
direcciones bajo la acción de una fuerza impulsora. Tal como se aprecia en la
figura de 16, en la que se ven actuar tanto potencial de energía (Tbs) como de
masa (W) en la propiedad de entalpía para la mezcla de aire seco y vapor de
agua.

Figura 16: Remoción del calor latente y sensible de un serpentín

Fuente: Handbook de Refrigeración Industrial de Stoecker Wilbert. Pág.207

15
Tomado de Heating, Ventilating, and Air Conditioning Analysis and Design, Fourth Edition. McQuiston
Faye C, chapter 14 Extended surface heat exchangers y de Handbook de Refrigeración de ASHRAE 81678,
capítulo 15 Sistemas de aire: componentes-ventiladores, serpentines, filtros y deshumidificadores.

56
Es de práctica común en refrigeración y el aire acondicionado tener en cuenta el
Número Lewis que relaciona para este fenómeno como se expresa a continuación:

1.6 CAPACIDAD DE UN SISTEMA DE REFRIGERACIÓN

La capacidad de un sistema de refrigeración corresponde a la cantidad de calor

que el refrigerante es capaz de sacar del medio refrigerado y es igual al producto
de la masa del refrigerante circulado por unidad de tiempo por la diferencia de
entalpias del refrigerante entre la salida y entrada del evaporador.

1.6.1 Efecto refrigerante16

El aumento de la entalpía del refrigerante en el evaporador se conoce como el

efecto refrigerante (ER). Debido a que representa asimismo la cantidad de calor
removido del medio que se debe enfriar por cada libra o kilogramo de refrigerante
que fluye. Esto se deduce de la ecuación de la energía. Esto es:

En donde:
Efecto de refrigeración
= Entalpía del refrigeración a la salida del evaporador

16
Tomado del libro Principios y Sistemas de Refrigeración. Termodinámica del ciclo de refrigeración por
compresión de vapor. Págs.86-87. Pita, Edward G.Editorial Limusa. 2006

57
 = Entalpía del refrigeración a la salida del evaporador

1.6.2 El subenfriamiento del refrigerante17

El líquido refrigerante condensado es usualmente subenfriado a una temperatura

más baja que la temperatura saturada correspondiente a la temperatura de
condensación del refrigerante. Cuando se presenta éste, el efecto refrigerante
aumenta con el subenfriamiento, y que se produce menor cantidad de gas de
vaporización. También, resulta en una potencia más baja requerida por el
compresor, por unidad de capacidad, y un mayor coeficiente de rendimiento.
Además, disminuye el desplazamiento requerido por el compresor. Desde el punto
de vista de la conservación de energía, el subenfriamiento puede fácilmente
reducir el consumo de energía entre un 5% y 15%. Sin embargo, los grados de
subenfriamiento dependen principalmente de la temperatura del refrigerante
durante la condensación y la construcción y capacidad del condensador.

Figura 17: Subenfriamiento en el refrigerante durante la condensación

Fuente: Tomado del libro Handbook de Aire Acondicionado y refrigeración de la ASHRAE. Capítulo
9. Refrigerantes, ciclos de refrigeración y sistemas de refrigeración. Pás. 26-27. Año 2001

17
Tomado del libro Handbook de Aire Acondicionado y refrigeración de la ASHRAE. Capítulo 9.
Refrigerantes, ciclos de refrigeración y sistemas de refrigeración. Pás. 26-27. Año 2001

58
La entalpía del líquido refrigerante , puede ser calculada:

Donde:
= Entalpía del refrigerante de líquido saturado a la temperatura de
condensación
Calor especifico del líquido refrigerante a presión constante
= Temperatura saturada del líquido refrigerante a la presión de condensado
= Temperatura del refrigerante líquido subenfriado

59
 2. MODELAMIENTO MATEMÁTICO PARA UNIDADES MANEJADORAS DE
AIRE

Para determinar los calores tanto sensible como latente y las condiciones de
salida del aire al paso por la unidad enfriadora, en función de los parámetros
geométricos de un evaporador (indicados en la figura 16) y las condiciones de los
fluidos a la entrada de él, como temperatura de bulbo seco y temperatura de bulbo
húmedo para el aire y la temperatura de entrada del fluido de trabajo (agua o
refrigerante). Se construyó un modelo para cada fluido de trabajo que tiene en
cuenta la evolución del aire a su paso por el serpentín del evaporador (ver figura
15), dicho de otro modo, tiene en cuenta si hay condensación del vapor de agua o
no en cada fila del serpentín, dichos modelos son:

Modelamiento matemático para unidad manejadora enfriada por agua

Modelamiento matematico para unidad manejadora enfriada por refrigerante

60
Figura 18: Condición curva del serpentín

Fuente: Handbook de Refrigeración Industrial de Stoecker Wilbert. Pág.202

2.1 MODELAMIENTO MATEMATICO PARA UNIDAD MANEJADORA

ENFRIADO POR AGUA

2.1.1 Parámetros conocidos

Conocidas las condiciones de entrada de los flujos: para el agua, temperatura[K]

y los galones por minuto; para el aire, presión [kPa], temperatura de bulbo seco[K],
temperatura de bulbo húmedo[K] y los pies cúbicos por minuto; como también en
cuanto al intercambiador, material y geometría- Ver figura 19-(diámetro interno,
diámetro externo, número de tubos por fila, número de filas, paso transversal,
paso longitudinal, longitud del tubo, número de aletas y espesor de aletas) se
obtiene como función de lo anterior, los calores sensible, latente, total y las
condiciones de salida de los flujos.

61
 Figura 19: Geometría de unidad manejadora enfriada por agua

Fuente: Autores

2.1.2 Análisis del fenómeno de transferencia de masa y calor

En las unidades manejadoras de aire acondicionado es muy usual que sobre su

superficie se deposite algo del contenido de humedad asociada a él, es por tanto
necesario considerar el fenómeno de deshumidificación dentro del análisis del
modelo matemático. Entonces, con las condiciones de entrada conocidas para los
fluidos: agua y aire húmedo, junto con la geometría de la unidad y teniendo
presente que la configuración para este tipo de serpentines con proceso de
enfriamiento y deshumidificación es de flujo cruzado se plantea el siguiente
análisis:

Siendo la distribución de agua en una fila vertical y el comportamiento de un tubo

similar al resto de los otros tubos, se plantea dividir el intercambiador por filas
verticales (Ver figura 20) y evaluar para cada fila las condiciones de salida de los

62
fluidos de trabajo y los calores: total, el latente y el sensible. Donde cada porción
de tubo está vinculada a la porción adyacente por la condiciones de frontera, que
para este caso son las propiedades, ya que las condiciones de salida de un tramo
me representan la entrada del otro pues debe existir continuidad en las
condiciones de los flujos.

El método de solución propuesto se fundamenta en el balance de energía y

transferencia de masa tomando como volumen de control una fila del serpentín, en
donde interviene una rata de transferencia de calor convectiva interna desde
superficie exterior del tubo hacia el agua de enfriamiento que debe ser igual a la
rata de transferencia de calor convectiva externa desde el aire hacia la superficie
exterior del tubo más la rata de transferencia de masa difusiva desde la superficie
exterior aleteada hacia el aire -Ver figura 21. El modo de solución planteado está
condensado globalmente en un algoritmo y detalladamente en un diagrama de
flujo, los cuales se muestran a continuación:

Figura 20: Evolución del aire al pasar por cada fila del serpentín de UMA
enfriada por agua

Fuente: Autores

63
Figura 21: Balance de energía en una aleta circular

Fuente: Autores

64
2.1.2.1 Algoritmo global del modelamiento de la unidad manejadora enfriada
por agua

65
2.1.2.2 Diagrama de flujo del modelamiento para evaporador enfriado por
agua En el diagrama de flujo que se muestra a continuación, aparece un
compilado hecho en REFPROP para calcular las propiedades tanto del agua,
como del aire, por ser mucho más preciso que el resultado de la interpolación
lineal que se podría efectuar teniendo en cuenta las tablas para dichas sustancias.

También para la herramienta computacional elaborada se hizo uso de un

programa en Matlab disponible gratuitamente en Internet que genera una tabla con
un máximo de cinco columnas e ilimitado número de filas , adicionalmente posee
la barra de herramientas, para posibilitarle la opción de adquirir los resultados y
almacenarlos para análisis posteriores evitando incurrir en posibles errores a la
hora de transcribir los datos arrojados de la herramienta computacional, de igual
modo, se uso como un instrumento para mostrar los resultados de la evolución del
aire a medida que recorre las filas del serpentín en nuestra herramienta
computacional 18.

18 Tomado de la página de MathWorks, programa titulado mltable cuyo autor es: Morris Maynard, y está
basado en el código de Gregory Gershanok. 2005

66
67
68
69
2.2 MODELAMIENTO MATEMATICO PARA EVAPORADOR ENFRIADO POR
REFRIGERANTE

2.2.1 Parámetros conocidos.

Conocidas las condiciones de entrada de los flujos: para el refrigerante,

temperatura de entrada del evaporador y a la presión de alta de descarga del
compresor[K]; para el aire, presión [kPa], temperatura de bulbo seco[K],
temperatura de bulbo húmedo[K] y los pies cúbicos por minuto; como también en

70
cuanto al intercambiador, material y geometría(diámetro interno, diámetro externo,
número de tubos por fila, número de filas, paso transversal, paso longitudinal,
longitud del tubo, número de aletas y espesor de aletas) mostrados en la figura 22,
se obtiene como función de lo anterior, los calores sensible, latente, total y las
condiciones de salida de los flujos.

Figura 22: Geometría de unidad manejadora enfriada por refrigerante

Fuentes: Autores

2.2.2 Análisis del fenómeno de transferencia de masa y calor

En las unidades manejadoras de aire acondicionado es muy usual que sobre su

superficie se deposite algo del contenido de humedad asociada a él, es por tanto
necesario considerar el fenómeno de deshumidificación dentro del análisis del
modelo matemático. Entonces, con las condiciones de entrada conocidas para los
fluidos: agua y aire húmedo, junto con la geometría de la unidad y teniendo
presente que la configuración para este tipo de serpentines con proceso de

71
enfriamiento y deshumidificación es de flujo cruzado se plantea el siguiente
análisis:

Siendo la distribución del refrigerante en una fila vertical y el comportamiento de

un tubo similar al resto de los otros tubos, se plantea dividir el intercambiador por
filas verticales (Ver figura 23) y evaluar para cada fila las condiciones de salida de
los fluidos de trabajo y los calores: total, el latente y el sensible. Donde cada
porción de tubo está vinculada a la porción adyacente por la condiciones de
frontera, que para este caso son las propiedades, ya que las condiciones de salida
de un tramo me representan la entrada del otro pues debe existir continuidad en
las condiciones de los flujos.

Figura 23: Evolución del aire al pasar por cada fila del serpentín de UMA
enfriada por agua

Fuente: Autores

72
El método de solución propuesto se fundamenta en el balance de energía y
transferencia de masa tomando como volumen de control una fila del serpentín, en
donde interviene una rata de transferencia de calor evaporativa-convectiva interna
desde superficie exterior del tubo hacia el refrigera-nte que debe ser igual a la rata
de transferencia de calor convectiva externa desde el aire hacia la superficie
exterior del tubo más la rata de transferencia de masa difusiva desde la superficie
exterior aleteada hacia el aire -Ver figura 21. El modo de solución planteado está
condensado globalmente en un algoritmo y detalladamente en un diagrama de
flujo, los cuales se muestran a continuación:

2.2.2.1 Algoritmo global del modelamiento de la unidad manejadora enfriada

por refrigerante

73
2.2.2.2 Algoritmo del modelamiento para evaporador enfriado por
refrigerante. Para desarrollar el siguiente proceso de cálculo se tuvo en cuenta
las siguientes consideraciones:

Para el Número de Bo, se presumió la condición de salida del aire húmedo para
establecer la rata de calor que se requiere para calcular este número, necesario
para hallar el coeficiente de transferencia de calor de evaporación con turbulencia.
Dado que la recomendación para la este calor oscila entre 5-70 W/m2.

De acuerdo al libro de Kakac de Boilers, Evaporators and Condenser, la calidad

del refrigerante para efectos de cálculos de transferencia se recomienda tomarla
como la media, es decir 0.5.

Adicionalmente, se recuerda el uso del compilado REFPROP y el programa de la

tabla descargada en Internet mencionada en el anterior modelamiento.

74
75
76
77
78
79
 3. HERRAMIENTA COMPUTACIONAL PARA ESTUDIAR LAS UNIDADES
MANEJADORAS EN AIRE ACONDICIONADO

La herramienta computacional inicia en la interfaz grafica de Matlab inicio. fig. Al

momento de dar run se despliega la siguiente interfaz.

Figura 24. Interfaz principal de la herramienta computacional

Fuente: Autores

Esta herramienta calcula analíticamente (basado en un modelo matemático de

transferencia de calor y mecánica fluidos) la carga sensible, carga latente, carga
total, temperatura y húmeda de salida del aire en una unidad enfriadora de tipo

80
banco aleteado. La fiabilidad de la herramienta se establece a través de la
comparación de los resultados con unidades ya constituidas, tales como:

Unidades enfriadoras por catalogos

Sistema Climatizador aire - agua Modelo YG de YORK
Unidad de Expansión Directa Fancoil Modelo R1MX de YORK

Unidad instalada experimental

Unidad Manejadora aire – agua Modelo UMEC LSL 114FCV1 UMA2 de
Paramo

Los datos a comparar son los resultados del cálculo analítico que elabora la
herramienta computacional los cuales son:

Temperatura del aire a la salida (TEM AIR OUT)

Carga sensible (THS)
Carga Latente (THL)
Carga Total (TH)
Humedad relativa (φ)

Con el fin de disminuir el trabajo de entrada de los datos físicos de las unidades
manejadoras a analizar se han introducido en la herramienta computacional las
unidades enfriadoras como:

a. Modelo YG Aire-Agua de York

b. Modelo R1MX de expansión directa

Para comparar los resultados y determinar el desempeño de las unidades

manejadoras lo realizaremos según la herramienta computacional siguiendo el
siguiente caso:

81
3.1 ANÁLISIS DE UNIDADES MANEJADORAS

El análisis para las unidades manejadoras puede ser por catálogos o

personalizado. En la interfaz que se encuentra a continuación se realizara el
análisis por Catálogos.

Figura 245 Selección del análisis de unidades manejadoras

Fuente: autores

Una vez de haber seleccionado el análisis de las unidades manejadoras y de

haber dado click en el botón siguiente se inicia una nueva interfaz con el fin de
seleccionar el medio refrigerante. A continuación se muestra la interfaz de
selección del medio refrigerante.

82
 Figura 256 Selección del medio refrigerante

Fuente. Autores

Realizada la selección del medio refrigerante se inicia el proceso de análisis de

unidades manejadoras por catalogo. Donde se debe seguir todos los pasos
establecido en este análisis. Para llegar a una respuesta correcta es necesario
seguir los pasos mostrado en la interfaz que se encuentra a continuación.

3.2 SELECCIÓN DE UNIDADES.

La selección de Unidades es representada en la interfaz a continuación.

83
Figura 267 Pasos a seguir para el análisis de las unidades manejadoras por
catálogos

Fuente: Autores

Al oprimir el botón de selección de unidades se despliega una nueva interfaz para

seleccionar unidades enfriadoras por catalogos. En la interfaz de selección es
necesario llevar unos pasos para realizar una correcta selección de unidad
manejadora. El primer paso es seleccionar el modelo que se quiere. En la interfaz
a continuación es representada.

84
 Figura 278 Selección de unidad enfriadora

Fuente. Autores

El segundo paso para realizar la correcta selección de unidades manejadora es

escoger el número de filas en la unidad que se quiere analizar. Este procedimiento
de seleccionar el número de filas es representado en la interfaz a continuación.

85
 Figura 289 selección del numero de filas

Fuente. Autores

Una vez hecho la selección del numero de filas en la unidad manejadora es

necesario de oprimir el botón seleccionar para poder garantizar una correcta
selección. Este procedimiento es ilustrado en la siguiente interfaz.

86
 Figura 30 seleccionar unidad

Fuente. Autores

Al terminar la selección de la unidad y se desea conocer los datos geométricos

que tiene el serpentin de la unidad manejadora se puede oprimir botón geometría,
y asi conocer los datos geométricos. A continuación se ilustra este procedimiento.

87
Figura 31 Ilustración para ver la geometría

Fuente: Autores

Figura 329 Geometria de la unidad

Fuente: Autores

88
Una vez visto la geometría de la unidad manejadora, es necesario oprimir el botón
atrás, el cual este botón se encuentra en la interfaz de geometría del serpentín. A
continuación se muestra la ilustración.

Figura 33 Ilustración para continuar en el proceso

Fuente. Autores

Finalizado el chequeo de la geometría de la unidad manejadora, vuelve a la

interfaz principal (selección de unidad enfriadora por catalogo), el cual es
necesario oprimir el botón aceptar para que se realiza una adecuada selección.
En la interfaz a continuación es ilustrado.

89
 Figura 34 Selección de la unidad enfriadora

Fuente Autores

Nuevamente regresamos a los pasos para el proceso de análisis de unidades

manejadoras por catálogos.

Hecha una correcta selección se requiere oprimir el boton evaluar unidad para
seguir avananzando en el analisis de unidades manejadoras. A continuacion se
representa el procedimiento correcto.

90
 Figura 35. Pasos a seguir para el análisis de las unidades manejadoras por
catálogos 2

Fuente: Autores

Un procedimiento realizado por el programa es válida la selección de la unidad

manejadora ya hecha. Por eso, es necesario realizar una correcta selección, sin
interrumpir el procedimiento.

3.3 COMPROBACION DE LA UNIDAD MANEJADORA

En el segundo paso que es la comprobación de la unidad manejadora se realiza

los cálculos analíticos hecho por el programa. Al oprimir el botón de comprobación

91
de la unidad manejadora que se muestra en la interfaz siguiente se abre una
nueva ventana.

Figura 36. Ilustración para la comprobación de la unidad manejadora

Fuente: Autores

En esta fase se comprueba la unidad manejadora seleccionada del catalogo con

los cálculos analíticos hechos por la herramienta computacional.

En la interfaz que se muestra a continuación no es posible alterar los datos de

entrada de la unidad, los cuales fueron seleccionado en el catalogo.

92
 Figura 37. Fase de comprobación para la unidad seleccionada

Fuente: Autores

Dentro de la interfaz de comprobación de la unidad manejadora por catalogo se

encuentran consignados los datos de catalogo, tales como:

Temperatura del aire a la salida: T air out

Calor sensible: TSH
Calor latente:TLH
Temperatura del agua a la salida de la unidad. T wat out
Humedad relativa del aire a la salida: HR
Numero de filas de la unidad seleccionada: NF

93
Todos estos datos son comparados con el calculo analítico que realiza el
programa. A continuación se mostrara las interfaces que realiza la herramienta
para realizar el cálculo.

Figura 38 Fase de comprobación para la unidad seleccionada 2

Fuente: Autores

94
 Figura 39 Fase de espera para la comprobación

Fuente: Autores

Una vez que la herramienta termine de hacer los cálculos, los datos analíticos son
colocados en la parte de resultados de la unidad con el fin de compararlos.

Es necesario recordar que los cálculos fueron hechos para la unidad que ya fue
seleccionado. Un caso particular de un ejemplo es mostrado a continuación. los
2cálculos fueron hecho para la unidad YG York modelo 2.

95
 Figura 40. Fase de comprobación para la unidad selecciona 3

Fuente: Autores

Cuando los cálculos hechos por la herramienta computacional terminen, se activa

un botón de evolución del aire mostrado en la siguiente interfaz.

96
 Figura 41. Fase de comprobación para la unidad selecciona 4

Fuente: Autores

El botón que se llama tabla de evolución del aire y se represento en la interfaz

anterior tiene como objetivo mostrar la evolución del aire a través del numero del
filas de la unidad manejadora seleccionada. Los datos que ofrece este botón son.

Calor sensible: Qs(Kw)

Calor latente= Ql(Kw)
Humedad relativa: Humedad (%)
Temperatura de salida del aire: T aire (°C)

En la parte de los calores latente y sensible la tabla muestra la suma total de los
calores al final de cada uno.

97
 Figura 42 Ilustración de la tabla de evolución del aire

Fuente: Autores

Una de las ventajas que tiene esta tabla o matriz es que tiene la opción de guardar
el reporte de los cálculos que se realizo en la unidad manejadora. Esto se
encuentra en la barra de tarea donde dice File, y ahí se tiene la opción para
guardar la tabla. A continuación se ilustra este procedimiento.

98
 Figura 43 Ilustración para guardar la tabla de evolución del aire

Fuente: Autores

Ya terminado el proceso de comprobación es necesario de oprimir el botón de

aceptar para la ejecución de la herramienta y entrar a la última fase que es
calcular.

99
 Figura 44 Fase de comprobación para la unidad selecciona 5

Fuente: Autores

Al oprimir el botón aceptar nuevamente volvemos a los pasos principales del

análisis de unidades manejadoras por catálogos. Para avanzar en este último
paso en necesario de oprimir el botón evaluar unidad. A continuación se visualiza
el ejemplo.

100
 Figura 45. Fase de proceso de análisis de unidades para calcular

Fuente: Autores

Una vez de oprimir el botón evaluar unidad entramos al siguiente paso, el cual es
calcular.

101
3.4 CALCULO PERSONALIZADO DE LA UNIDAD MANEJADORA

Figura 46 Fase de proceso de análisis de unidades para calcular 2

Fuente: Autores

En la fase final del proceso de análisis de unidades manejadoras por catálogos

Entramos a la interfaz calcular unidad manejadora que se muestra a continuación

y l nos brinda la forma de modificar las condiciones de entrada de la unidad. Y así
conocer la cargas de la unidad.

102
 Figura 47 Cálculo de la unidad manejadora

Fuente: Autores

Un dato adicional del proceso de calculo es que la herramienta grafica el proceso

psicométrico a las condiciones que se desean. A continuación se ilustra el proceso
psicométrico hecho por la herramienta.

103
 Figura 48. Ilustración del diagrama psicométrico en la unidad a calcular

Fuente: Autores

Después de haber realizado los cálculos es necesario cerrar el programa y si se

desea analizar otro modelo volver al primer paso que es selección de la unidad
manejadora

104
 4. ANÁLISIS DE LOS RESULTADOS DE LA HERRAMIENTA
COMPUTACIONAL

4.1 CORROBORACIÓN DE LOS RESULTADOS DE LA HERRAMIENTA CON

OTRAS FUENTES

Para los resultados obtenidos de la herramienta computacional, primero se realiza

su comparación contra los extraídos de las siguientes fuentes, tal como sigue:

Para enfriamiento por agua:

Catálogo Modelo YGH de York19

Mediciones tomadas de la UMA en el Edificio de la Facultad de Ciencias
Humanas- Primer piso Sala de Computo

Para enfriamiento por refrigerante:

Catálogo RMX de York20

4.1.1 Comparación de resultados del programa Vs. Catálogo Modelo YGH de

York para el Agua

Para corroborar los resultados de la herramienta se extrajeron varias tablas de

comparación con los resultados del catalogo como sigue a continuación:

19
Ver Anexo A. MANUAL TECNICO DE COMFORT PACKAGE FAN & COIL UNIT- YG
CAPACIDADES DE 1360 A 10200 m3/h
20
Ver Anexo B. CATALOGO DE YORK UNIDADES FANCOILS DE EXPANSIÓN DIRECTA
CAPACIDADES DE 12000 A 36000 BTU/H

105
4.1.1.1 Comparación de la carga sensible entre el catálogo de la UMA para
agua de enfriamiento vs. la herramienta computacional Para la comparación
se selecciono del catálogo del fabricante los modelos 02, 03, 09 y cada número de
filas (3, 4, 6 y 8 filas) para elaborar las tablas de comparaciones, teniendo en
cuenta las condiciones de entrada:

Para el aire húmedo:

Temperatura de bulbo seco a la entrada Tbs=26.7°C

Temperatura de bulbo húmedo a la entrada Tbh=19.4 °C
Presión atmosférica=101.325kPa

Para el agua:

Temperatura del agua= 6,7 ºC

GPM @ modelo (Ver anexo A)

Las tablas de comparación son las siguientes:

Tabla 1. Comparación de carga sensible del catálogo modelo 02 Vs Carga

sensible del programa
CATALOGOS YG YORK 02
NF total TSH Programa TSH
Error (%)
(filas) (Kw) Catalogo(Kw)
3 3.07716 3.706 16.97
4 4.39924 4.40613 0.16
6 5.59347 6.1694 9.34
8 6.82085 6.9167 1.39
ERROR PROMEDIO(%) 6.96139421
Fuente: Autores

106
 Figura 49. Geometría Modelo 2 YG del catálogo de York 21

Fuente: Autores

Tabla 2. Comparación de carga sensible del catalogo modelo 03 Vs Carga

sensible del programa
CATALOGOS YG YORK 03
NF total TSH Programa TSH
Error (%)
(filas) (Kw) Catalogo(Kw)
3 1.67494 1.8111 7.52
4 5.089 7.1679 29.00
6 9.4774 9.3498 1.36
8 9.47747 10.3606 8.52
ERROR PROMEDIO(%) 11.6024152
Fuente: Autores

21
La geometría entre el mismo modelo pero distinto número de filas, es igual a excepción de este parámetro
geométrico.

107
 Figura 50. Geometría Modelo 3 YG del catálogo de York 22

Fuente: Autores

Tabla 3. Comparación de carga sensible del catálogo modelo 09 Vs Carga

sensible del programa
CATALOGOS YG YORK 09
NF total TSH Programa TSH
Error (%)
(filas) (Kw) Catalogo(Kw)
3 16.5717 17.2842 4.12
4 19.8661 21.488 7.55
6 24.94 26.96 7.49
8 16.59 17.6063 5.77
ERROR PROMEDIO(%) 6.23378571
Fuente: Autores

22
La geometría entre el mismo modelo pero distinto número de filas, es igual a excepción de este parámetro
geométrico.

108
 Figura 51. Geometría Modelo 9 YG del catálogo de York 23

Fuente: Autores

4.1.1.2 Comparación de la carga latente entre el catálogo de la UMA para

agua de enfriamiento vs. la herramienta computacional. Para la comparación
se selecciono del catálogo del fabricante los modelos 02, 03, 07 y cada número de
filas (3, 4, 6 y 8 filas) para elaborar las tablas de comparaciones, teniendo en
cuenta las condiciones de entrada:

Para el aire húmedo:

Temperatura de bulbo seco a la entrada Tbs=26.7°C

Temperatura de bulbo húmedo a la entrada Tbh=19.4 °C

23
La geometría entre el mismo modelo pero distinto número de filas, es igual a excepción de este parámetro
geométrico.

109
 Presión atmosférica=101.325kPa

Para el agua:

Temperatura del agua= 6,7 ºC

GPM @ modelo (Ver anexo A)

Las tablas de comparación son las siguientes:

Tabla 4. Comparación de carga latente del catalogo modelo 02 Vs Carga

latente del programa24
CATALOGOS YG YORK 02
NF TLH TLH
Modelo Error (%)
(filas) Catalogo(Kw) Programa(Kw)

3 2 0.59 0.569246 3.51762712

4 3 1.8772 1.52355 18.8392286
6 4 2.2307 1.99986 10.3483212
8 5 3.4 3.793 11.5588235
ERROR PROMEDIO (%) 11.0660001
Fuente: Autores

24
Ver datos de geometría del modelo 2 en la figura n.

110
 Tabla 5. Comparación de carga latente del catalogo modelo 03 Vs Carga
latente del programa25

CATALOGOS YG YORK 03

NF TLH TLH
Modelo Error (%)
(filas) Catalogo(Kw) Programa(Kw)

3 2 5.771 4.76394 17.4503552

4 3 2.4831 3.2336 30.2243164
6 4 3.4773 3.95344 13.6928076
8 5 3.81 3.95344 3.7648294
ERROR PROMEDIO (%) 16.2830772
Fuente: Autores

Tabla 6. Comparación de carga latente del catalogo modelo 07 Vs carga

latente del programa

CATALOGOS YG YORK 07

NF TLH TLH Error

Modelo
(filas) Catalogo(Kw) Programa(Kw) (%)

3 2 5.89 10.19 4.2

4 3 2.4831 3.2336 17.97
6 4 4.29 4.341 1.16
8 5 3.81 5.445 17.95
ERROR PROMEDIO (%) 10.32
Fuente: Autores

25
Ver datos de geometría del modelo 3 en la figura n.

111
 Figura 52. Geometría Modelo 7 YG del catálogo de York 26

Fuente: Autores

4.1.1.3 Comparación de temperatura de bulbo seco del aire de salida entre el

catálogo de la UMA para agua de enfriamiento vs. la herramienta
computacional. Para la comparación se selecciono del catálogo del fabricante los
modelos 02, 03, 07 y cada número de filas (3, 4, 6 y 8 filas) para elaborar las
tablas de comparaciones, teniendo en cuenta las condiciones de entrada:
Para el aire húmedo:

Temperatura de bulbo seco a la entrada Tbs=26.7°C

Temperatura de bulbo húmedo a la entrada Tbh=19.4 °C
Presión atmosférica=101.325kPa

26
La geometría entre el mismo modelo pero distinto número de filas, es igual a excepción de este parámetro
geométrico.

112
Para el agua:

Temperatura del agua= 6,7 ºC

GPM @ modelo (Ver anexo A)

Las tablas de comparación son las siguientes:

Tabla 7. Comparación de temperatura de bulbo seco del aire de salida del

catalogo modelo 02 Vs temperatura de bulbo seco del aire de salida del
programa27

CATALOGOS YG YORK 02
NF
Modelo T Catalogo (ºC) T programa (ºC) Error (%)
(filas)
3 2 18.3 18.9 3.278668
4 3 15.6 16.1 3.205128
6 4 12.8 14.1 10.15625
8 5 11 12.5 13.63636
ERROR PROMEDIO (%) 7.56907
Fuente: Autores

Tabla 8. Comparación de temperatura de bulbo seco del aire de salida del

catalogo modelo 03 Vs temperatura de bulbo seco del aire de salida del
programa28
CATALOGOS YG YORK 03
NF T Catalogo T programa
Modelo Error (%)
(filas) (ºC) (ºC)
3 2 18 20.5178 13.9877778
4 3 15.9 16.18 1.76100628
6 4 12.6 13.3693 6.10555556
8 5 11.1 12.1 9.00900900
ERROR PROMEDIO(%) 7.7158370
Fuente: Autores

27
Ver datos de geometría del modelo 2 en la figura n.
28
Ver datos de geometría del modelo 3 en la figura n.

113
 Tabla 9. Comparación de temperatura de bulbo seco del aire de salida del
catalogo modelo 07 Vs temperatura de bulbo seco del aire de salida del
programa29
CATALOGOS YG YORK 07
NF
Modelo T Catalogo (ºC) T programa (ºC) Error (%)
(filas)
3 2 17.9 19.91 9.35
4 3 15.9 16.3 2.515
6 4 12.6 13.6 7.936
8 5 11.1 14.608 9.48
ERROR PROMEDIO(%) 7.3203
Fuente: Autores

4.1.1.4 Comparación de la humedad relativa del aire salida entre el catálogo

de la UMA para agua de enfriamiento vs. la herramienta computacional. Para
la comparación se selecciono del catálogo del fabricante los modelos 12 y 07 y
cada número de filas (3, 4, 6 y 8 filas) para elaborar las tablas de comparaciones,
teniendo en cuenta las condiciones de entrada:

Para el aire húmedo:

Temperatura de bulbo seco a la entrada Tbs=26.7°C

Temperatura de bulbo húmedo a la entrada Tbh=19.4 °C
Presión atmosférica=101.325kPa

Para el agua:

Temperatura del agua= 6,7 ºC

GPM @ modelo (Ver anexo A)

29
Ver datos de geometría del modelo 7 en la figura n.

114
Las tablas de comparación son las siguientes:

Tabla 10. Comparación de humedad relativa del aire de salida del catálogo
modelo 12 Vs humedad relativa del aire de salida del programa
CATALOGOS YG YORK 12
Humedad
NF total Humedad relativa
relativa Error (%)
(filas) Catalogo (%)
Programa (%)
3 66,7691 70,5 7.32
4 70,9246 81 10.68
6 77,9352 90 26.95
8 83,5238 96 23.12
ERROR PROMEDIO(%) 17.02
Fuente: Autores

Figura 53. Geometría Modelo 7 YG del catálogo de York 30

Fuente: Autores

30
La geometría entre el mismo modelo pero distinto número de filas, es igual a excepción de este parámetro
geométrico.

115
 Tabla 11. Comparación de humedad relativa del aire de salida del catalogo
modelo 07 Vs humedad relativa del aire de salida del programa31
CATALOGOS YG YORK 7
Humedad Humedad
NF total
relativa relativa Catalogo Error (%)
(filas)
Programa (%) (%)
3 66,2579 70 7.32
4 70,4144 79 10.68
6 77,4776 90 26.95
8 83,0507 96 23.12
ERROR PROMEDIO(%) 17.02
Fuente: Autores

4.1.2 Mediciones tomadas de la UMA en el Edificio de la Facultad de Ciencias

Humanas- Primer piso Sala de Computo.

Figura 54: Unidad Manejadora de Aire marca Paramo

Fuente: Autores

31
Ver datos de geometría del modelo 7 en la figura n.

116
Unidad en funcionamiento que se encuentra ubicada en la sala de cómputo del
edificio de ciencias humanas de la Universidad Industrial de Santander. La cual
presenta las siguientes características:

Modelo UMEC LSL 114FCV1 UMA2

Serie 2006A5032-91344
Motor 5 HP
Frecuencia 60 Hz

Geometría de la Unidad
Número de Filas= 4
Número de Tubos por Fila=20
Número de aletas/pulgada =15
Diámetro externo tubo=17 mm
Diámetro interno tubo= 15.5 mm
Tubos de cobre escalonados
Largo=1,76 m
Alto= 86 cm
Aletas de aluminio corrugadas
Profundidad= 12,2 cm
Paso transversal= 31 mm
Paso longitudinal= 35 mm
Espesor de Aleta= 0.15 mm
Número de Aletas totales= 1039

Para calcular las cargas del recinto y las cargas del aire exterior es necesario usar
las siguientes fórmulas:

117
 Temperatura del recinto a condicionar (°C)
Temperatura del aire impulsado en la unidad manejadora (°C)
Temperatura de Bucaramanga (°C)
Humedad absoluta del recinto a condicionar (gv/Kgas)
Humedad absoluta del aire impulsado por la unidad (gv/Kgas)
Humedad absoluta de Bucaramanga (gv/Kgas)
Calo sensible impuesto por el recinto
Calor sensible impuesto por el aire externo
Calor latente impuesto por el aire externo
Calor latente que hay en el recinto
Calor sensible total a extraer la unidad manejadora
Calor latente total a extraer la unidad manejadora
Calor total que hay en el recinto

El flujo de agua y la masa de condensado se obtiene mediante la ecuación:

Instrumentos con los que se tomaron las pruebas fueron:

118
 Medidor de temperatura y humedad

Modelo: Tecpel 320 Series

Mide: Temperatura & Humedad
Referencia. RS-232 Test link

Figura 5530: Medidor de T y

Fuente: Autores

Pistola infrarrojo
Mide: temperatura del agua en la unidad
Ajustada con Emisividad de 0.9
Marca. Fluke

Figura 56: Pistola Infrarrojo

Fuente: Autores

119
 Termo Anemómetro
Mide: velocidad del aire (m/s)
Modelo: ERASMUS EA-200

Figura 31: Anemómetro

Fuente: Autores

Duración de cada prueba: 30 minutos

Numero de prueba/día=5
Días de prueba=3

Para la comparación de los datos experimentales con los resultados de los

cálculos analíticos de la herramienta computacional es necesario tomar los datos
experimentales medidos con los instrumentos de medición (tablas 9, 11 y 14) de
cada prueba realizada según el día y introducirlo en la herramienta computacional
con el fin de evaluar los datos de la unidad experimental y compararlos.

A continuación se muestra cada prueba, el cual es correspondiente a cada día.

Muestra los datos tomados con los instrumentos de medición y fueron tabulados
en tablas, con el fin de comparar los resultados con la herramienta computacional.

120
PRUEBA DIA 1
Fecha: Abril 12 de 2010 Hora Inicio: 9:00 am Hora Final: 10:00 am
Tiempo/Dato=15 min

Tabla12. Primera tabla de datos experimentales

Aire Exterior de
Aire Entrada en UMA Aire Impulsión en UMA Aire en el recinto CFM
Bucaramanga
T φ T T φ T φ ω
(°C) (%) ω (g/Kg.as) (°C) φ (%) Tbh (°C) (°C) (%) ω (g/Kg.as) (°C) (%) (g/Kg.as) Imp Aire Ext
DATO 1 27.5 61 16.55 22 55 18.4 19.8 70 8.9 23.8 56 12.1 4000 441
DATO 2 27.5 61 16.55 20 56 18 12 76 10.5 15 57 14 4000 441
DATO 3 27.5 61 16.55 20 63 18 12 80 9 15 53 13 4000 441
DATO 4 27 60.5 15.92 20 64 18 13 82 9.5 16 60 13.5 4000 441
DATO 5 26.4 63 15.98 18.5 83 16 10 74 9.5 14 64 12.1 4000 441
Fuente: Autores

Tabla 13. Cálculos de los datos experimentales

RSH (Kw) RLH (Kw) OASH (Kw) OALH (Kw) TSH (Kw) TLH (Kw) TH (Kw) Gpm
DATO 1 20.079 10.560 2.0431 1.6191 22.1224 12.1793 34.302 20
DATO 2 15.060 11.550 6.9023 0.9278 21.9618 12.4781 34.440 20
DATO 3 15.060 13.200 6.9023 1.2916 21.9618 14.4919 36.454 17
DATO 4 15.060 13.200 6.0740 0.8805 21.1335 14.0808 35.214 15
DATO 5 20.079 8.580 6.8471 1.4117 26.9264 9.9919 36.918 20
Fuente: Autores

121
PRUEBA DIA 2
Fecha: Abril 14 de 2010 Hora Inicio: 2:00pm Hora Final: 4:30 pm
Tiempo/dato: 30 min

Tabla 14. Primera tabla de datos experimentales

Aire Exterior de Aire Entrada en Aire Impulsion en
Aire en el recinto CFM
Bucaramanga UMA UMA
T φ ω T φ Tbh T φ ω T φ ω Aire
(°C) (%) (g/Kg.as) (°C) (%) (°C) (°C) (%) (g/Kg.as) (°C) (%) (g/Kg.as) Impulsion Externo
DATO 6 27 61.2 16.11 21.8 88 20 15 74 8 20 55.5 14 4000 441
DATO 7 28 61 17.06 18 67.1 15 20 76.8 10.54 22.8 60 12.18 4000 441
DATO 8 27 61.3 16.13 21.8 58.5 16.3 20.1 76.9 11.05 23.4 63 13.27 4000 441
DATO 9 27 61 16.66 20 58 16.4 20.4 79.4 10.85 23.9 59 13.5 4000 441
DATO 10 27 61 16.05 19 63.4 17.3 19.5 80.1 10.81 23.1 63.4 14.3 4000 441
Fuente: Autores

Tabla 15. Cálculos de los datos experimentales

RSH (Kw) RLH (Kw) OASH (Kw) OALH (Kw) TSH (Kw) TLH (Kw) TH (Kw)
DATO 6 25.0992 19.80048 3.8652768 0.76768936 28.9644768 20.5681694 49.5326462
DATO 7 14.055552 5.4121312 2.87134848 1.77550904 16.9269005 7.18764024 24.1145407
DATO 8 16.565472 7.3261776 1.98785664 1.04056473 18.5533286 8.36674233 26.920071
DATO 9 17.56944 8.745212 1.71176544 1.14971487 19.2812054 9.89492687 29.1761323
DATO 10 18.071424 11.5172792 2.15351136 0.63670919 20.2249354 12.1539884 32.3789237
Fuente. Autores

122
Comparando las condiciones de salida del programa con las mediciones el
catálogo se obtuvo lo siguiente:

Para la prueba 1:
Se comparo con el programa utilizando la Presión atmosférica de
Bucaramanga=87.83 Kpa
Con sus respectivos datos a:
Tabla de datos experimentales prueba día 1, condiciones de entrada de aire a la
unidad manejadora. Ver tabla de datos experimentales.

Tabla 16. Tablas de comparación de carga latente de la unidad Paramo Vs

Carga latente del programa
TLH TLH
Toma de NF total
programa unidad Error (%)
Datos (filas)
(Kw) (Kw)
DATO 1 4 13,34 12,179 10.24
DATO 2 4 14.26 12,478 14.28
DATO 3 4 14.95 14,491 3.07
DATO 4 4 14.55 14,080 3.33
DATO 5 4 13.17 9,991 31.81
Error Promedio 12.54
Fuente. Autores

Tabla 17. Tablas de comparación de carga total de la unidad Paramo Vs

Carga total del programa
TH TH
Toma de NF
programa unidad Error (%)
Datos humedas
(Kw) (Kw)
DATO 1 4 30.83416 34,302 10.11
DATO 2 4 37.72 34,440 9.52
DATO 3 4 36,415 36,454 0.10
DATO 4 4 35.54 35,214 0.92
DATO 5 4 31.72 36,918 14.07
Error Promedio 6.94
Fuente. Autores

123
Tabla 18. Tablas de comparación para la temperatura de la unidad Paramo Vs
temperatura de salidad del programa

T T
Toma de NF total
programa unidad Error (%)
Datos (filas)
(°C)
DATO 1 4 9.15 19.8 53.70
DATO 2 4 9.21 12 23.20
DATO 3 4 8.98 12 25.10
DATO 4 4 9.2 13 31.60
DATO 5 4 9.05 10 9.50
Error Promedio 28.6
Fuente. Autores

Tabla 19. Tablas de comparación para la humded relativa de la unidad

Paramo Vs temperatura de salidad del programa
Humedad
Humedad
Toma de NF total relativa
relativa Error (%)
Datos (filas) programa
unidad (%)
(%)
DATO 1 4 100 70 42.80
DATO 2 4 100 76 31.50
DATO 3 4 100 80 25.00
DATO 4 4 100 82 21.90
DATO 5 4 100 74 35.10
Error Promedio 31.26
Fuente. Autores

Para la prueba 2:

Se comparo con el programa utilizando la Presión atmosférica de

Bucaramanga=87.83 Kpa

124
Con sus respectivos datos a:
Tabla de datos experimentales prueba día 2, condiciones de entrada de aire a la
unidad manejadora. Ver tabla de datos experimentales.

Tabla 20. Tablas de comparación de carga latente de la unidad Paramo Vs

Carga latente del programa
TLH TLH
Toma de NF total
programa unidad Error (%)
Datos (filas)
(Kw) (Kw)
DATO 6 4 21.06 20,568 2.39
DATO 7 4 6.66 7,187 7.33
DATO 8 4 7.11 8,366 15.01
DATO 9 4 9.47 9,894 4.28
DATO 10 4 13.74 12,153 13.05
Error Promedio (%) 8,412
Fuente. Autores

Tabla 21. Tablas de comparación de carga total de la unidad Paramo Vs

Carga total del programa
TH TH
Toma de NF
programa unidad Error (%)
Datos humedas
(Kw) (Kw)
DATO 6 4 47.87 49,532 3.35
DATO 7 4 21,894 24,114 9.20
DATO 8 4 29 26,920 7.72
DATO 9 4 29.09 29,176 0.29
DATO 10 4 33.32 32,378 2.90
Error Promedio 4,692
Fuente. Autores

125
Tabla 22. Tablas de comparación para la temperatura de la unidad Paramo Vs
temperatura de salidad del programa
T T
Toma de NF total
programa unidad Error (%)
Datos (filas)
(°C)
DATO 6 4 10.5 15 30.00
DATO 7 4 16.2 20 19.00
DATO 8 4 17.4 20.1 13.40
DATO 9 4 18 20.4 11.70
DATO 10 4 14.8 19.5 24.10
Error Promedio 19.64
Fuente. Autores

Tabla 23. Tablas de comparación para la humedad relativa de la unidad

Paramo Vs humedad relativa de salida del programa
Humedad
Humedad
Toma de NF total relativa
relativa Error (%)
Datos (filas) programa
unidad (%)
(%)
DATO 6 4 100 74 35.10
DATO 7 4 100 76.8 30.20
DATO 8 4 100 76.9 30.00
DATO 9 4 100 79.4 25.90
DATO 10 4 100 80.1 24.80
Error Promedio 29.2
Fuente. Autores

4.1.3 Comparación de resultados del programa Vs. Catálogo Modelo RMX de

York para el Refrigerante R-22

Para corroborar los resultados de la herramienta se extrajeron varias tablas de

comparación con los resultados del catalogo como sigue a continuación:

4.1.3.1 Comparación de la humedad relativa del aire salida entre el catálogo

de la UMA para vs. la herramienta computacional para refrigerante R-22 Las

126
 comparaciones se hicieron entre los datos ofrecidos por los fabricantes y los
obtenidos por el programa. Los datos obtenidos por el programa se realizo a dos
temperaturas de evaporación a 4.44°C y 7.22°C con los siguientes datos:

Temperatura de bulbo seco a la entrada Tbs=23.5°C

Temperatura de bulbo húmedo a la entrada Tbh=16.67 °C
Se asumió una temperatura de subenfriamento=2°C
Se asumió una temperatura de alta=35°C

Tabla 24. Comparación de carga latente del catalogo Vs Carga latente del
programa para Tev=4.44 ºC
CATALOGO R1MX - YORK R22
T CFM NF TLH TLH
Modelo evaporización entrada total catalogo programa Error (%)
R22 (ºC) (filas) (Kw) (Kw)
12 4.44 500 2 0.76 0.727 4.34210526
18 4.44 600 3 0.78 0.757 2.94871795
24 4.44 800 3 1.04 1.06 1.92307692
30 4.44 1000 3 1.32 1.29 2.27272727
36 4.44 1200 3 1.3 1.39 6.92307692
Error Promedio 3.68194087
Fuente. Autores

Tabla 25. Comparación de carga latente del catalogo Vs Carga latente del
programa para Tev=7.22 ºC
CATALOGO R1MX - YORK R22
T CFM NF TLH TLH
Modelo evaporización entrada total catalogo programa Error (%)
R22 (ºC) (filas) (Kw) (Kw)
12 7.22 500 2 0.28 0.262 6.42857143
18 7.22 600 3 0.27 0.245 9.25925926
24 7.22 800 3 0.32 0.309 3.4375
30 7.22 1000 3 0.46 0.49 6.52173913
36 7.22 1200 3 0.39 0.34 12.8205128
Error Promedio 7.69351653
Fuente. Autores

127
 Tabla 26. Comparación de Temperatura de salida del catalogo Vs
Temperatura de salida del programa Tev=4,44 ºC

CATALOGO R1MX - YORK R22

T CFM T T
entrada NF total
Modelo evaporización Catalogo programa Error (%)
(filas)
R22 (ºC) (ºC) (ºC)
12 4.44 500 2 13.89 14.5 4.39164867
18 4.44 600 3 13.89 13.37 3.7437005
24 4.44 800 3 14.44 13.089 9.35595568
30 4.44 1000 3 14.44 12.8 11.3573407
36 4.44 1200 3 15 13.74 8.4
Error Promedio 7.44972911
Fuente. Autores

Tabla 27. Comparación de Temperatura de salida del catalogo Vs

Temperatura de salida del programa Tev=7,22 ºC
CATALOGO R1MX - YORK R22

T CFM T
entrada NF total T Catalogo
Modelo evaporización programa Error (%)
(filas) (ºC)
R22 (ºC) (ºC)

12 7.22 500 2 15 15.83 5.53333333

18 7.22 600 3 15.56 14.86 4.49871465
24 7.22 800 3 15.56 14.65 5.84832905
30 7.22 1000 3 15.56 14.47 7.00514139
36 7.22 1200 3 16.11 15.11 6.20732464
Error Promedio 5.81856861
Fuente. Autores

128
 Tabla 28. Comparación de carga sensible del catalogo Vs Carga sensible del
programa con Tev=4,44 ºC

CATALOGO R1MX - YORK R22

T CFM TSH TSH

entrada NF total
Modelo evaporización catalogo programa Error (%)
(filas)
R22 (ºC) (Kw) (Kw)

12 4.44 500 2 2.74 2.38 13.13868613

18 4.44 600 3 3.16 3.107 1.67721519
24 4.44 800 3 4.15 4.27 2.891566265
30 4.44 1000 3 5.16 5.44 5.426356589
36 4.44 1200 3 5.68 5.95 4.753521127
Error Promedio 5.57746906
Fuente. Autores

Tabla 29. Comparación de carga sensible del catalogo Vs Carga sensible del
programa con Tev=7,22 ºC
CATALOGO R1MX - YORK R22
T CFM NF TSH TSH
Modelo evaporización entrada total catalogo programa Error (%)
R22 (ºC) (filas) (Kw) (Kw)
12 7.22 500 2 2.4 2.17 9.583333333
18 7.22 600 3 2.75 2.85 3.636363636
24 7.22 800 3 3.62 3.53 2.486187845
30 7.22 1000 3 4.57 4.603 0.722100656
36 7.22 1200 3 5.07 4.97 1.972386588
Error Promedio 3.680074412
Fuente. Autores

129
 Tabla 30. Comparación de humedad relativa del catalogo Vs humedad
relativa del programa con Tev=4,44 ºC

CATALOGO R1MX - YORK R22

T CFM Humedad Humedad

NF total
Modelo evaporización entrada catalogo programa Error (%)
(filas)
R22 (ºC) (%) (%)

12 4.44 500 2 77 69.89 9.23376623

18 4.44 600 3 84 79.6 5.23809524
18 4.44 800 3 78 80.45 3.14102564
30 4.44 1000 3 78 82.036 5.17435897
36 4.44 1200 3 82 78.27 4.54878049
Error Promedio 5.46720532
Fuente. Autores

Tabla 31. Comparación de humedad relativa del catalogo Vs humedad

relativa del programa con Tev=7,22 ºC
CATALOGO R1MX - YORK R22

T CFM Humedad Humedad

entrada NF total
Modelo evaporización catalogo programa Error (%)
(filas)
R22 (ºC) (%) (%)

12 7.22 500 2 79 70.79 10.3924051

18 7.22 600 3 77 78.54 2
24 7.22 800 3 79 79.59 0.74683544
30 7.22 1000 3 76 79,73 4.9
36 7.22 1200 3 63 77.78 23.4603175
Error Promedio 8.29991159
Fuente. Autores

130
4.2 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA
HERRAMIENTA TENIENDO EN CUENTA VARIACIONES EN LOS
PARÁMETROS CONOCIDOS

Siendo ya corroborado los resultados de la herramienta computacional tanto con

fuente teórica como experimental, se prosiguió con analizarlos. También se
efectúo un análisis cuando los parámetros conocidos cambian, de tal modo que
se observe el efecto que tiene ellos sobre las condiciones de salida del aire y los
calores extraídos.

4.2.1 Análisis del Efecto de las condiciones de entrada de los fluidos de

trabajo

Para el modelo 4 YG del catálogo de la Unidad Manejadora de Aire (UMA) de

York, con geometría según figura N. y con las siguientes condiciones de entrada:

Flujo de aire: 1600,99 CFM

Presión atmosférica: 101.325 kPa
Temperatura de bulbo seco de entrada del aire: 26,7
Temperatura de bulbo húmedo de entrada del aire: 19,4
Flujo del agua: 6,16 GPM
Temperatura de entrada del agua: 6,7

131
 Figura 58. Geometría Modelo 4 YG del catalogo de York con NF=3

Fuente. Autores

Se realizó las variaciones a cada ítem de ellos, tal como se muestran en las
siguientes tablas:

Tabla. Efecto de la variación del flujo de aire sobre las condiciones de salida del
aire y de los calores extraídos por la unidad

Tabla. Efecto de la variación del bulbo seco del aire sobre las condiciones de
salida del aire y de los calores extraídos por la unidad

Tabla. Efecto de la variación del bulbo húmedo de aire sobre las condiciones de
salida del aire y de los calores extraídos por la unidad

132
Tabla. Efecto de la variación del flujo del agua sobre las condiciones de salida del
aire y de los calores extraídos por la unidad

133
 Tabla 32. Efecto de la variación del flujo de aire sobre las condiciones de salida del aire y de los calores
extraídos por la unidad
DATOS DE ENTRADA PARA LA UMA MODELO 4
DATOS DE SALIDA DEL PROGRAMA
YG NF=3
ITEM AIRE AGUA CALOR CALOR CALOR
Ta_out
SENSIBLE LATENTE TOTAL
CFM TBS (ºC) TBH (ºC) GPM T (ºC) (ºC)
(kW) (kW) (kW)
1 800 26,7 19,4 6,16 6,7 4,76139 1,43652 6,19791 17,408 79,4938
2 900 26,7 19,4 6,16 6,7 4,9918 1,46355 6,45473 18,0287 77,3801
3 1100 26,7 19,4 6,16 6,7 5,37027 1,48087 6,85114 19,0445 74,0689
4 1400 26,7 19,4 6,16 6,7 5,85316 1,51187 7,36503 20,1248 70,5616
5 1600,99 26,7 19,4 6,16 6,7 6,09827 1,4866 7,58488 20,6968 68,83
6 1800 26,7 19,4 6,16 6,7 6,3336 1,48335 7,81695 21,1463 67,4347
7 2000 26,7 19,4 6,16 6,7 6,52024 1,43896 7,9592 21,545 66,2789
8 2200 26,7 19,4 6,16 6,7 6,69834 1,40406 8,1024 21,8785 65,3011
9 2400 26,7 19,4 6,16 6,7 6,87337 1,38319 8,25656 22,1595 64,6408
10 2500 26,7 19,4 6,16 6,7 6,98372 1,41164 8,39537 22,271 64,0748
Fuente. Autores

134
Observando la tabla anterior se puede apreciar lo siguiente:

Para el calor sensible, el aumento es progresivo es decir es un aumento continuo

pero no proporcionado entre uno y otro valor, que está alrededor de menos de 1
kW. En la siguiente figura se observa los puntos graficados de la tabla de estas
dos variables.

Figura 59. Flujo de aire vs Calor sensible

Fuente. Autores

Para el calor latente, a medida que aumenta el flujo de aire, disminuye el calor
latente, lo cual es coherente con el fenómeno registrado ya que a mayor flujo de
aire más turbulencia externa lo cual hace que el coeficiente convectivo de calor
aumente junto con el coeficiente difusivo del potencial de masa.

Para el calor total, a medida que aumenta el flujo de aire, aumenta el calor total, lo
cual concluye que siendo el calor total la suma del latente y el sensible, y
respectivamente, el primero disminuyendo y el segundo aumentando, el
comportamiento del total se ve influenciado mucho más por el sensible ya que hay
más proporción de este en comparación con el latente, sin embargo mientras más
caiga el latente menor va a ser el aumento del total.

135
 Figura 60. Flujo de aire vs Calor total

Fuente. Autores

Para la temperatura de bulbo seco a la salida, el aumento es progresivo es decir

es un aumento continuo pero no proporcionado entre uno y otro valor, lo cual se
ve reflejado en el resultado de calor sensible extraído.

Para la humedad relativa a la salida, a medida que aumenta el flujo de aire,

disminuye, lo cual se ve reflejado en el resultado de calor latente extraído.

136
 Tabla 33. Efecto de la variación del la temperatura del bulbo seco de aire sobre las condiciones de salida del
aire y de los calores extraídos por la unidad
DATOS DE ENTRADA PARA LA UMA MODELO 4
DATOS DE SALIDA DEL PROGRAMA
YG NF=3
ITEM AIRE AGUA CALOR CALOR CALOR
Ta_out
SENSIBLE LATENTE TOTAL
CFM TBS (ºC) TBH (ºC) GPM T (ºC) (ºC)
(kW) (kW) (kW)
1 1600,99 20 19,4 6,16 6,7 3,28999 4,3753 7,66529 17,1749 97,8072
2 1600,99 21 19,4 6,16 6,7 3,72432 3,95756 7,68187 17,6542 92,826
3 1600,99 22 19,4 6,16 6,7 4,12996 3,50399 7,63396 18,2357 88,0875
4 1600,99 24 19,4 6,16 6,7 7,62104 2,64735 7,62104 19,2839 79,3226
5 1600,99 26,7 19,4 6,16 6,7 6,09827 1,4866 7,58488 20,6968 68,83
6 1600,99 28 19,4 6,16 6,7 6,65204 0,949893 7,60103 21,3604 64,2661
7 1600,99 29 19,4 6,16 6,7 7,05148 0,505427 7,5569 21,8914 60,9616
8 1600,99 30 19,4 6,16 6,7 10,8688 0,204196 11,073 22,4097 57,4753
9 1600,99 31 19,4 6,16 6,7 15,2202 0,0461591 15,2663 22,9544 53,691
10 1600,99 33 19,4 6,16 6,7 19,2083 0 19,2083 24,1408 45,7631
Fuente. Autores

137
Observando la tabla anterior se puede apreciar lo siguiente:

Para el calor sensible, cuando se incrementa el bulbo seco paulatinamente crece

el calor sensible, sin embargo su incremento no es proporcionado entre uno y otro
grado de bulbo seco, siendo el delta de calor sensible más pequeño cuando las
temperaturas de bulbo seco y bulbo húmedo de entrada son próximas entre ellas,
es decir, cuando la Tbs=20 ºC y la Tbh=19.4 ºC.

Para el calor latente, cuando se incrementa el bulbo seco paulatinamente decrece

el calor latente, ya que la humedad absoluta de entrada se reduce lo cual implica
reducción del potencial de masa, esto se reafirma con los datos de la tabla en
donde se ve que a medida que se aleja la Tbs de la Tbh decrece el calor latente
hasta aproximarse a cero.

Para el calor total, cuando se incrementa el bulbo seco paulatinamente crece el

calor total, ya que la humedad absoluta de entrada se reduce lo cual implica
reducción del potencial de masa, esto se reafirma con los datos de la tabla en
donde se ve que a medida que se aleja la Tbs de la Tbh decrece el calor latente
hasta aproximarse a cero.

Para la temperatura de bulbo seco a la salida, el aumento es progresivo es decir

es un aumento continuo pero no proporcionado entre uno y otro valor, lo cual se
ve reflejado en el resultado de calor sensible extraído.

Para la humedad relativa a la salida, a medida que aumenta el flujo de aire,

disminuye, lo cual se ve reflejado en el resultado de calor latente extraído.

138
Tabla 34. Efecto de la variación del bulbo húmedo de aire sobre las condiciones de salida del aire y de los
calores extraídos por la unidad
DATOS DE ENTRADA PARA LA UMA MODELO 4
DATOS DE SALIDA DEL PROGRAMA
YG NF=3
ITEM AIRE AGUA CALOR CALOR CALOR
Ta_out
SENSIBLE LATENTE TOTAL
CFM TBS (ºC) TBH (ºC) GPM T (ºC) (ºC)
(kW) (kW) (kW)
1 1600,99 26,7 16 6,16 6,7 13,8298 0 13,8298 20,2963 47,4652
2 1600,99 26,7 17 6,16 6,7 14,2802 0 14,2802 20,2946 55,2462
3 1600,99 26,7 18 6,16 6,7 9,18051 0,197547 9,37806 20,3632 62,4413
4 1600,99 26,7 18,5 6,16 6,7 6,22081 0,58453 6,80534 20,4842 64,8997
5 1600,99 26,7 19 6,16 6,7 6,15405 1,08358 7,23763 20,6008 67,0777
6 1600,99 26,7 19,4 6,16 6,7 6,09827 1,4866 7,58488 20,6968 68,83
7 1600,99 26,7 20 6,16 6,7 5,99793 2,08057 8,0785 20,8565 71,4736
8 1600,99 26,7 21 6,16 6,7 5,84442 3,11786 8,96229 21,1138 75,9172
9 1600,99 26,7 21,5 6,16 6,7 5,75935 3,63979 9,39914 21,2512 78,1536
10 1600,99 26,7 22 6,16 6,7 5,68551 4,18668 9,87219 21,3799 80,4053
Fuente. Autores

139
Observando la tabla anterior se puede apreciar lo siguiente:

Para el calor sensible, cuando se incrementa el bulbo húmedo paulatinamente

decrece el calor sensible; el decremento es pronunciado cada vez más que
aumenta en un grado la temperatura de bulbo húmedo a la entrada.

Para el calor latente, cuando se incrementa el bulbo húmedo paulatinamente el

calor latente aumenta, lo cual implica ampliación del potencial de masa, esto se
reafirma con los datos de la tabla en donde se ve que a medida que acerca la Tbh
de la Tbs, el calor latente aumenta vertiginosamente, pero debe notarse que este
tipo de proceso no es recomendable para lograr el confort de las personas.

Para el calor total, cuando se incrementa la temperatura de bulbo húmedo a la

entrada del serpentín paulatinamente decrece el calor total, ya que se observa la
fuerte dependencia del calor total al calor sensible, pues el potencial de
transferencia de calor debido a la variación de temperatura es de mayor peso que
el de humedades absolutas.

Para la temperatura de bulbo seco a la salida, la disminución es progresiva es

decir es un decremento continuo pero no proporcionado entre uno y otro valor, lo
cual se ve reflejado en la reducción del calor sensible extraído.

Para la humedad relativa a la salida, a medida que aumenta el flujo de aire,

aumenta, lo cual se ve reflejado en la elevación del calor latente extraído.

140
Tabla 35. Efecto de la variación del flujo del agua sobre las condiciones de salida del aire y de los calores
extraídos por la unidad
DATOS DE ENTRADA PARA LA UMA MODELO 4
DATOS DE SALIDA DEL PROGRAMA
YG NF=3
ITEM AIRE AGUA CALOR CALOR CALOR
Ta_out
SENSIBLE LATENTE TOTAL
CFM TBS (ºC) TBH (ºC) GPM T (ºC) (ºC)
(kW) (kW) (kW)
1 800 26,7 19,4 4,8 6,7 5,61743 0,844966 6,4624 21,1149 68,9591
2 900 26,7 19,4 5 6,7 5,65211 0,892386 6,54449 21,0848 68,9478
3 1100 26,7 19,4 5,7 6,7 5,66962 0,916522 6,58614 21,0697 68,9412
4 1400 26,7 19,4 6 6,7 5,66962 0,916522 6,58614 21,0697 68,9412
5 1600,99 26,7 19,4 6,16 6,7 6,09827 1,4866 7,58488 20,6968 68,83
6 1800 26,7 19,4 7 6,7 7,13554 2,75538 989.092 19,7857 68,8333
7 2000 26,7 19,4 7,5 6,7 7,28386 2,92695 10,2108 19,6546 68,8517
8 2200 26,7 19,4 8 6,7 7,41508 3,07745 10,4925 19,5385 68,8691
9 2400 26,7 19,4 8,5 6,7 7,52955 3,20719 10,7367 19,437 68,8875
10 2500 26,7 19,4 9 6,7 7,61095 3,29882 10,9098 19,3648 68,9015
Fuente. Autores

141
Observando la tabla anterior se puede apreciar lo siguiente:

Para el calor sensible, cuando se incrementa el flujo de agua dentro del serpentín
paulatinamente crece el calor sensible; el crecimiento es mucho menor al
registrado en la variación de Tbs, esto se debe a que al aumentar la rata de flujo
del agua, el calor total a extraer es mayor por tanto se observa un crecimiento en
el calor sensible pero sin embargo, es mucho menor al observado en el calor
latente.

Para el calor latente, cuando se incrementa el flujo de agua dentro del serpentín
ocurre un crecimiento mucho mayor que el registrado en el calor sensible. Y es
debido a que el potencial de temperaturas afecta más directamente al calor
sensible mientras que este cambio afecta mucho más al latente sin ser igual el
incremento que se produce por variación de la temperatura de bulbo húmedo del
aire a la entrada del serpentín.

Para el calor total, siendo el calor total función de la rata de flujo del agua, un
incremento de la rata de flujo del agua representa un aumento proporcional en el
calor total a extraer.

Para la temperatura de bulbo seco a la salida, cuando la rata de flujo del agua
aumenta, la turbulencia interna aumenta haciendo que el coeficiente interno
convectivo crezca, en decremento de la caída de temperatura del agua y por ende
de las condiciones del aire a la salida.

Para la humedad relativa a la salida, ocurre un decremento cuando la rata de flujo

del agua aumenta, debido la turbulencia interna ya que el coeficiente interno
convectivo crece, en decremento de la caída de temperatura del agua y por ende
de las condiciones del aire a la salida.

142
4.3 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA
HERRAMIENTA TENIENDO EN CUENTA LA EVOLUCIÓN DEL AIRE AL PASO
POR LAS FILAS DEL SERPENTÍN

También se efectúo un análisis de la evolución del aire cuando el fluido de trabajo

es agua para conocer el comportamiento de las condiciones de salida del aire y
de los calores extraídos.

4.3.1 Análisis del Efecto de las condiciones de entrada de los fluidos de

trabajo

Para ilustrar lo anterior, se tomó como referencia el modelo 2 YG del catálogo de

York con número de filas tres (NF=3) cuya geometría se muestra en la figura n. y
con las siguientes condiciones de entrada:

Flujo de aire: 800,496 CFM

Presión atmosférica: 101.325 kPa
Temperatura de bulbo seco de entrada del aire: 26,7
Temperatura de bulbo húmedo de entrada del aire: 19,4
Flujo del agua: 3,52 GPM
Temperatura de entrada del agua: 6,7

Se calculo para diferentes CFM en el Modelo 2 YG NF=3, la evolución del aire por
filas tanto la temperatura de bubo seco, humedad relativa y el cambio de los
calores al paso por el serpentín enfriador, tal como se muestran en la siguiente
tabla:

Tabla. Variación de las condiciones de salida del aire y de los calores extraídos al
paso por cada fila

143
Figura 61. Geometría del modelo 2 YG del catalogo de York con NF=3

Fuente. Autores

144
 Tabla 36. Variación de las condiciones de salida del aire y de los calores extraídos al paso por cada fila
UNIDAD MANEJADORA MODELO 2 YG NF=3
400,496 800,496
NF/CFM
Qt Ql Ta_out Qt Ql Ta_out
1 1014,45 160,49 22,89 62,41 1089,30 76,77 24,44 58,01
2 988,46 326,04 19,96 70,60 1155,44 280,48 22,50 63,71
3 947,66 402,21 17,58 76,22 1211,53 432,13 20,78 68,22

Condición salida
2950,57 888,74 17,58 76,22 3456,27 789,38 20,78 68,22
de la UMA

1500,496 2800,496
NF/CFM
Qt Ql Ta_out Qt Ql Ta_out
1 1202,04 6,17 25,27 55,51 2121.64 0 25,70 54,14
2 1268,89 184,11 23,99 59,47 2122,71 0 24,78 57,19
3 1389,75 381,23 22,80 62,84 1267,62 69,46 24,02 59,76

Condición salida
3860,68 571,51 22,80 62,84 5511,97 69,46 24,02 59,76
de la UMA

Fuente. Autores

145
De los resultados obtenidos de la grafica anterior se puede concluir lo siguiente:

Que a medida que aumenta el número de filas en una unidad manejadora de

aire: el calor latente, la temperatura de bulbo seco y la humedad relativa de la
salida van aumentando paulatinamente.

También que a mayor CFM la evolución de la temperatura y el calor latente se

va reduciendo tanto la total de la salida del serpentín como la parcial por filas.

Para un aumento grande de las CFM, ejemplo: 2800,496 CFM el calor total es
muy grande (5511,97 W) en comparación con el de 400,496 CFM (2950,57 w).
Sin embargo, el calor latente es mayor para el segundo que para el primero
(888,74 W contra 69,46 W), si esta unidad con 2800,496 CFM se colocara en
Bucaramanga no funcionaria bien ya que el aire tiene humedad relativa mayor
a 70 lo cual hace que se requiere disponibilidad de extraer humedad para
lograr el confort humano en un recinto cerrado.

4.4 COMPARACIÓN Y ANÁLISIS DE LOS RESULTADOS DE LA

HERRAMIENTA TENIENDO EN CUENTA LA VARIACION DEL NUMERO DE
FILAS

En las siguientes tablas se realizo cambio a la geometría del modelo YG 02 y la

unidad de expansión directa R1MX 018. La geometría fue modificada al cambiar el
número de filas de la unidad. El cambio de número de filas se realizo manteniendo
el flujo del refrigerante (Agua o R-22) constante.

Condiciones de entrada para la unidad YG York 02

Temperatura de bulbo seco a la entrada Tbs=26.7°C
Temperatura de bulbo húmedo a la entrada Tbh=19.4 °C

146
 Presión atmosférica=101.325kpa

Tabla 37. Variación de la geometría para la unidad YG York 02

CATALOGOS YG YORK O2
Flujo
TSH T
de NF TLH TH Humedad
Programa programa
agua (filas) Programa(Kw) Programa(Kw) (%)
(Kw) (ºC)
(Kg/s)
3.52 3 3.0495 0.533 3.5825 20.606 70.28
3.52 4 5.267 0.595 5.85 19.5 75.1
3.52 6 6.08 0.922 7.002 17.52 82.68
3.52 8 7.161 1.6508 8.8118 15.47 88.37
Fuente. Autores

De la tabla anterior se observa:

Manteniendo el flujo constante y cambiando el número de filas al modelo 2 YG

del catalogo de York, ocurre un aumento en el calor sensible mucho más
grande que el que se registra con el calor latente, lo cual conlleva a ampliar el
calor total.

Condiciones de entrada para la unidad R1MX 018

Temperatura de bulbo seco a la entrada Tbs=23.5°C
Temperatura de bulbo húmedo a la entrada Tbh=16.67 °C
Se asumió una temperatura de subenfriamento=2°C
Y se asumió una temperatura de alta=35°C
Presión atmosférica=101.325 kPa

147
 Tabla 38. Variación de la geometría para R1MX 018

CATALOGO FANCOIL R1MX YORK 018

Flujo de TSH T
NF TLH TH Humedad
Refrigerante Programa programa
(filas) Programa(Kw) Programa(Kw) (%)
R22 (Kg/s) (Kw) (ºC)
0.02 3 3.022 1.527 4.549 12 76.64
0.02 4 3.56 1.884 5.444 10.13 83.7
0.02 6 4.278 2.415 6.693 7.71 92.43
0.02 8 4.7 2.76 7.46 6.66 97
Fuente. Autores

De la tabla anterior se observa:

Lo mismo que ocurre con el modelo del agua sólo que el incremento en el calor
latente es un tanto mayor para refrigerante R-22.

148
 CONCLUSIONES

La herramienta computacional permite comprobar y / o determinar la capacidad

de retirar calor y las condiciones de salida del aire (Temperatura, humedad) en
las unidades enfriadoras ya sea de expansión directa o enfriadas por agua
cuando operan por datos del fabricante o personalizada.
La comparación entre la herramienta computacional desarrollada en Matlab y
la información suministrada por los fabricante (Catálogos) permitieron
establecer en:
a) Unidades de expansión directa
El error que se presenta en los resultados obtenidos respecto de la
temperatura, carga sensible, carga latente de los diferentes modelos de
fancoils R1MX de York son los siguientes:
CATALOGO R1MX - YORK R22
T Error Error
error carga error carga
Modelo evaporización humedad Temperatura
latente (%) sensible (%)
R22 (ºC) relativa (%) (%)
12 4.44 4.34210526 13.13868613 9.23376623 4.39164867
12 7.22 6.42857143 9.583333333 10.3924051 5.53333333
18 4.44 2.94871795 1.67721519 5.23809524 3.7437005
18 7.22 9.25925926 3.636363636 2 4.49871465
24 4.44 1.92307692 2.891566265 3.14102564 9.35595568
24 7.22 3.4375 2.486187845 0.74683544 5.84832905
30 4.44 2.27272727 5.426356589 5.17435897 11.3573407
30 7.22 6.52173913 0.722100656 4.9 7.00514139
36 4.44 6.92307692 4.753521127 5.17435897 8.4
36 7.22 12.8205128 1.972386588 10.34567 6.20732464
Fuente: Autores
b) Unidades enfriada por agua
Los resultados obtenidos respecto de la temperatura, carga sensible, carga
latente evaluado por la herramienta para la unidad manejadora YG York
difieren con la del fabricante de la siguiente manera:
CATALOGO YG YORK
Error Error
Error carga Error carga
Modelo humedad Temperatura
latente (%) sensible (%)
relativa (%) (%)
7 10 10.1 8.11 9.35
3 6.33 6.44 8.98 6.1
9 5.54 5.567 6.778 3.72

149
 CATALOGO YG YORK
Error Error
Error carga Error carga
Modelo humedad Temperatura
latente (%) sensible (%)
relativa (%) (%)
30 6.99 6.23 5.66 7.00514139
12 8.33 4.04 5 8.4
Fuente: autores

La herramienta permite al usuario modificar los datos en cualquier momento

del proceso, lo que posibilita el análisis de múltiples situaciones, con variar uno
o más parámetros en un mismo problema. Esta herramienta resulta ser útil
debido a que se adapta a las necesidades del usuario para determinar un
análisis más exacto sobre el efecto de cada variable en el desempeño de las
unidades manejadoras.
La comparación de la herramienta computacional y a las pruebas tomadas de
la Unidad manejadora de aire para agua helada, marca Páramo ubicada en el
edificio de ciencias humanas sala de cómputo de la Universidad Industrial de
Santander se pudo concluir que:
a. La unidad evaluada según el criterio de la herramienta se encuentra con
una humedad relativa del 100%
b. Los errores promedios que difieren entre carga latente, carga total y
temperatura y humedad relativa son:
 Carga latente. (10%)
 Carga total (6%)
 Temperatura de salida (20%)
 Humedad relativa (29%)
Los porcentajes de error que se presentaron en los resultados obtenidos en la
herramienta computacional en relación a unidades reales están dentro de los
rangos aceptables. La HTRI (Heat Transfer Research Institute) considera
después de múltiples y cuidadosas investigaciones que los errores
presentados en experimentos de laboratorios están en el rango de 15 %. Ref.
Solution of Shell Side Flow Pressure Drop and Heat Transfer By Stream.
Analysis Method.No. 62 Vol. 65, J.W Pal and Jerry Tabor

150
 RECOMENDACIONES

Se recomienda incluir más catálogos de unidades manejadoras dentro de las

bases del programa para crear versiones futuras más completas del programa
obteniendo posiblemente un porcentaje de error más favorable.

Se recomienda probar otras correlaciones experimentales para calcular el

coeficiente de transferencia de ebullición-forzada en tubos horizontales para
refrigerante R-22.

Se recomienda probar otros modelos de transferencia de calor para la

superficie de la unidad enfriadora seca- húmeda al paso del aire.

Se recomienda leer el manual de usuario del funcionamiento del programa

antes de ingresar a evaluar unidades manejadoras.

151
 REFERENCIAS BIBLIOGRAFICAS

[1] CENGEL Yunus. Termodinámica. Cuarta Edición. Editorial McGraw Gill.

Capítulo 12 Mezclas de gases. Capítulo 13 Mezcla de gas vapor y
Acondicionamiento de aire.

[2] MARADEY Juan F. Termodinámica Aplicada. Ediciones UIS. Capítulo 2

Equilibrio de Fases

[3] CARRIER AIR CONDITIONING. Manual de aire acondicionado. Boixareu

Editores. Capítulo 8 Empleo del Diagrama Psicrometrico. 1980. Págs. 109-
125.

[4] QUADRI Néstor. Sistemas de Aire Acondicionado. Primera Edición. Capítulo

2 Confort Térmico. 2001

[5] CENGEL Yunus. Thermodynamics.Fifth edition. Chapter 14 Mixture of gas-

vapor and air conditioning.

[6] CENGEL Yunus. Transferencia de calor. Segunda edición. Capítulo 13

Intercambiadores de calor y Capítulo 3 Conducción de calor en estado
estable 2008

[7] ASHRAE. Handbook de Refrigeración. Capítulo 15 Sistemas de aire:

componentes-ventiladores, serpentines, filtros y deshumidificadores. 2002

[8] INCROPERA De Witt. Transferencia de calor. Capítulo 11 Intercambiadores

de calor. Cuarta Edición.

152
[9] GELVEZ A. Omar. Transferencia de calor. Capítulo 9 Intercambiadores de
calor. Págs. 250-255.

[10] MILLS A. F. Transferencia de calor. Capítulo 8 Intercambiadores de calor.

Capítulo 4 Fundamentos y correlaciones de la convección. Primera edición.

[11] FERNANDEZ D. Pedro. Transferencia de Calor, capítulo X Superficies

ampliadas de sección transversal variable.

[12] McQUISTON Faye C. Heating, Ventilating, and Air Conditioning Analysis and
Design. Chapter 14 Extended surface heat exchangers. Fourth Edition.

[13] KAKAC.Boilers, evaporators, and Condenser. Chapter 12 Evaporators and

condensers for refrigeration and air-conditioning systems. Editorial John
Wiley & Sons, Inc 1991.

[14] PITA Edward. Principios y Sistemas de Refrigeración. Termodinámica del

ciclo de refrigeración por compresión de vapor. Págs.86-87. Editorial Limusa.
2006

[15] ASHRAE. Handbook de Aire Acondicionado y refrigeración. Capítulo 9.

Refrigerantes, ciclos de refrigeración y sistemas de refrigeración. Pás. 26-27.
Año 2001

153
PAPERS:

 Pirompugd W, Wang.C.C., Wongwises S.A fully wet and fully dry tiny circular fin
method for heat and mass transfer characteristics for plain fin-and-tube heat
exchangers under dehumidifyng conditions, by W.. J.Heat Transfer 129
(9)(2007).

 Kandlikar S. G., A general correlation for saturated two-phase flow boiling heat
transfer inside horizontal and vertical tubes, by W.. J Heat Tranfer 112 (1990)

154
ANEXOS

155
ANEXO A. MANUAL TECNICO DE COMFORT PACKAGE FAN & COIL UNIT-
YG CAPACIDADES DE 1360 A 10200 m3/h

156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
ANEXO B. CATALOGO DE YORK UNIDADES FANCOILS DE EXPANSIÓN
DIRECTA CAPACIDADES DE 12000 A 36000 BTU/H

191
192
193
194
ANEXO C. CORRELACIÓN GENERAL PARA FLUJO DE DOS FASES
SATURADO PARA TRANSFERENCIA DE CALOR DENTRO DE TUBOS
VERTICALES U HORIZONTALES

195
196
197
198
199
200
201
202
203
204
205