SciELO - Scientific Electronic Library Online

 
vol.31 número1Propuesta Tecnológica para la Implementación de una Planta Fotovoltaica de 4.5 Kwp para la Comunidad Rural Altoandina San Francisco de Raymina, Ayacucho-PerúConductividad Térmica de Compuestos Tipo Sándwich Usados en la Industria de la Construcción índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

  • Não possue artigos citadosCitado por SciELO

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Tecnia

versão impressa ISSN 0375-7765versão On-line ISSN 2309-0413

Tecnia vol.31 no.1 Lima jan./jun. 2021

http://dx.doi.org/10.21754/tecnia.v21i1.1083 

Artículos originales

Simulación del Consumo Energético para la Calefacción de un Módulo de Vivienda Experimental Ubicado a 4500 Msnm

Simulation of Energy Consumption for the Heating of an Experimental Housing Module Located at 4500 Msnm

Pablo Ratache-Rojas1 
http://orcid.org/0000-0001-6648-2258

Cesar Rivasplata-Cabanillas2 
http://orcid.org/0000-0002-2129-1893

Manfred Horn-Mutschler1 
http://orcid.org/0000-0003-4808-658X

Mónica Gómez-León1 
http://orcid.org/0000-0003-0990-0593

1Universidad Nacional de Ingeniería, Facultad de Ciencias, Lima, Perú

2Empresa Solarsur EIRL, Tacna, Perú

RESUMEN

En el presente trabajo se simula la demanda energética en un módulo de vivienda experimental construido en el poblado de Imata, Perú, ubicado a 4500 msnm. Previo a esto se desarrollaron las mediciones de temperatura y los cálculos energéticos necesarios para validar el modelo de simulación con el programa EnergyPlus. La predicción de la simulación indica que se necesita en un año 3,54 MWh, con un máximo diario de 11 kWh en julio y un mínimo diario de 8,4 kWh en diciembre (promedios mensuales), para mantener una temperatura neutra de 12 °C y 15 °C respectivamente.

Palabras Clave: Simulación Energética; EnergyPlus; Temperatura Neutra; Demanda Energética Térmica.

ABSTRACT

In the present work, the energy demand is simulated for an experimental housing module built in the village of Imata, Peru, located at 4500 masl. Prior, the temperature measurements and the energy calculations necessary to validate the simulation model with the program EnergyPlus were performed. To maintain a neutral temperature of 12 °C in July and 15 °C in December, the simulation indicated that 3,54 MWh are needed yearly, with a daily maximum of 11 kWh in July and a daily minimum of 8,4 kWh in December (monthly averages).

Keywords: Energy Simulation; EnergyPlus; Neutral Temperature; Thermal Energy Demand

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.

Fig. 1 Módulos construidos en el centro poblado de Imata, para efectos comparativos [2]. 

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.

Fig. 2 Modelado de un módulo experimental de vivienda (MEV) sin invernadero adosado en 3D elaborado con SketchUp y ubicación de zonas térmicas. 

Fig. 3 a) Vista superior del MEV (Google, s.f.). Adaptado de Google Maps [6]. b) En morado se muestra edificaciones que proyectan sombra en el MEV. 

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.

Tabla 1 Valores de propiedades termo físicas de materiales usados en la simulación MEV. 

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].

Tabla 2 Propiedades ópticas de material usado para la claraboya en la simulación del MEV. 

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.

Fig. 4 Curva de temperatura horaria del aire al interior de la ZI del MEV, del 25/05/2019 al 10/06/2019. En azul la temperatura simulada por el programa EnergyPlus y en rojo la temperatura registrada por el sensor de temperatura. 

Fig. 5 Curva de temperatura horaria del aire al interior de la ZI del MEV, del 14/09/2019 al 7/10/2019. En azul la temperatura simulada por el programa EnergyPlus y en rojo la temperatura registrada por el sensor de temperatura. 

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.

Fig. 6 Comparación mensual de la energía solar promedio diario contra la energía requerida promedio diario para mantener la temperatura neutra en cada mes, calculada con el EnergyPlus. Se muestra en barras rojas la energía requerida en promedio diario para elevar la temperatura del aire al interior de la ZI para cada mes y en barra azul la energía solar promedio diario incidente en un metro cuadrado de plano inclinado 20° hacia el norte. 

