GEOSOL: UNA HERRAMIENTA COMPUTACIONAL PARA EL CLCULO DE COORDENADAS SOLARES Y LA ESTIMACIN DE IRRADIACIN SOLAR HORARIA Alejandro L.

Hernndez1 INENCO Instituto de Investigaciones en Energa No Convencional U.N.Sa. - CONICET Universidad Nacional de Salta Buenos Aires 177, 4400, Salta, Argentina FAX: 54-87-4255489, TE: 54-87-4255579 EMAIL: [email protected]
RESUMEN El conocimiento de la posicin del sol a lo largo del da es de suma importancia para todas las aplicaciones que aprovechan la energa solar. Por otra parte, si no se dispone de valores medidos de irradiacin solar es necesario contar con mtodos analticos que permitan estimarla. En este trabajo se presenta el programa GEOSOL V.2.0 para Windows que permite, a travs de sus clculos, cumplir con estos dos objetivos. Su salida grfica es til no slo desde el punto de vista prctico sino tambin educativo ya que representa la trayectoria del sol en 2D (carta solar) y 3D, posibilitando comparar las diferencias entre invierno y verano y entre distintos lugares del planeta. Posee tres mtodos de estimacin de irradiacin solar: de Page y de Hottel para das claros y de Liu-Jordan para das medios mensuales. Se presenta la validacin de los tres mtodos con mediciones realizadas en la Universidad Nacional de Salta, denotando un excelente acuerdo entre valores medidos y estimados. Palabras clave: radiacin solar, aplicaciones con energa solar. EL PROGRAMA GEOSOL V.2.0 GEOSOL es un programa de clculo y graficacin creado por el autor del presente trabajo. Corre bajo Windows y permite obtener, para cualquier lugar y da del ao, los siguientes valores numricos (de J. Duffie & W. Beckman, 1991): Horas solares de salida y puesta del sol sobre el horizonte. Duracin del da. Diferencia horaria entre Hora Solar y Hora Oficial. Declinacin solar para la fecha consignada. Angulo horario, altitud y azimut solares, hora por hora, desde la salida hasta la puesta del sol. Irradiacin solar directa, difusa y total en MJ/m2, hora por hora, sobre cualquier superficie especificada por el usuario.

GEOSOL genera, adems, un conjunto de grficos con los cuales se representa: La evolucin horaria de las irradiaciones solares directa, difusa y total mediante un grfico de barras. La trayectoria del sol sobre la bveda celeste en 2D (Carta Solar) y 3D (perspectiva isomtrica). Grfica de obstculos a la trayectoria del sol en 2D. La pantalla principal de GEOSOL est dividida en 5 campos segn se observa en la figura 1. En el de la derecha, de color gris, y sobre su sector superior se encuentra la ventana para el ingreso de Datos del Lugar: latitud, longitud, uso horario en horas respecto del tiempo medio de Greenwich (GMT), altitud sobre el nivel del mar, albedo del entorno y da de clculo. Si se desconoce el uso horario, presionando el botn Uso horario (GMT) se despliega un planisferio donde se observa que a la Argentina le corresponde el valor 3 (tres horas menos que la correspondiente a Greenwich). Al presionar el botn Albedo se despliega una pantalla con los valores de albedo correspondientes a los tipos de coberturas de suelo ms comunes. El n-simo da del ao se calcula ingresando la fecha correspondiente en la pantalla que se despliega al presionar sobre el botn Da de clculo. En la ventana denominada Datos del plano se deben ingresar la pendiente y el ngulo azimutal (en grados) de la superficie sobre la cual se ha de estimar la irradiacin solar. La convencin de signos empleada en GEOSOL es la que figura en la tabla 1.

Integrante del Proyecto Foncyt de la Agencia Nacional de Promocin Cientfica y Tcnica, FUNDALTES BID 802/OC AR PICT 2000 N 13-09991

