Energy and Power

p-ISSN: 2163-159X    e-ISSN: 2163-1603

2016;  6(2): 29-38

doi:10.5923/j.ep.20160602.01

 

General Polynomial for Optimizing the Tilt Angle of Flat Solar Energy Harvesters Based on ASHRAE Clear Sky Model in Mid and High Latitudes

Alsadi Samer Yasin1, Nassar Yasser Fathi2, Amer Khaled Ali2

1Electrical Engineering Department, Faculty of Engineering and Technology, Palestine Technical University-Kadoorie, Tulkarm, Palestine

2Mechanical Engineering Department, Engineering and Technology Faculty, Sebha University, Brack, Libya

Correspondence to: Nassar Yasser Fathi, Mechanical Engineering Department, Engineering and Technology Faculty, Sebha University, Brack, Libya.

Email:

Copyright © 2016 Scientific & Academic Publishing. All Rights Reserved.

This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

Abstract

The aim of this study was to identify the opportunity of optimising the collection of solar energy as far as it is available in order to increase its utilisation. The second aim of this study was to enhance the performance of solar energy harvesters that depend on it through the appropriate determination of optimum solar collector tilt angles. In this study, ASHRAE clear sky model is applied to calculate the solar radiation; and the fundamental solar energy equations were programmed to determine optimum tilt angles in any location on the earth. The optimum tilt angle was presented in a polynomial form to enable the most convenient use of the function of latitude angle on a Julian day in a monthly, seasonal and an annual manner. According to the comparison between the obtained results with those of the local measured authoritative data and those of NASA published data, it can be safely recommended that these polynomials be used especially in mid and high latitudes ( > 20°) in the two hemispheres. The presented study could serve as a reference for the domestic solar electrical and thermal applications.

Keywords: Optimum tilt angle, Solar energy

Cite this paper: Alsadi Samer Yasin, Nassar Yasser Fathi, Amer Khaled Ali, General Polynomial for Optimizing the Tilt Angle of Flat Solar Energy Harvesters Based on ASHRAE Clear Sky Model in Mid and High Latitudes, Energy and Power, Vol. 6 No. 2, 2016, pp. 29-38. doi: 10.5923/j.ep.20160602.01.

1. Introduction

Power supply affects a lot a country’s economic activity. Currently, fossil fuel is the most important part of power supply in the world; however, for the consideration of the issue of environment pollution, renewable energy -including solar energy, wind energy, and biomass energy-is being paid more and more attention. Solar energy is the best choice for most of the areas all over the world [1] as it can be utilized through solar collector or photovoltaic (PV) cell. Many researchers have been dealing with this subject, but they almost were able to determine the optimum angle for limited locations only [1-14]. In addition, a comprehensive research like this current study is not dealt with in the related literature; In fact; this was the motive to carry out this research which extends comprehensively to present formula for the optimum tilt angle of flat solar harvesters for any location and any time.
To maximize the collected energy, proper installing of the collector is really needed. One way is to install the solar collectors in the correct tilt and orientation angle, in which they would obtain the maximum insolation over a specific period of time. The tilt angle depends mainly on the position of the sun and, therefore, is different from one location to another all over the world; in regard to the best orientation angle, it is advised to be directed towards the equator. There are already plenty of investigations dealing with this subject to optimize solar power systems according to the correct tilt angle or orientation. Many authors have provided empirical models to calculate the optimum tilt angle by searching for the maximum total solar radiation on the collector surface. In reference to a specific period of time and purpose, daily, monthly, seasonally or yearly values have been calculated [15-17]. Elsayed in [18] also presented an analytical model based on long-term averaging of solar data. He outlined values of optimum tilt angles given in different literature and conducted that value of tilt angle that can be recommended.
In fact there is a wide range of tilt (±20) which is dependent on the applied model and the location. Some authors noted the existence of a correlation between the optimum tilt angle and the latitude. Frequently, it is recommended to apply the rule of thumb, in which the yearly optimum tilt angle is about (L: latitude) and a difference of tilt with about 10° would hardly affect the performance. Authors of [19] determined monthly optimum tilt angles for Izmir, Turkey; they found the optimum tilt angle to be equal to U throughout the year, while for summer and for winter was suggested, these values were suggested also by [20]. They advised to mount the solar collector at the monthly average tilt angle. During the last decade there have also been investigations by using simulation software. This software which was developed in order to simulate an entire solar power plant, takes the most influential parameters into account. The soft ware usually possesses a database of monthly mean global radiation data and different empirical models: for example authors of [21] applied the simulation software TRNSYS to calculate the optimum tilt angle for Cairo, Egypt. What they did was calculating the monthly mean of solar radiation data and compared it with the output power of solar cells. They stated the yearly optimum tilt angle to be

