Analysis of the Inﬂuence of Terrain Orientation on the Design of PV Facilities with Single-Axis Trackers

: This paper investigates how to optimally orient the photovoltaic solar trackers of an axis parallel to the terrain, applying the sky model of Hay–Davies. This problem has been widely studied. However, the number of studies that consider the orientation (inclination and azimuth of the terrain) is very limited. This paper provides an examination of incident solar irradiance that can be extended to terrain with variable orientation and in consideration of di ﬀ erent azimuths of the axis of rotation. Furthermore, a case study of the south of Spain is provided, considering di ﬀ erent inclination and orientation terrain values. The results obtained in this study indicate, as a novelty, that for lands that are not south facing, the rotation axis azimuth of solar trackers should be di ﬀ erent from zero and adjusted to the same direction as the land azimuth in order to maximize energy production. Annual energy production is sensitive to changes in the rotation axis azimuths of solar trackers (an inﬂuence of around 3% of annual energy production).


Introduction
At present, population growth, deficiencies in natural resources and global warming are producing challenges around the world. In recent years, various initiatives have been working on the development or replacement of the current power grid to obtain greater efficiency in the electrical system and a reduction in energy waste by reducing losses during distribution [1]. Such a system would contribute to the promotion of renewable resources and the minimisation of the environmental impact in the future. Therefore, there is a global requirement to preserve the environment and to improve the penetration of alternative energy resources [2]. Within the framework of renewable energies, solar energy is the form of energy that has shown the most remarkable growth in recent years thanks to the reduction of costs and to legal requirements in many countries [3,4]. Similarly, photovoltaic (PV) technology is one of the current technologies with a better future projection due to its simplicity and scalability and the continuous manufacture, operation and maintenance reduction costs it enables for solar panels [5].
However, the lack of the linearity in the solar energy received by solar panels, mainly caused by Earth-Sun relative displacement, is a disadvantage to consider. It is necessary to redefine solar Backtracking has become a technique used to increase energy production in photovoltaic plants in cloudy periods. The diffuse component of solar radiation is much more relevant in these cases and "approaches" that ignore this component (such as astronomical ones) are thus no longer optimal. This is not the only advantage of employing backtracking. There is also the resulting reduction in hot spots on the trackers, which increases their life span and reduces breakdowns. Backtracking also improves the ground cover ratio (GCR) of the installations, which has the effect of reducing the economic impact with regard to the soil required when building plants [27].
For these reasons, there is a niche in scientific work focused on the development of options utilising backtracking and on the optimisation of solutions responding to the challenges posed by solar panels which are not exclusively focused on classical tracking [9,26,28,29].
From this perspective, work based on tracking has been undertaken with the aim of optimising production throughout the movement of a panel path, including the backtracking periods [27]. To achieve this, the search is underway for a panel path that maintains the optimum production under the premise that the panel is isolated and therefore does not receive shade from any adjacent panel. Once this objective has been considered, a condition can be established to verify the possible shading between photovoltaic panels. After these two steps, the shaded orientations can be indicated against those that are not shaded. Thus, a panel path that optimises production, including the non-shading of modules, can be established.
Since it is difficult to find a completely horizontal location for photovoltaic plants, further study and deeper knowledge of photovoltaic trackers are required to offer optimisation guidelines for solar capture. The present study was undertaken to include: (a) mathematical modelling of solar capture on trackers; (b) programming that allows automation of the calculations; (c) simulations of case studies; and (d) analysis of optimal conclusions. In light of this, conclusions were drawn on orientations, land inclinations and distances between modules, etc., in terms of produced energy. Figure 2a,b represent angles χ (azimuth) and β (inclination) for a generic terrain. The angle γ (azimuth of the collector rotation axis) is also represented. This study, based on the Hay-Davies model, characterises solar radiation on trackers using angles χ, β and γ. The remainder of this paper is organised as follows. In Section 2 the model used, the mathematical optimisation bases, the backtracking method and the relevant software development Backtracking has become a technique used to increase energy production in photovoltaic plants in cloudy periods. The diffuse component of solar radiation is much more relevant in these cases and "approaches" that ignore this component (such as astronomical ones) are thus no longer optimal. This is not the only advantage of employing backtracking. There is also the resulting reduction in hot spots on the trackers, which increases their life span and reduces breakdowns. Backtracking also improves the ground cover ratio (GCR) of the installations, which has the effect of reducing the economic impact with regard to the soil required when building plants [27].
For these reasons, there is a niche in scientific work focused on the development of options utilising backtracking and on the optimisation of solutions responding to the challenges posed by solar panels which are not exclusively focused on classical tracking [9,26,28,29].
From this perspective, work based on tracking has been undertaken with the aim of optimising production throughout the movement of a panel path, including the backtracking periods [27]. To achieve this, the search is underway for a panel path that maintains the optimum production under the premise that the panel is isolated and therefore does not receive shade from any adjacent panel. Once this objective has been considered, a condition can be established to verify the possible shading between photovoltaic panels. After these two steps, the shaded orientations can be indicated against those that are not shaded. Thus, a panel path that optimises production, including the non-shading of modules, can be established.
Since it is difficult to find a completely horizontal location for photovoltaic plants, further study and deeper knowledge of photovoltaic trackers are required to offer optimisation guidelines for solar capture. The present study was undertaken to include: (a) mathematical modelling of solar capture on trackers; (b) programming that allows automation of the calculations; (c) simulations of case studies; and (d) analysis of optimal conclusions. In light of this, conclusions were drawn on orientations, land inclinations and distances between modules, etc., in terms of produced energy. Figure 2a,b represent angles χ (azimuth) and β (inclination) for a generic terrain. The angle γ (azimuth of the collector rotation axis) is also represented. This study, based on the Hay-Davies model, characterises solar radiation on trackers using angles χ, β and γ. Backtracking has become a technique used to increase energy production in photovoltaic plants in cloudy periods. The diffuse component of solar radiation is much more relevant in these cases and "approaches" that ignore this component (such as astronomical ones) are thus no longer optimal. This is not the only advantage of employing backtracking. There is also the resulting reduction in hot spots on the trackers, which increases their life span and reduces breakdowns. Backtracking also improves the ground cover ratio (GCR) of the installations, which has the effect of reducing the economic impact with regard to the soil required when building plants [27].
For these reasons, there is a niche in scientific work focused on the development of options utilising backtracking and on the optimisation of solutions responding to the challenges posed by solar panels which are not exclusively focused on classical tracking [9,26,28,29].
From this perspective, work based on tracking has been undertaken with the aim of optimising production throughout the movement of a panel path, including the backtracking periods [27]. To achieve this, the search is underway for a panel path that maintains the optimum production under the premise that the panel is isolated and therefore does not receive shade from any adjacent panel. Once this objective has been considered, a condition can be established to verify the possible shading between photovoltaic panels. After these two steps, the shaded orientations can be indicated against those that are not shaded. Thus, a panel path that optimises production, including the non-shading of modules, can be established.
Since it is difficult to find a completely horizontal location for photovoltaic plants, further study and deeper knowledge of photovoltaic trackers are required to offer optimisation guidelines for solar capture. The present study was undertaken to include: (a) mathematical modelling of solar capture on trackers; (b) programming that allows automation of the calculations; (c) simulations of case studies; and (d) analysis of optimal conclusions. In light of this, conclusions were drawn on orientations, land inclinations and distances between modules, etc., in terms of produced energy. Figure 2a,b represent angles χ (azimuth) and β (inclination) for a generic terrain. The angle γ (azimuth of the collector rotation axis) is also represented. This study, based on the Hay-Davies model, characterises solar radiation on trackers using angles χ, β and γ. The remainder of this paper is organised as follows. In Section 2 the model used, the mathematical optimisation bases, the backtracking method and the relevant software development The remainder of this paper is organised as follows. In Section 2 the model used, the mathematical optimisation bases, the backtracking method and the relevant software development are presented. Section 3 outlines how the methodology was tested in Córdoba, Spain, and the results are discussed. In Section 4, conclusions are drawn based on the restrictions shown and potential further research is indicated.

