Radiacin Solar Sara Moreno Tejera

8
3. Radiacin Solar

La radiacin solar es la radiacin electromagntica emitida por el sol con una
distribucin espectral muy similar a la producida por un cuerpo negro a 5777 K.
El sol hace llegar una irradiancia directa normal de 1367 W/m
2
a la superficie de la
atmsfera que recibe el nombre de constante solar. Desde la superficie de la Tierra,
dicho valor vara en funcin del momento del da, los componentes atmosfricos y la
latitud.
Se llama irradiancia, denotndola con la letra I, a la energa incidente por unidad de rea
y unidad de tiempo procedente del sol (W/m
2
). Su integracin en un intervalo de tiempo
determinado es la irradiacin (kJ/m
2
o Wh/m
2
), representada por la letra H.


3.1 Geometra Sol - Tierra

3.1.1 Los movimientos de la tierra

La Tierra, al igual que el resto de los planetas y de acuerdo con la primera ley de
Kepler, gira en torno al Sol describiendo una rbita elptica con dicho astro situado en
uno de sus focos. El eje mayor de esta elipse se denomina lnea de los psides, siendo sus
extremos el perihelio, punto de la rbita terrestre ms prximo al Sol, y el afelio, punto
ms alejado. El plano que contiene la rbita de la Tierra se denomina plano de la
eclptica.

La distancia media de la Tierra al Sol, r
o
, define una unidad de longitud
denominada unidad astronmica, UA, cuyo valor es

UA m 5979 . 1 = k
10
49 1 =
r
8
o
[3.1]

La excentricidad de la elipse, es decir, la razn entre la distancia del centro de la
elipse a un foco y el semieje mayor, es igual aproximadamente a 0.017. Al ser esta
excentricidad tan pequea, la rbita terrestre podra asimilarse a una circunferencia de
radio r
o
. No obstante, a lo largo del ao, la distancia Tierra-Sol, r, vara apreciablemente
entre un valor mnimo de 0.983 UA en el paso por el perihelio (aproximadamente el da
3 de enero), y un valor mximo de 1.017 UA en el paso por el afelio (aproximadamente
el 4 de julio). Estas variaciones originan una oscilacin anual de un 3 % en la cantidad
de radiacin solar recibida por la Tierra

El perodo orbital de la Tierra o intervalo de tiempo que transcurre entre dos
pasos consecutivos por el mismo punto de la rbita, define el denominado ao sidreo,
y tiene una duracin de 365.25636 das.





Radiacin Solar Sara Moreno Tejera

9
3.1.1.1. Movimiento de traslacin

El movimiento orbital de la Tierra no es uniforme. La velocidad lineal media es
de unos 29.8 km/s, siendo mxima en el perihelio y mnima en el afelio. Este
movimiento viene determinado por la segunda ley de Kepler o ley de las reas, segn la
cual el radio vector que une el Sol con un planeta, barre reas iguales en tiempos iguales
esto es, la velocidad areolar es constante. A efectos prcticos, interesa conocer la
distancia Tierra Sol un da cualquiera del ao. Las expresiones ms usadas son las
siguientes:
o bien,

0
1 0, 033cos(2 / 365) E n = + [3.3]
donde:

Siendo r la distancia Tierra Sol el da n del ao a partir del 1 de Enero, r
0
la
distancia media Tierra Sol, igual a una UA. Donde es el ngulo del da en radianes.


3.1.1.2. Movimiento de rotacin

Simultneamente al movimiento de traslacin descrito, la Tierra gira sobre un
eje que pasa por los polos, en un movimiento de rotacin que puede suponerse
perfectamente uniforme. El perodo de rotacin de la Tierra, es decir, el tiempo que
tarda en dar una vuelta completa sobre s misma, es ligeramente inferior a un da,
exactamente 23 h 56 min. 4.099 s. El plano normal al eje de rotacin que pasara por el
centro de gravedad de la Tierra se denomina plano ecuatorial o plano del ecuador.

El eje de rotacin de la Tierra no es perpendicular al plano de la eclptica. Esto
da origen a las estaciones del ao - invierno, primavera, verano y otoo - y a la desigual
duracin de los das y de las noches. Durante su movimiento de traslacin alrededor del
Sol, el eje de rotacin de la Tierra se desplaza paralelamente a si mismo, conservando
prcticamente invariable su direccin en el espacio y formando siempre el mismo
ngulo con la normal al plano de la eclptica. Este ngulo se denomina oblicuidad de
la eclptica y su valor es de 23 27'.

Esta oblicuidad permite definir sobre la rbita terrestre algunos puntos
singulares. Para ello, se ha de considerar un plano que pasando por el centro geomtrico
del Sol sea normal al eje de rotacin de la Tierra. La interseccin de este plano, paralelo
al del ecuador, con el plano de la eclptica recibe el nombre de lnea de los equinoccios
y determina sobre la rbita terrestre dos puntos conocidos como equinoccio de
E
0
= 1.000110 + 0.034221cos + 0.001280sen + 0.000719cos 2 + 0.000077sen 2

[3.2]


2
0
0

=
r
r
E [3.4]
( )
365
1 2
=
n
[3.5]
Radiacin Solar Sara Moreno Tejera

10
primavera y equinoccio de otoo. Anlogamente, la perpendicular trazada en el plano de
la eclptica por el centro del Sol a la lnea de los equinoccios se llama lnea de los
solsticios cuya interseccin con la rbita terrestre da lugar a otros dos puntos
denominados solsticio de invierno y solsticio de verano. El equinoccio de primavera,
tambin denominado punto vernal, se utiliza frecuentemente en astronoma como punto
de referencia, ya que su posicin puede determinarse con gran precisin.

