ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
QUASI-DYNAMIC TEST OF FLAT PLATE SOLAR COLLECTORS IN URUGUAY ACCORDING TO THE ISO 9806:2017 STANDARD
Juan Manuel Rodríguez Muñoz
Universidad de la República (UDELAR). Centro Universitario Litoral Norte. Departamento de Física del Litoral. Laboratorio de Energía Solar. Uruguay.
jrodrigue@fing.edu.uy
http://orcid.org/0000-0001-8523-0582
Andrés Monetta
Universidad de la República (UDELAR). Centro Universitario Litoral Norte. Departamento de Física del Litoral. Laboratorio de Energía Solar. Uruguay.
andres130402@gmail.com
http://orcid.org/0000-0003-1087-2569
Italo Carlos Bove Vanzulli
Universidad de la República (UDELAR). Facultad de Ingeniería. Instituto de Física. Laboratorio de Energía Solar. Uruguay.
italo@fing.edu.uy
http://orcid.org/0000-0001-8450-8183
Rodrigo Alonso-Suárez
Universidad de la República (UDELAR). Facultad de Ingeniería. Instituto de Física. Laboratorio de Energía Solar. Uruguay.
r.alonso.suarez@gmail.com
http://orcid.org/0000-0002-5216-5232
Recibido: 30/05/2020 y Aceptado: 03/08/2020
ENERLAC. Volumen IV. Número 2. Diciembre, 2020 (10 - 26)
ISSN: 2602-8042 (impreso) / 2631-2522 (digital)
QUASI-DYNAMIC TEST OF FLAT PLATE SOLAR COLLECTORS IN URUGUAY ACCORDING TO THE ISO 9806:2017 STANDARD
Juan Manuel Rodríguez Muñoz
Universidad de la República (UDELAR). Centro Universitario Litoral Norte. Departamento de Física del Litoral. Laboratorio de Energía Solar. Uruguay.
jrodrigue@fing.edu.uy
http://orcid.org/0000-0001-8523-0582
Andrés Monetta
Universidad de la República (UDELAR). Centro Universitario Litoral Norte. Departamento de Física del Litoral. Laboratorio de Energía Solar. Uruguay.
andres130402@gmail.com
http://orcid.org/0000-0003-1087-2569
Italo Carlos Bove Vanzulli
Universidad de la República (UDELAR). Facultad de Ingeniería. Instituto de Física. Laboratorio de Energía Solar. Uruguay.
italo@fing.edu.uy
http://orcid.org/0000-0001-8450-8183
Rodrigo Alonso-Suárez
Universidad de la República (UDELAR). Facultad de Ingeniería. Instituto de Física. Laboratorio de Energía Solar. Uruguay.
r.alonso.suarez@gmail.com
http://orcid.org/0000-0002-5216-5232
Recibido: 30/05/2020 y Aceptado: 03/08/2020
ENERLAC. Volumen IV. Número 2. Diciembre, 2020 (10 - 26)
ISSN: 2602-8042 (impreso) / 2631-2522 (digital)
Foto de Samyag Shah de Unsplash.
RESUMEN
Los ensayos de colectores solares en estado transitorio permiten aumentar la capacidad de ensayo de los laboratorios respecto a los ensayos en estado estacionario. Esto ha sido validado en laboratorios internacionales, especialmente en Europa. En Latinoamérica este método no ha tenido mucha difusión y a la fecha ningún laboratorio de esta región había validado (ni implementado) su uso en climas locales, por lo que los laboratorios de la región que están desarrollando capacidades de ensayo de colectores solares no consideran aún este método en forma comercial. En este trabajo se presentan los resultados de un ensayo en estado transitorio (cuasi dinámico) de un colector solar de placa plana de acuerdo a la norma ISO 9806:2017 en Uruguay, región de la Pampa Húmeda en Sudamérica. Se compara sus resultados con la metodología de ensayo tradicional en estado estacionario y se demuestra que es posible aplicarlo en la región incluso con la ventaja de ampliar la cantidad anual de ensayos, como se ha observado en Europa.
Palabras clave: Colectores Solares Térmicos, Ensayos de Desempeño, Modelos Dinámicos, Estado Estacionario, ISO 9806, Uruguay
Abstract
Transient-state efficiency tests for solar thermal collectors allow increasing the laboratories annual test capacity in comparison to steady-state efficiency tests. This characteristic has been evaluated by international laboratories, especially in Europe. In Latin America this method has not been widely disseminated and to date no regional laboratory had validated (or implemented) its utilization under the local climates, so laboratories in the region that are developing solar collector testing capabilities do not yet consider this method commercially. This article presents the results of a transient-state (quasi-dynamic) test of a flat plate solar collector according to ISO 9806:2017 standard in Uruguay, a part of the broader Pampa Húmeda region in South America. Test results are compared to those of the traditional steady-state methodology, demonstrating that the method can be applied in the region and has also the advantage of increasing the annual test capacity, as observed in Europe.
Keywords: Solar Thermal Collectors, Performance Test, Dynamic Model, Steady State, ISO 9806, Uruguay.
Los ensayos de colectores solares en estado transitorio permiten aumentar la capacidad de ensayo de los laboratorios respecto a los ensayos en estado estacionario. Esto ha sido validado en laboratorios internacionales, especialmente en Europa. En Latinoamérica este método no ha tenido mucha difusión y a la fecha ningún laboratorio de esta región había validado (ni implementado) su uso en climas locales, por lo que los laboratorios de la región que están desarrollando capacidades de ensayo de colectores solares no consideran aún este método en forma comercial. En este trabajo se presentan los resultados de un ensayo en estado transitorio (cuasi dinámico) de un colector solar de placa plana de acuerdo a la norma ISO 9806:2017 en Uruguay, región de la Pampa Húmeda en Sudamérica. Se compara sus resultados con la metodología de ensayo tradicional en estado estacionario y se demuestra que es posible aplicarlo en la región incluso con la ventaja de ampliar la cantidad anual de ensayos, como se ha observado en Europa.
Palabras clave: Colectores Solares Térmicos, Ensayos de Desempeño, Modelos Dinámicos, Estado Estacionario, ISO 9806, Uruguay
Abstract
Transient-state efficiency tests for solar thermal collectors allow increasing the laboratories annual test capacity in comparison to steady-state efficiency tests. This characteristic has been evaluated by international laboratories, especially in Europe. In Latin America this method has not been widely disseminated and to date no regional laboratory had validated (or implemented) its utilization under the local climates, so laboratories in the region that are developing solar collector testing capabilities do not yet consider this method commercially. This article presents the results of a transient-state (quasi-dynamic) test of a flat plate solar collector according to ISO 9806:2017 standard in Uruguay, a part of the broader Pampa Húmeda region in South America. Test results are compared to those of the traditional steady-state methodology, demonstrating that the method can be applied in the region and has also the advantage of increasing the annual test capacity, as observed in Europe.
Keywords: Solar Thermal Collectors, Performance Test, Dynamic Model, Steady State, ISO 9806, Uruguay.
INTRODUCCIÓN
El aumento sostenido del consumo final de energía mundial sumado al hecho de que el 85% proviene de combustibles fósiles (IEA, 2016a) genera gran preocupación respecto a las emisiones de CO2 y su efecto en el cambio climático. Las dos acciones más importantes para reducir estas emisiones son la incorporación de fuentes de energías renovables al uso final de energía y la reducción de la demanda a través de medidas de eficiencia energética. En 2016, el sector edilicio fue responsable del 30 % del consumo final de energía mundial (IEA, 2016a) mientras que en Uruguay este sector represento el 25 % (IEA, 2016b). Una parte importante de este consumo se debe al calentamiento de agua sanitaria (ACS) y a los sistemas de calefacción y refrigeración.
La energía solar térmica de baja temperatura se presenta como una candidata para sustituir el uso tradicional de energía para ACS y para los sistemas de acondicionamiento térmico en el sector edilicio. Entre 2007 y 2017, la capacidad instalada de colectores solares térmicos aumentó de forma sostenida pasando de 145 GW en 2007 a 472 GW en 2017 (REN21, 2018). No obstante, si bien la energía solar térmica ha crecido de forma notoria, la participación de esta tecnología en la matriz energética mundial sigue siendo pequeña. Para lograr una mayor participación se deben abordar distintas problemáticas, una de ellas es la calidad de los productos que se comercializan. En este sentido, varias organizaciones internacionales de reconocido prestigio, como la International Organization for Standarization (ISO) y la American Society of Heating, Refrigeration and Air-Conditioning Enginners (ASHRAE), han elaborado normativas que describen procedimientos para evaluar el desempeño térmico de colectores solares y la calidad de los materiales con los que se fabrican. El hecho que un país cuente con estas capacidades de ensayo no garantiza que la capacidad instalada de colectores solares aumente, dado que esto dependerá en definitiva de las políticas energéticas que adopte cada país. No obstante, dicha capacidad es necesaria para el desarrollo del sector, ya que difícilmente una política energética tendrá éxito si la calidad de los colectores comercializados no es la adecuada. Además, los ensayos de desempeño térmico constituyen un insumo básico para la evaluación económico financiera de los proyectos.
Tradicionalmente los ensayos de eficiencia de colectores solares (desempeño térmico) se han realizado en estado estacionario (Steady State Testing, SST). Para asegurar el estado estacionario se requiere mantener con gran estabilidad diversas variables durante el ensayo. Algunas de ellas son controlables, como la temperatura del fluido a la entrada del colector y el caudal, y otras son no controlables, como la irradiancia solar o las condiciones de viento sobre el colector. Es posible mitigar la acción del viento en los alrededores forzando viento de forma artificial a una velocidad determinada sobre el colector, elemento que es considerado por todos los laboratorios al exterior. El control del ángulo de incidencia es también una limitante, dado que el movimiento aparente del sol varía para cada momento del año. Por ende, dada todas las restricciones geométricas, no es sencillo conseguir los ángulos de incidencia necesarios para un ensayo bajo norma, incluso utilizando sistemas de seguimiento solar. Estas condicionantes, en especial, las relacionadas con irradiancia solar y ángulo de incidencia, limitan la cantidad de ensayos que un laboratorio puede realizar en un año.
Buscando solucionar esta situación, en la década de los 90’ se realizaron los primeros ensayos de eficiencia en estado transitorio (Quasi Dynamic Testing, QDT) para colectores solares (Muschaweck y Spirkl, 1993; Perers, 1993) y se comenzaron a evidenciar las ventajas de este método en sitios con nubosidad variable, donde las condiciones meteorológicas para el ensayo SST son difíciles de conseguir. Algunos laboratorios en Europa ubicados en sitios con estas características climáticas adoptaron esta metodología (Fischer et al., 2004; Osório y Carvalho, 2014; García de Jalón, 2011; Zambolin, 2010) y también en Estados Unidos (Rojas, 2008), pero en Latinoamérica, por distintos motivos, esta metodología no ha sido adoptada por la mayoría de los laboratorios, existiendo muy pocos antecedentes a la fecha (Kratzenberg, 2006). Dado que el método SST requiere condiciones exigentes de cielo despejado para su correcta implementación, el método QDT se presenta como una alternativa más factible de realizarse en climas con mayor presencia de nubosidad, donde los días de cielo despejados son escasos. Dadas sus características, la metodología QDT no sólo permite, sino que es un requisito de la misma, que durante el ensayo hayan ocurrido diferentes condiciones de cielo (cielo despejado, parcialmente nublado y nublado), siendo difícil de implementar en sitios con predominancias de cielos despejados (en climas áridos o desérticos, por ejemplo). Los laboratorios pueden optar por utilizar una u otra metodología según sean las condiciones meteorológicas locales (utilizando la misma infraestructura de ensayo).
El objetivo central de este artículo es demostrar la aplicabilidad del método QDT en el clima de la región conocida como la Pampa Húmeda, sudeste de Sudamérica, que abarca la región climática afín de Uruguay, las provincias de Buenos Aires, Entre Ríos, Santa Fe, Córdoba y la Pampa, y el sur de Brasil. Para ello se realizan dos cosas: (1) se presenta en este trabajo el ensayo de un colector solar de placa plana en Uruguay bajo las metodologías QDT y SST, siguiendo normativa ISO vigente, y se muestra que sus resultados de ensayo son equivalentes, y (2) se proporciona un análisis de la cantidad de ensayos por mes y por año de la metodología (QDT) en la región, encontrándose que esta metodología es capaz de aportar más del doble de ensayos anuales en la región en comparación con la metodología tradicional. El trabajo busca además darle mayor difusión a la metodología de ensayo QDT en Latinoamérica, describiendo una implementación del método realizada y validada por un laboratorio local.
Este artículo se organiza de la siguiente forma. En primer lugar se describe el modelado matemático de la norma ISO 9806:2017 para describir el comportamiento térmico de un colector solar con cubierta, incluyendo una descripción de las variables de ensayo relevantes. En la siguiente sección (p.17) se presenta brevemente la plataforma de ensayo, a saber, el Banco de Ensayos de Calentadores Solares (BECS) del Laboratorio de Energía Solar (LES) de la Udelar, y se describen las características comunes a ambos ensayos. Posteriormente, (p. 20) se describe el procedimiento de ensayo utilizado para ambas metodologías (SST y QDT) y la siguiente sección (p. 23) presenta y compara sus resultados. Se presenta al final de la misma el análisis de la cantidad de ensayos anuales que se podrían obtener con ambas metodologías, evidenciando las ventajas del método QDT en la región. Finalmente, se resumen las conclusiones (p. 25).
MODELADO TÉRMICO DEL COLECTOR
La norma ISO 9806:2017 (ISO, 2017) propone dos modelos matemáticos para describir el desempeño térmico de un colector solar con cubierta. El primer modelo es el más simple y se describe según:
El aumento sostenido del consumo final de energía mundial sumado al hecho de que el 85% proviene de combustibles fósiles (IEA, 2016a) genera gran preocupación respecto a las emisiones de CO2 y su efecto en el cambio climático. Las dos acciones más importantes para reducir estas emisiones son la incorporación de fuentes de energías renovables al uso final de energía y la reducción de la demanda a través de medidas de eficiencia energética. En 2016, el sector edilicio fue responsable del 30 % del consumo final de energía mundial (IEA, 2016a) mientras que en Uruguay este sector represento el 25 % (IEA, 2016b). Una parte importante de este consumo se debe al calentamiento de agua sanitaria (ACS) y a los sistemas de calefacción y refrigeración.
La energía solar térmica de baja temperatura se presenta como una candidata para sustituir el uso tradicional de energía para ACS y para los sistemas de acondicionamiento térmico en el sector edilicio. Entre 2007 y 2017, la capacidad instalada de colectores solares térmicos aumentó de forma sostenida pasando de 145 GW en 2007 a 472 GW en 2017 (REN21, 2018). No obstante, si bien la energía solar térmica ha crecido de forma notoria, la participación de esta tecnología en la matriz energética mundial sigue siendo pequeña. Para lograr una mayor participación se deben abordar distintas problemáticas, una de ellas es la calidad de los productos que se comercializan. En este sentido, varias organizaciones internacionales de reconocido prestigio, como la International Organization for Standarization (ISO) y la American Society of Heating, Refrigeration and Air-Conditioning Enginners (ASHRAE), han elaborado normativas que describen procedimientos para evaluar el desempeño térmico de colectores solares y la calidad de los materiales con los que se fabrican. El hecho que un país cuente con estas capacidades de ensayo no garantiza que la capacidad instalada de colectores solares aumente, dado que esto dependerá en definitiva de las políticas energéticas que adopte cada país. No obstante, dicha capacidad es necesaria para el desarrollo del sector, ya que difícilmente una política energética tendrá éxito si la calidad de los colectores comercializados no es la adecuada. Además, los ensayos de desempeño térmico constituyen un insumo básico para la evaluación económico financiera de los proyectos.
Tradicionalmente los ensayos de eficiencia de colectores solares (desempeño térmico) se han realizado en estado estacionario (Steady State Testing, SST). Para asegurar el estado estacionario se requiere mantener con gran estabilidad diversas variables durante el ensayo. Algunas de ellas son controlables, como la temperatura del fluido a la entrada del colector y el caudal, y otras son no controlables, como la irradiancia solar o las condiciones de viento sobre el colector. Es posible mitigar la acción del viento en los alrededores forzando viento de forma artificial a una velocidad determinada sobre el colector, elemento que es considerado por todos los laboratorios al exterior. El control del ángulo de incidencia es también una limitante, dado que el movimiento aparente del sol varía para cada momento del año. Por ende, dada todas las restricciones geométricas, no es sencillo conseguir los ángulos de incidencia necesarios para un ensayo bajo norma, incluso utilizando sistemas de seguimiento solar. Estas condicionantes, en especial, las relacionadas con irradiancia solar y ángulo de incidencia, limitan la cantidad de ensayos que un laboratorio puede realizar en un año.
Buscando solucionar esta situación, en la década de los 90’ se realizaron los primeros ensayos de eficiencia en estado transitorio (Quasi Dynamic Testing, QDT) para colectores solares (Muschaweck y Spirkl, 1993; Perers, 1993) y se comenzaron a evidenciar las ventajas de este método en sitios con nubosidad variable, donde las condiciones meteorológicas para el ensayo SST son difíciles de conseguir. Algunos laboratorios en Europa ubicados en sitios con estas características climáticas adoptaron esta metodología (Fischer et al., 2004; Osório y Carvalho, 2014; García de Jalón, 2011; Zambolin, 2010) y también en Estados Unidos (Rojas, 2008), pero en Latinoamérica, por distintos motivos, esta metodología no ha sido adoptada por la mayoría de los laboratorios, existiendo muy pocos antecedentes a la fecha (Kratzenberg, 2006). Dado que el método SST requiere condiciones exigentes de cielo despejado para su correcta implementación, el método QDT se presenta como una alternativa más factible de realizarse en climas con mayor presencia de nubosidad, donde los días de cielo despejados son escasos. Dadas sus características, la metodología QDT no sólo permite, sino que es un requisito de la misma, que durante el ensayo hayan ocurrido diferentes condiciones de cielo (cielo despejado, parcialmente nublado y nublado), siendo difícil de implementar en sitios con predominancias de cielos despejados (en climas áridos o desérticos, por ejemplo). Los laboratorios pueden optar por utilizar una u otra metodología según sean las condiciones meteorológicas locales (utilizando la misma infraestructura de ensayo).
El objetivo central de este artículo es demostrar la aplicabilidad del método QDT en el clima de la región conocida como la Pampa Húmeda, sudeste de Sudamérica, que abarca la región climática afín de Uruguay, las provincias de Buenos Aires, Entre Ríos, Santa Fe, Córdoba y la Pampa, y el sur de Brasil. Para ello se realizan dos cosas: (1) se presenta en este trabajo el ensayo de un colector solar de placa plana en Uruguay bajo las metodologías QDT y SST, siguiendo normativa ISO vigente, y se muestra que sus resultados de ensayo son equivalentes, y (2) se proporciona un análisis de la cantidad de ensayos por mes y por año de la metodología (QDT) en la región, encontrándose que esta metodología es capaz de aportar más del doble de ensayos anuales en la región en comparación con la metodología tradicional. El trabajo busca además darle mayor difusión a la metodología de ensayo QDT en Latinoamérica, describiendo una implementación del método realizada y validada por un laboratorio local.
Este artículo se organiza de la siguiente forma. En primer lugar se describe el modelado matemático de la norma ISO 9806:2017 para describir el comportamiento térmico de un colector solar con cubierta, incluyendo una descripción de las variables de ensayo relevantes. En la siguiente sección (p.17) se presenta brevemente la plataforma de ensayo, a saber, el Banco de Ensayos de Calentadores Solares (BECS) del Laboratorio de Energía Solar (LES) de la Udelar, y se describen las características comunes a ambos ensayos. Posteriormente, (p. 20) se describe el procedimiento de ensayo utilizado para ambas metodologías (SST y QDT) y la siguiente sección (p. 23) presenta y compara sus resultados. Se presenta al final de la misma el análisis de la cantidad de ensayos anuales que se podrían obtener con ambas metodologías, evidenciando las ventajas del método QDT en la región. Finalmente, se resumen las conclusiones (p. 25).
MODELADO TÉRMICO DEL COLECTOR
La norma ISO 9806:2017 (ISO, 2017) propone dos modelos matemáticos para describir el desempeño térmico de un colector solar con cubierta. El primer modelo es el más simple y se describe según:
El
segundo modelo surge de separar la irradiancia global en sus
subcomponentes directa y difusa, y remplaza el modificador de ángulo
considerando su efecto distinto sobre ambas componentes, de la
siguiente manera:
Lo que da lugar a:
En el caso del método SST se determinan los coeficientes η0,hem, a1, a2, Khem(θ)
y C de la Ecuación (1) en forma separada a través de tres ensayos
independientes: (i) Ensayo de rendimiento en donde se determinan los
coeficientes η0,hem, a1 y a2; (ii) Ensayo de modificador por ángulo de incidencia en donde se determina Khem(θ);
y (iii) Ensayo de capacidad térmica en donde se determina la constante
C. Cabe destacar que los primeros dos ensayos son en estado
estacionario pero el tercero no. En el caso del método QDT se
determinan los coeficientes η0,b, a1, a2, Kb(θ), Kd y C de la Ecuación (3) en un único ensayo que integra los tres ensayos independientes del método SST.
Para el modificador por ángulo de incidencia Khem(θ) se utiliza la siguiente relación:
Para el modificador por ángulo de incidencia Khem(θ) se utiliza la siguiente relación:
En el anexo B de la norma ISO 9806:2017 se describe como determinar η0,b, Kb(θ) y Kd a partir η0,hem y Khem(θ) y viceversa. Para calcular Kd primero se asume que Khem(θ) ≅ Kb(θ) y luego Kb(θ) es promediado y normalizado sobre el ángulo sólido visto por el colector, esto es:
La constante η0,b
se calcula a partir de la Ecuación (2) asumiendo incidencia normal y
una fracción difusa de 15% en el plano del colector. Para calcular η0,b,
se utilizó la mejora propuesta por Kovacs et al. (2011) de utilizar en
la Ecuación (6) la fracción difusa promedio en el plano del colector
durante el ensayo de rendimiento en estado estacionario. Entendemos que
este procedimiento es el más apropiado.
CARACTERÍSTICAS GENERALES DE LOS ENSAYOS
Las pruebas se realizaron en el Banco de Ensayos de Calentadores Solares (BECS) del Laboratorio de Energía Solar (LES, http://les.edu.uy/) de la Universidad de la República (Udelar) ubicado en la localidad de Salto (31° sur), Uruguay. En la Figura 1 se muestra la ubicación del LES dentro de la región climática afín de la pampa húmeda (verde claro). Dicha región se ubica al sureste de Sudamérica y está clasificada como Cfa (clima templado-cálido, húmedo, con veranos calurosos) según la clasificación actualizada de Köppen-Geiger (Peel et al., 2007). Las zonas identificadas en verde oscuro corresponden a la clasificación Cfb, muy similar a Cfa, con la excepción de que los veranos son más frescos por la influencia de la costa.
Se utilizó el mismo banco de ensayos para los dos métodos (SST y QDT). Cabe destacar que dicho banco fue diseñado por investigadores del propio laboratorio, en base a las instalaciones preexistentes del Centro Nacional de Energías Renovables (CENER) de España. Recientemente esta capacidad de ensayo participó de una inter-comparación de laboratorios de ensayo a nivel de Latinoamérica organizada por el PTB (Physikalisch-Technische Bundesanstalt), Instituto Metrológico de Alemania. En esta inter-comparación (de ensayos SST) se obtuvo una puntuación satisfactoria en todos los parámetros medidos (valor absoluto de Z-score menor 2). El BECS fue el único laboratorio de Latinoamérica que presentó además resultados de ensayo bajo la metodología QDT, obteniendo valores satisfactorios en todos los parámetros medidos, excepto en el parámetro Kd, en el que se obtuvo un valor absoluto de Z-score entre 2 y 3. El informe final de dicha actividad se encuentra en (Fisher, 2020).
Las pruebas se realizaron en el Banco de Ensayos de Calentadores Solares (BECS) del Laboratorio de Energía Solar (LES, http://les.edu.uy/) de la Universidad de la República (Udelar) ubicado en la localidad de Salto (31° sur), Uruguay. En la Figura 1 se muestra la ubicación del LES dentro de la región climática afín de la pampa húmeda (verde claro). Dicha región se ubica al sureste de Sudamérica y está clasificada como Cfa (clima templado-cálido, húmedo, con veranos calurosos) según la clasificación actualizada de Köppen-Geiger (Peel et al., 2007). Las zonas identificadas en verde oscuro corresponden a la clasificación Cfb, muy similar a Cfa, con la excepción de que los veranos son más frescos por la influencia de la costa.
Se utilizó el mismo banco de ensayos para los dos métodos (SST y QDT). Cabe destacar que dicho banco fue diseñado por investigadores del propio laboratorio, en base a las instalaciones preexistentes del Centro Nacional de Energías Renovables (CENER) de España. Recientemente esta capacidad de ensayo participó de una inter-comparación de laboratorios de ensayo a nivel de Latinoamérica organizada por el PTB (Physikalisch-Technische Bundesanstalt), Instituto Metrológico de Alemania. En esta inter-comparación (de ensayos SST) se obtuvo una puntuación satisfactoria en todos los parámetros medidos (valor absoluto de Z-score menor 2). El BECS fue el único laboratorio de Latinoamérica que presentó además resultados de ensayo bajo la metodología QDT, obteniendo valores satisfactorios en todos los parámetros medidos, excepto en el parámetro Kd, en el que se obtuvo un valor absoluto de Z-score entre 2 y 3. El informe final de dicha actividad se encuentra en (Fisher, 2020).
Figura 1. Ubicación del LES dentro de la región de la pampa húmeda (verde)
Fuente: Mapa de Peel et al. (2007) para la aplicación Google Earth
Figura 2. Instalación de ensayo
Fuente: Elaboración propia
Para este trabajo se ensayó un colector solar térmico de placa plana de 2.01 m2
de área bruta entre el 17 y el 25 de noviembre de 2019. Los datos se
procesaron de diferente forma según el método empleado, SST o QDT. El
colector se montó sobre un seguidor móvil como se muestra en la Figura
2. La inclinación horizontal de este seguidor se puede ajustar
manualmente entre 0° y 75° y el azimut se puede ajustar de forma manual
o automática con un paso de 2 minutos entre -90° y 90°. Se impuso una
velocidad de 3 m/s (promedio espacial) sobre el colector utilizando los
ventiladores axiales que se muestra en Figura 2. El caudal de ensayo se
fijó en 2.4 lts/min y durante todas las pruebas el ángulo de
inclinación del seguidor se fijó en 45°, el azimut se ajustó de forma
manual o automática según el ensayo. En todos los casos, antes del
periodo de medición, al colector se lo hizo pasar por un periodo de
acondicionamiento previo de 15 minutos a la temperatura de ensayo
correspondiente; dicho período no fue incluido en la identificación de
parámetros de ensayo en ninguna de las metodologías.
En la Figura 3 se muestra un esquema simplificado de la instalación termo-hidráulica del BECS. La misma cuenta con tres circuitos independiente: (1) circuito primario, en verde, (2) circuito de calentamiento, en rojo, y (3) circuito de enfriamiento, en azul. Los tramos de cañería en negro en el circuito primario corresponden a las cañerías flexibles que se pueden observar en la Figura 2. El control temperatura a la entrada del colector (circuito primario) se realiza en dos etapas. Primero el fluido caliente a la salida del colector debe ser enfriado, proceso que se realiza en el intercambiador de calor IC1. Luego el fluido es calentado en forma precisa a la temperatura de entrada establecida para el ensayo a través del intercambiador de calor IC2. El circuito de enfriamiento utiliza agua a 10 °C que proviene de una enfriadora eléctrica de agua (no mostrada en el diagrama) y el circuito de calentamiento utiliza agua caliente que proviene del termotanque TTE1. Los caudales en los tres circuitos son regulados mediante válvulas reguladores electroneumáticas (V51, V52 y V53) accionadas por controladores PID (indicados en línea punteada). Todo el sistema de control se desarrolló localmente utilizando un PLC marca Siemens modelo S7-1200. El diagrama indica además la ubicación de los sensores de temperatura de agua (T1, T2, T3, T4 y T5), el sensor de temperatura ambiente (Tamb), los sensores de irradiancia global horizontal (GHI) y en plano inclinado (GTI), el sensor de irradiancia difusa horizontal (DHI), la medida de viento (V) y el forzador de viento (GV).
En la Figura 3 se muestra un esquema simplificado de la instalación termo-hidráulica del BECS. La misma cuenta con tres circuitos independiente: (1) circuito primario, en verde, (2) circuito de calentamiento, en rojo, y (3) circuito de enfriamiento, en azul. Los tramos de cañería en negro en el circuito primario corresponden a las cañerías flexibles que se pueden observar en la Figura 2. El control temperatura a la entrada del colector (circuito primario) se realiza en dos etapas. Primero el fluido caliente a la salida del colector debe ser enfriado, proceso que se realiza en el intercambiador de calor IC1. Luego el fluido es calentado en forma precisa a la temperatura de entrada establecida para el ensayo a través del intercambiador de calor IC2. El circuito de enfriamiento utiliza agua a 10 °C que proviene de una enfriadora eléctrica de agua (no mostrada en el diagrama) y el circuito de calentamiento utiliza agua caliente que proviene del termotanque TTE1. Los caudales en los tres circuitos son regulados mediante válvulas reguladores electroneumáticas (V51, V52 y V53) accionadas por controladores PID (indicados en línea punteada). Todo el sistema de control se desarrolló localmente utilizando un PLC marca Siemens modelo S7-1200. El diagrama indica además la ubicación de los sensores de temperatura de agua (T1, T2, T3, T4 y T5), el sensor de temperatura ambiente (Tamb), los sensores de irradiancia global horizontal (GHI) y en plano inclinado (GTI), el sensor de irradiancia difusa horizontal (DHI), la medida de viento (V) y el forzador de viento (GV).
Figura 3. Esquema de instalación termo-hidráulica
Fuente: Elaboración propia
Para
medir la temperatura a la entrada y a la salida del colector (T1 y T2)
se utilizaron PT100 de 3 hilos con transmisores de 4-20 mA de la marca
Herten. Dichos sensores fueron calibrados bajo normativa ISO vigente en
el LES reportándose una incertidumbre estándar de 0.02 °C. La
temperatura ambiente (Tamb) se registró con un sensor PT1000 de 2 hilos
Honeywell también calibrado en el LES con una incertidumbre estándar de
0.02 °C. La medida de caudal (C1) se realizó con un caudalímetro
electromagnético de Endress & Hauser
con una incertidumbre estándar de 0.5 % de la medida. La velocidad del
viento paralelo al plano del colector (V) se midió con un anemómetro de
copa NGR con incertidumbre estándar de 0.25 m/s. La irradiancia global
en el plano del colector (GTI) se midió con un piranómetro Kipp & Zonen
CMP 10. La irradiancia directa Gb en el plano del colector se estimó de
la siguiente manera: por un lado, se midió la irradiancia global en
plano horizontal (GHI) con un piranómetro Kipp & Zonen CMP11 y la irradiancia difusa en plano horizontal (DHI) con un piranómetro Kipp & Zonen
CMP6 montado con una banda de sombra del mismo fabricante, luego, la
irradiancia en incidencia normal (DNI) se estimó por diferencia de
estas dos medidas, y, multiplicando la DNI por el coseno del ángulo de
incidencia en el plano del colector se obtuvo Gb. La medida de
irradiancia difusa en plano inclinado con la banda de sombra es
corregida con un factor provisto por el fabricante. Dichos piranómetros
son calibrados anualmente en el LES de acuerdo a la norma ISO 9747:1992
contra un estándar secundario Kipp & Zonen CMP22 que es enviado cada cinco años a calibrar contra la referencia radiométrica mundial en World Radiation Center en Davos, Suiza. Todas las medidas se registraron cada 10 segundos usando un adquisidor de datos Fischer Scientific DT85.
PROCEDIMIENTO
Condiciones requeridas
La Tabla 1 muestra las condiciones requeridas por cada metodología para cada variable, incluyendo la variabilidad permitida. Estas son las condiciones que deben cumplir las medidas registradas durante el ensayo para poder ser utilizadas para determinar los parámetros con cada metodología. Como se puede observar, la metodología SST presenta requerimientos más estrictos respecto a los valores requeridos y su variabilidad. Por otro lado, la metodología QDT requiere que se representen diferentes condiciones meteorológicas durante el ensayo, como se explica en la página 24.
PROCEDIMIENTO
Condiciones requeridas
La Tabla 1 muestra las condiciones requeridas por cada metodología para cada variable, incluyendo la variabilidad permitida. Estas son las condiciones que deben cumplir las medidas registradas durante el ensayo para poder ser utilizadas para determinar los parámetros con cada metodología. Como se puede observar, la metodología SST presenta requerimientos más estrictos respecto a los valores requeridos y su variabilidad. Por otro lado, la metodología QDT requiere que se representen diferentes condiciones meteorológicas durante el ensayo, como se explica en la página 24.
Tabla 1. Condiciones de ensayo y variabilidad permitida
Fuente: ISO 9806:2017
Ensayo en Estado Estacionario SST
Ensayo de Rendimiento: Para este ensayo se midió la energía útil por unidad de área para cuatro temperaturas de entrada diferentes, y la temperatura de ensayo más alta fue 85 ⁰C. A su vez, por cada temperatura de entrada, se obtuvieron cuatro puntos independientes dando así un total de 16 puntos. Cada punto consistió en un promedio sobre un período de medición 10 minutos. Para verificar el estado estacionario, para cada punto, se tomaron promedios en períodos sucesivos de 30 segundos y se verificaron que las variaciones de dichos promedios respecto al promedio 10 minutal se encontraran dentro de los límites dados por la Tabla 1. Finalmente, los parámetros η0,hem, a1 y a2 se determinaron por regresión lineal usando el siguiente modelo:
Ensayo de Rendimiento: Para este ensayo se midió la energía útil por unidad de área para cuatro temperaturas de entrada diferentes, y la temperatura de ensayo más alta fue 85 ⁰C. A su vez, por cada temperatura de entrada, se obtuvieron cuatro puntos independientes dando así un total de 16 puntos. Cada punto consistió en un promedio sobre un período de medición 10 minutos. Para verificar el estado estacionario, para cada punto, se tomaron promedios en períodos sucesivos de 30 segundos y se verificaron que las variaciones de dichos promedios respecto al promedio 10 minutal se encontraran dentro de los límites dados por la Tabla 1. Finalmente, los parámetros η0,hem, a1 y a2 se determinaron por regresión lineal usando el siguiente modelo:
En donde Q́ u’’ es la energía útil del colector por unidad de área bruta y se calcula según:
Se
definió como variable dependiente para la regresión lineal a la energía
útil por unidad de área y como variables independientes G, (Tm - Ta) y (Tm - Ta)2; los coeficientes a determinar son η0,hem, a1 y a2.
Para asegurar la condición de incidencia normal se realizó un
seguimiento acimutal durante el ensayo. La Ecuación (7) resulta de
imponer la condición de estado estacionario (dTm/dt = 0) y la condición de incidencia normal (Khem = 1) en la Ecuación (1).
Ensayo de Modificador por Ángulo de Incidencia: En este caso el acimut del seguidor se fijó en 0° (orientación norte) y se tomaron distintas medidas de energía útil para diferentes ángulos de incidencia entre 40° y 70° antes y después del mediodía solar. El procesamiento de datos se realizó de igual forma que en el ensayo de rendimiento. Se impuso, para cada periodo de medición, una temperatura de entrada tal que la temperatura media del fluido fuera lo más próxima posible a la temperatura ambiente (para disminuir las pérdidas térmicas). El modificador por ángulo de incidencia se determinó usando:
Ensayo de Modificador por Ángulo de Incidencia: En este caso el acimut del seguidor se fijó en 0° (orientación norte) y se tomaron distintas medidas de energía útil para diferentes ángulos de incidencia entre 40° y 70° antes y después del mediodía solar. El procesamiento de datos se realizó de igual forma que en el ensayo de rendimiento. Se impuso, para cada periodo de medición, una temperatura de entrada tal que la temperatura media del fluido fuera lo más próxima posible a la temperatura ambiente (para disminuir las pérdidas térmicas). El modificador por ángulo de incidencia se determinó usando:
La Ecuación (9) resulta de despejar de la Ecuación (1) el modificador por ángulo de incidencia Khem e imponer el estado estacionario (dTm/dt = 0). Finalmente, el parámetro b0 de la Ecuación (4) se determinó por regresión lineal usando como variable dependiente (Khem – 1) y como variable independiente (1/cos(θ) – 1).
Ensayo de Capacidad Térmica: Este ensayo se hizo de acuerdo a la sección 25.2 de la norma ISO 9806:2017 pero teniendo en consideración el ajuste de segundo orden en las pérdidas térmicas, es decir el coeficiente a2. Al inicio del ensayo se fijó una temperatura de entrada igual a la temperatura ambiente y luego se cubrió el colector con una manta reflectiva dejando que alcance el estado estacionario. Luego, la cubierta se removió y se dejó que el colector alcance un nuevo punto en estado estacionario, diferente al primero por el efecto de la irradiancia solar. La capacidad térmica efectiva se determinó integrado la Ecuación (1) sobre el periodo de tiempo entre los dos puntos de funcionamiento estacionarios, como se muestra en la Ecuación (10). Para esto se utilizaron los datos promediados cada 30 segundos.
Ensayo de Capacidad Térmica: Este ensayo se hizo de acuerdo a la sección 25.2 de la norma ISO 9806:2017 pero teniendo en consideración el ajuste de segundo orden en las pérdidas térmicas, es decir el coeficiente a2. Al inicio del ensayo se fijó una temperatura de entrada igual a la temperatura ambiente y luego se cubrió el colector con una manta reflectiva dejando que alcance el estado estacionario. Luego, la cubierta se removió y se dejó que el colector alcance un nuevo punto en estado estacionario, diferente al primero por el efecto de la irradiancia solar. La capacidad térmica efectiva se determinó integrado la Ecuación (1) sobre el periodo de tiempo entre los dos puntos de funcionamiento estacionarios, como se muestra en la Ecuación (10). Para esto se utilizaron los datos promediados cada 30 segundos.
En este caso, también se realizó un seguimiento acimutal y por este motivo en la Ecuación (10) se asume incidencia normal (Khem = 1).
Ensayo en estado Cuasi-Dinámico QDT
Para este ensayo se midió la energía útil por unidad de área para cuatro temperaturas de entrada diferentes. Los datos fueron promediados cada 5 minutos y para cada temperatura el periodo de medición fue de al menos 3 horas. Se verificó además que durante cada período de medición se cumplieran las condiciones de la Tabla 1, comparando los promedios cada 30 segundos con los promedios de todo el período de medición. Para la temperatura más baja se tomó una serie de datos realizando seguimiento acimutal, de modo de obtener medidas con bajo ángulo de incidencia, y otra serie de datos fijando el azimut del seguidor en 0° (orientación norte) para obtener ángulos de incidencia altos. Ambas series se tomaron bajo condiciones de cielo claro. Las temperaturas restantes se realizaron bajo distintas condiciones de cielo realizando seguimiento acimutal para trabajar con ángulos de incidencia pequeños (Khem = 1). Si bien esto no es requisito de la norma, permite un mayor desacople de las variables dependientes y mejora la identificación de parámetros. Finalmente, las constantes η0,b, a1, a2, b0, Kd y C se determinaron todas a la vez por regresión lineal utilizando el modelo siguiente:
Ensayo en estado Cuasi-Dinámico QDT
Para este ensayo se midió la energía útil por unidad de área para cuatro temperaturas de entrada diferentes. Los datos fueron promediados cada 5 minutos y para cada temperatura el periodo de medición fue de al menos 3 horas. Se verificó además que durante cada período de medición se cumplieran las condiciones de la Tabla 1, comparando los promedios cada 30 segundos con los promedios de todo el período de medición. Para la temperatura más baja se tomó una serie de datos realizando seguimiento acimutal, de modo de obtener medidas con bajo ángulo de incidencia, y otra serie de datos fijando el azimut del seguidor en 0° (orientación norte) para obtener ángulos de incidencia altos. Ambas series se tomaron bajo condiciones de cielo claro. Las temperaturas restantes se realizaron bajo distintas condiciones de cielo realizando seguimiento acimutal para trabajar con ángulos de incidencia pequeños (Khem = 1). Si bien esto no es requisito de la norma, permite un mayor desacople de las variables dependientes y mejora la identificación de parámetros. Finalmente, las constantes η0,b, a1, a2, b0, Kd y C se determinaron todas a la vez por regresión lineal utilizando el modelo siguiente:
En donde la derivada temporal de la temperatura media del fluido se estimó por diferencias hacia adelante según:
En
donde Δt = 300 s (correspondiente a 5 minutos). La Ecuación (11) surge
de remplazar la Ecuación (4) en la Ecuación (3) y reordenar algunos
términos. La variable dependiente en este caso es la energía útil por
unidad de área y las variables independientes Gb, (1/cos(θ)-1) Gb, Gd, (Tm-Ta), (Tm-Ta)2 y dTm/dt; los coeficientes a determinar son η0,b, η0,bb0, η0,bKd, a1, a2 y C/Ag.
RESULTADOS
Comparación de métodos de ensayo
En la Tabla 2 se presentan los resultados de los ensayos para ambos métodos. En dicha tabla se presentan los coeficientes de los modelos térmicos de las Ecuaciones (1) y (3), sus respectivas incertidumbres y el estadístico t para cada parámetro (cociente entre el valor de un parámetro y su incertidumbre). No se reportó la incertidumbre de los parámetros η0,b y Kd para el caso SST por tratarse de magnitudes derivadas (Ecuaciones (5) y (6)).
RESULTADOS
Comparación de métodos de ensayo
En la Tabla 2 se presentan los resultados de los ensayos para ambos métodos. En dicha tabla se presentan los coeficientes de los modelos térmicos de las Ecuaciones (1) y (3), sus respectivas incertidumbres y el estadístico t para cada parámetro (cociente entre el valor de un parámetro y su incertidumbre). No se reportó la incertidumbre de los parámetros η0,b y Kd para el caso SST por tratarse de magnitudes derivadas (Ecuaciones (5) y (6)).
Tabla 2. Resultados
Fuente: Elaboración propia
En
todos los casos se obtuvo un estadístico t mayor 3, que indica que los
parámetros son estadísticamente significativos. Los parámetros η0,b y b0
difieren en menos de 1% (cociente entre diferencia de los valores y el
valor promedio) lo cual da cuenta de una buena concordancia entre estos
parámetros. En el caso de los coeficientes a1 y a2 la diferencia es mayor, 3% y 18% respectivamente. Ahora, si se presta atención a la Tabla 2 se puede observar que el término a1 en el caso SST es mayor que en el QDT y el término a2
en el caso SST es menor que en el caso QDT, es decir, que la diferencia
tiende a compensarse. El signo negativo de la covarianza entre estos
dos parámetros, obtenido en los dos casos, confirma esta relación
inversa. Este comportamiento no es nuevo; lo mismo se observa en las
referencias (Fischer et al., 2004; Osório y Carvalho, 2014). Una forma
más adecuada de comparar los coeficientes de pérdidas es utilizando un
coeficiente de pérdida linealizado a una diferencia de temperatura dada
ΔT, eso es:
Para una diferencia de 50 K se tienen a(50 K) = 4.585 W/m2.K para el método SST y a(50 K) = 4.553 W/m2.K
para el caso QDT, lo que da una diferencia menor al 1%. Las mayores
diferencias se presentan en el modificador de ángulo para la
irradiación difusa Kd y para la capacidad térmica por unidad de área C/Ag, 5% y 12% respectivamente. Entendemos que la diferencia en Kd
puede deberse en parte al factor de corrección empleado para corregir
la media de difusa en plano horizontal por el efecto de la banda de
sombra. Este factor supone que la difusa es isotrópica, por lo tanto,
este factor se puede mejorar utilizando factores de corrección que
contemplen la anisotropía de la atmósfera (Batlle et al., 1995). Otro
causante de esta diferencia puede ser que las horas bajo condiciones de
cielo parcialmente y totalmente cubierto durante el ensayo realizado
estén sub-representadas, dado que una adecuada representación de estas
condiciones es importante para la determinación de este parámetro bajo
la metodología QDT. La diferencia en la capacidad térmica efectiva
puede deberse a dos motivos: (i) puede ser que la cantidad de horas en
condiciones de cielo variable (variable dTm/dt alta en valor
absoluto) estén sub-representadas, y (ii) que el intervalo de tiempo 5
minutos para promediar los datos sea demasiado grande y por este motivo
tal vez no se capturó correctamente el comportamiento dinámico del
colector. Todos estos aspectos son motivo de nuestro trabajo actual, en
particular, en lo que refiere a los factores de corrección de la
radiación difusa.
Finalmente, para estudiar el efecto combinado de las diferencias entre los parámetros, se calculó para cada caso la energía útil producida por el colector utilizando la Ecuación (11), asumiendo incidencia normal y estado estacionario, para distintas condiciones de temperatura y de cielo. Para las distintas condiciones de cielo se utilizaron las condiciones de reporte estándar (Standard Reporting Conditions, SRC) dadas por la norma ISO 9806:2017. Los resultados se muestran en la Tabla 3.
Finalmente, para estudiar el efecto combinado de las diferencias entre los parámetros, se calculó para cada caso la energía útil producida por el colector utilizando la Ecuación (11), asumiendo incidencia normal y estado estacionario, para distintas condiciones de temperatura y de cielo. Para las distintas condiciones de cielo se utilizaron las condiciones de reporte estándar (Standard Reporting Conditions, SRC) dadas por la norma ISO 9806:2017. Los resultados se muestran en la Tabla 3.
Tabla 3. Potencia útil (W) producida por el colector bajo distintas condiciones de cielo y temperatura
Fuente: Elaboración propia
Para
condiciones de cielo claro la diferencia en la energía útil es poco
significativa (menor al 1%). En condiciones de cielo parcialmente
nublado la diferencia comienza crecer a medida que crece la diferencia
de temperatura, si bien sigue siendo poco significativa. En condiciones
de cielo cubierto la diferencia comienza a ser significativa alcanzando
el 15 % a una diferencia de temperatura de 40 K. Entendemos que esa
diferencia se debe principalmente a la diferencia en el modificador por
ángulo de incidencia para la irradiancia difusa Kd, lo que es motivo actual de estudio.
Estimación de cantidad de ensayos anuales
Para mostrar la aplicabilidad del método QDT en el clima de la región, en la Tabla 4 se muestra un estimativo de la cantidad de ensayos posibles por mes en la localidad de Salto (Uruguay) tanto para seguir el método SST realizando un seguimiento acimutal como para seguir el método QDT con una orientación norte, en ambos casos, con una inclinación horizontal fija de 45°. Una hora se considera válida para un tipo de ensayo si cumple las condiciones presentadas en la Tabla 1. Para el caso SST se verificó además que las horas fueran en condiciones de cielo claro. El análisis se realizó utilizando el año meteorológico típico de la localidad de Salto (latitud = 31.28° S, longitud = 57.92° O, altitud = 56 m sobre el nivel del mar), Uruguay, a escala horaria (Alonso-Suárez et al., 2019). Se supuso una duración de los ensayos de 10 horas para el caso SST y de 15 horas para el caso QDT, tiempo estimado en el cual se logra medir las variables necesarias para cada ensayo si se cumplen las condiciones de la Tabla 1. Además, para el caso QDT se impuso que al menos un 30% de las horas fueran con cielo claro, debido a que es importante para este ensayo capturar cierta dispersión en las condiciones meteorológicas, y las condiciones de cielo despejado no deben ser sub-representadas. Para verificar esta condición se impuso que el índice de cielo claro fuera mayor a 0.8 y que la irradiación horaria en incidencia normal (DNI) fuera mayor a 200 Wh/m2. Para el cálculo del índice de cielo claro se utilizó la expresión dada por Pérez et al. (1990):
Estimación de cantidad de ensayos anuales
Para mostrar la aplicabilidad del método QDT en el clima de la región, en la Tabla 4 se muestra un estimativo de la cantidad de ensayos posibles por mes en la localidad de Salto (Uruguay) tanto para seguir el método SST realizando un seguimiento acimutal como para seguir el método QDT con una orientación norte, en ambos casos, con una inclinación horizontal fija de 45°. Una hora se considera válida para un tipo de ensayo si cumple las condiciones presentadas en la Tabla 1. Para el caso SST se verificó además que las horas fueran en condiciones de cielo claro. El análisis se realizó utilizando el año meteorológico típico de la localidad de Salto (latitud = 31.28° S, longitud = 57.92° O, altitud = 56 m sobre el nivel del mar), Uruguay, a escala horaria (Alonso-Suárez et al., 2019). Se supuso una duración de los ensayos de 10 horas para el caso SST y de 15 horas para el caso QDT, tiempo estimado en el cual se logra medir las variables necesarias para cada ensayo si se cumplen las condiciones de la Tabla 1. Además, para el caso QDT se impuso que al menos un 30% de las horas fueran con cielo claro, debido a que es importante para este ensayo capturar cierta dispersión en las condiciones meteorológicas, y las condiciones de cielo despejado no deben ser sub-representadas. Para verificar esta condición se impuso que el índice de cielo claro fuera mayor a 0.8 y que la irradiación horaria en incidencia normal (DNI) fuera mayor a 200 Wh/m2. Para el cálculo del índice de cielo claro se utilizó la expresión dada por Pérez et al. (1990):
En donde kt
es el índice de claridad horario usual (cociente entre la irradiación
horaria global en plano horizontal y la correspondiente irradiación en
un plano horizontal en el tope de la atmósfera) y m es la masa de aire (Duffie y Beckman, 2013).
Tabla 4. Comparativa de cantidad de ensayos posibles por mes
Fuente: Elaboración propia
En
la Tabla 4 se puede ver que la cantidad de ensayos posibles en un año
en el caso QDT es más del doble que la cantidad de ensayos en el caso
SST. Esta proporción coincide con lo estimado en otras localidades
(Rojas, 2008; Osório y Carvalho, 2014).
CONCLUSIONES
Se estudió la aplicabilidad del método cuasi dinámico en Uruguay y se determinó la cantidad de ensayos a nivel anual que podrían llegar a realizarse si se siguiera este método. Esta cantidad es mayor al doble de la cantidad de ensayos que podrían obtenerse con la metodología tradicional en estado estacionario, lo cual muestra la ventaja de la metodología QDT en Uruguay y en la región climática homogénea de la Pampa Húmeda del sureste de Sudamérica.
Por otra parte, se concluye que el método QDT se implementó correctamente para el colector estudiado y se puede utilizar en el clima de la región, en tanto, su comparación con los ensayos SST está dentro de lo esperado y reportado en otras publicaciones. La principal diferencia entre los parámetros se observó en Kd y C/Ag. Teniendo en cuenta que el parámetro Kd afecta la producción de energía útil del colector en condiciones de cielo parcialmente nublado y cubierto, y que además una porción significativa de las horas en el año en nuestro país corresponde a este tipo de condiciones (Alonso-Suárez et al., 2019), consideramos importante continuar estudiando los motivos de esta diferencia. La diferencia en el parámetro C/Ag es otro punto a estudiar a futuro. Si bien en muchos programas de simulación de colectores este término no es tenido en cuenta, restándole importancia (asumen estado estacionario), en el ensayo QDT este término es determinado de forma simultánea con los otros parámetros y, por lo tanto, una mala determinación de éste puede afectar negativamente a los otros.
El presente trabajo creemos es un paso significativo en el desarrollo de capacidades de ensayo de colectores solares en la región. No obstante, entendemos necesario continuar ensayando colectores de placa plana y en especial otros tipos de colectores, como los de tubos al vacío, para poner a punto la metodología QDT. Esperamos que el presente trabajo impulse el uso de la metodología de ensayos QDT no sólo en nuestra región, sino en el resto de Latinoamérica y que ésta pueda ser validada para su uso en otros climas o países.
CONCLUSIONES
Se estudió la aplicabilidad del método cuasi dinámico en Uruguay y se determinó la cantidad de ensayos a nivel anual que podrían llegar a realizarse si se siguiera este método. Esta cantidad es mayor al doble de la cantidad de ensayos que podrían obtenerse con la metodología tradicional en estado estacionario, lo cual muestra la ventaja de la metodología QDT en Uruguay y en la región climática homogénea de la Pampa Húmeda del sureste de Sudamérica.
Por otra parte, se concluye que el método QDT se implementó correctamente para el colector estudiado y se puede utilizar en el clima de la región, en tanto, su comparación con los ensayos SST está dentro de lo esperado y reportado en otras publicaciones. La principal diferencia entre los parámetros se observó en Kd y C/Ag. Teniendo en cuenta que el parámetro Kd afecta la producción de energía útil del colector en condiciones de cielo parcialmente nublado y cubierto, y que además una porción significativa de las horas en el año en nuestro país corresponde a este tipo de condiciones (Alonso-Suárez et al., 2019), consideramos importante continuar estudiando los motivos de esta diferencia. La diferencia en el parámetro C/Ag es otro punto a estudiar a futuro. Si bien en muchos programas de simulación de colectores este término no es tenido en cuenta, restándole importancia (asumen estado estacionario), en el ensayo QDT este término es determinado de forma simultánea con los otros parámetros y, por lo tanto, una mala determinación de éste puede afectar negativamente a los otros.
El presente trabajo creemos es un paso significativo en el desarrollo de capacidades de ensayo de colectores solares en la región. No obstante, entendemos necesario continuar ensayando colectores de placa plana y en especial otros tipos de colectores, como los de tubos al vacío, para poner a punto la metodología QDT. Esperamos que el presente trabajo impulse el uso de la metodología de ensayos QDT no sólo en nuestra región, sino en el resto de Latinoamérica y que ésta pueda ser validada para su uso en otros climas o países.
AGRADECIMIENTOS
Los autores desean agradecer al Ministerio de Industria, Energía y Minería (MIEM, Uruguay), en especial a su Dirección Nacional de Energía (DNE), al Fideicomiso Uruguayo de Ahorro y Eficiencia Energética (Fudaee, Uruguay) y a la Corporación Nacional para el Desarrollo (CND, Uruguay), por haber proporcionado apoyo financiero y logístico para el desarrollo del BECS y haber impulsado este proyecto con capacidades nacionales. Los autores agradecen también al PTB de Alemania por impulsar y financiar el inter-laboratorio de ensayos de eficiencia de colectores solares, lo que nos ha brindado certezas técnicas sobre la capacidad de ensayo local.
Los autores desean agradecer al Ministerio de Industria, Energía y Minería (MIEM, Uruguay), en especial a su Dirección Nacional de Energía (DNE), al Fideicomiso Uruguayo de Ahorro y Eficiencia Energética (Fudaee, Uruguay) y a la Corporación Nacional para el Desarrollo (CND, Uruguay), por haber proporcionado apoyo financiero y logístico para el desarrollo del BECS y haber impulsado este proyecto con capacidades nacionales. Los autores agradecen también al PTB de Alemania por impulsar y financiar el inter-laboratorio de ensayos de eficiencia de colectores solares, lo que nos ha brindado certezas técnicas sobre la capacidad de ensayo local.
NOMENCLATURA
Ag área bruta del colector (m2)
a factor de pérdidas térmicas linealizado (W/(m2.K))
a1 factor de pérdidas térmicas de primer orden (W/(m2.K))
a2 factor de pérdidas térmicas de segundo orden (W/(m2.K2))
b0 factor de ajuste del modificador de ángulo (adimensionado)
C capacidad térmica del colector (J/K)
Cp calor específico del fluido a la temperatura Tm (J/Kg.K)
fd fracción difusa en el plano del colector (adimensionado)
G irradiancia solar global en el plano del colector (W/m2)
Gb irradiancia solar directa en el plano del colector (W/m2)
Gd irradiancia solar difusa en el plano del colector (W/m2)
kt índice de claridad
kt' índice de cielo claro
Khem modificador por ángulo de incidencia de la irradiancia solar global (adimensionado)
Kb modificador por ángulo de incidencia de la irradiancia solar directa (adimensionado)
Kd modificador por ángulo de incidencia de la irradiancia solar difusa (adimensionado)
m masa de aire
ḿ caudal másico de fluido que circula a través del colector (kg/s)
Q́u'' energía útil producida por el colector por unidad de área (W/m2)
Ti temperatura del fluido a la entrada del colector (°C)
To temperatura del fluido a la salida del colector (°C)
Tm temperatura media del fluido en el colector (promedio de la entrada y la salida) (°C)
Ta temperatura ambiente (°C)
η0,hem rendimiento máximo del colector referido a irradiancia solar global (adimensionado)
η0,b rendimiento máximo del colector referido a irradiancia solar directa (adimensionado)
θ ángulo de incidencia en el plano del colector (°)
Ag área bruta del colector (m2)
a factor de pérdidas térmicas linealizado (W/(m2.K))
a1 factor de pérdidas térmicas de primer orden (W/(m2.K))
a2 factor de pérdidas térmicas de segundo orden (W/(m2.K2))
b0 factor de ajuste del modificador de ángulo (adimensionado)
C capacidad térmica del colector (J/K)
Cp calor específico del fluido a la temperatura Tm (J/Kg.K)
fd fracción difusa en el plano del colector (adimensionado)
G irradiancia solar global en el plano del colector (W/m2)
Gb irradiancia solar directa en el plano del colector (W/m2)
Gd irradiancia solar difusa en el plano del colector (W/m2)
kt índice de claridad
kt' índice de cielo claro
Khem modificador por ángulo de incidencia de la irradiancia solar global (adimensionado)
Kb modificador por ángulo de incidencia de la irradiancia solar directa (adimensionado)
Kd modificador por ángulo de incidencia de la irradiancia solar difusa (adimensionado)
m masa de aire
ḿ caudal másico de fluido que circula a través del colector (kg/s)
Q́u'' energía útil producida por el colector por unidad de área (W/m2)
Ti temperatura del fluido a la entrada del colector (°C)
To temperatura del fluido a la salida del colector (°C)
Tm temperatura media del fluido en el colector (promedio de la entrada y la salida) (°C)
Ta temperatura ambiente (°C)
η0,hem rendimiento máximo del colector referido a irradiancia solar global (adimensionado)
η0,b rendimiento máximo del colector referido a irradiancia solar directa (adimensionado)
θ ángulo de incidencia en el plano del colector (°)
REFERENCIAS
Alonso-Suárez R., Bidegain M., Abal G., Modernell P. (2019). Año meteorológico típico para aplicaciones de energía solar (AMTUes): series horarias típicas para 5 sitios del Uruguay. Ediciones Universitarias (UCUR) de la Universidad de la República (Udelar), ISBN 978-9974-0-1674-7.
Batlles, F.J., Olmo, F.J. y Alados-Armboledas, L. (1995). On shadowband correction methods for diffuse irradiance measurements. Solar Energy, 54(2): 105-114.
Duffie, J.A. y Beckman, W.A. (2013). Solar Engineering of Thermal Processes. Wiley and Sons, Hoboken, New Jersey, third edition, 2006.
Fischer, S., Heidemann, W., Müller-Steinhagen, H., Perers, B., Bergquist, P., Hellström, B. (2004). Collector test method under quasi-dynamic conditions according to the European Standard EN 12975-2. Solar Energy, 76(1-3):117-123.
Fisher, S. (2020). Expert report - quality infrastructure for energy efficiency and renewable energy in Latin America and the Caribbean. + Report #95309. Recuperado de http://les.edu.uy/papers/Expert_Report-RR-Test_R1.pdf
García de Jalón, A., Sallaberry, F., Olano, X., Mateu, E., Astiz R., Ezcurra, M., Ramírez, L. (2011). Comparison of thermal efficiency curves of solar collectors tested in outdoor conditions. Proceedings of ISES Solar World Congress, 28 august - 2 september, Kassel, Germany.
International Energy Agency (IEA). (2016a). Final energy consumption in the world. Recuperado de https://www.iea. org/Sankey/#?c=World&s=Final%20consumption
International Energy Agency (IEA). (2016b). Final energy consumption in Uruguay. Recuperado de https://www.iea. org/Sankey/#?c=Uruguay&s=Final20consumption
ISO 9806. (2017). Solar Energy – Solar thermal collectors – Test methods. Switzerland: International Organization of Standarization.
Kovács, P., Petterson, U., Persson, M., Perers, B., Fischer, S. (2011). Improving the accurancy in performance prediction for new collector deings. Proceedings of ISES Solar World Congress, 28 august - 2 september, Kassel, Germany.
Kratzenberg, M.G., Beyer, H.G., Colle, S. (2006). Uncertainty calculation applied to different regression methods in the quasi-dynamic collector test. Solar Energy, 80, 1453-1462.
Muschaweck, J. y Spirkl, W. (1993). Dynamic solar collector performance testing. Solar Energy Materials and Solar Cells, 30(2):95-105.
Osório, T. y Carvalho, M. (2014). Testing of solar thermal collectors under transient conditions. Solar Energy, 104, 71-81.
Peel, M.C., Finlayson, B.L. y Mcmahon, T.A. (2007). Update world mapo of the köppen-geiger climate classification. Hydrology and Earth System Sciences Discussions, 11, 1633- 1644.
Perers, B. (1993). Dynamic method for solar collector array testing and evaluation with standard database and simulation programas. Solar Energy, 50(6), 517-526.
Pérez, R., Ineichen, P., Seals, R. y Zelenka, A. (1990). Making full use of the clearness index for parameterizing hourly insolation conditions. Solar Energy 45(2), 111-114.
Renewable Energy Policy Newtwork for the 21st Century (REN21). (2018). Renewables 2018 global status report. Recuperado de http://www.ren21.net/gsr-2018/pages/ foreword/foreword/
Rojas, D., Beermann, S.A., Klein, S.A. y Reindl, D.T. (2008). Thermal performance testing of flat-plate collectors. Solar Energy, 82, 746-757.
Zambolin, E. y Del Col D. (2010). Experimental analysis of thermal performance of flat plate and evacuated tube solar collectors in stationary standard and daily conditions. Solar Energy, 84, 1382-1396.
Alonso-Suárez R., Bidegain M., Abal G., Modernell P. (2019). Año meteorológico típico para aplicaciones de energía solar (AMTUes): series horarias típicas para 5 sitios del Uruguay. Ediciones Universitarias (UCUR) de la Universidad de la República (Udelar), ISBN 978-9974-0-1674-7.
Batlles, F.J., Olmo, F.J. y Alados-Armboledas, L. (1995). On shadowband correction methods for diffuse irradiance measurements. Solar Energy, 54(2): 105-114.
Duffie, J.A. y Beckman, W.A. (2013). Solar Engineering of Thermal Processes. Wiley and Sons, Hoboken, New Jersey, third edition, 2006.
Fischer, S., Heidemann, W., Müller-Steinhagen, H., Perers, B., Bergquist, P., Hellström, B. (2004). Collector test method under quasi-dynamic conditions according to the European Standard EN 12975-2. Solar Energy, 76(1-3):117-123.
Fisher, S. (2020). Expert report - quality infrastructure for energy efficiency and renewable energy in Latin America and the Caribbean. + Report #95309. Recuperado de http://les.edu.uy/papers/Expert_Report-RR-Test_R1.pdf
García de Jalón, A., Sallaberry, F., Olano, X., Mateu, E., Astiz R., Ezcurra, M., Ramírez, L. (2011). Comparison of thermal efficiency curves of solar collectors tested in outdoor conditions. Proceedings of ISES Solar World Congress, 28 august - 2 september, Kassel, Germany.
International Energy Agency (IEA). (2016a). Final energy consumption in the world. Recuperado de https://www.iea. org/Sankey/#?c=World&s=Final%20consumption
International Energy Agency (IEA). (2016b). Final energy consumption in Uruguay. Recuperado de https://www.iea. org/Sankey/#?c=Uruguay&s=Final20consumption
ISO 9806. (2017). Solar Energy – Solar thermal collectors – Test methods. Switzerland: International Organization of Standarization.
Kovács, P., Petterson, U., Persson, M., Perers, B., Fischer, S. (2011). Improving the accurancy in performance prediction for new collector deings. Proceedings of ISES Solar World Congress, 28 august - 2 september, Kassel, Germany.
Kratzenberg, M.G., Beyer, H.G., Colle, S. (2006). Uncertainty calculation applied to different regression methods in the quasi-dynamic collector test. Solar Energy, 80, 1453-1462.
Muschaweck, J. y Spirkl, W. (1993). Dynamic solar collector performance testing. Solar Energy Materials and Solar Cells, 30(2):95-105.
Osório, T. y Carvalho, M. (2014). Testing of solar thermal collectors under transient conditions. Solar Energy, 104, 71-81.
Peel, M.C., Finlayson, B.L. y Mcmahon, T.A. (2007). Update world mapo of the köppen-geiger climate classification. Hydrology and Earth System Sciences Discussions, 11, 1633- 1644.
Perers, B. (1993). Dynamic method for solar collector array testing and evaluation with standard database and simulation programas. Solar Energy, 50(6), 517-526.
Pérez, R., Ineichen, P., Seals, R. y Zelenka, A. (1990). Making full use of the clearness index for parameterizing hourly insolation conditions. Solar Energy 45(2), 111-114.
Renewable Energy Policy Newtwork for the 21st Century (REN21). (2018). Renewables 2018 global status report. Recuperado de http://www.ren21.net/gsr-2018/pages/ foreword/foreword/
Rojas, D., Beermann, S.A., Klein, S.A. y Reindl, D.T. (2008). Thermal performance testing of flat-plate collectors. Solar Energy, 82, 746-757.
Zambolin, E. y Del Col D. (2010). Experimental analysis of thermal performance of flat plate and evacuated tube solar collectors in stationary standard and daily conditions. Solar Energy, 84, 1382-1396.