# An Algorithm to Determine the Optimum Tilt Angle of a Solar Panel from Global Horizontal Solar Radiation

Authors: Emanuele Calabrò
Journal: Journal of Renewable Energy
 (2013)

---

## Abstract

This paper proposes an algorithm to calculate the optimum tilt angle of solar panels by means of global horizontal solar radiation data, provided from Earth-based meteorological stations. This mathematical modeling is based on the maximization of the theoretical expression of the global solar irradiation impinging on an inclined surface, with respect to the slope and orientation of the panel and to the solar hour angle. A set of transcendent equations resulted, whose solutions give the optimum tilt and orientation of a solar panel. A simulation was carried out using global horizontal solar radiation data from the European Solar Radiation Atlas and some empirical models of diffuse solar radiation. The optimum tilt angle resulted was related to latitude by a linear regression with significant correlation coefficients. The standard error of the mean values resulted increased significantly with latitude, suggesting that unreliable values can be provided at high latitudes.


---

## Body

## 1. Introduction

Most countries in the world have realized the need for reduction of gases emission to contrast the adverse global climatic change, encouraging the use of renewable and sustainable sources of energy. Indeed, large quantities of carbon dioxide, nitrogen, and sulfur oxides are emitted in the world by conventional energy sources, which are released to the earth’s atmosphere contributing to climate change.Furthermore, the world will soon run out of its conventional energy resources because of the rapid depletion of fossil fuel reserves. This future scenario and the risks associated with CO2 emissions and global warming have increased the interest in renewable energy.The major renewable energy systems include photovoltaics (PVs), solar thermal, wind, biomass, hydroelectric, and geothermal. However, among various renewable energy sources, the photovoltaic technology for power generation is considered well-suited technology, particularly for distributed power generation.Solar panel is the energy conversion fundamental component of PV systems or solar collectors. Solar panels use light energy from the sun to generate electricity through the photovoltaic effect, whereas solar thermal systems generate heat.The amount of electrical power produced from PV systems is related to the amount of solar irradiation projecting on the modules. Hence, the global solar irradiation on tilted surfaces facing in different directions should be considered to estimate thermal and electrical power obtained in architectural planning.The literature provides that solar power supplied by the modules depends on many extrinsic factors, such as insolation levels, temperature, load conditions, and orientation of the panel.The solar radiation is also a function of the nature and extent of cloud cover and of the atmosphere’s water vapour content, because solar radiation entering Earth’s atmosphere is scattered by atmospheric gases, aerosols, and clouds.Indeed, meteorological parameters used as predictors include the amount and distribution of clouds or other observations such as the fractional sunshine and water content [1, 2].Aerosols can either absorb or scatter the radiation and alter the energy balance of Earth, especially under clear skies [3].These parameters, obviously, cannot be modified, whereas other variables can be changed to maximize the solar energy acquired by the panel. In fact, the design of a solar energy module involves complex tradeoffs due to the interaction of several factors such as the characteristics of the solar cells, power supply requirements and power management features of the embedded system, application behavior, inclination, and orientation of the panel.Hence, it is essential to understand and exploit these factors in order to maximize the energy efficiency of a solar module.In particular, for fixed absorber surfaces, solar energy gain is strictly related to the slope and azimuth angles of a solar panel.The global solar radiation for inclined surfaces can be calculated by the values of direct and diffuse solar radiation on the corresponding horizontal surface. Meteorological data from all parts of the world are needed to know the horizontal global solar radiation, and, for some regions, measured data may only be applied within a radius of about 50 km from weather stations. This circumstance leads to interpolate parameters between stations.Furthermore, spectral irradiance is usually not measured routinely, so that the energy instantaneous production of a solar plant is to rely on appropriate models.Otherwise, the users need accurate computations of the slope and orientation of solar panels in order to maximize the solar energy that is collected by fixed solar panels.Advantages of fixed solar panels are mainly related to their tolerance to misalignment, as approximately 20% of the incident solar radiation is diffuse light, available at any angle of misalignment with the direct sun.In contrast, the main advantage of tracking systems is to collect solar energy for the longest period of the day with the most accurate alignment as sun’s position shifts with the seasons. Indeed, daily solar energy collected was calculated to be 19%–24% higher by a solar PV panel with one axis east-west tracking system than by a fixed system [4].Nevertheless, since solar tracking systems have high operation and maintenance costs and are not always applicable, it is often convenient to set the solar collector at a fixed value of an optimum tilt angle [5].However, some algorithms for the minimization of the energy loss generated by the driving actuator were recently proposed [6, 7].Otherwise, the increase of diffuse solar radiation at low latitudes at some locations with respect to high latitudes makes fixed solar panels a competitive alternative to other energy sources mainly in those locations [8].Hence, the following main question regarding the design of fixed solar panels arises.What tilt and orientation of a solar panel have to be chosen to maximize the solar radiation?Solar radiation impinging on an inclined surface can be divided into direct, diffuse, and ground reflected radiation. Hence, a method to determine these quantities for each latitude has to be studied to answer to that fundamental question.The diffuse solar component is the most difficult quantity to determine, because the distribution of the sky-diffuse radiance strictly depends on the local condition of the sky.One way to estimate the diffuse solar component is to study the regression between the global and the diffuse solar radiation at locations where appropriate data are available, establishing models which may be used to predict the diffuse solar radiation.Liu and Jordan performed the first studies on this subject determining a relationship between daily diffuse and global radiation on a horizontal surface, assuming an isotropic diffusion of solar radiance in the sky [9].Erbs et al. used a database acquired from four US weather stations, composed of hourly direct normal radiation and global radiation, to develop an estimation of the diffuse fraction of hourly, daily, and monthly average global radiation [10].Moreover, [11, 12] used also two predictors for their correlations: the clearness index and the solar elevation.Garrison proposed a model to represent the dependency of the diffuse fraction on the surface albedo, atmospheric precipitable water, atmospheric turbidity, solar elevation, and global horizontal radiation [13].Reindl et al. considered two more significant predictors, the ambient temperature and relative humidity [14], reducing the standard error of Liu and Jordan-type models.On the other side, [15] developed an exponential model for the estimation of the direct normal beam radiation from the global radiation, named the Disc model. This model was then improved by [16].However, accurate mathematical modeling of global and diffuse solar radiation that was used for the simulation in this study is proposed in the following section.Various optimum tilt angle values were provided in the literature for fixed solar panels.For instance, Qiu and Riffat suggested the tilt angle of the solar collector set within the optimum tilt angle of ±10° as an acceptable practice [17].Other, computations led to the values ofβ=φ±20° [18], β=φ±8° [19], β=φ±5° [20], and β=φ±15° [21], where φ represents the geographical latitude and the signs + and − refer to winter and summer months, respectively.The disagreement among these values may be due to two main reasons:(1)
firstly, the different methods of calculation that were used for the determination of the optimum slope value of a solar panel;(2)
secondly, the different empirical models that were considered for the determination of diffuse solar radiance and its link with the amount of global solar radiation.The aim of this study was to propose an algorithm for the determination of the optimum tilt and orientation of a solar panel using a mathematical model based on the orientation of a generic surface with respect to the position of the sun in the sky. The other physical variables and the empirical model of diffuse solar radiation are considered successively in the algorithm, linking it to a specific location.


