D:\Desktop\Mech100102(3).xps TJER 2012, Vol. 9, No. 2, 27-35 1. Introduction The knowledge of the availability of solar radiation at any given location is important for the design and performance evaluation of solar drying equipment. The data on hourly, daily and monthly average solar radiation strike on horizontal and inclined surfaces are essential for designing solar drying equipment. Actual measurements of insolation are generally made in only a few places because of the cost of the instruments required, and the care needed for maintaining it. Simulations can be made on the basis of theoretical and statistical relationships developed for estimating ________________________________________ *Corresponding author’s e-mail: basunia@squ.edu.om insolation from various weather observations, most of which are usually made at many locations. The tem- perature and relative humidity have great influence on the open air natural sun drying of farm crops as it is commonly practiced in Bangladesh and in other devel- oping countries. In most tropical countries, including Bangladesh, the climate is characterized by hot and humid air during harvesting. This unfavorable climat- ic condition dictates the need for a more effective method of drying, grain. Because of the difficulties with open air sun drying particularly in wet weather, an alternative and efficient drying method is needed in tropical countries. Simulation of Solar Radiation Incident on Horizontal and Inclined Surfaces MA Basunia*a, H Yoshiob, and T Abec aDepartment of Soils, Water & Agricultural Engineering, College of Agriculture and Marine Sciences, Sultan Qaboos University, P.O. Box 34, Al-Khod 123, Muscat, Sultanate of Oman bDepartment of Biomechanical Systems, Ehime University, Matsuyama, 790-8566, Japan cDepartment of Biomechanical Systems, Ehime University, Matsuyama, 790-8566, Japan Received 2 January 2010; accepted 25 December 2011 Abstract: A computer model was developed to simulate the hourly, daily and monthly average of daily solar radia- tion on horizontal and inclined surfaces. The measured hourly and daily solar radiation was compared with simulat- ed radiation, and favourable agreement was observed for the measured and predicted values on clear days. The meas- ured and simulated monthly averages of total (diffuse and beam) daily solar radiation were compared and a reason- able agreement was observed for a number of stations in Japan. The simulation showed that during the rice harvest- ing season, September to October, there is a daily average of 14.7 MJ/m2 of solar irradiation on a horizontal surface in Matsuyama, Japan. There is a similar amount of solar radiation on a horizontal surface during the major rice har- vesting season, November to December, in Bangladesh. This radiation can be effectively utilized for drying rough rice and other farm crops. Keywords: Solar energy, Simulation, Geometric factor, Daily and hourly solar radiation 28 MS Basunia, H Yoshio and T Abe Bangladesh receives abundant solar radiation throughout the year, so a solar heated dryer can be used to solve the drying problem. Most of the farmers in Bangladesh operate on a small scale and cannot afford mechanically powered drying systems. An intermediate solution takes advantage of the ready availability of solar energy to utilize it for drying. It was reported (Bala and Zia-Uddin 1990) that solar dryers can be designed to enhance the effectiveness of solar drying in Bangladesh, and that there is potential for their adoption and application by small farmers. This would reduce grain loss and help to maintain the quality of dried grain. Solar drying systems must be properly designed in order to meet the particular dry- ing requirements of specific crops and to give a satis- factory performance with respect to energy require- ments (Steinfeld and Segal 1986). The designer should investigates both theoretically and practically, the physical needs of solar crop dryers and evaluate the performance of the solar dryers based on available solar energy. This paper deals with a simulation model to be used for the simulation of hourly, daily, and monthly average radiation incidents, both on horizon- tal and tilted surfaces. 2. Theoretical Considerations 2.1 Radiation Geometry The geometric relationships between a plane of any particular orientation relative to the earth at any time and the incoming beam solar radiation, that is, the position of the sun relative to that plane, can be described in terms of several angles (Benford and Bock 1939). Simplification of the geometric relation- ships relating to the angle of incidence of beam radia- tion on a surface to the other angles leads to the fol- lowing equation. (1) where is the angle of incidence the angle between the beam radiation on a surface and the normal to that surface; z is the angle of incidence; is the angle between the plane of the surface considered and the horizontal, 0 180º, s is the solar azimuth angle; is the surface azimuth angle, the deviation of the projection on a horizontal plane of the normal to the surface from the local meridian, with zero due south, -180º 180º. For the horizontal surface, the angle of incidence is the zenith angle of the sun z. For this condition, = 0, and Eqn. (1) becomes (2) where is the latitude, the angular location north or south of the equator; is the angular position of the sun at solar noon, called declination; is the hour angle, the angular displacement of the sun east or west of the local meridian due to rotation of the earth on its axis at 15º per hour, with morning negative and after- noon positive;. The relationship for the angle of incidence of sur- faces sloped due north or due south can be derived from the fact that surfaces with slope to the north or south have the same angular relationship to beam radi- ation as a horizontal surface at an artificial latitude of ( - ). Modifying Eqn. (2) yields (3) The declination can be computed from the fol- lowing equation (Duffie and Beckman 1991): (4) where n is the day of the year. Recommended average days for months and values of n by month are given by Klein (1977) and are shown in Table 1. Equation (2) can be solved for the sunset hour angle, s, when z = 90º (5) The number of day light hours is given by follow- ing equation: (6) 2.2 Geometric Factor (Rb) The ratio of direct (beam) radiation on the tilted sur- face to that on a horizontal surface is called geometric factor and it is usually expressed as: (7) Cos z and Cos are both determined from Eqns. (2) and (3). The optimum azimuth angle for the flat- plate collectors is usually zero in the Northern Hemisphere and 180º in the Southern Hemisphere. In this case Eqs. (2) and (3) can be used to determine Cos z and Cos , respectively, leading in the northern hemisphere, ( = 0), to (8) 2.3 The Theoretical (Extraterrestrial) Radiation on a Horizontal Surface The amount of extraterrestrial radiation or radia- tion above Earth's atmosphere at any time of year is given by the following equation (Duffie and Beckman 1991). 29 Simulation of Solar Radiation Incident on Horizontal and Inclined Surfaces (9) where En is the extraterrestrial radiation measured on the plane and normalized to the radiation on the nth day of the year, W/m2; and Sc is the solar constant (1367 W/m2). At any given time, the extraterrestrial radiation inci- dent on a horizontal plane outside of the atmosphere is the normal incident of solar radiation as given by Eqn. (9) and divided by Rb. (10) Substituting the value of Cos z from Eqn. (2) in Eqn. (10) yields E for a horizontal surface at any time between sunrise and sunset as: (11) (12) (13) The hourly extraterrestrial radiation strikes on a horizontal surface can be estimated from Eqn. (12) by integrating it between hour angles 1 and 2. The integration for an hour period results in the following equation, where Io is the hourly extraterrestrial radi- tion on horizontal surface, J/m2. (14) The hourly extraterrestrial radiation can also be determined from the following approximate equation, evaluating at the mid point of the hour. (15) The daily theoretical solar radiation, H0 striking on a horizontal surface is obtained by integrating Eqn. (12) over the period from sunrise to sunset. The inte- gration yields where H0 is the daily extraterrestrial radiation, J/m2, and s is the sun set hour angle in degrees. (16) 2.4 Estimation of Average Solar Radiation on a Horizontal Surface The original Angstrom-type regression equation relates monthly average daily radiation to radiation on a clear day, at a given location. The verage fraction of possible sunshine hours is given by (17) Table 1. Recommended average days for months and value of n (the day of the year) by months Source: Klein (1977) m d c m 30 MS Basunia, H Yoshio and T Abe where Hm is the monthly average daily radiation on a horizontal surface Hc is the average clear sky daily radiation for the location and month considered a, and b are empirical constants, nd is the monthly average daily hours of bright sunshine Nm is the monthly aver- age of the maximum possible daily hours of bright sunshine i.e., the day length represented by N of the average day of the month in Eqn. (6). Page (1964) and others have modified the method to base it on extraterrestrial radiation on a horizontal surface rather than on clear day radiation, using equa- tion: (18) where Hom is the extraterrestrial radiation for the given location, averaged over the period considered, and a1 and b1 are constants depending on location. Lof et al. (1966) developed sets of constants, a1 and b1 for vari- ous climate types and locations based on radiation data then available. For Poona, India, a1 = 0.30 and b1 = 0.51 through Eqn. (18). The ratio (Hm/Hom) is termed as the monthly average clearness index. 2.5 Estimation of Hourly Solar Radiation on a Horizontal Surface from Daily Data When hour-by-hour performance estimations for a solar system are needed, it is necessary to start with daily data and then to estimate hourly values from daily data. The ratio of hourly total radiation I to daily total radiation H on a horizontal surface is expressed as a function of day length. The ratio I/H is represent- ed by the following equation of (Collares-Pereira, Rabl 1979): (19a) where the parameters c1 and c2 are given by: (19b) (19c) The ratio of hourly diffuse radiation (Id ) to the daily diffuse radiation (Hd) on a horizontal surface is a func- tion of time and day length. Liu and Jordan (1962) reported that Id/Hd is the same as I0/H0, and is estimat- ed from the following equation. (20) 2.6 Beam and Diffuse Components of Hourly Global Radiation on a Horizontal Surface An hourly clearness index (kh) can be defined as (21) where I is the hourly global radiation on a horizontal surface, J/m2. The usual approach is to correlate Id/I, which expresses the fraction of the hourly radiation on a horizontal plane which is diffuse, with kh, the hourly clearness index. Id is the hourly diffuse component of global radiation (I). This correlation may not represent a particular hours very closely, but over a number of hour it adequately represents the diffuse fraction (Duffie and Beckman, 1991). Orgill and Hollands (1977) used data from Canadian stations (Erbs et al. 1982) have used data from United States and one Australian station, and (Reindl et al. 1990) used independent data from the United States and Europe. The three correlations are essentially identical, although they were derived from three sepa- rate data sources. The correlation between hourly dif- fuse and global radiation developed by (Erbs et al. 1982) is as follows: (22a) (22b) (22c) 2.7 Beam and Diffuse Components of Daily Global Radiation on a Horizontal Surface The daily clearness index, kd is defined as the ratio of a particular day’s radiation H on a horizontal sur- face, to the extraterrestrial radiation H0 for that day. Hd is the daily diffuse component of the total global radi- ation, H on that particular day, J/m2. Mathematically it is expressed as (23) The following statistical relationship was given by Collares-Pereira and Rabl (1979): (24a) 31 Simulation of Solar Radiation Incident on Horizontal and Inclined Surfaces (24b) (24c) (24d) 2.8 Beam and Diffuse Components of Monthly Global Radiation on a Horizontal Surface The monthly average clearness index (Km) is defined as the ratio of monthly average daily radiation on a horizontal surface, Hm to the monthly average daily extraterrestrial radiation H0m. Hdm is the diffuse component of the total global monthly average daily radiation Hm, J/m2. Mathematically it can be expressed as (25) The (Erbs et al. 1982) correlation between the frac- tion of the monthly average daily radiation Hdm/Hm which is diffuse, and the monthly average clearness index Km is represented by the following equations: (26a) (26b) 2.9 Hourly Radiation on a Tilted Surface The radiation on the tilted surface includes three components: beam, isotropic diffuse, and solar radia- tion diffusely reflected from the ground (Reindl et al. 1990). Considering all diffusion is isotropic, the result is an equation that gives the hourly radiation on a tilt- ed surface, It , in terms of parameters that can be deter- mined either theoretically or empirically as is seen in Eqn. 27. (27) where Ib is the beam component of I, Rb is the geomet- ric factor, Fvs is the view factor to the sky for a surface sloped at an angle from the horizontal, Fvg is the view factor of the surface to the ground, and r is the surrounding diffuse reflectance for the total solar radi- ation. Equation (27) can be written as (28) where It is total hourly radiation incident on a tilted surface and is the slope of the surface considered in degrees. 2.10 Daily Total Radiation on a Tilted Surface If the diffuse and ground reflected radiations are each assumed to be isotropic, then in a manner similar in Eqn. (27), for a given day, the radiation on tilted sur- face can be written as (29) where Rbm is the monthly average daily geometric fac- tor, or the ratio between the monthly average daily beam radiation on the tilted surface and that on a hor- izontal surface. For a surfaces with = 0, Rbm is expressed as (30) where m is the sunset hour angle for the sloped sur- face for the mean day of the month, which is given by (31) where "min" means the smaller value of the two items in the brackets. If we divide by the monthly average daily extraterrestrial radiation H0m and substitute K = H/H0m, then Eqn. (30) becomes (32) The ratio Hd/H is the daily fraction of the diffuse radiation and can be found from Eqn. (24) as a func- tion of kd. Therefore Eqn. (32) can be used to calculate the daily total radiation on a tilted surface. for 0.17 kd < 0.75 32 MS Basunia, H Yoshio and T Abe 2.11 Monthly Average Daily Radiation on a Tilted Surface For use in solar process design, the monthly aver- age daily radiation on a tilted surface is needed. The method described by Liu and Jordan (1963) which is extended by Klein (1977) has been widely used. If the diffuse and ground reflected radiations are each assumed to be isotropic, then in a manner similar to Eqn. (28), the monthly mean solar radiation on an unshaded tilted surface can be expressed as (33) where Htm is the monthly average daily total radiation on a sloped surface and Hdm/Hm is a function of Km. 3. Simulation Procedure for Hourly and Daily Radiation on Horizontal and Tilted Surfaces The basic equations to be used are Eqs. (28) and (29). The first steps are to obtain Id/I, Rb, Hd/H and Rbm. The ratios Id/I and Hd/H are functions of kh and kd, and can be obtained from Eqns. (22) and (24), respectively. For a particular day of the year and loca- tion where n and are known, can be calculated from Eqn. (4). For known values of and , sunset hour angle s and day light hours Nm or N are calculat- ed from Eqns. (5) and (6), respectively. For the calcu- lation of hourly radiation on a tilted surface, Rb is needed and can be computed from Eqn. (8). The cal- culation of Rbm requires sunset hour angle of the sloped collector. Equation (31) is used to calculate the sunset hour angle m. With known n and s, from Eqns. (15) and (16), the hourly extraterrestrial (theo- retical) radiation I0 and daily extraterrestrial (theoreti- cal) radiation, Ho, can be estimated. Then with assumed daily hours of bright sunshine (nd), the nd/Nm ratio is calculated to be used in Eqns. (18) and to find the value of daily radiation H on a horizontal surface. From the known value of H, the value of hourly radi- ation I on a horizontal surface is determined from Eqn. (19). So both kh and kd are known from Eqns. (21) and (23), respectively. From known kh and kd ratios, Id/I and Hd/H ratios can be calculated from Eqs. (22) and (24), respectively. Id value also can be calculated from a known Hd using Eqn. (20). Ib is the difference between I and Id. is considered equal to . Then Eqns. (28) and (29) are used to calculate the hourly and daily radiation incident on a tilted surface. 4. Simulation Procedure for Monthly Average Daily Radiation on Horizontal and Tilted Surfaces The basic equation to be used is Eqn. (33). The first steps are to obtain Hdm/Hm and Rbm. The ratio Hdm/Hm is a function of Km and it can be obtained from Eqn. (26). For the mean day of any month n is known (Table 1) and can be calculated from Eqn. (4). For known values and , sunset hour angle s and monthly average of the maximum possible daily hours Nm or N are calculated from Eqns. (5) and (6), respec- tively. Then with assumed monthly average daily hours of bright sunshine (nd), the nd/Nm ratio is calcu- lated to be used in Eqn. (18). With known n and s, the extraterrestrial (theoretical) radiation H0m is calculated from Eqn. (16). The monthly average daily radiation Hm on a horizontal surface is calculated from Eqn. (18). Then Km is obtained from Eqn. (25). Equation (26) is used to calculate the Hdm/Hm ratio from Km. The calculation of Rbm requires sunset hour angle m. Equation (31) is used to calculate sunset hour angle s. Finally, global radiation on an inclined surface is calculated using Eqn. (33). In this way monthly aver- age daily radiation incidences on horizontal and tilted surfaces are calculated for the 12 months of the year. A computer program was written in BASIC program- ming language. 5. Measured and Simulated Solar Radia- tion Data A pyranometer we used to measure the solar radia- tion incident on a horizontal surface for several days in October 2004 from morning to evening. The pyra- nometer was connected to a personal computer via a data logger (Green klt 77) to collect the radiation data at one minute intervals using BASIC programming language. Sensitivity of the pyranometer was 10.5 mV/(cal.cm2.min) (Moll-Goregiatic type MS-61). The model described was applied to generate total global hourly, daily, and monthly average radiation on a hor- izontal and tilted surfaces. The comparison between the measured total (beam and diffuse) and the simulat- ed total (beam and diffuse) solar radiation on a hori- zontal surface versus time for clear days are shown in Figures 1 and 2. There was favourable agreement between the measured and predicted values on a day with clear skies. For example in October 5, 2004 day, actual meas- ured insulation on a horizontal surface was 8.04 33 Simulation of Solar Radiation Incident on Horizontal and Inclined Surfaces MJ/m2 and simulated one was 9.5 MJ/m2. Similarly in October 24, 2004 day, actual measured insulation was 8.6 MJ/m2 and simulated one was 8.9 MJ/m2. So the hourly and daily radiations incident on horizontal sur- faces can be simulated from the theoretical relation- ships. The measured monthly average daily radiation inci- dents on a horizontal surface of a flat plate collector were compared for a number of stations in Japan. There was favorable agreement between the measured and simulated values (Figure 3). It is expected that simulated values will show the same level of accuracy with the measured data of other places. Simulations show that during the paddy har- vesting period from mid-September to mid-October, about 14.7 MJ/(m2.day) of solar energy is available on Figure 1. Total (beam and diffuse) measured and simulated solar radiation a horizontal surface versus time for a clear day on October 5, 2004 Figure 2. Total (beam and diffuse) measured and simulated solar radiation on a horizontal surface versus time for a clear day on October 24, 2004 34 MS Basunia, H Yoshio and T Abe a tilted surface at Matsuyama, Japan, and that an almost similar amount (14.1 MJ/(m2.day)) of solar energy is also available during the major rice harvest- ing season in Bangladesh from November-December. This energy can be utilized for drying rough rice and other farm crops. The model has been developed on the basis of available relations for hourly, daily, and monthly average radiation on horizontal and tilted sur- Figure 3. Variations of measured and simulated monthly average daily global solar radiation with time of the year at three stations of Japan J F M A M J J A S O N D Month 25 20 15 10 5 0 Table 2. Simulated monthly average daily global solar radiation incident on horizontal and inclined surfaces at Matsuyama, Japan and Bangladesh Note: Tilted angle is considered as equal to the latitude of respective places 35 Simulation of Solar Radiation Incident on Horizontal and Inclined Surfaces faces. Bangladesh lies in a narrow range of latitude (20º - 26º N) and, therefore, the same level of solar radiation is expected to be incident throughout the country. 6. Conclusions A simulation model was developed on the basis of best available relations for daily, hourly and monthly average radiation on horizontal and tilted surfaces. The simulated hourly and daily solar radiation values were compared with the measured values and a favourable agreement was observed between them. The measured monthly average daily solar radiation was compared with the simulated values and a reasonable agreement was observed for a number of stations in Japan. Considering the narrow range of latitude and similar terrain and weather conditions, the same level of solar radiation is expected throughout Bangladesh. The amount of solar energy available during major rice harvesting season in Matsuyama, Japan, and Bangladesh are almost the same and can be used for drying rough rice and other farm crops. References Bala, BK, Zia-Uddin ATM (1990), Simulation of solar drying of rough rice. Proclamation of the 1st World Renewable Energy Congress, Reading, UK 938-942. Benford F, Bock JE (1939), Transaction of the Illumination Engineering Society 34:200-203. Collares-Pereira M, Rabl A (1979), The average dis- tribution of solar radiation-correlation between diffuse and Hemispherical and between day and hourly insolation values. J. of Solar Energy 22:155-162. Duffie JA, Beckman WA (1991), Solar engineering of thermal process. A Wiley-International Publication. John Wiley and Sons, New York. Erbs DG, Klein SA, Duffie JA (1982), Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation. J. of Solar Energy 28:293-300. Klein SA (1977), Calculation of monthly average inso- lation on tilted surfaces. J. of Solar Energy 19: 325-330. Liu, BYH, Jordan RC (1962), Daily insolation on sur- faces tilted toward the equator. ASHRAE Journal 3(10):53-58. Liu, BYH, Jordan RC (1963), The long-term average performance of flat-plate solar energy collectors. J. of Solar Energy 7:53-59. Lof GOG, Duffie JA, Smith CO (1966), World distri- bution of solar radiation. Engineering experimen- tation station report 21, University of Wisconsin, Madison. Orgill JF, Hollands KGT (1977), Correlation equation for hourly diffuse radiation on a horizontal sur- face. J. of Solar Energy 19:357-363. Page JK (1964), The estimations of monthly mean val- ues of daily total short-wave radiation of vertical and inclined surfaces from sunshine records for latitudes 400N-400S. Proc. of the UN conference on new sources of Energy 4:378-385. Reindl DT, Beckman WA, Duffie JA (1990), Diffuse fraction correlation. J. of Solar Energy 45:1-6. Steinfeld A, Segal I (1986), A simulation model for solar thin layer drying process. J. of Drying Technology 4:536-554.