C:\Users\raoh\Desktop\Paper_3.xps The Journal of Engineering Research Vol 8 No 2 (2011) 19-31 1. Introduction Building transient thermal simulation software pro- grams such as TRNSYS require hourly weather parameters in standard Typical Meteorological Year format (TMY2, for example). Each line of the TMY2 record has a total of 70 parameters and indices that need to be specified based on measured data or calcu- lated from measured data or estimated. Some of the important elements that normally are not measured are _________________________________________ *Correseponding author’s e-mail: alrawahi@squ.edu.om the hourly terrestrial solar radiation fluxes on horizon- tal surfaces. Terrestrial solar radiation; direct beam, diffuse and global hourly solar fluxes are modelled and calculated based on average daily measured data. Other data such as dew point and wet bulb tempera- tures are also calculated based on the measured dry bulb and relative humidity data. The design of build- ing conventional HVAC systems requires cooling load calculations which are based on hourly weather data. Also, the design of solar cooling systems requires Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman N.Z. Al-Rawahi*, Y.H. Zurigat and N.A. Al-Azri Department of Mechanical and Industrial Engineering, College of Engineering, Sultan Qaboos University, P.O. Box 33, Postal Code 123, Al Khoud, Muscat, Oman Received 7 April 2011; accepted 20 September 2011 Abstract: In this paper, hourly terrestrial radiation: direct beam, diffuse and global solar radiation are modelled and calculated based on daily measured data for a horizontal surface. In addition, the same parameters were modelled for inclined surfaces. Most of the parameters modelled in this work represent a part of the input data required by building thermal simulation and solar energy systems software. Important trends of the solar radiation on tilted surfaces as a function of time and direction are being pre- sented and discussed. The comparison of some of the results with measured data from other sources shows good agreement. The effect of tilt angle and orientation on the incident solar radiation fluxes are presented along with optimum surface tilt angles and directions for maximum solar radiation collection in Muscat area. The results presented in this paper are quite useful for quick estimation of solar radiation for calculations of cooling load and solar collector performance. Also, the models and the computer code developed in this work form the backbone of any computer-aided building thermal design and solar sys- tems design calculations. Keywords: Solar radiation, Solar cooling, Solar collectors 20 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri cooling load calculations, collector and thermal ener- gy storage sizing, and performance evaluation under transient conditions. Thus, the need for accurate pre- dictions of hourly solar flux data for a given location is quite obvious. Solar radiation models are used to predict direct, diffused, and global radiation based on certain known variables. The models can be classified into either parametric or decomposition models. The models in the former category require the knowledge of atmos- pheric conditions (ASHRAE handbook, 1999; Iqbal 1978; Gueymard 1993). The decomposition models, on the other hand, use the known solar global radiation data (Liu and Jordan 1960; Lam and Li 1996; Badescu 2002). A good comparison between the two types of models can be found in (Wong and Chow 2008) where comparison between representative models from both categories with measured data was presented. Decomposition models use solar global radiation data on horizontal surface to predict the direct and diffuse radiation components. The use of these models is sup- ported by the availability of a large number of loca- tions that keep records of the global radiation on a hor- izontal surface. The decomposition models can be fur- ther classified according to the used method as empir- ical, analytical, numerical and statistical, and neural networks models (Liu and Jordan 1962; Lam and Li 1996; Tymvios et al. 2005; Zervas et al. 2008; Ulgen and Hepbasli 2009; Mellit et al. 2010). A number of models were shown to predict the global solar radiation. However, modelling the diffuse solar radiation has been a challenge for researchers. This difficulty in modelling the diffuse radiation com- ponent reflects on the ability to have accurate models for the solar radiation on tilted surfaces which is gen- erally calculated using total, direct, and diffuse solar radiation fluxes. The accuracy of models for solar radi- ation on tilted surfaces is highly affected by the accu- racy of the diffuse solar radiation models. The isotrop- ic model of Liu and Jordan (1962) and the anisotropic model of Perez et al. (1990) are among the well known models in this field and will be used in this work. In their study, (Padovan and Del Col 2010) com- pared the predictions of a number of models used to predict diffuse and total radiation on tilted surfaces including the models of (Liu and Jordan 1962; Perez et al. 1990). They found that for southern tilted surfaces, the models have similar accuracy but for the western tilted surfaces the anisotropic models, like that of Perez et al. (1990), perform better. This is because, in the northern hemisphere unlike the west-tilted sur- faces, the south-tilted surfaces receive more direct solar radiation and less diffuse radiation. They found that the error in the diffuse solar radiation increases the error in estimating radiation on tilted surfaces signifi- cantly. Similar conclusions were reported by Noorian et al. (2008). In the presented work, decomposition models for the prediction of hourly solar fluxes on horizontal and inclined surfaces have been implemented with detailed description presented in the next sections. These mod- els and the computer code developed in this work form the backbone of any computer-aided building thermal design and solar systems design calculations. 2. Estimation of Hourly Solar Flux Compo -nents The total hourly solar radiation falling on a horizon- tal surface, Ith in W.h/m2, could be estimated from the daily global radiation, H (in kJ/m2 per day) in accor- dance with (Duffie and Beckman 1984) as follows: (1) Here 3.6 is a conversion factor (3.6 W.h/m2 = kJ/m2), r is the ratio Ith / H and is given by (Collares- Pereira and Rable 1979): (2) (3) (4) W is the hour angle of the sun (in degrees) given by: (5) h is the time of the day (solar time) in hours and Ws is the sunset hour angle (in degrees), given by: (6) is the location latitude angle (in degrees) and is the declination angle of the sun (in degrees) given by: (7) day is the day of the year under consideration (1 to 365). It should be noted that while the equations involv- ing solar radiation terms use solar time, the different climate data obtained from meteorological stations are normally given in standard time (local time in the 21 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman weather station). In addition, the output from the solar radiation calculations needs to be expressed in stan- dard time as it is used in solar thermal systems simu- lations. It is noted, however, that many published work do not take this time difference into considera- tion. This could be due to the fact that the difference between the standard time and the solar time is some- times small. In this work, the difference between solar and standard timings was taken into account. The well established relation between the solar time tsol and standard time tstd is given by (Duffie and Beckman, 1984): (8) where Lst is the standard meridian for the local time zone, Lloc is the longitude of the weather station, and the parameter T is the equation of time (in minutes) defined as: (9) where B is given by: (10) The difference between solar time and standard time over the year for Seeb/Muscat area is shown in Fig. 1. Although the difference is small reaching a maximum of 24 minutes, the data reported by Seeb Meteorological Center based on local standard time The total solar radiation consists of the direct or beam radiation coming directly from the solar disc and the diffuse component scattered to the ground from the sky dome. The latter depends on the clarity of the sky and could be estimated from the correlation of (Collares-Pereira and Rabl 1979) which gives the daily average diffuse radiation, Hd, as: (11) Where KT is the clearness index for the day, defined as the ratio of the daily radiation on a horizontal sur- face to the daily extraterrestrial radiation (Ho) on that surface, that is: (12) and (13) Isc is the solar constant = 1367 W/m2. The hourly values of the diffuse solar radiation can be estimated from the equation by (Liu and Jordan 1960) which gives: (14) The direct (beam) component of the solar radiation on a horizontal surface is then obtained from: (15) The hourly value of the direct solar radiation on a Ibn (16) Where s is the solar zenith angle calculated by (17) Alawi (1995) for developing Typical Solar Radiation Data for Oman. 3. Estimation of Extraterrestrial Radiation Components The extraterrestrial radiation on a plane normal to the sun rays is a function of the distance between the Day 100 150 200 250 300 350 40050 Figure 1. Difference between solar time and stand ard time over the year for Seeb/Muscat area 22 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri earth and the sun and is given by Duffie and Beckman (1984): (18) Using the solar zenith angle, the extraterrestrial radiation on a horizontal plane can be calculated from: (19) Combining Eqs. (17) and (19), we get the extrater- restrial radiation on a horizontal plane at any hour angle, W as: (20) If the total extraterrestrial solar radiation energy incident on a horizontal plane, Ieh, over a certain peri- od of time is required, then Eq. (20) is to be integrat- ed from hour angles W1 to W2. This results in the fol- lowing relation: (21) Here, Ieh is in J/m2. To have Ieh in W.h/m2 Eq. (21) is to be divided by 3600. On the other hand, the total extraterrestrial radiation on a plane normal to the sun rays over a certain period of time within a certain day is given simply by multiplying Eq. (18) by that time period in seconds. This is because the normal total extraterrestrial radiation is constant over the day. 4. Estimation of Hourly Solar Flux on Inclined Surfaces The total radiation received by an inclined surface consists of beam, diffuse, and reflected radiation. The reflected radiation here is the radiation reflected from the surrounding ground. Different models for estimat- ing the total solar radiation on inclined surfaces had been proposed and evaluated. Such studies include the work of Noorian et al. (2008) who evaluated twelve models, the study by Nijmeh and Mamlook (2000), and that by Notton et al. (1996). Most of the present- ed models in these reviews proved to give good results. The majority of these models use similar terms for the beam and reflected radiation while differ in the method of calculating the diffuse radiation portion. The general form of total hourly solar radiation energy incident on an inclined surface Iti (in kJ/m2), is given by: (22) where ITd is the diffuse radiation on the tilted surface, is the angle the surface makes with the horizontal, r is the surrounding diffuse reflectance for the total solar radiation, and Rb is the hourly geometric factor given by: (23) Here is the incidence angle, the angle between the beam radiation on a surface and the normal to that surface, and it is given by: (24) where is the surface azimuth angle defined as the deviation of the projection on a horizontal plane of the normal to the surface from the local meridian, with zero due south, east negative, and west positive. Among the methods available for prediction of radi- ation fluxes on tilted surfaces is the method of (Liu and Jordan 1962). This is a simple isotropic model that had been tested by many investigators and found to give good results. In this method the diffuse radiation on the tilted surface is given by: (25) Another method for the estimation of the diffuse solar radiation is the anisotropic model of (Perez et al. 1990). Most of the reviewers agree that this model gives better results (Noorian et al. 2008; Padovan and Del Col 2010). The diffuse radiation on the tilted sur- face in the Perez et al. (1990) method is given by: (26) where the coefficients F1 and F2 are given by: (27) where is the sky brightness and F11, F12, F13, F21, F22, and F23 are given in Table 1 as functions of the SCeh 23 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman sky clearness index, , and the terms a and b are given by: (28) (29) The sky's brightness, gives an indication about the thickness of the clouds and is given by: (30) where m is the relative optical air mass (Kasten 1966). The sky clearness index, , represents the transition from totally overcast sky to a low turbidity clear sky. It is given by: (31) It is used to specify the sky clearness categories which are used in Table 1. The upper and lower bound for each category can be found in Perez et al. (1990). For example, a sky clearness index of = 1 represents an overcast sky, whereas a sky clearness index of = 8 represents a clear sky. 5. Comparison with Measured Data To validate the model presented in this work, a com- parison with hourly measured data from different sources has been conducted. Figures 2 and 3 show comparisons of predictions of the model developed in this work with data from Phoenix, Arizona for June 21st and with those from Seeb Meteorological Center, Table 1. Coefficients used in Eq. (27) for different sky clearness, , categories Figure 2. Comparison between predicted and me- asured hourly global radiation data of Figure 3. Comparison between predicted and mea- sured hourly global solar radiation data of Seeb/Muscat/Oman 1200 1000 800 600 400 200 0 5 7 9 11 13 15 17 19 21 Hours 1000 900 800 700 800 600 500 400 300 200 100 0 24 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri respectively. Clearly, the predictions are in good agreement with the measured data. The maximum dif- ference between predicted and measured data does not exceed 7%. This shows that the model is a reliable tool in predicting hourly solar flux data from measured daily averaged data. 6. Results Based on the above formulations a computer pro- gram has been developed to calculate the hourly solar radiation parameters on horizontal and inclined sur- faces. The daily average solar radiation flux data of Typical Meteorological Year for Seeb/Muscat (Zurigat et al. 2003) were used as input. Sample results that show the effect of tilt angle and orientation on the inci- dent solar radiation are presented. In addition, the opti- mum surface tilt angles and directions for maximum solar radiation collection in Muscat area is deduced. This information is useful when a simple solar collec- tor is to be installed without a tracking system. Figure 4 shows the hourly extraterrestrial, global, and diffuse radiation on a horizontal surface on January 21st and June 21st. The extraterrestrial and the global radiation increase from January to June by about 30% and 80%, respectively. The increase in the Figure 4. Terrestrial and extraterrestrial solar radiation in Seeb/Muscat area during January 21st and June 21st 25 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman Figure 5. Sun path diagram for Seeb / Muscat area Figure 6. Solar radiation in Seeb/Muscat area for planes with different tilt angles, , towards South (a), East (b), North (c), and West (d) directions for Januatry 21st 26 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri diffuse solar radiation is about 95%. These significant increases in the global and diffuse solar radiation indi- cate the matching between the availability of solar energy and the need for cooling during hot summer months. Thus, the potential use of solar thermal cool- ing where it is needed most is obvious. A number of simulations were carried out to inves- tigate the effect of the tilt angle and the tilt direction. To make sense of the results, the sun path diagram for Seeb/Muscat area is presented in Fig. 5. The diagram shows that the sun is in the southern half of the sky most of the year, especially during the midday period when the solar radiation is at its daily peak. However, the sun is in the northern hemisphere during most of June which represents the hottest month of the year. It can be concluded from the diagram that, for most part of the year, tilting solar collection surfaces towards the south will increase the incident solar radiation on these surfaces. This is clearly illustrated by the results pre- sented in this study. Figures 6 and 7 show the hourly solar radiation in Seeb/Muscat area for a flat surface at different tilt angles towards South, East, North, and West directions during January 21st and June 21st, respectively. With the exception of surfaces tilted towards the south, the solar radiation on horizontal and tilted surface during June 21st (see Fig. 7) is greater than that during January 21st. The general trend of the incident solar radiation for surfaces tilted towards the east and the west is almost identical for the two days with higher values for June 21st. On the other hand, as the tilt angle towards the North increases, a higher effect on the hourly solar radiation during January 21st is noticed. It is also noticed that the general trend of the incident solar radiation is almost the same with north-tilted sur- faces during January 21st, i.e. the solar radiation increases until it reaches a maximum at noon time then it decreases. While tilting surfaces towards the north has less effect on the solar radiation magnitude, the general trend of the solar changes as the tilt angle increases beyond 60o. For larger tilt angles, the increase of the solar radiation levels off before noon and starts to decrease to have a local minimum point at noon. This is clearly understood by looking at the solar path diagram (see Fig. 5). In June, before 12 noon, the sun moves from the northern half to the southern half of the sky which means that the north tilted surfaces will receive less direct solar radiation. The solar radia- tion on the inclined surface during this period consists of the diffuse solar radiation and the reflected solar radiation only. This shows the importance of having good models for these two components. Figure 7. Solar radiation in Seeb/Muscat area for planes with different tilting angles, , towards South (a), East (b), North (c), and West (d) directions for June 21st 27 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman For south-tilted surfaces, Figs. 6a and 7a show dif- ferent trends. Compared with June 21st, the solar radi- ation during January 21st is less sensitive to changes in the tilt angle. The maximum radiation is reached at a tilt angle of 40o at which the received solar radiation increases 40% above that received by a horizontal sur- face. In June 21st, the received solar radiation decreas- es significantly as the tilt angle increases. The trends shown in Figs. 6 and 7 indicate that strategies to con- trol the amount of received solar radiation on a surface depend on the surface orientation. The combined effect of tilt angle and surface orien- tation on the total incident solar radiation during the months of January and June is shown in Fig. 8. Orientation angles, , of 0o, 90o, 180o, and 270o corre- spond to South, West, North, and East directions, respectively. This figure also shows whether orienta- tions other than the main four (E, W, S, N) have desir- able features in terms of received solar radiation. While the figure shows different trends for the two months, there are no special features for other orienta- tions regarding the received solar radiation. During January the minimum received solar radiation is for surfaces tilted towards the North while the maximum received solar radiation is for surfaces tilted towards the South. During the month of June this behaviour is reversed. In January, the maximum received solar radiation is achieved with a tilt angle of about 40o towards the south. During June, the horizontal surface collects more solar radiation than other tilted surfaces Figure 8. Solar radiation in Seeb/Muscat area for surfaces of different inclination angles, , as a function of the oriengation angle, , for January (a) and June (b). Oriengation angle is measured with 0o due South and positive towards East 28 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri Figure 9. Average monthly incident solar radiation in Seeb/Muscat area for planes with different tilt angles angles, , toards South (a), East (b), North (c), and West (d) directions, respectively Figure 10. Average monthly solar radiation in Seeb/Muscat area for vertical planes with different orientations. Orientation angle is measured with 0o due South and positive towards East 29 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman except for surfaces with a tilt angle of 20o towards the north where the received solar radiation is slightly higher. It should be noted that the results presented here assume that the surrounding diffuse reflectance is isotropic, i.e. the same in all directions (Dry bare ground of r = 0.2). Looking at the broader picture, Fig. 9 shows the incident solar radiation on tilted surfaces for the whole year. As noticed previously, surfaces tilted towards the East and the West have similar trends and the trend is almost the same for other tilt angles. For surfaces tilt- ed towards the South, in the period between October to March, the effect of tilt angle on the received solar radiation is relatively small. However, the tilt angle has more significant effect between April and September. For surfaces tilted towards the North the received solar radiation is sensitive to the tilt angle throughout the year. Vertical surfaces ( = 90o) require further attention as they represent the major part of the building envel- op which contributes significantly to the cooling/heat- ing loads. Figure 10 shows the received solar radiation for vertical surfaces at different compass orientations. During summer the West- and East-facing vertical sur- faces receive higher solar radiation than the South- and North-facing ones. These results explain why the ori- entation of buildings in Muscat area should be East- West with the longest vertical surface facing south. During summer time, this orientation minimizes the solar gain and consequently, the cooling load. It is noticed also that South-facing vertical surfaces receive more solar radiation than North-facing ones for most part of the year and the difference is considerably large. Tables 2a and 2b present quite useful information regarding solar radiation in Seeb/Muscat area. For each month it shows the average solar radiation on horizontal surfaces, the tilt angle that results in maxi- mum received solar radiation and the corresponding solar radiation, the range of tilt angles that will main- tain the received solar radiation within 1% from the maximum value, and the solar radiation on vertical surfaces in the four main directions. The information about the received solar radiation on horizontal and vertical surfaces can be quite useful in estimating the building heat gain. Also, the information about the maximum received solar radiation and the tilt angle are useful for sun tracking solar collectors that are manually adjusted every month or so. One of the gen- eral rules of thumb in solar radiation is that during Table 2a. Summary of incident solar energy for tilted surfaces for the months of January through June using the model of (Perez et al. 1990) Table 2b. Summary of incident solar energy for tilted surfaces for the months July through December using the model of (Perez et al. 1990) 30 N.Z. Al-Rawahi, Y.H. Zurigat and N.A. Al-Azri summer the tilt angle that will give the maximum received solar radiation equals the location latitude minus 15o. The presented data is in agreement with this rule since the latitude of Seeb/Muscat area is 23.35o. Another rule of thumb is that during winter time the tilt angle that will give the maximum received solar radiation equals the location latitude plus 15o. This again is in agreement with the presented range of tilt angle for maximum received solar radiation. Figure 11 shows the average received solar radia- tion if the tilt angle is kept the same throughout the year. When the tilt angle of a solar collector is fixed, a tilt angle of about 25o will receive the maximum solar radiation. Figure 11 also shows the resulting reduction in the incident solar radiation, in percentage, as a function of tilt angle. For example, using a fixed horizontal surface ( = 90o) will result in about 10% reduction in total solar radiation incident on that sur- face. Conclusions In this paper, prediction of hourly terrestrial solar radiation on a horizontal surface: direct beam, diffuse and global from measured daily averaged global solar radiation on the same surface has been conducted. The same parameters were also modelled for inclined sur- faces. Thus, an important tool for building and solar thermal systems design, simulations, and performance assessments has been developed. The predicted hourly solar radiation incident on a horizontal surface was compared with hourly data measured locally at Seeb Meteorological Centre as well as elsewhere outside Oman. The effect of tilt angle and orientation on the inci- dent solar radiation fluxes are presented along with optimum surface tilt angles and directions for maxi- mum solar radiation collection in Muscat area. This information is useful when a simple solar collector is to be installed without a tracking system. In general, surfaces tilted towards the East and the West have sim- ilar trends and the trend is almost the same for other tilt angles. For surfaces tilted towards the South, in the period between October to March, the effect of tilt angle on the received solar radiation is relatively small. However, the tilt angle has more significant effect between April and September. For surfaces tilt- ed towards the North the received solar radiation is sensitive to the tilt angle throughout the year. During summer the West- and East-facing vertical surfaces receive higher solar radiation than the South- and North-facing ones. These results explain why the orientation of buildings in Muscat area should be East- West with the longest vertical surface facing south. In addition, a table that presents important monthly solar radiation information for horizontal, inclined and ver- tical surfaces has been presented. Assuming isotropic reflection of diffuse radiation, the results show that for January the maximum received solar radiation is Figure 11. Average daily incident solar radiation energy in Seeb / Muscat area for different tilt angles 31 Prediction of Hourly Solar Radiation on Horizontal and Inclined Surfaces for Muscat/Oman achieved with tilt angle of about 40o towards the south. For June, the horizontal surface collects more solar radiation than other tilted surfaces except for surfaces with a tilt angle of 20o towards the north where the received solar radiation is slightly higher. It is believed that the results presented in Tables 2a and 2b are quite useful for quick estimation of solar radiation for calcu- lations of cooling load and solar collector perform- ance. Acknowledgments The Research leading to these results has received Research Project Funding from the Research Council of the Sultanate of Oman, Research Agreement No. CR SQU 010 001. References Al-Hinai AH, Al-Alawi SM (1995), Typical solar radiation data for Oman. Applied Energy 52:153- 163 ASHRAE handbook (1999), HVAC applications, Atlanta (GA): ASHRAE. Collares-Pereira M, Rabl A (1979), The average dis- tribution of solar radiation - correlations between diffuse and hemispherical and between daily and hourly insolation values. Solar Energy 22:155- 164. Duffie JA, Beckman WA (1980), Solar engineering of thermal processes. John Wiley, New York. Gueymard C (1993), Critical analysis and perform- ance assessment of clear-sky solar irradiance models using Theoretical and Measured Data. Solar Energy 51(2):121-138. Iqbal M (1978), Estimation of the monthly average of the diffuse component of total insolation on a hor- izontal surface. Solar Energy 20(1):101-105. Kasten A (1966), A new table and approximate formu- la for relative optical air mass. Arch. Meteorol Geophys. Bioklimatol Ser. B 14:206-223. Lam JC, Li DHW (1996), Correlation between global solar radiation and its direct and diffuse compo- nents. Building and Environment 31(6):527-535. Liu BYH, Jordan RC (1962), Daily insolation on sur- faces tilted towards the equator. Trans ASHRAE 67:526-541. Mellit A, Eleuch H, Benghanem M, Elaoun C, Massi Pavan A (2010), An adaptive model for predicting of global, direct and diffuse hourly solar irradi- ance. Energy Conversion and Management 51:771-782. Nijmeh S, Mamlook R (2000), Testing of two models for computing global solar radiation on tilted sur- faces. Renewable Energy 20(1):75-81. Noorian AM, Moradi M, Kamali GA (2008), Evaluation of 12 models to estimate hourly diffuse irradiation on inclined surfaces. Renewable ener- gy 33(6):1406-1412. Notton G, Muselli M, Louche A (1996), Two estima- tion methods for monthly mean hourly total irradi- ation on tilted surfaces from monthly mean daily horizontal irradiation from solar radiation data of ajaccio, corsica. Solar Energy 57(2):141-153. Padovan A, Del Col D (2010), Measurement and modeling of solar irradiance components on hori- zontal and tilted planes. Solar Energy 84:2068- 2084. Perez R, Ineichen P, Seals R, Michalsky J, Stewart R (1990), Modelling daylight availability and irradi- ance components from direct and global irradi- ance. Solar Energy 44:271-289. Tymvios FS, Jacovides CP, Michaelides SC, Scouteli C (2005), Comparative study of angstrom's and artificial neural networks methodologies in esti- mating global solar radiation. Solar Energy 78:752-762. Ulgen K, Hepbasli A (2009), Diffuse solar radiation estimation models for turkey's big cities. Energy Conversion and Management 50:149-156. Zervas PL, Sarimveis H, Palyvos JA, Markatos NCG (2008), Prediction of daily global solar irradiance on horizontal surfaces based on neural-network techniques. Renewable Energy 33:1796-1803. Zurigat YH, Sawaqed N, Al-Hinai H, Jubran B (2003), Development of typical meteorological years for different climatic regions in Oman, Final Report, Petroleum Development Oman.