Abstract

Incident solar radiation on photovoltaic (PV) solar panels is not constant throughout the year. Besides dependence on the season, solar radiation is reliant on the location and weather conditions. For a given location on Earth, the best-fixed orientation of a PV panel can be determined by achieving the maximum incident solar irradiance throughout the year or for a predetermined period. In this paper, we use a sophisticated atmospheric radiative transfer model to calculate the direct and diffuse solar irradiation (radiant exposure) for the solar spectrum incident on PV solar panels to determine the best tilt angle of the panel in order to maximize absorption of solar radiation for selected periods. We used the Regula-Falsi numerical method to obtain the tilt angle at which the derivative of solar irradiation (concerning the tilt angle) approaches zero. Moreover, the spectral response of typical silicon cells is taken into account. These calculations were carried out in São Carlos (SP), a town in the southeast of Brazil. The best tilt angle was obtained for three selected periods. Additionally, we provide results for Southern latitudes ranging from 0° to −55° in steps of −5° for the meteorological seasons. We have shown that for each period, there is an increase in solar radiation absorption compared to the traditional installation angle based exclusively on the local latitude. These calculations can be extended to any location.

1. Introduction

Electricity, as a type of energy, has become historically the cheapest and most widely used and has emerged as an essential and strategic resource for humanity. Nevertheless, electricity has not yet reached the entire world population, often depriving communities of the benefits this valuable resource can provide. In Brazil, there is a predominance of hydroelectric generation [1], comprising 64.54%, while the photovoltaic (PV) generation is only 0.02%. This alternative power generation is often presented as the only viable solution in remote regions or distributed forms of power generation and small installations.

Photovoltaic solar tracking systems aim to maximize exposure of the panels to sunlight keeping them as close as possible in the orthogonal position to the incident rays. However, in some situations such as in remote locations or rural areas, tracking systems may be impractical and/or expensive and therefore are not always applicable [2]. Thus, in many applications, fixed solar panels are the best option. Tracking systems may not fully take into account the incidence of radiation scattered by the atmosphere and the surroundings. Certainly, the direct normal irradiance from the sun is a major component of global irradiance incident on the panels. However, on cloudy days, solar panels still capture the sun’s scattered radiation, indicating an important contribution. Figure 1 schematically shows the components of the total irradiance on a panel that is tilted at an angle to the horizontal.

There are various experimental and theoretical studies in the literature that determine ideal tilt angles of PV panels based on many assumptions to account for characteristics of the installation site [312]. An extensive literature review on the optimization of the tilt angle to maximize incident solar radiation has been provided elsewhere by Yadav and Chandel [13]. As pointed out by these authors [13], in most studies, direct and ground-reflected components of global irradiance were computed by using simplified models [13]. The diffuse component, by its nature, is more complicated and needs more sophisticated models [1416]. Moreover, most studies in the literature [13] are for Northern hemisphere sites and use measured solar irradiation and/or use hybrid experimental/theoretical methods. Few works in the literature [13] rely exclusively on theoretical methods such as the present study. Thus, this research is useful for sites where measured data are unavailable or unreliable. To the best of the authors’ knowledge, the literature is very scarce, if any, concerning the use of spectral responses of silicon cells in weighting solar irradiance to improve calculations. Indeed, none of the references cited in the aforementioned review [13] take into account the spectral responsivity of solar cells. Therefore, our work is situated in this area where there are limited available resources.

In view of this rationale, this paper is aimed at determining an optimum tilt angle for solar PV arrays in order to maximize incident solar irradiance on the array for specific periods, considering all solar irradiance components, i.e., direct normal, diffuse, and ground-reflected ones. In addition, the solar irradiance incident into the PV solar panel is weighted by the spectral response of typical silicon cells, strengthening the relevance of present work. We also use a very efficient and robust, yet simple, numerical method to determine the best tilt angle. This combination of methods is unique in the field. Furthermore, this is a strictly theoretical approach that does not rely on collected data, turning our approach universal.

2. Methods

In this section, the calculation and computer simulations will be detailed as follows.

2.1. Theoretical Foundations

A typical situation that this work aims to address is that of a flat solar panel, located in the Southern hemisphere, with a tilt angle to the horizontal and directed to the geographic North. In addition, the solar panel is exposed to atmospheric conditions described as clear sky or without clouds. Thus, the calculation of spectral irradiance and solar irradiation at the tilted surface is required. For this purpose, the solar spectrum needs to be simulated using an atmospheric radiative transfer model [1416].

