Tesis

Análisis del desempeño energético y desarrollo del coeficiente de
transferencia de calor y masa en un sistema de acondicionamiento de aire

residencial para clima tropical
Tesis para optar al título de doctor en ingeniería mecánica
Hermes Javier Ramírez León

Candidato a Doctor en Ingeniería Mecánica
Director
Antonio José Bula Silvera
PhD. en Ingeniería Mecánica
Programa de Doctorado en ingeniería mecánica

Programa de Ingeniería mecánica
División de Ingenierías
Universidad del norte
Barranquilla, 30 de agosto de 2019
Tesis doctoral
Análisis del desempeño energético y desarrollo del coeficiente de
transferencia de calor y masa en un sistema de acondicionamiento de aire
residencial para clima tropical
Realizada por: Hermes Javier Ramírez León

Dirigida por: Antonio José Bula Silvera PhD
Tribunal calificador:
Presidente: Mauricio Yilmer Carmona García, PhD

Jurado 1: Cesar Isaza Roldán, PhD
Jurado 2: Marco Sanjuán Mejía, PhD
Jurado 3: Ricardo Vásquez Padilla, PhD
Programa de Doctorado en ingeniería mecánica

Programa de Ingeniería mecánica
División de Ingenierías
Universidad del norte
Barranquilla, 30 de agosto de 2019
Agradecimientos
A Dios, gracias por tu bendición y fortaleza

Mi esposa Diana, mis padres y hermanos, gracias por su amor, paciencia y
compresión
Mi tutor Antonio Bula, mis amigos, compañeros y cada uno que intervino para
que este trabajo pudiera realizarse a todos muchas gracias. Por último, a la
Universidad del Norte por el contrato de beca No UN-OJ-2013-22062 y también
al Departamento Administrativo de Ciencia, Tecnología e Innovación
COLCIENCIAS por la Convocatoria No 727 de 2015; gracias por el apoyo
financiero y académico.
Contenido
Lista de figuras....................................................................................................................................................... 6
Lista de tablas ........................................................................................................................................................ 8
Nomenclatura ......................................................................................................................................................... 9
Resumen .............................................................................................................................................................. 11
Introducción ......................................................................................................................................................... 13
Justificación ......................................................................................................................................................... 13
Planteamiento del problema .................................................................................................................................. 24
Objetivo General .................................................................................................................................................. 24
Objetivos específicos ............................................................................................................................................ 24
1 Análisis energético del aire típico de un clima tropical sobre un aire acondicionado de uso residencial .......... 25
1.1 Características del aire en el clima tropical ........................................................................................... 25
1.2 Diseño de Experimentos ...................................................................................................................... 25
1.2.1 Factores ...................................................................................................................................... 26
1.2.2 Variables de respuesta ................................................................................................................. 26
1.2.3 Niveles factor .............................................................................................................................. 27
1.2.4 Factor bloque .............................................................................................................................. 27
1.3 Descripción del sistema........................................................................................................................ 27
1.4 Puesta en marcha y Resultados ............................................................................................................. 31
1.4.1 Calculo tipo................................................................................................................................. 33
1.5 Metodología de análisis de resultados ................................................................................................... 36
1.6 Análisis de resultados para el calor absorbido en el evaporador ............................................................. 37
1.7 Análisis de resultados para el trabajo eléctrico del compresor ............................................................... 38
1.8 Análisis de resultados para el coeficiente de desempeño del equipo ...................................................... 40
1.9 Conclusiones del capítulo..................................................................................................................... 44
2 Análisis del proceso de transferencia de calor entre el aire de clima tropical y el refrigerante R410a en el
interior de un evaporador ...................................................................................................................................... 45
2.1 Descripción matemática ....................................................................................................................... 46
2.2 Cálculos de coeficiente de transferencia de calor y masa en un evaporador de tubos aleteados ............... 49
2.3 Análisis de resultados para el coeficiente de transferencia de calor........................................................ 52
3 Análisis estadísticos a partir del resultado en el sistema experimental ............................................................ 55
3.1 Modelo de regresión ............................................................................................................................ 55
3.2 Trazabilidad del error........................................................................................................................... 61
3.3 Propuesta de estudio: Mapa de funcionamiento del equipo de aire acondicionado ................................. 64
4 Modelo dinámico del evaporador.................................................................................................................. 67
4.1 Formulación matemática ...................................................................................................................... 67
4.2 Modelación del aire húmedo ................................................................................................................ 71
4.3 Solución del modelo dinámico ............................................................................................................. 72
4.3.1 Variables dinámicas .................................................................................................................... 73
4.3.2 Condiciones iniciales ................................................................................................................... 73
4.4 Resultados ........................................................................................................................................... 74
Conclusiones ........................................................................................................................................................ 78
Trabajos futuros ................................................................................................................................................... 80
Anexos ................................................................................................................................................................. 82
I. Selección del diseño de experimento ............................................................................................................ 82
a. Factores ................................................................................................................................................... 82
b. Variables de respuesta.............................................................................................................................. 82
c. Niveles factor .......................................................................................................................................... 82
d. Otras consideraciones del diseño de experimento...................................................................................... 83
e. Factor Bloque .......................................................................................................................................... 83
f. Descripción del sistema............................................................................................................................ 83
II. Resultados y validación del diseño de experimentos ..................................................................................... 84
a. Análisis de resultados para el calor absorbido en el evaporador ................................................................. 84
b. Análisis de resultados para el trabajo eléctrico del compresor ................................................................... 85
c. Análisis de resultados para el coeficiente de desempeño del equipo .......................................................... 87
III. Prueba de idoneidad de los modelos de regresión...................................................................................... 89
a. Modelo 1 ................................................................................................................................................. 89
b. Modelo 2 ................................................................................................................................................. 90
c. Modelo 3 ................................................................................................................................................. 90
d. Modelo 4 ................................................................................................................................................. 90
e. Modelo 5 ................................................................................................................................................. 91
f. Modelo 6 ................................................................................................................................................. 91
g. Modelo 7 ................................................................................................................................................. 91
h. Modelo 8 ................................................................................................................................................. 92
i. Modelo 9 ................................................................................................................................................. 92
Lista de figuras
Figura 0-1 Demanda de unidades de aire acondicionado residencial y comercial a nivel mundial desde
2011-2016. Fuente [3] ........................................................................................................................... 13
Figura 0-2 Distribución de la población mundial por latitud y longitud. Fuente [4]. ............................... 14
Figura 0-3 índice trilema del consejo mundial de la energía. fuente [5] ................................................. 15
Figura 0-4 Intercambio de calor entre el aire y el refrigerante dentro del evaporador de un sistema DX de
alta eficiencia (Izq. Condición de alta temperatura. Der. Condición de baja temperatura). Fuente:
elaboración propia. ................................................................................................................................ 17
Figura 0-5 Esquema para un modelo dinámico de un intercambiador de calor. Fuente: [8] .................... 21
Figura 0-6 Esquema del modelo de límites movibles. Fuente: elaboración propia .................................. 22
Figura 1-1 Carta psicrométrica con los niveles de los factores (Zona naranja) y la zona de confort para
verano (zona azul). Fuente: elaboración propia ...................................................................................... 28
Figura 1-2 Izq. Sistema Experimental Der. Etiqueta energética. Fuente: elaboración propia .................. 28
Figura 1-3 Esquema del sistema experimental y ubicación de los sensores Fuente: Elaboración propia.. 29
Figura 1-4 Esquema del generador de vapor, resistencia eléctrica y cámara de mezcla para suministrar
condiciones de entrada del aire al evaporador. Fuente: elaboración propia.............................................. 30
Figura 1-5 Hoja de Cálculo – Izq. Adquisición de datos. Der. Propiedades termodinámicas y cálculos
energéticos. ........................................................................................................................................... 32
Figura 1-6 Comparación: calor cedido por el aire vs calor absorbido por el refrigerante. Fuente:
elaboración propia ................................................................................................................................. 33
Figura 1-7 Potencia del compresor termodinámico vs potencia eléctrica del compresor. Fuente:
elaboración propia ................................................................................................................................. 33
Figura 1-8 Comportamiento flujo refrigerante con los lazos de control del inverter y la temperatura del
aire a la salida en condiciones extremas de operación. Fuente: elaboración propia .................................. 38
Figura 1-9 Resultado ANOVA Diagrama de Pareto para el ANOVA. Fuente: elaboración propia ......... 39
Figura 1-10 Gráfico de interacciones para el consumo del compresor. Fuente: elaboración propia ......... 39
Figura 1-11 Superficie de respuestas para el consumo eléctrico del compresor en kW. Ejes: Temperatura
codificada del aire a la entrada del evaporador, Humedad absoluta codificada del aire a la entrada del
evaporador. Izq. Nivel bajo de velocidad, el Rango de velocidad para este nivel y máximo Reynolds
alcanzado. Der. Nivel alto de velocidad, rango de la velocidad para este nivel y máximo Reynolds
alcanzado. Fuente: elaboración propia ................................................................................................... 40
Figura 1-12 Diagrama de Pareto para el ANOVA del COP. Fuente: elaboración propia ........................ 41
Figura 1-13 Gráfico de efectos principales (izquierda) e interacciones (derecha) para el Coeficiente de
Desempeño. Fuente: elaboración propia ................................................................................................. 42
Figura 1-14 Superficie de respuestas para el coeficiente de desempeño del equipo. Ejes: Temperatura
codificada del aire a la entrada del evaporador, Humedad absoluta codificada del aire a la entrada del
evaporador Izq. Nivel bajo de velocidad, rango de velocidad para este nivel y máximo Reynolds
alcanzado Der. Nivel alto de velocidad, rango de velocidad para este nivel y máximo Reynolds alcanzado.
Fuente: elaboración propia ..................................................................................................................... 43
Figura 2-1 Análisis de condensación en presencia de un gas no condensable en la evaporador. Fuente:
[19] ....................................................................................................................................................... 46
Figura 2-2 Comparación del coeficiente de transferencia de calor obtenido por el modelo desarrollado y
el coeficiente presentado en la tabla 13. Fuente: elaboración propia ....................................................... 51
Figura 2-3 Comparación de resultados del coeficiente de transferencia de calor en la carta psicrométrica
con base a las condiciones de entrada del aire. Izquierda: método tradicional asumiendo Le=1. Derecha:
método presentado en este documento. Fuente: elaboración propia. ....................................................... 52
Figura 2-4 Izquierda: Coeficiente global de transferencia de calor y masa. Derecha: Coeficiente de
transferencia de calor. Ambas gráficas en función de las condiciones de entrada del aire y del consumo
eléctrico del compresor. Fuente: elaboración propia ............................................................................... 53
Figura 3-1 Gráfico de residuos estudentizados vs predicho para el factor j H para verificar independencia.
Fuente: elaboración propia ..................................................................................................................... 60
Figura 3-2 Gráfico de residuos estudentizados vs predicho para el factor jM para verificación de
independencia. Fuente: elaboración propia ............................................................................................. 60
Figura 3-3 Comparación del modelo de regresión vs datos experimentales. Fuente: elaboración propia. 61
Figura 3-4 Mapa de la capacidad de enfriamiento de la unidad de aire acondicionado tipo mini split en
una tabla psicrométrica a diferentes velocidades. Fuente: elaboración propia ......................................... 65
Figura 4-1 Esquema del modelo dinámico del evaporador y aire. Fuente: elaboración propia ................ 72
Figura 4-2 Modelo dinámico del evaporador en SIMULINK. Fuente: elaboración propia ...................... 72
Figura 4-3 Superior. T aire entrada y T aire salida vs tiempo de simulación. (Inferior) T aire salida vs w
aire salida vs tiempo de simulación Fuente: elaboración propia .............................................................. 75
Figura 4-4 Variables espacio estado del evaporador y condiciones psicrométrica del aire a la entrada del
evaporador vs tiempo de computación. Fuente: elaboración propia. ....................................................... 76
Figura 4-5 Superior. T aire entrada y T aire salida vs tiempo de simulación. (Inferior) T aire salida vs w
aire salida vs tiempo de simulación bajo perturbación Fuente: elaboración propia .................................. 77
Figura 4-6 Variables espacio estado del evaporador y condiciones psicrométrica del aire a la entrada del
evaporador vs tiempo de computación bajo perturbación. Fuente: elaboración propia. ............................ 77
Figura II-1 Grafica de probabilidad normal de los efectos estandarizados para el calor absorbido por el
evaporador. Fuente: elaboración propia.................................................................................................. 84
Figura II-2 Grafica de probabilidad normal de los efectos estandarizados para el consumo eléctrico del
compresor. Fuente: elaboración propia................................................................................................... 85
Figura II-3 Gráfica de probabilidad normal para validar normalidad (superior-izquierda), Gráfica de
residuos vs orden de ejecución para validación de independencia (superior-derecha) y Gráfica de residuos
vs predicho para validación de homocedasticidad (inferior). Fuente: elaboración propia ......................... 86
Figura II-4 Grafica de probabilidad normal de los efectos estandarizados para el coeficiente de
desempeño. Fuente: elaboración propia.................................................................................................. 87
Figura II-5 Gráfica de probabilidad normal para validar normalidad (superior-izquierda), Gráfica de
residuos vs predicho para validación de homocedasticidad (superior-derecha) y Gráfica de residuos vs
orden de ejecución para validación de independencia (inferior). Fuente: elaboración propia ................... 88
Lista de tablas
Tabla 0-1 Comparación entre categorías de modelos dinámicos para intercambiadores de calor. ............ 21
Tabla 1-1 Factores y niveles del diseño de experimento ......................................................................... 27
Tabla 1-2 Nivel de humedad relativa según niveles de los factores ......................................................... 27
Tabla 1-3 Descripción de los sensores instalados en el sistemas experimental ........................................ 30
Tabla 1-4 Descripción de la ubicación de los sensores utilizados en el sistema experimental .................. 31
Tabla 1-5 Orden de ejecución de la experimentación.............................................................................. 32
Tabla 1-6 Orden de experimentación y resultados obtenidos de las variables de respuestas (factores
codificados) ........................................................................................................................................... 36
Tabla 1-7 Resultado ANOVA para el calor cedido del aire a través del evaporador en función de los
efectos principales y sus interacciones. AB: Temperatura-Humedad, AC: Temperatura-Nivel velocidad,
BC: Humedad-Nivel velocidad. ............................................................................................................. 37
Tabla 1-8 Resultado ANOVA para el consumo eléctrico del compresor en función de las interacciones.
AB: Temperatura-Humedad, AC: Temperatura-Nivel velocidad, BC: Humedad-Nivel velocidad. .......... 38
Tabla 1-9 Resultado ANOVA para el coeficiente de desempeño en función del nivel de velocidad: C, y
las interacciones. AB: Temperatura-Humedad, AC: Temperatura-Nivel velocidad, BC: Humedad-Nivel
velocidad. .............................................................................................................................................. 41
Tabla 2-1 Resultado del coeficiente global de transferencia de calor y masa, coeficiente de transferencia
de calor y coeficiente de transferencia de masa para cada corrida experimental realizada. ...................... 50
Tabla 2-2 Resultado, Difusividad másica entre vapor de agua-aire y el número de Lewis por cada corrida
experimental .......................................................................................................................................... 51
Tabla 3-1 Modelos de regresión recomendados por STATGRAPHICS (c=coeficiente; m=intercepto) .... 56
Tabla 3-2 Resultados del cálculo de factor JH y JM en las corridas experimentales en orden ascendente
según Re. .............................................................................................................................................. 57
Tabla 3-3 Análisis de idoneidad de los modelos de regresión para jHT .................................................... 57
Tabla 3-4 Análisis de idoneidad de los modelos de regresión para jMT .................................................... 58
Tabla 3-5 Resultado del modelo de regresión para los factores j HT y jMT ................................................. 59
Tabla 3-6 ANOVA significancia del modelo para jHT ............................................................................. 59
Tabla 3-7 ANOVA significancia del modelo para jMT ............................................................................ 59
Tabla 3-8 Resultado del análisis de los parámetros del modelo jHT ......................................................... 59
Tabla 3-9 Resultado del análisis de los parámetros del modelo jMT ......................................................... 59
Tabla 3-10 Propagación del error en las variables medidas ..................................................................... 63
Tabla 3-11 Propagación del error en las variables calculadas y de respuesta ........................................... 63
Tabla 3-12 Porcentaje de la propagación del error en las variables ......................................................... 64
Tabla 3-13 Rango de variación de la temperatura, humedad y velocidad del aire a la entrada del
evaporador ............................................................................................................................................ 64
Tabla 4-1 Tabla de los elementos de la matriz Z(x,u) para el evaporador ................................................ 70
Tabla 4-2 Condiciones iniciales de las variables de entrada, dinámicas y de salida del modelo dinámico 74
Tabla 4-3 Condiciones de estado estable de las variables dinámicas de entrada ...................................... 75
Nomenclatura
Símbolo Descripción
A Área (m2)
C Coeficiente
COP Coeficiente de Desempeño
CP Calor específico a presión constante (J/kgK)
DX Expansión directa
F Factor multiplicativo bifásico
FCV Volumen de control finito
G Flujo másico por unidad de aire (kg/sm2)
GTC Coeficiente de global de transferencia de calor y masa (kg/s)
h Entalpía (kJ/kg)
HR Humedad relativa (%)
HTC Coeficiente de transferencia de calor (W/m2 °C)
I Corriente (A)
j Factor de Colburn
L Longitud (m)
Le Número de Lewis
LM Límites movibles
LMED (ΔhLM) Diferencia de entalpía media logarítmica (kJ/kg)
LMTD Diferencia de temperatura media logarítmica (°C)
ṁ Flujo másico (kg/s)
NTU Número de unidades de transferencia
P Presión (kPa)
p Perímetro (m)
Pr Número de Prandtl
Q Calor (kW)
q Calor por unidad de área (kW/m2)
Re Número de Reynolds
S Factor de supresión de ebullición nucleada
Sc Número de Schmidt
T Temperatura (°C)
t Tiempo
U Velocidad frontal (m/s)
V Voltaje (V), Volumen (m3)
W Trabajo (kW)
w Humedad absoluta (kgagua/kgaire) / (gagua/kgaire)
X Calidad del vapor saturado del refrigerante
z Dirección de flujo del refrigerante
Letra griega
α Coeficiente Convectivo de transferencia (W/m2 °C)/ (kg/s)
γ Fracción de vacío
η0 Eficiencia total del intercambiador (%)
λ Longitud normalizada
μ Relación entre densidades
ρ Densidad (kg/m3)
σ Incertidumbre (según variable a evaluar) / (%)
Subíndices
1 Zona bifásica
2 Zona vapor sobrecalentado
a Aire promedio
A Zona externa del evaporador/contacto con el aire
a,i Entrada del aire
a,s_x Salida del aire en para la zona x (1 o 2)
agua Agua
Aire Aire
Comp Compresor
CS Sección transversal
DC Diámetro del collar de aleta = Diámetro externo tubo*2(espesor aleta)
Evap Evaporador
f Fase líquida del refrigerante
fg Vaporización del agua
FZ Forster y Zuber
G Transferencia global de calor y masa
g Fase vapor del refrigerante
HT Transferencia de calor
i Zona interna del evaporador/contacto con el refrigerante
i,e Zona de entrada el refrigerante
i,o Zona de salida del refrigerante
in Entrada
Latente Calor latente
mid Zona intermedia del lado del refrigerante
MT Transferencia de masa
O Condición de superficie libre/máxima
Out Salida
pared Pared de la(s) tubería(s) del evaporador
ref Refrigerante
Sensible Calor sensible
Total Total
Resumen
Esta investigación se realizó con el propósito de aportar información sobre el funcionamiento de
los sistemas de acondicionamiento de aire tipo residencial en la región Caribe colombiana y en
general, en las zonas tropicales húmedas, analizando el comportamiento energético, los efectos de
transferencia de calor y transferencia de masa del equipo bajo las características del clima tropical
colombiano, caracterizado por su alta humedad y temperatura durante todo el año.
El estudio se proyectó en tres etapas: Una fase experimental, una fase evaluativa y por último una
fase de modelación dinámica. Para la primera parte, se construyó un sistema experimental que
permite analizar energéticamente, un aire acondicionado doméstico de eficiencia media con
inverter y condiciones normales de operación, este sistema consta de 21 sensores instalados a lo
largo del sistema para identificar variables termodinámicas del refrigerante, variables
psicrométricas del aire y variables energéticas del compresor. Con este sistema experimental, se
estableció un diseño experimental 23 tomando como factores controlables el nivel de velocidad del
ventilador en el evaporador, la temperatura y humedad del aire a la entrada del evaporador. Las
variables de respuesta fueron el calor cedido del aire hacia el evaporador, el consumo eléctrico del
compresor y el COP. Se observó, que el calor absorbido no es afectado por ninguno de los factores
principales (Temperatura, humedad y nivel de velocidad del aire), sin embargo, el consumo
eléctrico del compresor y el COP del equipo (por el lado externo) si es afectado por las
interacciones de temperatura-humedad del aire, temperatura-nivel de velocidad del ventilador,
humedad-nivel de velocidad del ventilador. De igual forma, se establece que si el aire ingresa al
evaporador a 28°C ,14 g/kg de humedad (condiciones típicas del aire tropical) manteniendo una
alta velocidad del ventilador, su eficiencia energética puede ser 10% menor que si el aire que
ingresa al evaporador a 23°C y 14 g/kg, y un nivel bajo de velocidad del ventilador.
Durante la fase evaluativa, se analiza los efectos de la humedad en el fenómeno de transferencia
de calor. Esta fase fue propuesta para desarrollar un coeficiente de transferencia de calor, en
evaporadores de tubos aleteados en climas tropicales, para esto, se estableció una metodología de
cálculo que incluye el número de Lewis en el proceso para considerar los efectos de condensación
en presencia de gas no condensable en este estudio. Los resultados se compararon con un
coeficiente de transferencia de calor sin considerar la condensación; se aprecia cuando existe una
variación de la humedad, el coeficiente de transferencia de calor varía entre 4 a 16 W/-(m2 °C) en
comparación con el coeficiente sin considerar la condensación. Además, se obtuvo un coeficiente
global de transferencia de calor y masa, donde se observó una relación entre este coeficiente con
la velocidad del aire, y el consumo energético: Cuando el coeficiente global de transferencia de
calor y masa (αG)se encuentra por debajo de 0,1 kg/(m2s) o por encima de 0,15 kg/(m2s) se
observan los mayores consumos energéticos del sistema, mientras que, a niveles intermedios, el
consumo disminuye en un 30%, básicamente porque el fenómeno de transferencia de calor mejora
gracias al condensado generado por el aire húmedo. Finalmente, para terminar esta fase, se
obtuvieron dos modelos de regresión para establecer un factor de Colburn para la transferencia de
calor y un factor de Colburn para la transferencia de masa, y así poder obtener un coeficiente de
transferencia de calor y un coeficiente de transferencia de masa, aplicado a un evaporador de tubos
aleteados para el proceso de enfriamiento y deshumidificación.
En la última fase se establece un modelo dinámico del evaporador del sistema experimental a
través del método de límites movibles donde se incluyen los efectos de la condensación del aire
en la metodología, este modelo se realizó a través de SIMULINK, MATLAB y con apoyos de la
librería REFPROP, con el propósito de identificar el fenómeno interno y externo del evaporador y
establecer una tendencia o posible estrategia de control. Este modelo permite establecer hasta
donde avanza el fenómeno de cambio de fase, teniendo en cuenta que el flujo másico de
refrigerante no es suministrado por un plenum, es decir, existe un punto donde el refrigerante se
evapora completamente y comienza a sobrecalentarse.
Con todo esto, se pretende desarrollar un campo de estudio que ha avanzado en otras regiones, que
requiere mayor visibilidad hacía la comunidad científica.
En el primer capítulo, se presenta la construcción y desarrollo de un sistema experimental que
simula la operación típica de un equipo mini split con inverter de consumo moderado de energía.
Se presenta el diseño de experimento, los factores principales, las condiciones experimentales, las
variables de respuestas y consideraciones utilizadas para tal evaluación.
El segundo capítulo, presenta un análisis de los efectos de la humedad en la transferencia de calor
del evaporador, se desarrolla un modelo de regresión para obtener los coeficientes de transferencia
de calor y de transferencia de masa del aire húmedo en el evaporador
El tercer capítulo describe el desarrollo de un modelo dinámico tipo “caja gris” del evaporador del
sistema de DX, trabajando con R-410A como refrigerante y analiza el fenómeno del cambio de
fase del refrigerante bajo las condiciones típica de operación en la región tropical.
Por último, se desarrollan las conclusiones y recomendaciones en el último capítulo.
Introducción
Justificación
Los sistemas de aire acondicionado ya no son de uso exclusivo de las industrias y de los edificios
comerciales, los costos unitarios de estos equipos han disminuido, y su demanda ha tenido un
crecimiento exponencial a lo largo de las últimas décadas [1]; aunque los usos de tecnologías
eficientes han mejorado el consumo energético; se estima que la demanda energética es mucho
mayor a la esperada. La agencia internacional de energía (IEA), informa que las ventas anuales de
estos equipos se han cuadriplicado desde 1990 hasta 2017 para alcanzar 135 millones de unidades
a nivel mundial, esto representa aproximadamente 2000 TWh (terawatts horas) de consumo
eléctrico al año. La agencia advierte que el uso excesivo de estos equipos ha causado un aumento
en la demanda de energía eléctrica y, por consiguiente, el aumento de las emisiones de CO2, bajo
este ritmo, es muy probable que entremos a una inminente “crisis del frio”[2]. Uno de los grandes
causantes de esta alta demanda energética, se debe a países de clima tropical cuyo
aprovechamiento energético en general es ineficiente y el consumo energético en sistemas de aire
acondicionado es muy elevado. La Figura 0-1, presenta la demanda a nivel mundial de los aires
acondicionado y su crecimiento por sector [3]. Se observa la demanda de los equipos orientada
hacia el nivel doméstico que a nivel comercial/industrial. Se observa también que las regiones de
mayor demanda de estos equipos corresponden a China, Asia, América del Norte y Japón.
Figura 0-1 Demanda de unidades de aire acondicionado residencial y comercial a nivel mundial desde 2011-2016.
Fuente [3]
La Figura 0-2 representa la distribución de la población mundial por longitud y latitud en el 2015
presentado por Datagraver [4]; el color representa el tipo de clima que presenta la región (cálido –
rojo, templado – verde, frío – azul) y los ejes horizontales y verticales representa la suma de la
población por latitud y longitud. Se puede observar que alrededor del 60% de la población mundial
vive en regiones tropicales; estas regiones se caracterizan por tener alta temperatura y alta humedad
además de tener una variación mínima anual de la temperatura y la humedad. Ciudades con alta
densidad poblacional como las ciudades de Mumbai (India); Karachi (Pakistán) Manila (Filipinas),
Dhaka (Bangladesh), Sao Paulo (Brasil), Lagos (Nigeria), Yakarta (Indonesia), Taipéi (Taiwán) y
Kinshasa (Congo) se encuentran ubicada en zona tropical, por lo tanto, se puede relacionar la
demanda de aire acondicionado con el alto consumo energético de estas unidades.
Figura 0-2 Distribución de la población mundial por latitud y longitud. Fuente [4].
Esta correlación es más evidente al identificar lo que se informa en el índice de trilema energético
de la WEC [5] (Figura 0-3) . Este indicador es presentado con el propósito de clasificar el uso
eficiente de los recursos energéticos de 125 países evaluados. Al analizar los países mencionados
anteriormente Brasil es el país de mejor posición (53 en el segundo Cuartil) y la república del
Congo el de peor posición (122 en el cuarto Cuartil). La mayoría de los países se encuentran en el
tercer cuartil, se puede decir entonces que todos estos países de alta densidad poblacional que se
ubican en zona tropical se pueden clasificar como ineficientes en su el uso energético, y es
necesario implementar lineamientos en estas regiones que puedan evitar esta crisis del frío [2].
Figura 0-3 índice trilema del consejo mundial de la energía. fuente [5]
En Colombia, a pesar de no tener una alta densidad poblacional en sus zonas tropicales, y tener un
índice trilema similar al de Brasil (posición 49) presenta estos mismos inconvenientes de uso
ineficiente energético; según JRAIA (siglas en inglés de Asociación de Industrias Japonesas de
Refrigeración y Aire Acondicionado), la demanda de aire acondicionados pasó de 176000
unidades en 2012 a 286000 unidades en 2017, un crecimiento del alrededor 60% en 5 años, de
estas unidades, 263000 en 2017 corresponde a unidades domésticas [3]. Según ACAIRE, el 60%
del consumo energético que se presenta en los edificios de la región tropical corresponde al uso de
sistemas de acondicionamiento de aire y es necesario implementar técnicas de ahorro energético
para reducir este consumo.[6]
La EIA estableció varias oportunidades estratégicas para ahorrar energía en sistemas de
acondicionamiento de aire; una de estas líneas de oportunidades es el obtener un aumento de la
eficiencia energética de los equipos de aire acondicionado. Las investigaciones de mejora en la
eficiencia energética de estos equipos pueden desarrollarse a través de varios enfoques; Chua et al
[7] agrupó estos enfoques en tres grandes grupos: 1) mejoras de los equipos principales (mejoras
tecnológicas), 2) mejoras innovadoras (mejoras en el diseño), 3) Implementación de los sistemas
de control (mejoras operacionales). Las investigaciones asociadas al grupo 1 corresponden a
trabajos que se realizan para incrementar la eficiencia a partir de mejoras en los dispositivos que
hacen parte del sistema de acondicionamiento de aire (evaporador, condensador, compresor, y
dispositivos de expansión). En este grupo pertenecen los compresores tipo scroll, tipo tornillo,
evaporadores de enfriamiento evaporativo, válvulas de expansión electrónica etc. Las
investigaciones asociadas al grupo 2 corresponden al uso de sistemas alternativos para la
realización de procesos de acondicionamiento de aire, como es el caso de utilización de energías
limpias para el funcionamiento de los equipos; cogeneración, desarrollo de distritos de
enfriamientos, entre otros. El último grupo corresponde a la implementación de estrategias de
control que permitan una óptima operación de los sistemas de acondicionamiento de air, de este
grupo cabe destacar la inclusión de los equipos inverter, los cuales ayudan a ajustar el desempeño
bajo ciertas condiciones de temperatura, regulando el consumo energético del compresor; como
consecuencia, un mayor ahorro de los equipos acondicionado en comparación a un equipo sin
inverter.
Al estudiar más a fondo lo presentado anteriormente, podemos destacar; los países mencionados
se localizan en zona de clima tropical (cuya temperatura media se encuentra entre 24 – 30°C y la
humedad relativa promedio está entre 66 y 96%), condiciones altísimas a lo largo del año. El
objetivo de un sistema de aire acondicionado es proveer condiciones adecuadas para el “confort
térmico humano” y estar satisfecho con el ambiente térmico que lo rodea desde un punto de vista
físico, psicológico y fisiológico [8]. En otras palabras, desarrollar la operación adecuada para que
el clima del espacio habitado no afecte la actividad realizada por las personas; esto se hace
controlando 3 de los siguientes 6 parámetros: temperatura, humedad, velocidad del aire,
temperatura media radiante, nivel de actividad y vestimenta. Todos los equipos de
acondicionamiento de aire actuales, de acuerdo a nuestro mejor conocimiento, se enfocan en
controlar los tres primeros parámetros, por consiguiente, es importante que un sistema de
acondicionamiento de aire, incluyendo equipos para uso residencial de expansión directa, tenga la
capacidad de mantener la temperatura y humedad del aire en el espacio bajo condiciones aceptadas
de confort. Considerando todo lo anterior, es coherente suponer, que los equipos de
acondicionamiento de aire tendrán un alto consumo energético (en comparación a un equipo que
trabaje en zona templada) esto es porque el compresor debe trabajar más; debe comprimir más
refrigerante para remover la carga térmica del espacio, y también debe aumentar su presión de alta
para que pueda trabajar correctamente el condensador debido a las condiciones ambientales,
causando una menor eficiencia en el sistema [7], [9], [10] este trabajo parte de este problema; la
necesidad de mejorar la eficiencia energética de los sistemas de aire acondicionamiento de aire en
climas tropicales. Buscando condiciones que permitan usar el refrigerante de la mejor manera para
enfriar y deshumidificar.
Para afrontar los efectos de clima tropical y ser más eficientes, los equipos de acondicionamiento
de aire tipo comercial e industrial (es decir de grandes superficies y capacidad como Chillers,
Sistemas Centrales, etc.) incluyen otros dispositivos, como es el caso de sistemas desecantes
(sólidos o líquidos) para remover la humedad del aire[11], [12], sistemas de deshumidificación
mecánica [13] o implementar estrategias de control más agresivas [10], [14]. Sin embargo, cuando
queremos implementar esto a sistemas más sencillos como los usado a nivel residencial, se
observan problemas asociados a limitaciones geométricas o aumento en el costo inicial del equipo
[15]. Teniendo en cuenta todo lo anterior, este trabajo se enfocará en el análisis energético de un
sistema de acondicionamiento de aire tipo residencial que opera en climas tropicales; con el
objetivo de encontrar oportunidades de mejoras energéticas y asegurando condiciones adecuada
de confort térmico
Se ha considerado en este trabajo realizar de un diseño de experimentos; hemos discutido
previamente que los equipos de acondicionamiento de aire por expansión directa (DX ) son los
sistemas más usados a nivel residencial, comparando estos equipos con los equipos de aire
centralizado, o equipos de acondicionamiento de aire tipo agua-aire, estos equipos tienen una
configuración más sencilla, pueden alcanzar un aceptable rendimiento energético y son fáciles de
instalar, operar y mantener, convirtiéndolos en sistemas adecuados para nivel residencial y para
edificios de mediana y pequeña capacidad[8], pero, debido a su simplicidad estos equipos no son
capaces de mantener una temperatura de control estable, además no es capaz de realizar el proceso
de deshumidificación y enfriamiento por separado [16]; incluir la tecnología desecante y/o de
control avanzado puede presentar problemas porque se debe incluir más dispositivos en el
evaporador, y esto ocasionaría aumento el costo inicial del equipo y aumentaría su complejidad
operativa. Otro inconveniente como se mencionó anteriormente es la separación entre eficiencia
energética y deshumidificación en equipos DX; la operación de equipos DX de alta eficiencia han
reportado inconvenientes en regiones de clima tropical. La Figura 0-4 representa el efecto de
intercambio de calor entre el aire y el refrigerante en el evaporador de un equipo residencial de
alta eficiencia; la parte derecha de la figura indica cuando el aire se encuentra a una alta
temperatura (28-35°C), el flujo másico del refrigerante (y su temperatura) es el adecuado para
reducir la temperatura del aire a condiciones de diseño (10-13°C); además, permite una
deshumidificación a lo largo del evaporador a través de un área efectiva de deshumidificación.
Cuando la temperatura del aire que ingresa al evaporador disminuye, el equipo de alta eficiencia
reduce gradualmente el consumo energético del compresor, causando una disminución del flujo
másico del refrigerante que va hacia el dispositivo de expansión. Esto se debe al lazo de control
del sistema, el cual consiste en regular el consumo energético del sistema midiendo la temperatura
del aire a la entrada, comparando con la temperatura de ajuste y controlando con un inverter el
trabajo eléctrico del compresor.
Figura 0-4 Intercambio de calor entre el aire y el refrigerante dentro del evaporador de un sistema DX de alta
eficiencia (Izq. Condición de alta temperatura. Der. Condición de baja temperatura). Fuente: elaboración propia.
Al disminuir el consumo energético, disminuye el trabajo compresión y aumenta la presión de baja
del compresor; esto genera un aumento de la temperatura del refrigerante en el evaporador como
se observa en la parte izquierda de la Figura 0-4. El proceso de transferencia de calor con fines de
remoción de calor sensible no se verá afectado significativamente por la disminución del flujo
másico del refrigerante, aunque el aire saldrá con una mayor temperatura seguirá cumpliendo el
objetivo de disminuir temperatura. Sin embargo, la condensación sí se verá afectada de manera
significativa; principalmente porque una mayor temperatura del refrigerante ocasionará que este
último se evapore más rápidamente debido a la disminución de la zona de cambio de fase en el
domo, lo que disminuye el área de zona bifásica ocupada por el refrigerante y disminuirá el área
efectiva de deshumidificación. Además, al disminuir el coeficiente de transferencia de calor,
disminuirá el coeficiente de trasferencia de masa. Al disminuir el coeficiente y área de
transferencia de masa afectará directamente el proceso de deshu midificación dentro del
evaporador. Xu et al, manifiesta otro inconveniente relacionado; las tendencias actuales de diseño
de evaporadores para sistemas residenciales están orientadas a mejorar los índices de eficiencia
energética, esto se consigue aumentando el área frontal y reducido la profundidad del evaporador,
con esto permiten que la temperatura del refrigerante dentro del evaporador sea más alta,
reduciendo el trabajo del compresor, sin embargo, esto ha deteriorado el proceso de
deshumidificación a nivel operativo, ya que presenta los mismos inconvenientes explicados
anteriormente, específicamente al no respetarse el factor de calor sensible (SHR), en el diseño del
evaporador por deshumidificación [15].
En consecuencia de todo lo anterior, se presenta una condición de “no confort” (baja temperatura,
alta humedad) y para solucionar esta condición, lo que comúnmente se realiza es disminuir la
temperatura de control de operación del equipo; la disminución permitirá que el equipo trabaje
más tiempo permitiendo mayor deshumidificación, sin embargo, conllevará a un mayor consumo
energético y una temperatura inadecuada en el sistema, este fenómeno se llama
“sobreenfriamiento”[8].
La tendencia actual de diseño de sistemas de acondicionamiento de aire tipo residencial, ha tenido
que sacrificar la condensación del vapor de agua presente en el aire para mejorar el consumo
energético; autores han expresado que la mejor opción para que un sistema de aire acondicionado
en climas tropicales pueda operar eficientemente es desacoplar los procesos de enfriamiento y
deshumidificación [15], esto ha ocasionado estudios orientados a diseño de control múltiple de
temperatura y humedad en el aire; como el estudio presentado por Deng y Li [9] donde desarrollan
una estrategia de control digital para controlar simultáneamente la temperatura y la humedad del
aire que ingresa a una unidad de aire acondicionado de expansión directa a partir del control del
factor de calor sensible y un algoritmo de computación numérica que permitió un ajuste del control
y la energía más eficiente que un sistema sencillo; uno de los inconvenientes de este control se
debe a su elevado costo de fabricación, lo que elevaría considerablemente el costo de la unidad de
aire acondicionado.
Otra investigación, realizada por Qi Qi et al [16], diseño una estrategia de control MIMO para
mantener simultáneamente la humedad y la temperatura del aire de un recinto; el control de
múltiples entradas y múltiples salidas (MIMO) se mantenía gracias a la manipulación de la
velocidad simultánea de la velocidad del compresor y la velocidad del ventilador en el evaporador.
Diseñado bajo la técnica Gaussiana Lineal Cuadrática (LGQ), los resultados demostraron que el
sistema de control puede manejar desde pequeñas hasta grandes perturbaciones de calor que
ingresa al espacio acondicionado. No obstante, el diseño del sistema de control fue realizado en
un modelo dinámico linealizado, esto es un inconveniente porque sus resultados se aplicarían
únicamente a una condición específica. Sin embargo, sus resultados presentan un aporte
significativo para la investigación de control de aire acondicionado en climas tropicales.
Por otra parte, hay estudios que demuestran que el ahorro energético también se puede alcanzar a
través de la optimización. Huh y Brandemuehl presentan los resultado de la optimización en un
sistema de control para un equipo de acondicionamiento de aire tipo DX con el fin de minimizar
el consumo energético conservando las condiciones de confort, especialmente manteniendo
control de humedad; los resultados demuestran que sí es posible obtener ahorros energéticos a
través de la optimización del sistema de control y evitando a su vez los problemas de
sobreenfriamiento [10]. Aunque cabe resaltar, que el estudio fue realizado con un equipo tipo
paquete con ductería, que permite un mejor control independiente, ya que utiliza un plenum para
control de aire (elementos que no se incluyen en un equipo de acondicionamiento de aire para uso
residencial), pero a su vez, demuestra que es posible identificar oportunidades de ahorro
simplemente a través de la optimización del control operacional del equipo.
Un segundo estudio planteado para este trabajo es la de analizar el efecto de transferencia de calor
y masa entre el aire de clima tropical y un evaporador DX de un sistema residencial. Para esto,
hemos realizado una búsqueda sobre coeficientes de transferencia de calor y transferencia de masa
en intercambiadores de placa y tubos de pequeña capacidad en climas tropicales. Se identifica
inicialmente que los trabajos realizados por Kays y London [17] McQuiston [18], Threlkeld [19]
son los más utilizados sin tener en cuenta las condiciones del clima, como lo presenta Qasem y
Zubair en [20]. Cabe destacar que los modelos de correlación calculado experimentalmente son
bajo ciertas condiciones de laboratorio y con una gran variedad de intercambiadores con diferentes
geometrías.
Al estudiar más a fondo se identifica los trabajos desarrollado por Threlkeld [19] y Wang [21]–
[23] donde incluye la deshumidificación en sus investigaciones de cálculo del coeficiente de
transferencia de calor al incluir la eficiencia de la aleta para aire húmedo. Wang incluye además
en su investigación la correlación de transferencia de masa. Recientemente la investigación ha sido
avanzada por Kim [24], [25], Yang [26], [27] y Zhang [28], en cada uno de estas investigaciones
se han basado teniendo en cuenta la metodología de Threlkeld para el cálculo del coeficiente de
correlación de Colburn (factor j) [19]. Los resultados indican controversias respecto a la influencia
de la humedad relativa en el cálculo de los coeficientes de transferencia de calor y masa; para
algunos autores, la humedad relativa no afecta [18], [24], para otros sí [21], [26]. No obstante,
todos los autores afirman que un intercambiador de calor con tubos aleteados húmedos (con agua
condensada en las paredes) mejora su coeficiente transferencia de calor en comparación con el
intercambiador con aletas secas (el aire no condensa el vapor de agua) en casos donde el flujo de
aire tiene un número de Reynolds bajo (Re<2000) [18], [22]. Analizando todo lo anterior, se infiere
que el estudio a realizar en este trabajo, no se considerará asumir un coeficiente de correlación
para el análisis de transferencia de calor y masa; deberá calcular y validar un coeficiente ajustado
de transferencia de calor y de masa que permita analizar los efectos de enfriamiento y
deshumidificación del aire húmedo en un sistema de acondicionamiento de aire tipo residencial.
Un último estudio a plantear en el trabajo es el desarrollo de un modelo dinámico; bajo la misma
línea de identificación en la operación de un sistema de acondicionamiento de aire, es necesario
investigar la dinámica del refrigerante dentro del evaporador del sistema. El motivo de analizar el
evaporador como un proceso dinámico se debe a la naturaleza del intercambio de calor y masa,
relacionado a su vez por las fluctuaciones de temperatura y flujo másico del refrigerante. Como se
mencionó anteriormente, un sistema de acondicionamiento de aire tipo DX es controlado a partir
de la temperatura de ingreso del aire al evaporador. Por este motivo, durante una operación es muy
común observar la variación de la temperatura del aire que ingresa al evaporador a lo largo de la
operación del equipo, por lo que el flujo másico y la temperatura en consecuencia se regula para
garantizar un trabajo conveniente de transferencia de calor. Si el equipo contiene un sistema de
control orientado al ahorro de energía como un inverter, el refrigerante y el evaporador tendrá un
desempeño y operación que no es posible evaluar desde una modelación estacionaria.
Otra consideración se debe a la naturaleza del refrigerante; en intercambiadores de calor de
múltiples fases, existe la complejidad del cambio de fase del refrigerante a lo largo del
intercambiador. El refrigerante cambia su estado durante la trayectoria de estos dispositivos, y el
punto donde deja de cambiar de estado no está definido físicamente; debido a esto, la complejidad,
y desarrollo del modelo es mucho mayor que en intercambiadores de calor con refrigerantes sin
cambio de fases (agua, gas, etc.). Es importante detallar estas regiones ya que la transferencia de
calor, la deshumidificación del aire, y la capacidad de enfriamiento del intercambiador se verá
afectado por las condiciones termodinámicas y la fase en la que el refrigerante se encuentra durante
su recorrido a lo largo del intercambiador; la capacidad del refrigerante para absorber el calor es
mucho mayor fase saturada que en fase gas sobrecalentada; a su vez, la condensación del agua en
el aire es mejor cuando el refrigerante se encuentra en fase saturada, debido que el refrigerante
permanece a una temperatura de saturación (que por lo general se encuentra por debajo de la
temperatura de rocío del aire húmedo); en cambio, si el refrigerante está en fase de gas
sobrecalentada es posible que se encuentre a una temperatura mayor a la temperatura de rocío del
aire y por lo tanto no ocurrirá condensación en el aire húmedo. A partir de la modelación dinámica,
es posible identificar el área ocupada por el refrigerante a lo largo del intercambiador de calor en
su fase de mezcla saturada y en su fase de gas sobrecalentado, permitiendo así establecer algunas
estrategias operacionales o incluso desarrollar estrategias de control avanzado.
Como última consideración, el evaporador trabaja bajo la naturaleza térmica del proceso, su
variable dominante es la temperatura y presión del refrigerante, y debido a esto, durante cualquier
cambio en las condiciones de entrada del evaporador, el tiempo de respuesta de estos dispositivos
tienden a ser más lentos que otros dispositivos que hacen parte del sistema de acondicionamiento
de aire, (como el compresor) y debe ser analizado como un modelo dinámico. La estructura de un
modelo dinámico se presenta en la Figura 0-5 y es necesario identificar estas variables:
 Variables de entradas: Estas son variables conocidas que actúan en el sistema, existen dos
