Noticias

May 04, 2023

Modelado detallado del colector solar de placa plana con acristalamiento al vacío

Fecha: 20 de marzo de 2023 Autores: Viacheslav Shemelin y Tomas Matuska Académico

Fecha: 20 de marzo de 2023

Autores: Viacheslav Shemelin y Tomas Matuska

Editor académico: Stoian Petrescu

Fuente: hindú | https://doi.org/10.1155/2017/1587592

Se presenta un análisis teórico de colectores solares de placa plana con acristalamiento al vacío. Se han investigado diferentes configuraciones del colector mediante un modelo teórico detallado basado en un balance combinado de energía externa e interna del absorbedor. Se han obtenido las características de rendimiento de las alternativas de colectores de placa plana al vacío. Posteriormente, se evaluaron las ganancias anuales de energía para una variante seleccionada y se compararon con colectores de tubo de vacío de última generación. Los resultados de la modelización indican que, en el caso de utilizar acristalamientos de vacío avanzados con revestimiento optimizado de baja emisividad (emisividad 0,20, transmitancia solar 0,85), es posible alcanzar parámetros de eficiencia similares o incluso mejores que los colectores de tubos de vacío. El diseño presentado en este documento puede considerarse prometedor para la extensión del rango de aplicabilidad de FPC y podría usarse en aplicaciones que requieren un nivel de temperatura bajo a medio.

El aprovechamiento térmico de la energía solar para la producción de calor a partir de la luz solar es uno de los métodos de transformación energética más antiguos. Esta tecnología se conoce y, a veces, incluso inconscientemente, se utiliza desde hace mucho tiempo. Ha sido redescubierto y utilizado de nuevo durante los últimos 45 años. Hoy en día, está listo para su aplicación, pero después de este corto tiempo de crecimiento, existe un gran potencial de desarrollo en este campo, especialmente en el campo de los colectores solares.

Actualmente, el tipo de colector solar más utilizado en Europa es el colector solar de placa plana (FPC). Estructura simple, alta eficiencia óptica, bajo costo y operación segura son sus principales características. Sin embargo, FPC generalmente está diseñado para un nivel de temperatura bajo entre 40 °C y 60 °C, que es principalmente el caso del sistema de agua caliente sanitaria. Cualquier cambio a un nivel de temperatura más alto podría traer la extensión del rango de aplicabilidad de FPC. Por lo tanto, se están realizando esfuerzos para mejorar el rendimiento de los colectores solares de placa plana. El rendimiento de un colector solar de placa plana está influenciado en gran medida por las pérdidas térmicas del absorbedor al ambiente a través de la cubierta transparente. Una forma de reducir esta pérdida de calor es reducir la transferencia de calor por convección natural en el espacio entre el absorbedor y la cubierta mediante su partición con el uso de paneles de vidrio adicionales, películas de plástico o materiales aislantes transparentes (TIM). Otra forma de reducir esta pérdida de calor es utilizar gas con una conductividad térmica más baja en lugar de aire o evacuar el espacio.

Veinberg BP y Veinberg VB [1] investigaron el uso de "mallas profundas y estrechas" como aislamiento solar transparente en forma de panal. Además, Hollands [2] presentó las características teóricas de rendimiento de un panal celular como un dispositivo de supresión de convección colocado entre el absorbedor y la cubierta de vidrio exterior del FPC. Tabor [3] presentó una breve descripción de la construcción de panal celular, indicando que un uso exitoso del aislamiento de panal requiere un material con mejores propiedades físicas y técnicas de fabricación. Más tarde, Rommel y Wagner [4] demostraron que el FPC que contiene capas de panal de policarbonato de 50–100 mm funciona bien con una temperatura de trabajo fluida entre 40 y 80 °C. Kesentini et al. [5] presentó un FPC con aislamiento de plástico transparente y sistema de protección contra sobrecalentamiento de bajo costo destinado al suministro de calor de 80 a 120 °C. También es posible alcanzar temperaturas de trabajo más altas de hasta 260 °C utilizando panales de vidrio, ya que las cubiertas de plástico son susceptibles de derretirse a temperaturas superiores a 120 °C.