2. Modelling of Solar Radiation on Tilt Surface

The moment the solar radiation, which is coming from the direction of the sun, reaches the earth’s surface without being significantly scattered is called direct normal irradiance (or beam irradiance). Some of the scattered sunlight is scattered back into space and some of it also reaches the surface of the earth. The scattered radiation reaching the earth’s surface is called diffuse radiation. Some radiation is also scattered off the earth’s surface and then re-scattered by the atmosphere to the observer. This is also part of the diffuse radiation which the observer can see. This amount can be significant in areas in which the ground is covered with snow or reflectors. The total solar radiation on a horizontal surface is called global irradiance and it is the sum of incident diffuse radiation plus the direct normal irradiance projected onto the horizontal surface. If the surface under study is tilted with respect to the horizontal, the total irradiance is the incident diffuse radiation plus the direct normal irradiance projected onto the tilted surface plus ground reflected irradiance that is incident on the tilted surface.
As recommended by ASHRAE (1985) and presented in all text books of solar energy [22, 23], hourly global radiation hourly beam radiation hourly diffuse radiation and hourly reflect radiation on the inclined surface on a clear day are calculated using the following expression:
(1)
(2)
(3)
(4)
Where: is the beam normal radiation A is the apparent solar-radiation constant in B is the atmospheric extinction coefficient, and C is the diffuse sky factor and there values are tabulated in table 1 for a widely range of latitudes
Table 1. ASHRAE values of clear sky model parameters A, B and C [23]
     
and are coefficients and is ground Albedo and was assumed to be 0.5. is the ratio between global solar energy on a horizontal surface and global solar energy on a tilt surface. Meanwhile, is the view factor between the sky and the tilt surface and is the view factor between the ground and the surface that tilted at angle from the horizontal. These coefficients evaluated from:
(5)
for isotropic diffuse:
(6)
(7)
where: is the solar incident angle and it is given by:
(8)
where: is latitude, the angular location north or south of the equator, north positive;
slope of the surface; is the surface azimuth angle, with zero due south, east negative, west positive, is the hour angle, morning negative, zero at noon and afternoon positive; and is solar declination angle, it can be found from:
(9)
where: is the Julian day. Meanwhile, is the solar zenith angle, which equal to:
(10)

3. Results and Discussions

A comprehensive computer program in FORTRAN has been created in order to calculate the hourly solar radiation incident on a surface when the tilt angle is changed by interval of The surface azimuth is either for south facing or for north facing. In regard to the daily optimization, the basis is the daily total hourly radiation, for the monthly optimization the basis is the monthly total daily and for annual optimization the basis is the annual total monthly.
It is well known that a unique surface tilt angle which exists for each time, and which corresponds to each latitude angle for a particular day through which the solar radiation is at peak value. Accordingly, this study enhances these two parameters as independent variables for the daily, monthly, seasonal and annual optimum angle prediction.
The seasonal average tilt angle was calculated by finding the average value of the tilt angle for each season, but the implementation of this condition requires the collector tilt to be changed four times a year. The process of adjusting the tilt angle to its monthly optimum values throughout the year does not seem to be practical, and at the same time rises the consideration of changing the tilt angle once seasonal.
The obtained results from the FORTRAN program treated by MATLAB in order to get out polynomials those fitted the daily, monthly, seasonal and annual optimum angles.
The general form of the polynomial that prescribed the optimum tilt angle is expressed in fifth order three dimension as:
(11)
The coefficients and the variables and are defined below in table 2.
Table 2. The definition of the coefficients P00 ...P05 and the independent variables x, y for the offered polynomial
The obtained results were plotted and tabulated for comparison with other measurements and data recommended by NASA.
Figure 1 presents the Daily optimum angles for various latitudes for a complete calendar year. The abscissa presents the Julian day The markets present the calculated results and the dash lines present the daily polynomial. Obviously, depending on figure 1, it can be stated that the polynomial is well fitted to the calculated data. The positive value depicted means that the surface was inclined towards the equator, where as the negative value means that the surface was inclined towards the North Pole. The slopes were within the range of for any location on the northern hemisphere.
Figure 1. Daily optimum angle as a function of the Julian day n and the latitude angle L. The marks refer to the calculated results and the dash lines refer to the polynomial Sopt = f(n, L)
With the adoption of the long-term solar radiation data, the optimum tilt angle of a surface by using the monthly total daily solar irradiation on diversely latitudes was simulated for 12 months. Figures 2 presents the optimum angles for monthly tracking from January to December. From figure 2, the step-like lines present the calculated results, and the dash curves illustrate the polynomial equation that fitted the calculated data. Results were tabulated in table 4, for many cities from other references beside NASA data available on the internet to ease the comparison. Table 3 shows a considerably agreement in results obtained by the offered polynomial with the local measurements specially for mid and high latitudes even in some times better than NASA's data. Unfortunately, for low latitudes the model needs to more arrangements.
Figure 2. Monthly optimum angle as a function of the Julian day n and the latitude angle L. The solid lines refer to the calculated results and the dash lines refer to the polynomial Sopt = f(n, L)
Table 3. Values of monthly optimum angles from many sources
Table 4. Description of Julian day for northern and southern hemispheres that used in polynomial (11)
     