tipos de variables de entradas: variables controladas por el experimento (flujo másico del
refrigerante) y variables no controladas (entalpía del aire a la entrada del condensador)
 Estructura del sistema y parámetros/propiedades: Son variables que nos ayudan a describir
físicamente el sistema (propiedades mecánicas, masa del evaporador, etc.)
 Variables de salida: Variable de respuesta o dependiente que describe la reacción del
sistema a partir de las variables de entrada. Por lo general la energía es una variable de
salida.
Figura 0-5 Esquema para un modelo dinámico de un intercambiador de calor. Fuente: [8]
Rasmussen [29], [30] clasificó los modelos dinámicos de intercambiadores de calor que funcionan
con refrigerante de cambio de fase en tres categorías: modelos de parámetros concentrados (LP),
modelos de límites movibles (MB) y modelos de volumen de control fijos (FCV), esta
categorización se realizó según el tiempo de modelo computacional y precisión obtenida como se
observa en la Tabla 0-1.
Tabla 0-1 Comparación entre categorías de modelos dinámicos para intercambiadores de calor.
Tipo de Modelo Carga computacional Precisión
LP Bajo Menor
MB Moderado Moderado
FCV Alto Mayor
Fuente: elaboración propia

Los modelos LP son los modelos computacionales más simple, y rápidos (eran lo más utilizados
hace más de 20 años, cuando los equipos no tenían mucha capacidad de carga computacional) sin
embargo imprecisos para el desarrollo de estrategias de control avanzadas [30]. Se caracterizan
porque las propiedades termodinámicas del refrigerante dependen únicamente del estado y no se
tienen en cuenta consideraciones asociadas al cambio de fase [31]. Por lo general evalúan la
modelación de un intercambiador de vapor a través de métodos de cálculo rápido como los
métodos LMTD y/o NTU; o a través de métodos detallados como el método de parámetros
distribuidos o el método de zona [30]. Estos modelos son ideales para describir y predecir la acción
que ocurre por fuera del evaporador.
Por otra parte, los modelos FCV desarrollan una mayor exactitud de sus resultados en comparación
de las otras categorías [30]; este tipo de modelos se caracteriza por utilizar las ecuaciones de
conservación de energía, masa y continuidad y aproximarlas a partir del método de diferencias
finitas a lo largo del evaporador [31]. El evaporador a su vez es dividido en volúmenes de control,
entre más volúmenes mayor precisión. Software de CFD y paquetes de simulación computacional
de flujos utilizan estos modelos, aprovechando la capacidad computacional y la utilización de
nodos computacionales. La desventaja de estos modelos es su alto costo computacional necesario
para garantizar independencia entre los valores predichos y la discretización [30].
Los modelos MB, también conocido como modelo de limites movibles, (Figura 0-6), son modelos
que permiten una mayor rapidez respecto a un modelo FCV y una mejor precisión en comparación
con los modelos LP. Los modelos MB se caracterizan por dividir el intercambiador de calor en
zonas (regiones) definidas por la fase del refrigerante (líquido, vapor, bifásico), luego en cada
región se evalúa bajo parámetros promedios por cada fase del fluido y con las ecuaciones de
conservación de energía, masa y momento [31]. Este tipo de modelo al ser más rápido que los
modelos FCV permite un buen ajuste para el desarrollo de estrategias de control o análisis de
intercambiadores [16], [31]–[36]. Este enfoque se llama así porque conforme las variables de
entradas son modificadas, también se modifica las regiones del refrigerante, así que el método
realiza los cálculos con base a la nueva región generada.
Figura 0-6 Esquema del modelo de límites movibles. Fuente: elaboración propia
Conclusiones Introducción
Teniendo en cuenta todo lo anterior, este trabajo presentará una investigación agrupado en la
sección de mejoras operacionales, para identificar mejoras en la eficiencia energética a partir de
estrategias de control y de identificación del proceso. Para esto, este trabajo se compone de tres
estudios que permitirán identificar el comportamiento energético y operativo de un equipo de
acondicionamiento de aire tipo residencial que trabaja en condiciones tropicales.
El primer estudio, corresponde a la caracterización operacional de un sistema de
acondicionamiento de aire tipo residencial bajo diferentes condiciones de entrada del aire a partir
de un diseño de experimento, este se hace con el fin de establecer una relación entre las condiciones
psicrométricas del aire a la entrada del evaporador, y la eficiencia energética del sistema DX.
El segundo estudio, busca analizar la relación entre coeficiente de transferencia de calor,
coeficiente de transferencia de masa y la operación del equipo, a fin de relacionar la
deshumidificación a través de la condensación en presencia de un gas no condensable (aire seco)
con la eficiencia del sistema y establecer una correlación entre estos.
El tercer estudio desea establecer un modelo dinámico del evaporador, a fin de investigar el
comportamiento del refrigerante, durante el proceso de intercambio de calor y de
deshumidificación del aire en el evaporador, a fin de identificar la región donde el refrigerante
cambia de fase y por consiguiente el proceso de intercambio de temperatura y transferencia de
masa sea modificado, ya que el flujo másico es finito.
Por último, se destaca que durante la búsqueda de la información, los estudios sobre sistemas de
acondicionamiento de aire en climas tropicales encontrados fueron desarrollados en la región
subtropical del continente asiático [7], [15], [37]–[44]; y reportes para clima templado durante la
estación de verano [10], [15], a nivel del país, los trabajos que se encontraron como el de Bula y
Méndez, [45] Fajardo et al [46] Balbis et al [47], se desarrollan en unidades de aire acondicionado
para clima tropical, sin embargo, estos estudios investigan el comportamiento de los sistema de
acondicionamiento de aire, sin considerar los efectos del aire húmedo propio de la región tropical
colombiana.; con este trabajo se pretende aportar información necesaria sobre la caracterización y
evaluación energética de un sistema de acondicionamiento de aire tipo residencial bajo el clima
tropical de la región.
Planteamiento del problema
¿Existe una oportunidad de identificar nuevas estrategias de ahorro energético en sistemas de aire
acondicionado tipo residencial, teniendo en cuenta la afectación directa entre el desempeño
energético y las condiciones ambientales del aire tropical colombiano caracterizado por su alta
temperatura y humedad?
Objetivo General
Estudiar la influencia del aire húmedo en el desempeño de un equipo de aire acondicionado tipo
residencial, con el propósito de establecer un punto de partida en el desarrollo de nuevas
oportunidades de ahorros energéticos.
Objetivos específicos
 Analizar energéticamente un sistema de acondicionamiento de aire bajo las condiciones