Svendsen y Jensen [6] y Svendsen [7] demostraron experimentalmente que la eficiencia de la FPC solar se puede mejorar significativamente llenando el espacio de aire entre el absorbedor y la cubierta con un aerogel de sílice monolítico y evacuando a 10 kPa. Duan [8] estudió la reducción de la pérdida de calor del lado frontal colocando la capa de aerogel entre la cubierta transparente y la placa absorbente mostrando un aumento del 21% en la eficiencia del colector con respecto al colector convencional. Estos estudios han demostrado que las pérdidas de calor por convección se reducen significativamente con el uso de TIM debido a la partición del espacio entre el absorbedor y la cubierta que restringe el transporte de calor por convección y, por lo tanto, se logró un mayor rendimiento de FPC. Los resultados de las pruebas fueron alentadores y se obtuvo un rendimiento comparable al de los colectores de tubos de vacío.

Sin embargo, la mayoría de los materiales aislantes transparentes disponibles todavía no son una buena opción para los colectores de placa plana de alta temperatura. O bien no pueden soportar niveles de temperatura elevados porque están hechos de plástico (principalmente materiales de nido de abeja y capilares), son higroscópicos y no pueden soportar la humedad dentro del colector (aerogeles, etc.) o son muy caros (capilares de vidrio).

Se sabe que el uso de un vacío moderado en colectores de placa plana reduce las pérdidas de calor superiores desde el trabajo de Eaton y Blum [9]. El concepto de un colector de placa plana de vacío se realizó comercialmente y está disponible en el mercado. Además del mayor rendimiento térmico, estos colectores tienen la ventaja de una vida útil más prolongada en comparación con los colectores sin vacío, ya que no se producen problemas de humedad ni condensación dentro de la carcasa. Las presiones interiores típicas, que se pueden mantener económicamente, se encuentran entre 1 y 10 kPa. Significa que, aunque se suprimen las pérdidas por convección, la conducción de calor del gas permanece completamente desarrollada.

Además, Benz y Beikircher [10] construyeron un colector prototipo basado en el colector de placa plana disponible comercialmente. Para implementar una alta eficiencia térmica en el rango de temperatura media, se han reducido las pérdidas térmicas del absorbedor utilizando un absorbedor selectivo de baja emisividad, un relleno de criptón a baja presión (5 kPa) en la carcasa del colector. El colector prototipo ha sido probado dinámicamente y ha mostrado eficiencias muy altas de más del 60% a 100°C. Posteriormente, Benvenuti [11] presentó un FPC, que es capaz de alcanzar los 300°C. Eso ha sido posible gracias al vacío ultraalto (1,33x10⁻⁷Pa) mantenido por una bomba getter alimentada por el sol. En cuanto a los estudios más recientes, Moss y Shire [12] indican una mejora del 25 % para un FPC convencional al 60–65 % para FPC de vacío cuando se opera a 140 °C por encima de la temperatura ambiente. Más recientemente, Shire et al. [13] destaca que el colector FPC de vacío podría proporcionar calor hasta 200°C con una eficiencia superior al 50%.

Este artículo presenta la idea de sustituir el acristalamiento simple, que se utiliza en la mayoría de los colectores solares de placa plana, por un acristalamiento plano al vacío que, por un lado, presentará un bajo nivel de pérdidas de calor (revestimiento de baja emisividad, alto vacío). ) y, por otro lado, demostrará una alta transmitancia de energía solar. Los colectores solares de placa plana con baja pérdida de calor (al nivel de los colectores de tubos de vacío) y con una eficiencia óptica suficientemente alta podrían usarse de manera efectiva para la integración en envolventes de edificios (residenciales, industriales), que están ampliamente disponibles.

El acristalamiento al vacío consta de dos láminas de vidrio selladas entre sí alrededor de la periferia. Las láminas de vidrio están soportadas por un conjunto de pilares dispuestos en un patrón de rejilla cuadrada regular, y el espacio entre las láminas se evacua a una presión inferior a 0,1 Pa, eliminando así de manera efectiva tanto la conducción como la convección gaseosa. Mantener la presión por debajo de 0,1 Pa durante un período de vida útil previsto de 30 años representa un importante desafío de ingeniería. Colmillo et al. [14] llevaron a cabo pruebas extremas de ciclos térmicos.

