CALCULOS DE LA CINEMATICA SOLAR

1.- Cálculo de la irradiancia solar sobre una superficie

Consideremos por ahora el cálculo de la radiación directa sobre una superficie.

El flujo de energía sobre una superficie determinada depende no sólo de la irradiancia
que exista, sino también de la orientación que tenga la superficie en cuestión respecto de
la dirección de propagación de la radiación. La irradiancia será máxima sobre un plano
que esté en posición normal a la dirección de propagación de la radiación (es decir, que
le dé "de frente" ). En cambio, será mínima si el plano es paralelo a la radiación (es
decir "de canto")

La radiación extraterrestre que incide sobre la Tierra está sujeta a las variaciones
geométricas y a las condiciones físicas del propio Sol.

Por otro lado, la órbita que describe la Tierra alrededor del Sol no es circular,
sino cuasielíptica. La pequeña excentricidad de la órbita hace que, alrededor del 4 de
enero, cuando la Tierra se encuentra en el perihelio (mínima distancia al Sol -sí:
mínima-) la radiación solar extraterrestre sea máxima. Por otro lado, alrededor del 1 de
julio -seis meses después- la Tierra se encuentra en el afelio (máxima distancia al Sol) y
entonces la radiación solar extraterrestre es mínima.

Para el estudio de equipos que hacen uso de la potencialidad solar es necesario

encontrar los valores de la radiación por hora para ello se utiliza el método de Duffie –
Beckman para encontrar el comportamiento de la radiación solar a lo largo del día, para
ello es necesario tomar en cuenta los siguientes elementos necesarios para la
compresión del procedimiento.

 (Phi) Simboliza la latitud geográfica del lugar que toma valores positivos hacia el
norte de la línea ecuatorial y negativos hacia el sur ( - 90º  º ).

La ciudad de loja tiene una latitud Sur de 4º por lo que el valor de  = - 4º

Z


t
N
W

s
S  E

Figura 1. Ángulos formados por el sol con la horizontal

 (Delta) Declinación solar. La declinación es máxima en el solsticio de verano (23º
27') y mínima en el solsticio de invierno (-23º 27'). Puede calcularse
aproximadamente mediante la ecuación:
284 +n
(
δ=23 . 45 sen 360
365 )
(1.1)

 (Beta) Pendiente del plano sobre el cual incide la radiación, esto es, su inclinación
respecto de la horizontal.
Si se tratase de un destilador podría =10º
 (Gamma) Acimut de la superficie, ie, el ángulo que forma la normal a la superficie,
respecto del meridiano local.

El valor de  se lo considera cero por que el destilador estará frente al sur

 (Omega) Angulo horario, esto es, el desplazamiento angular del Sol, al este o al oeste
del meridiano local, debido a la rotación de la Tierra a la velocidad angular de 15 o
por hora (360º en 24 horas). Negativo por la mañana y positivo por la tarde

Tabla 1. Ecuaciones para convertir el día del mes, en el número de día del año.

“n” para el i-esimo

Mes
día del mes
Enero i
Febrero 31 + i
Marzo 59 + i
Abril 90 + i
Mayo 120 + i
Junio 151 + i
Julio 181 + i
Agosto 212 + i
Septiembre 243 + i
Octubre 273 + i
Noviembre 304 + i
Diciembre 334 + i

1.1 Determinación de la Declinación Solar por ejemplo para el día 10 del mes de
Junio () y energía incidente.

La declinación solar es cero el día que comienza la primavera y el día que