The global solar irradiance (W/m2) for a given time and place is calculated in terms of the spectral irradiance:

In this expression, is the solar spectral irradiance as a function of time and vector represents the location’s spatial coordinates, given by geographic coordinates (latitude, longitude, and altitude). In Equation (1) and onwards, refers to the wavelength. Figure 2 shows the angles involved in the calculation of solar radiation on a tilted surface.

The direct normal irradiance, denoted by , is obtained by: in which the angular function is the angle between the sun’s incident rays and the panel’s surface (see Figure 2) as a function of the location, the day of the year, and the selected time of day. The irradiance is given by Equation (1). The diffuse () and reflected () components are considered isotropic. To calculate the global solar irradiation (J/m2) on a tilted surface over a period, the three components of irradiance should be added, i.e.,

The integrals in Equations (1) and (3) are numerically calculated by a five-point Gauss-Legendre quadrature method. This scheme is very robust and provides accurate results by knowing only five points of the integrand. In order to increase precision, the integration intervals are subdivided and the integrals are calculated in each subinterval. The results of each subinterval are finally added up.

The challenge for these calculations is to obtain the spectral irradiance at the panel surface in the desired location for the day of the year and time of day. Furthermore, typical atmospheric data such as ozone concentration, water vapor column, and other information may be considered.

Figure 3 shows typical global solar spectral irradiance at the top of the atmosphere (top curve) and at sea level (bottom curve). Note the attenuation resulting from absorption bands of gases and water vapor at sea level.

Experimental data on solar irradiance are not available for all locations. Thus, irradiance must be obtained using an atmospheric model. Some solar spectrum modeling at the Earth’s surface exists in the literature for semiquantitative applications [17, 18]. In this paper, we propose using a well-known model in the literature, called the Simple Model of the Atmospheric Radiative Transfer of Sunshine, Version 2 (SMARTS2) [1416]. This model uses solar spectral irradiance data collected by satellites as a primary data source. This model was chosen due to its accuracy, reliability, and availability of open-source codes, allowing it to be adapted to local conditions. This model proved to be precise and accurate when compared with proprietary and sophisticated models [19], as well as other methods in the literature [20, 21]. The curves presented in Figure 3 were obtained using the SMARTS2 spectral radiation model [1416].

Typically, PV panels absorb in the wavelength () ranging from 300 nm to 1200 nm. Figure 4 shows the spectral response of a typical silicon solar cell [22]. The dots represent experimental data adapted from Ref. [22]. The solid line represents a polynomial adjustment as further explained in Section 3.2.

At wavelengths below 300 nm (ultraviolet), the glass absorbs most of the UV radiation and the electromagnetic response of the cells is very low. At intermediate wavelengths, the cell absorption increases with the wavelength and approaches the ideal performance. The response vanishes at long wavelengths (infrared) because silicon is a semiconductor with an indirect energy gap at room temperature [23], which is equivalent to a wavelength of 1107 nm. Therefore, all calculations to be carried out will use, when needed, the absorption range of silicon cells, i.e., from 300 nm to 1200 nm. Furthermore, the curve shown in Figure 4 should be used as a weighting factor of the spectral irradiance because the spectral response of a solar cell is not uniform for all wavelengths. Considering this, Equation (1) should be rewritten as follows: where is the spectral response of the photovoltaic cell, shown in Figure 4, i.e., acts as a weighting factor for the solar spectral irradiance . Because function is obtained experimentally, Equation (4) should be interpolated with the points presented in Figure 4 so that intermediate values are achieved in the integration. Thus, a high-order polynomial interpolation in multiple sectors should be used to ensure high accuracy with respect to the experimental points. Equation (4) can be related to the short-circuit current density ().

In order to determine the best solar panel tilt angle to obtain the maximum power throughout the year, one could simply discretize and scan the tilt variable , calculate the daily solar irradiation for each tilt angle for the entire year, and thus inspect the value of that provides the maximum annual solar irradiation. The main negative aspect of this approach is the unnecessary amount of calculations, overloading the computer resources. This approach can be avoided by adopting a more efficient procedure, which is determining the value of for which the derivative of the solar irradiation approaches zero. Therefore, one needs to calculate not only the solar irradiation given by Equation (3) but also its first derivative. In this study, a numerical fourth-order derivation method is used, i.e., the local truncation error is of the order of the step-size to the fourth power. This method uses Lagrange’s formula for five-point derivation [24].