Se encontró que la conductancia térmica en la región central había aumentado en un 10 %, a partir de lo cual se determinó que la presión de vacío dentro del espacio evacuado había aumentado desde el nivel insignificante de menos de 0,1 Pa a 0,16 Pa. Más tarde, Koebel et al. [15] investigó las posibles fuentes de aumento de la presión y concluyó que, teniendo en cuenta las condiciones ideales del proceso, debería ser posible mantener la presión total por debajo de 0,1 Pa después de 30 años. Es fundamental evitar cualquier tipo de fuga o microfisura porque la calidad del vacío está directamente ligada al rendimiento del aislamiento térmico.

Tres mecanismos diferentes de transferencia de calor contribuyen al coeficiente total de transferencia de calor hg₁-g₂ del acristalamiento: conducción térmica a través de un gas residual, conducción térmica a través de espaciadores y transferencia de calor por radiación entre las dos láminas en el acristalamiento al vacío. El coeficiente de transferencia de calor total hg₁-g₂ entre las láminas de vidrio de un acristalamiento al vacío se puede aproximar simplemente sumando los coeficientes de transferencia de calor individuales como [16, 17]

donde P es la presión interna, σ es la constante de Stefan-Boltzmann (5,67x10⁻⁸W/m² x K⁴), Tmean es el promedio de las temperaturas T₁ y T₂ de las láminas de vidrio, λ es la conductividad térmica de los pilares de vidrio, r es la radio de los pilares de vidrio, d es la distancia entre los pilares, y la emitancia efectiva, εeff, se escribe convencionalmente de la siguiente manera:

La ecuación (1) es válida para el espacio con presión inferior a 0,1 Pa, es decir, alto vacío.

El acristalamiento al vacío ya aparece en la cartera de proveedores de ventanas para uso en edificios (ver Figura 1). El acristalamiento al vacío comercialmente disponible con una configuración de 3-0,2-3 mm tiene una transmitancia de la radiación solar de τ=62 % y un valor U de transmitancia en el centro del acristalamiento de 1,1 W/m² x K. El vacío proporciona un alto nivel de aislamiento térmico que reduce la pérdida de calor por conducción y convección. Sin embargo, la pérdida de calor aún puede ocurrir debido a la radiación. El uso de un revestimiento de baja emisividad en la superficie exterior del panel de vidrio interior minimiza las pérdidas de calor por radiación.

Los recubrimientos low-e de última generación se han desarrollado casi exclusivamente para la arquitectura. Para mantener el confort térmico y visual en los edificios, se utilizan principalmente sistemas de revestimiento a base de plata, que pueden proporcionar una emisividad extremadamente baja (menos de 0,03) y una alta transmitancia visible (hasta 0,90). Sin embargo, la transmitancia solar (dentro de todo el rango del espectro solar) rara vez es superior a 0,60. Una transmisión de energía solar baja, provocada por la reflectancia del revestimiento de baja emisividad para la radiación del infrarrojo cercano (NIR) en el espectro solar, no es adecuada para su uso en colectores solares. Sin embargo, se pueden conseguir valores de hasta 0,80 y la correspondiente emisividad superior (entre 0,15 y 0,20) utilizando capas de plata muy finas, que se han desarrollado en los últimos años para triple acristalamiento, o con óxidos metálicos (Figura 2) [18] . También es posible utilizar los revestimientos antirreflectantes de vidrio externo en ambas superficies y, por lo tanto, reducir la reflexión en las dos superficies límite de aire y vidrio.

Para evaluar el potencial de la aplicación del acristalamiento al vacío en el diseño de colectores solares térmicos, se han llevado a cabo simulaciones detalladas para tres variantes diferentes del acristalamiento de la cubierta del colector. La variante de referencia (REF) es un simple vidrio solar bajo en hierro. La segunda variante (VG1) tiene un acristalamiento al vacío basado en dos vidrios bajos en hierro sin ningún tipo de recubrimiento. La última variante (VG2) es un acristalamiento de vacío avanzado con un revestimiento de baja emisividad en la superficie exterior del vidrio interior (posición 3). Las propiedades ópticas del recubrimiento son emisividad IR de 0,2 y transmitancia solar de 0,85. Las configuraciones de los acristalamientos de colector considerados se muestran gráficamente en la Figura 3. Los parámetros de los acristalamientos de cubierta utilizados para el estudio comparativo se enumeran en la Tabla 1.

Tabla 1 Comparación de las propiedades físicas de los acristalamientos de cobertura.

3.1. Descripción del modelo

