Prediction of Symmetrical and Asymmetrical of Diurnal Global Solar Irradiance Distribution—New Approach ()
1. Introduction
The prediction of the diurnal global solar radiation q(t) W/m2 is needed as one important input parameter to study theoretically the design and performance estimation of solar systems for solar energy exploitation, for example, the performance and efficiency of a solar cell, flat plate collector, water heating and treatment, pool heating, space heating, solar cookers, (Heating, Ventilation and air conditioning) (HVAC) technological systems. It is also required to study the production of electricity using molten salt technologies in which the liquid salt is pumped through panels in a solar collector for further heating.
Analysis of solar radiation measurements has aroused the interest of many investigators. As an example, different distributions for solar irradiance with different fitting degrees are given [1] - [6]. A lot of experimental meteorological data for many locations are published [7] - [14]. Trials to introduce governing formulae are given [5] [10] [11] [15].
The received solar energy is a function of several variables [2] such as the nature and extent of cloud cover, the aerosol and other atmospheric constituents such as O2, N2, CO2, O3, dust etc.
Such a function depends also on other parameters such as the sunshine hours, the solar declination angle, the latitude, the altitude and the relative humidity [2]. As a result of these challenges, it is not always possible to predict theoretically the actual shape of such a function to get accurate values of the received irradiance for a given location. Different trials are given by different authors with different degrees of fitting accuracy [1] [2] [3] [7] [16] [17] [18]. Most of such trials are either semi empirical or incomplete to form a closed system or they are difficult to be integrated.
The need for more accurate trials with better fitting degrees is still required. El-Adawi et al. [2] introduce a power expression for such a function, the parameters of which were determined through the least fitting technique. The given expression is not easy to be integrated.
Good fitting with published experimental meteorological data is obtained with maximum relative error 11%. Other trials expressed the required distribution in the form of polynomial [in(t − tmax)] [16] with relative maximum error 16%, or in
[18] with maximum relative error 15% or as polynomial in (t/td) with a correction factor
[4] with maximum relative error 15%.
The present trial represents a new approach to introduce a suggested formula based on well-established solar data such as the length of the solar day “td” in hours, which is well defined in [15], and is also expressed through the maximum value of the daily solar irradiance qmax W/m2. The expression for (td) is well defined in literatures on meteorological basis [15].
To get a closed system, the value of qmax is suggested in terms of the extraterrestrial solar constant adjusted for the variation of the distance between the sun and the earth along the time of year [2] [3]. Thus the introduced distribution is not a semi-empirical one. This is an advantage of the present trial. Moreover, it can be easily integrated and thus it is feasible for practical applications.
A comparative study between the experimental meteorological published data of the received global solar irradiance in different locations [7] [9] [12] [13] and that computed using the present suggested model is given. The relative errors are indicated.
2. Theory
The experimental measured meteorological values of the global solar irradiance q(t), W/m2 received on a horizontal surface as measured by different authors [1] [2] [7] revealed a symmetrical distribution about a maximum average value qmax acquired at midday time (tmax between sunrise tr and sunset (ts) i.e.
.
This symmetrical behavior is shown to be true for the whole solar year [7].
Moreover, the behavior of this function for different locations reveals its universal character [1] [2] [8] [9] [17] [19].
However, some authors [4] discussed the case of asymmetrical distribution for which qmax occurs at “tmax” shifted from the midday times
i.e.
This case will be considered in the present trial.
In the present trial the suggested model to predict the function q(t) W/m2 is given in the form:
(1)
Shifted time scale is considered for which tr = 0. This distribution satisfies the following conditions:
i) At
(2)
ii) At
(3)
iii) At
(4)
This gives:
(5)
iv) At
(6)
This gives:
(7)
For symmetrical distribution,
.
This gives:
(8)
Finally, one gets for symmetrical distribution the following expression:
(9)
For asymmetrical distribution:
(10)
For symmetrical distribution the total daily solar energy received per unit area of a horizontal surface is given as:
(11)
Authors of different trials obtained for the same integral the following values:
[1]
[4]
[16]
[17]
While for asymmetrical distribution the obtained value is:
[4].
This shows that the daily totals of the global solar irradiance on a horizontal surface depends on the degree of symmetry about the point t = tmax, at which the received solar irradiance attains its maximum value [4]
It is worth to note that the length of the day “td” can be expressed in terms of the latitude L and the solar declination “δ” [15] as follows:
(12)
where,
(13)
and “n” is the day number of the year starting from 1 January i.e., (1 ≤ n ≤ 365).
To get a closed system of equations, the physical quantity qmax is suggested on theoretical basis to be in the form [2] [3] [5] [10] [11] :
(14)
where,
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 [11] :
(15)
And, s = 1353 W/m2 [10] is the solar constant. It is worth to note that qmax is computed [3] according to Equation (14) for Jeddah and Makkah.
The estimated value qmax (Jeddah) = 856.8 W/m2 while the experimental value is 915 W/m2 and qmax (Makkah) = 878 W/m2, while the experimental value is 938 W/m2. The obtained relative error is 6% [3].
And “α“ is a correction factor. It is well defined in the literatures [2] [10] Its value is given by the relation:
(16)
Its value is estimated [3] to be 0.65 and 0.63 for Jeddah and for Makkah respectively.
This value depends on the optical thickness, the reflectivity of the underlying terrain and depends also on the solar zenith angle [2] [3].
3. Computations
The function q(t) is computed according to Equation (9) for different locations. The obtained results are compared with the corresponding meteorological published data as illustrative examples.
The measure of fitting is taken as
This step is summarized as follows:
1) The considered data for Barcelona (Spain) (41˚23'N, 2˚7'E) January 1973 [9] are given in Table 1 and are illustrated graphically (as shown in Figure 1).
2) The data for Hong Kong (China) (22˚19'N, 114˚10'E) December 1978 [7] are given in Table 2 and are illustrated graphically (as shown in Figure 2).
Table 1. Comparison between experimental [9] and computed values (Equation (9)) for the incident solar irradiance (W/m2) for Barcelona (Spain) (Latitude 41˚23'N, Longitude 2˚7'E), (June 1973).
*tmax = 8 hr, td = 16 hr, qmax = 710 W/m2.
Figure 1. Comparison between experimental [9] and computed values (Equation (9)) for the incident solar irradiance (W/m2) for Barcelona (Spain).
Table 2. Comparison between experimental [7] and computed values (Equation (9)) for solar irradiance (W/m2) Hong Kong. China. December (1978) [22˚19'N, 114˚10'E].
tr = 6.5, ts = 17.5, tmax = 12, td = 11, qmax = 594.5 W/m2.
Figure 2. Comparison between experimental [7] and computed values (Equation (9)) for solar irradiance (W/m2) Hong Kong (China).
3) The data for Jeddah (Saudi Arabia)
(21˚37'N 40˚25'E) April 1982 [12] are given in Table 3 and are illustrated graphically (as shown in Figure 3).
4) The data for Makkah (Saudi Arabia) (38.5˚E, 21.5˚N) March 1983 [13] are given in Table 4 and are illustrated graphically (as shown in Figure 4).
The relative errors obtained according to our model are compared with the corresponding published meteorological values obtained for the same locations and at the same local day time as shown in the corresponding tables.
It is revealed that our model is promising and gives better fitting relative to some other trials as [4] [17] irrespective of the extreme points. The model itself gives good fitting for the intermediate points. This is the more effective one. The
Table 3. Comparison between experimental [12] and computed values (Equation (9)) for solar irradiance (W/m2) Jeddah Located at [21˚37'N, 40˚25'E] (April 1982).
[tr = 6.8, ts = 18.66, tmax = 12.37, td = 12.58, qmax = 915 W/m2].
Figure 3. Comparison between experimental [12] and computed values (Equation (9)) for solar irradiance (W/m2) Jeddah (Saudi Arabia).
Table 4. Comparison between experimental [13] and computed values (Equation (9)) for the incident solar irradiance (W/m2) for Makkah (38.5˚E, 21.5˚N), (March 1983).
tmax = 6.5 hr, td = 13 hr, qmax = 938.
Figure 4. Comparison between experimental [13] and computed values (Equation (9)) for the incident solar irradiance (W/m2) for Makkah (Saudi Arabia).
obtained relative errors are: 12% for Hong Kong and 7% for Barcelona, Jeddah and Makkah respectively.
4. Conclusions
1) The introduced trial to predict the daily global solar irradiance for clear days is promising. It gives good fitting
when compared with the corresponding measured data.
2) The introduced formula can easily be integrated along the local day time to get the total energy received per day per unit area.
This is of vital importance for technological applications.
3) The symmetrical and asymmetrical distributions are considered.
4) The given distribution is expressed in terms of a well-defined parameter which is the length of the solar day td [15]. It is also expressed in term of the maximum value of the solar irradiance qmax attained during the considered day.
5) The latter is suggested in the present study to be expressed in terms of a modified solar constant. Thus the formula is based totally on pure theoretical arguments.
Nomenclatures
t, Time variable (hr.).
tr, sunrise time (hr.).
ts, sunset time (hr.).
td = (ts − tr), The length of the day (hr.).
tmax, The midtime between sunrise and sunset.
q(t), W/m2 Solar irradiance.
L, latitude.
Greek Symbols
, Solar declination (defined in the text).