Once the derivative is calculated, its root can be sought. The root of a function can be found by various numerical methods, for instance, the bisection method, secant method, Newton-Raphson method, among others [25]. There are other powerful methods, such as Ridders’ and Brent’s methods [25], but they incur significant computational resources. Other algorithms, such as bisection and secant, are computationally simple and inexpensive, but may not always work [25]. The chosen method for this work is a compromise between speediness and robustness, i.e., the method is simple and numerically effective, reaching a solution with very few iterations. The method used is called Regula-Falsi [25], also known as the False Position method. This method has been used successfully in atomic [2636] and solid-state physics [37, 38], in which the quantum mechanical equations call for very efficient numerical methods for minimizations. Thus, using this technique adds a novelty to the PV field.

2.2. Computational Details

This work required intensive use of a computer system for floating-point operations in double precision (REAL8) aimed at minimizing rounding and propagation errors. The developed codes and free SMARTS2 [1416] software that models the solar spectrum in the Earth’s atmosphere use the FORTRAN programming language, which is still widely adopted in the scientific community, although nowadays, there are other programming languages that offer more flexibility. Nonetheless, the SMARTS2 [1416] codes are only available in FORTRAN.

The SMARTS2 model requires input information about site location (geographical coordinates, altitude, azimuth, and tilt angle), date, and time. Furthermore, the model characterizes the local atmosphere (ozone concentration, aerosols, turbidity, and soil reflectance). Concerning the atmosphere, for Southern latitude calculations, we selected the SMARTS2 built-in “Tropical” standard atmosphere, which has typical gas concentrations without pollutants. In addition, the local environment was assumed as an urban area with a ground made of concrete. An assumption of a clear sky was made. Among a variety of output results, the data for spectral irradiances (direct, diffuse, and ground reflected) are returned to the main program to proceed with calculations.

3. Results

This section presents the results for solar irradiance and solar irradiation calculations from for the city of São Carlos (SP), Brazil, with geographic coordinates of latitude -22°0103 S, longitude -47°5327 W, and an average altitude of 854 m. The calculations were carried out for all days of the year and, for each day, the incident radiation throughout the day was considered, i.e., from dawn to sunset. These results were used to determine the best tilt angle of a PV solar panel to optimize the absorption of solar radiation. This section is divided into subsections to present the results more clearly.

3.1. Photoperiod

Daylight length is called photoperiod and comprises the elapsed time from sunrise to sunset. Therefore, the sunrise and sunset times were calculated and the photoperiod was obtained for each day. The calculations of these times were made using the SMARTS2 [1416] program by determining the time of day for which the sun elevation was zero, i.e., when the sun is on the horizon, more specifically, when half of the solar disk was below the horizon and the other half above. Figure 5 shows the photoperiod for each day of the year.

Figure 5 shows that daily sunlight lasts over 13 hours on summer days while it lasts slightly over 11 hours in the wintertime. The minimum period of sunlight extent occurs for day number 172, which corresponds to June 21. This day will be considered as half of the year and will be mentioned later as such.

3.2. Spectral Response

The first step was to fit the experimental points of the spectral response of a silicon solar cell, shown earlier in Figure 4. The fitting was done using polynomials up to the sixth order in four different spectral bands. These ranges are as follows (in nm): (1)(2)(3)(4)

These ranges were chosen to ensure minimum deviation in the fittings. For each of the above ranges, the fitting was of the sixth-order polynomial [25], each having, therefore, its own coefficients. Thus, the function obtained for the spectral response of the solar panel can be safely used instead of the experimental data, which are not available for all wavelength values. This polynomial [] was used to calculate the solar irradiance given by Equation (4).

3.3. Spectral Solar Irradiance

The spectral irradiance was calculated using the SMARTS2 [1416] codes for each value of and An example can be seen in Figure 6, which shows the global solar irradiance (top curve) at ground level, in São Carlos. This curve can be compared with Figure 3 and it should be noted that the model accurately reproduces the solar spectrum observed at ground level.

The bottom curve in Figure 6 shows the global solar spectral irradiance weighted by the spectral response of the solar panel, i.e., this curve is the product of the solar irradiance by the cell’s response for each wavelength. This weighted function should be integrated (compare Equation (4)) to give the total irradiance.

3.4. Total Solar Irradiance