Para analizar el rendimiento térmico de FPC con variantes consideradas de la cubierta transparente, se ha utilizado un modelo teórico detallado del colector de placa plana. El modelo detallado se originó a partir de la herramienta de diseño KOLEKTOR 2.0 [19] desarrollada originalmente como el programa Visual Basic.

Para comparar el rendimiento de los acristalamientos dados, se ha considerado que el colector de placa plana consiste en un absorbente colocado en la caja aislada cubierta con una cubierta transparente dada. Existe un espacio de aire entre el absorbedor y la cubierta y entre el absorbedor y su aislamiento posterior, ambos definidos por el espesor y la pendiente. El absorbedor está diseñado como un arpa con tuberías de distribución y subida (definidas por longitud, distancia y diámetro). La cubierta transparente se considera con conductancia térmica dependiente de la temperatura según la ecuación

donde hg₀, hg₁ y hg₂ son coeficientes para estructuras de cubierta dadas y Tg₁-g₂ (°C) es la temperatura media de acristalamiento. Las capas de aislamiento térmico se consideran de manera similar como dependientes de la temperatura.

El modelo detallado del colector solar de placa plana permite realizar un cálculo detallado de la transferencia de calor en el colector solar. El flujo de energía desde la superficie del absorbedor al ambiente y desde la superficie del absorbedor a un líquido de transferencia de calor, junto con una distribución de temperatura en el colector, se calculan en los bucles de iteración. Un colector solar se puede especificar mediante una serie de parámetros detallados, propiedades ópticas del acristalamiento y absorbente, y propiedades termofísicas de los componentes principales de un colector solar (marco, absorbente y cubierta transparente) en el modelo.

3.2. Ecuaciones Básicas

El modelo matemático para el colector solar de líquido de placa plana resuelve los balances de transferencia de calor unidimensionales. Hottel y Woertz [20], Hottel y Whillier [21] y Bliss [22] desarrollaron las suposiciones más simples: se desprecian las capacidades térmicas y se considera un valor único del coeficiente de pérdida de calor total del colector. Con base en estos supuestos y considerando que la transferencia de calor es principalmente unidimensional y predominante en la dirección normal al absorbedor, Duffie y Beckman [23] desarrollaron un modelo simplificado (con la analogía eléctrica) para caracterizar el colector solar en estado estacionario. condiciones. El modelo resuelve el balance de energía del colector solar en condiciones de estado estacionario de acuerdo con el principio de la ecuación de Hottel-Whillier para la producción térmica utilizable:

En esta ecuación, Aabs es el área del absorbedor, FR es el factor de eliminación de calor del colector, τ es la transmitancia solar de la cubierta del colector, α es la absorbancia solar del absorbedor, Gt es la radiación solar total, U es la pérdida de calor total coeficiente del colector, Tin es la temperatura del fluido de entrada y Tamb es la temperatura ambiente.

Los planos principales del colector son la superficie exterior de la cubierta (f₂), la superficie interior de la cubierta (f₁), el absorbedor (abs), la superficie interior del aislamiento posterior (b₁), la superficie exterior del marco posterior (b₂), la superficie interior del aislamiento del borde (e₁), y superficie exterior del marco de borde (e₂). Se determina una temperatura superficial para cada plano del colector durante el procedimiento de cálculo. Los principales planos colectores se describen esquemáticamente en la Figura 4.

El modelo matemático en general consta de dos partes: balance de energía externa del absorbedor (transferencia de calor desde la superficie del absorbedor al entorno ambiental) y balance de energía interna del absorbedor (transferencia de calor desde la superficie del absorbedor al fluido de transferencia de calor). Tanto el balance de energía externo como el interno son mutuamente dependientes. El coeficiente total de pérdida de calor del colector U como salida principal del balance externo es una de las entradas para el balance interno. Por otro lado, la temperatura media del absorbedor Tabs como una de las salidas del balance interno se usa como entrada necesaria para el balance externo. Se introdujo un ciclo de iteración para transferir los resultados del balance externo al balance interno inicial y los resultados del balance interno se colocan en el balance externo. El bucle itera siempre que la diferencia entre las temperaturas del absorbedor calculadas en dos pasos de iteración adyacentes sea mayor que el mínimo requerido (consulte la Figura 5).

3.3. Validación Experimental