For the southern hemisphere another arrangement must be taken to estimate the monthly optimum angles; one can use the same polynomial with the same latitude but the months will be reversed. That means December in the northern hemisphere will be June in the southern hemisphere. The months must be shifted every 6 months after the ordinary calendar. Table 4 described the Julian day for both hemispheres for monthly tilt angle optimization, and the optimum tilt angle for southern hemisphere is negative sing of the optimum tilt angle for northern hemisphere for the same latitude.
(12)
Sometimes the cause of the difference in some values for some sites is a result of the large step of tilt angle 10°, or because the researchers had not included angles greater than 60° or less than 10°. Anyhow, it can be obviously stated from table 3 the agreement in values obtained from the mathematical model proposed in this study compared with the those results presented in other previous studies.
The study of seasonal optimum tilt angle for maximising energy collection by the absorber has been carried out. In literatures, there are two classifications for seasons: Heating and cooling group interesting with two seasons heating season and cooling season. Another classification may be adopted for solar energy group dependent on the position of the sun in the sky, they classified seasons to four seasons in 12 months. This study adopted the second classification. Each season consists of 91 days, therefore, winter (6th November to 4th February), spring (5th February to 5th May), summer (6th May to 5th August), and autumn (6th August to 5th November).
These seasonal is valid for northern hemisphere. In the southern hemisphere the situation is reverse the summer season (6th November to 4th February), autumn (5th February to 5th May), winter (6th May to 5th August), and spring (6th August to 5th November). Figure 3 shows the optimum seasonal tilt angle based on the total daily solar radiation reaching on a tilt surface. The step-like lines present the calculated optimum tilt angle for seasonal optimization. The simulated results have agreed to the local available measurements with similar pattern to reflect the trend at the site. From the figure, one polynomial equation were proposed to depict the trend of seasonal optimum as functions of the latitude and the day. The equation shows a good fit to the output as it obvious from the figure 4, in where the dash curves present the polynomial. A comparative results are tabulated in table 5. Unfortunately, there are no enough results to include in the table; moreover, some researches classified the season into only two seasons cold and hot seasons for heating and cooling of buildings purposed.
Figure 3. Seasonal optimum angle as a function of the Julian day n and the latitude angle L. The solid lines refer to the calculated results and the dash lines refer to the polynomial Sopt = f(n, L)
Table 5. Values of seasonal optimum angles from many sources
     
The annual optimum tilt angle has been computed based on the estimated annual total monthly solar radiation. The results illustrated graphically in figure 4 and tabulated in table 6 for comparative purpose. Results showed that the optimum values were almost positive for the northern hemisphere and ranged facing to the south (equator). The value was similar and agreed with the optimum slope presented by other researchers and for those of NASA, specially for high latitudes even in some times better than NASA's data. Additionally, the computed results were fairly consistent with the general rule that the yearly optimal tilt angle was about the latitude of the location facing to the equator, specially for middle latitudes. For the southern hemisphere the optimum tilt angle will be negative sign of the northern hemisphere at the same latitude.
Figure 4. Annual optimum angle as a function of the latitude angle L. The marks refer to the calculated results and the dash lines refer to the polynomial Sopt = f( L)
Table 6. Values of annual optimum fixed angles from many sources and NASA
     