## 2. Mathematical Modeling of Global and Diffuse Solar Radiation

The literature provides several methods to estimate the global and diffuse solar radiation using climatologic parameters.A typical empirical model is the regression equation of the Angstrom type [22]
(1)HgHo=a+b(nN),
where Hg is the monthly average of daily global solar radiation impinging on a horizontal surface at a location, Ho is the monthly mean of daily radiation on a horizontal surface in the absence of atmosphere, n is the monthly mean of daily number hour of observed sunshine, N is the monthly mean value of day length of the interested location, and “a” and “b” are the regression constants, determined from climatological data. The ratio n/N is often called the “possible sunshine hour percentage.”Regression coefficients “a” and “b” can be obtained from some relationships as proposed by [23, 24].A statistically worldwide 5.3% average decrease of global solar radiation, with largest decline between 45° and 30°N, was found by analyzing the data collected in 45 actinometric stations during the years 1958, 1965, 1975, and 1985 [25]. A significant decrease in mean yearly global solar radiation between the years 1964 and 1990 under completely overcast skies was found in some locations in Germany [26].The observed changes in clear-sky radiation could be related to the recovery by the volcanic eruptions effects, submicron aerosol particles with simultaneous reduction of aerosol mass concentration, and increasing absorption by urban aerosol.Hence, global horizontal solar irradiation data should be updated and have to be preferred with respect to empirical expression such as (1).Otherwise, daily extraterrestrial radiation on a horizontal surface, namedHo, can be computed for the day n from the following equation [27]:
(2)Ho=86400*Gscπ×(1+0.033cos(2πn365))cosφcosδ*(±(1-tan2φtan2δ))+cos-1(-tanφtanδ)*sinφsinδ,
where Gsc is the solar constant (1367 W/m2), n is the number of the day, δ is the solar declination, and φ is the geographical latitude. The hour angle of sunrise ωs has been expressed as a function of φ and δ by Cooper’s equation.The ratio of solar radiation at the surface of the Earth to extraterrestrial radiation is called the clearness index, with the monthly average clearness index,KT, defined as
(3)KT=H_gHo,
where H_g is the monthly average of daily solar radiation on a horizontal surface.The global horizontal solar irradiationHg, provided from meteorological stations, includes the horizontal direct beam irradiation Hb and the horizontal diffuse sky irradiation Hd.Some models can convert global horizontal irradiation to direct beam irradiation and diffuse sky irradiation on the horizontal plane by means of empirical relationships.Heliosat is an algorithm which was developed to estimate ground level global horizontal irradiance by using Meteosat satellites images taken in the visible band.Some predictive models for estimating global solar irradiation have to be used to obtain the solar radiation on a tilted surface, so that the relationship between the global solar irradiation on horizontal planes and that on tilted planes can be evaluated [9, 28–31].Unfortunately, no theoretic relationship between the horizontal sky diffuse irradiationHd and the horizontal global solar irradiation Hg can be determined rigorously. Indeed, a double integral like equation (5) of   [32] should be solved, which cannot be computed even in the simplest case of uniform and isotropic sky diffuse solar radiation, because the distribution of solar irradiance through the sky is difficult to be represented adequately.The direct and diffuse components of solar radiation can be estimated using empirical relationshipsHd=f(Hg) by means of the clearness index KT.Some of these algorithms, requiring the direct normal and diffuse radiation on a horizontal surface as input, can provide very different estimated results in different locations [33].Indeed, the most considerable cause of error in the computation of the optimum tilt and orientation of a solar panel depends on the model of diffuse solar radiation which is used.Muneer [34] recommends the model proposed by [35] for the desert and tropical locations:
(4)HdHg=1.35-1.61KT.
For temperate climates and locations out with the tropics, equation (5) given by [36] may be used:
(5)HdHg=1.00-1.13KT.
The daily horizontal diffuse irradiation Hd and direct beam irradiation Hb=Hg-Hd can be obtained by (4) and (5), as well.


## 3. The Algorithm to Maximize the Global Solar Radiation on a Fixed Sloped Surface

### 3.1. The Case of Isotropic Diffusion