Astronomical Bases and Irradiance Model
Knowing the concrete position of the sun at each moment of the day and on each day of the year is fundamental for determining the direction of solar beams. It is possible to determine the position of the sun accurately by using astronomical geometry. For this purpose, a solar vector was defined, expressed with respect to a terrestrial reference system ( Figure 2) and formed by the axes Ox, towards the west; Oy, towards the south; and Oz in the zenithal direction. The position was determined by Equations (1) and (2): where δ represents the solar declination, ϕ the latitude and Ωt the hour angle. This angle was calculated as the Earth's rotation speed, 2π/24 rad/h, for the hours since solar noon [30].
In PV technology, all irradiance components (direct, diffuse and reflected irradiance) are usable. Traditionally, the astronomical model has been used to determine the position of the sun. Applying the astronomical model to solar tracking means that the angle formed between the direct solar rays and the normal angle to the collector's surface θ must be as low as possible. With astronomical tracking, the value of the direct irradiance component is maximised, which is appropriate for applications focused on this component (such as concentration technologies). Therefore, this type of tracking is not the most suitable for the purposes of this study. For instance, on cloudy days, when the solar disk is not visible and direct radiation does not reach the collectors, those located on a fixed horizontal position collect more energy than those using astronomical tracking. Therefore, to study the influence of diffuse and reflected components on solar tracking, several sky models have been proposed.
The Hay-Davies method [31] was used for the calculations in this work because it adequately describes the anisotropy of radiation [32] in addition to being notable for its simplicity, as opposed to other more complex models such as those of Muneer [33] and Perez [34], and for obtaining high quality results [35][36][37]. This model establishes that a determined fraction of diffuse irradiance is directed according to the direction of the solar disk. Models such as that of Hay-Davies [31] describe the irradiance affecting a solar panel by considering its three components, direct, diffuse and reflected irradiance, but take into account the fact that the diffuse component has a preferential direction in which the diffuse radiation is greater. According to Hay [31], the mathematical expression of the model is given by Equation (3), where the first term indicates the direct irradiance I B , the second term refers to the diffuse irradiance I D and the third term corresponds to the reflected component, taking into account the visible soil fraction and the albedo, ρ. I B and I D are determined by the Collares-Pereira model [38].
where ξ is the inclination angle of the collector and θ Z is the solar zenith angle, that is, the angle between the solar vector and the zenith axis.