del clima tropical.
 Evaluar los efectos de la humedad y condensación del aire en los fenómenos de
transferencia de calor, masa, y desempeño energético de un sistema de aire acondicionado
 Obtener de coeficientes de transferencia de calor y masa en un evaporador de tubos
aleteados bajo sus condiciones operacionales y las condiciones climáticas del aire húmedo
de regiones tropicales.
 Diseñar un modelo dinámico del evaporador de la unidad de aire acondicionado para
analizar el efecto del refrigerante y sus parámetros termodinámicos cuando es sometido al
aire húmedo típico de regiones tropicales.
1 Análisis energético del aire típico de un clima tropical sobre un
aire acondicionado de uso residencial
Como se expresó previamente, la tendencia en diseño de evaporador para equipos de aire
acondicionado de uso residencial está orientadas para disminuir el consumo energético
aumentando el área frontal del evaporador (priorizando el enfriamiento) y disminuyendo las filas
del evaporador, disminuyendo el proceso de la deshumidificación. A niveles de clima templado
(o estacional), este tipo de diseño es funcional; pero en una región de clima tropical, la
deshumidificación es necesaria para garantizar condiciones de confort dentro del espacio a
acondicionar. Con base a lo anterior, se desea evaluar la influencia del clima tropical en el proceso
de enfriamiento y deshumidificación de un sistema de aire acondicionado a nivel residencial. Para
esto, se propuso realizar un diseño experimental, en el cual se modifica las condiciones
psicrométrica del aire (temperatura y humedad) que ingresa a un sistema de acondicionamiento de
aire. Con estas modificaciones se busca analizar los parámetros energéticos del sistema, su
eficiencia energética y el fenómeno de transferencia de calor/masa que se presenta en el sistema
experimental, con el propósito de caracterizar y establecer oportunidades de mejora en el
enfriamiento y deshumidificación del aire de clima tropical a nivel residencial.
1.1 Características del aire en el clima tropical

El aire de clima tropical corresponde a cualquier región del mundo que se encuentran ubicadas
entre el trópico de cáncer (23° latitud norte) y el trópico de capricornio (23° latitud sur). Se
caracteriza también porque la temperatura media del aire es superior a 18°C y su humedad relativa
también es elevada [48]. El clima tropical a su vez se subdivide en tres regiones (Ecuatorial,
Monzónico, y de sabana) [49], en el caso de Barranquilla se encuentra en una región tropical de
sabana, una región que se caracteriza por mantener un ciclo de temperatura muy alto a lo largo del
año. Según el IDEAM el clima de Barranquilla y sus alrededores tiene una precipitación al año de
822 mm, la temperatura promedio es de 26.9 °C y con una temperatura máxima entre 31 y 32 °C;
la temperatura mínima se presenta en las madrugadas con una oscilación entre 23 y 25 °C. Por otra
parte, la humedad relativa oscila entre 75-85% a lo largo del año.[50].
Para esta investigación, las condiciones psicrométricas del aire a la entrada del evaporador deben
ser similares a la de un local residencial ubicado en un clima tropical, como referencia se utilizará
las indicadas en la ciudad de Barranquilla, además, a fin de caracterizar adecuadamente la
temperatura y la humedad del aire en clima tropical, se considera utilizar otras condiciones
(niveles) de temperatura y humedad con el fin de analizar individualmente la influencia de estos
parámetros y eliminar cualquier tipo de inferencia en el estudio.
1.2 Diseño de Experimentos

Con base a todo lo mencionado en la sección anterior, se propuso realizar un diseño de
experimento, orientado hacia las condiciones del aire a la entrada del evaporador. Uno de los
grandes retos a enfrentarse durante el diseño de experimento en sistemas de aire acondicionado,
es la selección del tipo de diseño a considerar. A continuación, se presentan los detalles más
importantes del diseño de experimento. Para profundizar en las características deseadas en el
diseño de experimentos y sus justificaciones remitirse al Anexo “Selección del diseño de
experimento”.
1.2.1 Factores
Los factores a considerar en el diseño son la temperatura y la humedad del aire a la entrada del
evaporador, el último factor a considerar es el nivel de velocidad del ventilador del evaporador.
Teniendo en cuenta lo explicado en el anexo y las recomendaciones de [51], [52], un diseño
factorial 2k ; es el escogido para este estudio.
1.2.2 Variables de respuesta
Las variables de respuestas para este estudio fueron tres: el flujo de calor que el evaporador absorbe
del aire (Qevap), la energía requerida para que se realice la operación del equipo (Wcomp), y
finalmente el coeficiente de desempeño (COP) [8]. El COP es calculado como se presenta en la
ecuación (1-1).
QEvap
COP  (1-1)
Wcomp
QEvap, y Wcomp, se obtiene a partir de las ecuaciones (1-2) y (1-3).
 aire haire,in  haire,out 
QEvap  m (1-2)
V I (1-3)
Wcomp 
1000
Para el cálculo del flujo másico del aire, se usa la ecuación (1-4), a partir de la medición de
velocidad frontal del aire, la densidad y el área de sección transversal. La entalpía del aire se
obtiene calculando a partir de la relación temperatura de bulbo seco y humedad absoluta como se
presenta en la ecuación (1-5):
 air    U air  ACS
m (1-4)
hair  1.006  Taire  waire  2501  1.86  Taire  (1-5)

Los resultados de la ecuación (1-2) y ecuación (1-3) son comparados con el calor absorbido,
ecuación (1-6), y el trabajo del compresor, ecuación (1-7), realizado desde el lado del refrigerante
para verificar su validez.
QEvap _ ref  m ref href ,out  href ,in Evaporador (1-6)
WComp _ ref  m ref href ,out  href ,in Compresor (1-7)

El COP que se presenta en la ecuación (1-1), indica la relación entre las energías que ingresan al
sistema DX. Generalmente el COP representa la relación interna entre el calor del evaporador y el
trabajo del compresor, es decir, la variación de entalpía del refrigerante en el evaporador y en el
compresor como se muestra en la ecuación (1-8).
QEvap _ ref
COP  (1-8)
Wcomp _ ref
Por último, para calcular las propiedades psicrométricas del aire se utilizaron las ecuaciones que
se presentan en el manual de “Fundamentals” de ASHRAE [8]. Para la entalpía y otras propiedades
del refrigerante son obtenidas a partir del software de propiedades termodinámicas REFPROP
basándose en las condiciones de temperatura y presión de las mediciones del sistema.
1.2.3 Niveles factor
Los niveles para cada factor experimental se presentan en la Tabla 1-1 y Figura 1-1, una
explicación más detallada se observa en el anexo I. Como se observa en la Figura 1-1 solo uno de
los niveles experimentales (zona naranja) se encuentran en la zona de confort recomendada para
climas cálidos (zona azul), el resto representa condiciones que fácilmente podemos encontrar el
aire del recinto a acondicionar.
Tabla 1-1 Factores y niveles del diseño de experimento
Factores y niveles
Código Nivel
Código Rango
Factor Factor Unidades Tipo
Bajo Alto Bajo Alto Aceptación
A: Temperatura °C Continuo 23 28 -1 1 +0.5°C
-1 1 + 1% HR (ver tabla
B: Humedad Absoluta g/kg Continuo 11 14
4)
Velocidad del 1 2
C: Nivel Categórico Bajo Alto
ventilador

Teniendo en cuenta que el sistema experimental utiliza sensores de humedad relativa, los niveles
de humedad absoluta fueron evaluados bajo condiciones de humedad relativa como se observa en
la Tabla 1-2. Ver anexo donde se profundiza sobre el tema.
Tabla 1-2 Nivel de humedad relativa según niveles de los factores
Humedad relativa según temperatura y humedad absoluta

Humedad absoluta: (-1) Humedad absoluta: (1)
Temperatura (-1) 62% 75%
Temperatura (1) 46% 55%

1.2.4 Factor bloque
Se incluye como factor de bloque el día de la experimentación; se realizaban 2 corridas
experimentales por día, dando como resultados 4 bloques (ver Anexo I).
1.3 Descripción del sistema.
Para desarrollar esta investigación, se construyó un sistema experimental. En la Figura 1-2 se

presenta el sistema utilizado para este experimento se encuentra localizado en la Universidad del
Norte - Barranquilla, en el laboratorio de Termofluidos. del equipo funciona con un voltaje
nominal de 115V monofásico y está etiquetado como eficiencia D bajo el Reglamento Técnico de
Etiquetado (RETIQ), esto quiere decir, que la eficiencia energética de este equipo es de 3.14 W/W.
Figura 1-1 Carta psicrométrica con los niveles de los factores (Zona naranja) y la zona de confort para verano (zona
azul). Fuente: elaboración propia
El sistema consiste principalmente de un aire acondicionado comercial tipo mini split con inverter,
3.5 kW (1TR) de capacidad térmica y un consumo energético de 133.98 kWh/mes; la unidad opera
con refrigerante R410A, el compresor Además de la unidad de aire acondicionado, el sistema
contiene una serie de sensores instalados a lo largo del equipo para medir las variables del
refrigerante, aire y consumo eléctrico de la unidad. Este sistema fue diseñado para esta tesis
doctoral, a partir de un presupuesto otorgado por COLCIENCIAS en la convocatoria 727 de 2015.
Figura 1-2 Izq. Sistema Experimental Der. Etiqueta energética. Fuente: elaboración propia
El esquema de este sistema experimental se presenta en la Figura 1-3. Este sistema fue construido
con el objetivo de analizar el comportamiento de una unidad de aire acondicionado tipo mini split
DX bajo condiciones típicas de operación. Para mantener una condición de operación real, el
sistema no sufrió ninguna modificación para controlar su operación,
El sistema experimental contiene 21 sensores distribuidos a lo largo de todo el sistema, que
permiten analizar las propiedades termodinámicas del refrigerante, las propiedades psicrométrica
del aire y otras condiciones; está compuesto de 10 sensores de temperatura, 4 sensores de presión,
3 sensores de humedad relativa, 3 sensores de corriente eléctrica y 1 sensor de flujo volumétrico,
El flujo másico del refrigerante es obtenido a partir de las mediciones del flujo volumétrico y la
densidad del refrigerante como se presenta en la ecuación (1-9);
 ref  T 7, p2  F1
m (1-9)
La Tabla 1-3 presenta las especificaciones de los sensores; la nomenclatura y ubicación se
presentan en la Tabla 1-4 . Todos los sensores envían la información a un sistema de adquisición
de datos el cual envía la información ya sea en línea o a través de una memoria externa SD. Este
sistema experimental no busca evaluar la unidad de aire acondicionado, por lo tanto, no se
construye un sistema experimental tipo calorímetro o tipo método de entalpía como se presenta en
[53], con este equipo se desea evaluar el comportamiento energético, y el fenómeno de
transferencia de calor entre el refrigerante y el aire húmedo durante el proceso de intercambio de
calor en la unidad mini split . Sin embargo, durante la toma de mediciones en el sistema
experimental se siguieron los protocolos de pruebas como se establecen en las normas ISO 5151
[53]y AHRI 210 [54] para garantizar estabilidad del procedimiento.
Figura 1-3 Esquema del sistema experimental y ubicación de los sensores Fuente: Elaboración propia
Las condiciones de entrada del aire en el evaporador fueron obtenidos a través de una cámara de
aire donde se mezclaba, el aire del recinto, junto con un sistema generador de calor y vapor; esto
se conseguía gracias a un generador de vapor y una resistencia eléctrica Figura 1-4. Por último, un
anemómetro de hilo caliente fue utilizado para la caracterización de la velocidad del aire a la salida
del evaporador.
Esta medición no fue incluida en los sensores, debido que debía tomarse muchas mediciones en
varias posiciones, a lo largo del evaporador y a lo ancho; las características y precisiones se pueden
observar en Tabla 1-3.
Figura 1-4 Esquema del generador de vapor, resistencia eléctrica y cámara de mezcla para suministrar condiciones
de entrada del aire al evaporador. Fuente: elaboración propia
Tabla 1-3 Descripción de los sensores instalados en el sistemas experimental
Características de los sensores

Sensor Modelo Unidad Rango de Precisión Cant
trabajo
Temperatura + 0.5°C para
Termopar sumergible digital -55°C a
°C -10°C a 10
Maxim Integrated - DS18B20 +125°C
+85°C
Presión 0 a 500
Transductor de presión psi para
psig + 1 psi 4
Full Gauge – SB69 -25°C a
80°C
Humedad Sensor capacitivo de humedad 0 a 100%
Relativa tipo “, thermoset polymer para
% + 1.5% Error 3
capacitive sensing” Marca 0°C a
Honeywell HIH 4000 40°C
Corriente Sensor de corriente 1,5% Error a
A 0 a 5A 3
Allegro ACS712 25°C
Flujo Sensor de flujo tipo 200 a
+ 1.0% Full
microturbina para líquidos mL/min 2000 1
scale
McMillan 102 mL/min
Velocidad 0.15 a 30 1
Termo higro anemómetro + 3% de
m/s
KIMO VT 210 con sonda de hilo m/s lectura
-20 a
caliente + 0.03 m/s
80°C

El funcionamiento del sistema experimental consiste en lo siguiente: según sea el caso a considerar
(evaluar) se utilizaba uno o los dos equipos generadores de temperatura/humedad, para ajustar las
condiciones de temperatura del aire a la entrada del evaporador. La unidad de aire acondicionado
es iniciada y se procede a esperar que el sistema alcance su estado estable. Esto podría ocurrir entre
30 y 45 minutos; al alcanzar su estado estable[54], se inicia el protocolo de adquisición de datos.
Esta fase contaba con la duración de 1 hora y 45 minutos [53], durante el ese tiempo, se verificaba
que las condiciones de entrada del aire al evaporador se mantenían, que la temperatura de salida
del aire no disminuya por debajo de los 14°C y además se caracterizaba la velocidad del aire a la
salida del evaporador. Al finalizar el tiempo, se detiene la toma de datos y el equipo pasa a la
siguiente condición climática de ingreso.
Tabla 1-4 Descripción de la ubicación de los sensores utilizados en el sistema experimental
Ubicación de los sensores

Nomenclatura Descripción Ubicación
T1 Temperatura del refrigerante Entrada del compresor
T2 Temperatura del refrigerante Salida del compresor
T3 Temperatura superficial Compresor
T4 Temperatura del aire Entrada condensador
T5 Temperatura del aire Salida condensador
T6 Temperatura del refrigerante Entrada condensador
T7 Temperatura del refrigerante Salida condensador
T8 Temperatura superficial Tubo capilar
T9 Temperatura del aire Entrada evaporador
T10 Temperatura del aire Salida evaporador
H1 Humedad relativa del aire Cámara de aire
H2 Humedad relativa del aire Salida evaporador
H3 Humedad relativa del aire Entrada evaporador
P1 Presión del refrigerante Entrada del compresor
P2 Presión del refrigerante Salida del compresor
P3 Presión del refrigerante Entrada del compresor
P4 Presión del refrigerante Salida del compresor
I1 Corriente Entrada a la unidad externa
I2 Corriente Ventilador del condensador
I3 Corriente Ventilador del evaporador
F1 flujo volumétrico del refrigerante Salida condensador
1.4 Puesta en marcha y Resultados

El orden de las corridas se presenta en la Tabla 1-5. Para garantizar los resultados experimentales,
inicialmente se evaluó la normalización de los datos medidos calculando la desviación estándar de
las mediciones, con el propósito de verificar que los datos se encontraban entre un rango menor a
3 veces su desviación estándar. Luego se verificó a través balances de energía termodinámico los
resultados obtenidos entre: 1) Calor cedido por el aire vs el calor absorbido por el refrigerante, 2)
Trabajo eléctrico consumido por el compresor, vs trabajo termodinámico realizado por el
compresor, y 3) Coeficiente de desempeño del equipo.
Tabla 1-5 Orden de ejecución de la experimentación.
Orden aleatorio de experimentación

BLOQUE Temperatura Humedad Velocidad ventilador
(Días) °C Relativa % Nivel
1 28 46 Bajo
1 23 75 Alto
2 23 62 Alto
2 28 55 Bajo
3 23 75 Bajo
3 28 46 Alto
4 23 62 Bajo
4 28 55 Alto

Cada corrida experimental representa aproximadamente 2 horas de experimentación. Al finalizar
los datos son procesados en una hoja de cálculo donde se evalúan las propiedades termodinámicas
(Figura 1-5), y luego se calculan las variables de respuestas a partir de las ecuaciones (1-1)-(1-3).
Para verificar los resultados, se utilizan las ecuaciones (1-6) y (1-7) y se compara el porcentaje de
error. Podemos observar en la Figura 1-6, que el calor absorbido por el refrigerante es mayor que
el calor cedido por el aire esto es correcto ya que el refrigerante también absorbe calor de otros
medios. De igual manera en la Figura 1-7 se observa la comparación entre la potencia eléctrica del
compresor y el trabajo termodinámico del mismo calculado a partir de sus propiedades
termodinámicas; estos resultados son coherentes porque el consumo eléctrico del compresor debe
ser mayor al trabajo termodinámico del compresor considerando las irreversibilidades del proceso
de compresión [13].
Figura 1-5 Hoja de Cálculo – Izq. Adquisición de datos. Der. Propiedades termodinámicas y cálculos energéticos.
Fuente: elaboración propia.
Figura 1-6 Comparación: calor cedido por el aire vs calor absorbido por el refrigerante. Fuente: elaboración propia
Figura 1-7 Potencia del compresor termodinámico vs potencia eléctrica del compresor. Fuente: elaboración propia
1.4.1 Calculo tipo
Para obtener las variables de respuestas, el procedimiento es el siguiente:
1. Para el Qevap
a. Se calculó la entalpía del aire a la entrada y a la salida del evaporador a partir de las
condiciones psicrométricas obtenidas en las mediciones
b. Se calculó el flujo másico del aire a partir de la medición de la temperatura y velocidad
promedio del aire durante las medicioens
c. Se aplica la ecuación (1-2)
2. Para el Wcomp
a. Se mide la corriente eléctrica del compresor y de los ventiladores
b. El voltaje es constante durante todo el proceso
c. Se aplica la ecuación (1-3)
3. Para el COP
a. Se aplica la ecuación (1-1)
El cálculo tipo se presenta a continuación
Datos Corrida 1
Aire entrada evaporador:
kgagua
Tin  T 9  28.17C; HRin  HR3  49.16  win  0.01177 kgaire
Aire salida evaporador:
kgagua
Tout  T10  12.47C; HRout  HR2  wout  0.00803 kgaire
Velocidad promedio frontal aire:
U min  1.070 m

m
U  1.940 s   s
m
U max  2.530 s
A partir de estos datos obtenidos se calcula la entalpía del aire a la entrada y salida, densidad del aire y
flujo másico del aire:
kJ
hair _ in  1.006  28.17  0.01177  2501  1.86  28.17  58.398
kgaire
kJ
hair _ out  1.006  12.47  0.00803 2501  1.86  12.47  32.82
kgaire
  1.23 kgaire
m3
ACS  0.059m 2
m air    U air  ACS  mair  1.23  1.940  0.059  0.14148 kgaire s
Con el flujo másico y la diferencia de entalpía se calculará el calor absorbido en el evaporador:

 aire haire,in  haire,out   0.14148 53.398  32.82  3.619kW
QEvap  m
Se busca a partir de los datos de la corrida, la corriente consumida por el compresor para calcular
el trabajo del compresor:
I1  4.17 A I 2  0.122 A I 3  0.170 A
I  I1  I 2  I 3   I  3.878 A
V  120V
V  I 120  3.878
Wcomp    0.466kW
1000 1000
Por último, el COP se calcula con la relación QEvap y WComp
QEvap 3.619
COP    7.763
Wcomp 0.466
Se verifica los resultados de QEvap y WComp con los resultados obtenidos en el ciclo interno del
sistema de aire acondicionado, inicialmente se toman las medidas en diferentes puntos del
sistema:
Salida Refrigerante condensador:
Tref _ out  T 7  24.7C pref _ out  P3  240.80 psig f  F1  1091.04 mL min
Salida refrigerante evaporador
Tref _ out  T1  10.74C pref _ out  P1  79.71psig
Salida refrigerante compresor
Tref _ out  T 2  55.12C pref _ out  P2  262.57 psig
A partir de estas mediciones, con la ayuda de REFPROP, se procede a calcular las propiedades
termodinámicas del refrigerante en base a la presión y temperatura indicada
T 7  REFPROPT7, P3  1061.28 kg

m3
href _ in _ Evap  REFPROPT7, P3  239.30 kJ kg
href _ out _ Evap  REFPROPT1, P1  436.47 kJ kg
href _ out _ Comp  REFPROPT2, P2  458.530 kJ kg
Con las propiedades calculadas se procede a determinar QEvap_ref; Wcomp_ref, y COP_ref

1091.04  1061.28 kg
m ref  f  T 7   0.0193
60000000 s
QEvap _ ref  mref href ,out  href ,in Evaporador  0.0193  436.47  239.30  3.81kW

WComp _ ref  m ref href ,out  href ,in Compresor  0.0193  458.530  436.47  0.4259kW
QEvap _ ref
COPref   8.945
WComp _ ref
La Tabla 1-6 presenta el orden de las corridas experimentales y los resultados para cada una de las
variables de respuestas. Siendo un diseño experimental 23 con 4 bloques; los tamaños de bloque
2k-p son 2; es decir, que se ejecutaron 2 corridas por cada bloque. A su vez, como se explica en
Montgomery [51], los efectos confundidos con los bloques son las interacciones de segundo orden
(AB, AC y BC) .
Tabla 1-6 Orden de experimentación y resultados obtenidos de las variables de respuestas (factores codificados)
Orden de experimentación y variables de respuestas

Velocidad
Temperatura Humedad QEvap WComp
BLOQUE ventilador COP
°C g/kg (kW) (kW)
Nivel
1 1 -1 1 3,619 0,46621 7,763
1 -1 1 2 3,561 0,47597 7,481
2 -1 -1 2 3,579 0,51839 6,904
2 1 1 1 3,613 0,46564 7,76
3 -1 1 1 3,726 0,46081 8,086
3 1 -1 2 3,655 0,47469 7,7
4 -1 -1 1 3,805 0,61982 6,14
4 1 1 2 3,372 0,6719 5,018
1.5 Metodología de análisis de resultados
Para analizar los resultados del diseño de experimento según variable de respuesta, se realiza el
siguiente procedimiento:
1) Siendo un diseño factorial sin réplica es necesario examinar la gráfica mitad normal de los
efectos estandarizados para evaluar los efectos que son significantes[51]
2) Una vez excluido los efectos que no son significantes, se estudia la tabla de Análisis de
Varianza (ANOVA); donde se verifica si los efectos que quedaron en el modelo tengan un
valor p < 0.05, para establecer que los factores (y/o interacciones) son significativas para
la variable de respuesta.
3) Se verifica los estadísticos R2 y R2 Ajustada de la ANOVA
4) Se verifica la validación de los supuestos del experimento a través de las gráficas de
probabilidad normal (normalidad), residuos vs número de corrida (independencia) y
residuos vs predicho (homocedasticidad)
5) Finalmente se analizan las gráficas de efectos principales, gráfica de interacción y
finalmente superficie de respuesta.
En esta sección se presentarán únicamente los pasos 2 y 5, el resto de análisis se presentan para
cada una de las variables de respuestas en el Anexo II “Resultados y validación del diseño de
experimento”
1.6 Análisis de resultados para el calor absorbido en el evaporador
La Tabla 1-7 presenta la tabla ANOVA1, correspondiente a la significancia entre los factores
principales y sus interacciones respecto al calor cedido por el aire en el evaporador, la obtención
del ANOVA se observa en el Anexo II
Tabla 1-7 Resultado ANOVA para el calor cedido del aire a través del evaporador en función de los efectos
principales y sus interacciones. AB: Temperatura-Humedad, AC: Temperatura-Nivel velocidad, BC: Humedad-
Nivel velocidad.
Análisis de Varianza para Q_Evap

Fuente Suma de Cuadrados Gl Cuadrado Medio Razón-F Valor-P
A:Temperatura_aire_in 0,021218 1 0,021218 1,49 0,4374
B:Humedad_aire_in 0,0186245 1 0,0186245 1,30 0,4579
C:Nivel_velocidad 0,044402 1 0,044402 3,11 0,3284
AB+bloque 0,010404 1 0,010404 0,73 0,5502
AC+bloque 0,0043245 1 0,0043245 0,30 0,6797
bloques 0,005832 1 0,005832 0,41 0,6380
Error total 0,0142805 1 0,0142805
Total (corr.) 0,11329 7
Como se puede observar, ninguno de los factores principales, ni sus interacciones son
significativas para esta variable. Al realizar un análisis detallado del sistema experimental, se
observó que la temperatura del aire a la salida del evaporador es regulada por el lazo de control
del inverter; regulando el flujo másico del refrigerante bajo cualquier condición de entrada del aire
para desarrollar un efecto refrigerante adecuado. Lo anterior es válido siempre y cuando la
temperatura del aire a la salida del evaporador sea mayor a la temperatura de control; en caso que
la temperatura disminuya por debajo de este nivel, un segundo lazo de control actúa en el sistema.
En este lazo, el inverter envía la señal al compresor para disminuir el flujo másico y poder aumentar
la temperatura de salida del aire; este lazo de control evita el congelamiento de las tuberías del
evaporador (que a su vez podría generar problemas en el compresor). Al alcanzar la temperatura
1
Análisis de Varianza, una prueba de hipótesis para probar la igualdad de medias de los grupos a través de la
comparación de varianzas estadísticas.
de salida deseada, el lazo es desactivado y vuelve a actuar el primer lazo de control. En caso que
el sistema disminuya la temperatura una vez más, vuelve a actuar este segundo lazo de control.
Esto fue un impedimento en la realización de puntos extremos más alejados durante la
experimentación como fue explicado previamente y como se observa en la Figura 1-8.
Figura 1-8 Comportamiento flujo refrigerante con los lazos de control del inverter y la temperatura del aire a la salida
en condiciones extremas de operación. Fuente: elaboración propia
1.7 Análisis de resultados para el trabajo eléctrico del compresor