Figura 1: pantalla principal de GEOSOL. Hemisferio Sur Norte Angulo horario solar 0 al Norte, positivo a la maana 0 al Sur, positivo a la maana Azimut solar 0 al Sur, positivo al Este 0 al Sur, positivo al Este Azimut del plano 0 al Sur, positivo al Este 0 al Sur, positivo al Este

Tabla 1: convencin de signos para los ngulos solares y el azimut del plano colector. Debajo de la ventana correspondiente a los datos del plano se encuentra la que permite seleccionar el mtodo de estimacin de irradiacin solar a emplear. GEOSOL incluye tres mtodos de clculo: 1) 2) 3) Mtodo de da claro de Page. Mtodo de da claro de Hotel. Mtodo de Liu-Jordan para da medio mensual.

El primero, aplicable solamente a los das sin nubosidad, fue desarrollado por el grupo del Prof. J. K. Page de la Universidad de Sheffield en la dcada del 70. Para su uso es necesario dar el contenido de humedad del aire (como presin parcial de vapor) y un parmetro ajustable que depende del lugar, la turbidez atmosfrica, la cual toma en cuenta la atenuacin de la radiacin solar debida al polvo y dems aerosoles presentes en la atmsfera. En cuanto al modelo de cielo considera una distribucin anisotrpica de la radiacin difusa dividindola en dos componentes: istropa de fondo y circumsolar la que, desde el punto de vista geomtrico, se trata como directa. El modelo anisotrpico conduce a valores ms realistas de radiacin difusa que el modelo isotrpico. Cuando se selecciona este mtodo se despliega la ventana que se muestra en la figura 2. Dado que la presin parcial de vapor puede calcularse a partir de la temperatura de bulbo seco y de la humedad relativa del aire, se escogieron estas variables como entradas al programa ya que, generalmente, es ms fcil contar con valores medidos de humedad y temperatura que de presin parcial de vapor en la zona donde se pretende estimar la radiacin solar o en zonas prximas a ella. Por otra parte, GEOSOL calcula la turbidez atmosfrica mediante una correlacin que depende solamente de la altitud del lugar. Este valor es el que aparece en la casilla titulada Turbidez atmosfrica y slo debe cambiarse si se conoce un valor apropiado cuando la atmsfera del lugar presenta altos contenidos de aerosoles debidos a fenmenos anmalos locales (incendios forestales, erupciones volcnicas, densas nieblas, etc.). O bien puede emplearse a la turbidez como parmetro de ajuste si se dispone de valores de radiacin solar medidos in situ.

Figura 2: variables del mtodo de Page.

El segundo mtodo utiliza para la estimacin de la irradiancia solar directa sobre una superficie normal al haz la ecuacin emprica de Hottel (1976), vlida hasta los 2.500 m.s.n.m con una visibilidad de 23 Km. Para la estimacin de la radiacin difusa, GEOSOL emplea tambin en este caso el modelo de cielo anisotrpico de Page. Cuando se selecciona este mtodo se despliega la ventana Tipo de Clima que se muestra en la figura 3. La ecuacin de Hottel incluye en su formulacin tres coeficientes cuyos valores dependen de la altitud del lugar y de cuatro posibles tipos de clima. En aquellas regiones del planeta en que no existe una marcada diferencia entre inviernos y veranos (latitudes entre el Ecuador y los Trpicos) puede seleccionarse el tipo de clima Tropical. Entre los trpicos y los polos la diferenciacin entre invierno y verano aumenta sustancialmente y deben emplearse las opciones Verano de Latitud media e Invierno de Latitud media segn la poca del ao. Para latitudes prximas a los crculos polares puede emplearse el tipo de clima Verano Subrtico. Sin embargo, debe recordarse que slo en el hemisferio Norte existen tierras habitadas a ambos lados del crculo polar. Evidentemente con este mtodo no hay tanta flexibilidad como en el anterior a la hora de intentar el ajuste de valores de radiacin medidos. Una vez seleccionado el tipo de clima se procede a efectuar los clculos, como en el caso anterior, clickeando sobre el botn Calcular.

Figura 3: climas para el mtodo de Hottel.