The calculation of the spectral irradiance allows one to calculate its integral to obtain the total solar irradiance, in accordance with Equation (4). The global solar irradiance consists of direct, diffuse, and reflected irradiances. For instance, for the 172nd day of the year, which is June 21, the respective components of the global solar irradiances are shown in Figure 7. Note that, as expected, the main component is the direct irradiance on the solar panel. However, the diffuse and reflected components are not negligible. Collectively, they contribute about 20% of the overall irradiance.

The solar spectral irradiance weighted by the silicon cell’s response was used for the total solar irradiance calculation (see Figure 6 (bottom curve)).

3.5. Solar Irradiation

This section presents the essential results of this study. The calculation of the solar irradiation (Equation (3)) for different inclinations of the solar panel shows the superficial energy density absorbed by the plate with the specified geometry. Calculations were made considering three scenarios that maximize solar irradiation (sometimes referred to as radiant exposure) into the panel depending on its tilt angle. The first scenario is the most general case in which the global solar irradiation throughout the year is considered, henceforth referred to as Period I. The second case (Period II) considers a six-month period centered in the day with less solar irradiation, i.e., June 21, which is day 172 of the year and its midpoint, as mentioned in section III A. Thus, the period considered for the calculation of this case spans three months before day 172 and three months after day 172, ranging from day 80 to day 264, summing up 184 days between March 21 (day 80) to September 21 (day 264). This six-month period refers to the autumn and winter seasons in the Southern hemisphere. Finally, the third case (Period III) is concentrated in three months centered on day 172, i.e., from day 127 to day 217, corresponding to May 7 to August 5, respectively. These three scenarios are good choices to understand the dynamics of PV tilt angles throughout the year. Period I is the overall condition while Periods II and III are good in cases that are not connected with the grid (subdimensioned and/or with battery banks).

In order to generalize our approach, we provide results for Southern latitudes ranging from 0° to −55° in steps of −5° considering four different periods. These are the standard meteorological seasons [39], which are more appropriate than the astronomical seasons (for this kind of calculation) because meteorological seasons are based on the annual temperature cycle.

As a first example, given a 30° tilt for the solar panel, the solar irradiation for each day of the year can be seen in Figure 8. This figure shows the direct, diffuse, and reflected components, which are added together to obtain the global solar irradiation.

As expected, the solar irradiation is smaller in the central part of the graph, which corresponds to half of the year, namely, to the winter period (Southern hemisphere). The two maxima correspond to late summer and early spring. More precisely, they refer to the southern equinoxes. Figure 8 is a solar irradiation plot for each day of the year. Adding the values together for each day of the year, the annual superficial energy density is obtained. This amount is the total that the solar panel will be submitted during the year for the considered geometric condition. Thus, for a 30° tilt of the PV panel, the total annual solar irradiation is 4758.1 MJ/m2. Recall that this value takes into account the spectral response of the PV panel, shown in Figure 4.

3.5.1. Maximization by Calculating the First Derivative

For all considered scenarios (Periods I to III for São Carlos and meteorological seasons for Southern latitudes ranging from 0° to −55°), the tilt angle of the solar panel is determined by the maximum solar irradiation within the period of each scenario.

The Regula-Falsi method for obtaining the root of the derivative of the solar irradiation has proven to be very numerically efficient because the root was obtained with an accuracy of eight decimal digits after only six iterations of the algorithm. Evidently, there is no need for such precision in the tilt angle. One decimal digit is adequate for the physical installation of the panel. Convergence to one decimal digit would require only three iterations of the algorithm. Considering that the first two iterations are initial attempts required by the method, just one additional iteration would be sufficient. Therefore, the Regula-Falsi method proved to be extremely effective in determining the roots. The first two iterations are arbitrary depending on the desired range. In this paper, we have chosen and to start the iterations.

Table 1 shows the optimized tilt angle of a PV panel to maximize the absorption of solar radiation in the calculated periods.

Table 1 indicates the correct positioning for PV panels, for the town of São Carlos (SP) in Brazil, in order to obtain the maximum radiation absorption in the periods of the three considered scenarios. Note that the panel should always be directed to the geographic North and should not be determined by a compass. Table 1 also shows the total solar irradiation in each considered period and the gain in relation to the positioning of the panel at 22°, which is the default angle for São Carlos. It should be mentioned that this default angle is equal to the latitude installation site. Although the percentage gain for Period I may look small, we emphasize that usually PV panels are rarely tilted at 22° (São Carlos) and not always directed to the North as they are placed on top of roofs that are not designed for this purpose. Indeed, roof inclinations in São Carlos are commonly designed in the 2°–20° range, depending on the roof covering.