Se analizan los resultados del diseño de experimentos para el trabajo del compresor; en el anexo
II, sección b se presentan los resultados para validar los efectos significativos para esta variable de
respuesta. La Tabla 1-8 presenta el ANOVA donde se indica que el factor bloque y las
interacciones Temperatura-Humedad del aire, Temperatura del aire-Nivel de velocidad del
ventilador y Humedad del aire-Nivel velocidad del aire describen en gran proporción la
variabilidad del compresor
Tabla 1-8 Resultado ANOVA para el consumo eléctrico del compresor en función de las interacciones. AB:
Temperatura-Humedad, AC: Temperatura-Nivel velocidad, BC: Humedad-Nivel velocidad.
Análisis de Varianza para WComp

AB+bloque 0,0198075 1 0,0198075 27,40 0,0064
AC+bloque 0,0113259 1 0,0113259 15,67 0,0167
BC+bloque 0,0123536 1 0,0123536 17,09 0,0144
Error total 0,0028914 4 0,00072285
Total (corr.) 0,0463783 7
Como se presentó en la Tabla 1-8, las interacciones (sin los efectos principales) afectan
estadísticamente el consumo eléctrico del compresor, en la Figura 1-9 se observa que la interacción
Temperatura-Humedad afecta significativamente, y de forma positiva, en el consumo del
compresor; las otras interacciones también afectan de forma positiva este consumo, pero a menor
medida. La gráfica de interacciones (Figura 1-10) nos permite analizar más a fondo como es
afectado el consumo con estas variables. El compresor alcanza un mayor consumo eléctrico
cuando las interacciones AB, AC y BC se encuentran en nivel bajo, o en nivel alto. Es decir, que
el consumo energético de la unidad de aire acondicionado aumenta cuando este recibe una
corriente de aire con un nivel alto de temperatura y de humedad o si recibe una corriente de aire
con una temperatura y humedad baja. Otros casos donde se pueden presentar un mayor consumo,
ocurre si el equipo opera con un nivel alto de velocidad del ventilador Y el aire a la entrada que
ingresa al evaporador, se encuentre a alta temperatura o alta humedad. Por último, se presenta un
mayor consumo si el equipo opera con un nivel bajo de velocidad en el ventilador y el aire ingresa
con una temperatura y/o humedad baja.
Diagrama de Pareto Estandarizada para W_comp
+
-
AB+bloque
BC+bloque
AC+bloque
0 1 2 3 4 5 6
Efecto estandarizado
Figura 1-9 Resultado ANOVA Diagrama de Pareto para el ANOVA. Fuente: elaboración propia
Gráfica de Interacción para W_comp
0,58
- +
0,56 - + - +
0,54
W_comp
0,52
0,5
0,48 + - + -
+ -
0,46
-1,0 1,0 -1,0 1,0 -1,0 1,0
AB AC BC
Figura 1-10 Gráfico de interacciones para el consumo del compresor. Fuente: elaboración propia
Estos resultados también son respaldados por las gráficas de superficies de respuestas que se
presentan en la Figura 1-11 donde se observa una fuerte interacción entre la temperatura del aire a
la entrada del evaporador, cuando el ventilador tiene un nivel bajo de velocidad (figura 16
izquierda); pero a su vez se observa una interacción fuerte de la humedad cuando el ventilador
tiene un nivel alto de velocidad (figura 16 derecha). En el caso del nivel bajo de velocidad en el
ventilador es entendible, ya que una menor velocidad indica una mayor permanencia del aire en el
evaporador, permitiendo un mejor intercambio de calor y una transferencia de masa más estable;
si las condiciones psicrométricas del aire son muy altas, requerirá un mayor consumo energético
el proceso de transferencia de calor/masa.
Por otra parte, a una mayor velocidad del ventilador, la humedad absoluta aire afecta más el
proceso de transferencia de calor; a mayor velocidad, mayor número de Reynolds, esto disminuye
la capa límite de agua condensada que se forma en el evaporador, causando un goteo excesivo en
el evaporador. Este goteo afecta la transferencia de calor, causando una mayor resistencia [8][18];
por consiguiente, una mayor necesidad de flujo de refrigerante y un mayor consumo energético.
Superficie de Respuesta Estimada Superficie de Respuesta Estimada
Rango Velocidad del aire: 1.5 -2.0 m/s Re_Max = 1800 Rango Velocidad del aire: 2.5 - 3 m/s
Re_Max = 1200 Nivel_velocidad=1,0
Nivel_velocidad=-1,0
0,7 0,7
0,66 0,66
0,62
W_comp
0,62
W_comp
0,58 0,58
0,54 0,54
0,5 0,5
-1 1
-0,6 0,46 0,6
0,46 -0,2 0,2
0,2 -1 -0,2
1 0,6 0,2 0,6 -0,6 -0,2 -0,6
-0,2 -0,6 -1 1 Humedad_aire_in 0,2 0,6 -1 Humedad_aire_in
1
Temperatura_aire_in Temperatura_aire_in
Figura 1-11 Superficie de respuestas para el consumo eléctrico del compresor en kW. Ejes: Temperatura codificada
del aire a la entrada del evaporador, Humedad absoluta codificada del aire a la entrada del evaporador. Izq. Nivel bajo
de velocidad, el Rango de velocidad para este nivel y máximo Reynolds alcanzado. Der. Nivel alto de velocidad, rango
de la velocidad para este nivel y máximo Reynolds alcanzado. Fuente: elaboración propia
Todo lo analizado anteriormente indica que, para un mejor consumo energético del compresor, el
cual es gobernado actualmente por el inverter, es necesario establecer estrategias de control donde
se incluya la medición de la humedad del aire y la manipulación de la velocidad del ventilador. Ya
que, los sistemas actuales orientadas al ahorro energético no están tomando en cuenta la
deshumidificación y cuando estos sistemas de acondicionamiento de aire tipo residencial se
enfrentan a estas condiciones climáticas el consumo energético se verá afectado.
1.8 Análisis de resultados para el coeficiente de desempeño del equipo

La Tabla 1-9 presenta el ANOVA para el COP; la validación de este diseño se observa en la sección
c del anexo II; se indica que el factor bloque, las interacciones, el factor principal “Nivel de
velocidad” son estadísticamente significativas y afectan significativamente el COP del equipo.
Tabla 1-9 Resultado ANOVA para el coeficiente de desempeño en función del nivel de velocidad: C, y las
interacciones. AB: Temperatura-Humedad, AC: Temperatura-Nivel velocidad, BC: Humedad-Nivel velocidad.
Análisis de Varianza para COP

C:Nivel_velocidad 0,875164 1 0,875164 11,18 0,0443
AB+bloque 3,39041 1 3,39041 43,30 0,0071
AC+bloque 1,09816 1 1,09816 14,02 0,0332
BC+bloque 2,04829 1 2,04829 26,16 0,0145
Error total 0,234905 3 0,0783018
Total (corr.) 7,64693 7
Fuente: Elaboración propia
Los resultados presentados por la Tabla 1-9 y la Figura 1-12, indica que el efecto más significativo
es la interacción Aire-Humedad; seguido de la interacción humedad-nivel velocidad; temperatura-
nivel velocidad y por último el efecto principal nivel de velocidad. Todos estos efectos ejercen una
influencia inversamente proporcional al COP, es decir, al aumentar el nivel de un factor
significativo (o todos los factores) el desempeño del equipo disminuye significativamente.
Diagrama de Pareto Estandarizada para COP
+
AB+bloque -
BC+bloque
AC+bloque
C:Nivel_velocidad
0 2 4 6 8
Efecto estandarizado
Figura 1-12 Diagrama de Pareto para el ANOVA del COP. Fuente: elaboración propia
La Figura 1-13 izquierda, presenta la gráfica de los efectos principales, se puede identificar que el
desempeño de la unidad se reduce en 9.5% únicamente por aumentar el nivel de velocidad del
ventilador. En la gráfica de interacción (Figura 1-13 derecha) se observa que la unidad obtiene los
mejores niveles de COP cuando el nivel de velocidad del ventilador es bajo, y los niveles de
humedad son altos (interacción BC). Otro nivel del alto desempeño ocurre cuando la velocidad del
ventilador es baja, y los niveles de temperatura sea alto (interacción AC). Finalmente, para
interacción AB, se presentan un mejor aprovechamiento del consumo energético de la unidad en
dos casos 1) nivel de temperatura bajo, nivel de humedad alto, y 2) nivel de temperatura alto, nivel
de humedad bajo.
Gráfica de Efectos Principales para COP Gráfica de Interacción para COP
7,5 8
7,43725 -
+ - -
7,7
7,3
7,4
+
COP
COP
7,1 7,1 +
-
-
6,8
6,9
6,5 - +
6,77575 +
6,7 +
6,2
1 2 -1,0 1,0 -1,0 1,0 -1,0 1,0
Nivel_velocidad AB AC BC
Figura 1-13 Gráfico de efectos principales (izquierda) e interacciones (derecha) para el Coeficiente de Desempeño.
Cuando se analiza la superficie de respuesta estimada (Figura 1-14) se aprecia que el nivel de
velocidad afecta significativamente en el coeficiente de desempeño; cuando el ventilador del
evaporador trabaja en un nivel bajo de velocidad, se obtiene COP en el rango de 5.8-8.2; mientras
que a un nivel alto de velocidad, el COP se mantiene en el rango de 5.4 y 7.6 Puede entenderse a
partir de la superficie de respuesta para un nivel bajo de velocidad (Figura 1-14 izquierda), ), que
la variación del COP es mínima cuando la temperatura del aire se encuentra en un nivel alto, sin
importar los niveles de humedad. A su vez, se observa a niveles altos de humedad que, sin importar
la temperatura, el COP no tiene una variación significativa. Sin embargo, al disminuir algunos de
estos factores (o en simultáneo), el desempeño energético comienza a variar significativamente
hasta alcanzar el menor valor; cuando la humedad y temperatura hayan disminuido.
Para la superficie de respuesta para un nivel alto de velocidad (Figura 1-14 derecha) además de
observar una disminución significativa del coeficiente de desempeño, se observa un
comportamiento distinto al expresado en el nivel de velocidad bajo. Se observa niveles estables de
COP cuando la temperatura y la humedad se encuentran en niveles bajos, pero una disminución
significativa del COP cuando uno (o los dos) factores incrementan su valor.
Rango de Velocidad del aire: 2.5 - 3 m/s
Superficie de Respuesta Estimada Superficie de Respuesta Estimada

Re_Max = 1200 Rango de Velocidad del aire : 1.5 - 2 m/s Re_Max = 1800 Rango de Velocidad del aire: 2.5 - 3 m/s
Nivel_velocidad=-1,0 Nivel_velocidad=1,0
9 9
8 8
COP
COP
7 7
6 6
1
0,6 5 -1
5 0,2 -0,6
-0,2 -0,2
1 0,6 0,2 -0,6 -1 -0,6 0,2
-0,2 -0,6 -1 -1 Temperatura_aire_in -0,2 0,2 0,6
0,6 1 1 Temperatura_aire_in
Humedad_aire_in Humedad_aire_in
Figura 1-14 Superficie de respuestas para el coeficiente de desempeño del equipo. Ejes: Temperatura codificada del
aire a la entrada del evaporador, Humedad absoluta codificada del aire a la entrada del evaporador Izq. Nivel bajo de
velocidad, rango de velocidad para este nivel y máximo Reynolds alcanzado Der. Nivel alto de velocidad, rango de
velocidad para este nivel y máximo Reynolds alcanzado. Fuente: elaboración propia
Este comportamiento está asociado con lo explicado en la sección de consumo eléctrico. A niveles
bajo de velocidad del ventilador, el proceso de intercambio de calor es más estable en comparación
a niveles de alta velocidad; cuando el ventilador del evaporador trabaja a niveles bajo de velocidad,
y la temperatura del aire es alta, el control del inverter permite que el consumo energético sea casi
igual. En caso de una alta humedad, el fenómeno de condensación del agua en el evaporador es un
fenómeno de condensación en presencia de un gas no condensable; este fenómeno a un nivel de
baja velocidad del aire genera una película de agua alrededor de las placas y tubos del evaporador
causando una mejora en la transferencia de calor en el evaporador [8], [55], mejorar la
transferencia de calor permite una disminución de refrigerante y mejora el coeficiente de
desempeño.
Al disminuir la temperatura del aire a la entrada, el sistema de control puede tener dificultades para
regular la temperatura a la salida, e incluso esta pueda salir por debajo (o cerca) de la temperatura
de control, causando el descontrol del equipo explicado en el análisis del consumo del compresor.
A su vez, un aire menos húmedo no permite formar la película de agua de manera uniforme, esto
puede disminuir el coeficiente de transferencia de calor y en consecuencia una disminución del
COP.
El caso de la temperatura y humedad cuando el equipo trabaja con un nivel alto de velocidad del
ventilador; podemos explicarlo a partir de la ecuación (2); el equipo tiene una limitación física: la
capacidad de enfriamiento a la cual fue diseñada; si el equipo aumenta el nivel de velocidad del
ventilador, el aire, seguirá teniendo la misma capacidad de enfriamiento, es debido a esto, que
deberá compensarse con la disminución de la variación de entalpía y esta puede ser compensada
ya sea con un nivel bajo de temperatura, o un nivel bajo de humedad absoluta. Esto es válido,
siempre que el aire a la salida del evaporador no alcance (o se encuentre por debajo) de la
temperatura de control. Además, un bajo nivel de temperatura y humedad un nivel de velocidad
alta del ventilador en el evaporador, genera una menor capacidad de enfriamiento en el equipo,
por consiguiente, un menor flujo másico de refrigerante y en conclusión un mayor COP.
Aire a un nivel alto de temperatura no sobrepasará la capacidad de enfriamiento, simplemente
trabajará hasta una temperatura de aire a la salida del evaporador mucho mayor de la esperada; por
lo tanto, si el equipo se encuentra a un nivel alto de velocidad del ventilador y recibe aire a alta
temperatura y/o alta humedad, requerirá de más refrigerante para realizar su trabajo, en conclusión,
más consumo energético y un menor COP. Como se explicó anteriormente, el aumento de la
velocidad del ventilador, aumenta el flujo del aire, causando una disminución de la capa límite del
agua condensada, formando goteo en el aire que afecta el proceso de transferencia de calor
significativamente.
1.9 Conclusiones del capítulo

Para este capítulo tenemos las siguientes conclusiones:
 La temperatura, humedad y velocidad del ventilador son factores significativos en el
desempeño energético de un sistema de aire acondicionado tipo residencial que trabaja bajo
condición de clima tropical. Su significancia proviene de las interacciones, lo cual es
completamente lógico al ser un proceso de enfriamiento con fines de deshumidificación de un
vapor condensable en presencia de un gas no condensable.
 Si el espacio a acondicionar se encuentre en una condición de alta temperatura y humedad, es
recomendable que el equipo trabaje inicialmente a baja velocidad del ventilador, garantizando
un mayor tiempo de permanencia del aire en el evaporador, para poder deshumidificar más el
aire y tener un mejor desempeño energético.
 El mejor desempeño energético en este experimento se encontró cuando el nivel de velocidad
del ventilador es bajo (nivel bajo), la temperatura es baja (23°C) y la humedad absoluta es alta
(14 g/kg). Estas condiciones de temperatura y humedad pueden presentarse típicamente en
climas tropicales.
 Un ajuste del nivel de velocidad del ventilador permite el ahorro de hasta el 10% en consumo
energético del sistema de acondicionamiento de aire.
 Se observa que los factores significativos obtenidos para el COP difieren de los factores
significativos obtenidos para el Wcomp, esto ocurre porque es posible que el Qevap sea
descrito con alguno de los 3 factores y/o sus interacciones presentadas, pero con un modelo
diferente al analizado en la sección 1.5. Se recomienda la ejecución de un diseño experimental
diferente a fin de identificar el correcto modelo que describa el Qevap y poder encontrar una
relación con el modelo que se describe en el COP.
2 Análisis del proceso de transferencia de calor entre el aire de clima
tropical y el refrigerante R410a en el interior de un evaporador
El proceso de transferencia de calor del aire de clima tropical es un fenómeno muy particular.
Siendo un aire con alta humedad, el aire no puede analizarse simplemente como aire seco, o
considerarlo despreciable. Cuando este aire ingresa al evaporador El aire disminuirá su
temperatura hasta que el aire húmedo alcance la temperatura de rocío, cuando esto ocurra,
además de disminuir su temperatura a menor ritmo en condiciones de saturación, parte del
vapor de agua que contiene es condensado, separándose de la corriente de aire y adhiriéndose
a la tubería del evaporador. El aire húmedo saldrá del evaporador a menor temperatura y
humedad gracias al proceso de intercambio de calor con el refrigerante. El fenómeno de
deshumidificación del aire húmedo en el evaporador se conoce como condensación en
presencia de un gas no condensable que para este caso es aire[55]. Este fenómeno es
completamente diferente al de condensación de una sustancia pura, ya que existe una interface
liquido gas que generan un efecto completamente distinto al proceso de condensación pura. En
la Figura 2-1 se observa una corriente de aire húmedo, conformado por una mezcla de aire
seco y vapor de agua, que fluye a través de una superficie que se encuentra a una temperatura
TW; entonces esta temperatura es lo suficientemente baja para alcanzar la temperatura de
saturación del vapor de agua a la presión parcial del aire húmedo (TA)El fenómeno de
condensación ocurrirá, el agua condensada se adhiere a la superficie formando una capa límite
de agua condensada. Al coexistir el agua en estado líquido y vapor de agua, ocurrirá
transferencia de masa por difusión entre el aire húmedo y el agua condensada, en consecuencia,
la humedad absoluta (proporción agua-aire) en la superficie del agua (WA) es mayor con
respecto a la humedad absoluta en la corriente de aire húmedo (WO). Además de lo anterior, el
aire seco formará una capa límite de velocidad con velocidad cero en la superficie del
condensado (U=0) y aumentará gradualmente hasta alcanzar la velocidad de superficie libre
(UO); por último, la capa de condesado genera a su vez una resistencia térmica, causando que
una diferencia entre la temperatura de la pared (Tw) la temperatura del agua condensada (Ta)
y la corriente de aire húmedo (To). En resumen, en climas tropicales, durante el enfriamiento
del aire en el evaporador, se formarán tres capas límites en el proceso (velocidad, temperatura,
humedad). Teniendo en cuenta lo anterior, podemos asegurar que las condiciones
psicrométricas del aire húmedo, y su velocidad afectan el fenómeno de transferencia de calor.
Figura 2-1 Análisis de condensación en presencia de un gas no condensable en la evaporador. Fuente: [19]
Además de la transferencia de calor, la deshumidificación también se verá afectada por el aire de

la región tropical. Antes que nada, hay que indicar que existe semejanza entre la transferencia de
calor y transferencia de masa; además, el mecanismo de transferencia de masa puede ser estimado
de igual forma como el fenómeno para la transferencia de calor [15], como se demuestran en la
analogía del factor J de Chiton Y Colburn; si el coeficiente de transferencia de calor es alto, sería
correcto afirmar que existe también un alto coeficiente de transferencia de masa. Sin embargo, hay
que tener en cuenta que el refrigerante que fluye dentro del evaporador tiene una limitación física
respecto a la capacidad de absorber calor, y un aumento en el coeficiente de transferencia de calor
no influye en la tasa de calor absorbido,, influye en el área superficial de transferencia de calor, es
decir al aumentar el coeficiente de transferencia de calor, el área superficial necesaria para realizar
la transferencia de calor será menor porque el refrigerante alcanzará más rápidamente la fase de
gas sobrecalentado. Al llevar este análisis a la transferencia de masa, el aumento del coeficiente
de transferencia de masa no afecta la diferencia de concentración, entre el aire húmedo y el agua
condensada; reducirá el área de condensación en el evaporador, por consiguiente, el aire realizará
la condensación del agua en un área menor y conlleva a un proceso de deshumidificación mucho
menor de lo esperado.
En esta sección, se investiga el proceso de transferencia de calor entre el aire húmedo y el
refrigerante R410A; de igual modo, se investigará la transferencia de masa entre el aire húmedo y
el agua condensada que se encuentra dentro de las paredes del evaporador. Esta investigación se
realiza con el fin de analizar el efecto de la transferencia de calor y masa del aire húmedo sobre el
proceso energético de un aire acondicionado tipo residencial. Para esto, se utilizará los resultados
obtenidos del sistema experimental de la sección anterior y se analizará los fenómeno de
transferencia de calor/masa en el evaporador; partiendo del método de reducción de Threlkeld [19]
2.1 Descripción matemática

En un intercambiador de calor, la tasa total de transferencia de calor que existe es un promedio
matemático entre el calor cedido por el aire con el calor absorbido por el refrigerante:
QAire  Qref (2-1)
QTotal 
2
Este calor promedio puede ser utilizado para evaluar el diseño de intercambiadores, por lo general
se aplica los conceptos de diferencia de temperatura media logarítmica (LMTD) o el método de
efectividad (ε-NTU) y a partir de los conceptos teóricos de transferencia de calor [18], se evalúa
el intercambiador. Esta metodología desarrollada por Threlkeld [19], es la más desarrollada a nivel
de investigación como se observa en [21], [22], [24], [56].
Ahora bien, uno de los inconvenientes del uso de esta metodología, es el desarrollo de una iteración
con base a un parámetro ficticio para “ajustar” las variables del proceso, necesario a su vez para
recalcular el coeficiente de transferencia de calor. Además, no toma en cuenta el proceso de
condensación que ocurre en el evaporador, aunque cabe resaltar que la metodología incluye
conceptos de aire húmedo, humedad absoluta y eficiencia húmeda de la aleta.
Para poder asociar el fenómeno de condensación en presencia de un gas no condensable, debemos
recordar que el proceso de transferencia de calor que ocurre en aire húmedo consiste no solamente
en el proceso de cambio de temperatura (calor sensible), también se debe incluir el proceso de
cambio de fase del agua condensada (calor latente); por lo tanto, la energía por unidad de área que
existe en el proceso de enfriamiento del aire es:
qAire  qSensible  qLatente (2-2)

Donde qSensible es el flujo de calor por unidad de área gracias a la diferencia de temperatura entre
el aire húmedo y la interface; qLatente es el flujo de calor por unidad de área generado por la
condensación del agua las ecuaciones (2-3)y (2-4) presentan como obtener estos parámetros:
qSensible  CpaireTa  Tpared  (2-3)
qLatente  m
 aguahfg (2-4)
Donde ṁaire es el flujo másico del aire húmedo CPaire es el calor específico en el aire húmedo, T A
y T pared son temperaturas promedio del aire en el evaporador y de la superficie de la pared
respectivamente. ṁagua es el flujo másico de agua condesada en la superficie y hfg es la entalpía de
evaporización del agua bajo la presión parcial correspondiente. El calor latente también puede
expresarse como la diferencia de humedad absoluta del aire a la temperatura promedio, con
respecto a la humedad de saturación del aire a la temperatura de la pared:
 airehfg wa  wpared 
qLatente  m (2-5)
qAire  CPaireTa  Tpared   hfg wa  wpared  (2-6)

Transformando la ecuación, podemos analizar la transferencia de calor en el aire húmedo con
fenómeno de transferencia de calor/masa:
CPaire   HT 
 HT Ta  Tpared   hfg wa  wpared 
qAire  (2-7)
 HT  CPaire 
El término que acompaña la diferencia de humedad se puede expresar también con el número de
Lewis:
2  HT 2  (2-8)
Le 3    MT Le 3  HT
CPaire MT CPaire
Por lo tanto,
CPaire 
 HT Ta  Tpared   MT Le 3 hfg wa  wpared 
2
qAire  (2-9)
 HT  
El término αMT es el coeficiente de transferencia de masa. El número de Lewis, por lo general es
tomado como una unidad, para simplificar el procedimiento [18], [19], [21], [23], sin embargo,
teniendo en cuenta que este número es importante en el fenómeno de condensación en presencia
de un gas no condensable, se calculará para cada proceso. El número de Lewis se calcula a partir
de la relación entre la difusión térmica y másica; también se puede expresar en términos del número
de Schmidt (Sc) y el número de Prandtl (Pr):
DifusiónTermica Sc (2-10)
Le  
DifusiónMásica Pr
El número de Schmidt es calculado utilizando la ecuación
 (2-11)
Sc 
DV
Donde μ, ρ son la viscosidad dinámica y la densidad del aire respectivamente. La variable D, es
la difusividad másica del vapor de agua en el aire húmedo y es calculada a partir de esta ecuación
0.926  T 2.5  (2-12)

DV   
p  T  245 
Las unidades para la presión p y la temperatura T para calcular la difuividad másica DV son kPa y
K respectivamente. Además, los parámetros Pr, μ, ρ son calculados a partir del software
REFPROP en función de p y T, esta última es la temperatura promedio del aire en el evaporador.
A continuación, recordemos que el calor cedido por el aire, consiste en la diferencia de entalpía
entre la superficie y el aire, considerando que el aire cambia constantemente su entalpía, no es
posible realizar únicamente una diferencia aritmética, McQuiston y Threlkeld [18], [19] proponen
la utilización de una diferencia de entalpía media logarítmica (LMED) para garantizar que la
transferencia de calor sea lo más precisa posible:
qAire  AAG hLM  AAG

h  hpared,1   haire, 2  hpared, 2 
aire,1
 haire,1  hpared,1  (2-13)


ln  

 haire, 2  hpared, 2 