Por su parte, el mtodo conocido como de Liu Jordan (Liu B. y Jordan R., 1960) permite estimar las irradiaciones horarias directa y difusa sobre cualquier superficie a partir de la irradiacin solar global diaria media mensual sobre superficie horizontal. Este dato es el promedio de todos los valores diarios de irradiacin solar global registrados durante un mismo mes del ao, promediados adems a lo largo de varios aos. Su valor, nico para todo el mes, representa la cantidad de energa solar promedio que incide por cada m2 de superficie horizontal a lo largo de un da, se mide en MJoule/m2 y puede obtenerse consultado las tablas editadas por el Servicio Meteorolgico Nacional o los registros existentes en los aeropuertos ms cercanos a la localidad de inters. Dado que deriva de un promedio entre das que pueden ser claros, seminublados o nublados, su valor es ms conservativo que el correspondiente a los dos mtodos anteriores. Cuando se selecciona este mtodo se despliega la ventana que se muestra en la figura 4. Luego de ingresar el valor de la irradiacin solar diaria, GEOSOL permite seleccionar entre dos modelos de cielo para el tratamiento de la radiacin difusa: anisotrpico e isotrpico. En el primer caso utiliza el modelo de Reindl et al. (1990) basado en la hiptesis de que toda la radiacin difusa puede ser representada por tres componentes: istropa de fondo, circumsolar y difusa de horizonte. En el segundo caso toda la radiacin difusa proviene en forma istropa desde la bveda celeste. Sobre superficies inclinadas, el modelo anisotrpico conduce a valores de radiacin estimada sensiblemente superiores (alrededor del 7 %) respecto de los obtenidos con el modelo isotrpico. Figura 4: ventana desplegada al seleccionar el mtodo de Liu Jordan. Una vez seleccionado el mtodo de clculo, clickeando sobre el botn Calcular correspondiente se actualizan los otros cuatro campos de la pantalla principal y se generan dos archivos en cdigo ASCII: resultados.dat con los valores de las irradianciones horarias y angulos.dat con las coordenadas solares. Segn se observa en la figura 1, en el sector superior izquierdo de la pantalla figuran los valores de las irradiaciones horarias directa, difusa y total, hora por hora, desde la salida hasta la puesta del sol. Presenta, adems, el valor de irradiacin total diaria sobre la superficie escogida. En el sector superior central se grafica la evolucin temporal de los valores horarios donde, con color rojo, se representa a la componente directa y con turquesa a la difusa. El mximo valor de irradiacin solar disponible corresponde a su valor normal al haz fuera de la atmsfera (Radiacin Extraterrestre) y ha sido representado en color azul para su comparacin con los valores obtenidos a nivel del suelo. El sector inferior izquierdo presenta los valores de las horas solares de salida y puesta del sol, la duracin del da, la declinacin solar y los ngulos horario (ws), altitud y azimut solares, hora por hora. En el sector inferior derecho se encuentra una ventana con la diferencia entre la hora solar y la oficial para ese da. Finalmente, en el sector inferior central se representa, mediante un grfico 3D en perspectiva isomtrica, la trayectoria seguida por el sol a lo largo del da. Esta forma de visualizar el movimiento del sol es muy interesante desde el punto de vista educativo ya que permite comparar rpidamente las marcadas diferencias existentes entre el invierno y el verano y entre distintos puntos del planeta. En la figura 5 se representa la trayectoria del sol el da 17 de enero en 4 latitudes distintas. El caso (a) corresponde al invierno en el hemisferio Norte por lo que el sol se mueve sobre un plano inclinado al Sur mientras que (c) y (d) corresponden al verano en el hemisferio Sur por lo que el sol se mueve sobre un plano inclinado al Norte. Si se clickea con el mouse sobre el grfico 3D, se despliega la pantalla con la Carta Solar mostrada en la figura 6. Esta carta es una representacin 2D de la esfera celeste vista desde el cenit, en la cual se grafica la trayectoria del sol mediante sus coordenadas angulares de altitud (crculos de color azul) y azimut (radios de color rojo). El puntero del mouse se mueve exclusivamente sobre la curva de la trayectoria (color verde) a fin de facilitar la determinacin de la posicin del sol en cualquier instante del da, la que aparece en la ventana titulada Valores instantneos sobre el sector inferior izquierdo.