Cuando la Tierra est situada sobre los equinoccios, lo que ocurre
aproximadamente los das 21 de marzo y 22 de septiembre, el Sol se encuentra sobre el
plano del ecuador y, por consiguiente la duracin del da natural de luz es igual a la de
la noche en todos los lugares de la Tierra. Hacia el da 21 de junio, fecha del solsticio de
verano, el Sol ilumina ms directamente el hemisferio norte de la Tierra, donde el Sol
alcanza su mxima inclinacin respecto al plano del ecuador. Seis meses ms tarde,
hacia el 22 de diciembre, fecha del solsticio de invierno, la situacin anterior se ha
invertido, siendo el hemisferio Sur, el ms iluminado. En los solsticios la diferencia
entre el da natural de luz y la noche es mxima.

Esta rotacin es la causa de la visibilidad o no (da, noche) del Sol desde un
determinado lugar de la Tierra y, por tanto, proporciona la definicin del tiempo solar
en relacin con la posicin angular relativa observador Sol.

3.1.2 Posicin de un observador sobre la superficie terrestre

La forma real de la superficie de la Tierra se asemeja, con bastante
aproximacin, a un elipsoide de revolucin de dos ejes, denominado elipsoide de
referencia (de radio ecuatorial 6378.16 km y radio polar 6356.77 km), con centro en el
de gravedad de la Tierra y con el eje menor de la elipse generatriz coincidente con el eje
de rotacin de la Tierra. Dicho eje corta al elipsoide en dos puntos, denominados polo
Norte y polo Sur.

A efectos prcticos se puede tomar la Tierra como una esfera cuyo plano
perpendicular al eje determina una circunferencia llamada ecuador terrestre. Los planos
que pasan por los polos y que, por consiguiente, son normales al plano del ecuador, se
denominan planos meridianos y determinan sobre la esfera terrestre unas
circunferencias conocidas como meridianos terrestres. Finalmente, los planos paralelos
al ecuador cortan a la esfera terrestre formando circunferencias denominadas paralelos.

Un punto cualquiera sobre la superficie terrestre queda determinado por sus
coordenadas geogrficas que son:

La longitud, , definida como el angulo diedro entre el plano meridiano cero y el
plano meridiano que pasa por dicho punto. Se cuenta de 0 a 180, positiamente
hacia el Lste y negatiamente hacia el Oeste.

La latitud, , que se deine como el angulo ormado por la normal a la esera
terrestre por dicho punto con el plano del ecuador. Se cuenta de 0 a 90, a partir
del ecuador, positiamente hacia el Norte y negatiamente hacia el Sur.

Radiacin Solar Sara Moreno Tejera

11

Figura 3.1 Coordenadas geogrficas

En la figura 2, se pueden observar ambos ngulos sobre la superficie de la
Tierra, representada en forma intencionadamente exagerada de elipsoide. En la realidad
es mucho ms esfrica.

La inclinacin del eje de rotacin de la Tierra respecto al plano de la eclptica
determina sobre el elipsoide unos paralelos denominados trpicos y crculos polares. En
el solsticio de verano, durante el medioda, los rayos solares caen verticalmente sobre
los lugares situados en el trpico de Cancer, mientras que en el crculo polar antartico el
Sol est en el horizonte. Otro tanto sucede, en el solsticio de invierno, con el trpico de
Capricornio y el crculo polar rtico.

3.1.3 La esfera celeste. Sistemas de referencia

En astronoma, la posicin de un astro se determina de ordinario mediante
coordenadas esfricas. Sin embargo, la distancia a los astros es generalmente
desconocida en principio, por lo que se considera su posicin proyectada sobre una
esfera de radio unidad (o si se prefiere de radio infinitamente grande) y centro arbitrario
(casi siempre se piensa en el centro de la Tierra), denominada esfera celeste.

Se puede suponer la esfera celeste centrada en el lugar de observacin, sobre la
superficie de la Tierra, o en el centro de sta. Las direcciones trazadas a un astro desde
ambos puntos forman un ngulo, conocido como paralaje, que en el caso de la posicin
solar, puede despreciarse.

Para determinar la posicin de un astro sobre la esfera celeste, dada por una
pareja de ngulos denominados coordenadas celestes, es necesario establecer
previamente algn sistema de referencia. Esto exige adoptar un circulo mximo
fundamental, determinado por un plano o un eje fundamental normal al mismo, que
Radiacin Solar Sara Moreno Tejera

12
pase por el centro de la esfera; un punto arbitrario sobre dicho circulo que se toma como
origen y un sentido positivo para ambas coordenadas.


Figura 3.2 Coordenadas celestes horizontales

El primer sistema de referencia a considerar es el sistema de coordenadas
celestes horizontales, que representa geomtricamente el aspecto del cielo visto por un
observador situado en un lugar determinado de la superficie terrestre. Utiliza como
plano fundamental el horizonte, definido como el plano tangente a la esfera terrestre en
el lugar de observacin. Tambin puede ser entendido como el plano normal a la
vertical de observador (lnea que une a este con el centro de la Tierra) y que, para un
observador en un punto cualquiera de nuestro planeta, viene determinado por la lnea
que marca la plomada. Este plano determina sobre la esfera celeste un crculo mximo
denominado horizonte celeste y la recta normal al mismo (vertical del lugar o del
observador), trazada por el centro de la esfera, corta a esta en dos puntos: el cenit, sobre
el observador, y el nadir, oculto por la Tierra, bajo el observador. Los crculos mximos
que pasan por el cenit y el nadir se denominan crculos verticales y son, por
consiguiente, normales al horizonte.

Las coordenadas celestes horizontales de un astro (vase figura 3) son:

El acimut, , arco de horizonte celeste comprendido entre el punto cardinal Sur y
el punto donde el crculo ertical que pasa por el astro corta al horizonte. Se cuenta
a partir del Sur, de 0 a 180, positiamente hacia el Oeste y negatiamente hacia
el Lste.