Global solar radiation on a tilted surfaceIT consists of daily direct solar radiation Ib, diffuse solar radiation Id, and ground reflected radiation Ir.Daily solar radiation on a tilted surface for a given month can be estimated as follows [9]:
(6)IT=Ib+Id+Ir.
The daily direct radiation on a tilted surface Ib can be obtained by means of Rb, the ratio of the average daily direct radiation on a tilted plane to that on a horizontal plane and the parameters to it correlated [37]:
(7)Ib=cosθcosθz=HbRb,(8)Rb=cosβ-sinβ(sinδcosφcosγ-sinβ-cosδsinφcosωcosγ-cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1,
where β is the slope of the panel as to the horizontal plane, γ is the azimuth, ω is the solar hour angle, φ is the latitude, and θ and θz are the solar incidence angle on the considered plane and the solar zenith angle, respectively.If we consider a uniform and isotropic distribution of diffuse solar radiation over the sky hemisphere,Id would be easily obtained from the simple approximation of [9]:
(9)Id(iso)=Hd(1+cosβ)2.
The evaluation of the ground-reflected diffuse radiation depends on Ir. Most studies consider that the ground reflection process is ideally isotropic, in which a specific case Ir can be simplified as follows:
(10)Ir=Hgρ(1-cosβ)2,
where ρ represents the diffuse reflectance of the ground (also called ground albedo).Finally, the daily global solar irradiation on slopesIT can be expressed as the sum of (7), (9), and (10).The global solar radiation incident on a sloped surface depends on the position of the sun along its daily trajectory (represented by the solar angleω) and on the orientation of the panel (represented by the slope β and the azimuth γ).The other quantities appearing in (6), (7), (8), (9), and (10) depend on the local conditions and can be considered fixed parameters.The algorithm proposed here is based on the assumption that the daily solar irradiation impinging on a collecting surface is maximum with respect to the anglesβ, γ, and ω, where β is the slope of the panel as to the horizontal plane, γ is the azimuth, and ω is the solar hour angle, respectively.This physical condition for the maximization of the solar radiation acquired by a solar panel can be represented by the mathematical expressions [8]:
(11)∂IT∂β=0,∂IT∂γ=0,∂IT∂ω=0.
Their application to the isotropic model of [9] provides the following expressions [8]:
(12)∂IT∂β=sinβ(Ho+Hd2-Hgρ2)-cosβHo×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1=0,∂IT∂γ=sinγ(cosδsinφcosω-sinδcosφ)-cosγcosδsinω=0,∂IT∂ω=sinωsinδcosγ-cosδcosφsinγ-cosωsinδsinφsinγ=0.
This set of equations can be solved with respect to the angles β, γ, and ω [8]:
(13)β=tan-1{(sinδsinφ+cosδcosφcosω)(Ho+Hd2-Hgρ2)]-1Ho×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)×[(Ho+Hd2-Hgρ2)×(sinδsinφ+cosδcosφcosω)(Ho+Hd2-Hgρ2)]-1},(14)γ=tan-1{cosδsinω(cosδsinφcosω-sinδcosφ)},(15)ω=sin-1{(sin2δcos2γ+sin2δsin2φsin2γ)-1(sinδcosδsinγcosγcosφ±[sin4δ(sin2φsin2γcos2γ+sin4φsin4γ)-sin2δcos2δsin2φcos2φsin4γ]1/2)×(sin2δcos2γ+sin2δsin2φsin2γ)-1}.
The transcendent equations (13), (14), and (15) can be solved by iterative methods with respect to the angles β, γ, and ω, which provide, respectively, the optimum tilt and orientation of the solar panel and the angular position of the sun in the sky where the maximization of the solar radiation on the panel occurs.


### 3.2. The Case of Anisotropic Diffusion