Optimisation of Collector Orientation
The aim of this study was to establish guidelines for designing facilities with solar trackers on a rotation axis → e . This axis was considered as parallel to the terrain; no restrictions were initially imposed on it. The general results are valid for any disposition.
In this study there were two degrees of freedom for the orientation of a surface, the azimuth χ and inclination β. In the case of single-axis trackers, there would only be the elevation. The orientation of a panel was characterised by a single vector → n , normal to the surface of the panel. The objective was therefore the calculation of → n such that that the solar capture was the maximum. Expanding on the treatment of the vectors, the cosines of the angles θ, θ Z and ξ that appear in the expressions of the different models of the irradiance are expressed as the scalar products of Substituting the previous expressions into the Hay-Davies model [31] gives the following: To optimise the irradiance n , the Lagrange multiplier method was used with the following restrictions defined: The differential is described by: by matching the brackets to zero.
By naming three vectors → u , → n and → e appear. These vectors fulfil the relation: So, they are linearly dependent and therefore coplanar.
As → n is perpendicular to → e , it can be expressed more easily. Replacing the vectors ⃗ and ⃗ with their expressions results in the equation that defines the vector ⃗: This result makes it possible to calculate the optimal direction of the normal vector for any irradiance model in any direction of the rotation axis and for different positions of the sun. By replacing the values of the partial derivatives in Equation (14) with those obtained from Equation (7), we get: Utilising Equation (16), the values of ⃗ can be obtained at any given time for the optimum irradiance incidence.
For example, for the simplest case of a direct irradiance model which only takes into account the first term of Equation (3), and with the rotation axis of the solar panels assumed to be horizontal to the ground, a normal vector like the following results, where the vector ⃗ is: being the inclination of the terrain, the azimuth of the terrain and the azimuth of the rotation axis collector.