La altura solar, , arco de dicho crculo ertical comprendido entre el astro y el
horizonte. Se mide a partir del horizonte de 0 a 90, positiamente hacia el cenit y
negatiamente hacia el nadir.

En lugar de la altura se emplea frecuentemente el ngulo cenital,
z
, que es el
arco de crculo vertical entre el cenit y el astro. Se cuenta a partir del cenit, de 0 a 180,
con lo que la relacin entre la altura solar y el ngulo cenital viene dada por:

Radiacin Solar Sara Moreno Tejera

13

z
= 90 [3.5]

Otro sistema de reerencia que se puede utilizar para la esera celeste, es el denominado
sistema de coordenadas celestes horarias basado en el moimiento de rotacin de la 1ierra.
Como eje undamental se adopta el eje de rotacin, que deine sobre dicha esera los polos
celestes, Norte y Sur, y el ecuador celeste, normal al mismo, como crculo maximo
undamental. Se llaman crculos horarios o meridianos celestes los crculos maximos
secundarios que pasan por los polos y paralelos celestes los crculos menores paralelos al
ecuador. Ll semicrculo horario que partiendo del polo Norte celeste, pasa por el cenit, corta al
horizonte en el punto cardinal Sur y inalmente, pasa por el polo Sur celeste, recibe el nombre
de meridiano superior del lugar.

Las coordenadas celestes horarias de un astro son, como puede verse en la figura
4:



Figura 3.3 Coordenadas celestes horarias

El ngulo horario, , que es, por deinicin, el arco del ecuador celeste ,o el
angulo diedro correspondiente, comprendido entre el meridiano superior del lugar y
el circulo horario que pasa por el astro. Se mide sobre el ecuador a partir del punto
de interseccin entre el meridiano superior del lugar y aquel, de 0 a 180,
positiamente hacia el Oeste y negatiamente hacia el Lste.

La declinacin, , es el arco del crculo horario que pasa por el astro, comprendido
entre el ecuador y ste. Se cuenta a partir del ecuador de 0 a 90, positiamente
hacia el polo Norte y negatiamente hacia el polo Sur.

El ngulo que forman los planos del horizonte y del ecuador depende de la
latitud geogrfica del observador. En efecto, la altura del polo Norte celeste sobre el
horizonte es igual a la latitud geogrfica del observador, afirmacin que sigue siendo
vlida en el hemisferio Sur, donde el polo Norte est por debajo del horizonte y la
Radiacin Solar Sara Moreno Tejera

14
latitud es negativa. As pues, el ngulo formado por horizonte y ecuador es el
complemento de la latitud geogrfica del observador.

Las coordenadas celestes horizontales de un astro vienen dadas, en funcin de
sus coordenadas celestes horarias y de la latitud del observador, por las siguientes
relaciones trigonomtricas, obtenidas de las relaciones de Bessel del tringulo esfrico
formado por el polo (Norte en el hemisferio Norte y Sur en el hemisferio Sur), el cenit y
el astro.

La expresin ms simple usada para calcular el acimut viene dada por la expresin
3.6:

cos cos
sen sen sen
cos
-
= [3.6]


Los dos sistemas de coordenadas enumerados, horizontales y horarios, dependen
del lugar en que est situado el observador sobre la superficie terrestre, por lo que se
denominan por este motivo coordenadas celestes locales. En efecto, las coordenadas
horizontales, acimut y altura, de un mismo astro, en un determinado instante, son
diferentes para dos observadores situados en distintos lugares de la Tierra. Otro tanto
ocurre con el ngulo horario del astro, salvo que ambos observadores estn situados
sobre el mismo meridiano. Por el contrario, todos los observadores de la Tierra,
cualquiera que sea su lugar de observacin, miden la misma declinacin de un
determinado astro.

3.1.4 Movimiento aparente del Sol sobre el horizonte

Las coordenadas horizontales de un astro varan continuamente con el tiempo.
Como consecuencia del movimiento de rotacin de la Tierra, los astros describen sobre
la esfera celeste crculos paralelos al ecuador que recorren diariamente en sentido
retrgrado (de Este a Oeste) y sobre el que pueden destacarse algunos puntos singulares.
Se denomina orto de un astro el momento de su cruce con el horizonte hacindose
visible por el Este, y ocaso el momento en que cruza el horizonte por el Oeste
ocultndose. Asimismo, el instante de paso de un astro por el meridiano superior del
lugar ( = 0), en el que alcanza su altura mxima sobre el horizonte, se conoce como
culminacin superior del astro.

Para un astro carente de movimiento aparente propio (por ejemplo, una estrella
lejana), este movimiento diario se realiza exactamente por un paralelo celeste por lo
que, como el movimiento de rotacin de la Tierra es uniforme, se tiene que el intervalo
de tiempo transcurrido entre dos culminaciones superiores consecutivas de tal astro, que
coincide con el perodo de rotacin de la Tierra es constante y se denomina da sidreo y
constituye una unidad fundamental de tiempo en astronoma. Su duracin es de 23 h 56
min. 4.091 s.

Sin embargo, el Sol, debido al movimiento orbital de la Tierra, est animado de
un movimiento propio aparente con respecto a las estrellas fijas, que realiza en sentido
directo (de Oeste a Este). Este movimiento se realiza sobre la eclptica, que est
inclinada un ngulo ~ 23 2' sobre el ecuador. Por tanto, la declinacin solar ara a lo largo
Radiacin Solar Sara Moreno Tejera

15
del ano, para el hemiserio Norte, entre un maximo de ~ - c, en el solsticio de erano y un
mnimo de o ~ - c, en el solsticio de inierno, y al res para el Sur, siendo o ~ 0 durante los
equinoccios en ambos casos.

