Received 10 December 2015; accepted 17 April 2016; published 20 April 2016
The social and economic development of a nation is assessed on the basis of energy consumption per capita. In advanced countries 7.35 KW per person is required to maintain the standard of living and comfort. For Europe it is 5.02 KW per person, 1.30 KW per person in the developing countries and only 0.228 KW per person in Bangladesh. Energy consumption per capita over the world is 2.37 KW  . So, the countries poor in energy are in an extremely unhappy position. The world population will reach approximately 12 billion by the year 2100. If the current trend in technological progress and innovation continues, the demand for energy then will be five times greater than what it is now. Thus it is obvious that the present exhaustible forms of conventional energy resources, such as wood, cool, natural gas and oil will not be able to meet this every increasing energy demand. Besides the fast rate of consumption of conventional energy resources, the other very important factor responsible for energy crisis is the odd distribution and consumption of these resources. This is to be mentioned here that one-third of the world population consumes about four fifths of the total world energy while a very meager share is left to the majority of the population living in developing or underdeveloped countries.
The present day energy crisis has therefore resulted in the search for alternate energy resources in order to cope with the drastically changing energy picture of the world. With the rapid industrialization and fast growth rate of energy consumption, the world will soon run out of its conventional fuel. This problem demands an urgent and serious attention, especially by the inhabitants of developing and underdeveloped countries, not only to the husbanding of present resources, but also to the short and long-term development of future energy resources. This is important because the social and economic development of these countries heavily leans on the consumption and harnessing of their energy resources. Out of the available alternate energy resources, solar energy has emerged as one of the most powerful, clean, safe and virtually inexhaustible sources of energy; world can depend on, in the years to come.
Solar energy is quite simply the energy produced directly by the sun and collected elsewhere, typically on the earth. The geographical location of Bangladesh lies between 20˚34' and 26˚38' North latitude and 88˚01' to 92˚41' East longitude, and is endowed with abundant sunshine for at least eight months of the year  . It has a total land area of 147,570 square km. On an average 5 KW/m2 per day solar energy falling on this land over the year, in certain regions the bright sunny days may be more. Even if one percent of this land is used to harness solar energy for electricity generation at an overall efficiency of 12%, 7200 MWh/year electricity is generated  . This is a redundant amount of electricity.
The availability of solar energy in a particular location depends basically on the astronomical, geographic and climatic factors, with greatest influences on atmospheric transmissivity given water vapor concentrations, aerosols and clouds   . Bangladesh lies between 15˚N and 35˚N known as the most favorable belt embracing the regions endowed with the most suitable conditions for solar energy applications. In this respect, Bangladesh is fairly rich as far as solar radiation is concerned.
2. Theoretical Analysis of Solar Radiation
In this study, the value of solar radiation has been estimated on the basis of meteorological parameters of the last 30 years collected by Bangladesh Meteorological Department with the help of different empirical formulae. Solar radiations (global, diffuse and direct solar radiation) over 10 stations (districts) of Bangladesh has been predicted by using the value of monthly averaged of different parameters e.g. the monthly averaged value of maximum temperature, minimum temperature, humidity and sunshine hours for the period of 1980-2009 years. The analysis has been done by using MS Excel, Origin, SPSS.
2.1. Estimation of Monthly Averaged Daily Radiation (Global, Diffuse and Beam) on Horizontal Surface
Generally to estimate the monthly averaged daily global (G) and diffuse (D) radiation on a horizontal surface, following Garg and Garg, M. Hussain  correlation has been used.
where, Go is the extra-terrestrial radiation, N is the maximum duration for which the Campbell-Stakos sunshine recorder is active, n is the sunshine duration and w is the atmospheric water content.
The total solar radiation on a horizontal surface is equal to the sum of diffuse and beam radiation on a horizontal surface. Therefore we can find out beam radiation (I) by subtraction method
Diffuse solar radiation represents the shortwave energy of solar origin scattered downwards by gas molecules, aerosols, water vapor and clouds in the atmosphere. The variability in the amount and type of cloud has a dominant role on the value of diffuse radiation. The solar scientist gave a theoretical relation between diffuse solar radiation (D) and the clearness index (KT) which is given by
where, c and d are the correlation coefficients. The values of c and d are obtained from monthly averaged data of all months of the year for a location. The term clearness index KT is defined as the ratio of monthly averaged radiation on a horizontal surface to the monthly averaged daily extraterrestrial radiation i.e.
An angstrom like correlation using sunshine hour data
has been used to compute diffuse radiation in terms of global radiation for different locations of Bangladesh. Here c' and d' are the constant and regression coefficient.
In case of direct solar radiation, a page likes regression using sunshine hour data
and another angstrom like regression using sunshine hour data
have been used to estimate monthly average direct solar radiation over Bangladesh. Here e, f, , are the regression coefficients to be determined.
2.2. Estimation of Monthly Averaged Daily Solar Radiation (Global, Diffuse and Beam) on Tilted Surface
The method developed by Liu and Jordan  of estimating the average daily radiation for each month on a tilted south-facing surface is described in Klien’s review paper  . It is also extended to calculate the radiation for off-south orientation (east or west of south). The average daily radiation on a tilted surface is given as
where, R is the ratio of the daily average radiation on a tilted surface to that on a horizontal surface for each month. R can be estimated by considering individually the direct, diffuse and reflected components of radiation incident on the tilted surface. Assuming diffuse and reflected radiation to be isotropic, Liu and Jordan  suggest R as equal to
where, D is the monthly averaged daily diffused radiation, Rb is the ratio of the average beam radiation on tilted surface to that on the horizontal surface, s is the tilt angle from the horizontal and ρ is the ground reflectance, usually taken to be 0.2 when there is no snow.
From a surface facing directly southwards
where is the sunset hour angle for tilted surface given by
where is the sunrise hour angle in radian and is given by
where, is the sun’s declination which is the angular distance of the sun’s rays north or south of the equator. Also
where n is the day of the year.
3. Results and Discussion
3.1. Estimation of Solar Radiation on Horizontal Surface
Calculation of global, diffuse and beam radiation for each station has been evaluated by using S. P. Garg and H. N. Garg formula  , M. Hussain formula  and subtraction method. Monthly averaged daily global, diffuse and beam radiation in different stations of Bangladesh are shown in Figures 1-3.
Figure 1. Illustrating the variation of monthly averaged global (G) solar radiation over selected stations.
Figure 2. Illustrating the variation of monthly averaged diffuse (D) solar radiation over selected stations.
Figure 3. Illustrating the variation of monthly averaged beam (I) solar radiation over selected stations.
The dependence of global solar radiation on maximum temperature is intensified by sunshine hour. From all the curves, it is clear that the first peak in the global solar radiation occurs during summer (April/May). In first period; both sunshine hour and temperature are high. But the second peak occurs in August (autumn) which is not so prominent, due to high temperature but short sun shining period. Again, during winter (November/ December), though there is enough sun shining periods but the temperature is low. Therefore, it results in low global solar radiation. The maximum value of global radiation is 6.34 Kwh/m2 day in Chittagong for the month of April when the bright sunshine hour, average temperature and humidity are (8.1 hours), (28˚C) and (77%) respectively, contrary to the minimum value of global radiation is 3.74 Kwh/m2 day in Chittagong for the month of November when the bright sunshine hour, average temperature and humidity are (5.2 hours), (25˚C) and (76%) respectively. The average maximum value of global radiation has been found in Cox’s Bazar 5.22 Kwh/ m2 day and the minimum value is 4.75 Kwh/m2 day in Mymensingh.
The diffuse solar radiation depends on relative humidity and atmospheric water content. It increases with the decrease of sun shining hour and increase of atmospheric water content. Therefore, the diffuse radiation is maximum in June/July (rainy) and minimum in December/January (winter). The maximum value of diffuse radiation is 3.83 Kwh/m2 day in Sylhet for the month of June when the bright sunshine hour, average temperature and humidity are (3.77 hours), (28.07˚C) and (86.1%) respectively, and the minimum value is 1.37 Kwh/m2 day in Rajshahi for the month of December when the bright sunshine hour, average temperature and humidity are (8.05 hours), (19.01˚C) and (77.6%) respectively.
The direct solar radiation is directly related to sunshine duration. The direct solar radiation is, therefore, maximum in March/April (summer) and minimum in July (rainy). The peak value of direct radiation is 3.64 Kwh/m2 day in Chittagong for the month of March when the bright sunshine hour, average temperature and humidity are (8.53 hours), (25.60˚C) and (73.28%) respectively, and the minimum value is 0.92 Kwh/m2 day in Sylhet for the month of July when the bright sunshine hour, average temperature and humidity are (3.64 hours), (28.40˚C) and (86.9%) respectively.
3.2. Monthly Variation of Clearness Index
From Figure 4 of clearness index at different stations, we can predict four solar seasons for Bangladesh. These seasons are as follows:
Season-1: December, January, February; Season-2: March, April, May; Season-3: June, July, August; Season-4: September, October, November.
3.3. Determination of Regression Coefficient and Constant of Different Stations
Table 1 shows the values of the regression coefficients and constants a, b, c, d, c′, d′, e, f, e′, f′ for ten stations of Bangladesh. It is evident that, the values of the coefficient “a” vary from 0.232 to 0.257, while the coefficient “b” varies from 0.487 to 0.573. Using the average values of the coefficients a and constant b for ten stations, a approximate equation recommended for the estimation of monthly averaged global solar radiation for any area
Figure 4. Monthly variation of clearness index.
Table 1. The values of regression coefficient and constants for different stations are shown in following table.
over Bangladesh, even where the radiation data is missing or unavailable, is:
which implies that about 24% of extraterrestrial radiation over Bangladesh filtrates the atmosphere into on a fully cloudy day (n = 0) and about 76% on a clear sky day respectively.
The maximum and minimum values of other regression coefficients and constants c, d, c′, d′, e, f, e′ and f′ are (1.547 and 1.37), (−1.616 and −1.937), (0.389 and 0.377), (−0.205 and −0.244), (−0.37 and −0.547), (1.983 and 1.616), (−0.12 and −0.156) and (0.817 and 0.692) respectively.
Using the average values of the coefficients and constants in Equations (15-18), the linear equations recommended for the estimation of monthly averaged diffuse and direct solar radiation all over Bangladesh are:
3.4. Calculation of Monthly Averaged Daily Solar Radiation on Tilted Surface
In the estimation of solar radiation on a tilted surface the angles of inclination have been chosen from 0˚ to 90˚ with an interval of 5˚. Monthly averaged daily solar radiation on tilted surface for all the months of the year for different stations is shown in Figure 5.
Figure 5. Illustrating the variation of monthly averaged daily total solar radiation for various inclination angles of different stations of Bangladesh.
Divergent value of solar radiation is found due to various inclination angles for different month. An angle for which solar radiation becomes highest is called as optimal tilt angle which has been discussed hereafter. It is apparent that the highest value of annual total solar radiation is found for 20˚ and 25˚ inclination angles in different stations of Bangladesh and the lowest for 90˚ inclination angle in all the stations of Bangladesh.
3.5. Conversion Factor
Conversion factor R implies the ratio of radiation on tilted surface GT to radiation on horizontal surface. The value of R depends upon the value of GT not on G. If the value of GT increases R will also increase and if decreases R will decrease because the value of G is always fixed. The value of GT varies from month to month and angle to angle as well as location to location. The highest value of R is found on optimal tilt surface for every distinct month. The value of R is fluctuated highly with an increase of different angles when applied consecutively. A typical conversion factor for Sylhet station is presented below by Table 2 and Figure 6 and Figure 7.
3.6. Optimal Tilt Angle
Optimal tilt angle refers to such angle for which the value of solar radiation on that surface becomes highest. Figure 8 shows optimal tilt angle for every concerned station for every month, and Figure 9 shows optimal tilt angle for few stations of Bangladesh for the year round.
In January, the optimal tilt angle is 40˚ for all the stations except Sylhet (45˚). In February, the optimal tilt angle is 35˚ for the stations Dhaka, Rajshahi, Sylhet, Rangpur, Mymensingh and 30˚ for the stations Chittagong,
Figure 6. Conversion factor R, for various inclina- tion angles for all month of the year.
Figure 7. Conversion factor R, for various months of the year.
Figure 8. Monthly optimal tilt angle for different stations.
Table 2. Values of conversion factor R for various inclination angles in Sylhet station.
Khulna, Barisal, Comilla and Cox’s Bazar. In March, the optimal tilt angle is 25˚ excluding Chittagong (20˚). In April, the optimal tilt angle is 10˚. In May, the optimal tilt angle is 5˚ except two stations (Chittagong, Cox’s Bazar) where their optimal tilt angle is 0˚. In June & July, the optimal tilt angle is 0˚ for all the stations. In August, the optimal tilt angle is 5˚ for all the stations except one station (Rangpur) where its optimal tilt angle is 10˚. In September, the optimal tilt angle is 15˚. In October, the optimal tilt angle is 30˚ for all the stations except 3 stations (Chittagong, Barisal, Cox’s Bazar), where their optimal tilt angle is 25˚. In November, the optimal tilt angle is 40˚ for all the stations except Chittagong (35˚). In December, the optimal tilt angle is 45˚ for all the stations.
The optimal tilt angle is mainly depends on Sun’s declination and latitude. When the sun’s declination increases from December solastic to June solastic, then the value of optimal tilt angle decreases. After June solastic the declination decreases and the value of optimal tilt angle increases. In December the optimal tilt angle is 45˚ when the declination is −21.16˚, while the optimal tilt angle is 0˚ for June & July when the declination is from +23.17˚ to +21.37˚. For the same month, the optimal tilt angle varies with latitude. The value of optimal tilt angle also depends on weather condition. That is why, the value of ground reflected radiation rate is higher in fine days than cloudy days. So the value of optimal tilt angle is higher in fine days than cloudy days.
3.7. Monthly Averaged Daily Beam, Diffuse, Ground Reflected and Total Solar Radiation on Tilted Surface in Dhaka
Total solar radiation on tilted surface consists of ground reflected, diffuse and beam radiation. The value of ground reflected radiation increases with the increase of the inclination angle. For this reason the lowest value of ground reflected radiation is found for 0˚ angle and the highest for 90˚ angle for the whole year and all stations as aforesaid in Bangladesh. The value of diffuse radiation always decreases with the increase of inclination angle owing to the value of scattered and reflected radiation on tilted surface decreases. That is why; the highest value of diffuse radiation is seen for 0˚ angle and the lowest for 90˚ inclination angle, for all month in the year and all stations in Bangladesh. The value of beam radiation relies upon various inclination angles and the highest value of beam radiation is found on the optimal tilt surface. Figure 10 shows the variation of different radiation for various inclination angles in Dhaka station.
Figure 9. Yearly optimal tilt angle for different stations.
Figure 10. Variation of ground reflected, diffuse, beam and total solar radiation for various inclination angles in Dhaka for month January (a), April (b), July (c) and October (d).
3.8. Comparison of Solar Radiation on Horizontal Surface and Tilted Surface at Dhaka Station
In January, we have found the value of solar radiation is 5.058 Kwh/m2 day for optimal tilt angle 40˚. On the other hand, the value of solar radiation is 4.149 Kwh/m2 day on horizontal surface. Therefore, we get 22% gain on optimal tilt surface with respect to horizontal surface. For the other months solar radiation on optimal tilt surface and horizontal surface, and gain also are shown in Table 3. The average value of solar radiation is 5.337 Kwh/m2 day for the optimal tilt angle 20.833˚ whereas the average value of solar radiation is 4.952 Kwh/m2 day on the horizontal surface. So, the average gain is 8.71% on optimal tilt angle with respect to horizontal surface for the whole year.
In January, we have found the value of solar radiation is 4.831 Kwh/m2 day on fixed 20.83˚ tilt surface. On the other hand, the value of solar radiation is 4.149 Kwh/m2 day on horizontal surface. Therefore, we get 16.43% gain on fixed tilt surface with respect to horizontal surface. For the other months solar radiation on fixed tilt surface and horizontal surface, and gain also are shown in Table 4. The average value of solar radiation is 5.213 Kwh/m2 day on fixed 20.83˚ tilt surface whereas the average value of solar radiation is 4.952 Kwh/m2 day on the horizontal surface. So, the average gain is 5.977% on fixed tilt surface with respect to horizontal surface for the whole year. Variable optimal tilt surface gives another interesting result in terms of solar radiation. Solar radiation on variable optimal tilt surface is equal or higher than that of solar radiation on horizontal surface or solar radiation on fixed 20.83˚ tilt surface. Figure 11 shows the maximum value of solar radiation on variable optimal tilt surface.
Table 3. Extra Solar radiation gain on variable optimal tilt surface from horizontal surface of Dhaka station.
Table 4. Extra solar radiation gain on fixed 20.83˚ tilt surface from horizontal surface of Dhaka station.
Figure 11. Comparison in terms of irradiance per day among horizontal surface, fixed tilt surface, variable optimal tilt surface.
Energy is very important part for every country and every person. But energy crisis is vital problem in the world at present. Our Bangladesh also faces mostly energy problem now. At this energy situation solar energy can be mitigated our energy problem. The research work reported in this paper will help regard state of knowledge relating to solar radiation on horizontal and tilted surface and estimate the solar energy potential for practical and efficient utilization in Bangladesh. That is why it will encourage the interested personnel dealing with promotion of solar energy throughout the country to install economically viable solar device by taking into account above- mentioned observation under which increased solar energy can be found. Ultimately it may mitigate energy crisis and thereby reduce sufferings of the general people as well as lessen investment crisis in Bangladesh.
Authors thank Dr. Md. Nurul Islam, Dept. of Physics, University of Chittagong for helpful discussions and Mr. Md. Abdul Mannan, Meteorologist, Bangladesh Meteorological Dept. for the cooperation in collecting data.
 Codato, G., Oliveira, A.P., Soares, J., Escobedo, J.F., Gomes, E.N. and Pai, A.D. (2008) Global and Diffuse Solar Irradiances in Urban and Rural Areas in Southeast Brazil. Theorical Applied Climatology, 93, 57-73.
 Escobedo, J.F., Gomes, E.N., Oliveira, A.P. and Soares, J. (2009) Modeling Hourly and Daily Fractions of UV, PAT and NIR to Global Solar Radiation under Various Sky Conditions at Botucatu, Brazil. Applied Energy, 86, 299-309.
 Klein, S.A. (1977) Calculation of Monthly Average Insolations on Tilted Surfaces. Solar Energy, 19, 325-329.
 Nahar, S., Golam Sarwar, A.T.M. and Chowdhury, S.A. (2009) A Theoretical Analysis of Optimizing Solar Irradiance: Bangladesh Perspective. 1st International Conference on the Developments in Renewable Energy Technology, Dhaka, 17-19 December 2009, 161-164.