3.5.2. Maximization by Direct Calculation

The calculation of maximization by direct calculation is redundant as the best tilt angle was determined in the previous section. However, the direct calculation is aimed at demonstrating the consistency of the results and ensuring that the searching method for the root of the derivative actually determines the point at which the total solar irradiation is maximum, i.e., its derivative with respect to the tilt angle of the solar panel is zero. The direct calculation is nothing more than calculating the annual solar irradiation to fixed intervals of the panel tilt angle. This allows the plot of the solar irradiation incident on the panel as a function of the tilt angle for the considered scenarios. Figure 9 shows the calculation for Period I as an example.

Note that the direct component reaches its maximum just before the global maximum, showing that the diffuse and reflected components have no negligible contributions in the overall calculation. To have a clearer understanding, Figure 10 shows the global solar irradiation incident in a solar panel as a function of inclination for the three considered periods.

The shift of the peak towards larger tilt angles is due to the maximum solar elevation of the days in the year that becomes smaller as the winter period approaches, causing the solar panel to be increasingly tilted in order to increase the radiation absorption.

3.5.3. Comparison and Generalization

It can be observed in Figure 10 that, indeed, the Regula-Falsi method was effective in determining the tilt angle to maximize solar irradiation. Visual inspection also shows the accuracy of the method. Comparing the data in Table 1 with Figure 10, the consistency of the results can be observed, i.e., the values shown in Table 1 match the global maximum points of the respective curves.

Figure 11 shows the daily solar irradiation incident on the solar panel for each of the periods. For each period, the calculation was performed using the best solar irradiation panel tilt angle listed in Table 1, calculated using the proposed method. For shorter periods, the curves move upward, indicating that the solar irradiation is maximized as the tilt angle increases, and thus increases the absorption of solar radiation.

The curves in Figure 11 are the global solar irradiation, that is, the sum of the direct, diffuse, and reflected components. Figure 11 also shows the global irradiation for standard inclination, i.e., the latitude tilt angle, for the sake of comparison.

In order to demonstrate the capabilities and versatility of this approach, Figure 12 shows results for Southern latitudes ranging from 0° to −55° in steps of −5° for the meteorological seasons.

4. Conclusions

The irradiance and solar irradiation incident on PV panels were calculated, and the best panel tilt angle was determined from the maximum absorption of solar radiation in the chosen periods.

For irradiance and solar irradiation calculations, we proposed to use a diffusion and atmospheric radiative transfer model, which is well known and validated in the literature. Moreover, the solar irradiation is weighted by the response curve of the silicon PV panels. Subsequently, to determine the tilt angle of PV panels that retrieves the highest possible incident solar irradiation, we proposed the Regula-Falsi numerical method in order to find the point at which the derivative of the solar irradiation with respect to the tilt angle approaches zero. These three aspects are novelties to the field. This combination of methods is unique in the field. Furthermore, this is a strictly theoretical approach that does not rely on collected data, turning our approach universal.

As an application, the calculations were carried out for the city of São Carlos. The optimal tilt angle of a PV panel was obtained by the proposed method for three selected periods. Additionally, optimal tilt angles were provided for Southern latitudes and periods corresponding to the meteorological seasons.

For each considered period, this work showed an increase in solar radiation absorption compared to the general installation angle, which at best is equal to the latitude of the location. The highest gain was for Period III (winter season), which amounts to over 14% compared to the standard latitude tilt angle. At the level of GJ/m2, this improvement is significant.

Finally, the computer codes developed are able to run autonomously, with minimum user intervention. The data set for each location of interest is provided to the main program as input files that can be edited by the user. Thus, these calculations can be extended to anywhere in the Southern or Northern hemispheres. Moreover, different desired periods of investigation are achievable. For instance, daily, weekly, or monthly values of optimal tilt angles can be provided for photovoltaic arrays.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Aside from the SMARTS2 code, MM developed the other codes. In the remaining of this work, the authors have worked and contributed equally.

Acknowledgments

The authors are grateful to FAPESP (grant numbers 2013/08038-7 and 2014/16938-0), which financially supported our research. This study was financed in part by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES) (Finance Code 001 (MM)).