El modelo ha sido validado experimentalmente en el marco de pruebas de colectores solares de acuerdo con la norma europea EN ISO 9806 en el Laboratorio Solar acreditado operado por el Centro Universitario para Edificios Eficientes en Energía, Universidad Técnica Checa en Praga. Los colectores solares térmicos han sido probados para obtener una salida térmica de estado estacionario en condiciones de funcionamiento constantes de temperatura de entrada (±0,05 K) y tasa de flujo másico (±0,002 %) del fluido de transferencia de calor (agua) que ingresa al colector y en condiciones climáticas constantes de energía solar. irradiación (±1,4%) y temperatura ambiente (±0,05 K).

La eficiencia instantánea se ha calculado a partir de la producción térmica del colector en relación con la entrada de radiación solar total (incidente en el área de referencia del colector: área bruta). Los puntos de datos experimentales de la eficiencia del colector solar se combinan con barras de incertidumbre uniformes en los gráficos. La incertidumbre expandida de eficiencia y la diferencia de temperatura reducida se evaluaron para datos experimentales de incertidumbres de tipo A (estadísticas) y tipo B (instrumentales) considerando el factor de cobertura k = 2 con un nivel de confianza del 95% (distribución normal).

El cálculo teórico de la característica de eficiencia por parte del modelo está sujeto a la incertidumbre de los parámetros reales del colector que se utilizan como entradas para el modelo. Si bien los parámetros geométricos están fácilmente disponibles con un alto grado de confianza, el número de parámetros que definen las propiedades de las partes del colector se encuentra incierto dentro de un rango estrecho (p. ej., parámetros ópticos de absorción y acristalamiento, en su mayoría ±2%), rango medio (p. ej., conductividad de la capa de aislamiento depende de su temperatura y densidad, ±10 %), y un rango bastante amplio (p. ej., la emitancia de la parte posterior del absorbedor y la emitancia de la capa de aislamiento o el marco del colector, >10 %). Por lo tanto, los resultados del cálculo teórico podrían presentarse como dos curvas delimitadoras donde los valores de eficiencia del colector se pueden encontrar en la realidad.

El modelo matemático ha sido validado en el campo de los colectores solares atmosféricos de placa plana (colectores solares de primera calidad con absorbedor de cobre soldado con láser de última generación recubierto con un recubrimiento selectivo de alto rendimiento y acristalamiento solar como cubierta transparente). Se han utilizado cuatro colectores solares diferentes para la validación detallada del modelo. La mayoría de los parámetros del colector solar térmico (p. ej., la conductividad térmica del aislamiento, la transmitancia solar del acristalamiento y la emisividad del absorbedor) se han medido experimentalmente para reducir el rango de incertidumbre. El modelo también ha sido probado en el caso de varios valores de pendiente, caudal másico, velocidad del viento y radiación incidente. Se puede encontrar más información sobre la validación del modelo en Shemelin y Matuska [24]. La Figura 6 muestra puntos de eficiencia medidos experimentalmente y características de eficiencia modeladas teóricamente. Es evidente a partir de los resultados que las características de eficiencia simuladas se ajustan relativamente bien a las medidas, lo que da confianza sobre el modelo desarrollado.

Se modelaron cuatro configuraciones de FPC con dimensiones 1 × 2 m. Todas las variantes tienen un absorbedor común soldado con láser de cobre de última generación con un recubrimiento selectivo de alto rendimiento con una absortividad de 0,95 y una emisividad de 0,05. El espesor del absorbedor es de 0,2 mm y la distancia entre tubos es de 100 mm. El aislamiento térmico trasero y el aislamiento térmico de borde con conductividad térmica de 0,04 W/m·K tienen un espesor de 50 y 20 mm, respectivamente. El espesor del entrehierro entre absorbedor y acristalamiento es de 30 mm (excepto variante VC4, explicada a continuación). La diferencia entre las variantes consideradas está solo en el acristalamiento de la cubierta.

La primera configuración de FPC (RC) tiene un acristalamiento de referencia (REF) a modo de tapa transparente con los parámetros que se muestran en la Tabla 1. La segunda configuración VC1 tiene un acristalamiento al vacío VG1 en lugar del acristalamiento de referencia. La tercera configuración VC2 considera la variante VG2 como acristalamiento de cobertura. La última configuración VC3 tiene una configuración principalmente diferente. Esta variante no tiene espacio entre el absorbedor y la tapa. Aquí, el absorbedor está unido al VG2 mediante un gel de silicona altamente transparente y permanentemente flexible para reducir el grosor del colector (alternativa de colector "delgado" adecuada para la integración de la envolvente del edificio). Las configuraciones consideradas de los colectores FPC RC, VC1, VC2 y VC3 se muestran en la Figura 7.