comienza el otoño (equinoccios). En estos días, el Sol se encuentra en el ecuador
celeste, que es la extrapolación del ecuador terrestre en la bóveda celeste. La
declinación es máxima en el solsticio de verano (23o 27') y mínima en el solsticio de
invierno (-23o 27'). La declinación solar puede calcularse aproximadamente mediante la
ecuación

De acuerdo con la tabla 1 el día 10 de Junio es el día 161 del año, entonces: n = 161

284+n
δ=23 , 45 ⋅sen 360
( ( 365 ))
(1)

284+161
δ=23 , 45 ⋅sen 360
( ( 365 ))
 = 23,012 º
1.2 Angulo horario de la salida del sol (Ws)

−1
W s=cos (−tg φ⋅tg δ )
(2)
−1
W s=cos (−tg ( − 4 )⋅tg (23 , 012))
−1
W s=cos ( 0, 029699 )
Ws = 88,298º

1.3 Determinación de la radiación solar extraterrestre sobre una superficie normal

al sol (Ho)

)) ( )
24⋅3600⋅Gsc π ⋅W s
Ho =
π (
1+ 0, 033 cos
360⋅(n )
365( ⋅ cos φ⋅cos δ ⋅senW s +
180 (
senφ ⋅senδ )
(3)
Donde:
Gsc es la Constante Solar y tiene un valor de 1367 W/m2
Luego tenemos que:

Ho =
( 24 )⋅( 3 600 )⋅( 1 367 )
π (
1 + 0, 033 cos
360⋅( 161 )
365 ( ⋅¿ ))
¿

(
¿ cos ( − 4 )⋅cos ( 23 ,012 )⋅sen ( 88 ,298 ) + (180π⋅( 88 , 298 ) )⋅sen ( − 4 )⋅ sen ( 23 , 012 ))
7
Ho=3 .191 x 10 [ J /m2]
 3
Ho=31 ,91 x 10 [ kJ /m2]

1.4 Cálculo de la radiación total diaria por unidad de área de una

superficie horizontal (H')

H ' = K T⋅Ho
(4)
Donde:
KT = Índice de Nubosidad o Factor de Turbiedad de la atmósfera y tiene en cuenta la
transparencia de la atmósfera por efectos de nubes, polvo, aerosoles, ozono, etc.).

12204
Kt=
Ho

Kt = 0.382

H '=0,382 ⋅ 31 ,91 x 10 3

H’= 1,22 x 10 4 [kJ/m2]

1.5 Cálculo de la Radiación Difusa Diaria que cae sobre una superficie
horizontal (Hd’)

H d ' = f (K T )⋅H '

(5)

Donde f(KT) viene dada por la siguiente expresión:

f (K T ) = 1,3903−4, 0273⋅( K T ) + 5, 5315⋅( K T )2 −3, 1080⋅( K T )3

(6)

2 3
f (K T )=1, 3903−4, 0273⋅(0,349 )+ 5,5315⋅(0, 349) −3, 1080⋅( 0, 349 )

f (K T )=0, 485

Ahora bien
 −4
H d '=0, 485 ⋅ (1, 22 x 10 )

Hd’ = 5.923 x 103 [kJ/m2]

1.6 Determinación de la Radiación Solar Total Horaria por unidad de área que cae
sobre una superficie horizontal (H)

H=K⋅H '

(7)

Siendo:

K, la relación entre la Radiación Solar Horaria y la Radiación Solar Total Diaria y se

calcula por la siguiente expresión tomada de Duffie (1991):

( cosW−cos W s )

( )
π
K= ( )
24
⋅( a+b⋅cosW )⋅
( π⋅W s )
senW s − ( 180 ) ⋅cosW s

Donde

a = 0,409 + 0,5016 sen (Ws – 60)

a = 0,409 + 0,5016 sen ( 91.702 – 60)
a = 0.647
b = 0,6609 – 0,4767 sen (Ws – 60)
b = 0,6609 – 0,4767 sen ( 91.702 – 60 )
b = 0.435

W = Angulo Horario que toma valor 0 para el mediodía solar y varía 15 cada hora.

W = 75º, W = 60º, W = 45º, W = 30º, W = 15º, W = 0º

Para W = 75º

( cos 75−cos 88, 298 )

( )
π
K= ( )
24
⋅( 0,647+0,435⋅cos 75 )⋅
sen 88 ,298−
180(
( π⋅88,298 )
)
⋅cos 88,298

K= 0,024
 Como H = K.H’

H = 0,024 (1,22 x 104)

H = 291.401

Se realiza el mismo procedimiento para calcular los valores de k para los demás
ángulos y luego se calcula la Radiación Solar Total horaria por unidad de área sobre una
superficie horizontal (H ).

Tabla 2. Radiación Solar Total Horaria por unidad de Área

Horas del día W K H (W/m2 h)

6-7
75o 0,024 291,401
17 - 18
7-8
60o 0,056 680,77
16 - 17
8–9
45o 0,089 1,038 x 103
15 – 16
9 – 10
30o 0,117 1,434 x 103
14 - 15
10 – 11
15o 0,137 1,673 x 103
13 – 14
11 – 12
0o 0,144 1,758 x 103
12 - 13

1.7 Determinación de la Radiación Difusa Horaria por unidad de área que cae
sobre una superficie horizontal (Hd).

En la Página 91 del Duffie (1991) aparece la ecuación de la relación entre la

radiación difusa horaria y la radiación difusa diaria sobre una superficie horizontal
como función de la longitud del día. De esta ecuación se obtienen los valores de K d y
se calcula la radiación difusa horaria.
H d =K d⋅H d ' (8)

Los valores de Kd fueron encontrados por la siguiente formula

 ( cosW−cosW s )

( )
π
K d=
24 ( )
⋅
( π⋅W s )
senW s− ( ) 180
⋅cosW s

(9)

( cos75 º−cos 88,298 º )

( )
π
K d= ( )
24
⋅
sen 88,298º− (
( π⋅88 ,298 º )
180 )
⋅cos88 ,298º

Kd = 0,031
3
H d =0,031⋅ 5,923 x 10

Hd = 186.231

Con cada uno de los valores horarios de K d se calculan los valores de H d los
cuales se muestran en la tabla 3

Tabla 3. Radiación Difusa Horaria por unidad de área

Horas del día W Kd Hd (W/m2 )

6-7
75o 0.031 186,231
17 - 18
7-8
60o 0.065 382,267
16 - 17
8–9
45o 0.093 550,607
15 - 16
9 – 10
30o 0.115 679,779
14 - 15
10 – 11
15o 0.128 760,98
13 – 14
11 – 12
0o 0.133 788,676
12 - 13

1.8 Determinación del valor del Angulo de Incidencia Horario de la Radiación

Solar Directa sobre el calentador solar ( t).

Estos valores se pueden calcular teniendo en cuenta los términos ( = 4°,  = 10°,  =
23,012º. El Angulo de Incidencia de la Radiación Solar sobre el calentador solar, viene
dado por:
 cosΘ= ( sen δ⋅senφ⋅cos β−sen δ⋅cos δ⋅sen β⋅cos γ+cosδ⋅cos φ⋅cos S⋅cos W + )
( +cos δ⋅sen φ⋅sen β⋅cos γ⋅cos W +cos δ⋅sen β⋅sen γ⋅senW )
(10)

En la mayoría de las aplicaciones de esta ecuación, la única variable a lo largo

del día es el ángulo horario W. Todo lo demás suele ser constante en una fecha y lugar
determinados. Debido a esto, la ecuación tiene un gran número de constantes que, al
aplicarla, simplifican mucho la expresión.

Teniendo en consideración que el ángulo cenital  es 1. Luego aplicando

identidades trigonométricas la expresión nos queda de la siguiente manera:

cos θt =[ cos( φ+ β )⋅cos δ⋅cos W +( sen( φ + β )⋅sen δ ) ]

(11.a)

θt =arccos[cos (φ+ β )⋅cos δ⋅cos W +(sen (φ+ β )⋅sen δ)]

(11.b)

Donde:

 = Angulo de Inclinación de la superficie captadora.

 = Latitud de Loja
 = Declinación Solar para el día medio del mes de Junio
W = Angulo Horario que toma valor 0 para el mediodía solar y varía 15 cada hora.
Los resultados tanto de t como de Cos t son los que se muestran en la tabla
siguiente:

Tabla 4. Angulo de Incidencia Horario

Horas del día W Cos  t t

6-7
75o 0,265 74,638º
17 - 18
7-8
60o 0,486 60,898º
16 - 17
8–9
45o 0,677 47,427º
15 - 16
9 – 10
30o 0,822 34,67º
14 - 15
10 – 11
15o 0,914 23,912º
13 – 14
11 – 12
0o 0,945 19,012º
12 - 13

1.9 Determinación del Angulo Cenital ( z).

 Los valores del ángulo cenital se calculan a partir de la ecuación (10) donde se lo
calculan sobre una superficie horizontal con lo que  = 0 entonces sen = 0 y cos= 1
con lo que desaparecen algunos términos de la expresión. Estos valores se calculan a
partir de la expresión teniendo en cuenta que la cubierta está orientada al norte entonces
el acimut  = 0, por lo que el Angulo Cenital se calcula según la siguiente formula
tomada de la página 16 del Duffie

cos θ z =cos φ⋅cos δ⋅cos W +sin δ⋅sin φ

(12)

Los valores del ángulo Cenital por cada hora son los que se muestran en la tabla 5.

Tabla 5. Angulo Cenital

Horas del
W Cos  z z
día
6-7
75o 0,21 77,856º
17 - 18
7-8
60o 0,432 64,417º
16 - 17
8–9
45o 0,622 51,539º
15 - 16
9 – 10
30o 0,768 39,834º
14 - 15
10 – 11
15o 0,86 30,725º
13 – 14
11 – 12
0o 0,891 27,012º
12 - 13
1.10 Cálculo de la Relación de la Radiación Solar Directa sobre una superficie
inclinada y la Radiación Solar Directa sobre una superficie horizontal (Rb).
Este coeficiente nos da la relación entre la radiación que incide sobre una
superficie inclinada y la que incide sobre una superficie horizontal como se observa en
la figura 2

Gb Gbt
Gbn Gbn 
z


Figura 2. Radiación directa en superficies horizontal e inclinada

Y viene dada por la siguiente expresión:

Ht H n cos θ t
Rb = =
H H n cos θ z (13)

cos θ t
Rb =
cos θ z
(14)

De esta ecuación se obtienen los siguientes valores:

Tabla 6. Valores Rb para cada hora del Día

 Horas del día W Rb
6-7
75o 1,262
17 - 18
7-8
60o 1,125
16 - 17
8–9
45o 1,088
15 - 16
9 – 10
30o 1,07
14 - 15
10 – 11
15o 1,063
13 – 14
11 – 12
0o 1,061
12 - 13

1.11 Determinación de los valores de la Radiación Solar Total Horaria que cae
sobre la superficie del equipo solar (Ht).

La Radiación Solar Total Horaria se calcula según Duffie (1991) página 126 por
la siguiente expresión.

Ht=( H−H d )⋅ Rb + H d ( 1 + 2cos S ) + H⋅ρ ( 1−cos

t
2 )
S

(15)

Siendo t el Albedo o reflexividad de la tierra (t = 0.2) (valor usual).

Tabla 7. Valores de Radiación Solar Total Horaria

Horas del Ht Ht
día (kJ/m2 h) (kcal/m2 h)

6-7
318,324 76,15
17 - 18
7-8
716,886 171.33
16 - 17
8–9
1,128 x 103 269,68
15 - 16
9 – 10
1,485 x 103 355,263
14 - 15
10 – 11
1,728 x 103 413,001
13 – 14
11 – 12
1,815 x 103 433,79
12 - 13