La composicin de ste con el movimiento diario causado por la rotacin
terrestre, determina el movimiento aparente del Sol sobre el horizonte. El tiempo
transcurrido entre dos culminaciones superiores del Sol define el da solar verdadero, que
es un poco mayor que el da sidreo como consecuencia de que la Tierra tiene que
avanzar en su movimiento sobre la eclptica hasta enfrentar al Sol.


3.1.5 Tiempo solar y tiempo oficial. La ecuacin del tiempo

El tiempo solar verdadero, TSV, en un lugar dado es, por definicin, el ngulo
horario del centro del Sol expresado en horas. Cada hora equivale a 15, luego:

15

= TSV [3.7]

Por tanto, el da solar verdadero est dividido en 24 horas de tiempo solar
verdadero, que empieza a contarse a partir del medioda.

El tiempo local aparente, TLA, es, por definicin, el tiempo solar verdadero
aumentado en 12 horas, por lo que empieza a contarse a partir de medianoche.

TLA = TSV + 12 [3.8]

El da solar verdadero, sin embargo, tiene una duracin desigual a lo largo del
ao. El movimiento aparente del Sol sobre el horizonte es la composicin de dos
movimientos, uno diario, debido al movimiento de rotacin de la Tierra, que se puede
suponer perfectamente uniforme y otro anual, a causa del movimiento orbital terrestre,
que no es uniforme, debido a la excentricidad de esta rbita. As pues, el tiempo solar
verdadero, que viene determinado por la observacin del Sol, tiene una naturaleza ms
geomtrica que fsica, ya que no es una funcin lineal del tiempo, y no sirve por
consiguiente, para establecer una escala de tiempo uniforme.

Corrigiendo todas las irregularidades del da solar verdadero se obtiene lo que se
denomina da solar medio, de duracin constante, que es el utilizado normalmente en la
vida diaria. El tiempo solar medio, TSM, en un lugar dado es, por definicin, el tiempo
solar verdadero corregido de todas sus irregularidades. Definiremos como ecuacin del
tiempo, E
t
, la suma algebraica de todas las correcciones que es preciso deducir del
tiempo solar verdadero para despejarlo de todas sus irregularidades. Se tiene pues, por
definicin:

TSV = TSM + E
t
[3.9]

La figura 3.4 representa grficamente la ecuacin del tiempo a lo largo del ao.
Alcanza un valor mximo de 16 min. 24 s el da 4 de noviembre, un valor mnimo de
-14 min. 22 s el da 11 de febrero y se anula cuatro veces al ao, los das 16 de abril, 15
Radiacin Solar Sara Moreno Tejera

16
de junio, 2 de septiembre y 25 de diciembre. Estas fechas y valores son tan slo
aproximadas y oscilan ligeramente de un ao a otro.

-20
-15
-10
-5
0
5
10
15
20
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
Da Juliano
E
c
u
a
c
i

n

d
e
l

T
i
e
m
p
o

[
M
i
n
u
t
o
s
]

Figura 3.4 Representacin grfica de la ecuacin del tiempo


El tiempo civil, TC, tambin denominado tiempo local medio, es, por definicin,
el tiempo solar medio aumentado en 12 horas. Se empieza a contar, por consiguiente, a
medianoche, instante en el que se cambia la fecha del da.

TC = TSM + 12 [3.10]

As pues, dos lugares de la Tierra situados en meridianos distintos tendrn un
tiempo civil diferente, siendo esta diferencia proporcional a su diferencia en longitud
geogrfica. Para sincronizar estos dos lugares se utiliza el tiempo universal, TU, definido
como el tiempo civil del meridiano 0 o de Greenwich
1
, que no debe confundirse con el
tiempo solar medio del meridiano 0, TMG, ya que este ltimo se cuenta a partir del
medioda. Por lo tanto, para un lugar de longitud geogrfica (positiva hacia el Este), el
tiempo civil en horas viene dado por:

TC = TU +
15

[3.11]

El carcter local del tiempo civil impide su utilizacin oficial en la vida prctica.
Para obviar este inconveniente y extender a todo el mundo el tiempo universal, se ha
dividido la superficie terrestre en 24 husos horarios, con una amplitud de 15 (1 hora) de
longitud cada uno, numerados de 0 a 23 hacia el Este, a partir del huso horario 0, que se
toma como referencia. Este huso tiene como meridiano central el 0, extendindose, por
tanto, entre las longitudes +7.5 y -7.5. Por convencin arbitraria dentro de cada huso

1
No siempre ue ese meridiano el que se emple como reerencia. lubo pocas histricas en que lo ue el que pasaba por
Salamanca y, con posterioridad el que pasaba por San lernando ,Cadiz, porque all se encontraba ,y se encuentra todaa,
el Obseratorio Astronmico de mayor importancia en su poca.
Radiacin Solar Sara Moreno Tejera

17
rige el tiempo local estndar, TLE, definido como el TU aumentado en un nmero entero
de horas igual al nmero del huso. As pues, el tiempo local estndar es el tiempo civil
del meridiano central de dicho huso. Para una determinada localidad, de longitud
geogrfica , la diferencia entre el tiempo civil y el tiempo local estndar, conocida
como correccin de longitud, viene dado por

TC TLE =
( )
15
S

[3.12]

donde
S
es la longitud del meridiano central del huso correspondiente. La
correccin de longitud es positiva si el lugar est al Este de dicho meridiano y negativa
en caso contrario.

El tiempo oficial, TO, que es el que marcan los relojes que usamos en la vida
ordinaria, depende de la normativa legal de cada pas, pero, normalmente, difiere un
nmero entero de horas del correspondiente al huso horario. En este caso

TO = TLE + AO [3.13]

siendo AO el adelanto oficial sobre el huso horario, que puede variar a lo largo
del ao (en Espaa se les suele denominar, horarios de invierno y verano y toman los
valores 1 y 2 horas respectivamente). Por consiguiente, de las definiciones dadas se
tiene que tiempo solar verdadero y tiempo oficial estn relacionados por

