Prediction of the daily global solar irradiance received on a horizontal surface -New approach

International Refereed Journal of Engineering and Science (IRJES)
ISSN (Online) 2319-183X, (Print) 2319-1821
Volume 4, Issue 6 (June 2015), PP.85-93
www.irjes.com 85 | Page
Prediction of the daily global solar irradiance received on a
horizontal surface -New approach
M.K. El-Adawi 1,*
, S.A. Shalaby1
, S.E.-S.Abd El-Ghany2
, F.E.Salman2
and M.A.Attallah2
1
Physics Department –Faculty of Education, Ain Shams University, Heliopolis, Cairo, Egypt:
2
Physics Department –Faculty of Science, Benha University, Benha, Egypt
Abstract:- A general distribution function q (t) W/m2
for clear days is suggested to predict the global solar
irradiance incident on a horizontal surface as a function of the local day time and the maximum daily
value of the incident solar irradiance qmax. Moreover, the authors suggested a form to express qmaxin terms of
the solar constant. The value of the solar constant is also adjusted to its variation with the distance between earth
and the sun.
This is done to get a complete theoretical expression for q (t) rather than semi- empirical formulae introduced in
previous trials. Comparison between the computed values according to the introduced theoretical formula of
q (t) and the corresponding published experimental data recorded in Barcelona(Spain) , Hong Kong (China),
Cairo(Egypt), Makkah ,Jeddah(Saudi Arabia) is made.Good fitting is obtained except for few extreme points.
Keywords:- Prediction of solar irradiance, Daily global solar radiation, Solar distribution function for clear
days, Computation on theoretical basis of qmax..
I. INTRODUCTION
The prediction of the global solar radiation incident on a horizontal surface is important in solar
energy exploitation. For example, the performance of a solar cell is a function of the received incident solar
irradiance( defined as : the radiant energy received by unit area of a surface normal to the pencil of solar rays
per unit time W/m2
) . Several attempts have been made to predict the daily values of the global solar
irradiance q(t) W/m2
, where '' t '' is the local day time in with different degrees of fitting accuracy . (El-
Adawiet.al.,1986) introduced a power expression in terms of │t-tr│ with parameters tr(sunrise), td (the length of
the solar day), ts (sunset). . Good fitting was obtained between experimental and corresponding calculated values
with maximum error ≈11%. In such a trial the authors suggested an expression to estimate the value of qmaxin
terms of a solar constant. (Shalaby, 1994) expressed q(t) as a polynomial in │t-t0│with a certain correction
factor F(t) and parameters tr,ts,t0,td . Where t0=(ts-tr)/2=td/2. Computations revealed that the relative error did not
exceed 16%. (El- Adawiet.al.andNuaim, 2001)introduced an expression of q(t) as a polynomial in( t/td ) with a
correction factor sin(π t/td )with maximum error 15% . (El- Adawiet.al. ,2002)expressed q(t) as a polynomial in(
t/td ) with a correction factor in terms of t0,td. Maximum relative error 12 % was obtained. The evaluation of the
qmax in terms of the solar constant is also suggested in(El- Adawi, 2002). Noticing that the suggestion to estimate
the value of qmax in terms of the solar constant in the previous trials was not tested through any concrete
computations. While semi-empirical formulae are introduced in which qmax is taken as an experimental input.
The present trial represents an approach to the problem of evaluation of the average global solar
irradiance incident on a horizontal surface .The suggested formula is an exponential function in time. The aim of
the trial is to get better fitting with the experimental data. Moreover, qmax in terms of the solar constant is
estimated theoretically in order to have a closed system.
II. THE SUGGESTED MODEL
The experimental measurements for the considered distribution q(t) ,W/m2
for clear sky(El-Adawi et
al.,1986; Shalaby , 1994; Quraishee, 1969 ; Munroe, 1980 ; Tiwari, 1997 ; Leung, 1980 ; Villarrubia et al.
1980; El-Bar, 1983; EL-Gendi, 1983; Katsoulis and Parachistopoulos, 1978; Khogali and Ramadan, 1982;
Harty et al.,1999; Cooper, 1969)show a symmetrical distribution that passes through a maximum value qmax at
the midday time '' t0=td/2 ''between sunrise '' tr'' and sunset '' ts'' in hours .In the present trial ,the distribution
function q(t) ,W/m2
for the hourly daily global irradiance received on a horizontal surface is such that, it satisfies
the following conditions :
at t= tr q(tr)=0 (1)
at t= ts q(ts)=0 (2)