Backtracking
This result makes it possible to calculate the optimal direction of the normal vector for any irradiance model in any direction of the rotation axis and for different positions of the sun. By replacing the values of the partial derivatives in Equation (14) with those obtained from Equation (7), we get: Utilising Equation (16), the values of → n can be obtained at any given time for the optimum irradiance incidence.
For example, for the simplest case of a direct irradiance model which only takes into account the first term of Equation (3), and with the rotation axis of the solar panels assumed to be horizontal to the ground, a normal vector like the following results, where the vector β being the inclination of the terrain, χ the azimuth of the terrain and γ the azimuth of the rotation axis collector.

Backtracking
In the previous section, the methodology used to determine the orientation of the solar tracker in order to obtain the maximum solar capture was detailed. This section describes the method used to prevent a solar tracker from shading an adjacent tracker. Only two factors needed to be considered for this method: the position of the sun and the orientation of the panels. To study intershading it is convenient to use an auxiliary reference system where the representation plane is perpendicular to the tracker axis. For this purpose, the vector → q was defined as: where → q is a vector perpendicular to the ground normal vector convenient to use an auxiliary reference system where the representation plane is perpendicular to the tracker axis. For this purpose, the vector ⃗was defined as: where ⃗ is a vector perpendicular to the ground normal vector ⃗, as well as to the rotation axis vector ⃗. The solar vector ⃗and the plane normal vector in the optimal position ⃗ were named ′ ⃗ and ′ ⃗ in the reference system composed by ⃗, ⃗, ⃗ (Figure 4).  convenient to use an auxiliary reference system where the representation plane is perpendicular to the tracker axis. For this purpose, the vector ⃗was defined as: where ⃗ is a vector perpendicular to the ground normal vector ⃗, as well as to the rotation axis vector ⃗. The solar vector ⃗and the plane normal vector in the optimal position ⃗ were named ′ ⃗ and ′ ⃗ in the reference system composed by ⃗, ⃗, ⃗ (Figure 4).  In order to obtain the vector ⃗, both terms were multiplied by the vector ′ ⃗ : In order to obtain the vector → a , both terms were multiplied by the vector → n : Substituting the value of τ from Equation (19) into Equation (18) results in: Vector → a can be obtained according to the following expression: As mentioned above, the examination of vector → a enables the determination of whether shading will occur. If the module of the vector is less than the width of the panel, shading will occur.
To avoid this situation, backtracking was used to look for a new orientation of the panels that prevented shading. The new orientation was intended to involve a minimum deviation from the optimal capture orientation. The new inclination angle α was that which caused the shading of a panel to fall on the edge of the adjacent panel ( Figure 6).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 16 Vector ⃗can be obtained according to the following expression: As mentioned above, the examination of vector ⃗ enables the determination of whether shading will occur. If the module of the vector is less than the width of the panel, shading will occur.
To avoid this situation, backtracking was used to look for a new orientation of the panels that prevented shading. The new orientation was intended to involve a minimum deviation from the optimal capture orientation. The new inclination angle was that which caused the shading of a panel to fall on the edge of the adjacent panel ( Figure 6).

This inclination angle
was calculated based on a triangle formed by the line that joins the center of the two panels , the projection of the centre of the panel on the line of the adjacent ′ ⃗ and the distance that corresponds to the module of vector ⃗ in Figure 5, which in Figure 6 is equivalent to the width of the panels, ℎ.
The angle is the one between the solar vector ′ ⃗ and the axis: The value of the angle can already be calculated as: Therefore, the inclination angle of the solar panel can be calculated as follows: The normal vector that leads to the determination of the orientation of the panel in order to avoid shading is:  Figure 5, which in Figure 6 is equivalent to the width of the panels, h.
The angle ζ is the one between the solar vector → s and the X axis: The value of the angle δ can already be calculated as: Therefore, the inclination angle of the solar panel can be calculated as follows: The normal vector that leads to the determination of the orientation of the panel in order to avoid shading is:

Software Applications for Analysis
To calculate the irradiance received by each solar panel, a function was implemented in Visual Basic for Applications (VBA) for Excel in which the calculations described above were developed. As seen in Figure 7, the function was based on several input parameters (latitude, Julian day, geometrical values, irradiance model considered) and provided an irradiance value every three minutes. Some restrictions, such as whether it as daytime, the maximum radiation related to the used model and eventual shading with backtracking response were assumed.