TSV = TO AO +
( )
12
15
+

t
S
E

[3.14]

donde al igual que en las relaciones anteriores, el tiempo est expresado en
horas, lo mismo que E
t
y en grados (positiva hacia el Este).

El tiempo solar verdadero se cuenta, como se ha visto, a partir del medioda. Sin
embargo, en la bibliografa ms comn es frecuente contarlo a partir de medianoche. Es
decir, se toma como tiempo solar verdadero el tiempo local aparente, TLA.

Finalmente, habra que decir que es frecuente en el estudio de la radiacin solar
agregar sta en perodos horarios. En este caso se entiende por hora 1 la comprendida
entre medianoche verdadera y las 1 horas TLA, por hora 2 la comprendida entre las 1
horas y las 2 horas TLA, y as sucesivamente. La hora 24 sera entonces la comprendida
entre las 23 horas TLA y la medianoche verdadera del da siguiente.

3.1.6 Clculo de la posicin solar. Ecuaciones aproximadas

La duracin del ao trpico es de 365.24219 das y no de un nmero entero de
das. Para obviar este inconveniente sin producir un corrimiento secular apreciable de
las estaciones a lo largo del ao, el calendario gregoriano, de uso casi universal en la
actualidad, y en el que son bisiestos los aos mltiplos de cuatro, excepto los mltiplos
de cien que no son divisibles por cuatrocientos como 1700, 1800, 1900, 2100, etc.,
utiliza aos de duracin desigual. Por consiguiente, los instantes de paso del Sol por
Radiacin Solar Sara Moreno Tejera

18
equinoccios, solsticios, perihelio y afelio, as como las variables utilizadas para
determinar la distancia y posicin del Sol en un instante determinado (distancia Tierra-
Sol, declinacin solar, ecuacin del tiempo, etc.) varan ligeramente de un ao a otro.

La utilizacin del calendario gregoriano, tanto en el clculo de la posicin y
distancia solar como en el anlisis estadstico de los datos de radiacin solar, resulta
poco adecuada, por lo que en su lugar suelen emplearse otros mtodos.

Un mtodo habitual consiste en utilizar nicamente aos de 365 das,
prescindiendo del da 29 de febrero de los aos bisiestos, que es ignorado. Este mtodo
conlleva un pequeo error, que sin embargo es aceptable en la mayora de las
aplicaciones, y tiene la ventaja de considerar todos los aos idnticos, en lo que a
posicin y distancia solar se refiere, simplificando en gran medida los clculos
necesarios. Cada da de este ao se identifica mediante un nmero denominado da del
ao, n, que varia entre 1 para el da 1 de enero y 365 para el da 31 de diciembre.