a) Lat.: 25 N

b) Lat.: 0 (Ecuador)

c) Lat.: 25 S

d) Lat.: 69 S

Figura 5: trayectoria del sol el da 17 de Enero en 4 latitudes distintas: 25 N, Ecuador, 25 S y 69 S.

Figura 6: pantalla con la carta solar y los valores de las coordenadas solares instantneos y hora por hora.

La carta solar puede imprimirse clickeando sobre el botn Imprimir Carta con lo cual se dispone de una til herramienta de campo cuando se deben realizar instalaciones solares en lugares remotos donde resulta imposible utilizar una PC. Presionando sobre el botn Dibujar Obstc. se despliega la pantalla del graficador de obstculos que se observa en la figura 7. Cuando se elige el lugar de emplazamiento de algn artefacto solar (colectores, cocinas, paneles fotovoltaicos, etc.) resulta conveniente saber si los objetos alrededor de l (casas, rboles, instalaciones, etc.) han de sombrearlo en alguna poca del ao. Si se miden los ngulos de azimut y altitud de cada vrtice de los obstculos presentes en el lugar y se los grafica juntamente con la trayectoria del sol, se puede analizar en qu momentos del da o del ao el sol queda apantallado. El graficador de obstculos de GEOSOL permite realizar esta tarea ya que, cuando se lo invoca, aparece el grfico 2D de la figura 7 slo con la trayectoria del sol para la fecha elegida y, con las herramientas de dibujo que contiene, se procede a dibujar los obstculos de la misma forma que se dibuja con el Paint de Windows. A fin de facilitar la tarea, en las ventanas ubicadas arriba a la derecha figuran los valores del azimut y la altitud del punto sealado por el puntero el mouse.

Figura 7: pantalla del graficador de obstculos en 2D con un ejemplo correspondiente al da 17 de Enero en Salta. Una vez finalizado el dibujo se lo puede guardar en un archivo con extensin .grf que slo contiene las figuras de color azul. De esta forma, se puede cambiar la fecha de clculo para ver qu ocurre en otra poca del ao pero el dibujo de los obstculos (que aparecer cuando se cargue el archivo .grf) seguir siendo el mismo. En el ejemplo de la figura se observa que, para el da 17 de Enero en Salta, el edificio ubicado al Sudeste apantalla el sol durante las 3 primeras horas de la maana, mientras que el rbol ubicado al Sudoeste lo hace entre las 15:30 y las 18:00, hora solar. VALIDACIN DE LOS VALORES DE IRRADIACIN TOTAL HORARIA ESTIMADOS CON GEOSOL En la figura 8 se comparan los resultados arrojados por los tres mtodos de estimacin de GEOSOL con valores medidos sobre superficie horizontal, un da 7 de Julio, en la Universidad Nacional de Salta. Debido a que no se efectuaron medidas diferenciadas de irradiancia directa y difusa, la comparacin corresponde a los valores globales. De acuerdo a los datos medidos se emple una temperatura media de 6 C y una humedad relativa media del 68% para el mtodo de Page, resultando una presin de vapor de 6,43 mb y una turbidez atmosfrica de 0,092 (correspondiente a 1.230 m.s.n.m.). Para le mtodo de Hottel se seleccion el tipo de clima Invierno de Latitud media y para el de Liu-Jordan se ingres una irradiacin global diaria sobre superficie horizontal de 15,9 MJ/m2, valor obtenido a partir de las mediciones. Segn se observa en la figura, el acuerdo entre los valores medidos y estimados por los tres mtodos es significativo. Se aprecia cierta asimetra en las mediciones despus del medioda solar, lo que puede responder a condiciones atmosfricas particulares (cristales de hielo en la alta atmsfera, etc.). Los tres mtodos arrojaron valores levemente inferiores a los medidos pero el ajuste es suficientemente bueno para realizar predicciones, por lo menos en das claros.