Software Applications for Analysis
To calculate the irradiance received by each solar panel, a function was implemented in Visual Basic for Applications (VBA) for Excel in which the calculations described above were developed. As seen in Figure 7, the function was based on several input parameters (latitude, Julian day, geometrical values, irradiance model considered) and provided an irradiance value every three minutes. Some restrictions, such as whether it as daytime, the maximum radiation related to the used model and eventual shading with backtracking response were assumed.

Results and Discussion
The practical results provided the value of the estimated annual radiation using the Hay-Davies model [31] for various land configurations. All the calculations were made for the province of Córdoba, Spain, therefore assuming the same geographical values (latitude = 37.75492 • N; longitude = 5.04548 • W) and the same climatic values throughout the year, varying only the values of inclination β and azimuth χ for a fixed panel width h = 3 m and separation between panel lines d = 6 m.
Terrains with negative inclinations were not considered since PV plants are not commonly installed in shaded areas.
For facilities located on terrain with a null azimuth (oriented south), independently of its inclination even if it is null, the optimal orientation of the axis of rotation is southward (γ = 0 • ).
For each value of the inclination of the terrain β, it was also observed that the maximum value of irradiance was obtained when the azimuth of the terrain χ was 0 • . That is, the optimal situation is when the terrain is oriented south and the axis of rotation is directed in the same direction. If the terrain has another orientation, the value of the optimal angle towards which the axis of rotation should be oriented would move in that same direction. This fact is more evident at higher values of the terrain inclination β, although in no case will the azimuth axis γ be as high as that of the terrain χ.
The results for a case study with a β = 15 • inclination are presented in Figure 8. Figure 8a shows the variation of radiation with regard to γ when the azimuth of the terrain χ = 0 • . However, when the terrain was not oriented to the south, as shown in Figure 8b in which the terrain was turned χ = 30 • , the maximum radiation was obtained by orienting the axis with an azimuth γ of 6 • towards the south. The losses with regard to the maximum radiation obtained in a terrain are shown in Figure 8c, where it is again verified that, by orienting the axis at 6 • , these losses can be minimised. Finally, Figure 8d shows the optimal orientation to direct the axis for the different values of the terrain azimuth χ. In this way it is possible to estimate the orientation that achieves optimal radiation considering the azimuth of the terrain χ.
As already mentioned, for each inclination value the maximum irradiance value was obtained when the terrain had an angle χ = 0 • . In addition, as the inclination of the terrain increased, the radiation affecting the panels was greater, as seen in Figure 9. The explanation can be found in the latitude for which the calculations were made and in the fact that the direction of the sun is not perpendicular to the horizontal plane. The inclination of the terrain influenced the radiation, increasing until a value of β of around 21 • , where the maximum was reached.