Fig. 7 Temperaturas simuladas horaria del aire de la ZI del MEV para el mes de julio. En verde, las temperaturas simuladas bajo los parámetros ingresados para la verificación térmica. En celeste, las temperaturas simuladas con el funcionamiento de un SCS, suministrando diariamente 4,6 kWh de calor al MEV, en forma uniforme de 192 W, durante las 24 horas. 

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.

AGRADECIMIENTOS

Este trabajo se llevó a cabo bajo el auspicio del Consejo Nacional de Ciencia, Tecnología e Innovación del Perú (CONCYTEC), por medio del Convenio N° 024-2016-Fondecyt. Se agradece al Servicio Nacional de Meteorología e Hidrología del Perú (SENAMHI) y a la empresa SOLCAST por facilitar los datos meteorológicos necesarios para la creación del archivo de clima usado en la simulación.

REFERENCIAS

[1] J. Molina, G. Lefebvre, M. Horn, y M. Gómez, “Diseño de un módulo experimental bioclimático obtenido a partir del análisis de simulaciones térmicas para el centro poblado de Imata (4519 m s.n.m.) ubicado en Arequipa, Perú”, Inf. tecnol., vol. 31, no. 2, pp. 173-185, abr. 2020. [ Links ]

[2] J. Molina, M. Ponce, M. Horn y M. Gómez, “Towards a sustainable bioclimatic approach for the Peruvian high Andean rural area: Evaluation of the thermal contribution of a greenhouse attached to a dwelling”, Proc. ISES Solar World Congress 2019 IEA SHC International Conference on Solar Heating and Cooling for Buildings and Industry 2019, pp. 355-364, nov. 2019. [ Links ]

[3] SENAMHI (2003, Enero). ATLAS DE ENERGÍA SOLAR DEL PERÚ[Online]. Available: https://hdl.handle.net/20.500.12542/343. [ Links ]

[4] SENAMHI (2021, Febrero). Estaciones Automáticas [Online]. Available: https://web2.senamhi.gob.pe/?p=estaciones--automaticas. [ Links ]

[5] D. B. Crawley et al., “EnergyPlus: creating a new generation building energy simulation program”, Energy Build., vol. 33, no.4, pp. 319-331, abr. 2001. [ Links ]

[6] Google Maps (2020, Dic. 12). GoogleMaps Imata [Online]. Available: https://www.google.com.pe/maps/place/15%C2%B050'23.5%22S+71%C2%B005'34.8%22W/@-15.8398611,-71.094028,377m/data=!3m2!1e3!4b1!4m6!3m5!1s0x0:0x0!7e2!8m2!3d-15.839872!4d-71.0929901?hl=es-419&authuser=0.Links ]

[7] Ministerio de Vivienda, Construcción y Saneamiento (2014, May. 13). Confort Térmico y Lumínico con Eficiencia Energética. (Norma EM. 110) [Online]. Available: http://www3.vivienda.gob.pe/dnc/archivos/Estudios_Normalizacion/Normalizacion/normas/DS-006-2014-VIVIENDA.pdfLinks ]

[8] M. A. Eben Saleh, “Adobe as a thermal regulating material”, Sol. Wind. Tech., vol. 7, no. 4, pp. 407-416, 1990. [ Links ]

[9] Rotulos Peru (2019). Policarbonato Ficha Técnica [Online]. Available: https://rotulosperu.com/wp-content/uploads/2019/07/Ficha-Tecnica-policarbonato-Glanze.pdfLinks ]

[10] M. A. Humphreys, Field Studies of Thermal Comfort Compared and Applied (BRE-CP 76/75), Garston: Building Research Establishment, 1975. [ Links ]

[11] M. A. Humphreys, Comfortable indoor temperatures related to the outdoor air temperature, BRE Publication Draft PD 117/76, Garston: Building Research Establishment, 1976. [ Links ]

Recibido: 28 de Diciembre de 2020; Aprobado: 18 de Junio de 2021

Autor para la correspondencia: Mónica Gómez-León, Email mgomez@uni.edu.pe

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons