RADIACIÓN Y ENERGÍA SOLAR

LIC. FERNANDO NOLLAS

1
TRABAJO PRÁCTICO DE RADIACIÓN
SOLAR PARA EL 05-05-2022

ENVIAR POR MAIL A

[email protected] 2
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

3
ESPECTRO ELECTROMAGNÉTICO

4
FUENTE DE RADIACIÓN ELECTROMAGNÉTICA

https://misistemasolar.com/el-sol-y-su-temperatura/
5
RADIACIÓN DE CUERPO NEGRO

LEY DE PLANCK

LEY DE WIEN

LEY DE STEFAN-BOLTZMANN

6
7
LA TIERRA COMO CUERPO NEGRO

8
 ROTACIÓN DIFERENCIAL DEL SOL

http://fundacionalfaomegamc.blogspot.com/2016/01/que-es-la-rotacion-diferencial-del-sol.html

9
Ciclo solar

10
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

11
CAPAS DE LA ATMÓSFERA Y EL TOPE
EL SOL DESDE LA TIERRA

13
14
 El movimiento de revolución de Su apartamiento respecto a una
la Tierra alrededor del Sol se circunferencia es del 1.7%. La
realiza en una órbita elíptica de Tierra esta mas cerca del Sol en
baja excentricidad enero y mas lejos en junio.

Este cambio en la distancia

produce una variación de la
I=Ics (1+0.033 cos[360°(N-2)/365] )
irradiancia solar ( I ) que alcanza
el tope de la atmósfera terrestre Ics =constante solar
y puede ser aproximada por I=irradiancia

15
POSICIÓN SOBRE LA SUPERFICIE TERRESTRE

16
 LATITUDES IMPORTANTES

https://www.aristasur.com/contenido/sistema-de-coordenadas-geograficas-longitud-y-latitud

17
ÁNGULO DE DECLINACIÓN

18
 DECLINACIÓN

https://personales.unican.es/gonzalmi/ssolar/articulos/movapsol.html

19
 DECLINACIÓN
δ=23.45°

DECLINACIÓN

δ= -23.45°

20
ESTA VARIACIÓN PERIÓDICA
EN LA DECLINACIÓN PUEDE δ=23.45°sen(360°(N–81)/365)
SER APROXIMADA POR
N=DIA JULIANO

21
¿POR QUÉ LOS DIAS TIENEN 24 HS?

¿EL SOL SIEMPRE ESTA EN LA MISMA POSICION A LA MISMA HORA?

22
 La rotación diaria de la Tierra sobre su eje hace que el Sol parezca dar una vuelta a su
alrededor en un día

Se define

El día solar VERDADERO es el intervalo de tiempo (generalmente no 24 hs)

transcurrido entre dos pasajes sucesivos del Sol sobre un observador estacionario en
la Tierra
-Los días solares no son uniformes en longitud debido a la elipticidad de
la órbita terrestre alrededor del sol. Esto genera que la velocidad de
traslación de la Tierra crezca cerca del sol(perihelio) y disminuya lejos del
sol(afelio).
-Se define entonces el día solar medio el cual es el promedio anual del día
solar y se divide en 24 hs iguales.
-Se define el tiempo solar como la medida de tiempo basada en el día
solar (o verdadero). Puede ser medido con un reloj de sol y toma un valor
máximo de ts=12hs en lo que se denomina “mediodía solar” ( el punto de
máxima altura solar del día ) 23
 ECUACIÓN DEL TIEMPO

Se define “Ecuación del Tiempo” a la diferencia entre el día solar aparente y el día
solar medio.
Esta diferencia es función del día del año y debe ser tenida en cuenta para corregir la
hora reloj y conocer la hora solar.

E [min]=9.87 sen(2B) - 7.53 cos(B)-1.5 sen(B)

Esta función se suele
aproximar mediante
con B=360°(N-81)/365
N=dia juliano 24
 HUSOS HORARIOS

4 3

https://www.diferenciador.com/husos-horarios/
25
 HORA SOLAR
Como en una rotación completa,
Cada hora el Sol parece rotar
ejecutada en aproximadamente 24
360°/24=15° ( o 1° cada 4 minutos)
hs, la Tierra gira un ángulo de 360°
Para determinar la hora solar aparente, ts, deben considerarse tres correcciones
muy importantes.

1) Ecuación del tiempo

2) Corrección por diferencia en Longitud del observador (Lobs) y longitud del
uso horario oficial(Lof) Hora del uso donde se encuentra la capital
3) Corrección por cambio de horario legislado (Δ) Normalmente para
ahorro de energía

La hora solar se calcula entonces con la siguiente ecuación:

ts(hs)= tof + Δ + 4min/° (Lof – Lobs) + E

26
 HORA SOLAR

Una simplificación es considerar que actualmente estamos en hora UTC-3 y que, para ese
huso horario, la longitud oficial es 45°.

ts(hs)= tof + 4min/° (45° – Lobs) + E

Para hoy (N=88) E= - 7,21 min y hora local 18:30 hs

ts= 18,5 + (4 * ( 45 – 68,83 ))/60 - 7,21/60 = 16,79 = 16:47 hs

27
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

28
 RADIACIÓN SOBRE UNA SUPERFICIE

Aspectos importantes a la hora de hablar sobre la incidencia de radiación

ÁNGULO DE INCIDENCIA

29
http://www.geo.mtu.edu/KeweenawGeoheritage/Lake/Temperature.html

30
 Θ

LEY DE LAMBERT

31
 ÁNGULO SOLAR HORARIO (w)

El ángulo horario es una variable angular asociada con el tiempo.

-Es positivo durante la mañana

-Toma el valor nulo al mediodía solar Los signos cambian para el HN
-Es negativo durante la tarde

w=15 °/h (12- ts)

EJEMPLO: Calcular el ángulo horario para la ciudad de Buenos Aires para el día 3 de febrero a las
13 hs ( hora oficial) conociendo que LongBs As= 58.38° O

Se obtiene que ts= 11.88 para este caso w= 1.87°

Para la ciudad de Mendoza el ángulo horario sería w=12.3°

Eso quiere decir que en Mendoza falta, “angularmente” más que en BS AS

para llegar al mediodía solar
32
ÁNGULOS DE LA SUPERFICIE (NO SOLAR)

β= inclinación de la superf.
γ= ángulo acimutal de la superf, orientación.

33
 RELACIONES ENTRE LOS ÁNGULOS SOLARES

Las diferentes magnitudes angulares nombradas hasta ahora se relacionan bajo la

siguiente expresión:

cos(θ)= sen(δ) sen(φ) cos(β) – sen(δ) cos(φ) sen(β) cos(γ) +

+ cos(δ) cos(φ) cos(β) cos(w) + cos(δ) sen(φ) sen(β) cos(γ) cos(w) +
+ cos(δ) sen(β) sen(γ) sen(w)

φ = latitud
δ= declinación
θ= ángulo de incidencia (entre normal de la superf. y
dirección de incidencia)
w= ángulo horario
β= inclinación de la superf.
γ= ángulo acimutal de la superf.

34
 CASOS PARTICULARES

cos(θ)= sen(δ) sen(φ) cos(β) – sen(δ) cos(φ) sen(β) cos(γ) + cos(δ) cos(φ) cos(β) cos(w) +
cos(δ) sen(φ) sen(β) cos(γ) cos(w) + cos(δ) sen(β) sen(γ) sen(w)

1) Para planos en el HS orientados hacia el norte e inclinados un ángulo β se tiene que

γ=0° y la expresión queda

cos(θ)= cos(φ - β) cos(δ) cos(w) + sen(φ - β) sen(δ)

2) Para el caso de un plano horizontal (β=0) se obtiene para el ángulo cenital la siguiente
expresión

cos(θz)= sen(δ) sen(φ) + cos(δ) cos(φ) cos(w)

35
 ÁNGULO HORARO DE SALIDA / PUESTA DEL SOL

Cuando se desea calcular la duración del día ello se logra determinando el

ángulo horario de salida del sol, lo que corresponde a un ángulo de
incidencia de 90°

Para el caso de un plano cos(90)=0= sen(δ) sen(φ) + cos(δ) cos(φ) cos(ws)

horizontal esa condición
lleva a obtener la siguiente
expresión para el ángulo
cos(ws)= -tg(δ)tg(φ)
horario de salida del sol:
ws= cos-1( -tg(δ)tg(φ))

Este ángulo se puede utilizar en situaciones en las que se desea conocer la hora de salida o
puesta del sol para un sitio especifico en un momento del año específico.
También lo vamos a usar a la hora de estudiar la radiación solar que hay en un día
completo, o sea desde el amanecer hasta el anochecer.
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACIÓN SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

37
 DIFERENCIA ENTRE POTENCIA Y ENERGIA

ENERGÍA ES LA POTENCIA INTEGRADA EN UN DETERMINADO PERIODO DE TIEMPO

IRRADIANCIA SOLAR POTENCIA INSTANTANEA (W/m2)

IRRADIACIÓN SOLAR ENERGÍA EN UN DETERMINADO

PERÍODO DE TIEMPO (J/m2)

MINUTO,HORA,DIA,ETC
RADIACIÓN SOLAR EXTRATERRESTRE

Como se comentó antes, la energía que llega del Sol al tope

de la atmósfera varía en pequeña proporción debido a la
excentricidad de la órbita terrestre.

Si se quiere calcular la IRRADIANCIA solar que incide en el

tope de la atmósfera normalmente sobre una superficie
horizontal paralela a la superficie sobre un punto de la
Tierra, se usa la siguiente expresión:

I0 =Ics (1+0.033 cos[360°(N-2)/365] ) cos(θz) cos(θz)= sen(δ) sen(φ) + cos(δ) cos(φ) cos(w)

Ics =1366 W/m2

39
 IRRADIACIÓN DIARIA

Si se quisiera conocer la irradiación solar recibida al tope de la atmósfera durante

todo un día habría que integrar entre los ángulos de salida y entrada del Sol y daría la
siguiente ecuación:

H0(J/m2. día) = (24/π) 3600 Ics [1+ 0.033cos (360° (N-2)/365] *

* [cos(φ)cos(δ)sen(ws) + ws (π/180°) sen(φ)sen(δ)

Radiación extraterrestre horaria

Resulta también de interés calcular la irradiación extraterrestre sobre superficie

horizontal para un período de una hora. Integrando la ecuación de la irradiacia
extraterrestre entre los ángulos horarios w1 y w2 que limitan la hora de interés se
obtiene:

12  3600  r   (w 2 − w1 )
2

Io = GSC  0   sensen + cos  cos ( senw 2 − senw1 )
 r  180 

(Io en J/m2·hora)

Otra forma de calcular la irradiación horaria es calcular la irradiancia

correspondiente al punto medio de la hora y multiplicar este valor por 3600 s.

La diferencia entre las irradiaciones horarias calculadas por estos dos métodos será
un poco más grande en las horas próximas a la salida y puesta del sol pero aún así
será pequeña.
41
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

42
 EQUILIBRIO RADIATIVO EN LA ATMÓSFERA

Nube de pequeño espesor. El color blanco Nube de varios kilómetros de espesor. Al

bajo la nube indica la fuerte presencia de pie de la nube no llega radiación. Los
todas las longitudes de onda del espectro distintos tonos de gris indican el espesor de
visible la nube 43
Equilibrio radiativo de radiación solar, atmósfera y
superficie terrestre

Fuente: https://www.weather.gov/jetstream/energy
Atenuación de la radiación solar por la atmósfera

Al llegar la radiación a la atmósfera sufre dos fenómenos muy importantes al entrar en

contacto con ella. Ellos son:
1) Dispersión o scattering: es causada por la interacción de la radiación con las moléculas de
aire, agua (vapor y gotas) y de polvo.
2) Absorción atmosférica: es causada por la absorción de radiación solar de moléculas de
ozono, agua y dióxido de carbono. Ocurre para longitudes de onda del U.V., visible e I.R.

45
 ABSORCIÓN SELECTIVA

Los principales absorbedores

atmosféricos son el vapor de agua,
el dióxido de carbono y las
moléculas de ozono; son moléculas
triatómicas y producen muchas
bandas de absorción desde el
ultravioleta hasta el infrarrojo
lejano. El CO2, vapor de agua, N2O,
metano y otros son considerados
gases de efecto invernadero.

Para resumir, se puede decir que el

fenómeno de absorción selectiva es
particularmente marcado en el
ultravioleta, debido a la presencia
del ozono, y en el infrarrojo, por las
fuertes bandas de absorción del
vapor de agua y de otros gases. 46
 PROCESOS DE DISPERSIÓN

DISPERSIÓN (“SCATTERING”)

Este fenómeno es el más importante en la parte visible del espectro

solar. Deben distinguirse dos casos según el tamaño de las partículas
dispersoras:
a) si la longitud de onda de la luz es mucho mayor que el tamaño de
las partículas (dispersión de Rayleigh), que es lo que ocurre con las moléculas
de los gases atmosféricos, se observa que la distribución espacial de la
radiación dispersada es simétrica y mayor en la dirección de avance y de
retroceso que en la transversal, y que la cantidad de radiación dispersada es
inversamente proporcional a la cuarta potencia de la longitud de onda (esto
permite explicar el color azul del cielo limpio al mediodía).

47
DISPERSION DE RAYLEIGH

48
 Dispersión de MIE

b) si la longitud de onda de la luz es del orden o menor que el tamaño

de las partículas (dispersión de Mie), que es el caso de la interacción
con los aerosoles, la distribución es más complicada y no es simétrica
(presenta simetría axial en la dirección de avance de la radiación):
mayor cantidad de energía es dispersada en la dirección de avance que
en la de retroceso. En la atmósfera, el polvo y el agua tienden a
aparecer en grandes moléculas debido a la agregación de moléculas de
agua y a la condensación de agua en partículas de polvo.

49
 DISPERSION DE MIE

La dependencia de la dispersión respecto a la longitud de onda es pequeña, y

disminuye conforme aumenta el tamaño de las partículas (lo que explica el color
blanco del cielo en presencia de nubes).

50
 RADIACIÓN SUPERFICIAL
CON Y SIN NUBES
CIELO CON NUBOSIDAD

CIELO DESPEJADO

CIELO TOTALMENTE NUBLADO 51

CON NUBES DENSAS
52
 Impacto de la cobertura nubosa en la
radiación

Attenuation
UV Index

Clear Sky Model

reference

E. Wolfram , J. Salvador, F. Orte, D. Bulnes, R. D’Elia, M. Antón, L. Alados-Arboledas And E. Quel-Study Of Cloud Enhanced Surface UV Radiation At
Atmospheric Observatory Of Southern Patagonia, Río Gallegos, Argentina.. AIP CONFERENCE PROCEEDINGS. Melville, New York: AIP CONFERENCE
PROCEEDINGS. 2013 vol.1531 n°. p907 - 910. 53
Impacto de la cobertura nubosa en la
radiación

Enhancement
UV Index

Clear Sky Model

reference

UTC (hs)
 Impacto de la cobertura nubosa en la
radiación

Measurements ; minutal time rate

Model : represent clear sky condition
-Orte F, Wolfram E., Salvador J., D’Elia R., Quiroga J., Quel E., Mizuno A. Attenuation By Clouds Cover Abnormal UV Radiation Impact for Low
Stratospheric Ozone Situations. International Radiation Conference AIP Conference Proceedings 1810, 110009 (2017); doi:
http://dx.doi.org/10.1063/1.4975571
 Influencia de la cobertura nubosa en la Radiación Solar
Día Despejado
Día Despejado 25/11/2017
25/11/2017

Modelo Iqbal
Inputs: Año, mes, día ,elevación solar, HR, temperatura, presión, alfa, beta, single scattering
albedo, factor de corrección, ozono, albedo.

Modelo Paramétrico Iqbal Model C, 1983

 Influencia de la cobertura nubosa en la Radiación Solar

12:20

Δt1=1min → ΔI1=863,0W/m2
incremento ~ 18%

12:19
12:19 Δt
 Influencia de la cobertura nubosa en la Radiación
Resultados Solar
Diciembre 2017

Días de incremento
58
 Latitud Cobertura nubosa
Altura sobre el nivel del mar

¿DE QUÉ DEPENDE LA RADIACION QUE ALCANZA

UNA SUPERFICIE EN UN DETERMINADO SITIO?

Albedo superficial
Columna total de agua
Aerosoles Obstáculos

59
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

60
 Medición de la radiación global

Piranómetros termoeléctricos

Los piranómetros termoeléctricos

tienen como elemento sensible una
termopila de Möll, denominándose
así a un conjunto de pares
termoeléctricos o termocuplas en
serie (vale la pena recordar que una
termocupla es un par de metales
determinados de manera tal que, al
hallarse sus junturas a diferentes
temperaturas, se genera una
diferencia de potencial entre ellas
que es proporcional al salto
térmico). Así, la salida eléctrica es
proporcional a la irradiancia.

61
 CARACTERISTICAS DE UN PIRANOMETRO

RANGO ESPECTRAL
SENSIBILIDAD
TIEMPO DE RESPUESTA
NO-LINEALIDAD
NO-ESTABILIDAD
DEPENDENCIA DE LA TEMPERATURA
ERROR DE INCLINACION
RESPUESTA DIRECCIONAL (LEY DE COSENO)
INCERTIDUMBRE

NORMA ISO 9060 :Specification and classification of instruments for

measuring hemispherical solar and direct solar radiation 62
 Medición de la radiación global

Solarímetros fotovoltaicos

Los solarímetros fotovoltaicos poseen

como elemento sensible una oblea de
silicio que al recibir un flujo radiante
genera una corriente eléctrica
proporcional la que, al circular por una
resistencia, da origen a una diferencia
de potencial; esta respuesta no es
plana sino selectiva, dependiendo de la
longitud de onda de la radiación, lo
cual limita su uso a la determinación
de valores diarios de la irradiación
global una vez que se ha comparado
estadísticamente su comportamiento
con el de un piranómetro
termoeléctrico y determinado su
constante media.

63
Medición de la irradiación global

64
 Medición de la radiación directa normal
La irradiancia solar directa es la procedente de manera directa del sol y
de la zona circunsolar (3° aprox alrededor del sol) sin modificar su
dirección al atravesar la atmósfera.

PIRHELIOMETRO

Fuente: https://www.powerfromthesun.net/Book/chapter02/chapter02.html

65
 Usos de la radiación directa

http://www.sc.ehu.es/sbweb/energias-renovables/temas/termoelectrica/revision/revision.html

66
https://www.reutersevents.com/renewables/csp-es/el-crecimiento-de-la-energia-solar-en-india-liderado-por-la-fv-pone-la-mirada-en-los
 Medición de la radiación difusa
La irradiancia solar difusa es la que resulta de la dispersión en
el aire y las nubes y se mide obstruyendo la radiación solar
directa en un sensor. De ésta manera el sensor recibe
únicamente la radiación solar proveniente de la atmósfera.

67
En días despejados representa un 10-15% de la radiación total incidente pero
en días con nubosidad la dispersión aumenta disminuyendo la radiación
global y aumentando la difusa.

QUE OCURRE UN DIA NUBLADO?

68
 COMPONENTES DE LA RADIACIÓN

RADIACIÓN SOLAR GLOBAL RADIACIÓN SOLAR DIFUSA RADIACIÓN SOLAR DIRECTA

EN IRRADIANCIA GLOBAL= DIFUSA + DIRECTA x cos (SZA)

ENERGÍA GLOBAL= DIFUSA + DIRECTA

69
 ESTIMACIONES INDIRECTAS
Estimación de la irradiación global a través de
heliógrafos

Heliógrafos
Los heliógrafos tienen por objeto medir
la duración de la insolación solar
(heliofanía efectiva). El principio de
funcionamiento de este equipo se basa
esencialmente en la convergencia de
los rayos solares sobre una faja de
papel que provoca que ésta se queme,
lo que ocurre cuando los valores de la
radiación solar directa superan un
cierto umbral que presenta cierta
variabilidad, dependiendo de la
ubicación geográfica del instrumento,
del clima y del tipo de faja utilizada.

70
 RELACIONES ESTABLECIDAS ENTRE LA IRRADIACIÓN GLOBAL Y OTRAS
VARIABLES METEOROLÓGICAS DISPONIBLES

La amplitud térmica diaria ΔT ( = Tmax – Tmin ) es otro de los parámetros que han
utilizado algunos autores para relacionar la irradiación solar de un lugar con parámetros
meteorológicos locales

El Atlas de Energía Solar del Perú,

publicado en 2003, fue realizado
utilizando mediciones de
estaciones piranométricas y
completado con información de
amplitud térmica según la relación:

71
Raichijk en 2005 utilizó una expresión lineal de la raíz cuadrada del promedio de
amplitud térmica diaria propuesta por Hargreaves ( 1985).

Utilizaron estaciones distribuidas a lo largo del país y

luego realizaron cartas donde mostraron la
distribución espacial del promedio mensual de
irradiación diaria. Obtuvieron un error cuadrático
medio relativo (RMSE) entre 9% y 14,2% y un desvío
medio relativo (MBE) entre -1,2% y -7,1%.

VOLVER
72
ESTIMACIÓN DE IRRADIACIÓN SOLAR A PARTIR DE INFORMACIÓN SATELITAL

73
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

74
 RADIACIÓN SOBRE SUPERFICIES INCLINADAS
Como ya hemos visto, variar el ángulo de la superficie puede proporcionarnos mas o menos
radiación sobre la misma porque varía el angulo de incidencia

LA MAYOR PARTE DE LA INFORMACIÓN DISPONIBLE SE ENCUENTRA EN LA HORIZONTAL

así que partir de ahora vamos a estudiar métodos para estimar la irradiación sobre un
plano inclinado partiendo de la base de que se conoce la radiación sobre plano
horizontal. 75
ACLARACIÓN: VAMOS A OBTENER EXPRESIONES PARA RADIACION INCIDENTE EN
SUPERFICIES ORIENTADAS HACIA EL ECUADOR (=0 Ó =π SEGÚN LA BIBLIOGRAFÍA)

La irradiación e irradiancia que inciden sobre un plano inclinado se

pueden separar en sus componentes:

Donde:

Gbi= es la radiación directa que incide en el plano

Gdi= es la radiación difusa que incide en el plano
Gri= es la radiación proveniente de superficies reflectoras

76
 Componente DIRECTA:

Se define un factor Rb como el cociente entre la proyección de la irradiancia incidente

sobre la normal a un plano inclinado un ángulo β y la proyectada sobre la normal a un
plano horizontal

Donde s=1 para el HS (ya que φ<0) y s=-1 para el HN

77
 Componente REFLECTIVA:

La irradiancia proveniente de reflexión en superficies cercanas se puede expresar como

Donde para la superficie j-esima, Ij es la radiación en incidencia normal, Aj el área y Fc-j

es el factor de vista con que la superficie j ve al colector. j representa la reflectancia
global en todas las longitudes de onda

PARA SIMPLIFICAR SUPONDREMOS QUE SOLO SE REFLEJA EN EL SUELO EL CUAL SE

CONSIDERA HORIZONTAL ENTONCES QUEDA:

78
Componente DIFUSA:

Existen dos enfoques para esta componente:

MODELO ISOTRÓPICO MODELO ANISOTRÓPICO

Directa

Cielo

Difusa

β Suelo

79
MODELO ISOTRÓPICO

Directa

Cielo

Difusa

β Suelo

 1 + cos  
I di = I dh  
 2 

80
Modelo difuso anisotrópico (más realista que el anterior)
Este modelo divide a la irradiación difusa en tres componentes: difusa isótropa de
fondo, difusa circunsolar y difusa de horizonte. Es conocido como el modelo Hay,
Davies, Klucher, Reindl (autores del mismo). La radiación difusa depende en este
modelo de dos factores, Ai y f, que definen las fracciones circunsolar y de fondo:

  1 + cos    3    
I d , p = I d  Ai Rb + (1 − Ai )  1 + f sen   
  2   2  
I I Indice de anisotropía
Ai = bn = b
I on I o
f = Ib / I Factor de modulación

81
 LA RADIACIÓN TOTAL QUE RECIBE UNA PLANO INCLINADO SERÁ ENTONCES

ISOTRÓPICO

 1 + cos    1 − cos  
I p = I b Rb + I d   + I  
 2   2 

ANISOTRÓPICO

 1 + cos   3     1 − cos  
I p = (I b + I d Ai ) Rb + I d (1 − Ai ) 1 + f sen   + I  
 2   2   2 

82
 INCLINACIÓN ÓPTIMA
Para sacar el mayor provecho en invierno conviene inclinar mas la superficie para que la
incidencia se mas normal a la misma.

https://www.tutiendaenergetica.es/blog/63_rendimiento-instalacion-solar-invierno

83
 INCLINACIÓN ÓPTIMA

https://tipos-de-energia.blogspot.com/2006/02/energia-solar-trmica.html

CONCLUSIÓN: SE ESTUDIA CON DATOS 84

 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

85
 MODELOS DE CIELO CLARO PARA ESTIMAR RADIACION SOLAR

https://www.sciencedirect.com/science/article/pii/S1364032119302187

86
Estimación de irradiación solar de día claro. Método de Page

Este método, aplicable solamente a los días sin nubosidad, fue desarrollado por el
grupo del Prof. J. K. Page de la Universidad de Sheffield en la década del ´70. Para
su uso es necesario dar el contenido de humedad del aire (como presión parcial de
vapor) y un parámetro ajustable que depende del lugar, la turbidez atmosférica ta,
la cual toma en cuenta la atenuación de la radiación solar debida al polvo y demás
aerosoles presentes en la atmósfera.
Primero se calcula la irradiancia directa teórica sobre superficie normal al haz
bajo una atmósfera sin aerosoles, después de interactuar con el aire, ozono,
vapor de agua y anhídrico carbónico. La irradiancia solar directa a nivel del suelo
resulta entonces:

Gbn = Gon .tairema . taguamw . tozonomo . tCO2mg . Exp(-ta . ma) (ma: masa de aire)

La irradiancia directa de cielo claro de Page sobre plano horizontal es:

Gbc = Gbn.cos z

La irradiancia directa de cielo claro de Page sobre plano inclinado es:

Gbp = Gbn .cos  = Gbc.Rb
87
En cuanto al modelo de cielo considera una distribución anisotrópica de la radiación
difusa dividiéndola en dos componentes: isótropa de fondo, Gdf, y difusa circunsolar
Gdc la que, desde el punto de vista geométrico, se trata como directa. La irradiancia
difusa sobre plano horizontal es entonces:

Gd = Gdf + Gdc = Gdf + Gdcn .cos z , con Gdcn: difusa circunsolar normal al haz

Gdf = c.b0 – b1.Gbc y Gdcn = c.c0 – c1.Gbc

Donde: 7 7

 
Gon
c= , bj = b ji (cos z ) i
y cj = c ji (cos z )i
1367 i =0 i =0

i b0i b1i c0i c1i

0 2,000 272,000 1,000 452,400
1 331,965 371,162 536,917 1.549,010
2 -658,223 -9.332,020 -802,612 -12.294,800
3 4.356,270 45.822,100 3.836,530 36.647,200
4 -15.563,000 -108.407,000 -12.557,000 -60.846,000

5 26.253,800 137.668,000 20.257,800 59.732,800

6 -20.505,800 -89.661,600 -15.276,300 -32.313,200
7 6.037,430 23.460,600 4.358,420 7.381,700
88
La irradiancia solar total de Page sobre plano horizontal está dada por:
Gth = Gbc + Gdcn.cos z + Gdf

Directa

Cielo

Difusa Difusa
Circumsolar
Isótropa de
fondo

Difusa de Horizonte
b Suelo
Difusa reflejada en el suelo

La irradiancia solar total sobre plano inclinado está dada por:

Gtp = Gbc.Rb + Gdcn.cos  + Gdf .(1 + cos b) / 2 + r.Gth. (1 – cos b) / 2
89
Estimación de irradiación solar de día claro. Método de Hottel

Hottel desarrolló en 1976 un método para estimar la irradiancia directa transmitida a

través de atmósferas claras (sin nubes) el cual tiene en cuenta el ángulo cenital z del
sol y la altura del lugar (hasta 2,5 km) para cuatro tipos de clima.
La transmitancia atmosférica para radiación directa normal al haz, b, es Gbn/Gon y
está dada por:
b = a0 + a1.exp(-k/cos z)

Las constantes a0, a1 y k para atmósferas con visibilidad de 23 km se encuentra a

partir de a0*, a1* y k*, los cuales están dados para altitudes menores de 2,5 km por:
a0* = 0,4237 – 0,00821.(6 - H)2 , H: altura del lugar en km sobre el nivel del mar
a1* = 0,5055 + 0,00595.(6,5 - H)2
k* = 0,2711 + 0,01858.(2,5 - H)2

Los cambios climáticos son tenidos en cuenta mediante 3 factores de corrección :

r0 = a0/a0* r1 = a1/a1* y rk = k/k*
90
Tipo de clima r0 r1 rk
Tropical 0,95 0,98 1,02
Verano de 0,97 0,99 1,02
latitud media
Verano 0,99 0,99 1,01
subártico
Invierno de 1,03 1,01 1,00
latitud media
La irradiancia directa de cielo claro normal al haz es entonces:
Gbnc = Gon.p

La irradiancia directa de cielo claro sobre plano horizontal es:

Gbc = Gon.p.cos z

La irradiancia directa de cielo claro sobre plano inclinado es:

Gbcp = Gbc.Rb

91
Estimación de irradiación solar horaria a partir de datos diarios.
Método de Liu -Jordan

IRRADIACION DIARIA IRRADIACION HORARIA

(O MEDIA DIARIA MENSUAL)

DIFUSA GLOBAL

Cuando se hacen estudios de performance de artefactos solares, hora por hora,

puede ocurrir que sólo se disponga de datos de irradiación solar media diaria, H, o
diaria media mensual, Ĥ, sobre superficie horizontal.

El método de Liu Jordan permite, mediante curvas obtenidas experimentalmente,

obtener la irradiación solar horaria, I, para cada hora de sol y para planos de
cualquier inclinación a partir de valores de irradiación solar diaria (H o Ĥ) sobre
plano horizontal.

Primeramente se debe determinar la radiación difusa diaria Hd sobre plano

horizontal mediante el uso de correlaciones empíricas que dependen de un
parámetro llamado Indice de Claridad, KT definido mediante: KT= H/ Hext
92
 ÍNDICE DE CLARIDAD

kt para valores horarios Kt para valores diarios Kt para valores mensuales

La correlación de Collares Pereira y Rabl para datos diarios es:

0,99 para KT  0,17


H d  1,188 − 2,272 KT + 9,473 KT − 21,865 KT + 14,648 KT para 0,17  KT  0,75
2 3 4

=
H  − 0,54 KT + 0,632 para 0,75  KT  0,80
 0 ,2 para K T  0,80

93
 RELACIÓN ENTRE LOS VALORES PROMEDIO DE LAS COMPONENTES
HORARIAS SOBRE UN PLANO HORIZONTAL Y SUS PROMEDIOS DIARIOS

Los estudios realizados por Liu y Jordan en 1960 junto con los realizados por Collares-Pereira
y Rabl (1979) dieron como resultado las siguientes expresiones. Si H es el valor medio de la
irradiación global diaria, la correspondiente fracción horaria puede calcularse como:

rt =
Hh 
= (a + b cos ) cos −cos s Donde ws es el promedio
H 24 sen s − s
cos s mensual del ángulo horario de
180 salida del sol y w es el ángulo
horario del centro del intervalo
horario medido respecto del
a = 0,409 + 0,5016sen( s − 60) mediodía solar
b = 0,6609 − 0,4767 sen( s − 60)

Componente difusa

H hd  cos − cos s
rd = =
Hd 24 sen −   s cos
s s
180
94
IRRADIACIÓN MENSUAL

El modelo isotrópico de Liu y Jordan (1961) establece la siguiente relación, donde Hβ es el

valor medio de la irradiación total recibida sobre el plano inclinando

H Hd Hd  1 + cos    1 − cos  
= (1 − ) Rb +   +  
H H H  2   2 

