Study of the Dependence of Solar Radiation Regarding Design Variables in Photovoltaic Solar Installations with Optimal Dual-Axis Tracking

: Solar tracking is an efﬁcient strategy to increase the radiative capture of photovoltaic collectors. Within the multiple efforts made in recent decades to improve the production of these facilities, various works have studied solutions to optimize the number of rotation axes (single or dual rotation axes), the degree of collector coverage, the distances between trackers, the geometric arrangement of trackers or the minimization of shading between collectors. However, although in this type of installation it is common to ﬁnd collectors with geometric shapes other than rectangles, no studies on the inﬂuence of the shape of the collectors on the radiative incidence are found in the literature. In this connection, the present work systematically addresses the study of incident solar radiation in photovoltaic installations with dual-axis trackers with collectors of different geometric shapes. By means of the exhaustive study, the conclusion is drawn that, for dual-axis photovoltaic installations with an optimal tracking strategy, the main variables that inﬂuence the annual radiative incidence are the spacing between collectors, the coverage ratio (GCR), and the collector surface, while the type of arrangement of collectors and the shape of these do not show predictive values. the main conclusions of the present study are presented.


Introduction
There is no doubt about the important role that energy plays in our societies. Its implications go beyond mere technical aspects. Social structuring, economic development, and the environment, among others, are aspects that make up the complex implication of energy in the global agenda [1]. In this context, renewable energies have been experiencing sustained growth in recent years. Thus, the International Energy Agency (IEA) foresees a record increase of 218 GW in the year 2021 in the net capacity of renewable electricity installed in the world, in an average scenario, that could reach up to 266 GW, in an accelerated scenario [2]. In fact, during the first quarter of 2020, renewable energies were the only source of electricity whose demand increased despite the 2.5% decrease in global electricity demand caused by the blockades implemented by different governments to curb the spread of COVID-19 [3].
Within the field of renewables, solar energy, in general, and photovoltaic (PV), in particular, are candidates to satisfy a large part of the global energy demand in the coming years due to their abundance and competitiveness [4]. In fact, in the IEA predictions for the year 2021, solar PV accounts for 54% of the growth in the world's installed net renewable electricity capacity [2]. The remarkable progress that the technology associated with the sensors, recalculating the position of the sun. Finally, intelligent driver systems incorporate artificial intelligence techniques to control the movement of collectors [24].
Another of the determining factors in the performance of a PV plant is the shading of its modules, since these shadows not only imply a lower incident irradiance [25] but also give rise to the appearance of hot spots that bring with them overheating and losses in energy production [26]. In this regard, the behavior of PV modules when they are partially shaded has been widely debated. Several simulation models have been used to find a configuration less susceptible to shadow problems of solar cells [27][28][29][30]. Specifically, Díaz-Dorado et al. [31] have analyzed the effect of shading in a PV tracker with partially shaded astronomical tracking based on the exact arrangement of shaded cells and modules. Other authors have analyzed the energy cost in the production of PV plants as a function of the connections between the cells and the modules [32] of the ground cover ratio (GCR) of the plant that depends on the variables of their design [33] or of the tracking strategy [26,34].
One possible solution to alleviate the effect of shadows in PV plants with solar tracking is back tracking [35], which consists of modifying the orientation of the collecting surfaces in shading situations between panels in order to eliminate such shadows. Another possible solution is to modify the geometry of the collecting surfaces. Thus, although most of the PV panels found in PV plants connected to grids are rectangular, there are already some installations with dual-axis solar tracking where the collectors have other geometric shapes, such as those developed by the Deger Ibérica company in Tarragona (Spain), Ontario (Canada), or Estonia, with 15.6 kWp, 24 MWp, and 100 kWp installed, respectively [36]. However, no previous works have been found in the literature aimed at characterizing the geometric shape of the collectors or their degree of modularity in terms of optimizing the performance of a photovoltaic installation in the event of the possible incidence of shading.
In this context, in the present work, the annual radiative uptake has been quantified in a wide set of PV installations with dual-axis monitoring with different geometries in which the shape of the collectors and the design parameters have been systematically varied. For each case, it has been assumed that the collectors follow the optimal solar tracking strategy proposed by Fernández-Ahumada et al. [37,38]. Unlike traditional solar tracking methods that search for the position of the sun at each moment using astronomical models [39][40][41], according to this strategy, the collecting planes are oriented at each instant of time towards the direction of space in which the irradiance is maximum, except during the moments when such orientation implies the shading of one collector over another. When this occurs, normally at the beginning and end of the day, the collectors are oriented in the direction in which, without causing shading, the incident irradiance on the collectors is maximum. In this way, the production of each installation, which is considered proportional to annual solar radiation, is calculated under the hypothesis that each installation will follow an optimal tracking strategy adapted to its own geometry. Therefore, this study aims to advance the characterization of the electrical behavior of shaded solar trackers, which is an issue where the scientific community has made a considerable effort, simulating tracking strategies for an improvement in photovoltaic production [12,14,42].
To achieve this objective, after this introduction in which the scientific advances made by the scientific community in the field of solar tracking in PV plants are presented, the following section describes the methodology followed in this study to simulate annual solar radiation incident on collectors with a dual-axis tracking strategy that optimizes radiative uptake while avoiding shading between collectors. Similarly, the methodology established for the study of the influence of the design parameters of a PV plant on this annual solar irradiance is explained. Based on this, Section 3 presents the results when applying the methodology described to an existing PV plant ("El Molino", Córdoba, Spain) with its design parameters systematically modified. Similarly, an adjustment model is proposed that represents in a simplified manner the dependence of the annual solar irradiance with respect to the design variables studied, and these dependencies are quantified. Finally, in Section 4, the main conclusions of the present study are presented.