Donde αG es el coeficiente global de transferencia de calor y masa del aire húmedo. Lo desarrollado
por Threlkeld y McQuiston [18], [19] difiere de este coeficiente por lo mencionado del número de
Lewis. Finalmente, al sustituir, qaire, en la ecuación (2-9) y considerando el trabajo presentado por
Xia et al [57] se tiene que el coeficiente global de transferencia de calor y masa del aire húmedo
por el área del evaporador por el lado del aire (AA) es:
1
AAG 
  
  
  
 b C paire  1  (2-14)
     
 Ai ref 0 AA HT 
2
Le 3  1 
 1 q 2

  sensible
 Le 3 
  qlatente 
Donde b es un coeficiente lineal, producto de una relación entre la entalpía del aire saturado y la
temperatura del aire saturado [19]. Ai el área del evaporador por el lado del refrigerante, αref el
coeficiente de transferencia de calor del lado del refrigerante, η0 la eficiencia total del
intercambiador (en función de la eficiencia de la aleta)
El coeficiente de transferencia de calor del refrigerante es calculado usando la ecuación
i   FZ S   f F (2-15)
Siendo αFZ el coeficiente de ebullición nucleada obtenida por la correlación de Forster y Zuber[29],
[31], αf el coeficiente de transferencia de calor para el refrigerante en fase líquida, y S y F factores
de correlación en función de la calidad del refrigerante a la entrada del evaporador [31]. Se omite
la transferencia de calor del refrigerante en fase gaseosa porque se estima que la transferencia de
calor cuando el refrigerante se encuentra en este estado es mínima en comparación cuando está en
estado bifásico.
2.2 Cálculos de coeficiente de transferencia de calor y masa en un evaporador de

tubos aleteados
A partir de las ecuaciones (2-1), (2-13), (2-14), (2-15), de la metodología LMTD (adaptada para
LMED), de la metodología para el cálculo de eficiencia para climas húmedos se analiza los
resultados obtenidos en el diseño experimental para calcular el coeficiente global de transferencia
de calor y masa del aire en el evaporador. Para realizar esto, un proceso iterativo se debe realizar
de la siguiente forma:
a) Se evalúan las condiciones geométricas termodinámicas y psicrométricas del evaporador,
aire y refrigerante, para obtener áreas, αi, Cp, Le, qsensible y qlatente
b) Se asume un coeficiente αHT
c) Se calcula la eficiencia de la aleta y seguido a esto la eficiencia total del evaporador η0,
para este estudio, el cálculo se realizó a partir de la metodología presentada por Wang [21]
d) Se calcula el coeficiente global de transferencia de calor y masa del aire húmedo ecuación
(2-14)
e) Se calcula el calor cedido por el aire con la ecuación (2-13)
f) Se compara el calor obtenido con el calor calculado con la ecuación (2-1), en caso que el
resultado no sea satisfactorio, repetir desde el paso b)
En la Tabla 2-1 se observa inicialmente por cada corrida experimental, el número de Reynolds
calculado. Luego, el coeficiente de transferencia de calor (αHT) del aire utilizando la metodología
presentada por Wang [21], [23]. Seguido de esto, se observa los resultados obtenidos del
coeficiente de transferencia de calor (αHT), coeficiente de transferencia de masa (αMT) y el
coeficiente global de transferencia de calor y masa (αG). Estas tres últimas variables fueron
calculadas con la metodología explicada previamente.
Tabla 2-1 Resultado del coeficiente global de transferencia de calor y masa, coeficiente de transferencia de calor y
coeficiente de transferencia de masa para cada corrida experimental realizada.
Coeficiente de transferencia de calor

Coeficiente Coeficiente Coeficiente
Coeficiente
αHT αMT αG
αHT 2 2
[W/(m K)] [kg/(m s)] [kg/ (m2s)]
Corrida Re [W/(m2 K)]
Fuente: Fuente: Fuente:
Fuente: [21],
Elaboración Elaboración Elaboración
[23]
propia propia propia
1 985,034 29,05 43,93 0,0494 0,1223
2 1151,353 35,51 24,190 0,0269 0,1161
3 1677,192 25,38 24,360 0,032 0,1247
4 1105,271 27,21 26,11 0,027 0,1045
5 1243,820 30,14 26,300 0,029 0,1123
6 1819,337 20,14 19,330 0,028 0,1127
7 1018,052 48,26 48,610 0,0545 0,1480
8 1151,511 22,14 28,020 0,0313 0,1011
Al comparar los resultados del coeficiente de transferencia de calor de la metodología propuesta

por Wang y la presentada en este trabajo, de la Tabla 2-1 se pueden observar diferencia entre los
coeficientes de transferencia de calor; esto se debe principalmente al incluir el número de Lewis
dentro de la metodología, en lugar de asumir Le=1. En la Tabla 2-2 se puede observar que la
variación del número de Lewis por corrida es menor que 1, por esta razón, se afecta el coeficiente
de transferencia de calor (por lo tanto, el coeficiente de transferencia de masa) siendo menor en la
mayoría de las corridas experimentales.
Tabla 2-2 Resultado, Difusividad másica entre vapor de agua-aire y el número de Lewis por cada corrida
experimental
Difusividad másica vapor agua-aire y Número de Lewis

Temperatura DV
Corrida Le
Aire entrada (°C) (mm2/s)
1 28,176 25,041 0,8301
2 23,446 24,818 0,842
3 22,822 24,678 0,8412
4 27,385 25,217 0,8371
5 23,717 24,753 0,8396
6 28,116 25,397 0,8387
7 23,549 24,407 0,8319
8 27,641 25,280 0,8378

La Figura 2-2 presenta la comparación entre resultados; se puede identificar que no existe una
variación en la tendencia de los coeficientes.
Figura 2-2 Comparación del coeficiente de transferencia de calor obtenido por el modelo desarrollado y el coeficiente
presentado en la tabla 13. Fuente: elaboración propia
Desde un punto actual de diseño, no parece significativo la inclusión del número de Lewis para el
cálculo del coeficiente de transferencia de calor, sin embargo, al analizar los resultados de ambos
métodos en una carta psicrométrica, con base a las condiciones de entrada del aire húmedo al
evaporador, se puede observar una diferencia entre la mayoría de los resultados. En la Figura 2-3
se observa las condiciones psicrométrica del aire húmedo a la entrada del aire, los colores
representan el coeficiente de transferencia de calor calculado (azul= el menor coeficiente, rojo= el
mayor coeficiente), y el tamaño de los puntos representa la velocidad frontal del aire húmedo
(mayor tamaño, mayor velocidad).
Figura 2-3 Comparación de resultados del coeficiente de transferencia de calor en la carta psicrométrica con base a
las condiciones de entrada del aire. Izquierda: método tradicional asumiendo Le=1. Derecha: método presentado en
este documento. Fuente: elaboración propia.
2.3 Análisis de resultados para el coeficiente de transferencia de calor

Si analizamos ambos métodos se presenta coherencia entre los resultados y tendencias asociada al
fenómeno de transferencia de calor y masa:
1. Al aumentar la velocidad y/o temperatura del aire húmedo, disminuye el coeficiente de
transferencia de calor
2. Al aumentar la humedad mientras que el aire fluya a una velocidad entre 2.5-3.5 m/s
(establecida en este documento como rango alto) el coeficiente de transferencia de calor
aumenta
3. Si el aire a la entrada del evaporador aumenta su humedad, y además el aire fluya a una
velocidad entre 1,5 -2.2 m/s (rango bajo) la transferencia de calor disminuirá
La teoría de condensación en presencia de un gas no condensable, indica que la mayor
concentración de humedad se encuentra en las paredes, un aumento de la humedad en el aire puede
mejorar proceso de transferencia de calor, ya que aumenta la capa límite del agua y este a su vez
absorbe más calor, del aire y desde el refrigerante, siempre y cuando la velocidad del aire permita
la formación de la capa límite. Si la velocidad del aire aumenta, se producirá desprendimiento de
agua condensada en las paredes de la tubería, esto provoca una disminución de la capa límite del
agua condensado, mientras que el goteo generado por el agua condensada causaría una resistencia
extra en el aire, ocasionando una disminución de la transferencia de calor [8].
En comparación entre los métodos presentados, se observa un mejor comportamiento del modelo
desarrollado cuando los efectos de la humedad son más significativos; por lo presentado en el
capítulo anterior, los efectos de la humedad predominan en la variable de respuesta cuando la
temperatura del aire trabaja en altos niveles. Se observa que la transferencia de calor cuando trabaja
a un nivel de alta temperatura, alta velocidad y baja humedad, el coeficiente de transferencia de
calor es similar, sin embargo, a otros niveles de humedad y velocidad el coeficiente varía. De igual
forma al analizar el comportamiento de los modelos cuando los efectos de la temperatura
predominan, se observa un comportamiento muy similar, aunque existe una mejor sensibilidad del
proceso en el modelo actual, es decir, al no considerar variación del número de Lewis se presenta
una variación mayor del coeficiente de transferencia de calor en comparación de si considerarlo.
La Figura 2-4 presenta las condiciones psicrométricas del aire a la entrada del evaporador en cada
corrida experimental, en esta ocasión, en la figura del lado izquierdo, los colores representan el
coeficiente combinado de transferencia de calor y masa, mientras que la figura derecha los colores
indican el valor del coeficiente de transferencia de calor. Por último, los tamaños de los puntos
representan el consumo eléctrico del compresor (a mayor círculo, mayor consumo).
Analizando el comportamiento de la transferencia global, se observa que no existe un cambio
significativo, ni el consumo ni el valor del coeficiente, cuando las condiciones el aire trabaja a baja
temperatura - alta humedad; o cuando trabaja a alta temperatura – baja humedad. Las afectaciones
de este coeficiente corresponden básicamente en las otras dos condiciones operacionales (alta
humedad - alta temperatura o baja humedad - baja temperatura), ya que se puede observar que
alcanza su alto valor cuando las condiciones de la temperatura y humedad aire son bajas, sin
embargo, alcanza sus rangos más bajos cuando la temperatura y humedad del aire trabaja en altas
condiciones (28°C y 14 g/kg). Teniendo en cuenta el consumo eléctrico para estas condiciones, se
obtiene un alto consumo eléctrico cuando se tiene un coeficiente por debajo de 0,1 kg/ (m2 s) o un
coeficiente mayor a 0,15 kg/ (m2 s).
Figura 2-4 Izquierda: Coeficiente global de transferencia de calor y masa. Derecha: Coeficiente de transferencia de
calor. Ambas gráficas en función de las condiciones de entrada del aire y del consumo eléctrico del compresor. Fuente:
elaboración propia
Con base a la Figura 2-4, los puntos extremos del coeficiente global de transferencia de calor y
masa (menor y mayor denominación) se alcanzan cuando se trabaja a un nivel bajo de velocidad a
temperatura y humedad baja, y cuando el nivel de temperatura y humedad es alta el coeficiente
global alcanza un nivel máximo cuando el nivel de la velocidad del aire es alta, esto expresa la
relación entre la temperatura, humedad y velocidad del aire con respecto al consumo eléctrico del
compresor. Teniendo en cuenta lo anterior, se confirma lo expresado en el capítulo anterior: la
operación adecuada del sistema se obtiene cuando existe interacción de la humedad y temperatura.
Esto también se puede observar en la figura derecha de la Figura 2-4, el consumo es menor en los
niveles de baja temperatura – alta humedad; alta temperatura – baja humedad, sin embargo el
coeficiente es más estable cuando trabaja en el nivel de baja temperatura – alta humedad y esto se
refleja en un alto COP de la unidad del sistema.

Analizando para clima de zonas tropicales, si el aire acondicionado recibe aire a alta temperatura
y humedad, el coeficiente de transferencia de calor (y el global) sería muy bajo (20 W/m2 °C y 0.1
kg/m2s respectivamente), y básicamente el proceso de enfriamiento será dominado por los efectos
de la humedad. Se ha identificado que poder controlar la velocidad del ventilador en el evaporador,
permitirá controlar un poco más el proceso de transferencia de calor, asegurando por lo menos que
el equipo no sea castigado energéticamente por las condiciones operacionales. Esto se pudo
observar durante el capítulo 1; cuando el aire del espacio acondicionado se encuentre en zona de
sobreenfriamiento (alta humedad, baja temperatura) lo recomendable es mantener un bajo nivel de
velocidad para garantizar una operación óptima del equipo; sin embargo, cuando el aire del espacio
disminuya su humedad entre 10-11 g/kg humedad, es necesario aumentar la velocidad de giro del
ventilador con el fin de mantener el coeficiente de desempeño del equipo en niveles aceptables de
consumo energético.
Otras conclusiones que se obtienen de este capítulo se expresan a continuación:
 Se obtuvo un coeficiente global de transferencia de calor y masa que considera los efectos de
condensación del aire húmedo en el proceso.
 El coeficiente de transferencia de calor es afectado por la humedad del aire, este coeficiente
aumenta si la temperatura en el aire a la entrada del evaporador disminuye; pero la variación
de este coeficiente con respecto a la humedad depende de la velocidad del aire. A bajos niveles
de velocidad, un aumento de la humedad disminuye el coeficiente, pero a un alto nivel de
velocidad, el coeficiente aumenta si la humedad absoluta aumenta.
 Cuando los modelos de transferencia de calor asumen un número de Lewis = 1, se presenta
una variación mayor del coeficiente de transferencia de calor respecto a los resultados
encontrados en esta investigación, donde el número de Lewis utilizado es el real, el cual se
encuentra alrededor de 0.8.
 El aire de clima tropical que ingresa al evaporador del aire acondicionado a alta temperatura y
humedad, afectará el coeficiente global de transferencia de masa y calor, pues la variable tendrá
un valor muy bajo, requiriendo en el diseño del evaporador, que el área de transferencia de
calor sea alta, pero, básicamente el proceso de enfriamiento será dominado por los efectos de
la humedad.
 Una reducción de la velocidad del ventilador, permitirá controlar el proceso de transferencia
de calor, asegurando por lo menos que el equipo no consuma energía excesiva por las
condiciones operacionales.
 Si el aire que ingresa se encuentra en zona de sobreenfriamiento (alta humedad, baja
temperatura) es recomendable mantener un bajo nivel de velocidad para garantizar una
operación adecuada del equipo. Se debe tener cuidado al alcanzar niveles bajos de temperatura
y humedad, a fin de evitar consumo excesivo de energía eléctrica en el equipo.
3 Análisis estadísticos a partir del resultado en el sistema

experimental
En este capítulo, se desarrolla tres análisis estadísticos a partir de los resultados experimentales en
la influencia del aire húmedo sobre un sistema de acondicionamiento de aire DX tipo residencial.
Inicialmente se presenta un modelo de regresión para los coeficientes de transferencia de calor y
de masa; esta se desarrolla a partir de los resultados obtenidos en el capítulo anterior, la obtención
de este coeficiente permitirá realizar simulación del sistema para el desarrollo de futuras
investigaciones. Siguiendo a esto, se presenta el análisis de incertidumbre de los resultados, a fin
de evaluar la trazabilidad del error del diseño experimental. Por último, a partir de los resultados
experimentales se plantea una propuesta teórica sobre estudios de sistemas de acondicionamiento
de aire DX para clima tropical.
3.1 Modelo de regresión

Como se expresó al inicio de este trabajo, no existe un conceso general sobre la influencia de la
humedad relativa en el coeficiente de transferencia de calor, por lo tanto, a partir de los datos
experimentales del coeficiente de transferencia de calor y coeficiente de transferencia de masa, se
desarrolla un modelo de regresión para obtener un factor de colburn j, para el sistema de aire
acondicionado experimental con el propósito de caracterizar los coeficientes de transferencia de
calor y masa en aire de región húmeda a través del sistema experimental. Inicialmente los
coeficientes de trasferencia de calor y masa deberán expresarse en función de la analogía de
Colburn, la ecuación (3-1) y ecuación (3-2) presentan las relaciones de colburn con los coeficientes
de transferencia de calor y masa respectivamente.
 HT 2
(3-1)
j HT  Pr 3
C PaireGaire
 MT 2 (3-2)
j MT  Sc 3
U aire
Donde Gaire es el flujo másico por unidad de área del aire, Pr el número de Prandtl y Sc el número
de Schmidt, un número que relaciona el número de Prandtl (Pr) con el número de Lewis (Le). Por
lo general un modelo de factor j es de forma
j  c Re m (3-3)
Como se observan en [21], [23], [24], [53], [54], [56], y considerando que este estudio está
enmarcado bajo las condiciones operacionales del equipo como la capacidad de enfriamiento, la
potencia consumida, la geometría definida del evaporador y las velocidades que se pueden
desarrollar. Teniendo en cuenta que para obtener la estructura de la ecuación (3-3), se utilizaron
hasta 40 evaporadores de diferentes geometrías, la cantidad de información obtenida por ellos es
más precisa en comparación con las 8 corridas experimentales presentadas en este trabajo. Por
consiguiente, es muy probable que los modelos de regresión tengan una estructura diferente. Es
debido a esto, que se analizaron varios modelos de regresión con el fin de obtener un modelo
idóneo al comportamiento descrito. La Tabla 3-1 presenta 10 modelos de regresión que fueron
seleccionados para describir los factores de Colburn con el Número de Reynolds. Estos modelos
de regresión presentados fueron elegidos a partir de las recomendaciones del software
STATGRAPHICS, no obstante, el análisis de idoneidad, necesario para validar los modelos, fue
desarrollado a partir del paquete “Análisis de regresión” de Excel. Luego de validar, se
seleccionaron dos modelos de regresión, cuyos resultados se presentan en este capítulo. Los
análisis de idoneidad para los modelos de regresión no seleccionados se presentarán en el Anexo
1.
Tabla 3-1 Modelos de regresión recomendados por STATGRAPHICS (c=coeficiente; m=intercepto)
Forma general de los modelos de regresión
Modelo No Estructura
1
1 j
c Re
c
2 j
Re
2
3  c 
j  
 Re 
1
4 j
c Re
1
5 j
c Re 2
6 j  c Re m
7 j  ec Re
8 j  ec Re
j  c ln Re
2
9
Re 2
10 j

c Re 2  m 
Para la validación del modelo de regresión, se realiza un análisis de idoneidad, para esto, tres
pruebas estadísticas son realizadas: (1) test de ANOVA del modelo para evaluar si hay
significancia entre la variable de respuesta (factor de Colburn) y los parámetros de regresión (Re.),
(2) prueba de significancia de los coeficientes de los modelos para identificar su contribución
individual y (3) Análisis gráfico de los residuos con el propósito de garantizar que no existan
estructura u homocedasticidad en los residuos. Si algunas de estas pruebas fallan, el modelo no
cumple con la idoneidad y debe ser descartado
La Tabla 3-2 presenta los resultados de los factores j por cada experimento; los resultados de los
modelos de regresión y su prueba de idoneidad para el factor jHT y jMT, se presentan en la Tabla
3-3.y Tabla 3-4 respectivamente. Se observa, en ambas tablas que todos los modelos cumplen con
las 2 primeras pruebas estadísticas, sin embargo, solo 2 cumplen con el análisis de idoneidad en
su totalidad, (modelo 2 y modelo 10). Después de un análisis detallado de la gráfica de residuos,
se selecciona el modelo 10, como el aquel que describe idóneamente el Re con el factor jH y con
el factor jM ya que tiene mayor dispersión de sus residuos.
Tabla 3-2 Resultados del cálculo de factor JH y JM en las corridas experimentales en orden ascendente según Re.
Factor j en función de Re
Corrida Re_DC jHT_factor jMT_factor
1 985,03 0,0094 0,0115
7 1018,05 0,0102 0,0126
4 1105,27 0,0050 0,0060
2 1151,35 0,0044 0,0054
8 1151,51 0,0051 0,0061
5 1243,82 0,0045 0,0054
3 1677,19 0,0032 0,0032
6 1819,34 0,0023 0,0023
Tabla 3-3 Análisis de idoneidad de los modelos de regresión para jHT
Análisis de idoneidad
Prueba 2:
Prueba 1:
Modelo Significancia Prueba 3:
Estructura ANOVA
No Valor-p Gráfica residuos
Valor-p
m c
1 No cumple –
1 j 2,4837E-05 #N/A 8,02949E-06
Ajuste lineal mal calculado
c Re
c Cumple,
2 j 1,4684E-04 #N/A 6,26093E-05
Pero no escogido
Re
2
 c  No cumple - Ajuste lineal no
3 j   3,0111E-07 #N/A 4,75348E-08
adecuado
 Re 
1 No cumple - Ajuste lineal mal
4 j 2,2136E-04 #N/A 0,000100433
c Re calculado
1 No cumple - Ajuste no lineal con
5 j 1,8945E-06 #N/A 4,04444E-07
heterocedasticidad
c Re 2
No cumple - Ajuste lineal no
6 j  c Re m 1,9129E-03 0,016409075 0,001912945
adecuado
No cumple - Ajuste lineal mal
7 j  ec Re 2,4180E-06 #N/A 5,37156E-07
calculado
8 j  ec Re 2,9584E-09 #N/A 2,17088E-10
adecuado
2
9 5,7233E-05 #N/A 2,11011E-05
calculado
Re 2
j
10

c Re 2  m  2,2949E-04 1,104E-05 0,00022949 Cumple

Tabla 3-4 Análisis de idoneidad de los modelos de regresión para jMT
Análisis de idoneidad
Prueba 2:
Prueba 1:
Modelo Significancia Prueba 3:
Estructura ANOVA
No Valor-p Gráfica residuos
Valor-p
m c
1 j
c Re 1,0827E-04 #N/A 4,40607E-05 calculado
c Cumple,
2 j
Re 2,4354E-04 #N/A 0,000112089 Pero no escogido
2
 c  No cumple - Ajuste lineal no
3 j   4,1129E-07 #N/A 6,83492E-08 adecuado
 Re 
4 j 7,2450E-04 #N/A 0,000391146 calculado
c Re
1 No cumple - Ajuste no lineal con
5 j 5,5772E-07 #N/A 9,74438E-08 heterocedasticidad
c Re 2
6 j  c Re m 5,0651E-04 0,003268351 0,000506506 adecuado

7 j  ec Re 1,1612E-06 #N/A 2,28831E-07 calculado

c Re
8 j e 7,5848E-10 #N/A 4,43808E-11 adecuado

2
9 1,0355E-04 #N/A 4,1853E-05 calculado
Re 2
j
10

c Re 2  m  1,8357E-04 1,39755E-05 0,000183572 Cumple

La Tabla 3-5 presenta el modelo de regresión para jHT y jMT, con los coeficientes obtenidos en el
análisis de regresión. La Tabla 3-6 y Tabla 3-7 y presentan el ANOVA de los modelos
seleccionados para el jHT y jMT respectivamente, como se puede observar, en ambos el valor p es
menor de 0.05; lo que indica que las variables de regresión y la variable de respuesta tienen una
relación significativa y el modelo escogido puede predecir significativamente el resultado. Una
vez comprobada la significancia del modelo con la ANOVA, se analiza por independiente los
parámetros del modelo y su contribución al modelo, la Tabla 3-8 y Tabla 3-9 presentan la prueba
de los coeficientes para el intercepto (m) y el parámetro Re2(c) para el jHT y jMT respectivamente.
Como se observa, el valor-p una vez más es menor que el umbral de aceptación (0.05) por
consiguiente también se valida que los coeficientes de regresión impactan significativamente en el
modelo.
Tabla 3-5 Resultado del modelo de regresión para los factores jHT y jMT
Modelo de regresión del factor j y sus correlaciones

Factor Modelo de regresión
Re 2
jHT
504.563 Re 2  397972890
Re 2
jMT
506.321Re 2  432943428
Tabla 3-6 ANOVA significancia del modelo para jHT
Análisis de Varianza factor jHT

Modelo 75896,1629 1 75896,1629 61,2784 0,00023
Residuo 7431,28091 6 1238,54682
Total 83327,4438 7
Tabla 3-7 ANOVA significancia del modelo para j MT
Análisis de Varianza factor jMT

Modelo 89820,43326 1 89820,43326 66,4144 0,00018
Residuo 8114,540062 6 1352,423344
Total 97934,97332 7
Tabla 3-8 Resultado del análisis de los parámetros del modelo j HT
Análisis de los parámetros del modelo jHT

Error I.C I.C.
Parámetro Coeficiente. Valor-p
típico inferior 95% Superior 95%
c 506,321 39,56 1,398E-05 409,517 603,12
m -432943428 53125118,1 0,00018 -562935909,4 -302950947
Tabla 3-9 Resultado del análisis de los parámetros del modelo j MT
Análisis de los parámetros del modelo jMT

Error I.C I.C.
Parámetro Coeficiente. Valor-p
típico inferior 95% Superior 95%
c 506,321078 39,5617318 1,39755E-05 409,5170072 603,125148
m -432943428 53125118,1 0,000183572 -562935909,4 -302950947
Finalmente, la Figura 3-1 y Figura 3-2figura presentan los gráficos de residuales para jHT y jMT
respectivamente, las figuras demuestran que no existe algún patrón que sugiera estructura u
homocedasticidad, o en su defecto que no se consideraron factores que son importantes para el
modelo de regresión, finalmente con esto hemos validado que el modelo de regresión es idóneo
para predecir y relacionar el Número Re con los factores jHT y jMT respectivamente.
Figura 3-1 Gráfico de residuos estudentizados vs predicho para el factor jH para verificar independencia. Fuente:
elaboración propia
Figura 3-2 Gráfico de residuos estudentizados vs predicho para el factor jM para verificación de independencia.
Estos modelos de regresión presentados, pueden ser utilizados para estimar el comportamiento de
transferencia de calor y masa del aire de un sistema de aire acondicionado convencional con tubos
aleteados durante el proceso de enfriamiento y deshumidificación; ajustando y prediciendo el
coeficiente de transferencia de calor y de masa, adecuadamente bajo las condiciones a las cuales
el sistema experimental. La Figura 3-3 presenta el ajuste de los modelos a las variables
experimentales, aunque se observa una desviación de los datos experimentales y el modelo,
estadísticamente se ha comprobado, que el ajuste es válido para ambos casos.
Figura 3-3 Comparación del modelo de regresión vs datos experimentales. Fuente: elaboración propia.
3.2 Trazabilidad del error.