Results and Discussion
The practical results provided the value of the estimated annual radiation using the Hay-Davies model [31] for various land configurations. All the calculations were made for the province of Córdoba, Spain, therefore assuming the same geographical values (latitude = 37.75492° N; longitude = 5.04548° W) and the same climatic values throughout the year, varying only the values of inclination and azimuth for a fixed panel width ℎ = 3 and separation between panel lines = 6 . For each inclination (from 0° to 45°, 1° increment) and terrain azimuth (from −60° to 60°, 5° increment), the rotation axis was oriented with a wide range of azimuth γ (from −20° to 20°, 2° increment).
Terrains with negative inclinations were not considered since PV plants are not commonly installed in shaded areas.
For facilities located on terrain with a null azimuth (oriented south), independently of its inclination even if it is null, the optimal orientation of the axis of rotation is southward (γ = 0°).
For each value of the inclination of the terrain , it was also observed that the maximum value of irradiance was obtained when the azimuth of the terrain was 0°. That is, the optimal situation is when the terrain is oriented south and the axis of rotation is directed in the same direction. If the terrain has another orientation, the value of the optimal angle towards which the axis of rotation should be oriented would move in that same direction. This fact is more evident at higher values of the terrain inclination , although in no case will the azimuth axis γ be as high as that of the terrain .
The results for a case study with a = 15° inclination are presented in Figure 8. Figure 8a shows the variation of radiation with regard to γ when the azimuth of the terrain = 0°. However, when the terrain was not oriented to the south, as shown in Figure 8b in which the terrain was turned = 30°, the maximum radiation was obtained by orienting the axis with an azimuth γ of 6° towards the south. The losses with regard to the maximum radiation obtained in a terrain are shown in Figure 8c, where it is again verified that, by orienting the axis at 6°, these losses can be minimised. Finally, Figure 8d shows the optimal orientation to direct the axis for the different values of the terrain azimuth . In this way it is possible to estimate the orientation that achieves optimal radiation considering the azimuth of the terrain . (c) (d) Figure 8. Results for a case study for inclination 15°: (a) maximum radiation; (b) radiation variation for 30° azimuth; (c) maximum radiation loss for 30° azimuth; (d) axis azimuth for optimal terrain radiation.
As already mentioned, for each inclination value the maximum irradiance value was obtained when the terrain had an angle = 0°. In addition, as the inclination of the terrain increased, the radiation affecting the panels was greater, as seen in Figure 9. The explanation can be found in the latitude for which the calculations were made and in the fact that the direction of the sun is not perpendicular to the horizontal plane. The inclination of the terrain influenced the radiation, increasing until a value of of around 21°, where the maximum was reached. The difference observed between the radiation value for a particular axis azimuth and the maximum capture obtained in the optimum direction was greater as the inclination of the terrain increased. Thus, in low-inclination terrains the difference was low and while for steeper inclinations the losses were considerable. Figure 10 shows that, in the specific case of = 20°, the losses caused by a bad choice of rotation axis direction can reach values above 7%. (c) (d) Figure 8. Results for a case study for inclination 15°: (a) maximum radiation; (b) radiation variation for 30° azimuth; (c) maximum radiation loss for 30° azimuth; (d) axis azimuth for optimal terrain radiation.
As already mentioned, for each inclination value the maximum irradiance value was obtained when the terrain had an angle = 0°. In addition, as the inclination of the terrain increased, the radiation affecting the panels was greater, as seen in Figure 9. The explanation can be found in the latitude for which the calculations were made and in the fact that the direction of the sun is not perpendicular to the horizontal plane. The inclination of the terrain influenced the radiation, increasing until a value of of around 21°, where the maximum was reached. The difference observed between the radiation value for a particular axis azimuth and the maximum capture obtained in the optimum direction was greater as the inclination of the terrain increased. Thus, in low-inclination terrains the difference was low and while for steeper inclinations the losses were considerable. Figure 10 shows that, in the specific case of = 20°, the losses caused by a bad choice of rotation axis direction can reach values above 7%. The difference observed between the radiation value for a particular axis azimuth and the maximum capture obtained in the optimum direction was greater as the inclination of the terrain increased. Thus, in low-inclination terrains the difference was low and while for steeper inclinations the losses were considerable. Figure 10 shows that, in the specific case of β = 20 • , the losses caused by a bad choice of rotation axis direction can reach values above 7%.  Figure 11 shows the losses for the conditions γ = 0° and = 20°, considering an azimuth range from −60° to 60°. For terrains that are also oriented towards the south, or close to it, the losses will be zero or very low, but as the azimuth of the terrain increases, the losses become considerable.  Figure 11 shows the losses for the conditions γ = 0 • and β = 20 • , considering an azimuth range χ from −60 • to 60 • . For terrains that are also oriented towards the south, or close to it, the losses will be zero or very low, but as the azimuth of the terrain increases, the losses become considerable.  Figure 11 shows the losses for the conditions γ = 0° and = 20°, considering an azimuth range from −60° to 60°. For terrains that are also oriented towards the south, or close to it, the losses will be zero or very low, but as the azimuth of the terrain increases, the losses become considerable. Figure 11. Losses with 0° axis azimuth and 20° inclination.
Thus, it is possible to establish a relationship between the sensitivity of the collector axis azimuth γ and the inclination and azimuth of the terrain . It is worth highlighting the novelty of this result since no references have been found in the literature suggesting this relationship between variables.
Another key factor when configuring the layout of the collectors was the separation between them. This fact had an influence since the greater the panel separation is, the fewer shading effects are produced and the fewer losses due to a decrease in production. On the other hand, it was not possible to separate the rows of collectors as much as desired, since the costs associated with a higher terrain occupation increase rapidly [21]. Figure 12 shows that, depending on this parameter, the irradiance demonstrated asymptotic behaviour. This was because the shadow effects of one collector on another, when separated by a sufficient distance, were no longer perceived, thus resulting in an irradiance value like that of an isolated collector without adjacent collectors. Thus, it is possible to establish a relationship between the sensitivity of the collector axis azimuth γ and the inclination β and azimuth of the terrain χ. It is worth highlighting the novelty of this result since no references have been found in the literature suggesting this relationship between variables.
Another key factor when configuring the layout of the collectors was the separation between them. This fact had an influence since the greater the panel separation is, the fewer shading effects are produced and the fewer losses due to a decrease in production. On the other hand, it was not possible to separate the rows of collectors as much as desired, since the costs associated with a higher terrain occupation increase rapidly [21]. Figure 12 shows that, depending on this parameter, the irradiance demonstrated asymptotic behaviour. This was because the shadow effects of one collector on another, when separated by a sufficient distance, were no longer perceived, thus resulting in an irradiance value like that of an isolated collector without adjacent collectors.