Vector Treatment of the Solar Position and the Estimation of the Solar Irradiance in the Celestial Sphere
In accordance with the above, this work presents a study of the influence of design variables on the performance of a PV plant with dual-axis tracking and optimal tracking strategy. This tracking strategy, which is described in the following section, will determine the orientation of the solar collectors at all times, which allows for optimizing radiative capture while avoiding shadows between collectors. In the present work, to determine at each instant of time the orientation of these collectors, as well as the position of the Sun in the celestial sphere, vector notation is used. Figure 1 shows the coordinate system used, with the Ox axis oriented towards the west, the Oy axis towards the south, and the Oz axis towards the zenith. In this equation, → i , → j , and → k are the unit vectors on the axes Ox, Oy, and Oz, respectively; the unit vector that points towards the Sun, called the solar vector, → s , will be expressed by Equation (1), in which ϕ is the latitude of the place, δ the declination, and t the solar hour.
Sci. 2021, 11, x FOR PEER REVIEW 5 of 24 The orientation of the collector plane at each instant of time can be represented by the unit vector perpendicular to it, ⃗⃗, or by the pair of unit vectors, ⃗⃗ and ⃗, contained in the collector plane, of which ⃗⃗ is horizontal and ⃗ is parallel to the maximum slope direction of the collector plane ( Figure 1). Thus, γ and β being the angles representing the azimuth and inclination of the collectors, respectively, the vectors ⃗⃗, , ⃗⃗⃗⃗ and ⃗ will be given by Equations (4)-(6), respectively.
For the characterization of the incident solar irradiance on the collector planes, the model of Perez et al. [43] has been used. According to this model, the global solar irradiance, I, on an inclined plane at an angle β is the sum of the direct, diffuse, and reflected The orientation of the collector plane at each instant of time can be represented by the unit vector perpendicular to it, in the collector plane, of which → u is horizontal and → v is parallel to the maximum slope direction of the collector plane ( Figure 1). Thus, γ and β being the angles representing the azimuth and inclination of the collectors, respectively, the vectors → n , → u, and → v will be given by Equations (4)-(6), respectively.
For the characterization of the incident solar irradiance on the collector planes, the model of Perez et al. [43] has been used. According to this model, the global solar irradiance, I, on an inclined plane at an angle β is the sum of the direct, diffuse, and reflected irradiance on an inclined plane, which, in turn, depends on the direct irradiance, I B , and diffuse, I D , on a horizontal surface, according to Equation (7), where θ is the angle between the normal vector and the collector plane → n and the solar vector → s , θ z is the zenith angle, ρ is the albedo of the surface of incidence of the irradiance before being reflected towards the collector, a and b are parameters given by Equations (8) and (9), and F 1 and F 2 are the weighting factors for the decomposition of the inclined diffuse radiation, the second addition, in the three subcomponents considered by the authors: isotropic diffuse, circumsolar diffuse, and diffuse from the horizon (Perez et al., 1990).
Substituting the vector expressions (10)- (12) in Equation (7), we obtain Equation (13), in which the dependence of I on → n is made explicit, The production of a PV plant will be greater to the extent that the incident irradiance on its collectors is greater. In this regard, Fernández-Ahumada et al. [37] showed that the incident irradiance on the collectors is maximum when they are oriented in such a way that their normal vector verifies Equation (14). In accordance with this result, the authors propose a new tracking strategy for PV plants with dual-axis trackers in which they are oriented according to Equation (14), and consequently, radio capture and energy production are optimized.
Appl. Sci. 2021, 11, 3917 6 of 23 Therefore, when applying this tracking strategy to the model of Perez et al. (1990), the partial derivatives given by expressions (15) and (16) are obtained. Substituting these results in Equation (14), Equation (17) is obtained, which represents the direction of the normal vector to the collecting planes corresponding to the optimal tracking strategy proposed by Fernández-Ahumada et al. [37] for PV plants with dual-axis tracking.
Given the difficulty of solving the variable consequently, → n will be given by Equation (19). Given the nature of the problem, in this work, it has been considered that j = 25 is a sufficiently high value since the correct convergence of this value has been verified for all the practical cases studied.