Como se explicó en el capítulo 1, cada medición experimental consistía en una ejecución de 1 hora
y 45 minutos de toma de datos, en promedio de 330 datos de temperatura, velocidad, humedad
presión y voltaje se obtuvieron por cada medición. El COP es obtenido a partir del calor cedido
por el aire, y el consumo eléctrico del compresor como se expresa en la ecuación (1-1)-(1-3). Las
estimaciones de incertidumbre de las variables de respuestas se expresan a partir de la ecuación
(3-4) y ecuación (3-5):
2 2 2 2
 R2  1 R   1 R   1 R   1 R  (3-4)
  1     2     3   ...    n 
R  R x1   R x2   R x3   R xn 
2 2 2 2 2
 R  R   R   R   R  (3-5)
  1     2     3   ...    n 
R  x2   x2   x3  x
 n 
Siendo R la variable de respuesta, σR la incertidumbre en la respuesta x1, x2, x3, … xn las variables
independientes y σ1, σ2, σ3, …σn las incertidumbres de las variables independientes, es decir las
variables físicas de medición que contienen su propia incertidumbre. Esta ecuación es conocida
como la ley de propagación [58]. La ecuación (3-4) es utilizada cuando la variable respuesta es
presentada en forma de productos entre las variables independientes, la ecuación (3-5) es utilizada
cuando la variable respuesta es obtenida a partir de las sumas de las variables independientes.
Aplicando la ecuación (3-4) en la ecuación (1-1) se obtiene la ecuación (3-6)
2 2
 1   1  (3-6)
 COP    QEvap    WComp 
 QEvap   WComp 
Para poder obtener la incertidumbre del calor cedido en el aire y el trabajo del compresor se utiliza
la ecuación (3-4) y ecuación (3-5) respectivamente; y se obtienen la ecuación (3-7) y ecuación
(3-8)
2 2 2
 1    hair.in    hair.out  (3-7)
Q     m air     
 h h


Evap
m

 air h  h
  air,in air,out   air,in air,out 
W 
V   I1 2  V   I 2 2  V   I 3 2 (3-8)
Comp
1000
Las incertidumbres de las entalpías del aire dependen de las mediciones de la temperatura y
humedad absoluta como se observa en la ecuación (3-9) y la incertidumbre del flujo másico
depende de las mediciones de la velocidad del aire, las mediciones geométricas y de la densidad
del aire ecuación (3-10).
h 
air
1.006 1.86w    25011.86T  
air Tair
2
air wair
2 (3-9)
m  air
U A     A    U   
air cs air
2
air cs U air
2
air air Acs
2 (3-10)
La densidad del aire depende de las mediciones de temperatura y de la humedad como se observa
en la ecuación (3-11)
1
 air  (3-11)
(1  1.607858 wT 10, H 2 )
0.287042 (T10  273.15) 
101.325
Y su incertidumbre se presenta en la ecuación (3-12)
2
  
2

   
 2 
1  
 Tair   
1

 (3-12)
wair 
(1  1.6079  w )
air

 0.287  (Tair  273.15) 
2 air  
  0.287  (T  273.15)  (1  1.6079  wair ) 2
 
   
  101.325    
air
101.325 
La humedad absoluta es calculada con base a la relación psicrométrica presentada en la ecuación
(3-13). Los coeficientes (C9…C13) se presentan en el manual Fundamentals de ASHRAE[8], la
temperatura del aire T está en grados Kelvin y la humedad relativa en porcentaje (%). Debido a lo
extenso de la ecuación, se omite intencionalmente la ecuación de la incertidumbre de la humedad
absoluta (σwin – σwout). Por último, la incertidumbre asociada por la medición de área es presentada
en la ecuación (3-14):
 C  
0.623 exp 8  C9  C10T  C11T 2  C12T 3  C13 log T   HR / 100
wair   T   (3-13)
 C  
101325  exp 8  C9  C10T  C11T 2  C12T 3  C13 log T   HR / 100
 T  
L2 L1 2  L1 L2 2 (3-14)
A  CS
Siendo σL la incertidumbre de la cinta métrica usada para la medición =1mm

En Tabla 3-10, Tabla 3-11 y
Tabla 3-12, se observan los resultados de la incertidumbre de las mediciones; se puede identificar
que los errores asociados directamente a la medición obtienen un valor muy ajustado, esto nos
indica que el error asociado a la instrumentación es muy leve. El error de los datos aumenta
considerablemente conforme se van calculando variables como la entalpía y el flujo másico. Al
final, se observa que el COP contiene una incertidumbre asociada a cada medición realizada entre
el 8 y 15%; podemos considerar que es un error aceptable teniendo en cuenta la cantidad de
fórmulas desarrolladas para poder obtener el resultado. Se observa un error del 6% para las
mediciones de la potencia, también asociada a la fórmula utilizada para obtenerla, no de las
mediciones, ya que no se tomaron medidas directas de la potencia.
Tabla 3-10 Propagación del error en las variables medidas
Incertidumbre variables independientes

Corrida σT_in (°C) σT_out (°C) σwin (kg/kg) σwout (kg/kg) σU (m/2)
1 ± 0,02 ± 0,04 ±5,58,E-15 ±9,78,E-16 ± 0,06
2 ± 0,02 ± 0,03 ±3,40,E-15 ±1,24,E-15 ± 0,07
3 ± 0,02 ± 0,04 ±3,15,E-15 ±1,10,E-15 ± 0,10
4 ± 0,02 ± 0,03 ±5,15,E-15 ±1,38,E-15 ± 0,07
5 ± 0,02 ± 0,04 ±3,49,E-15 ±1,11,E-15 ± 0,07
6 ± 0,02 ± 0,03 ±5,54,E-15 ±1,61,E-15 ± 0,11
7 ± 0,02 ± 0,05 ±3,42,E-15 ±6,93,E-16 ± 0,06
8 ± 0,02 ± 0,03 ±5,28,E-15 ±1,45,E-15 ± 0,07

Tabla 3-11 Propagación del error en las variables calculadas y de respuesta
Incertidumbre variables de respuestas y calculadas

Σρ σmair σhin σhout σQ σW
Corrida σCOP
(kg/m3) (kg/s) (kJ/kg) (kJ/kg) (kW) (kW)
1 ±2,105,E-04 ± 0,01 ± 1,22 ± 1,41 ± 0,30 ± 0,03 ± 0,79
2 ±2,043,E-04 ± 0,01 ± 1,50 ± 1,42 ± 0,37 ± 0,03 ± 0,91
3 ±2,057,E-04 ± 0,01 ± 1,28 ± 1,45 ± 0,50 ± 0,03 ± 1,05
4 ±2,033,E-04 ± 0,01 ± 1,36 ± 1,40 ± 0,34 ± 0,03 ± 0,88
5 ±2,078,E-04 ± 0,01 ± 1,40 ± 1,44 ± 0,39 ± 0,03 ± 0,98
6 ±2,011,E-04 ± 0,01 ± 1,17 ± 1,40 ± 0,52 ± 0,03 ± 1,19
7 ±2,215,E-04 ± 0,01 ± 1,30 ± 1,49 ± 0,32 ± 0,02 ± 0,55
8 ±2,062,E-04 ± 0,01 ± 1,32 ± 1,41 ± 0,35 ± 0,02 ± 0,54

Tabla 3-12 Porcentaje de la propagación del error en las variables
Porcentaje de incertidumbre en variables independientes, de respuestas y calculadas

Corrida σT_in σT_out σwin σwout σU σρ σmair σhin σhout σQ σW σCOP
1 0,06% 0,32% 4,74E-13% 1,218E-13% 3,0% 0,02% 3,79% 2,09% 4,28% 8,21% 6,10% 10,23%
2 0,09% 0,24% 2,40E-13% 1,344E-13% 3,0% 0,02% 4,09% 2,52% 3,74% 10,42% 6,18% 12,11%
3 0,10% 0,28% 2,85E-13% 1,227E-13% 3,0% 0,02% 4,95% 2,50% 4,01% 13,88% 6,08% 15,15%
4 0,07% 0,21% 3,81E-13% 1,455E-13% 3,0% 0,02% 4,02% 2,19% 3,55% 9,49% 6,14% 11,30%
5 0,09% 0,28% 2,7E-13% 1,257E-13% 3,0% 0,02% 4,25% 2,47% 4,01% 10,50% 6,01% 12,10%
6 0,06% 0,18% 5,04E-13% 1,595E-13% 3,0% 0,02% 5,19% 2,08% 3,29% 14,14% 6,28% 15,47%
7 0,09% 0,56% 2,95E-13% 1,006E-13% 3,0% 0,02% 3,83% 2,45% 5,57% 8,45% 2,88% 8,92%
8 0,07% 0,20% 4,06E-13% 1,480-13% 3,0% 0,02% 4,11% 2,16% 3,46% 10,37% 3,08% 10,82%
3.3 Propuesta de estudio: Mapa de funcionamiento del equipo de aire

acondicionado
Los resultados presentados en el diseño de experimento también se pueden comprobar al analizar
la capacidad de enfriamiento del sistema de aire acondicionado, es decir a partir de la Ecuación
(1-2), construir un mapa de funcionamiento considerando diferentes condiciones psicrométricas
(Tin, win), y velocidad del aire (Uin) a la entrada del evaporador, pero manteniendo el aire a un nivel
de salida (Tout, wout). Esta salida se fijó en 13°C y 90% HR, una condición muy utilizada en los
cálculos de psicrometría [18]. Los rangos de variación de entrada del aire al evaporador se
presentan en la Tabla 3-13. Otras consideraciones son explicadas a continuación:
Tabla 3-13 Rango de variación de la temperatura, humedad y velocidad del aire a la entrada del evaporador
Rango de variación de los parámetros de entrada del aire

Parámetro Rango
Temperatura (°C) 10-40
Humedad absoluta (g/kg) 0-30
Velocidad del aire (m/s) 0.5-4.0

 Por limitación física, si el equipo necesitaba más de 4.2 kW de energía para alcanzar la
condición de salida, se omitía el resultado. Si esto se presenta en la vida real (que el aire
necesite una cantidad de energía mayor a la capacidad nominal para alcanzar la temperatura
de suministro), el equipo nunca alcanzará esta condición, y tendrá que conformarse con una
condición de mayor temperatura; debido a esto, es que se omite este resultado.
 Si el resultado obtenido era una capacidad de enfriamiento negativo, se omite el resultado
 Solo se puede trabajar por debajo de la condición del punto de roció, la cual se ve representada
en la curva negra que se observa en cada una de los mapas.
Los 9 mapas que se observan en la Figura 3-4 representan una carta psicrométrica con su respectiva
línea de saturación. Las zonas sombreadas son condiciones de entrada del aire, donde el equipo es
capaz de alcanzar el nivel de salida (13°C y 90%HR), cada zona de colores representa una
capacidad, que va desde 0,1 kW hasta 3.5 kW. Las zonas en blanco representan zonas de
penalización, ya sea porque el equipo no puede alcanzar el nivel de salida (zona derecha), porque
el equipo no tenía la necesidad de hacer proceso de enfriamiento (zona izquierda), y/o se encuentra
en una zona imposible de psicrometría (toda zona por encima de la línea de saturación). Por cada
mapa se presentan los cálculos de QEvap y bajo un nivel de velocidad promedio del aire a la entrada.
Esta presentación se realizó con el objetivo de analizar con más detalle el fenómeno que se observó
en el diseño de experimento. Como se puede observar a lo largo de los mapas, la velocidad del aire
es un factor muy significativo en el mapa de funcionamiento, el aumento de la velocidad, genera
una disminución del área de acción en el equipo de aire acondicionado; si la unidad trabaja con
una velocidad baja (1,5-2.2 m/s), y el aire ingresa al evaporador con una condición de 23°C, 12
g/kg el equipo trabajará más eficientemente en comparación si su velocidad trabajara con un alto
rango (2,5-3 m/s) ya que a esos niveles, se penalizaría la condición de diseño, por lo tanto el equipo
reducirá su capacidad de enfriamiento más para alcanzar una condición de salida, en una reducción
del COP, ambas situaciones se observan en la Figura 1-11 y Figura 1-14.
Figura 3-4 Mapa de la capacidad de enfriamiento de la unidad de aire acondicionado tipo mini split en una tabla
psicrométrica a diferentes velocidades. Fuente: elaboración propia
Si observamos otra condición más conservadora (20°C, 10 g/kg a la entrada del evaporador) se
observa si el equipo opera a altas velocidades que su mapa de funcionamiento indica que es capaz
de rechazar entre 1,8 y 2,5 kW, esto representa un menor flujo de refrigerante, y gracias al inverter,
el compresor trabajará menos, reflejando esto en un mejor COP del sistema. Por otra parte, si el
equipo opera con una velocidad baja, el mapa indica un rechazo entre 0,5 y 1 kW, un rango muy
bajo de rechazo de calor que podría reflejarse en una penalización por parte del inverter, ya que
este no será capaz de controlar la temperatura a la salida (fluctuará entre la temperatura de control
y una menor) y en consecuencia un mayor consumo del compresor y una disminución del COP,
estos dos casos también se pueden observar en Figura 1-11 y Figura 1-14.

Para esta sección podemos destacar las siguientes conclusiones:
 Se obtuvo un factor de Colburn para la transferencia de calor y uno para la transferencia de
masa en evaporadores de aire acondicionado tipo tubos aleteados, que trabajan en condiciones
de clima tropical en el proceso de enfriamiento y deshumidificación.
 Estos coeficientes fueron obtenidos teniendo en cuentas las condiciones y limitaciones
operacionales del sistema de aire acondicionado. A diferencia de otros factores que se analizan
bajo sistemas experimentales modificados y usando agua fría como refrigerante en lugar de un
refrigerante con cambio de fase.
 Los resultados de los coeficientes se ajustan a lo establecido por las teorías de transferencia de
calor: A mayor velocidad, disminuye la transferencia de calor y de masa [8], [18], [19], [23]
 Se realizó un análisis de trazabilidad del error para el cálculo del Q Evap, WComp y el COP, se
observa que el COP contiene una incertidumbre asociada a cada medición realizada entre el 8
y 15%; podemos considerar que es un error aceptable teniendo en cuenta la cantidad de
fórmulas desarrolladas para poder obtener el resultado.
 Se observa un error del 6% para las mediciones de la potencia, también asociada a la fórmula
utilizada para obtenerla, no de las mediciones, ya que no se tomaron medidas directas de la
potencia.
 El uso del mapa de funcionamiento podría influir en el estudio de sistemas de
acondicionamiento de aire tipo residencial para condiciones de clima tropical. Con este
esquema es posible diseñar estrategias de control orientadas al uso racional de la energía y al
confort térmico.
4 Modelo dinámico del evaporador
Resultados previos, identifican que la humedad es un factor importante en el proceso de
acondicionamiento de aire; hemos visto que el cambio del mismo puede afectar el consumo y el
rendimiento operacional de la unidad. Por otra parte, hemos observado que la transferencia de
calor del aire hacia el refrigerante también es afectada por la humedad y otros factores, como la
velocidad y la temperatura del aire al ingreso del evaporador. Las limitaciones del equipo para la
deshumidificación del aire se deben al área de transferencia de calor/masa, y esta área depende de
la capacidad del refrigerante de absorber calor durante su fase de mezcla saturada. Por
consiguiente, es necesario identificar físicamente el área ocupada por el refrigerante en estado
saturado dentro del evaporador, y así, estudiar más a fondo el proceso de deshumidificación desde
la interacción aire-refrigerante, y por esto, se plantea desarrollar un modelo dinámico del
evaporador.
Además de lo mencionado, el modelo dinámico nos permitirá proponer nuevas estrategias de
control asociadas a sistemas de aire acondicionado DX tipo residencial y así, poder alcanzar un
mejor desempeño de estos sistemas bajo condiciones de clima tropical.
El modelo dinámico del evaporador no puede ser un modelo de tipo “caja negra” ya que es
necesario establecer parámetros internos del equipo, además que el fenómeno del intercambiador
de calor es de naturaleza no lineal; tampoco puede ser un modelo tipo “caja blanca” porque
desarrollar los principios de conservación de masa, energía y momento supondría un costo de
procesamiento de cálculo muy elevado y por consiguiente costoso. La solución utilizada
generalmente para implementar estrategias de control en aire acondicionado por lo general es un
modelo de “caja gris”, este tipo de modelo trabaja con principios básicos pero el modelo es
ajustado experimentalmente. Mucho de los modelos dinámicos de aire acondicionado que se han
investigado para ajuste de control, tuning, diagnósticos de detección de fallas, etc. utilizan este
tipo de modelo [36].
4.1 Formulación matemática

Para desarrollar un modelo dinámico del evaporador a partir del método de límites movibles, en el
cual este estudio se basa, se tiene en consideración la metodología presentada por Rasmussen [30],
[32], [36] teniendo en cuenta las siguientes consideraciones:
 El evaporador es modelado como un tubo largo circular cuya longitud, es mucho mayor
que el diámetro de la tubería.
 El refrigerante que fluyen en el interior de la tubería es modelador solo en dirección
longitudinal
 La conducción de calor en el refrigerante es despreciable
 El espesor de la tubería se asume como despreciable; el perímetro del tubo interno y externo
son iguales y la resistencia por conducción a través de las paredes también es despreciable
 La caída de presión a lo largo del evaporador es considerada despreciable. Por lo tanto, la
presión en este evaporador se asume espacialmente uniforme.
 El aire fluye en dirección normal a la tubería del refrigerante, además se considerará los
efectos de la humedad en el proceso de transferencia de energía.
El estudio presentado por Rasmussen presenta el aire fluyendo en dirección paralela a la superficie,
y considera el aire como una frontera donde todo el calor que cede el aire se debe a la diferencia
de temperatura. Para este estudio, la última consideración fue modificada, ya que queremos
asegurar que el modelo simule una condición semejante a lo que ocurre realmente en el evaporador
en climas tropicales.
La ecuación (4-1), (4-2) y (4-3), son aquellas que gobiernan la modelación de límites movibles,
representan la conservación de masa, conservación de momento y de energía respectivamente.
A m  (4-1)
 0
t z
Ah  PA m h (4-2)
  pi Tpared  Ti 
t z
ACS C pTpared  (4-3)
 p i Ti  Tpared   p HT Ta  Tpared   p MT h fg wa  wpared 
t
La ecuación (4-3) fue modificada para este modelo, a fin de tener en consideración el proceso de
deshumidificación del aire. Estas ecuaciones se implementan en el refrigerante y en la pared del
tubo por cada región. Para poder eliminar la dependencia espacial, las ecuaciones son integradas
a lo largo de la tubería, permitiendo que el modelo se simplifique a una ecuación 1D. La
integración de esta ecuación parcial diferencial se realiza usando la regla de Leibniz (ecuación
(4-4)) y usando como límites de la integración, la longitud de cada región:
f z, t 
z2 ( t )
d 2
z (t )  d z t , t  d z t , t  (4-4)
z (t ) t dz    f  z , t dz   f z2 t , t  2  f z1 t , t  1
1
dt  z1 (t )  dt dt
Al aplicar la integración finalmente se obtiene ecuaciones diferenciales ordinarias las cuales son
organizadas en forma de matriz para obtener una respuesta en forma de “espacio de estado”
El evaporador está bajo condiciones de operación normal, desarrollará dos regiones del
refrigerante: región bifásica y región sobrecalentada, así que en cada región se desarrollan las
ecuaciones (4-1) - (4-3) y aplicando la ecuación (4-4) dando como resultado seis ecuaciones
diferenciales ordinarias (4-5)-(4-10):
 d f d 
 1     g   ACS L1Pi   f   g 1   ACS L1  m i,e  m mid (4-5)
 dPi dPi 
 d  1  d  dhg  1   
 2   2   ACS L2 Pi   2  A L h     A L  m
 mid  m
 i,s (4-6)
 dPi  2  dh2  dPi  2  h2  CS 2 i , s 1 g CS 1
 h2   Pi   Pi 
 d  f h f  d  h  
 1     g g   1 ACS L1Pi   f h f   g hg 1   ACS L1 (4-7)
 dPi dPi 
 m i ,e hi ,e  m mid hmid   i ,1 Ai 1 Tpared,1  Ti ,1 
 d   1  dhg    
 2
      2  h   1  dhg    1 A L P
 dPi   2  dPi  h2  2  2  dPi  2  CS 1 i
h2   Pi 

(4-8)
    1 
   2  h2   2   ACS L2 hi ,s   g hg   2 h2 ACS L1
 h2 Pi   2 
 
 m mid hmid  m i ,s hi ,s   i , 2 Ai 2 Tpared, 2  Ti , 2 
C p V  Tpared,1   i,1 Ai Ti,1  Tpared,1    HT AA Ta  Tpared,1 
pared
(4-9)
  MT AA h fg wa  w pared,1 
  T pared, 2  T pared,1   
C p V  pared T pared, 2  
 
L2
 L1    i , 2 Ai Ti , 2  T pared, 2 
(4-10)
   
  HT AA Ta  T pared, 2    MT AA h fg wa  w pared, 2 
Lo primero que se observan en las ecuaciones son las derivadas ordinarias (y parciales) de las
propiedades termodinámica del refrigerante. La obtención de estas derivadas puede realizarse ya
sea a través de diferencias finitas, por relaciones termodinámicas h  u  P   o también a partir
de software de propiedades termodinámicas (EES, COOLPROP, REFPROP). Los subíndices 1,2
son utilizadas para identificar las regiones que está dividido el modelo (1: zona bifásica, 2: zona
gaseosa o sobrecalentado); los subíndices i, a, indican el tipo de fluidos (i: refrigerante, a: aire).
Los subíndices f, g, i,e, mid, i,s indica la fase y ubicación del refrigerante (líquido saturado, vapor
saturado, entrada, intermedio y salida respectivamente). La variable λ es la longitud normalizada
de cada región, 1  L1 LT  y la variable  es la fracción de vacío promedio, esta propiedad
desarrollada por Wedekind [59] es utilizada para modelar la zona bifásica en intercambiadores de
calor. La fracción de vacío γ0 es la relación entre el volumen de vapor y el volumen total de una
zona puntual que se está estudiando/evaluando (No se debe confundir con la calidad del vapor
“X”, ya que esta última es la relación entre masa de vapor y masa total).
La fracción de vacío promedio es la integración de la fracción de vacío a lo largo de la zona bifásica

del refrigerante en el evaporador, desde que ingresa hasta que sale (si no se evapora completamente
el refrigerante) o cuando el refrigerante alcance la fase de vapor saturado (X=1). Existen diferentes
estudios donde desarrollan correlación para obtener la fracción de vacío y la fracción promedio
como se resume en [60], también esta propiedad puede ser considerada dinámica, sin embargo,
teniendo en cuenta que el sistema experimental es un equipo de pequeña capacidad térmica, la
respuesta dinámica del equipo suele ser muy pequeño y la fracción de vacío dinámico tiende a ser
más pequeño[30], por lo tanto es razonable considerar esta propiedad como invariante. La ecuación
(4-11) y ecuación (4-12) presenta la fracción de vacío local y fracción de vacío promedio utilizados
en esta tesis:
2
X  v  3 (4-11)
0   S   
X  1  X S  f 
 
 1   S   S  X 1  S   (4-12)
      2
ln  
 1  S   1  X 1  S    S  X 1  S   1  X 1  S 
Por último, se observa que las ecuaciones (4-5)-(4-10) contiene 5 derivadas dinámicas
 
L1 , Pi , hi ,s , Tpared,1 , Tpared, 2 las cuales deben resolverse y para esto, se combina algebraicamente las
ecuaciones para eliminar la variable m
 mid , ya que esta variable es muy difícil de determinar; por
lo tanto se organiza las ecuaciones en una matriz 5x5 de forma Z xx  f x, u  como se observa
en la ecuación (4-13); la matriz Z(x) contiene las propiedades del refrigerante, ẋ el vector de las
derivadas dinámicas, mientras que el vector f(x,u) contiene las variables de entrada y estructura
del sistema.
 z11 z12 0 0 0   L1   m i ,e hi ,e  hg   i ,1 Ai 1 Tpared,1  Ti ,1  

z
 21 z22 z23 0  
0   Pi
 
  m i,s gh  h i,s    A  T
i , 2 i 2 pared, 2  Ti,2  

 z31 z32 z33 0 0  hi , s    m i ,e  m i , s  (4-13)

0 0 0 z44
   

0  Tpared,1   i Ai Tpared,1  Ti ,1   AA  HT Ta  Tpared,1    MT h fg wa  wpared,1    
 z51 0 0 0  
z55  Tpared, 2  i Ai Tpared, 2  Ti , 2   AA  HT Ta  Tpared, 2    MT h fg wa  wpared, 2  
La tabla presenta los elementos de la matriz Z

Tabla 4-1 Tabla de los elementos de la matriz Z(x,u) para el evaporador
Elementos de la matriz Z y ecuación a resolver

Elemento Fórmula
z11  hf f  hg  1   Acs
 d  f h f  d f   d  h  d  
z12   h 1     
g   dP g g
hg   1 Acs L1
g
 dPi dPi   i dPi  

z21 2 hg  h2 Acs
    1     dh   dhg 
z22


2   
  2  g
h2  hg    
 2  
1 Acs L2
 Pi h2   2  h2 Pi  dPi   2  dPi 
  
 1     
z23
 
  
2 h  h    2  Acs L2
 2  h2 Pi 
2 g
  2 
  
z31  g  2    f   g 1    Acs
    1     dh   d g   
z32


2     2  g  L   d f 
1     

  L1  Acs
 Pi h2   2  h2 Pi  dPi 
2 

   dPi   dPi    
 1   2  
z33   ACS L2
 2  h2 Pi 
z44 C V 
p pared
 Tpared, 2  Tpared,1 
z51 C V 
p 
pared 
L

 2 