R = [cos( + s ) cos( ) sen( ws ' ) + ( / 180) ws ' sen( + s ) sen( )]

/[cos( ) cos( ) sen( ws ) + ( / 180) ws sen( ) sen( )]

Donde la barra indica valores medios mensuales y ws’ es el ángulo horario de salida
del sol que ve el plano, el cual es el menor entre

arccos(−tg ( −  )tg ( )
ws’=
arccos(−tg ( )tg ( ))

Esto lo usamos si tuviéramos que calcular la irradiación diaria media mensual sobre un
colector inclinado 95
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

96
 FRECUENCIAS TEMPORALES Y PROMEDIOS

MINUTAL HORARIO DIARIO MENSUAL

COMPONENTE HORARIO GLOBAL COMPONENTE HORARIO DIFUSA ECUACIONES

COMPONENTE DIARIA GLOBAL COMPONENTES HORARIOS DIFUSA Y GLOBAL

COMPONENTE DIARIO MEDIO MENSUAL GLOBAL COMPONENTE DIARIA MEDIA MENSUAL DIFUSA
ECUACIONES

97
 ÍNDICE
1. EMISIÓN SOLAR
2. COORDENADAS SOLARES---------> HORA SOLAR
3. ANGULO DE INCIDENCIA
4. RADIACION SOLAR EXTRATERRESTRE
5. LA ATMÓSFERA Y SU INTERACCION CON LA RADIACION SOLAR
6. SENSORES DE RADIACIÓN SOLAR. MEDICIONES DIRECTAS Y ESTIMACIONES
INDIRECTAS
7. RADIACIÓN SOLAR SOBRE SUPERFICIES INCLINADAS
8. MODELOS
9. FRECUENCIA TEMPORAL
10. BASES DE DATOS

98
ATLAS DE RADIACIÓN SOLAR
GROSSI GALLEGOS

99
 GEOSOL

AZIMUTH 180° PARA QUE ESTE MIRANDO AL NORTE

100
 SWERA

http://en.openei.org/apps/SWERA/
https://openei.org/datasets/dataset/swera

CELDAS DE 40 Km
101
102
¿¿Cómo se ve la información???

PARA MENDOZA

103
 COMPARACIÓN PARA LA CIUDAD DE MENDOZA

MEDIDO
8 SWERA
IRRADIACION MEDIA MENSUAL (kW/m2)

6
PUEDE OCURRIR QUE LA
5
RADIACIÓN SOLAR NO SEA
MÍNIMA EN JUNIO
4

2
IL

O
FE RO

EM E
ZO

IO

E
O

O BRE
BR

BR

BR
AY

PT ST
LI
R
ER

N
AR
E

AB

JU
JU

TU
SE GO
M

M
M
EN

BR

IE
IE

VI
A

IC
O

D
N

104
 MAPAS SWERA PARA NORTE DE ARGENTINA

Fuente:
https://www.researchgate.net/publication/258511376_TRAZADO_DE_MAPAS_MEDIOS_ANUALES_DE_ENERGIA_SOLAR_GLOBAL_DIRECTA_DIF
USA_Y_TILT_USANDO_LA_BASE_DE_DATOS_DE_SWERA_CASO_DE_ESTUDIO_PROVINCIAS_DE_SALTA_Y_JUJUY

105
 OTRAS BASES DE DATOS

http://www.soda-pro.com/ https://meteonorm.com

https://www.ecmwf.int/en/forecasts/datasets/ https://gmao.gsfc.nasa.gov/reanalysis/MERRA-
reanalysis-datasets/era5 2/

106
CONCLUSIÓN: EL USO DE MODELOS O DATOS OBTENIDOS A TRAVÉS DE ESTIMACIONES
CONLLEVA UN ERROR CONSIDERABLE POR LO QUE NO HAY NADA MEJOR QUE CONTAR
CON MEDICIONES CONFIABLES REALIZADAS IN-SITU.

A SU VEZ LO QUE ANDA BIEN PARA UNA UBICACIÓN NO SIEMPRE RESULTA IGUAL PARA
TODAS.

PLANILLA

107
 RELACIONES ESTABLECIDAS ENTRE LA IRRADIACIÓN GLOBAL Y LA
COBERTURA DE NUBES (NUBOSIDAD)

Black en 1956 determinó una relación que involucra la irradiación global y la

cobertura de nubes

H/H0= 0.803 - 0.340 C – 0.458 C2 con C<0.8

VOLVER

108
 COMPONENTES HORARIOS

ERBS (1982):

1.0 − 0.09 kT para kT  0.22


0.9511− 0.1604kT + 4.388 kT
2
Id
= kd = 
− + para 0.22  kT  0.8
3 4
I  16.638 k T 12. 336 k T
0.165 para kT  0.8


Reindl y otros (1990) plantearon una correlación similar:

VOLVER
kd

Ridley, Boland y Lauret (2010) desarrollaron un modelo donde consideran los

parámetros mas determinantes para la radiación solar:

kd=
109
Ejemplo: usando el modelo isotrópico, estimar las componentes de la irradiación horaria
y la total para un superficie inclinada 60° hacia el sur ( HN) a una latitud de 40° N para la
hora solar 9 a 10 am el 20 de febrero. La irradiación global fue de I= 1,04 MJ/m2 y ρ=0,6

Para esa hora la irradiación extraterrestre es 2,34 MJ/m2 y entonces el índice de claridad
horario kt será

kt= 1,04/2,34 = 0,445

Según la correlación de Erbs Id/I=0,766 1.0 − 0.09 kT para kT  0.22


0.9511− 0.1604kT + 4.388 kT
2
Id
= kd = 
− 16.638 kT + 12.336 kT para 0.22  kT  0.8
3 4
Id= 0,766*1,04=0,796 MJ/m2 I 
0.165 para kT  0.8

Ib= (1-0,766)*1,04=0,244 MJ/m2

El ángulo horario para el punto medio de la hora es 37,5° ( o -37,5° dependiendo de

como se defina) pero no afecta porque esta dentro de un cos()

110
 COMPONENTES DIARIOS A HORARIOS

MÉTODO DE LIU – JORDAN VISTO ANTERIORMENTE

COMPONENTE DIARIO A DIARIO

Si se cuenta con los datos de IRRADIACIÓN GLOBAL DIARIA entonces se puede estimar
la irradiación difusa DIARIA mediante la correlación de ERBS (1982):

También se puede utilizar

la correlación de Collares-
Pereyra vista
anteriormente

111
 COMPONENTES MENSUALES

Si se cuenta con los datos IRRADIACIÓN GLOBAL DIARIA MEDIA MENSUAL entonces se
puede estimar la irradiación difusa DIARIA media mensual mediante la correlación de
ERBS (1982):

La fracción de irradiación difusa diaria media mensual según Collares Pereira y

Rabl está dada por la siguiente correlación.

__
 
Fd = __ = 0,775 + 0,00606 ( s − 90 ) − 0,505 + 0,00455 ( s − 90 )cos115 K T − 103 
__
Hd
H  

Donde ws es el ángulo horario de puesta del sol para el día medio del mes bajo
consideración.

VOLVER 112