Method to Avoid Inter-Shading of Collectors
According to the studies found in the bibliographic review, the shading of solar panels negatively affects their energy production. Therefore, to optimize the performance of a PV plant with dual-axis solar trackers, as those discussed in this study, it is necessary to accompany an adequate tracking strategy with a procedure that prevents the collectors from shading each other. In this regard, Fernández-Ahumada et al. [38] presented a dichotomous criterion to quickly determine whether or not there is inter-shading between collectors. The method is based on the fact that, in an installation with collectors of the same geometric shape oriented towards the same direction of the celestial sphere (characterized by its normal vector, → n ), the shadow cast by a collector Π i on the plane ψ that contains the reference collector Π 0 is a polygon Π i with the same shape as the contour of the collector (Figure 2).

Method to Avoid Inter-Shading of Collectors
According to the studies found in the bibliographic review, the shading of solar panels negatively affects their energy production. Therefore, to optimize the performance of a PV plant with dual-axis solar trackers, as those discussed in this study, it is necessary to accompany an adequate tracking strategy with a procedure that prevents the collectors from shading each other.
In this regard, Fernández-Ahumada et al. [38] presented a dichotomous criterion to quickly determine whether or not there is inter-shading between collectors. The method is based on the fact that, in an installation with collectors of the same geometric shape oriented towards the same direction of the celestial sphere (characterized by its normal vector, ⃗⃗), the shadow cast by a collector Π on the plane ψ that contains the reference collector Π 0 is a polygon Π′ with the same shape as the contour of the collector ( Figure  2). Hence, vector ⃗ is the geometric displacement vector from polygon Π 0 to Π′ . This vector, as shown by the authors, is obtained through the Equation (20), although it must Hence, vector → d i is the geometric displacement vector from polygon Π 0 to Π i . This vector, as shown by the authors, is obtained through the Equation (20), although it must be considered that, for the collector Π i to be able to shade the collector Π 0 , the conditions expressed in Equations (21)-(23) must be fulfilled.
From expression (20), given that → d i is a vector included in the collector plane ψ, Equations (24) and (25) allow for determination of the components of this vector in the OXY reference system contained in the plane ψ and defined by the vectors Knowing → d i and its components d Xi and d Yi , the dichotomous criterion proposed by Fernández-Ahumada et al. [38], based on Minkowski's algebra [44][45][46] and whose validity is demonstrated in Appendix A, it is established that there will be an intersection of the polygons Π 0 and Π i . Therefore, mutual shading between the panels Π 0 y Π i , if, when representing vector → d i from the origin of coordinates, the end of the vector is included within the plane curve Σ, which is obtained as the envelope of the family of all polygons that can be drawn on the plane ψ by translating the polygon Π 0 with the condition that its perimeter is in contact with the coordinate origin of the reference system contained in ψ ( Figure 3).
OXY reference system contained in the plane ψ and defined by the vectors ⃗⃗ y ⃗.
Knowing ⃗ and its components and , the dichotomous criterion proposed by Fernández-Ahumada et al. [38], based on Minkowski's algebra [44][45][46] and whose validity is demonstrated in Appendix A, it is established that there will be an intersection of the polygons Π 0 and Π′ . Therefore, mutual shading between the panels Π 0 y Π , if, when representing vector ⃗ from the origin of coordinates, the end of the vector is included within the plane curve Σ, which is obtained as the envelope of the family of all polygons that can be drawn on the plane ψ by translating the polygon Π 0 with the condition that its perimeter is in contact with the coordinate origin of the reference system contained in ψ ( Figure 3). It is important to note that, since photovoltaic collectors are built by annexing rectangular photovoltaic modules, the resulting geometric shape for the collector is, in general, a closed polygon in which all sides of the perimeter are contained in two directions perpendicular to each other, as shown by the orange line in Figure 4. This fact leads to the determination of the envelope Σ being simplified and able to be performed analytically. Furthermore, this envelope (blue trace in Figure 4) will also be a closed polygon with the It is important to note that, since photovoltaic collectors are built by annexing rectangular photovoltaic modules, the resulting geometric shape for the collector is, in general, a closed polygon in which all sides of the perimeter are contained in two directions perpendicular to each other, as shown by the orange line in Figure 4. This fact leads to the determination of the envelope Σ being simplified and able to be performed analytically. Furthermore, this envelope (blue trace in Figure 4) will also be a closed polygon with the sides parallel and perpendicular to the directions of the perimeter of the collectors ( Figure 4). Therefore, the envelope is mathematically described by the coordinates X env (k) and Y env (k) of each vertex k of the polygon (1 < k < N env , N env being the total number of vertices of Σ). In Appendix B, the analytical obtaining of the enclosure for photovoltaic collectors is detailed in a generic way.
According to Figure 4, the end of vector → d i would be located inside the envelope Σ and, therefore, there would be shading between the panels Π 0 y Π i if there is some vertex k for which the mathematical condition given by Equation (26) is fulfilled.
According to Figure 4, the end of vector ⃗ would be located inside the envelope and, therefore, there would be shading between the panels Π 0 y Π if there is some v tex k for which the mathematical condition given by Equation (26) is fulfilled.
(2 In accordance with the above, Figure 5 shows, by means of a flow diagram, the int shading (⃗⃗) function, especially designed to implement the method described in this s tion to determine the existence of shadow between modules in PV installations with du axis tracking and regular distribution of collectors. This Boolean function depends on t vector ⃗⃗⃗⃗, Equation (19), returning a value "true" if ⃗⃗ implies inter-shading and "fals otherwise. As can be seen in Figure 5, only possible shading is studied in the referen collector, considered representative of the set, which is totally surrounded by tracke In accordance with the above, Figure 5 shows, by means of a flow diagram, the intershading → n function, especially designed to implement the method described in this section to determine the existence of shadow between modules in PV installations with dual-axis tracking and regular distribution of collectors. This Boolean function depends on the vector → n , Equation (19), returning a value "true" if → n implies inter-shading and "false" otherwise. As can be seen in Figure 5, only possible shading is studied in the reference collector, considered representative of the set, which is totally surrounded by N t trackers. In this way, in the case that the value of ⃗⃗, obtained from Equation (19), does not imply shadows between collectors, ℎ (⃗⃗) = , the optimal tracking strategy will propose the orientation of the collectors according to the address of ⃗⃗. On the contrary, if the value of ⃗⃗, obtained from Equation (19), implies shadows between collectors, ℎ (⃗⃗) = , the orientation that does not imply shadows and that max- In this way, in the case that the value of → n , obtained from Equation (19), does not imply shadows between collectors, intershadows → n = False, the optimal tracking strategy will propose the orientation of the collectors according to the address of → n . On the contrary, if the value of → n , obtained from Equation (19), implies shadows between collectors, intershadows → n = True, the orientation that does not imply shadows and that maximizes the incident irradiance must be sought. In that case, back tracking is used following the process indicated in Figure 6.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 11 of 24 Figure 6. Flow diagram to determine the orientation n ⃗⃗⃗⃗of maximum irradiance uptake without any inter-shading between collectors.

Calculation Scheme of Intercepted Solar Radiation
With the methodology described for solar tracking, the incident solar irradiance on solar collectors has been simulated for the 12 representative days, according to Klein [47]. From this irradiance, the incident radiation on the collectors on each representative day m has been calculated using Equation (27), where the integral has been approximated by discretizing the sum in time intervals of three minutes.

Calculation Scheme of Intercepted Solar Radiation
With the methodology described for solar tracking, the incident solar irradiance on solar collectors has been simulated for the 12 representative days, according to Klein [47]. From this irradiance, the incident radiation on the collectors on each representative day m has been calculated using Equation (27), where the integral has been approximated by discretizing the sum in time intervals of three minutes.
Once the daily radiation has been calculated for the representative days proposed by Klein [47] for each of the 12 months of the year (m = 1,2, . . . , 12), the annual global radiation H year is calculated according to Equation (28), where N m is the number of days in the month m.

Cases Analyzed
In this work, the effect of the shape of the collectors on the capture of annual radiation in a PV plant with dual-axis tracking has been studied. To generate the study scenarios, the design of "El Molino", a photovoltaic installation located in Córdoba

Cases Analyzed
In this work, the effect of the shape of the collectors on the capture of annual radiation in a PV plant with dual-axis tracking has been studied. To generate the study scenarios, the design of "El Molino", a photovoltaic installation located in Córdoba (latitude = 37.75492° N; longitude = 5.04548° W) was used as a starting point. It is an installation with dual-axis trackers (with azimuth and elevation movement) and rectangular collectors 8 m wide and 5 m high. Each collector is made up of 25 photovoltaic modules 1m high and 1.6m wide (Figure 7a). The trackers are arranged in a regular grid on horizontal ground, initially separated by a distance =20 m in the EW direction and =14 m in the NS direction (Figure 7b). Modification of the collector shape. The possibility of introducing cuts at the vertices has been considered to study the collectors with the shape indicated in Figure 8. Letting • Modification of the collector shape. The possibility of introducing cuts at the vertices has been considered to study the collectors with the shape indicated in Figure 8.  The crossing of all the possibilities generated = · · = 25000 different combinations of geometric designs for the PV plant "El Molino" with dual-axis tracking and its collectors. For each of these designs, the annual incident radiation on the collectors was obtained by the method set forth, Equation (28), considering the 12 characteristic days proposed by Klein [47], and the monthly mean radiation values [48] set forth in Table 1.   The crossing of all the possibilities generated = · · = 25000 different combinations of geometric designs for the PV plant "El Molino" with dual-axis tracking and its collectors. For each of these designs, the annual incident radiation on the collectors was obtained by the method set forth, Equation (28), considering the 12 characteristic days proposed by Klein [47], and the monthly mean radiation values [48] set forth in Table 1.  The crossing of all the possibilities generated N T = N G ·N D ·N S = 25,000 different combinations of geometric designs for the PV plant "El Molino" with dual-axis tracking and its collectors. For each of these designs, the annual incident radiation on the collectors was obtained by the method set forth, Equation (28), considering the 12 characteristic days proposed by Klein [47], and the monthly mean radiation values [48] set forth in Table 1.

Results
This section describes the incident radiation values on collectors and the results of the study of the influence of the design variables on the solar incidence on the PV plant designs considered. The synoptic values obtained for the annual incident radiation in collectors are summarized in Table 2. Figure 10 shows the distribution of values obtained depending on the membership intervals. In the set of values obtained, it is observed that certain forms of the collector offer the same value. These are the cases in which the set of cuts A, B, C and D of Figure 11a are permuted as shown in Figure 11b-d,f-i. A detailed analysis of the procedure followed makes it possible to verify that the shape of the solar collector does not influence directly but rather through its envelope. As a consequence, the collector shapes of the first column (11 (a), 11 (b), 11 (c), 11 (d)) give rise to identical annual radiation results since they all give rise to the same envelope Σ (represented in Figure 11e). On the other hand, the coincidence of radiation outcomes in the results between the collector shapes of the first (Figure 11a-d) and the second column (Figure 11f-i) should be understood as a consequence of the symmetry with respect to the NS axis of the studied configurations and of the symmetry with respect to this plane in the positions of the sun between the hours of the morning and afternoon.  In the set of values obtained, it is observed that certain forms of the collector offer same value. These are the cases in which the set of cuts A, B, C and D of Figure 11a permuted as shown in Figure 11b-d,f-i. A detailed analysis of the procedure follow makes it possible to verify that the shape of the solar collector does not influence direc but rather through its envelope. As a consequence, the collector shapes of the first colum (11 (a), 11 (b), 11 (c), 11 (d)) give rise to identical annual radiation results since they give rise to the same envelope Σ (represented in Figure 11e). On the other hand, the co cidence of radiation outcomes in the results between the collector shapes of the first (F ure 11a-d) and the second column (Figure 11f-i) should be understood as a consequen of the symmetry with respect to the NS axis of the studied configurations and of the sy metry with respect to this plane in the positions of the sun between the hours of the mo ing and afternoon. The dependence of the values obtained with respect to the considered design v bles was studied using an approximation and simplified model, in which the calcu annual radiation was expressed as a function of the following explanatory variables  Collector surface ( 2 ).  Distance between trackers in east-west direction, ( )  Distance between trackers in north-south direction, ( ) Figure 11. Set of collector shapes that generate identical annual incident radiation results. The forms (a-d) generate the same envelope Σ (e) and the set of symmetric figures with respect to a vertical axis (f-i) generate a symmetric envelope Σ' of Σ (j).
The dependence of the values obtained with respect to the considered design variables was studied using an approximation and simplified model, in which the calculated annual radiation was expressed as a function of the following explanatory variables:

•
Collector surface S col m 2 . • Distance between trackers in east-west direction, D EW (m) • Distance between trackers in north-south direction, D NS (m) • Discriminatory variable of the type of configuration T (T = 1 for staggered configurations and T = 0 for regular grids).
In this sense, it is worth highlighting that the model described and used (for each of the 25,000 cases) can be considered a mathematical function of the variables proposed. However, the complexity of the model and the need to aggregate results on the different representative days made it difficult to know the weight or influence of each variable on the final results. Thus, to overcome this difficulty, this paper proposes to replace this complex function with a mathematical function of simple expression reproducing the result of the complex model with the least possible error. The reader should assume that this is not a statistical problem but rather a problem of adjustment or approximation of a simple expression function to a complex function, so that statistical methods are not applicable. To address the fit, the set of simple variables was extended with composite variables obtained as products and ratios of simple variables. The proposed function (29) was selected from the set of fits to linear functions of composite variables. It lacks a clear physical meaning, but it allows for reproduction of the results of the model with an average relative error equal to 2.561 kWh/m 2 year; therefore, it is considered suitable for the study of the relative weight of the variables.
Equation (29) shows the mathematical expression of this model for which the parameters a, b, c, d, e, f, g, h, l, r, and w have been obtained by the least squares method (Table 3), with an adjustment coefficient R 2 = 0.993. Table 4 shows the synoptic values of the estimation errors of the model ε and ε rel , given by Equations (30) and (31), where H est year is the annual solar irradiance on the solar collectors estimated according to Equation (28) and H adj year is the one approximated by the model (29). The low value obtained for the mean square error means that the equation can be considered valid for the study of dependence of the annual solar irradiance with respect to the variables S col , D EW , D NS , and T.   It is important to note that the model adjusted in Equation (30) does not consider the geometric shape of the collectors itself since, given that the mean error, ε, of the proposed model is 3.3 kWh/m 2 year (Table 4), the geometric shape would not have an explanatory capacity superior to this ε. This reasoning allowed us to conclude the little influence of the shape of the collectors in facilities that follow the tracking/ back-tracking policy considered in this work.
For a better interpretation of the adjusted Equation (29), we can consider it as the addition of four terms separated by parentheses in Equation (32).
The first term only depends on the collector surface S col . Given that b < 0, it was found that, as the collectors were larger, the lower the annual incident radiation. This effect was due to the greater possibility of inter-shading as the collectors had more surface area. With the cases studied, the variation interval of the final result due exclusively to this term was 24 kWh/m 2 year.
The second term marks the importance of the geometric design of the plant given by D EW and D NS , regardless of whether it is a staggered or grid configuration. In the group of cases studied, the variation interval was 116 kWh/m 2 year.
The third term is interpreted as a function of the ground cover ratio (GCR) parameter defined by Equation (33).
According to this definition, the third term considered in the model given by Equation (29) can be rewritten obtaining the expression (34). Thus, this term showed that the variation in this term was 104 kWh/m 2 year for the GCR values considered in the set of cases studied.
Finally, given that the parameter w was negative, the term {wT} implied a small difference of 1.6 kWh/m 2 year to the detriment of the installations that were arranged in a staggered pattern compared to those with a regular grid.

Conclusions
This work presents a novel methodology for the productive study of PV solar collectors mounted on dual-axis trackers. The study, applied to multiple cases, generated as variations with respect to the design adopted in an existing PV Plant ("El Molino", located in Córdoba) allows for identification of the design variables that fundamentally influence the annual incident irradiation on the solar collectors and, therefore, on the energy production of the PV plant. Specifically, by systematically varying the geometry of the collectors, the distance between them, and their spatial distribution, 25,000 case studies were simulated. For all of them, the annual incident solar radiation on the solar collectors was calculated, using for this the irradiance estimation model of Perez [43] and assuming that they were governed by a tracking strategy that optimized radiative capture while avoiding intershading between collectors [37,38]. Although a comprehensive methodology was used to study the case, it was difficult with the data set to understand the influence of each variable on the final result. Therefore, a simple function was fitted, which accurately reproduced the result of the complex model ε_rel = 0.00115. From the irradiance data obtained for the different PV plant designs, a simple mathematical model has been obtained, Equation (29) with a high level of adjustment R 2 = 0.993 that represents the dependence of the annual solar irradiance on the PV plant with respect to design variables such as collector surface, S col , NS and EW distances between collectors D NS and D EW , and spatial distribution of the collectors: regular or staggered grid. Thus, the proposed model has made it possible to identify that the main variables that influence the annual incidence of irradiance are, in order of greatest influence, the geometric design of the plant as a function of the distances between its collectors, (D NS , D EW ), the GCR of the installation and the surface of its collectors (S col ), which can lead to variations in the electrical production of the PV plant of up to 116 kWh/m 2 year, 104 kWh/m 2 year, and 24 kWh/m 2 year. However, for practical purposes, with regard to the spatial distribution of the collectors, it is not appropriate to assume a better behaviour of the regular grid arrangement with respect to the staggered shape since the margin of 1.6 kWh/m 2 year for this variable falls within the uncertainty margin of the radiation prediction models. Similarly, the geometric shape of the collectors does not exert a significant influence on the irradiance received since it gives rise to variations of 3.3026 kWh/m 2 year, which, therefore, are lower than the uncertainty margin of the estimation model itself. Future works will study the influence of the shape in confluence with the orientation of the terrain.   To determine if there is inter-shading between collectors, in this work, a dichotomous criterion is used, based on Minkowski's algebra [44][45][46], according to which there is intershading between two collectors Π 0 y Π i if and only if the displacement vector → d i from collector Π 0 to Π i , which is the projection of collector Π i on the plane ψ that contains Π 0 and can be considered a translation of Π 0 , is included within the plane curve Σ, which is obtained as the envelope of the family of all polygons that can be drawn in the plane ψ when translating the polygon Π 0 , with the condition that its perimeter is in contact with the origin of the reference system contained in ψ.
Thus, this appendix demonstrates this dichotomous criterion for two generic surfaces Π and its translation Π T .
"A polygon Π and its translation Π T intersect if and only if the translation vector The demonstration has two steps: Since → OP is included in Σ ( Figure A1a), there will be at least one polygon Π that, passing through the origin of coordinates O, encloses point P ( Figure A1b). If this polygon is translated following → OP, it is evident that point O is translated to P ( Figure A1c). Since this point P, the boundary of Π T , is interior to Π, there will be an intersection between Π and Π T . shading between two collectors 0 y if and only if the displacement vector ⃗ from collector 0 to Π′ i , which is the projection of collector on the plane that contains 0 and can be considered a translation of 0 , is included within the plane curve , which is obtained as the envelope of the family of all polygons that can be drawn in the plane when translating the polygon 0 , with the condition that its perimeter is in contact with the origin of the reference system contained in .
Thus, this appendix demonstrates this dichotomous criterion for two generic surfaces: and its translation .
"A polygon and its translation intersect if and only if the translation vector ⃗⃗⃗⃗⃗⃗⃗ of with respect to is included in the envelope ". The demonstration has two steps:


Step 1: If ⃗⃗⃗⃗⃗⃗⃗ is included in , and intersect.
Since ⃗⃗⃗⃗⃗⃗⃗ is included in (Figure A1a), there will be at least one polygon that, passing through the origin of coordinates O, encloses point P ( Figure A1b). If this polygon is translated following ⃗⃗⃗⃗⃗⃗ , it is evident that point O is translated to P ( Figure A1c). Since this point P, the boundary of , is interior to , there will be an intersection between and . Figure A1. Demonstration sequence: If ⃗⃗⃗⃗⃗⃗ is included in , and intersect.


Step 2: If and intersect, ⃗⃗⃗⃗⃗⃗ is included in .
If and intersect, there will exist, at least. a point P, of the perimeter of that is included in . By calling O the origin of P, it is possible to define the vector ⃗⃗⃗⃗⃗⃗⃗ ( Figure  A2a). If the system of coordinate axes passing through O is considered and the envelope is drawn, it is evident that the point P is included in (Figure A2b). OP is included in Σ, Π and Π T intersect.
If Π and Π T intersect, there will exist, at least. a point P, of the perimeter of Π T that is included in Π. By calling O the origin of P, it is possible to define the vector → OP ( Figure A2a). If the system of coordinate axes passing through O is considered and the envelope Σ is drawn, it is evident that the point P is included in Σ ( Figure A2b).


Step 2: If and intersect, ⃗⃗⃗⃗⃗⃗ is included in .
If and intersect, there will exist, at least. a point P, of the perimeter of that is included in . By calling O the origin of P, it is possible to define the vector ⃗⃗⃗⃗⃗⃗⃗ ( Figure  A2a). If the system of coordinate axes passing through O is considered and the envelope is drawn, it is evident that the point P is included in (Figure A2b).

Appendix B. Obtaining the Generic Form σ for Collectors with Cuts
In this appendix, the procedure for obtaining the flat curve Σ associated with a collecting surface Π with a generic shape is explained in detail so that it is valid for all the collector shapes considered in this work.
Letting a plane ψ and a reference system OXY be in this plane, we considered a polygon Π contained in ψ with the shape of a rectangle with sides L x and L y with the vertices trimmed, the coordinates of these vertices being: (X ul , Y ul ), (X ur , Y ur ), (X dl , Y dl ), and (X dr , Y dr ) ( Figure A3).

Appendix B. Obtaining the Generic Form σ for Collectors with Cuts
In this appendix, the procedure for obtaining the flat curve Σ associated with a collecting surface Π with a generic shape is explained in detail so that it is valid for all the collector shapes considered in this work.
Letting a plane and a reference system OXY be in this plane, we considered a polygon Π contained in with the shape of a rectangle with sides and with the vertices trimmed, the coordinates of these vertices being: ( , ), ( , ), ( , ), and ( , ) ( Figure A3). The analysis of all the translations allows us to obtain the coordinates of the perimeter of Σ . In general, Σ is a polygon with 28 vertices that depend on the dimensions reflected in Figure A3 according to expression (A1).
In this way, these coordinates can be structured in two arrays of 28 elements: Xenv (28) and Yenv (28). Next, the VBA programming code that allows obtaining each of these coordinates is shown. The analysis of all the translations allows us to obtain the coordinates of the perimeter of Σ. In general, Σ is a polygon with 28 vertices that depend on the dimensions reflected in Figure A3 according to expression (A1). Σ = Σ L x , L y , X ul , Y ul , X ur , Y ur , X dl , Y dl , X dr , Y dr (A1) In this way, these coordinates can be structured in two arrays of 28 elements Xenv (28) and Yenv (28). Next, the VBA programming code that allows obtaining each of these coordinates is shown.
Xenv (1)  However, when there is some kind of symmetry or when any of the dimensions of the cutouts is null, some coordinates contained in the Xenv () and Yenv () arrays will be repeated. In these cases, the number of vertices of Σ will be less than 28. After eliminating the possible repeated points, it is mathematically verified that there are four combinations or symmetries that generate the same polygon Σ Equation (A2).
Σ L x , L y , X ul , Y ul , X ur , Y ur , X dl , Y dl , X dr , Y dr = Σ L x , L y , X dr , Y dr , X ur , Y ur , X dl , Y dl , X ul , Y ul = Σ L x , L y , X ul , Y ul , X dl , Y dl , X ur , Y ur , X dr , Y dr = Σ L x , L y , X dr , Y dr , X ur , Y ur , X dl , Y dl , X ul , Y ul (A2) In this way, it is possible to obtain the envelope Σ to the family of polygons Π on which the dichotomous criterion of inter-shading between collectors will be based.