z55 C V 
p pared

Además de modelar el evaporador, también se analizará el comportamiento del aire a través del
evaporador, esto con el propósito de calcular las condiciones de salida del mismo, bajo las
variables de salida del modelo del evaporador.
4.2 Modelación del aire húmedo

En esta parte, un modelo estático puede describir fácilmente el comportamiento del aire, ya que la
simplicidad del modelo dinámico genera que el aire no tenga acumulación o variación dinámica.
Para describir las condiciones de salida del aire, como la temperatura y la humedad absoluta, se
implementa el método ε-NTU bajo la modelación presentada por Corberán et al para aire húmedo
[61] las ecuaciones (4-14) y (4-16) presentan la obtención de los parámetros por zona “x”:

 MT PLx
  MT x 
 PL
ha, s _ x  ha,ie m aireb  hpared, x 1  e m aireb  (4-14)

 
 
 PL
 HT x   HT x 
 PL
Ta, s _ x  Ta,ie aire p ,aire  Tpared, x 1  e aire p ,aire 

m C m C (4-15)
 
 
Qaire  Qsensible
wa.s _ x  wa,i  (4-16)
m aireh fg
Finalmente, como el aire fluye perpendicular a la dirección del refrigerante, las condiciones de
salida por región son diferentes, y la salida final del aire en el evaporador corresponde a su
promedio ponderado.
4.3 Solución del modelo dinámico
La Figura 4-1, presenta el procedimiento del modelo dinámico, las condiciones de entrada del aire
son ingresados a través de MATLAB; calculando algunas propiedades y estructura del modelo;
para el cálculo de propiedades del refrigerante se incluye una rutina de REFPROP dentro del
programa.
Figura 4-1 Esquema del modelo dinámico del evaporador y aire. Fuente: elaboración propia
El modelo dinámico se realiza a través de un programa de SIMULINK, manteniendo

comunicación con MATLAB para el cálculo de los elementos de la matriz Z y del vector
independiente f. La matriz del evaporador se resuelve y luego ingresa a un segundo programa en
MATLAB para la solución del modelo del aire. Finalmente, la salida del modelo dinámico ingresa
una vez más a la entrada del modelo para iniciar de nuevo el modelo en un nuevo tiempo, con
nuevas condiciones de entrada.
En la Figura 4-2 se observa el programa realizado en SIMULINK, donde se observa las variables
de entradas, las rutinas y las variables de salida para esta modelación.
Figura 4-2 Modelo dinámico del evaporador en SIMULINK. Fuente: elaboración propia
Para resolver el modelo dinámico se utilizó el solver ode23t, la elección de este solver es debido a
su idoneidad para resolver ecuaciones rígidas, es decir, sistemas de ecuaciones donde algunos
componentes de la ecuación cambian rápidamente, junto con componentes que cambian
lentamente[62], este tipo de ecuaciones no son compatibles para su solución implementando el
método Runge-Kutta, ya que requieren que el paso del tiempo sea muy pequeño. El ODE23t es un
método implícito en la solución de ecuaciones diferenciales ordinarias que utiliza el método
trapezoidal (basado en la regla trapezoidal de integración) para resolver este sistema de ecuaciones
rígidas[62].
4.3.1 Variables dinámicas
Para la formulación de este modelo dinámico, todas variables utilizadas se presentan a
continuación:
 Variables de entradas (señal utilizada para evaluar):
o Flujo másico del refrigerante (Sinusoidal)
o Presión y temperatura del refrigerante a la salida del condensador (Constante)
o Temperatura y Humedad relativa del aire a la entrada del evaporador (Rampa)
o Velocidad frontal del aire (Sinusoidal)
 Estructura del sistema y parámetros/propiedades:
o Geometría del evaporador (área total, área interna, área externa)
o Masa del evaporador
o Propiedades térmicas del evaporador
 Variables de salida: (Esta a su vez se divide en variables del evaporador y del aire)
o Longitud del evaporador ocupada por la zona 1
o Presión del refrigerante en el evaporador
o Entalpía del refrigerante en la zona 2 del evaporador
o Temperatura de la pared del tubo en la zona 1
o Temperatura de la pared del tubo en la zona 2
o Temperatura y humedad del aire a la salida del evaporador
4.3.2 Condiciones iniciales
Para solucionar la ecuación (4-13) es necesario determinar condiciones iniciales que permitan
describir el modelo dinámico y obtener unas variables a las salidas coherentes, en este caso, son 5
parámetros dinámicos: Longitud zona 1, Presión evaporador, entalpía refrigerante zona 2,
Temperatura pared tubería zona 1 y temperatura pared tubería zona 2. Las condiciones iniciales se
escogen a partir de las variables geométricas del evaporador, además de los parámetros que se
observaron en 1 de las 8 corridas experimentales, la Tabla 4-2 presenta las condiciones iniciales
geométricas, intermedias y de salida escogida para este modelo dinámico. Para validar estas
condiciones, primero se ejecuta el modelo en estado estable y se verifica que las variables
dinámicas del modelo no existan variación en el tiempo. Segundo, con los resultados en estado
estable se realiza la modelación del aire húmedo. Finalmente se compara los resultados del modelo
dinámico con los datos experimentales. Para este trabajo se escoge la corrida experimental 5 para
la ejecución validación del modelo; se escoge esta corrida experimental, ya que es entre todos los
resultados, la del mejor COP registrado.
Tabla 4-2 Condiciones iniciales de las variables de entrada, dinámicas y de salida del modelo dinámico
Condiciones iniciales de las variables del modelo

Parámetro Descripción
Corrida Experimental 5
Nivel de velocidad del ventilador Nivel 1
Condiciones experimentales
Temperatura entrada 23
Humedad Relativa entrada 75
LTotal (m) 17.98
Ai (m2) 0.372
Condiciones geométricas y AA (m2) 7.8098
propiedades del evaporador ACS (m2) 3.421×10-5
CP (kJ/kgK) 0.6252
(ρV)Evaporador (kg) 6.5389
href,in (kJ/kg) 238,549
ṁref (kg/s) 0,019813
Condiciones iniciales
Ta,i(°C) 23,72
de entrada
Ua,i (m/s) 2,45
HRin (%) 70,34
L1 (m) 15 < L1 < 17
Pref,Evap (kPa) 610<Pref,Evap<680
Condiciones iniciales:
href,out (kJ/kg) 420 < href,out <450
variables dinámicas
Tpared,1 (°C) 0< Tpared,1 <19
T pared,2 (°C) Tpared,1 < T pared,2 < Ta,i
Pref (kPa) 624.97
Condiciones de salida Tref,out (°C) 9.89
experimentales Ta,o(°C) 13,469
wa,o(kg/kg) 0,0088
4.4 Resultados
La Tabla 4-3 presenta los resultados de las variables dinámicas y de salida del modelo en estado
estable, como se puede observar, las condiciones de salida del aire son similares a las obtenidas
durante el diseño experimental con una diferencia de 5% en la temperatura de salida del aire en el
modelo en relación al resultado experimental. Por otro lado, la temperatura a la salida del
refrigerante en el modelo (4.7°C) varía de los resultados experimentales (9.89°C), esto se debe a
la variación entre la presión obtenida en la simulación y la corrida experimental cuya diferencia es
de 6.2%, teniendo en cuenta que hay una diferencia de 1.3% en los resultados de la entalpía (href,out
= 436.24 kJ/kg experimental). es posible asumir que la variación de la temperatura se presente por
la presión. Esta variación puede ser ajustado si se incluye el modelo completo del sistema de
acondicionamiento de aire.
Tabla 4-3 Condiciones de estado estable de las variables dinámicas de entrada
Condiciones variables dinámicas Condiciones de salida

estado estable estado estable
Pref Tref,out Ta,o wa,o
L1 (m) Pref,Evap (kPa) href,out (kJ/kg) Tpared,1 (°C) T pared,2 (°C) (kPa) (°C) (°C) (kg/kg)
16.67 664.26 430.28 11.26 15.37 664.26 4.7 12.8 0.0088
El modelo trabaja y se ajusta adecuadamente a la dinámica del proceso como se observa en las
figuras, la Figura 4-3 presenta las variables dinámicas en estado estable temperatura y humedad
absoluta del aire a la salida del evaporador. La Figura 4-4 presenta el comportamiento de las
variables de estado del evaporador en estado estable. El evaporador en estado dinámico, se presenta
en la Figura 4-5 y Figura 4-6; la simulación inicia bajo condiciones de estado estable, luego en un
t = 500 la temperatura de entrada del aire disminuye, luego en t = 750 disminuye la humedad
absoluta del aire a la entrada del evaporador. Como se puede observar, la dinámica del sistema
responde adecuadamente a los cambios de condiciones de entrada del aire, a excepción de la
variable dinámica Temperatura de la pared en la zona 2, el cual inicialmente tiene una pequeña
elevación antes de realizar su respuesta natural (disminuir su temperatura) esto se debe a la
naturaleza de la solución del problema y debería disminuir en el momento que sea realice una
depuración al programa.
Figura 4-3 Superior. T aire entrada y T aire salida vs tiempo de simulación. (Inferior) T aire salida vs w aire salida
vs tiempo de simulación Fuente: elaboración propia
Figura 4-4 Variables espacio estado del evaporador y condiciones psicrométrica del aire a la entrada del evaporador
vs tiempo de computación. Fuente: elaboración propia.
Se destaca en los resultados que la zona 1, aumenta su longitud y reduce la presión, conforme
disminuye la temperatura del aire a la entrada del evaporados; y con la disminución de la humedad
se presenta una disminución para la longitud y un aumento de presión. Se recomienda desarrollar
un análisis más detallado de la humedad teniendo en cuenta resultados experimentales con
mediciones intermedias en el equipo para tener datos de avance de las propiedades en posición y
tiempo.
Figura 4-5 Superior. T aire entrada y T aire salida vs tiempo de simulación. (Inferior) T aire salida vs w aire salida
vs tiempo de simulación bajo perturbación Fuente: elaboración propia
Figura 4-6 Variables espacio estado del evaporador y condiciones psicrométrica del aire a la entrada del evaporador
vs tiempo de computación bajo perturbación. Fuente: elaboración propia.
Se desarrolló un modelo dinámico del evaporador a través del método de límites movibles, con el
propósito de desarrollar una estrategia de control o un análisis de optimización del proceso.
 El modelo dinámico presentado incluye los efectos de la humedad en el análisis energético
del evaporador, análisis que no se tiene en cuenta en la mayoría de los modelos dinámicos
desarrollados a través del método LM.
 El modelo presenta un ajuste adecuado en cuanto a las variables de salida del aire, pero con
fallas en las variables de salida del refrigerante. El modelo debe ser ajustado teniendo en
cuenta condiciones experimentales que consideren mediciones intermedias para poder tener
una mejor aproximación del fenómeno a nivel local.
 Se observa que la zona bifásica y un comportamiento esperado, sin embargo, es necesario la
inclusión de los otros dispositivos para un mejor análisis de los resultados.
Conclusiones
Esta tesis se realizó con el objetivo de estudiar el aire de un clima tropical y su influencia sobre
los equipos de acondicionamiento tipo residencial, a partir de sus características energéticas, de
transferencia de calor y de masa. A fin de conocer la relación actual aire/equipo e identificar
oportunidades de ahorro energético y de optimización del sistema. A partir de todo los resultados
y los análisis se presentan las siguientes conclusiones:
 Las características del aire en la región tropical influyen sobre las variables energéticas del
sistema de acondicionamiento de aire tipo residencial, y el control de estos equipos no debería
orientarse únicamente hacia la temperatura; pues ambas variables tienen influencia la una sobre
la otra. Esto se aprecia al momento de evaluar el factor de transferencia de calor y el de masa,
ya que debieron ser calculados de manera iterativa. Aún más, la velocidad del aire que ingresa
al evaporador influencia a ambas.
 Es posible ser eficiente y alcanzar las condiciones de confort siempre y cuando se establezca
una velocidad de aire adecuada a partir de las condiciones psicrométrica del aire a la entrada.
Para esto es necesario desarrollar una estrategia de control de la velocidad del aire; que junto
con la estrategia actual de control de temperatura se pueda orientar el control de un sistema de
acondicionamiento de aire hacia el desempeño eficiente y condiciones adecuadas de confort
 La transferencia de calor del aire húmedo hacia el evaporador varía a partir de las condiciones
psicrométricas del aire cuando ingresa al evaporador. A bajas velocidades, la transferencia de
calor disminuye si aumenta la humedad relativa. A altas velocidades la transferencia de calor
se incrementa si aumenta la humedad relativa. Los resultados analizados sugieren que la
velocidad del aire dentro del evaporador influye en la capa de interface agua-vapor de agua;
cuando el gradiente de temperatura entre el agua y el aire es muy pequeño, la capa límite es
dominada por la capa límite de la velocidad, un incremento de la velocidad ocasiona un
aumento de la capa límite y la interface, esto puede influye en la capa límite térmica,
disminuyendo en el fenómeno de transferencia de calor, ocasionando que la transferencia sea
ineficiente. En cambio, cuando existe un gradiente de temperatura alto, la capa límite es
controlada por la capa límite térmica; un aumento de la velocidad aumentará la capa límite de
la velocidad, y este se ajustará a las condiciones de transferencia de calor permitirá que el
fenómeno sea mejor para enfriar el aire.
 La eficiencia energética de un sistema de acondicionamiento de aire tipo residencial depende
de la relación temperatura-humedad-velocidad vs coeficiente global de transferencia de calor
y masa.
o Cuando estos tres factores se encuentren por debajo de 23°C-10 g/kg y un rango de 1.5-
2.0 m/s, el coeficiente global de transferencia de calor y masa es alto (0.15 kg/m2s) y
debido a esto, el calor es absorbido rápidamente y por consiguiente existe un exceso de
energía consumida por el compresor, dando como resultado un COP bajo; para mejorar
la eficiencia energética es necesario aumentar la velocidad del aire, para disminuir este
coeficiente.
o Cuando estos tres factores se encuentren en una temperatura de 28°C, 14 g/kg y un
rango de 2.5-3 m/s el coeficiente global es el mínimo, esto causará que exista un gran
consumo energético, debido al fenómeno de condensación del vapor de agua en
presencia del aire como gas no condensable, la velocidad desprende la capa de la
interface influyendo en el proceso de transferencia de calor y en consecuencia, el
compresor deberá trabajar más para poder absorber el calor del aire dando como
resultado un COP bajo; se recomienda en estos casos disminuir la velocidad del aire,
para aumentar ligeramente el coeficiente global de transferencia de calor y masa y
mejorar el desempeño energético
Es necesario estudiar más a fondo este coeficiente global de transferencia de calor y masa, a
fin de establecer estrategias para ahorrar energía, ya sea a nivel de control avanzado, como el
nivel de mejoras del evaporador diseñado con orientación hacia la transferencia global y el
ahorro energético.
 La deshumidificación del aire está ligada al coeficiente de transferencia de masa del vapor de
agua, sin embargo, este coeficiente no influye significativamente en el consumo y/o eficiencia
energética del sistema de acondicionamiento de aire. Aunque cabe destacar, que durante el
mejor punto de operación energética, el coeficiente global de transferencia de calor y masa no
sufrió una variación significativa, como en los casos donde existe un mayor consumo
energético (ver Figura 2-4). Esto parece indicar, que cuando el coeficiente global de
transferencia de calor y masa es independiente de las condiciones psicrométrica, se convierte
en un indicador de que el consumo energético del sistema es adecuado. Ya que la diferencia
de entalpía entre el aire y el refrigerante podrían ser óptimo. Es recomendable profundizar un
poco más sobre esta relación.
 Los coeficientes de transferencia de calor y de masa calculados permitirán realizar diferentes
análisis al proceso de enfriamiento y deshumidificación del aire en el evaporador, a fin de
desarrollar otros estudios relacionados al tema.
 El modelo dinámico nos permite simular la interacción aire/refrigerante dentro del evaporador,
esto no era posible realizarlo a través de un procedimiento de cálculo estacionario, ya que
existe una inercia térmica en el evaporador. A su vez permite desarrollar futuros trabajo de
estrategias de control avanzado.
Trabajos futuros
Un fin de este trabajo fue establecer un punto de partida en el desarrollo de nuevas oportunidades
de ahorros energéticos; por lo observado a lo largo de esta investigación, se identificó que la
regulación de la velocidad del ventilador puede mejorar en el desempeño energético de equipos
tipo residencial, siempre bajo ciertas condiciones de temperatura y humedad absoluta. Es por ello,
que una oportunidad de ahorro energético puede ser el desarrollo e implementación de una
estrategia de control a través de la medición de la temperatura y humedad del aire y la regulación
de la velocidad del compresor y del ventilador en el evaporador. Otra estrategia de control a tener
en consideración es el desarrollo experimental del mapa de funcionamiento de un sistema de aire
acondicionado tipo residencial; con este mapa, es muy probable aportar en el desarrollo de sistemas
de control inteligente a través de la regulación de las variables internas del evaporador y desarrollar
nuevas estrategias de control hacia la adecuada operación y una alta eficiencia energética.
Uno de los fenómenos que se encuentra en este documento y en ASHRAE [8] es la disminución
del coeficiente de transferencia de calor debido a la velocidad a la cual ingresa el aire al
evaporador, sabemos que esto se debe al desprendimiento de la capa límite del agua condensada y
la formación de gotas en el aire que generan una disminución del proceso de transferencia. Realizar
un análisis más profundo de este fenómeno, específicamente sobre la velocidad efectiva de
transferencia de calor, podría arrojar en una estrategia de ahorro energético, orientada al control
de la velocidad del aire que ingresa al evaporador.
En este estudio se consideró los rangos de variación entre los niveles de los efectos principales
eran lo suficientemente corto para no asumir curvatura; sin embargo, se observa en los resultados
que las superficies de respuestas (COP, Wcomp) no tienen una tendencia lineal definida y por
consiguiente podría existir curvatura. Se recomienda establecer un nuevo diseño de experimento
con puntos centrales para identificar si existe curvatura en las variables de respuestas; esta vez
incluyendo como factor principal la humedad relativa en lugar de la humedad absoluta. A su vez,
seleccionar un factor de bloque que permita confundir la interacción ABC en lugar de las
interacciones AB, AC y BC
Se deberá profundizar más en el coeficiente global de transferencia de calor, a modo de analizar

su afectación en el proceso de la unidad de aire acondicionado al incluir otras variables como la
geometría del evaporador (Diámetro, espesor aleta, paso aleta, etc.) y variables adimensionales
como el número de Jakob (Ja).
También se recomienda la implementación de los coeficientes de transferencia de calor y masa en
la aplicación de técnicas como el diseño de evaporadores a fin de realizar un estudio de
optimización del evaporador.
El modelo dinámico del evaporador desarrollado en este trabajo está orientado para simular
estrategias de control; no obstante, este modelo podría ser utilizado para estudiar más a fondo los
efectos del tamaño en el evaporador (geometría, forma, configuración, etc) en el desempeño del
proceso de enfriamiento y deshumidificación del aire. Para esto, es necesario modificar el modelo
dinámico e incluir las variables geométricas que se desean evaluar (Diámetro int/ext, tipo de aleta,
material, paso longitudinal/transversal, etc.) y así, ejecutar un análisis numérico para predecir y
calcular las características operacionales del evaporador, con el fin de encontrar una configuración
de evaporadores que puedan mejorar el proceso de deshumidificación y manteniendo el buen
desempeño energético en el sistema de aire acondicionado tipo residencial
El desarrollo de nuevos evaporadores orientado hacia el ahorro energético perjudica el proceso de
deshumidificación, no obstante, así como se planteó desarrollar estrategias de control para el
mejoramiento del desempeño energético a través de la medición de temperatura y humedad,
Anexos
I. Selección del diseño de experimento
a. Factores
La naturaleza del aire húmedo es una mezcla de aire seco y vapor de agua (humedad) [8]; y por lo
tanto, cualquier proceso psicrómetro donde el aire húmedo es sometido, se afectan tanto la
temperatura como la humedad. Por consiguiente, dos factores a considerar en el diseño será la
temperatura y la humedad del aire a la entrada del evaporador. Otro factor a considerar dentro del
proceso psicrométrico es el flujo másico del aire que fluye en el evaporador; modificar el flujo
másico del aire modifica la capacidad de absorción de calor Qevap en la unidad y por consiguiente
el COP. Sin embargo, este factor no puede ser controlable directamente, el factor que si puede ser
controlable es nivel de velocidad del ventilador del evaporador, una variación del nivel de
velocidad del ventilador afectará la velocidad del aire y por efecto, el flujo másico.
b. Variables de respuesta
Las variables de respuestas, para este estudio son los parámetros energéticos de la unidad de
acondicionamiento de aire. Ciertamente se está evaluando el aire y sus efectos en el evaporador,
pero considerar las condiciones psicrométricas del aire a la salida del evaporador, puede conllevar
a un error de análisis. Después de todo, en un equipo de acondicionamiento de aire tipo residencial
existen factores incontrolables por parte del usuario que influyen directamente en el proceso de
enfriamiento y deshumidificación del aire (sistemas de control, regulación del inverter, etc.), por
consiguiente, considerar la salida del aire en el evaporador como variable de respuesta afectará los
resultados. En el caso de los parámetros energéticos, son variables que dependen directamente de
los factores escogidos para este diseño de experimento. Para este estudio, se consideró tres
variables de respuestas: el flujo de calor que el evaporador absorbe del aire (Qevap), la energía
requerida para que se realice la operación del equipo (Wcomp), y finalmente el coeficiente de
desempeño (COP) [8]. El nombre de la variable utilizada para representar la eficiencia de energía
externa en un DX se denomina relación de eficiencia energética (EER), sin embargo, referirse a
este este indicador infiere en que las pruebas fueron realizadas bajo una norma o estándar de
pruebas de desempeño (AHRI, ASHRAE, ISO, etc.) y bajo condiciones controladas de laboratorio.
Desarrollar el experimento bajo esas condiciones inclinaría el diseño hacia el equipo de aire
acondicionado y dejará el aire de clima tropical como un segundo plano de análisis.
c. Niveles factor
La temperatura es un factor continuo; el nivel bajo representa la mínima condición de temperatura
donde la unidad operaba estable, el nivel alto representa el clima típico de la ciudad. La humedad
absoluta, también es factor continuo; el nivel alto, representa condiciones típicas de clima de la
ciudad; el nivel bajo representa un caso de operación extrema. Por último, el nivel de velocidad
del ventilador en el evaporador, es un factor categórico, ya que su manipulación es representada
en el nivel de velocidad deseada a partir del control de la unidad de aire acondicionado. La unidad
contiene 4 niveles de velocidad (Auto, bajo, medio, alto) y para este experimento, solo se consideró
dos niveles (bajo y alto).
Se observa en Tabla 1-3 y Tabla 1-4 que el sistema experimental utiliza sensores de humedad
relativa, cabe aclarar entonces, que los niveles de humedad absoluta representado en Tabla 1-1
fueron evaluados bajo condiciones de humedad relativa; esta última es dependiente además de la
temperatura del aire; por consiguiente, en cada corrida experimental, para obtener un nivel de
humedad absoluta, se debía alcanzar un nivel de humedad relativa según el nivel de temperatura
del aire; estos niveles de humedad relativa se presentan en la Tabla 1-2. Por ejemplo, si se quería
trabajar con un nivel de alta humedad absoluta y un nivel de baja temperatura, la humedad relativa
del diseño debía ser 75% de humedad relativa.
d. Otras consideraciones del diseño de experimento

Otras consideraciones que se tuvieron en cuenta en el diseño de experimento se presentan a
continuación:
 Se observó que, al disminuir la temperatura del aire a la entrada del evaporador por debajo de
23°C, el aire a la salida del evaporador salía entre 8-10°C y el inverter actuaba para aumentar
esta temperatura. Esto causaba una condición extrema de operación difícil de alcanzar un
estado estable de operación y los resultados obtenidos no eran fiables.
 Los niveles de temperatura, y humedad relativa, a la entrada son condiciones muy difíciles de
mantener en estado permanente durante más de una hora de medición, por consiguiente,
durante cada experimentación, las variables se desplazaba en un rango el cual se puede
observar en la Tabla 1-1.
e. Factor Bloque
La ejecución de una corrida experimental requiere de mucho tiempo de preparación, ejecución
y post ejecución; es por esto que se consideró analizar un factor de bloque en la
experimentación. Este factor de bloque es el día de la experimentación; se realizaban 2 corridas
experimentales por día.
f. Descripción del sistema.

El esquema de este sistema experimental se presenta en la Figura 1-3, el control del equipo es a
través de 2 métodos:
1) Control de temperatura ajustada por el usuario: La regulación de temperatura por el usuario,
se realiza a través del modo “Set Temp” que presenta en el control remoto de la unidad; se
indica la temperatura deseada a la entrada del aire al evaporador, y el equipo opera y toma
medición de la temperatura hasta alcanzar el punto deseado; conforme se acerca a esta
temperatura, el inverter regula el flujo másico del refrigerante a través del voltaje y corriente
enviada al compresor.
2) Control de seguridad ajustada por el inverter. Este segundo control de temperatura no puede
ser controlado por el usuario, es programado por el fabricante en el inverter; consiste en un
sensor de temperatura (diferente al usado en el de control del usuario) que detecta la
temperatura del aire a la salida del evaporador, si la temperatura del aire se encuentra por
debajo de la temperatura de control del inverter, este último pasa a controlar la regulación del
flujo másico del refrigerante disminuyendo la caída de presión en el compresor, pero
aumentando su flujo, esto causa que la temperatura de salida del aire se eleve, alejando la
operación de la temperatura de control; este control opera con el objetivo de evitar problemas
de funcionamiento dentro del equipo como congelamiento de las tuberías, líquido en el
evaporador y/o compresor.
II. Resultados y validación del diseño de experimentos