Prediction of the daily global solar irradiance received on a horizontal surface -New approach
at t= t0 q(t0)=qmax. (3)
att= t0 0
|
)(
0



tt
t
tq (4)
The suggested distribution is in the form:
q(t)= qmax
)t-(t)(
)t-(t
r
2
0
tts
e 

(5)
This form satisfies the above mentioned conditions( equations 1-4).
The length of the solar day is given as (Duffieeand Beckman 1974; Cooper, 1969):
td=
15
12
  tantancos 1

(6)
where:
Latitude  , is the angle made by the radial line joining the given location to the center of the earth with its
projection on the equatorial plane.
The solar declination angle  , is the angle between the line joining the centers of the sun and the earth and its
projection on the equatorial plane , is given as:
 =23.45
365
n284
360sin 




 
(7)
(1≤n≤365) starting from 1 January.
Where n is the day of the year.
Following (El-Adawi et al., 1986; El-Adawi, 2002; Wieder, 1982; Kondratyev, 1969 )we can suggest the
estimation of qmax in terms of the solar constant to be in the form:
qmax= s (8)
where:
sis the extraterrestrial solar constant adjusted for the variation of the distance between the sun and the earth
and along the time of the year, given as (Tiwari, 1997)
s = 












 
365
n360
cos0.033+1s (9)
Where:
s=1353(Rai, 1989) is the solar constant .
α=
𝛶+−𝛶−
1+𝐺 𝛶++𝐴−𝐵𝑅
𝑒 𝛶− 𝜏
0 (10)
2
1
2
]4)[(
2
1
)(
2
1
BDACAC 
 (11)
A=
0
0
2
2


(12)
B= 0 (13)
C= 02  (14)
D=
0
0
2

(15)
0 =

 s
where
s
 and  are the optical thickness due to scattering and total optical thickness (scattering
and absorption) .
G=-[
𝛶−+𝐴−𝐵𝑅
𝛶++𝐴−𝐵𝑅
]𝑒(𝛶−−𝛶+)𝜏
(16)
R is the reflectivity of the underlying terrain ,
0 = Zcos (17)

LDZ  (18)
Z is the solar zenith angle and D is the solar codeclinations which is the complementary angle of the
declination, L is the observer colatitudes which is the complementary angle of the latitude.
III. COMPUTATION
In the following qmax is computed according to equation (8) for Jeddah and Makkah only because the
day of the year is available to the authors for these two cities .Computations are made at R=0.4 (Hrens ,
2006), ω0=0.5 (El-Adawiet al.,1986).Moreover,td =18.5hr for Jeddah and td =18.66 hr for Makkah. The obtained
c-omputed values are as follows :
qmax(Jeddah)=856.8W/m2
, while the experimental value is 915 W/m2
with relative error 6% ,and
qmax(Makkah)=878W/m2
,while the experimental value is 938 W/m2
with relative error 6%.
These two values of qmax that are computed theoretically are inserted in equation (5) to fit the corresponding
published experimental data for Jeddahand Makkah. The obtained results are illustrated in figures (1) and(2)
respectively. We considered the published experimental values of qmax for Hong Kong ,Barcelona and Egypt in
fitting process to compute the corresponding q(t) for such cities, since the day number is not available for the
authors. The computed values of q(t) are compared with the corresponding experimental values. The obtained
results for Hong Kong, Barcelona and Egypt are illustrated graphically in Figures(3-11).
As a measure to the degree of fitting , the percentage relative error defined as
cal
cal
q
qq exp
is computed for
each point for the corresponding curve .
IV. RESULTS AND DISCUSSION
The comparison between the published experimental data and the corresponding calculated values
obtained using eq.(5) are illustrated graphically in Figures(1-11).The results show that, bad fitting is obtained
usually at the extreme points of the distribution near sunrise and sunset hours where the level of solar isolation is
not effective for practical exploitation of solar energy .
This situation is attributed to the fact, that the model is oriented to clear sky. While several parameters
affect such distribution function such as sunshine hours ,the declination angle, the latitude ,the altitude ,the
relative humidity (Quraishee, 1969; lin and Jordan, 1960; Sharma and Pal , 1965).It is also a function of such
variables as the nature and the extent of cloud cover , the aerosol and water vapor content of the atmosphere
(Munroe, 1980) .As a result some experimental points give large relative error such as the point indicated in
figure (4) with relative error 81%.Other than the extreme points the obtained relative errors are within 16% .This
indicates that the suggested formula to calculate the daily hourly solar irradiance is acceptable.
V. CONCLUSIONS
The comparison between the predicted global solar radiation q (t) W/m2
incident on a horizontal
surface and the corresponding experimental data shows in general satisfactorily agreement for points other than
the extreme points(~16%). The advantages of the introduced model is that q(t) is given in terms of well
established parameters such as the length of the solar day td and t0=td/2 for symmetrical distribution. Moreover
qmaxcan be estimated in terms of the solar constant adjusted for the variation of the distance between the sun and
the earth along the time of the year .This means that the trial represents a closed theoretical system to evaluate
the incident daily global solar radiation.
REFERENCE
[1]. El-Adawi,M.K.,El-Niklawi, M.M., Kutub, A.A., Al-Barakati,G.G.,1986.
[2]. Estimation of the hourly solar irradiance on a horizontal surface. Solar Energy. 36, 129-134.
[3]. Shalaby,S.A. ,1994. Estimation of diurnal global solar irradiance on a horizontal surface – new
approach.Indian J.Phys.68,503- 513.
[4]. El-Adawi,M.K., AL-Naim,I.A.,2001. Prediction of Diurnal Global solar Irradiance on A Horizontal
Surface-General Approach.proceedings, the first Saudi science conference ,KFUPM,16,63-78.
[5]. El-Adawi ,M.K. ,2002.New approach to modeling a flat plate collector :the Fourier transform technique
.Renewable energy.26,489-506.
[6]. Quraishee,M.M.,1969. Global solar radiation measurements at Kabul , Afghanistan. Solar Energy.12
,387-390 .
[7]. Munroe,M.M. ,1980.Estimation of totals of irradiance on a horizontal surface from U.K. average
meteorological data.Solar Energy.24, 235-238.

[8]. Wieder,S.,1982.AnIntroduction to Solar Energy for Scientist and Engineers. JohnWiley,New York
,p.58.
[9]. Kondratyev,K.Y.,1969. Radiation inatmosphere, Academic New York,p.203.
[10]. El- Adawi,M.K. ,AL-Naim,I.A.,2010. The temperature variation of a solar cell in relation to its
performance.Environmental Science and Engineering.4,56-59.
[11]. Duffiee,J.A. ,Beckman, W.A. ,1974 .Solar Energy Thermal Process, New York, Wiley –interscience.
[12]. Tiwari,G.N.,1997.Solar Thermal Engineering Systems.Narosa publishing House , London .
[13]. Rai,G.D.,1989.Solar Energy Utilization.Khannapublisher,Delhi, p35.
[14]. Hrens,A.,2006.An Introduction to Weather.Climate and the Environment. USA.
[15]. Leung,C.T. ,1980.The fluctuation of solar irradiance in Hong Kong .Solar Energy. 25, 485 -494.
[16]. Villarrubia ,M ., Coronas ,A.,Llorens,M.,1980. Solar radiation incident on tilted flat surface in
Barcelona Spain. Solar Energy. 25,259-263.
[17]. Data on The Hourly Daily Global Solar Irradiance on a horizontal Surface , Cairo 1980
(General Organization for Housing , Building and planning Research Center , Dokki , Cairo ,
Egypt).
[18]. El-Bar,O.,1983.Energy balance over slant surface. M.SC. Thesis , Faculty of Meteorology and
Environment Science, King Abdul Aziz University , Jeddah, Saudi Arabia.
[19]. EL-Gendi,M.H. ,1983. Theoretical and experimental study on Solar flat plate collector report, Physics
Department ,Faculty of Applied Science and engineering ,UMM Al QURA university ,Saudi Arabia.
[20]. Katsoulis, B.D., Parachistopoulos, C.E., 1978.Analysis of Solar radiation measurements at Athens
observatory and estimates of solar radiation in Greece.Solar Energy.21,217-226.
[21]. Khogali, A. , Ramadan, M.RI., 1982 .Optimum operating conditions of a solar cell panel and prediction
of solar radiation in Sanaa ,Yemen.Solar cells.5,173-181.
[22]. Harty, L.E. , Martinez-Lozano, J.A. ,Utrillas , M.P. ,Tena, F. , Pedros, R.,1999.The optimization of
the angle of inclination of a solar collector to maximize the incident solar radiation .Renewable
Energy.17,291-309.
[23]. Cooper , P.I. , 1969.The absorption of solar radiation in solar stills. Solar energy. 12,333-346.
[24]. Lin, B., Jordan, R.C . , 1960. The interrelationships and characteristic of direct ,diffuse ,and total
radiation .Solar Energy. 4,,1-19.
[25]. Sharma, M.R., Pal , R.S., 1965.Total , direct and diffuse Solar radiation in the tropics .Solar Energy.
9,183-192.
Global solar radiation intensity (W/m2
) received on a horizontal surface in Jeddah on (9/4/1982) located at 40̊
25` E, 41̊ 23 `N (El-Bar, 1983)
Fig.1.Experimental and calculated data for Jeddah April
(The maximum relative error equals to 12%)
0
200
400
600
800
1000
0 2 4 6 8 10 12 14
GlobelSolarIrradiance(W/m2)
Shifted time(=t-tr),h
Jeddah April
qexp qcal

) received on a horizontal surface
In Makkah on (9/4/1982) April located at 38.5 °E, 21.5°N (EL-Gendi, 1983)
Fig.2. Experimental and calculated data for Makkah, April
) received on a horizontal surface in Hong Kong on December (1979)
located at 114̊ 10`E, 22̊ 19` N (Leung, 1980)
Fig.3. Experimental and calculated data for Hong Kong, December
0
200
400
600
800
1000
0 2 4 6 8 10 12 14
time,h
Makkah
q exp qcal
0
200
400
600
800
1000
0 2 4 6 8 10 12
Hong Kong December
qexp qcal

) received on a horizontal surface in Hong Kong on November (1979)
Fig. 4.Experimental and calculated data for Hong Kong, November
) received on a horizontal surface in Hong Kong on January (1979)
Fig.5.Experimental and calculated data for Hong Kong, January
0
200
400
600
800
1000
0 2 4 6 8 10 12
Hong Kong november
qexp qcal
0
200
400
600
800
1000
0 2 4 6 8 10 12 14
Hong Kong January
qexp qcal

) received on a horizontal surface in Hong Kong on April (1979) located
at114° 10 `E, 22° 19`N (Leung, 1980)
Fig.6.Experimental and calculated data for Hong Kong, April
) received on a horizontal surface in Barcelona on December (1973-1975)
located at 2° 7` E, 41̊ 23` N (Villarrubia et al., 1980)
Fig. 7.Experimental and calculated data for Barcelona, December
0
200
400
600
800
1000
0 2 4 6 8 10 12 14 16
Hong Kong April
qexp qcal
0
200
400
600
800
1000
0 2 4 6 8 10 12
Barcelona Decmber
qexp qcal

) received on a horizontal surface in Egypt
(Cairo) on July (1980) located at 23̊ 58` N (Cairo, 1980)
Fig.8.Experimental and calculated data for Egypt, July
) received on a horizontal surface in Egypt(Cairo) onMarch (1980) located
at 23̊ 58` N (Cairo- Al-Ahram, 1980)
Fig. 9. Experimental and calculated data for Egypt, March
0
200
400
600
800
1000
1200
0 2 4 6 8 10 12 14 16
GlobalSolarIrradiance(W/m2)
Egypt july
qexp qcal
0
200
400
600
800
1000
0 2 4 6 8 10 12 14
Egypt March
qcal qexp

) received on a horizontal surface in Egypt (Cairo) on June (1980) located
at 23̊ 58` N (Cairo- Al-Ahram, 1980)
Fig. 10. Experimental and calculated data for Egypt, June
) received on a horizontal surface in Egypt (Cairo) on September (1980)
located at 23̊ 58` N (Cairo- Al-Ahram, 1980)
Fig. 11. Experimental and calculated data for Egypt, September
0
200
400
600
800
1000
1200
0 2 4 6 8 10 12 14
Egypt june
qexp qcal
0
200
400
600
800
1000
0 2 4 6 8 10 12 14
Egypt September
qcal qexp

Prediction of the daily global solar irradiance received on a horizontal surface -New approach

More Related Content

What's hot (20)

Similar to Prediction of the daily global solar irradiance received on a horizontal surface -New approach (20)

More from irjes (20)

Recently uploaded (20)

Prediction of the daily global solar irradiance received on a horizontal surface -New approach