The diffuse solar radiance influences the performance of most of the solar energy technologies. In fact, the scattering atmospheric processes redistribute solar energy out of the direct beam into the diffuse radiation. The PV systems, like flat-plate collectors, have the peculiarity to use both direct and diffuse forms of radiation, whereas solar concentrator systems can use only direct radiation.Unfortunately, the planetary network of diffuse solar radiation measurements stations is poor, whereas global solar radiation data are available for many locations.Nevertheless, empirical correlations between diffuse and global ratios and diffuse and direct solar radiation can be used to reduce errors in computing hemispherical radiation from estimates of direct and diffuse radiation, especially regarding nonclear sky conditions.Besides the isotropic model of diffuse solar radiation described in [9], several models have been proposed to represent the anisotropy of the diffuse component Id by means of empirical relations which should modify the expressions (12) and the relative solutions (13), (14), and (15).The assumption that the diffuse radiation originates entirely from the solar disk gives the relation(16)Id(disk)=Hdcosθcosθz,
which is the modeling opposite to the isotropic one of [9].Le Quere [38] and Hay and Davies [39] proposed a model for the diffuse solar component as a combination of the two components, that is, (9) and (16),
(17)Id=F*Id(iso)+(1-F)*Id(disk),
where (1-F) expresses the anisotropy degree.The valueF=0.8 was suggested by [38], while [39] assumed the ratio of terrestrial direct radiation to extraterrestrial radiation as the degree of anisotropy (1-F).Klucher [40] proposed a model in which the isotropic component Id (iso) is multiplied by two factors that represent both circumsolar and horizon brightening:
(18)Id=Id(iso)(1+Ksin3(β2))(1+Kcos2θsin3θz),
where the parameter K expresses the degree of anisotropy as a modulating function of the amount of direct radiation received by the surface, K=1-Hd/Hg (this model reduces to the Liu-Jordan isotropic model if the ratio of diffuse to global radiation Hd/Hg is close to the unity).Further models have been proposed in the literature to perform the anisotropy of diffuse solar irradiance, by means of some coefficients derived from statistical analyses of empirical data for specific locations [41–43].Nevertheless, the anisotropic models of [38–40] have been taken into account here, because of their easy applicability to every location. In addition, they are extensively validated models that convert hemispherical data on a horizontal surface to hemispherical data on a tilted surface, computing diffuse solar radiation.The mathematical conditions (11) applied to the models of [38, 39] provided the following solutions:
(19)β=tan-1{[(Ho+FHd2+(1-F)Hd-Hgρ2)(Ho+(1-F)Hd)×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)(Ho+FHd2+(1-F)Hd-Hgρ2)]×[(Ho+FHd2+(1-F)Hd-Hgρ2)×(sinδsinφ+cosδcosφcosω)(Ho+FHd2+(1-F)Hd-Hgρ2)]-1},(20)γ=tan-1{[(cosδsinφcosω-sinδcosφ)-1(Ho+(1-F)Hd)cosδsinω]×[(Ho+(1-F)Hd)×(cosδsinφcosω-sinδcosφ)]-1},(21)ω=sin-1{[×sin2φcos2φsin4γ]1/2(Ho+(1-F)Hd)sinδcosδsinγcosγcosφ±[[(Ho+(1-F)Hd)sinδ]4*(sin2φsin2γcos2γ+sin4φsin4γ)-[(Ho+(1-F)Hd)sinδcosδ]2×sin2φcos2φsin4γ]1/2]×[[(Ho+(1-F)Hd)sinδ]2×(cos2γ+sin2φsin2γ)]-1×sin2φcos2φsin4γ]1/2}.
Equations (19), (20), and (21) can be solved by iterative methods with respect to the angles β, γ, and ω. The mathematical conditions (11) applied to the model of [40], represented by (18), provided the following expressions: (22)∂IT∂β=-Hosinβ-Hocosβ[×(sinδsinφ+cosδcosφcosω)-1(sinδcosφcosγ-cosδsinφcosωcosγ-cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1]-Hdsinβ(1+Fsin3(β/2))2+3sin2(β2)HdF(1+cosβ)2-HdFsinβq1q2q32+3Fsin2(β2)q1q2HdF(1+cosβ)2+Hgρsinβ2+HdF(1+cosβ)(q1)*[cosβ(cosδsinφcosωcosγ+cosδsinγsinω-sinδcosφcosγ)-sinβ(sinδsinφ+cosδcosφcosω)]q2q3=0,(23)∂IT∂γ=sinβ(Hosinδcosφsinγ-Hocosδsinφcosωsinγ+Hocosδcosγsinω)×(sinδsinφ+cosδcosφcosω)-1+(1+Fsin3(β2))HdF(1+cosβ)q2(q1)×[sinβ(sinδcosφsinγ+cosδcosγsinω-cosδsinφcosωsinγ)(]=0,(24)∂IT∂ω=Hosinβ(cosδcosωsinγ-cosδsinφsinωcosγ)sinδsinφ+cosδcosφcosω+(1+Fsin3(β2))HdF(1+cosβ)q2q1×[sinβ(cosωcosδsinγ-sinωcosδsinφcosγ)-cosβcosδcosφsinω]+(1+Fsin3(β2))HdF(1+cosβ)(34)q1×[1-(sinδsinφ+cosδcosφcosω)2]1/2×(sinωcosδcosφ)=0,where
(25)q1=[cosβ(sinδsinφ+cosδcosφcosω)+sinβ(cosδsinφcosωcosγ+cosδsinωsinγ-sinδcosφcosγ)]2,q2=[1-(sinδsinφ+cosδcosφcosω)2]3/2,q3=[1+Fsin3(β2)].
This set of equation can be solved only using iterative methods.


## 3.1. The Case of Isotropic Diffusion

Global solar radiation on a tilted surfaceIT consists of daily direct solar radiation Ib, diffuse solar radiation Id, and ground reflected radiation Ir.Daily solar radiation on a tilted surface for a given month can be estimated as follows [9]:
(6)IT=Ib+Id+Ir.
The daily direct radiation on a tilted surface Ib can be obtained by means of Rb, the ratio of the average daily direct radiation on a tilted plane to that on a horizontal plane and the parameters to it correlated [37]:
(7)Ib=cosθcosθz=HbRb,(8)Rb=cosβ-sinβ(sinδcosφcosγ-sinβ-cosδsinφcosωcosγ-cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1,
where β is the slope of the panel as to the horizontal plane, γ is the azimuth, ω is the solar hour angle, φ is the latitude, and θ and θz are the solar incidence angle on the considered plane and the solar zenith angle, respectively.If we consider a uniform and isotropic distribution of diffuse solar radiation over the sky hemisphere,Id would be easily obtained from the simple approximation of [9]:
(9)Id(iso)=Hd(1+cosβ)2.
The evaluation of the ground-reflected diffuse radiation depends on Ir. Most studies consider that the ground reflection process is ideally isotropic, in which a specific case Ir can be simplified as follows:
(10)Ir=Hgρ(1-cosβ)2,
where ρ represents the diffuse reflectance of the ground (also called ground albedo).Finally, the daily global solar irradiation on slopesIT can be expressed as the sum of (7), (9), and (10).The global solar radiation incident on a sloped surface depends on the position of the sun along its daily trajectory (represented by the solar angleω) and on the orientation of the panel (represented by the slope β and the azimuth γ).The other quantities appearing in (6), (7), (8), (9), and (10) depend on the local conditions and can be considered fixed parameters.The algorithm proposed here is based on the assumption that the daily solar irradiation impinging on a collecting surface is maximum with respect to the anglesβ, γ, and ω, where β is the slope of the panel as to the horizontal plane, γ is the azimuth, and ω is the solar hour angle, respectively.This physical condition for the maximization of the solar radiation acquired by a solar panel can be represented by the mathematical expressions [8]:
(11)∂IT∂β=0,∂IT∂γ=0,∂IT∂ω=0.
Their application to the isotropic model of [9] provides the following expressions [8]:
(12)∂IT∂β=sinβ(Ho+Hd2-Hgρ2)-cosβHo×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1=0,∂IT∂γ=sinγ(cosδsinφcosω-sinδcosφ)-cosγcosδsinω=0,∂IT∂ω=sinωsinδcosγ-cosδcosφsinγ-cosωsinδsinφsinγ=0.
This set of equations can be solved with respect to the angles β, γ, and ω [8]:
(13)β=tan-1{(sinδsinφ+cosδcosφcosω)(Ho+Hd2-Hgρ2)]-1Ho×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)×[(Ho+Hd2-Hgρ2)×(sinδsinφ+cosδcosφcosω)(Ho+Hd2-Hgρ2)]-1},(14)γ=tan-1{cosδsinω(cosδsinφcosω-sinδcosφ)},(15)ω=sin-1{(sin2δcos2γ+sin2δsin2φsin2γ)-1(sinδcosδsinγcosγcosφ±[sin4δ(sin2φsin2γcos2γ+sin4φsin4γ)-sin2δcos2δsin2φcos2φsin4γ]1/2)×(sin2δcos2γ+sin2δsin2φsin2γ)-1}.
The transcendent equations (13), (14), and (15) can be solved by iterative methods with respect to the angles β, γ, and ω, which provide, respectively, the optimum tilt and orientation of the solar panel and the angular position of the sun in the sky where the maximization of the solar radiation on the panel occurs.


## 3.2. The Case of Anisotropic Diffusion

The diffuse solar radiance influences the performance of most of the solar energy technologies. In fact, the scattering atmospheric processes redistribute solar energy out of the direct beam into the diffuse radiation. The PV systems, like flat-plate collectors, have the peculiarity to use both direct and diffuse forms of radiation, whereas solar concentrator systems can use only direct radiation.Unfortunately, the planetary network of diffuse solar radiation measurements stations is poor, whereas global solar radiation data are available for many locations.Nevertheless, empirical correlations between diffuse and global ratios and diffuse and direct solar radiation can be used to reduce errors in computing hemispherical radiation from estimates of direct and diffuse radiation, especially regarding nonclear sky conditions.Besides the isotropic model of diffuse solar radiation described in [9], several models have been proposed to represent the anisotropy of the diffuse component Id by means of empirical relations which should modify the expressions (12) and the relative solutions (13), (14), and (15).The assumption that the diffuse radiation originates entirely from the solar disk gives the relation(16)Id(disk)=Hdcosθcosθz,
which is the modeling opposite to the isotropic one of [9].Le Quere [38] and Hay and Davies [39] proposed a model for the diffuse solar component as a combination of the two components, that is, (9) and (16),
(17)Id=F*Id(iso)+(1-F)*Id(disk),
where (1-F) expresses the anisotropy degree.The valueF=0.8 was suggested by [38], while [39] assumed the ratio of terrestrial direct radiation to extraterrestrial radiation as the degree of anisotropy (1-F).Klucher [40] proposed a model in which the isotropic component Id (iso) is multiplied by two factors that represent both circumsolar and horizon brightening:
(18)Id=Id(iso)(1+Ksin3(β2))(1+Kcos2θsin3θz),
where the parameter K expresses the degree of anisotropy as a modulating function of the amount of direct radiation received by the surface, K=1-Hd/Hg (this model reduces to the Liu-Jordan isotropic model if the ratio of diffuse to global radiation Hd/Hg is close to the unity).Further models have been proposed in the literature to perform the anisotropy of diffuse solar irradiance, by means of some coefficients derived from statistical analyses of empirical data for specific locations [41–43].Nevertheless, the anisotropic models of [38–40] have been taken into account here, because of their easy applicability to every location. In addition, they are extensively validated models that convert hemispherical data on a horizontal surface to hemispherical data on a tilted surface, computing diffuse solar radiation.The mathematical conditions (11) applied to the models of [38, 39] provided the following solutions:
(19)β=tan-1{[(Ho+FHd2+(1-F)Hd-Hgρ2)(Ho+(1-F)Hd)×(cosδsinφcosωcosγ-sinδcosφcosγ+cosδsinγsinω)(Ho+FHd2+(1-F)Hd-Hgρ2)]×[(Ho+FHd2+(1-F)Hd-Hgρ2)×(sinδsinφ+cosδcosφcosω)(Ho+FHd2+(1-F)Hd-Hgρ2)]-1},(20)γ=tan-1{[(cosδsinφcosω-sinδcosφ)-1(Ho+(1-F)Hd)cosδsinω]×[(Ho+(1-F)Hd)×(cosδsinφcosω-sinδcosφ)]-1},(21)ω=sin-1{[×sin2φcos2φsin4γ]1/2(Ho+(1-F)Hd)sinδcosδsinγcosγcosφ±[[(Ho+(1-F)Hd)sinδ]4*(sin2φsin2γcos2γ+sin4φsin4γ)-[(Ho+(1-F)Hd)sinδcosδ]2×sin2φcos2φsin4γ]1/2]×[[(Ho+(1-F)Hd)sinδ]2×(cos2γ+sin2φsin2γ)]-1×sin2φcos2φsin4γ]1/2}.
Equations (19), (20), and (21) can be solved by iterative methods with respect to the angles β, γ, and ω. The mathematical conditions (11) applied to the model of [40], represented by (18), provided the following expressions: (22)∂IT∂β=-Hosinβ-Hocosβ[×(sinδsinφ+cosδcosφcosω)-1(sinδcosφcosγ-cosδsinφcosωcosγ-cosδsinγsinω)×(sinδsinφ+cosδcosφcosω)-1]-Hdsinβ(1+Fsin3(β/2))2+3sin2(β2)HdF(1+cosβ)2-HdFsinβq1q2q32+3Fsin2(β2)q1q2HdF(1+cosβ)2+Hgρsinβ2+HdF(1+cosβ)(q1)*[cosβ(cosδsinφcosωcosγ+cosδsinγsinω-sinδcosφcosγ)-sinβ(sinδsinφ+cosδcosφcosω)]q2q3=0,(23)∂IT∂γ=sinβ(Hosinδcosφsinγ-Hocosδsinφcosωsinγ+Hocosδcosγsinω)×(sinδsinφ+cosδcosφcosω)-1+(1+Fsin3(β2))HdF(1+cosβ)q2(q1)×[sinβ(sinδcosφsinγ+cosδcosγsinω-cosδsinφcosωsinγ)(]=0,(24)∂IT∂ω=Hosinβ(cosδcosωsinγ-cosδsinφsinωcosγ)sinδsinφ+cosδcosφcosω+(1+Fsin3(β2))HdF(1+cosβ)q2q1×[sinβ(cosωcosδsinγ-sinωcosδsinφcosγ)-cosβcosδcosφsinω]+(1+Fsin3(β2))HdF(1+cosβ)(34)q1×[1-(sinδsinφ+cosδcosφcosω)2]1/2×(sinωcosδcosφ)=0,where
(25)q1=[cosβ(sinδsinφ+cosδcosφcosω)+sinβ(cosδsinφcosωcosγ+cosδsinωsinγ-sinδcosφcosγ)]2,q2=[1-(sinδsinφ+cosδcosφcosω)2]3/2,q3=[1+Fsin3(β2)].
This set of equation can be solved only using iterative methods.


## 4. Applying the Maximization Algorithm to a Data Set of Global Horizontal Solar Radiation

The data collected at the European Solar Radiation Atlas (available at the Internet site HelioClim) were used to test the algorithm proposed here.In particular, the global horizontal solar radiationHg relative to Tripoli, Athens, Rome, Paris, London, and Stockholm, located from 32° to 59° north latitude, was considered by averaging the values on a monthly basis from 1985 to 1989 and was used for this simulation.The calculations ofβ, γ, and ω were carried out at every month because their time variation cannot be neglected [44].The solar declination valuesδ were computed at the middle of every month by Cooper’s equation, which was also used to calculate the sunset hour angle ωs [37]. Daily extraterrestrial radiation on a horizontal surface, Ho, was computed at the same day by [27]; equation (5) from [39] was used as a relationship between Hd and Hg. Ground albedo was fixed at the typical value ρ=0.20.The isotropic model of [9] and the anisotropic models of [38–40] were used to calculate the diffuse solar component Hd.The solutions of (14), (15), (20), (21), (23), and (24) gave the values γ=0° and ω=0°, confirming that the optimum orientation of a solar panel is toward South, where in fact the solar hour angle at solar noon must be zero, as it is defined.Nevertheless, other models of the solar diffuse component could provide different values.The optimum tilt angle values, provided by the solutions of (13), (19), (22), were reported on Tables 1, 2, 3, 4, 5, and 6.Table 1
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of Tripoli (latitude +32°57′N; longitude +13°12′E).
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
50.7
12.44
55.3
13.05
53.3
12.73
57.0
18.54
1.17
1.26
d
2 = −13.289
43.6
17.13
47.5
17.96
46.0
17.59
48.7
21.56
0.95
0.87
d
3 = −2.819
33.7
22.54
37.0
23.61
35.8
23.18
37.4
24.06
0.72
0.28
d
4 = 9.415
21.9
27.08
24.3
28.24
23.5
27.79
23.8
26.01
0.45
0.42
d
5 = 18.792
12.9
29.26
14.3
30.16
13.8
29.83
13.6
27.56
0.25
0.50
d
6 = 23.314
8.6
30.92
9.5
31.53
9.2
31.34
8.9
29.31
0.17
0.44
d
7 = 21.517
10.6
31.92
11.4
32.55
11.2
32.39
10.9
30.00
0.15
0.50
d
8 = 13.784
18.4
30.71
19.8
31.58
19.4
31.35
19.4
28.77
0.25
0.55
d
9 = 2.217
29.5
25.74
31.9
26.77
31.2
26.43
32.0
25.90
0.50
0.22
d
10 =  −9.599
40.9
20.16
44.0
20.99
43.0
20.69
44.9
23.73
0.75
0.69
d
11 = −19.148
48.5
13.26
53.2
13.94
51.1
13.58
54.8
18.79
1.18
1.13
d
12 = −23.335
51.8
10.64
57.2
11.21
54.5
10.87
58.9
16.54
1.35
1.23Table 2
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of Athens (latitude +37°58′N; longitude +23°43′E).
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
50.1
6.58
60.5
7.42
53.5
6.73
60.2
9.89
2.22
0.66
d
2 = −13.289
42.3
9.42
52.7
10.66
45.8
9.68
52.2
12.16
2.19
0.54
d
3 = −2.819
33.0
13.16
42.4
14.94
36.2
13.59
40.9
14.65
1.87
0.37
d
4 = 9.415
24.7
20.47
30.2
22.59
27.1
21.24
28.8
20.13
1.03
0.47
d
5 = 18.792
17.1
24.75
20.5
26.65
18.7
25.56
18.9
23.23
0.61
0.62
d
6 = 23.314
13.6
28.03
15.8
29.58
14.8
28.82
14.6
26.03
0.39
0.66
d
7 = 21.517
15.7
29.27
17.7
30.82
16.9
30.15
16.8
26.97
0.36
0.73
d
8 = 13.784
22.9
27.03
25.9
28.83
24.7
28.03
25.3
25.33
0.55
0.65
d
9 = 2.217
33.6
21.73
37.8
23.49
36.0
22.62
38.2
22.39
0.90
0.31
d
10 = −9.599
43.2
14.21
49.5
15.64
46.2
14.77
50.6
17.48
1.45
0.62
d
11 = −19.148
48.2
7.38
58.5
8.32
51.6
7.56
58.2
10.63
2.19
0.65
d
12 = −23.335
52.3
5.92
62.5
6.67
55.6
6.06
62.2
9.35
2.18
0.69Table 3
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of Rome (latitude +41°54′N; longitude +12°27′E).
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
57.5
6.93
64.6
7.93
60.4
7.21
65.8
11.47
1.66
0.91
d
2 = −13.289
49.5
10.56
56.9
12.04
52.7
11.02
58.1
14.38
1.71
0.74
d
3 = −2.819
40.1
16.05
46.7
18.15
43.1
16.82
47.4
18.12
1.47
0.44
d
4 = 9.415
28
20.02
34.3
22.67
30.7
20.95
33.3
19.77
1.22
0.56
d
5 = 18.792
20.2
24.07
24.6
26.66
22.2
25.07
23
22.46
0.79
0.76
d
6 = 23.314
16.8
27.14
20
29.45
18.4
28.17
18.5
24.90
0.57
0.83
d
7 = 21.517
19
28.82
21.9
31.08
20.7
29.98
20.9
26.21
0.52
0.90
d
8 = 13.784
26.5
26.82
30.1
29.19
28.6
28.07
29.7
24.99
0.70
0.78
d
9 = 2.217
36.5
19.81
41.8
22.06
39.2
20.80
42.3
20.75
1.15
0.40
d
10 = −9.599
49
15.29
53.6
16.83
51.6
16.04
55.3
19.35
1.17
0.76
d
11 = −19.148
55.7
8.01
62.6
9.13
58.6
8.35
63.9
12.56
1.63
0.90
d
12 = −23.335
61.5
7.05
66.7
7.94
64
7.38
68.4
12.76
1.31
1.15Table 4
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of Paris (latitude +48°51′N; longitude +2°20′E).
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
56.7
3.98
68.1
4.94
59.9
4.10
66.5
6.54
2.33
0.51
d
2 = −13.289
46.9
6.05
60.4
7.53
50.4
6.25
58
8.21
2.74
0.44
d
3 = −2.819
38.3
10.30
50.2
12.66
41.7
10.73
48.3
11.94
2.41
0.47
d
4 = 9.415
28.7
15.74
38.1
18.95
31.7
16.47
35.7
15.92
1.80
0.64
d
5 = 18.792
20.4
17.66
28.4
21.00
22.8
18.36
24.5
16.94
1.45
0.76
d
6 = 23.314
18.3
22.16
23.9
25.46
20.4
23.11
21.2
20.59
1.00
0.88
d
7 = 21.517
21.1
24.83
25.9
28.13
23.3
26.05
24.2
22.76
0.87
0.97
d
8 = 13.784
27.8
21.73
33.9
25.06
30.5
22.92
32.9
20.72
1.17
0.80
d
9 = 2.217
39
17.48
45.7
20.24
42
18.50
46.2
18.67
1.46
0.49
d
10 = −9.599
48.7
10.22
57.1
12.14
52
10.75
57.9
13.32
1.88
0.60
d
11 = −19.148
58.1
5.93
66.3
7.12
61.1
6.21
66.9
9.67
1.83
0.73
d
12 = −23.335
62.6
4.31
70.3
5.22
65.3
4.50
70.6
8.02
1.70
0.74Table 5
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of London (latitude +51°32′N; longitude 0°05′W).
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
61.8
2.24
73.8
3.10
64.7
2.32
70.6
4.00
2.36
0.36
d
2 = −13.289
52.1
4.39
66.1
5.95
55.5
4.56
63
6.26
2.80
0.41
d
3 = −2.819
41.4
7.79
55.9
10.46
44.9
8.14
52.6
9.26
2.90
0.52
d
4 = 9.415
32
13.31
44
17.36
35.3
14.02
40.8
13.70
2.33
0.80
d
5 = 18.792
25.8
18.88
34.6
23.65
28.6
19.96
31.7
17.86
1.45
1.09
d
6 = 23.314
20.7
17.87
29.9
22.55
23.1
18.71
25.2
16.84
1.69
1.07
d
7 = 21.517
21.7
16.61
31.7
21.20
24.2
17.37
26.7
15.83
1.84
1.03
d
8 = 13.784
28.2
14.54
39.6
18.85
31.2
15.29
35.8
14.44
2.17
0.90
d
9 = 2.217
39.1
11.02
51.1
14.37
42.6
11.63
49.3
12.24
2.43
1.13
d
10 = −9.599
54.9
8.65
62.9
10.84
58
9.23
63.8
11.86
1.82
1.13
d
11 = −19.148
63.1
3.68
72.0
4.84
65.9
3.87
71.3
6.53
1.85
1.06
d
12 = −23.335
64.2
1.73
75.8
2.41
66.8
1.78
72.2
3.31
2.26
0.32Table 6
Optimum tilt angle values and monthly mean of the daily global solar radiation values computed by means of the algorithm proposed in Section2 applied to the monthly averages of the daily global horizontal solar irradiation data acquired from 1985 to 1989 at the site of Stockholm (latitude +59°17′N; longitude +18°03′E). The monthly mean of the daily global solar radiation was not acquired in the periods relative to the solar declinations d1, d11, and d12.
Solar declin. (deg)
Isotropic model
Le Quere’s model
Hay and Davies’ model
Klucher’s model
SEM values
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
β (deg)
E (MJ/m2)
d
1 = −21.269
n.c.
n.c.
n.c.
n.c.
n.c.
n.c.
d
2 = −13.289
67.1
5.18
73.2
6.99
69.5
5.65
73.8
8.79
1.37
0.70
d
3 = −2.819
54.7
10.25
63.1
13.77
57.9
11.15
63.8
12.56
1.88
0.68
d
4 = 9.415
43.1
18.64
51.3
24.32
46.4
20.37
51.8
18.32
1.80
1.19
d
5 = 18.792
33.0
21.99
41.8
29.10
36.2
23.88
40.7
19.95
1.77
1.70
d
6 = 23.314
29.5
25.39
37.3
32.92
32.4
27.54
35.9
22.19
1.52
1.95
d
7 = 21.517
29.8
21.84
39.0
29.12
32.9
23.59
37.0
19.68
1.78
1.75
d
8 = 13.784
33.8
14.17
46.5
20.12
37.2
15.15
43.4
13.98
2.50
1.25
d
9 = 2.217
46.3
10.25
58.0
14.45
49.8
11.02
56.9
11.52
2.44
0.79
d
10 = −9.599
60.5
5.65
69.5
7.93
63.5
6.09
69.3
8.25
1.92
0.56
d
11 = −19.148
n.c.
n.c.
n.c.
n.c.
n.c.
n.c.
d
12 = −23.335
n.c.
n.c.
n.c.
n.c.
n.c.
n.c.The corresponding monthly average of the daily global solar radiation (measured in MJ/m2) collected by a surface inclined at each β value was indicated in Tables 1–6, as well. Such a value can be considered representative of the conversion efficiency of a PV module, because of a linear relationship with solar radiation rate [45].The results led to the following considerations.(1)
The optimum tilt angle values of solar panels determined for winter months were confirmed to be very different from the values recommended for summer months.(2)
The valueγ=0°  was confirmed to be the optimum orientation value for a solar panel.(3)
The disagreement among the models of diffuse solar irradiance resulted was relevant for the winter months (represented by the solar declinationsd1–3 and d10–12, corresponding to the first three and the last three rows of Tables 1–6) and resulted in increase with increasing latitude.Furthermore, the optimum tilt angle values were averaged for each monthly solar declination and were reported as a function of latitude in Table7. The average tilt angles resulted were significantly related to latitudes. A correlation between these two sets of variables was studied, applying a linear relationship to predict the value of the optimum tilt angle for a given geographical latitude.Table 7
The averageβ of the optimum tilt angle values computed using the diffuse solar radiation models, relative to the latitude values used above. A linear regression with 95% confidence interval provided the coefficients a1 and a2 and the correlation coefficients ranging from 0.944 to 0.993, reported in the last column.
β (l1)
β (l2)
β (l3)
β (l4)
β (l5)
β (l6)
Coefficients
Correlationr
a
1
a
2
54.07
56.08
62.07
62.80
67.72
—
31.33
0.68
0.953
46.45
48.25
54.30
53.93
59.17
70.90
16.25
0.86
0.944
35.97
38.13
44.32
44.62
48.70
59.87
6.80
0.84
0.957
23.37
27.70
31.57
33.55
38.03
48.15
-
6.07
0.87
0.975
13.65
18.80
22.50
24.03
30.17
37.92
-
14.95
0.87
0.978
9.05
14.70
18.42
20.95
24.73
33.77
-
19.27
0.87
0.985
11.03
16.78
20.62
23.62
26.08
34.67
-
15.65
0.83
0.988
19.25
24.70
28.72
31.27
33.70
40.22
-
4.23
0.75
0.989
31.15
36.40
39.95
43.22
45.53
52.75
6.42
0.77
0.993
43.20
47.37
52.37
53.93
59.90
65.70
15.84
0.83
0.983
51.90
54.13
60.20
63.10
68.07
—
23.61
0.84
0.977
55.60
58.15
65.15
67.20
69.75
—
30.56
0.76
0.968In particular, the regression line was performed by the method of least squares, to make the sum of the squares of the differences between the ordinates of the points and those on the straight line as small as possible. The correlation coefficient “r” of the variables was calculated, defined as their covariance (that measures how the two variables are linearly related) divided by the product of their individual standard deviations. In addition, a confidence interval was considered to know how accurate is the regression.Linear regressions with 95% confidence interval provided correlation coefficients ranging from 0.944 to 0.993, reported in the last column of Table7 (larger correlation values were found for summer months). Two typical fits for two solar declinations are shown in Figures 1(a) and 1(b).Average values of optimum tilt angles, computed using diffuse solar radiation models, as a function of latitude, relative to the solar declinationd9 (a) and to the solar declination d3 (b). Linear regressions with 95% confidence interval provided the correlation coefficients r=0.993 and r=0.957, respectively.
(a)
(b)This result shows that the optimum tilt angle of a solar panel could be easily obtained as a function of geographical latitude by means of the coefficientsa1 and a2 of a linear regression β=a1+a2*φ. The coefficients a1 and a2, obtained by the regression, depend on the solar declination that was used.Furthermore, the application of the algorithm showed that the differences among the tilt angle values, computed using various models of diffuse solar irradiance, increase with increasing of the geographical latitude, suggesting that further empirical correlations between diffuse and global ratios and diffuse and direct solar radiation should be investigated, especially at high latitudes.Finally, the standard error of the mean (SEM) of the tilt angles and of the monthly average of the daily global solar radiation values was computed and reported in the last column of Tables1–6, to quantify the disagreement among the various models of diffuse solar irradiance used here.The average values of these SEM values were computed separately for summer months and for winter months and reported in Table8.Table 8
The average standard error of the mean (SEM) of the tilt angles and of the monthly average of the daily global solar radiation values of diffuse solar radiation, computed separately for summer months (solar declinationsd4–d9) and for winter months (solar declinations d1–d3 and d10–d12) as a function of latitude (The correlation coefficient r relative to winter months was computed excluding the values relative to the latitude of Stockholm).
Latitude (deg)
Average values of SEM values
Summer months
Winter months
β
E
β
E
32.95
0.30
0.44
1.01
0.91
37.96
0.64
0.57
2.01
0.59
41.90
0.82
0.71
1.50
0.82
48.85
1.26
0.75
2.15
0.58
51.53
2.02
0.92
2.33
0.46
59.28
1.96
1.44
n.c.
n.c.
Correlation coefficientr
0.956
0.935
0.835
−0.789A significant correlation between the averages of these SEM values and the latitudes was found for summer months. In fact, a linear regression with 95% confidence interval provided the correlation coefficientsr=0.956 and r=0.935 for the tilt angles and the monthly average of the daily global solar radiation values, respectively. These linear regressions were represented in Figures 2(a) and 2(b).The standard error of the mean (SEM) of the tilt angle values (a) and of the monthly average of the daily global solar radiation values (b), computed for the summer months (solar declinationsd4–d9) as a function of latitude. The linear fits provided the correlation coefficients r=0.956 and r=0.935, respectively.
(a)
(b)Regarding the winter months, lower correlations were found (see Table8). Nevertheless, the values obtained for the winter months can be considered less reliable because the decrease in intensity of direct and diffuse solar radiation during winter months produces an increase in the relative error for estimating the tilt angles and the monthly average of the daily global solar radiation values.


## 5. Conclusions

The method of maximization of global solar radiation proposed here led to take into account a set of equations, which can be solved with respect to the slopeβ and the azimuth γ of a collecting surface and to the solar hour angle ω. The proposed algorithm was applied to some models of diffuse solar irradiance, using global horizontal solar radiation data from the European Solar Radiation Atlas. The relative equations provided a set of angles β, γ, and ω, as a function of the other physical parameters. The better solar panel orientation toward South and the large time dependence of the optimum tilt angle were confirmed, suggesting the use of semifixed panels on building structures.The solutions of the set of equations provided tilt angles and monthly average of daily global solar radiation values strictly related to the model of diffuse solar irradiation which was used.The computed tilt angles resulted were significantly related to the latitude value, because linear regressions with 95% confidence interval provided correlation coefficients ranging from 0.944 to 0.993. This result allows to estimate the optimum tilt angle of a solar panel as a function of latitude by means of a linear regression, whose coefficients depend on global horizontal solar radiation data.Nevertheless, relevant differences among the monthly average of the daily global solar radiation​​and the tilt angle values, obtained using different models of diffuse solar radiation, were highlighted by the simulations performed separately for summer and winter months.The standard error of the mean of such values resulted in significant increase with increasing latitude for summer months, leading to the conclusion that values unreliable for the optimum tilt angle of a solar panel can be provided by the models of diffuse solar radiation at high latitudes.Hence, further research is needed to better estimate the diffuse solar radiation.



---

*Source: 307547-2013-04-09.xml*