Identificado el diseño de experimentos como un diseño factorial 2 3 con 4 bloques, y el
procedimiento para la ejecución del sistema experimental; se plantea el orden de las corridas
experimentales con ayuda del software STATGRAPHICS. La ejecución del experimento se realizó
como se establece en la norma AHRI para pruebas de equipos unitarios en el apartado 3.3 “Test
procedure for steady-state wet coil” [54]. Debido a lo extenso de las corridas experimentales; y,
teniendo en cuenta este tipo de experimentación fue una fase inicial, este diseño de experimento
fue realizado sin réplicas a expensa de asumir las consecuencias de tomar esta decisión. En total
se realizaron 8 corridas experimentales, y teniendo en cuenta los 4 bloques (día de
experimentación), las interacciones de 2 factores (AB, AC, BC) los grados de libertad del error es
1.
a. Análisis de resultados para el calor absorbido en el evaporador

Para el calor absorbido en el evaporador; como se observa en la Figura II-1; se puede identificar
que las interacciones ABC y BC siguen la línea de distribución normal, con media cero como se
indica en [51] es decir que sus efectos son insignificantes para el calor absorbido en el evaporador
y por lo tanto estos factores son excluidos del análisis.
Gráfica Mitad-Normal para Q_Evap
2
C:Nivel_velocidad
Desviaciones estándar
1,6
1,2 A:Temperatura_aire_in
B:Humedad_aire_in
0,8
ABC
BC+bloque
0,4
AB+bloque
AC+bloque
0
0 0,03 0,06 0,09 0,12 0,15
Efectos estandarizados
Figura II-1 Grafica de probabilidad normal de los efectos estandarizados para el calor absorbido por el evaporador.
La Tabla 1-7 presenta la tabla ANOVA, correspondiente a la significancia entre los factores
principales y sus interacciones respecto al calor cedido por el aire en el evaporador. Como se puede
observar, en la tabla el estadístico “valor-p” de los factores e interacciones se encuentra muy
alejado del umbral de aceptación para este diseño de experimento, el cual se estableció en un 95%
de confianza (p=0,05). Esto quiere decir, que ninguno de los factores principales, ni sus
interacciones son significativas para esta variable y por consiguiente no se continúa con el análisis.
b. Análisis de resultados para el trabajo eléctrico del compresor

Se analizan los resultados del diseño de experimentos para el trabajo del compresor; la Figura II-2
presenta la gráfica normal de los efectos estandarizados y se observa que las interacciones
(Temperatura-Humedad; Temperatura-Nivel de velocidad; Humedad-Nivel de velocidad), los
cuales están confundidos con los factor bloques, no son normalmente distribuidos y por esto, son
considerados como significativos; el resto de factores son excluidos del análisis.
La Tabla 1-8 presenta el ANOVA obtenida al descartar los efectos no significativos del
experimento. Se observa que las interacciones confundidas con los bloques son significativas para
esta variable de respuesta. Los estadísticos R2 (93,7656 %) indican que el modelo utilizado en el
análisis de varianza explica en gran parte la variabilidad del trabajo eléctrico del compresor.
Teniendo en cuenta que es un modelo con más de una variable independiente el estadístico a
analizar debe ser el R2ajustada (89,0898 %); este valor indica que el factor bloque y las interacciones
Temperatura-Humedad del aire, Temperatura del aire-Nivel de velocidad del ventilador y
Humedad del aire-Nivel velocidad del aire describen en gran proporción la variabilidad del
compresor
La validación de los supuestos del experimento se presenta la Figura II-3, se observan en cada uno
de los gráficos que se cumplen los supuestos de normalidad (comportamiento lineal ajustado a una
tendencia normal), independencia (no presentan tendencia alguna entre las corridas) y
homocedasticidad.
Gráfica Mitad-Normal para W_comp
2
AB+bloque
1,6
1,2 BC+bloque
AC+bloque
0,8
C:Nivel_velocidad
ABC
0,4
B:Humedad_aire_in
A:Temperatura_aire_in
0
0 0,02 0,04 0,06 0,08 0,1
Figura II-2 Grafica de probabilidad normal de los efectos estandarizados para el consumo eléctrico del compresor.
Para este último, a pesar de tener un patrón un poco irregular en la gráfica de residuos vs predichos,
el orden de magnitud en los residuos es muy pequeño, esto indica que la varianza es constante para
considerar heterocedasticidad en el diseño de experimento.
Gráfica de Probabilidad Normal para Residuos Gráfica de Residuos para W_comp
(X 0,001)
99,9 30
99
20
95
10
80
porcentaje
residuo
50 0
20 -10
5
-20
1
0,1 -30
-27 -17 -7 3 13 23 33 0 2 4 6 8
(X 0,001) número de ejecución
residuos
Gráfica de Residuos para W_comp
(X 0,001)
30
20
10
residuo
-10
-20
-30
0,46 0,5 0,54 0,58 0,62 0,66 0,7
predichos
Figura II-3 Gráfica de probabilidad normal para validar normalidad (superior-izquierda), Gráfica de residuos vs
orden de ejecución para validación de independencia (superior-derecha) y Gráfica de residuos vs predicho para
validación de homocedasticidad (inferior). Fuente: elaboración propia
c. Análisis de resultados para el coeficiente de desempeño del equipo
El último estudio que se presenta es el realizado al coeficiente de desempeño de la unidad. La
Figura II-4 presenta una distribución normal casi en todos los factores, excepto la interacción AB,
y el nivel de velocidad del ventilador. A su vez, debido que las interacciones AC y BC están
confundidas con el bloque, estos no deben salir del experimento, ya que el factor bloque aparecer
en la tabla ANOVA
Gráfica Mitad-Normal para COP
2
AB+bloque
1,6
1,2 BC+bloque
AC+bloque
0,8
C:Nivel_velocidad
ABC
0,4
A:Temperatura_aire_in
B:Humedad_aire_in
0
0 0,3 0,6 0,9 1,2 1,5
Figura II-4 Grafica de probabilidad normal de los efectos estandarizados para el coeficiente de desempeño. Fuente:
elaboración propia
La Tabla 1-9 presenta el ANOVA para el COP; los estadísticos para esta tabla es: R2 (96,9281 %)
y R2ajustada (92,8322 %) por lo tanto, el modelo y las variables independientes del ANOVA
describen en gran medida la variabilidad del COP. Por último, la Figura II-5 presenta la validez
gráfica de los supuestos del diseño del experimento, se observa un ajuste adecuado en la
normalidad (figura superior izquierda), es decir, comportamiento lineal se ajusta a una tendencia
normal, la homocedasticidad (figura superior derecha) se cumple ya que no hay aumento de los
residuos a lo largo de lo predichos ni tampoco tendencias definidas. Por ultimo en la gráfica de
independencia (figura inferior) no presentan tendencia alguna entre el orden de las corridas.
Gráfica de Probabilidad Normal para Residuos Gráfica de Residuos para CO P
99,9 0,25
99
0,15
95
80
porcentaje
0,05
residuo
50
-0,05
20
5
-0,15
1
0,1 -0,25
-0,25 -0,15 -0,05 0,05 0,15 0,25 5 6 7 8 9
residuos predichos
Gráfica de Residuos para CO P
0,25
0,15
0,05
residuo
-0,05
-0,15
-0,25
0 2 4 6 8
número de ejecución
Figura II-5 Gráfica de probabilidad normal para validar normalidad (superior-izquierda), Gráfica de residuos vs
predicho para validación de homocedasticidad (superior-derecha) y Gráfica de residuos vs orden de ejecución para
validación de independencia (inferior). Fuente: elaboración propia
III. Prueba de idoneidad de los modelos de regresión
En esta sección se presentarán el análisis de idoneidad obtenidas en cada modelo de regresión para
ambos factores de respuestas, realizado a partir del paquete “Análisis de Regresión de Excel”. Se
presenta en el modelo 1 la tabla ANOVA, la prueba de coeficientes y la gráfica de residuales,
finalmente se incluirá una observación para indicar porqué se rechaza este modelo.
a. Modelo 1
ANOVA
ANÁLISIS DE VARIANZA jHT ANÁLISIS DE VARIANZA jMT
GdL SC CM F Valor P GdL SC CM F Valor P

Regresión 1 462974,591 462974,591 134,312538 2,4837E-05 1 389770,759 389770,759 80,1972946 0,00010827
Residuos 7 24128,9621 3446,99459 7 34021,0394 4860,14848
Total 8 487103,553 8 423791,799
Coeficientes
Prueba de coeficientes j MT
Prueba de coeficientes jHT
Coeficientes Probabilidad Coeficientes Probabilidad
Intercepción 0 #N/A 0 #N/A
Re_Dc 0,18484575 8,0295E-06 0,16960378 4,4061E-05
Gráfica de residuos (izquierda jHT, derecha jMT)

Re_Dc Gráfico de los residuales
150
Re_Dc Gráfico de los residuales
150
100
100
50
50
Re siduos
Re siduos
0
0
900,00 1000,00 1100,00 1200,00 1300,00 1400,00 1500,00 1600,00 1700,00 1800,00 1900,00 900,00 1000,00 1100,00 1200,00 1300,00 1400,00 1500,00 1600,00 1700,00 1800,00 1900,00
-50
-50
-100
-150
-100 Re_Dc
Re_Dc
Se rechaza porque existe una estructura lineal en el residuo, es un ajuste lineal mal calculado
Considerando que todos los modelos tienen un valor p =0.05 para el ANOVA y los coeficientes,
solo se presentará las gráficas de los residuales para los siguientes modelos ya que estos son los
que rechazan el modelo. La gráfica a la izquierda representa el resultado de los residuales para el
modelo de jHT y el de la derecha para el modelo jMT.
b. Modelo 2
En este caso, no existe una estructura definida, sin embargo, el modelo 10 presenta un gráfico de
residual más aleatorio y por esto no es escogido este modelo.
c. Modelo 3
Este modelo es rechazado porque la estructura sugiere que existe una relación no lineal entre las
variables, y no la escogida en este modelo.
d. Modelo 4
e. Modelo 5
variables, y existe una aparente heterocedasticidad
f. Modelo 6
variables, y no la escogida en este modelo.
g. Modelo 7
h. Modelo 8
variables, y no la escogida en este modelo
i. Modelo 9
Bibliografía
[1] T. Harford, “How air conditioning changed the world - BBC News.” [Online]. Available:
https://www.bbc.com/news/business-39735802. [Accessed: 08-Feb-2019].
[2] I. International Energy Agency, “The Future of Cooling Opportunities for energy-efficient
air conditioning Together Secure Sustainable,” 2018.
[3] JRAIA, “World Air Conditioner Demand by Region The Japan Refrigeration and Air
Conditioning Industry Association,” Tokyo, 2018.
[4] DATA VAN DE DAG, “World population distribution by latitude and longitude - 2015 |
Datagraver,” 2016. [Online]. Available: https://www.datagraver.com/case/world-
population-distribution-by-latitude-and-longitude-2015. [Accessed: 11-Apr-2019].
[5] “World Energy Trilemma Index | 2016 BENCHMARKING THE SUSTAINABILITY OF
NATIONAL ENERGY SYSTEMS In Partnership with OLIVER WYMAN,” London -
England, 2016.
[6] ACAIRE, “Marzo 2016 3,” Comunicado ACAIRE, 2016. [Online]. Available:
http://spiroductos.com/2016/03/. [Accessed: 16-Apr-2019].
[7] K. J. Chua, S. K. Chou, W. M. Yang, and J. Yan, “Achieving better energy-efficient air
conditioning – A review of technologies and strategies,” Appl. Energy, vol. 104, pp. 87–
104, Apr. 2013.
[8] A. F. Handbook, “American society of heating, refrigerating and air-conditioning
engineers,” Inc. Atlanta, GA, USA, 2009.
[9] Z. Li and S. Deng, “A DDC-based capacity controller of a direct expansion (DX) air
conditioning (A/C) unit for simultaneous indoor air temperature and humidity control -
Part I: Control algorithms and preliminary controllability tests,” Int. J. Refrig., vol. 30, no.
1, pp. 113–123, 2007.
[10] J. H. Huh and M. J. Brandemuehl, “Optimization of air-conditioning system operating
strategies for hot and humid climates,” Energy Build., vol. 40, no. 7, pp. 1202–1213, 2008.
[11] A. Saeed and A. Al-Alili, “A review on desiccant coated heat exchangers,” Sci. Technol.
Built Environ., vol. 23, no. 1, pp. 136–150, 2017.
[12] C. X. Jia, Y. J. Dai, J. Y. Wu, and R. Z. Wang, “Analysis on a hybrid desiccant air-
conditioning system,” Appl. Therm. Eng., vol. 26, no. 17–18, pp. 2393–2400, Dec. 2006.
[13] ASHRAE, Handbook HVAC Systems and Equipment. 2012.
[14] Z. Li and S. Deng, “A DDC-based capacity controller of a direct expansion (DX) air
conditioning (A/C) unit for simultaneous indoor air temperature and humidity control -
Part II: Further development of the controller to improve control sensitivity,” Int. J.
Refrig., vol. 30, no. 1, pp. 124–133, 2007.
[15] X. Xu, Z. Zhong, S. Deng, and X. Zhang, “A review on temperature and humidity control
methods focusing on air-conditioning equipment and control algorithms applied in small-
to-medium-sized buildings,” Energy Build., vol. 162, pp. 163–176, 2018.
[16] Q. Qi and S. Deng, “Multivariable control of indoor air temperature and humidity in a
direct expansion (DX) air conditioning (A/C) system,” Build. Environ., vol. 44, no. 8, pp.
1659–1667, 2009.
[17] W. M. Kays and A. L. London, “Compact Heat Exchangers (Part III),” in Compact Heat
Exchangers , 1984.
[18] F. C. McQuiston and J. D. Parker, “Heating, ventilating, and air conditioning: analysis and
design,” 1982.
[19] T. H. Kuehn, J. W. (viaf)3764742 Ramsey, and J. L. T. environmental engineering.
Threlkeld, Thermal environmental engineering, 3rd ed. Upper Saddle River, 1998.
[20] N. A. A. Qasem and S. M. Zubair, “Compact and microchannel heat exchangers: A
comprehensive review of air-side friction factor and heat transfer correlations,” Energy
Conversion and Management. 2018.
[21] C. Wang, Y. Hsieh, and Y. Lin, “Performance of Plate Finned Tube Heat Exchangers
Under Dehumidifying Conditions,” J. Heat Transfer, 1997.
[22] C. C. Wang, W. S. Lee, and W. J. Sheu, “A new approach to correlate the frictional
performance of fin-and-tube heat exchangers in wet conditions,” Heat Transf. Eng., vol.
23, no. 2, pp. 15–21, 2002.
[23] C. C. Wang and K. Y. Chi, “Heat transfer and friction characteristics of plain fin-and-tube
heat exchangers, part I: New experimental data,” Int. J. Heat Mass Transf., vol. 43, no. 15,
pp. 2681–2691, Aug. 2000.
[24] M. H. Kim and C. W. Bullard, “Air-side performance of brazed aluminum heat
exchangers under dehumidifying conditions,” in International Journal of Refrigeration,
2002.
[25] N. H. Kim and T. Kim, “An experimental investigation on the airside performance of fin-
and-tube heat exchangers having slit fins under wet condition,” J. Mech. Sci. Technol.,
vol. 29, no. 11, pp. 5011–5019, 2015.
[26] Y. Zou, H. Li, K. Tang, and P. Hrnjak, “Round-Tube and Microchannel Heat Exchanger
Modeling at Wet Air Condition,” in International Refrigeration and Air Conditioning
Conference, 2016, pp. 1–10.
[27] Y. Zou, H. Li, and P. Hrnjak, “Simulation of a R410A Residential Air-Conditioning
System with Round-Tube and / or Microchannel Evaporators and Condensers under both
Dry and Wet Air Conditions,” Int. Compress. Eng. Refrig. Air Cond. High Perform. Build.
Conf., pp. 1–10, 2016.
[28] G. Zhang, B. Wang, X. Li, W. Shi, and Y. Cao, “Review of experimentation and modeling
of heat and mass transfer performance of fin-and-tube heat exchangers with
dehumidification,” Appl. Therm. Eng., vol. 146, pp. 701–717, Jan. 2019.
[29] B. P. Rasmussen, “Dynamic modeling for vapor compression systems-Part I: Literature
review,” HVAC R Res., vol. 18, no. 5, pp. 934–955, 2012.
[30] B. P. Rasmussen, C. Price, J. Koeln, B. Keating, and A. Alleyne, “HVAC system
modeling and control: Vapor compression system modeling and control,” in Advances in
Industrial Control, 2018.
[31] V. M. SÁNCHEZ, “Desarrollo de un Modelo Dinámico Reducido de un Sistema de
Refrigeración Industrial basado en el Ciclo de Compresión de Vapor,” 2012.
[32] B. P. Rasmussen and A. G. Alleyne, “Control-Oriented Modeling of Transcritical Vapor
Compression Systems,” J. Dyn. Syst. Meas. Control, vol. 126, no. 1, p. 54, 2004.
[33] Y. Yao, W. Wang, and M. Huang, “A state-space dynamic model for vapor compression
refrigeration system based on moving-boundary formulation,” Int. J. Refrig., vol. 60, pp.
1–16, 2015.
[34] D. Shiming, “A dynamic mathematical model of a direct expansion ( DX ),” Build.
Environ., vol. 35, pp. 603–613, 2000.
[35] Q. Qi and S. Deng, “Multivariable control-oriented modeling of a direct expansion (DX)
air conditioning (A/C) system,” Int. J. Refrig., vol. 31, no. 5, pp. 841–849, 2008.
[36] B. P. Rasmussen and B. Shenoy, “Dynamic modeling for vapor compression systems-Part
II: Simulation tutorial,” HVAC and R Research. 2012.
[37] Y. H. Yau, “Application of a heat pipe heat exchanger to dehumidification enhancement
in a HVAC system for tropical climates-a baseline performance characteristics study,” Int.
J. Therm. Sci., vol. 46, no. 2, pp. 164–171, 2007.
[38] Y. Yao, K. Yang, M. Huang, and L. Wang, “A state-space model for dynamic response of
indoor air temperature and humidity,” Build. Environ., vol. 64, pp. 26–37, 2013.
[39] M. Saraireh, “Heat Transfer and Condensation of Water Vapour From Humid Air in
Compact Heat Exchangers,” no. February, p. 151, 2012.
[40] X. Xiangguo, D. Shiming, H. Xiaohong, and Z. Xuejun, “A control-oriented semi-
physical model for a direct expansion air conditioning system,” Build. Serv. Eng. Res.
Technol., vol. 35, no. 6, pp. 585–599, 2014.
[41] H. Pangborn and A. G. Alleyne, “Dynamic Modeling of Heat Exchangers With Humidity
and Condensation,” in ASME 2015 Dynamic Systems and Control Conference, 2016, p.
V003T39A001.
[42] S. Y. Liang, M. Liu, T. N. Wong, and G. K. Nathan, “Analytical study of evaporator coil
in humid environment,” Appl. Therm. Eng., vol. 19, no. 11, pp. 1129–1145, 1999.
[43] Z. Li, X. Xu, S. Deng, and D. Pan, “A novel proportional-derivative (PD) law based fuzzy
logic principles assisted controller for simultaneously controlling indoor temperature and
humidity using a direct expansion (DX) air conditioning (A/C) system,” Int. J. Refrig.,
vol. 57, pp. 239–256, 2015.
[44] M. Kumar, I. N. Kar, and A. Ray, “State space based modeling and performance
evaluation of an Air-Conditioning system,” HVAC R Res., vol. 14, no. 5, pp. 797–816,
2008.
[45] A. J. Bula Silvera and A. F. Méndez Díaz, “Ahorro de energía eléctrica por
implementación de filtros desecantes en sistemas de acondicionamiento de aire por
compresión de vapor,” Rev. Fac. Ing., vol. 31, pp. 84–92, 2004.
[46] J. G. Fajardo Cuadro, B. Sarria López, and M. Álvarez Guerra Plasencia, “Estudio
exergético del espacio climatizado en una embarcación fluvial.,” Ing. Eléctrica, vol. 36,
no. 2, pp. 127–135, 2015.
[47] M. B. Universidad de Pamplona (Colò mbia), F. G. Reina, J. J. C. Eras, and V. S. Santos,
Revista colombiana de tecnologías de avanzada., vol. 2, no. 32. Universidad de
Pamplona, 2018.
[48] M. Sanchez, “El clima tropical.” [Online]. Available:
https://www.meteorologiaenred.com/el-clima-tropical.html. [Accessed: 03-Aug-2019].
[49] “Clasificación climática de Köppen - Meteo Navarra.” [Online]. Available:
http://meteo.navarra.es/definiciones/koppen.cfm. [Accessed: 03-Aug-2019].
[50] IDEAM, “Carácterísticas Climatológicas De Ciudades Principales Y Municipios
Turísticos,” 2015.
[51] D. C. Montgomery, “Design and analysis of experiments 5ed,” Quality and Reliability
Engineering International. 2001.
[52] M. Tanco, E. Viles, L. Ilzarbe, and M. J. Alvarez, “Implementation of Design of
Experiments projects in industry,” Appl. Stoch. Model. Bus. Ind., vol. 25, no. 4, pp. 478–
505, Jul. 2009.
[53] International Standard, “Non-ducted air conditioners and heat pumps — Testing and
rating for performance,” Iso 5151, vol. 2010, 2010.
[54] ANSI/AHRI, “Standard, Performance Rating of Unitary Air-Conditioning & Air-Source
Heat Pump Equipment,” ANSI/AHRI Stand. 210/240, vol. 2, no. December, 2012.
[55] V. P. Carey, Liquid-Vapor Phase-Change Phenomena. CRC Press, 2018.
[56] S. Wongwises and Y. Chokeman, “Effect of fin pitch and number of tube rows on the air
side performance of herringbone wavy fin and tube heat exchangers,” Energy Convers.
Manag., 2005.
[57] L. Xia, M. Y. Chan, S. M. Deng, and X. G. Xu, “A modified logarithmic mean enthalpy
difference (LMED) method for evaluating the total heat transfer rate of a wet cooling coil
under both unit and non-unit Lewis Factors,” Int. J. Therm. Sci., vol. 48, no. 11, pp. 2159–
2164, Nov. 2009.
[58] J. P. Holman, Experimental methods for engineers. 2012.
[59] G. L. Wedekind, B. L. Bhatt, and B. T. Beck, “A system mean void fraction model for
predicting various transient phenomena associated with two-phase evaporating and
condensing flows,” Int. J. Multiph. Flow, 1978.
[60] E. Data and B. Iii, “Chapter 17 Void Fractions in Two-Phase Flows,” Wolveri ne Tube
Inc, 2005.
[61] J. M. Corberán, P. Fernández de Córdoba, S. Ortuno, V. Ferri, and P. Montes, “Detailed
modelling of evaporators and condensers,” Proc. 8th Int. Refrig. Air Cond. Conf., 2000.
[62] T. O. F. Contents, “MATH2071 : LAB # 5 : Implicit ODE Methods,” no. x, pp. 1–8.

Tesis

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Tesis

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis

Cargado por

Copyright:

Formatos disponibles

Análisis del desempeño energético y desarrollo del coeficiente de

transferencia de calor y masa en un sistema de acondicionamiento de aire

Tesis para optar al título de doctor en ingeniería mecánica

Hermes Javier Ramírez León

Programa de Doctorado en ingeniería mecánica

Realizada por: Hermes Javier Ramírez León

Presidente: Mauricio Yilmer Carmona García, PhD

Programa de Doctorado en ingeniería mecánica

A Dios, gracias por tu bendición y fortaleza

Fuente: elaboración propia

1.1 Características del aire en el clima tropical

1.2 Diseño de Experimentos

hair  1.006  Taire  waire  2501  1.86  Taire  (1-5)

WComp _ ref  m ref href ,out  href ,in Compresor (1-7)

Fuente: elaboración propia

Humedad relativa según temperatura y humedad absoluta

Fuente: elaboración propia

1.3 Descripción del sistema.

Para desarrollar esta investigación, se construyó un sistema experimental. En la Figura 1-2 se

Tabla 1-3 Descripción de los sensores instalados en el sistemas experimental

Características de los sensores

Fuente: elaboración propia

Ubicación de los sensores

Fuente: elaboración propia

1.4 Puesta en marcha y Resultados

Orden aleatorio de experimentación

Fuente: elaboración propia

U min  1.070 m

Con el flujo másico y la diferencia de entalpía se calculará el calor absorbido en el evaporador:

Tref _ out  T 7  24.7C pref _ out  P3  240.80 psig f  F1  1091.04 mL min

Salida refrigerante evaporador

Tref _ out  T1  10.74C pref _ out  P1  79.71psig

Salida refrigerante compresor

Tref _ out  T 2  55.12C pref _ out  P2  262.57 psig

T 7  REFPROPT7, P3  1061.28 kg

href _ out _ Evap  REFPROPT1, P1  436.47 kJ kg

href _ out _ Comp  REFPROPT2, P2  458.530 kJ kg

Con las propiedades calculadas se procede a determinar QEvap_ref; Wcomp_ref, y COP_ref

Orden de experimentación y variables de respuestas

Fuente: elaboración propia

1.5 Metodología de análisis de resultados

1.6 Análisis de resultados para el calor absorbido en el evaporador

Análisis de Varianza para Q_Evap

Fuente: elaboración propia

1.7 Análisis de resultados para el trabajo eléctrico del compresor

Análisis de Varianza para WComp

Fuente: elaboración propia

Diagrama de Pareto Estandarizada para W_comp

Gráfica de Interacción para W_comp

1.8 Análisis de resultados para el coeficiente de desempeño del equipo

Análisis de Varianza para COP

Fuente: Elaboración propia

Superficie de Respuesta Estimada Superficie de Respuesta Estimada

1.9 Conclusiones del capítulo

Además de la transferencia de calor, la deshumidificación también se verá afectada por el aire de

2.1 Descripción matemática

qAire  qSensible  qLatente (2-2)

qSensible  CpaireTa  Tpared  (2-3)

qAire  CPaireTa  Tpared   hfg wa  wpared  (2-6)