800 700 600 Irradiancia (W/m2) 500 400 300 200 100 0 4 6 8 10 12 14 16 18 Hora S olar Valores M edidos Estimados por Page Estimados por Hottel Estimados por Liu-Jordan

Figura 8: comparacin de valores de irradiancia global medidos y estimados con GEOSOL. CONCLUSIONES En este trabajo se present la versin 2.0 para Windows del programa GEOSOL cuya primera versin, que data del ao 1991, corra bajo DOS. Las mejoras introducidas comprenden principalmente a la parte grfica y la inclusin del mtodo de estimacin de irradiacin solar de Page. Se subsanaron adems algunos errores existentes en el clculo del ngulo azimutal del sol y se incluyeron dos modelos de cielo anisotrpico para el clculo de la radiacin difusa. Los excelentes resultados obtenidos al validar los valores calculados por el programa con los medidos en la Universidad Nacional de Salta evidencian su potencialidad como herramienta de estimacin cuando no se dispone de valores de radiacin solar medidos en el lugar donde se ha de realizar algn tipo de aprovechamiento solar. La utilizacin de uno u otro mtodo depender de los datos de que se disponga. As, si slo se conoce la ubicacin y altura del lugar deber emplearse el mtodo de Hottel. Si adems se dispone de valores de temperatura y humedad relativa media se podr utilizar el mtodo de Page y si slo se conoce el valor de la irradiacin diaria media mensual sobre superficie horizontal deber emplearse el mtodo de LiuJordan. Con cualquiera de los tres mtodos se pueden realizar estudios para determinar la pendiente ptima de una superficie colectora de acuerdo al tipo de aprovechamiento (estacional o anual) que se desee realizar con ella. La carta solar y el graficador de obstculos son dos valiosas herramientas a la hora de instalar colectores solares o paneles fotovoltaicos ya que proporcionan toda la informacin necesaria sobre el posicionamiento del sol a lo largo del da y su posible apantallamiento, en cualquier poca del ao, con los obstculos circundantes al predio elegido para la instalacin. El grfico 3D de la trayectoria del sol es muy interesante desde el punto de vista educativo ya permite comparar cualitativamente las diferencias existentes entre el invierno y el verano y entre distintos puntos del planeta. Estos tres grficos y todos los datos numricos calculados con GEOSOL pueden imprimirse para su almacenado en papel. REFERENCIAS Duffie, J. A. & Beckman, W. A., 1991, Solar Engineering of Thermal Processes, Whiley Interscience Publications, USA. Hottel, H. C., (1976), A Simple Model for Estimating the Transmittance of Direct Solar Radiation Through Clear Atmospheres, Solar Energy 18, 129. Liu, B. y Jordan, R., (1960), The Interrelationship and Characteristic Distribution of Direct, Diffuse and Total Solar Radiation, Solar Energy, 4(3), 1. Reindl, D., Beckman, W. and Duffie, J., (1990), Evaluation of Hourly Tilted Surface Radiation Models, Solar Energy, 45, 9. ABSTRACT The knowledge of the solar position along the day is of supreme importance for all the applications that take advantage of the solar energy. On the other hand, if we doesn't have measured values of solar irradiation it is necessary to have analytic methods that allow to estimate it. In this work the code GEOSOL V.2.0 for Windows who allows, through its calculations, to fulfill these two objectives is presented. Its graphic output is useful not only from a practical point of view but also educational since it represents the trajectory of the sun in 2D (solar chart) and 3D, facilitating to compare the differences between winter and summer and between different places of the planet. It possesses three methods to estimate solar irradiation: Page and Hottel methods for clear days and Liu-Jordan method for means monthly days. The validation of three methods with mensurations carried out in the National University of Salta is presented, denoting an excellent agreement between measured and estimated values.