Conclusions
In this paper, the solar capture on horizontal one-axis trackers was mathematically modelled, thus establishing the rotation axis orientation for any irradiance model. Backtracking was used to avoid the shading of some panels on adjacent ones, modifying their orientation. The results obtained, combined with those in the literature, amount to an advance in the knowledge on photovoltaic trackers with axes parallel to the terrain, allowing an optimal design for terrains without the horizontality described in previous works.
A program and a subroutine were implemented in VBA to automatically calculate the radiation affecting a solar tracker for different terrain configurations and for which it was possible to vary the characteristics of the installation, namely the rotation axis orientation and the ratio dimension of the

Conclusions
In this paper, the solar capture on horizontal one-axis trackers was mathematically modelled, thus establishing the rotation axis orientation for any irradiance model. Backtracking was used to avoid the shading of some panels on adjacent ones, modifying their orientation. The results obtained, combined with those in the literature, amount to an advance in the knowledge on photovoltaic trackers with axes parallel to the terrain, allowing an optimal design for terrains without the horizontality described in previous works.
A program and a subroutine were implemented in VBA to automatically calculate the radiation affecting a solar tracker for different terrain configurations and for which it was possible to vary the characteristics of the installation, namely the rotation axis orientation and the ratio dimension of the panels h and the separation d between adjacent rows. Simultaneously, the program and subroutine provided the optimal azimuth γ value for the rotation axis to be suitably directed and the loss values for each possible maximum radiation value for each inclination β of the terrain.
The calculations of the amounts of radiation affecting solar panels, assuming fixed width and distance between lines of panels and varying terrain inclinations β and azimuth χ, were made in Córdoba, Spain.
For any given value regarding the terrain inclination β, it was verified that a higher value of radiation was obtained for terrain that was oriented south (χ = 0 • ). The optimum direction for the rotation axis collector γ should also be southwards.
An important contribution of this research is the relationship established between the azimuth γ of the collector axis and the inclination β and azimuth of the terrain. For terrain that is not south facing, the rotation axis azimuth γ of the solar trackers should be different from zero and varies in the same direction as the land azimuth χ. If the axis of rotation is positioned in a different direction to the optimum, losses in potential energy production arise. These losses vary according to the terrain inclination β. The greater the difference between the optimum azimuth γ and the terrain inclination β, the greater the production losses will be.
Another important conclusion drawn from this research concerns the collector layout on the power plant surface. As the distance d between rows of collectors increases, a greater amount of radiation is generated. Such growth is asymptotic.
Overall, some progress has been achieved towards the objective of deepening knowledge concerning photovoltaic trackers. Consequently, guidelines were proposed for solar capture optimisation in photovoltaic plants. This work opens up further lines of research into the geometric layout optimisation of solar trackers involving more complex models, such as those of Perez [34] and Muneer [33], and real databases.
As a result of the aforementioned conclusions, it is considered that this study can be the basis for further work, such as the study of one-axis solar trackers located on terrains with irregular topographies or of the collector distribution system, considering potential uses of the land.

Conflicts of Interest:
The authors declare no conflict of interest.