4. Conclusions

We have successfully applied the ASHRAE clear sky model to create comprehensive three dimension polynomials for daily, monthly, seasonal and annual optimum tilt angle for diversely latitudes. According to the comparison between the obtained results with local measured authoritative data and NASA published data, it can safely recommended to use of these polynomials especially in mid and high latitudes in the two hemispheres.

References

[1]  Feng-Jiao Liu, Tian-Pau Chang, Optimizing the Tilt Angle of Solar Collector under Clear Sky by Particle Swarm Optimization Method, Journal of Applied Science and Engineering, Vol. 18, No. 2, 2015, pp. 167-172.
[2]  T. Khatib, A. Mohamed, M. Mahmoud, K. Sopian, Optimization of the Tilt Angle of Solar Panels for Malaysia, Energy Sources, Part A: Recovery, Utilization, and Environmental Effects , Taylor & Francis Group, LLC, 37, 2015, pp. 606–613.
[3]  M. Benghanem, Optimization of tilt angle for solar panel: Case study for Madinah, Saudi Arabia, Applied Energy 88, 2011, pp. 1427–1433.
[4]  Jasmina Skerlić, Jasna Radulović, Danijela Nikolić, and Milorad Bojić, Maximizing performances of variable tilt flat-plate solar collectors forBelgrade (Serbia), J. Renewable Sustainable Energy 5, 2013, doi: 10.1063/1.4819254.
[5]  Belkacem Ihya, Abdellah Mechaqrane, Rachid Tadili, Mohammed Najib Bargach, Optimal tilt angles for solar collectors facing south at Fez city (Morocco), Journal of Natural Sciences Research, Vol.4, No.10, 2014, pp. 120-127.
[6]  Felix A. Uba and Emmanuel A. Sarsah, Optimization of tilt angle for solar collectors in WA, Ghana, Advances in Applied Science Research, 4(4), 2013, pp. 108-114.
[7]  M. Jamil Ahmad, G.N. Tiwari, Optimization of Tilt Angle for Solar Collector to Receive Maximum Radiation, The Open Renewable Energy Journal, 2, 2009, pp. 19-24.
[8]  Dhrupad Sarm, Parimal Bakul Barua, Diganta Hatibaruah, Optimization of Glazing Cover Parameters of a Solar Flat Plate Collector (FPC), International Journal of Engineering Trends and Technology (IJETT) – V. 14 No 2, 2014, pp. 74-80.
[9]  Ying-Pin Chang1, Der-An Wang, Optimization of Tilt Angle for Photovoltaic Modules Based on the Neural-Genetic Algorithm, Journal of Nan Kai, Vol.7, No.1, 2010, pp. 57-70.
[10]  Jayashri Vajpai Navaneet Maru Model based Optimization of Tilt angle for Solar PV Panels in Jodhpur, International Journal of Computer Applications, 2014, pp. 33-39.
[11]  Qais Azzam Khasawneh, Qatada Abdullah Damra, Omaymah Husni Bany Salman, Determining the Optimum Tilt Angle for Solar Applications in Northern Jordan, Jordan Journal of Mechanical and Industrial Engineering, V. 9 No. 3, 2015, pp. 187 - 193.
[12]  Arslanoglu A., Optimization of tilt angle for solar collectors: A case study for Bursa, Turkey, International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol:10, No:5, 2016, pp. 517-520.
[13]  Hasan M.S., Optimization of surface orientation angles to receive maximum solar radiation at Sebha city in Libya, M.SC. thesis, Sebelas Maret University, Surakarta, 2015.
[14]  Nassar Y.F., Simulation of solar tracking systems, Energy & Life Journal, No 21, March 2005, pp. 81-90.
[15]  EL-Kassa by, M.M., Monthly and daily optimum tilt angle for south facing solar collectors; theoretical model, experimental and empirical correlations. Solar & Wind Technology 5, 1988, pp. 589–596.
[16]  Soulayman, S.S., On the optimum tilt of solar absorber plates. Renewable Energy 1, 1991, pp. 551–554.
[17]  Ibrahim, D., Optimum tilt angle for solar collectors used in Cyprus. Renewable Energy 6, 1995, pp. 813–819.
[18]  Elsayed, M.M., Optimum orientation of absorber plates. Solar Energy 42, 1998, pp. 89–102.
[19]  Gunerhan, H., Hepbasli, A., Determination of the optimum tilt angle of solar collectors for building applications. Building and Environment 42, 2007, pp. 779–783.
[20]  H. P. Garg, Treatise on solar Energy, vol. 1, Fundamentals of solar energy, New York, Willy, 1983.
[21]  Hussein, H.E.-G.H.M.S., Ahmad, G.E., Performance evaluation of photovoltaic modules at different tilt angles and orientations. Energy Conversation and Management 45, 2004, pp. 2441–2452.
[22]  Duffie J. A. and Beckman W. A., Solar Engineering of Thermal Processes, 3ed., John Wiley & Sons, Inc., 2006.
[23]  Nassar Y. F., Solar energy engineering active applications, Sebha University, Libya, 2006. "Arabic language".
[24]  Hassane Darhmaoui, Driss Lahjouji, Latitude Based Model for Tilt Angle Optimization for Solar Collectors in the Mediterranean Region, Science Direct, Energy Procedia 42, 2013, pp. 426 – 435.
[25]  Khai Mun Ng, Nor Mariah Adam, Othman Inayatullah, Mohd Zainal Abidin Ab. Kadir, Assessment of solar radiation on diversely oriented surfaces and optimum tilts for solar absorbers in Malaysian tropical latitude, Int. J. Energy Environ. Eng., 2014, 5:75.
[26]  Yueren Wang, Chenfei Bi, Wen Lou, Yueren Wang, Chenfei Bi, Wen Lou, Collector solar hot water installation angle optimization in severe cold area, International Conference on Applied Science and Engineering Innovation (ASEI 2015).
[27]  Ioannis M. Michaelides, Computer simulation and optimisation of solar heating systems for Cyprus, PhD thesis, Faculty of Engineering and Science, Westminster Research http://www.westminster.ac.uk /west minster research, 1993.
[28]  Y. N. Udoakah1, N. I. Okpura, Determination of optimal tilt angle for maximum solar insolation for pv systems in Enugu-southern Nigeria, Nigerian journal of technology, vol. 34 no. 4, 2015, pp. 838 – 843.
[29]  Igor Shesho, Done Tashevski, Developing simulation application using graphical programming language for optimization of solar collector position, Mechanical Engineering – Scientific Journal, Vol. 31, No. 1–2, 2013, pp. 63-75.
[30]  Z. Bouhssine1, M. Faraji1, M. Najam1 and M. El Alami, Numerical Investigations of an Integrated Phase-Change-Material Solar Collector, Tech Science Press FDMP, vol.11, no.1, 2015, pp. 63-85.
[31]  Adnan Shariah, M-Ali Al-Akhras, I.A. Al-Omari, Optimizing the tilt angle of solar collectors, Renewable Energy 26, 2002, pp. 587–598.
[32]  Shaoning Wang, Bo Hong, Optimum Design of Tilt Angle and Horizontal Direction of Solar Collectors under Obstacle’s Shadow for Building Applications, Journal of Building Construction and Planning Research, 3, 2015, pp. 60-67.
[33]  Jože Rakovec, Klemen Zakšek, Kristijan Brecl, Orientation and Tilt Dependence of a Fixed PV Array Energy Yield Based on Measurements of Solar Energy and Ground Albedo – a Case Study of Slovenia, Energy Management Systems, 2011, pp.145-160. www.intechopen.com.
[34]  S. Armstrong, W.G. Hurley, A new methodology to optimise solar energy extraction under cloudy conditions, Renewable Energy 35, 2010, pp. 780–787.
[35]  Ahmad Qasaimeh, Solar Energy Optimization through Seasons: Case Study in Jordan, Smart Grid and Renewable Energy, 3, 2012, pp. 275-281. www.SciRP.org/journal/sgre.
[36]  Fawaz Sultan, Firas Aziz Ali, Tariq Khalid A. Razaq, Tilt Angle Optimization of Solar Collectors for Maximum Radiation in Three Iraqi cities, International Journal of Engineering and Industries (IJEI),Volume3, no 4, December 2012, pp. 99-107.
[37]  NASA surface meteorology and solar energy - available tables. https://eosweb.larc.nasa.gov/cgi-bin/see/grid.cgi?