Suponiendo que durante un determinado da del ao la declinacin solar es
constante (la mxima variacin en un da ocurre durante los equinoccios y es menor de
30'), se tendra que la trayectoria recorrida por el Sol ese da es exactamente un arco de
paralelo celeste. La posicin solar en un instante determinado vendr dada, en
coordenadas horarias, por el ngulo horario c
i
y la declinacin o del centro del Sol. El
ngulo horario (expresado en grados) es:

c = 15 t [3.15]

donde t est expresado en horas de tiempo solar verdadero y la declinacin solar
se calcula mediante la expresin trigonomtrica propuesta por Bourges, que ofrece un
error mximo de 12`, y que iene dada por:

t t t
t t t
cos 0.0201 sen 0.1712 - cos 0.3656 ...
... sen 0.1149 cos 0.7580 - sen 23.2567 0.3723


+ +
+ + + =
[3.16]

donde o est expresado en grados as como c
t
puede calcularse mediante:

( ) 436 . 79
365
360
= n
t
[3.17]

La ecuacin del tiempo expresada en minutos, puede calcularse por medio de la
siguiente expresin, debida a Spencer, y presenta un error mximo de unos 35
segundos:

( ) ( ) [
( ) ] ) 2 04089 . 0 2 cos 014615 . 0 ...
... 032077 . 0 cos 001868 . 0 000075 . 0
60
2 . 229

+ =
( sen
sen E
t
[3.18]


en la que I ya se ha definido en la expresin 3.4.


Radiacin Solar Sara Moreno Tejera

19
3.1.7 Posicin del Sol relativa a una superficie plana

La posicin relativa del Sol con respecto a una superficie se establece en funcin
de la posicin solar en la esfera celeste y de la posicin de la superficie sobre el plano
del horizonte, la cual viene caracterizada por los dos ngulos:

El ngulo acimutal, y, tambin llamado orientacin, que forma la proyeccin
de la normal a la superficie sobre el plano horizontal con el plano meridiano
del lugar. Se cuenta a partir del punto cardinal Sur, de 0 a 180,
positiamente hacia el Oeste y negatiamente hacia el Lste.

La inclinacin o pendiente, , de la superficie, que es el ngulo formado por
dicha superficie con el plano horizontal. Se mide de 0 a 180. Si > 90 la
superficie estara enfrentada al suelo.

La irradiacin solar recibida por una superficie est influenciada por el ngulo
de incidencia, i, que forman los rayos del Sol con la normal a la superficie. Este ngulo,
para una superficie arbitrariamente orientada e inclinada, viene dado en funcin de las
coordinas horizontales del Sol, por:

cos i = cos sen o + sen cos o cos ( , y ) [3.19]

Otra expresin ms compleja, pero que determina directamente el ngulo de
incidencia en funcin de las coordinas horarias del Sol y de la latitud geogrfica.
Utilizada para en la implementacin informtica. Viene expresada por:

cos i = sen sen o cos o sen c + sen cos o cos o sen v cos c -
- sen cos o sen o cos v + cos sen o sen v + cos cos o cos v cos c

[3.20]

Figura 3.5 ngulo de incidencia de los rayos solares sobre una superficie
Radiacin Solar Sara Moreno Tejera

20

Si la superficie est orientada al Sur ( o = 0 ) y es horizontal ( = 0 ), el
ngulo de incidencia sera el ngulo cenital, y se tendra la siguiente igualdad:

cos i = cos 0
z
= sen o = sen v sen o + cos v cos o cos c [3.21]

3.1.7 Radiacin solar extraterrestre

La radiacin solar extraterrestre determina un valor mximo terico de la energa
solar disponible, por lo que es ampliamente utilizada como referencia en el estudio de la
radiacin solar. Como patrn bsico se toma la constante solar, I
CS
, definida como la
cantidad total de energa procedente del Sol, en todas las longitudes de onda, por unidad
de tiempo y por unidad de rea de una superficie normal a los rayos solares y a la
distancia media entre la Tierra y el Sol. Medidas recientes de la constante solar (segn
el WRC de Davos, Suiza) dan un valor de:

I
CS
= 1367 W m
-2
[3.22]

La irradiacin extraterrestre incidente sobre una superficie cualquiera en un
instante dado, viene establecida por:

I
0
= I
CS
E
0
cos i [3.23]

Esta magnitud es utilizada en el caso de superficies horizontales, luego la
expresin ms til se tiene cuando el ngulo de incidencia viene dado por el ngulo
acimutal, por lo que la irradiancia extraterrestre horizontal es:

I
0
= I
CS
E
0
cos 0
z
= I
CS
E
0
sen o [3.24]

donde cos 0
z
o sen o se determinan con la expresin 3.21.

Y expresando la posicin solar en coordenadas horarias (ngulo horario y
declinacin ), queda

) cos cos cos (
0 0
+ = sen sen E I I
CS
[3.25]

La irradiacin, o cantidad total de energa recibida durante un determinado perodo de
tiempo [t1, t2], se obtiene integrando la irradiancia para ese perodo. Para integrar se
tomarn nicamente intervalos de integracin dentro de un determinado da del ao. De
esta forma, tanto como E
0
pueden suponerse constantes durante el perodo [t1, t2] y,
por consiguiente, I
0
depende nicamente del ngulo horario . En este caso


= =
2
1
0
2
1
0 0
12

d I dt I H
t
t
[3.26]

Ya que ngulo horario, en radianes, y tiempo solar verdadero, en horas, estn
relacionados por

Radiacin Solar Sara Moreno Tejera

21

=
12
2
24
t [3.27]

Sustituyendo e integrando, resulta

)) ( cos cos ) ((
12
1 2 1 2 0
0

sen sen sen sen E I H

CS
h
+ = [3.27]

Donde los ngulos horarios inicial
1
y final
2
del intervalo considerado estn
expresados en radianes. Esta expresin permite calcular la irradiacin extraterrestre
sobre una superficie horizontal para un determinado da del ao n, caracterizado por y
E
0
, y para un determinado perodo de tiempo [t
1
, t
2
] en ese da.

La irradiacin extraterrestre horizontal horaria se obtiene integrando sobre un perodo
de una hora. As pues, si
1
es el ngulo horario a mitad de dicho perodo, haciendo

24
2

+ =
i

24
1

=
i
[3.28]

Y sustituyendo se obtiene que

) cos cos cos
24
24
(
0
0
i CS
h
sen sen sen E I H

+ = [3.29]

como

1 9972 . 0
24
24
=

sen

Se puede tomar, con un error mnimo

) cos cos cos (
0
0
i CS
h
sen sen E I H + = [3.29]

Anlogamente, la irradiacin extraterrestre horizontal diaria se obtiene integrando
paratodo el da, entre el orto y el ocaso solar. As pues, haciendo

s
=
2

s
=
1


resulta

) cos cos cos
180
(
24
0
0
i
S
CS
d
sen sen E I H

+ = [3.30]

donde
s
es el ngulo horario del ocaso solar, en grados, para el da considerado.


Radiacin Solar Sara Moreno Tejera

22
3.2 Componentes de la irradiacin solar a nivel de suelo

La interaccin de la radiacin solar con la atmsfera, hace que se distingan dos
componentes: radiacin solar directa y radiacin solar difusa, de manera que el clculo de
la radiacin global incidente sobre una superficie inclinada a nivel de suelo, implica la
determinacin de tres magnitudes, como se observa en la figura 3.6, y que son:

- Radiacin directa desde el disco solar
- Radiacin difusa procedente del cielo
- Radiacin difusa reflejada desde el suelo


Figura 3.6. Componentes de la radiacin global incidentes sobre una superficie inclinada


La tarea de calcular la irradiacin solar directa sobre una superficie inclinada, es en
parte un problema geomtrico y en parte un problema fsico determinado por la interaccin
de la radiacin solar con la atmsfera. Para intervalos de tiempo cortos, del orden de una
hora o menos, se puede suponer constante la altura solar sin introducir un error elevado,
por lo que se puede determinar la irradiacin directa sobre la superficie en funcin de la
irradiacin directa normal y del ngulo de incidencia de sta con la superficie.

El clculo de la irradiacin difusa incidente en superficies inclinadas es, si no
hacemos ciertas simplificaciones y modelizaciones, de elevada complejidad debido a la
gran variabilidad tanto en el tiempo como en su distribucin espacial sobre el cielo, de la
radiacin difusa.

Existen en la bibliografa diversos modelos de distribucin de radiacin difusa. Se
basan tanto en distribuciones istropas sencillas, como en otras ms complejas, anistropas
y multiparamtricas de la radiacin difusa, y que se estudiarn detalladamente en secciones
posteriores.

Radiacin Solar Sara Moreno Tejera

23
Como resultado de la carencia de datos medidos en los que basar algoritmos ms
perfeccionados, es usual aceptar la simplificacin de obtener la irradiacin difusa reflejada
desde el suelo asumiendo de modo general que la reflexin es de naturaleza istropa,
semejante a la producida desde una superficie horizontal infinita frente a la superficie en
cuestin, si la radiacin global est compuesta principalmente de radiacin difusa o el
suelo es un reflector perfecto.

Se observa pues, que la determinacin de la irradiacin global incidente sobre una
superficie inclinada, es compleja y pasa primeramente por la determinacin de sus tres
componentes, de manera independiente.

3.3 Modelos de Radiacin

La cuantificacin de la radiacin solar incidente sobre cualquier superficie, tiene
gran importancia para el diseo adecuado de instalaciones de aprovechamiento de la
energa solar.

Aunque existen algunas estaciones que registran datos de irradiacin global
horizontal as como de horas de sol, no puede decirse lo mismo del registro de valores de
irradiacin difusa o directa. Por lo que es necesario, desarrollar algoritmos que estimen
alguna de estas componentes a partir de los valores de irradiacin global sobre superficie
horizontal.

Algunos investigadores con el fin de solucionar este problema, han desarrollado
relaciones entre la irradiacin difusa y global horizontal, para diversos perodos de tiempo
y basada en coeficientes y expresiones sencillas.

Previamente se han de definir algunos de estos coeficientes, en los que se basarn
las relaciones que se van a presentar, como son:

- K
t
: cociente entre el valor de la irradiacin global horizontal y la irradiacin
extraterrestre horizontal, en perodos horarios, diarios, etc.
- K
d
: cociente entre el valor de la irradiacin difusa horizontal y la irradiacin global
horizontal, en perodos horarios, diarios, etc.
- K
do
: cociente entre el valor de la irradiacin difusa horizontal y la irradiacin
extraterrestre horizontal, en perodos horarios, diarios, etc.


3.2.2 Relaciones de tipo horario

Entre las relaciones K
t
- K
d
horarias ms conocidas y a su vez de uso ms
extendido, cabe destacar las de Orgill - Hollands y M. Iqbal, por ser las que mejores
estimaciones proporcionan.
Asimismo, se presentan otras relaciones de tipo K
t
- K
do
, desarrolladas para Sevilla
y Madrid basadas en medidas de valores horarios de irradiacin global horizontal, directa
normal y difusa horizontal, registrados desde 1984 en la Escuela Tcnica Superior de
Ingenieros Industriales.

Radiacin Solar Sara Moreno Tejera

24

3.2.2.1 Algoritmo de Orgill-Hollands

La relacin entre K
d
y K
t
horarios viene expresada por medio de tres ecuaciones:

[3.25]

La base de datos original corresponde a cuatro aos de medidas de valores horarios
de radiacin difusa y global horizontales en Toronto (Canad) en el perodo Septiembre-
1967 / Agosto-1971. A partir de estas medidas, se obtuvieron los parmetros K
t
y K
d
ya
definidos. Dividiendo el eje de abscisas (eje K
t
) en intervalos de tamao 0.05, para cada
uno de ellos se promediaron los valores de K
d
all contenidos, con lo que se obtenan los
pares de puntos (K
T
, K
D
) representativos de cada intervalo y donde:

K
T
: punto medio del intervalo de K
t
considerado.
K
D
: valor medio de los puntos K
d
contenidos en el intervalo.

El rango 0 K
t
0.35, comprende el 32.4 % de los datos medidos. Los valores
dentro de este intervalo son representativos de das extremadamente cubiertos, con
alrededor de un 90 % de radiacin total de tipo difusa. A causa del alto porcentaje de sta,
los valores medidos de radiacin global tienden a ser pequeos, y por consiguiente
afectados por la sensibilidad y precisin de los instrumentos de medida.

En el rango 0.35 < K
t
0.75, se incluyen el 62 % de los datos medidos.

Para los valores de K
t
> 0.75, se asigna valor constante a K
d
. Esto es debido a que
este intervalo es representativo de das claros con alguna pequea nubosidad, aunque el sol
en s mismo no est sombreado por nubes. En este caso, la radiacin se refleja en ellas en
cantidades sustanciales, actuando stas como concentradores de radiacin difusa, que no
tiene la direccin de la radiacin directa. Debido al pequeo nmero de datos en este
intervalo, as como a lo impredecible de esa reflexin, se aconseja usar un valor constante,
igual a 0.177.

Estas expresiones se recomiendan para determinar la radiacin difusa horaria en
superficies horizontales, esperndose con ellas una suficiente precisin en latitudes
comprendidas entre 43 N y 54 N.

3.2.2.2. Relaciones de M. Iqbal

El resultado de la relacin de Iqbal fue presentado de una manera grfica en
funcin de K
do
, K
t
y la altura media del sol correspondiente a una hora solar determinada.
Estas relaciones, modelizadas analticamente para utilizarlas en clculo por ordenador,
pueden expresarse en la forma:



d t t
d t t
d t
K
= 1.0 - 0.249
K
0
K
0.35
K
= 1.577 - 1.84
K
0.35
K
0.75
K
= 0.177 0.75 <
K
si
si
si



Radiacin Solar Sara Moreno Tejera

25
Altura solar de 10 grados

[3.26]

donde
D = 0 si K
t
0.275
D = 1 si K
t
> 0.275


Altura solar de 20 grados

[3.27]

donde
D = 0 si K
t
0.325
D = 1 si K
t
> 0.325


Altura solar de 30 grados

donde
D = 0 si K
t
0.325
D = 1 si K
t
> 0.325


Altura solar de 40 grados

donde
D = 0 si K
t
0.35
D = 1 si K
t
> 0.35




do t t t
2
t
3
t
4
K
= 0.935
K
(1- D)+ D(-0.58+6.43
K
- 16.75
K
+17.87
K
- 6.42
K
)



do t t K
= 0.91
K
(1- D)+ D(0.258+0.067 (553
K
- 236)) cos



do t t K
= 0.89
K
(1- D)+D(0.253+0.083 (553
K
- 243)) cos


[3.28]

do t t K
= 0.87
K
(1- D)+D(0.254+0.083 (553
K
- 253)) cos


[3.29]

Radiacin Solar Sara Moreno Tejera

26
3.2.2.3. Relaciones obtenidas con datos horarios de Sevilla. Algoritmos de la
media, mediana y moda estadsticas

Con datos medidos durante cuatro aos en la Escuela Tcnica Superior de
Ingenieros Industriales de Sevilla, se han obtenidos diversas relaciones entre K
t
y K
do
, que
no son funcin de la altura solar.

Se ha utilizado un mtodo inspirado en el algoritmo de Orgill-Hollands
desarrollado para Toronto, con la diferencia de de emplear el parmetro K
do
y calcular para
cada intervalo de K
t
, adems del valor medio de aquel, los valores medianos y modales.
Para su implementacin analtica, se realiz un ajuste en serie de Fourier en la forma:


donde P(x) es un desarrollo en serie de Fourier, con uno o dos armnicos como
mximo y que viene dado por:

[3.31]

Finalmente se obtuvieron tres expresiones denominadas segn la variable
estadstica empleada, algoritmos de la media, mediana y moda, cuyas formulaciones son:


Algoritmo de la Media.
[3.32]

donde
D = 0 si K
t
0.275
D = 1 si K
t
> 0.275


Algoritmo de la Mediana


Y = (1- D) a x + D P(x)


[3.30]

P(x) =
A
+
A
x +
B
sen x +
A
2 x +
B
sen2 x =
=
A
+
C
( x - ) +
C
(2 x - )
0 1 1 2 2
0 1
1
2
2
cos cos
cos cos




do t t K
= 0.93
K
(1- D)+ D(0.2562+0.076 (550
K
- 242.6)) cos



do t t K
= 0.93
K
(1- D)+ D(0.2479+0.0954 (550
K
- 235.6)) cos


Radiacin Solar Sara Moreno Tejera

27
[3.33]

donde
D = 0 si K
t
0.275
D = 1 si K
t
> 0.275


Algoritmo de la Moda

[3.34]

donde
D = 0 si K
t
0.225
D = 1 si K
t
> 0.225


Estos algoritmos pueden utilizarse para valores de K
t
comprendidos entre 0 y 0,8.
De ellos, los denominados Mediana y Moda consiguen estimaciones levemente mejores
que la Media. Incluso, cuando se comparan con el algoritmo de Orgill - Hollands, obtienen
resultados ms parecidos a las medidas reales que ste, para los datos medidos en Sevilla,
aunque es significativo que el algoritmo de Orgill - Hollands ha obtenido diferencias en la
estimacin no demasiado importantes en relacin a los otros tres. Asimismo,
comparndolas con las relaciones de M. Iqbal, las obtenidas para Sevilla indican que, salvo
para la altura de 10, el considerar la altura del sol como parmetro en la correlacin, no es
significativo.

3.2.2. Relaciones de tipo diario

Para valores diarios, son varias las relaciones Kd-Kt que se pueden destacar de la
bibliografa. Aqu se presentan las relaciones obtenidas por Collares, Ruth-Chant, Muneer ,
Liu-Jordan y Frutos- Ruiz.
Salvo la relacin de Liu-Jordan (por otra parte la primera y que dio origen a las dems), las
dems proporcionan resultados semejantes en el intervalo de Kt comprendido entre 0.25 y
0.75, que a su vez constituye el de mayor porcentaje de datos registrados.
Las relaciones vienen expresadas de la siguiente forma:

Collares
[3.35]

do t t K
= 0.93
K
(1- D)+D(0.2205+0.0899 (520
K
- 215)) cos



d t
d t t
2
t
3
t
4
t
d t t
d t
K
= 0.99 ,
K
0.17
K
= 1.188 - 2.272
K
+9.473
K
- 21.865
K
+14.648
K
, 0.17 <
K
0.75
K
= - 0.54
K
+0.632 , 0.75
K
<0.8
K
= 0.2 , 0.8
K
si
si
si
si

Radiacin Solar Sara Moreno Tejera

28

Ruth-Chant

[3.36]

Muneer
[3.37]

Liu-Jordan
[3.38]

Frutos-Ruz
[3.39]



Por ltimo, incluimos los algoritmos que se obtuvieron con los datos recogidos por la
estacin radiomtrica de la Escuela Tcnica Superior de Ingenieros:





K
0.8 , 0.1323 =
K
0.8
K
< 0.1 ,
K
2.848 +
K K
1.15 0.91 =
K
0.1
K
, 0.98 =
K
t d
t t t t d
t d

si
si 936 . 4 4
si
3 2


K
0.77 , 0.16 =
K
0.77
K
< 0.17 ,
K
2 +
K
3.62 -
K
0.47 + 1.024 =
K
0.2
K
, 0.98 =
K
t d
t
3
t
2
t t d
t d

si
si
si


K
0.8 , 0.1319 =
K
0.
K
< 0.17 ,
K
2.631 +
K
3.522 -
K
0.0435 + 1.0045 =
K
0.2
K
, 0.94 =
K
t d
t
3
t
2
t t d
t d

si
8 si
si


d t
d t t
2
t
3
t
d t
K
= 0.924 ,
K
0.22
K
= 0.68+2.47
K
-6.955
K
+3.53
K
, 0.17 <
K
0.75
K
= 0.112 , 0.75
K
si
si
si

Radiacin Solar Sara Moreno Tejera

29

Para el perodo 1984-2000:

kd =1 si kt 0.22
kd = 1.390975-1.69912kt si kt >0.22

[3.40]


Para el perodo 2000-2003:

kd =1 si kt 0.22
kd = 1.384304-1.700599kt si kt >0.22

[3.41]