El gráfico de la Figura 8 muestra las curvas de eficiencia de las variantes de colectores solares consideradas. Con respecto a la norma EN ISO 9806, la eficiencia del colector η se basa en el área bruta del colector AG. Los cálculos se han realizado con el uso del modelo detallado teórico descrito anteriormente de FPC. El gráfico muestra la diferencia en la calidad de la energía de las variantes de FPC comparadas. La baja pendiente de la curva de eficiencia del VC2 se debe a que el colector VC2 tiene dos revestimientos de baja emisividad: el primero dentro del acristalamiento de vacío y el segundo en la superficie absorbente. Otras variantes RC, VC1 y VC3 solo tienen un revestimiento de baja emisividad. La alternativa de colector solar VC3 tiene una eficiencia térmica más baja que la alternativa de referencia REF. A pesar de la capa de vacío, la mayor emitancia del revestimiento de vidrio de baja emisividad (0,2) en la variante VC3 en lugar del revestimiento absorbente con una emitancia de 0,05 en la variante REF, lleva la pérdida de calor superior total a un valor similar, pero los parámetros ópticos de la configuración VC3 son más bajos (menor eficiencia de pérdida cero η₀). Los coeficientes resultantes de las características de eficiencia η₀, a₁ y a₂ se enumeran en la Tabla 2.

Tabla 2 Resumen de los resultados de la simulación del colector.

Por otro lado, la variante de colector solar VC3 tiene el grosor más bajo entre las variantes comparadas: solo 60 mm. Tal grosor brinda más posibilidades de integración de FPC en la envolvente del edificio debido al diseño delgado y compacto. Otras variantes RC, VC1 y VC2 tienen espesores de 87, 91 y 91 mm, respectivamente.

La Figura 9 presenta la comparación de las características de eficiencia de la variante de colector de placa plana VC2 y los colectores de tubo de vacío (con/sin reflector, absorbedor cilíndrico/plano) en relación con el área bruta de un colector. El rendimiento térmico de la variante de colector solar de placa plana VC2 es comparable al de los colectores de tubo de vacío. Además, la variante de colector solar VC2 muestra una eficiencia significativamente mayor que la mayoría de los colectores de tubos de vacío (VT).

Para obtener una imagen completa, el rendimiento anual de los colectores VC2 y VT se ha modelado utilizando el software ScenoCalc [25] para temperaturas de funcionamiento constantes de 25, 50, 75 y 100 °C y las condiciones climáticas de Würzburg. El rendimiento y las características ópticas de los colectores comparados se utilizaron como datos de entrada. Los resultados se muestran en la Tabla 3. Los resultados del modelado confirmaron que la variante de colector solar de placa plana VC2 tiene mayores ganancias de energía solar que la mayoría de los colectores de tubo de vacío hasta una temperatura de funcionamiento de 100 °C.

Tabla 3 Ganancia anual calculada del colector solar con respecto al área bruta del colector.

Se han investigado teóricamente diferentes diseños de colectores solares de placa plana basados ​​en un acristalamiento de vacío plano utilizando el modelo matemático detallado para mostrar el potencial de la aplicación del acristalamiento de vacío en colectores solares de placa plana. La variante seleccionada VC2 se ha comparado con colectores de tubo de vacío de última generación mediante la simulación anual de la producción de calor del colector en ScenoCalc. Los resultados han demostrado que existe un potencial significativo para aumentar la eficiencia de los colectores solares de placa plana mediante el uso de acristalamiento de vacío de alto rendimiento como cubierta transparente. La variante de colector VC2 muestra un mayor rendimiento que la mayoría de los colectores de tubos de vacío hasta una temperatura de funcionamiento de 100 °C.

Los autores declaran que no existe ningún conflicto de intereses con respecto a la publicación de este artículo.

Este trabajo ha sido apoyado por el Ministerio de Educación, Juventud y Deportes dentro del Programa Nacional de Sustentabilidad I (NPU I), Proyecto no. LO1605, por el Centro Universitario de Edificaciones Energéticamente Eficientes-Fase de Sostenibilidad, y por el proyecto SGS16/212/OHK2/3T/12—Modelado, control y diseño de instalaciones de ingeniería ambiental.