INTRODUCCIÓN
Trabajos previos han reportado el diseño y la construcción de dos módulos experimentales de vivienda (MEV), de dimensiones 2,94 m x 3,95 m (medidas interiores), con el análisis de diferentes técnicas constructivas de aislamiento térmico, incluyendo técnicas de climatización pasiva [1], [2]. Los MEV están construidos en el Perú, en el centro poblado de Imata, en el distrito de San Antonio de Chuca, de la provincia Caylloma, región Arequipa, a 4500 msnm con una latitud de 15° 50’ Sur y longitud de 71° 51’ Oeste, donde la temperatura mínima registrada para el año 2019 fue -14,4 °C [3], con una energía solar incidente diaria promedio anual de 5,5-6,0 kWh/m2 [4].
En el presente trabajo se evalúa la incorporación de un sistema de calefacción solar activo (SCSA) al MEV sin invernadero adosado, que se muestra en la Fig. 1, que provee la energía térmica necesaria para alcanzar una “temperatura neutra” (también denominada temperatura de confort), que brinde salud a los ocupantes del MEV.
METODOLOGÍA
Primero, se presenta un análisis de verificación térmica, consistente en comparar las predicciones de temperatura del modelo desarrollado con el programa EnergyPlus con los registros de temperatura medidos en el MEV, para asegurar que se implemente correctamente un modelo de transferencia de calor en EnergyPlus. Luego se calcula la temperatura neutra que se desea mantener al interior del MEV. Por último, se calcula, con EnergyPlus, la energía requerida para mantener la temperatura neutra al interior del MEV.
SIMULACIÓN CON ENERGYPLUS
El Departamento de Energía de los Estados Unidos (Department of Energy, DOE) ha desarrollado y publicado un programa relativamente nuevo de simulación energética de edificios completos, EnergyPlus [5]. El programa fusionó los programas BLAST (Building Loads Analysis and System Thermodynamics) y DOE-2 para utilizar sus mejores capacidades y características. EnergyPlus está diseñado para facilitar el desarrollo de módulos e interfaces de usuario.
Para ejecutar la simulación es necesario conocer de antemano la condición medioambiental y la situación del estado del MEV, tales como:
El diseño arquitectónico; conocer distribuciones y especificaciones geométricas de la construcción (pisos, paredes, techo, ventanas y puertas).
Las zonas térmicas; modelado en 3D del MEV con el uso del software SketchUp para el presente trabajo, de acuerdo a la cantidad de ambientes que se desean analizar. Una zona térmica representa un volumen de aire a temperatura uniforme la cual se delimita por superficies del modelado. El módulo se ha representado con cuatro zonas térmicas, como se ve en la Fig. 2: Zona Ático Este (ZAE), Zona Ático Claraboya (ZAC), Zona Ático Oeste (ZAO) y Zona Interior (ZI).
Ubicación y datos meteorológicos; se requiere conocer la ubicación geográfica del lugar para identificar posibles sombras, como se puede ver en la Fig. 3, y datos meteorológicos de la zona (Se usaron datos de SENAMHI y del programa SOLCAST).
Materiales constructivos y sus propiedades; conocer los diferentes materiales utilizados en la construcción de la edificación, sus propiedades termofísicas, como densidad, conductividad térmica y calor específico, y propiedades ópticas como, reflectancia, transmitancia, y absortancia. Estos datos están indicados en las Tab. I y Tab. II.
Variables de salida; definir parámetros que serán resultado de la simulación, por ejemplo, temperatura y humedad relativa al interior de las zonas térmicas, ganancia y pérdida de energía por convección en superficies interior y exterior, radiación solar incidente sobre superficie exterior e interior, entre otras.
Las superficies en morado de la Fig. 2 y Fig. 3b, restringen la incidencia solar directa con la proyección de sus sombras en el transcurso del sol para la ubicación del MEV, indicado en color rojo. Por ello es importante también establecer la orientación adecuada. En la Fig. 2 se aprecia la orientación del módulo: N21°30’W.
Nombre | Espesor (m) | Conductividad térmica (W/m-K) | Densidad (kg/m3) | Calor específico (J/kg-K) | Absortancia térmica | Absortancia solar | Absortancia visible |
Tierra1 | 0,25 | 0,52 | 2050 | 1840 | 0,88 | 0,7 | 0,7 |
Roca1 | 0,15 | 3,5 | 2750 | 880 | 0,9 | 0,7 | 0,7 |
Madera1 | 0,02 | 0,12 | 510 | 1380 | 0,9 | 0,6 | 0,7 |
EPS1 | 0,05 | 0,033 | 30 | 1700 | 0,9 | 0,7 | 0,7 |
Triplay1 | 0,004 | 0,14 | 560 | 1400 | 0,9 | 0,6 | 0,7 |
Calamina1 | 0,0005 | 110 | 7837 | 500 | 0,03 | 0,09 | 0,7 |
Adobe2 | 0,4 | 0,5 | 2000 | 950 | 0,92 | 0,65 | 0,7 |
Yeso1 | 0,02 | 0,4 | 900 | 1000 | 0,9 | 0,5 | 0,7 |
Fibrocemento1 | 0,004 | 0,18 | 1500 | 1550 | 0,9 | 0,7 | 0,7 |
Madera1 | 0,04 | 0,12 | 510 | 1380 | 0,9 | 0,6 | 0,7 |
Nota: 1 Valores obtenidos de la norma técnica peruana EM110. [7]. 2 Valor obtenido de [8].
Nombre | U-Factor (W/m2-K) | Coeficiente de Ganancias de Calor Solar | Transmitancia Visible |
Policarbonato Alveolar Transparente | 3,8 | 0,75 | 0,8 |
Nota: Valores obtenidos del catálogo GLANZE [9].
En la simulación se usaron las propiedades físicas y ópticas por defecto dadas en el EnergyPlus para la ventana (espesor: 3 mm, conductividad térmica: 0,9 W/m-K, área: 0,7 m2). En la Tab. II se puede ver las propiedades ópticas y físicas de la claraboya usadas en la simulación.
A cada zona térmica se le aplica la configuración de una infiltración de aire por hora. Esto consiste en renovar todo el aire al interior de cada zona por el aire del exterior.
Como se puede ver en la Fig. 3a las zonas son adyacentes entre sí. Se plantea la configuración de un intercambio de aire por hora entre zonas: ZI con ZAC, ZI con ZAO, ZI con ZAE y 0,1 intercambios de aire entre las zonas: ZAO con ZAC y ZAE con ZAC. Este intercambio de aire consiste en intercambiar en igual cantidad de aire cada dos zonas, lo que afecta el balance de energía para ambas zonas.
COMPARACIÓN DE TEMPERATURA
Los resultados de la simulación, con las propiedades físicas y condiciones de modelado ingresadas (ver Tab. I, Tab II, Fig 2 y Fig. 3), se comparan con las temperaturas del aire interior horaria registrados en el MEV, para las fechas 25/05/2019 al 11/06/2019 y del 14/09/2019 al 18/10/2019. El equipo usado para el registro fue un Data Logger HOBO U12-012.
Estos resultados indican que la predicción de la simulación, bajo los parámetros ingresados, demuestran una adecuada correlación entre la data experimental y la simulada.
Esta comparación de temperaturas permite extender la predicción para los demás periodos del año 2019 bajo los mismos parámetros ingresados.
TEMPERATURA NEUTRA
En 1975, Humphreys [10] encontró por estudios de campo la temperatura con un mínimo estrés térmico o temperatura neutra (Tn). Esta temperatura fue reportada en niveles medidos de temperatura de globo o temperatura del aire interior (𝑇𝑖). Las temperaturas experimentadas por los encuestados fueron entre 17 °C y 30 °C:
Posteriormente en 1976 Humphreys [11], sustituye la temperatura del aire interior por la temperatura media mensual exterior (𝑇𝑚) con resultados similares:
Según la data meteorológica de la estación IMATA [4], la 𝑇𝑚 para el año 2019 fue de 3,4 °C. Extrapolando la Ec. 2 de Humphreys se tiene que al reemplazar 𝑇𝑚 en la ecuación, se obtiene para Imata 𝑇𝑛 = 13,7 °C.
CÁLCULO DE CARGA DE CALEFACCIÓN
Las cargas de calefacción fueron calculadas de manera dinámica por el EnergyPlus. El promedio diario mensual de estas cargas dinámicas da la energía requerida en promedio diario para elevar la temperatura desde sus condiciones base. Una vez propuesta la temperatura neutra de 13,8 °C, según el modelo adaptativo de Humphreys, se calcula mediante simulación la energía requerida para mantener esta temperatura neutra al interior del MEV. La predicción de la simulación con los parámetros ingresados, indica que durante un año (los datos corresponden a 2019) se necesita 3,54 MWh, correspondiente, en promedio anual, a una demanda diaria de 9,7 kWh.
Como se observa en la Fig. 6, los requerimientos de energía varían según las condiciones climáticas, resultando un requerimiento diario (promedios mensuales) máximo de 11 kWh en Julio y un mínimo de 8,4 kWh por día en diciembre.
Cabe indicar que se tiene un alto requerimiento de energía, resultado de las grandes pérdidas de calor por la puerta, de calamina metálica, y por la ventana, con vidrio simple y marco de fierro, tal como es tradición de uso en las viviendas en esta región. Se aprecian grandes variaciones de la temperatura, resultado de las pérdidas por la alta conductividad térmica (Ver Tab. I) de puerta de calamina-EPS y ventana de vidrio simple. Por lo tanto, primeramente, se debe colocar una puerta y una ventana de buen aislamiento térmico, antes de pensar en la instalación de un sistema de calefacción solar activo.
En la Fig. 7 se muestra las simulaciones de las temperaturas horarias en el MEV (curva verde), con los cambios de puerta de madera de 5 cm de espesor y ventana con doble vidrio para el mes de julio (mes más frío), y si se incluiría un Sistema de Calefacción Solar (SCS), con 2 m2 de colector solar sobre un plano inclinado 20° hacia el norte y una eficiencia total de 35%, típica para un SCS (curva celeste).
El SCS produce un incremento de 4,7 °C de la temperatura en la ZI del MEV para el mes de julio.
CONCLUSIONES
Con el programa EnergyPlus, se realizaron cálculos dinámicos del comportamiento térmico de una edificación, se pudo simular la temperatura en el interior de un módulo de vivienda en Imata (Arequipa), a 4500 msnm, en buena concordancia con las temperaturas medidas.
Esto permitió calcular la demanda de energía que se requiere para mantener en la ZI del módulo una temperatura neutra de 12 °C, para el mes de julio: la demanda energética diaria varía de 8,4 kWh en diciembre y 11 kWh en julio, el mes más frío.
Esta alta demanda de energía indica que se debe primeramente reducir las pérdidas de calor del módulo, para reducir el requerimiento energético.