The Effect of Collector Azimuth on Inter-Row Shading in Photovoltaic Fields—A Comprehensive Point of View

: Rooftop buildings, mountain terrains, rivers, and lake beds may be oriented at different azimuths other than due south (in the northern hemisphere), on which photovoltaic (PV) systems may be installed. The present study presents a comprehensive point of view on the relation between the inter-row spacing and the inter-row shading in multiple-row collectors of PV ﬁelds; it analytically formulates the inter-row spacing and the associated shading losses for collectors deployed at any azimuth on horizontal, sloped-planed, and oriented toward the south and north, and is supported by the numerical results. Several criteria for the inter-row spacing are applied in the study; among them is a newly introduced criterion appropriate for the collectors deployed at a given azimuth. This study indicates that shorter spacing results in increased shading; minimum shading is obtained for collectors oriented toward the equator; inter-row spacing for maintenance purposes is used in limited dimensions of the PV ﬁelds and, usually, in optimal PV system designs; less shading occurs for systems deployed on sloped planes oriented toward the south; less shading is obtained for collectors based on the new criterion, where collectors are deployed at azimuth equal to the solar azimuth, and larger spacing is required for collectors deployed on northern slopes.


Introduction
The inter-row spacing in multiple-row deployments of the PV collectors affects the shading losses and, hence, affects the performance and economics of the photovoltaic (PV) systems.The inter-row shading reduces the generated electric energy of the system, and the means to mitigate the shading losses is an objective for PV system designs.Increasing the spacing results in a reduction in these losses, however, at the expense of the available land area.The inter-row spacing depends on the dimension of the available PV plant land, limited or unlimited.For limited land (e.g., rooftops), the PV collector rows are usually dense-deployed, with smaller inter-row spacing, as the objective of the PV system design is to meet some desired electric output power as an example, thus settling for a larger percentage of shading losses.Dense deployment of the collector rows may also result from using constraint optimization methods for limited land, depending on the objective of the system design, whether energetic or economic [1][2][3][4][5].The optimization process results, among others, in optimal inter-row spacing and optimal inclination of the collectors.The common approaches to determine the inter-row spacing are as follows: minimum required distance for maintenance purposes of the PV system and no shading on collectors on the winter solstice at solar noon or at different times.A few references mention the inter-row spacing in the topic of the articles [6,7], although shading in PV fields is dealt with in many very early studies [8][9][10][11].Analytical expressions were developed in [12] for the spacing between PV rows on horizontal and sloping grounds.Minimum spacing between the rows to avoid shadings is presented in [13].A vector analysis method for the row spacing in PV systems on horizontal and non-horizontal planes is presented in [14].Most studies on inter-row spacing and shading of PV systems deal with the deployment of the PV collectors on horizontal and inclined planes, including flat rooftops, that are oriented toward the equator.The effect of the collector azimuth on the energy yield was analyzed on monitored installed systems [15].Mountain terrains, riverbeds, rooftop buildings, lakes, and ponds may have slopes oriented at any azimuth on which a PV system may be installed.An analytical approach to investigate the effect of the collector azimuth on the inter-row spacing and on the inter-row shading in photovoltaic fields was not investigated thoroughly in the previous studies.The present study investigates the effect of the azimuth angle on the inter-row spacing and shading for four criteria: maintenance spacing, spacing for no shading on solstice day, spacing for optimal deployment of a PV system, and spacing for newly introduced criteria-no shading on collectors at collector azimuth equal to solar azimuth.The new criteria may result in a reduction in shading losses.The effect of collector azimuth on inter-row spacing and shading in photovoltaic fields deployed on horizontal and inclined planes oriented toward the south and the north are investigated.

Methods and Materials
The angles involved in solar calculations of multiple collector rows of PV systems are indicated in Figure 1, where β is the collector inclination angle with respect to horizontal; γ c , γ s are the collector and sun azimuth angles with respect to south, and α is the sun altitude angle.The calculations of the shadow height H s and length L s on the second and on the subsequent collector rows, deployed on a horizontal plane and oriented at any azimuth angle γ c , are based on [16], respectively: L s = L − (D + H cos β) sin β|sin(γ s − γ c )|/ tan α cos β + sin β cos(γ s − γ c )/ tan α (2) where H and L are the collector width and length, respectively, and D is the inter-row spacing.
Energies 2023, 16, x FOR PEER REVIEW 3 of 17 The sun altitude angle  is related to the collector site latitude  , to the solar dec- lination angle  , and to the solar time angle  by [9]: sin sin sin cos cos cos In winter solstice at solar noon, Four spacing criteria are applied to estimate the inter-row spacing and shading: spacing for maintenance purposes; spacing for "no shading" on winter solstice day at solar noon; spacing based on optimal deployment of PV systems in a solar field; and spacing based on new criteria for "no shading" for γ s = γ c .The shadow height depends on the parameters H, D, β, γ c , γ s , α.No shading on collectors; H s = 0 (see Equation ( 1)) occurs for collectors oriented toward the equator (γ c = 0 • ) at solar noon (γ s = 0).Equation (1) then reduces to The sun altitude angle α is related to the collector site latitude φ, to the solar declination angle δ, and to the solar time angle ω by [9]: In winter solstice at solar noon, δ = δ 0 = −23.45• , and ω = 0; Equation (4) becomes By substituting Equation ( 5) into Equation (3) the inter-row spacing is obtained.
The percentage of the annual shading losses of the PV system is defined by the difference between the direct beam radiation of the first and the second collector rows, relative to the global radiation on the first collector row: The shading losses in this study are based on hourly solar radiation data (direct beam and diffuse radiation, Israel, Meteorological Service-IMS) for Tel Aviv; latitude φ = 32 • 6 N; longitude 34 • 51 E. Four criteria are examined for the effect of the azimuth angle on the inter-row spacing and shading: maintenance spacing; spacing for no shading on solstice day; spacing for optimal deployment of a PV system; and spacing for no shading on collectors at collector azimuth equal to solar azimuth.This study deals with collectors deployed on horizontal and sloped planes.
To indicate the effect of the azimuth angle on the inter-row spacing and shading in a totally different geographical location, this study investigated the deployment of collectors at Lindenberg, Germany; latitude φ = 52 • 12 N; longitude 14 • 7 E, based on solar radiation data from Meteorologisches Observatorium Lindenberg, Richard Aßmann Observatorium, PANGAEA, Monitoring station no.12, https://dataportals.pangaea.de/bsrn/?q=LR1300(accessed on 1 January 2019); ten minutes of solar radiation data, the year 2006.The assumed lower limit for the inter-row spacing, in this study, is dictated by the maintenance requirement of D = 0.8 m.The assumed range for the collector azimuth angles γ c is from 0 • to 85 • , relative to the south for latitude φ = 32 • 6 N and from 0 • to 70 • for latitude φ = 52 • 12 N.For comparison, the assumed collector inclination angle for both latitudes is β = 20 • .The value of the upper limit of the inter-row spacing D depends on the differently applied criteria for the inter-row spacing.

Collectors Deployed on Horizontal Planes
This section deals with the inter-row spacing and shading of the PV collectors deployed on horizontal planes for different criteria for the inter-row pacing at different collector azimuth angles.

Maintenance Criteria for Inter-Row Spacing
The maintenance inter-row spacing approach is used in limited areas, for example, on rooftops.The PV collector rows are densely deployed, with smaller inter-row spacing to meet some desired electric output power of the available land area, thus settling for a larger percent of shading losses.It is assumed a maintenance spacing of D = 0.8 m.

No-Shading Criteria for Inter-Row Spacing on Winter Solstice at Solar Noon
The spacing D is determined by Equation (6).The shadow height and length are given in Equations ( 1) and ( 2), respectively.The inter-row spacing is D = 1.05 m for H = 2.12 m, β = 20 • ; the declination angle at solar noon is δ 0 = −23.45• at latitude φ = 32 • 6 N.

γ s = γ c Criterion for Inter-Row Spacing on day December 21
Collectors oriented toward the southwest at angle γ c , with respect to the south, are shown in Figure 2, where the direction of γ c is normal to the collector.The criteria for which the inter-row spacing is defined is when γ s = γ c .The sun azimuth angle γ s is related to the sun declination angle δ, to the sun hour angle ω, and to the sun altitude angle α by [1]).
Substituting into Equation ( 9) and solving for D, we obtain For a given H, β, the inter-row spacing depends on the sun altitude angle α at γ s .
sin tan( ) For a given , H  , the inter-row spacing depends on the sun altitude angle  at s  .An optimal design of a Photovoltaics system may be formulated by an objective function and a set of constraints.The optimal design may be based on energy or economic criteria, and the constraints are bounding conditions imposed on the field and on the collector parameters.In this section, the maximum annual incident energy on the collectors is the objective function.
The optimization problem is formulated as follows: Maximize yearly energy E : with respect to X ( , , , An optimal design of a Photovoltaics system may be formulated by an objective function and a set of constraints.The optimal design may be based on energy or economic criteria, and the constraints are bounding conditions imposed on the field and on the collector parameters.In this section, the maximum annual incident energy on the collectors is the objective function. The optimization problem is formulated as follows: Maximize yearly energy E: C(X) with respect to X (D, H, K, β) where K (integer) is the number of the collector rows in the PV field, and W is the field width.

Collectors Deployed on Sloped Planes Facing the South
Multiple collector rows deployed on an inclined plane ε, facing the south, are shown in Figure 3.The second collector is raised by an amount B, relative to the first collector, thus resulting in shading that takes place only on part H 1 of the collector.
Energies 2023, 16, x FOR PEER REVIEW 6 of 17 where K (integer) is the number of the collector rows in the PV field, and W is the field width.

Collectors Deployed on Sloped Planes Facing the South
Multiple collector rows deployed on an inclined plane  , facing the south, are shown in Figure 3.The second collector is raised by an amount B , relative to the first collector, thus resulting in shading that takes place only on part 1 H of the collector.The inter-row spacing for collectors mounted on inclined planes facing the south is developed in [18] where Inserting Equation (18) into Equation (17), we obtain The inter-row spacing for collectors mounted on inclined planes facing the south is developed in [18] where Inserting Equation (18) into Equation (17), we obtain By substituting Equation (19) into Equation ( 16), we, thus, obtain the inter-row spacing for collectors deployed on south facing slope ε: The shadow height H s on the second and on the subsequent collector rows is given in Equation (1), where H is replaced by H 1 (see [16]): γ s = γ c Criterion for Inter-Row Spacing on day December 21 The collector height of the second and subsequent collectors subject to shading is H 1 (see Figure 3)); therefore, based on Equation ( 13), the spacing is given by where H 1 is given in Equation (17).By solving for D south ε , we, thus, obtain: The collector height of the second and subsequent collectors subject to shading is 1 H (see Figure 3)); therefore, based on Equation ( 13), the spacing is given by where 1 H is given in Equation (17).By solving for south D  , we, thus, obtain:

Collectors Deployed on Sloped Planes Facing the North
A deployment of two collectors on an inclined plane  facing north is shown in ) ; therefore, the inter-row spacing, see Equation ( 13), is given by By substituting Equations ( 5) and (25) into Equation ( 24) and solving for D north ε , we obtain γ s = γ c Criterion for Inter-Row Spacing on day December 21 Based on Figure 4, the second and subsequent collectors are subject to shading for height where B = (D north ε + H cos β) tan ε; therefore, the inter-row spacing, see Equation ( 13), is given by where H * is given in Equation ( 27).Solving for D n • rth ε thus obtaining:

Results
The results of the effect of the azimuth angle on the inter-row spacing and shading pertain to the deployment of collectors on horizontal and sloped planes applying four spacing criteria.The deployments are at latitudes φ = 32 • 6 N and φ = 52 • 1 N.

Results
The results of the effect of the azimuth angle on the inter-row spacing and shading pertain to the deployment of collectors on horizontal and sloped planes applying four spacing criteria.The deployments are at latitudes   6) and ( 7)).

No-Shading Criteria for Inter-Row Spacing and Shading on Winter Solstice at Solar Noon,
The spacing D is determined by Equation ( 6).The shadow height and length are given in Equations ( 1) and ( 2), respectively.The inter-row spacing is It is also interesting to compare the shading losses on a monthly basis.Figure 6 shows that the time for "no shading" decreases as the collector azimuth increases.For   6) and ( 7)).It is also interesting to compare the shading losses on a monthly basis.Figure 6 shows that the time for "no shading" decreases as the collector azimuth increases.For γ c = 0 • , no shading occurs for seven months; for γ c = 10 • -for five months; for γ c = 20 • -for four months; for γ c = 30 • -for one month; and for larger azimuths, shading takes place for the entire year.13) is depicted in Figure 7.The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 8 for a constant spacing   6)).

γ s = γ c Criterion for Inter-Row Spacing on day December 21
The results for the γ s = γ c criterion are depicted for latitudes φ = 32 • 6 N and φ = 52 The variation of the sun elevation angle for the criteria γ c = γ s on winter solstice day needed for Equation ( 13) is depicted in Figure 7.   13) is depicted in Figure 7.The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 8 for a constant spacing  The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 8 for a constant spacing D = 1.05 m (in blue-Equation ( 6)) and for the spacing criterion γ c = γ s (in orange-Equation ( 13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 8 reveals that the inter-row spacing may be equally determined by Equation ( 6) or by Equation ( 13 6) and ( 13)).6), ( 7) and ( 13)).6) and ( 13)).6) and ( 13)).6), ( 7) and ( 13)).6), ( 7) and ( 13)).

Latitude
The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 10 for a constant spacing D = 2.836 m (in blue-Equation ( 6)) and for the spacing criterion γ c = γ s (in orange-Equation ( 13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 11 describes the shading losses percentage as a function of the collector azimuth for a constant spacing D = 2.386 m (in blue) and for a spacing criterion γ c = γ s (in orange).Less shading losses may be obtained by applying the criteria γ c = γ s .6) and ( 13)).6), ( 7) and ( 13)).6) and ( 13)).6) and ( 13)).6), ( 7) and ( 13)).

Inter-Row Spacing and Shading Losses in Optimal PV Systems Designs
The optimal design is performed at two locations   6), ( 7) and ( 13)).

Inter-Row Spacing and Shading Losses in Optimal PV Systems Designs
The optimal design is performed at two locations φ = 32 • 6 N (Tel Aviv, Israel) and φ = 52 • 12 N (Lindenberg, Germany).
The inter-row spacing D is bounded by a lower limit of 0.8 m for maintenance purposes.The results for maximum yearly incident energy, the optimal number of the collector rows K, and the optimal collector inclination angle β, for different azimuth angles of the collector γ c are depicted in Figure 12 for the latitude φ = 32 • 6 N. The optimization process tends to a minimum inter-row spacing of 0.8 m for all azimuth angles of the collectors in order to increase the number of collector rows in the given field and, hence, to obtain high incident energies.The relation between the optimal inclination angle β (in orange) and the azimuth angle γ c , as depicted in Figure 12, is of a stepwise shape.Therefore, the relationship between shading losses (in blue) and the parameters β and γ c is complicated.For example, shading losses of 1.91% correspond to β = 20 • , γ c = 40 • ; 2.68% losses correspond to β = 25 • , γ c = 20 • , and 4.40% losses correspond to β = 30 • , γ c = 10 • .

Inter-Row Spacing and Shading Losses in Optimal PV Systems Designs
The optimal design is performed at two locations   14) and ( 15)).

Latitude φ = 52 • 12 N
The results for maximum yearly incident energy, the optimal number of the collector rows K, and the optimal collector inclination angle β, for different azimuth angles of the collector γ c are depicted in Figure 13 for the latitude φ = 52 • 12 N.
The field and collector parameters at Lindenberg are the same as for Tel Aviv: W max = 100 m, L max = 40 m, H max = 2.12 m, 0.8 m ≤ D ≤ 5.0 m, 0.0001 ≤ β ≤ 90 • , 2 ≤ K ≤ 100.The optimization process tends to a minimum inter-row spacing of 0.8 m for all azimuth angles of the collectors in order to increase the number of collector rows in the given field.The field and collector parameters at Lindenberg are the same as for Tel Aviv: . The optimization process tends to a minimum inter-row spacing of 0.8 m for all azimuth angles of the collectors in order to increase the number of collector rows in the given field.14) and ( 15)).for collectors deployed on inclined planes facing the south is shorter than for the collectors deployed on horizontal planes [18]; therefore, the inter-row spacing in this section is determined by the spacing for maintenance purposes, e.g., 0.8 D  m. Figure 14 14) and ( 15)).for collectors deployed on inclined planes facing the south is shorter than for the collectors deployed on horizontal planes [18]; therefore, the inter-row spacing in this section is determined by the spacing for maintenance purposes, e.g., D = 0.8 m. Figure 14 6) and ( 23)).6), ( 7) and ( 23)).for collectors deployed on inclined planes facing north is longer than for collectors deployed on horizontal planers [18]; therefore, the inter-row spacing, in this section, is determined by Equation ( 6) for winter solstice at solar noon,   6), ( 7) and ( 29)).for collectors deployed on inclined planes facing north is longer than for collectors deployed on horizontal planers [18]; therefore, the inter-row spacing, in this section, is determined by Equation ( 6 6), ( 7) and ( 23)).for collectors deployed on inclined planes facing north is longer than for collectors deployed on horizontal planers [18]; therefore, the inter-row spacing, in this section, is determined by Equation (6) for winter solstice at solar noon,   6), ( 7) and ( 29)).6), ( 7) and ( 29)).The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 18 for a constant spacing D = 1.05 m, for spacing criterion γ c = γ s (for ε = 0 • ), and for ε = ±5 • and ε = ±10 • , respectively, for given parameters H = 2.12 m, β = 20 • , L = 40 m, φ = 32 • 6 N, where the (+) pertains to a southern plane and the (−) pertains to a northern plane.The inter-row spacing (Figure 18) and the annual shading percentage (Figure 19) are drawn up to the azimuth γ c = 40 • for the collectors because the required spacing becomes impractical even for smaller azimuths.Smaller inter-row spacing is required for collectors deployed on southern slopes and larger for northern slopes, with respect to horizontal planes.Figure 19 describes the annual shading losses percentage as a function of the collector azimuth γ c for a constant spacing D = 1.05 m (Equation ( 6)), for spacing criterion γ c = γ s (Equation ( 13)), and for ε = ±5 • and ε = ±10 • , respectively.The figure clearly emphasizes the reduction in the shading losses by applying the criteria γ c = γ s , and larger southern slopes +ε result in less shading.Collectors may be deployed on northern slopes with acceptable inter-row spacing on account of increasing shading.Collectors may be deployed on northern slopes with acceptable inter-row spacing on account of increasing shading.6) and ( 7)).

Discussion and Conclusions
The effect of collector azimuth on inter-row spacing and on inter-row shading in photovoltaic fields was not investigated analytically in previous studies, to the best of our knowledge.The present study presents a comprehensive point of view of the relation between the inter-row spacing and the inter-row shading in multiple-row collectors of the PV fields for any azimuthal deployment of the collector rows.PV collectors may be deployed at any azimuth to exploit available structures and land to promote the utilization of renewable energy.Rooftop buildings, mountain terrains, rivers, lakes, and water pond beds may be oriented at different azimuths other than to the equator, on which PV systems Collectors may be deployed on northern slopes with acceptable inter-row spacing on account of increasing shading.6) and ( 7)).

Discussion and Conclusions
The effect of collector azimuth on inter-row spacing and on inter-row shading in photovoltaic fields was not investigated analytically in previous studies, to the best of our knowledge.The present study presents a comprehensive point of view of the relation between the inter-row spacing and the inter-row shading in multiple-row collectors of the PV fields for any azimuthal deployment of the collector rows.PV collectors may be deployed at any azimuth to exploit available structures and land to promote the utilization of renewable energy.Rooftop buildings, mountain terrains, rivers, lakes, and water pond beds may be oriented at different azimuths other than to the equator, on which PV systems  6) and ( 7)).

Discussion and Conclusions
The effect of collector azimuth on inter-row spacing and on inter-row shading in photovoltaic fields was not investigated analytically in previous studies, to the best of our knowledge.The present study presents a comprehensive point of view of the relation between the inter-row spacing and the inter-row shading in multiple-row collectors of the PV fields for any azimuthal deployment of the collector rows.PV collectors may be deployed at any azimuth to exploit available structures and land to promote the utilization of renewable energy.Rooftop buildings, mountain terrains, rivers, lakes, and water pond beds may be oriented at different azimuths other than to the equator, on which PV systems may be installed.The inter-row spacing plays an important role in the design of PV systems affecting the generated electric energy.The inter-row spacing determines the number of collector rows in the given dimensions of the field and, hence, determines the shading and masking losses and the electric output energy of the PV system.The common approach to determine the inter-row spacing is based on the "winter solstice rule" ("shading angle"), where the inter-row spacing is determined by the shadow length on day December 21 at solar noon.This general rule is used in practice and in theoretical studies of PV system designs.For limited land (e.g., rooftops), the inter-row spacing is usually dictated by the required maintenance distance, and hence, the PV collector rows are usually densely deployed.The approach of optimal deployment of collectors in a given field may lead to the inter-row spacing that is dictated by the objective function and by the constraints on the problem parameters.A new criterion for the inter-row spacing is introduced in the present study, namely, no shading on collectors for collector azimuth equal to solar azimuth; γ c = γ s on the winter solstice.Applying the new criteria decreases the shading losses.The present study deals with the inter-row spacing and inter-row shading in the PV fields as affected by the collector azimuth angle for horizontal and sloped planes oriented toward the south and the north.Generally, shorter spacing results in increased shading; minimum shading is obtained for collectors oriented toward the equator; inter-row spacing for maintenance purposes is used in limited land of the PV fields and, usually, also in optimal PV system designs; less shading occurs for systems deployed on inclined planes oriented toward the equator; less shading is obtained for the collectors facing the south and oriented at azimuth γ c by applying the criterion γ c = γ s , and larger spacing is required for collectors deployed on the northern slopes than on the southern slopes.
Conclusions for the PV collectors of width H = 2.12 m, deployed with an inclination angle of β = 20 • with respect to horizontal at latitude φ = 32 • 6 N (Tel Aviv, Israel), have the following characteristics:

•
For collectors deployed on horizontal planes, the inter-row distance D may be determined nearly equally by the "winter solstice rule" (for example, D = 1.05 m) or by the criterion γ c = γ s , up to collector azimuth γ c = 20 • ; however, the inter-row spacing, based on γ c = γ s criterion, then increases rapidly for γ c > 40 • , resulting in large spacing; however, less shading losses are obtained; • Reduction in shading losses is obtained by applying the new γ c = γ s criterion for collector deployment with an azimuth = 0 • ; see Figure 9, for example;

•
The inter-row spacing for collectors deployed on south-facing slopes is usually determined by the distance required for maintenance; for example, D = 0.8 m;

•
Collectors may be deployed on northern slopes with acceptable spacing on account of increasing shading; • Optimal PV system designs for maximum annual incident solar energy tend to result in minimum inter-row spacing (for example D = 0.8 m-maintenance distance) for all azimuth angles of the collectors in order to increase the number of collector rows in the given solar field, and hence, to obtain high incident energies; • At the latitude φ = 52 • 12 N, it was assumed that an inclination angle β = 20 • , and collector width H = 2.12 m for comparison with the results for φ = 32 • 6 N. A low percentage of shading losses is obtained for γ c = γ s criterion due to large inter-row spacing.

2 year
shading losses [%] = 100 f irst row beam radiation kWh/m 2 year − sec ond row beam radiation kWh/m 2 year f irst row global radiation kWh/m

Figure 3 .
Figure 3. Collectors deployed on inclined plane  , facing the south.

Figure 3 .
Figure 3. Collectors deployed on inclined plane ε facing the south.

2. 3 .
Collectors Deployed on Sloped Planes Facing the North A deployment of two collectors on an inclined plane ε facing north is shown in Figure 4. From the figure, we obtain tan α γ s=0 = (H sin β + B)/D north ε (24) where B = (H cos β + D north ) tan ε (25) Energies 2023, 16, x FOR PEER REVIEW 7 of 17

Figure 4 .
Figure 4. From the figure, we obtain

Figure 4 .
Figure 4. Collectors deployed on inclined plane  facing north.
the annual shading losses percentage with collector azimuth angle c Figure 5 (orange line).

Figure 5 .
Figure 5. Variation of the annual shading losses percentage with collector azimuth c  , variation of the percentage of the shading losses, Equation(7), with the collectors' azimuth angle c  is shown in Figure5for figure shows that the shading losses increase with the increase in the collector azimuth angle c  .Shorter inter-row spacing, 0.8 D  m (in orange) results in higher shading losses.
no shading occurs for seven months; for 10 o c   -for five months; for month; and for larger azimuths, shading takes place for the entire year.
Criteria for Inter-Row Spacing and Shading on Winter Solstice at Solar Noon, φ = 32 • 6 N The spacing D is determined by Equation (6).The shadow height and length are given in Equations (1) and (2), respectively.The inter-row spacing is D = 1.05 m for H = 2.12 m, β = 20 • , declination angle at solar noon δ 0 = −23.45• , at latitude φ = 32 • 6 N. The variation of the percentage of the shading losses, Equation (7), with the collectors' azimuth angle γ c is shown in Figure 5 for D = 1.05 m (blue line) and L = 40 m.The figure shows that the shading losses increase with the increase in the collector azimuth angle γ c .Shorter inter-row spacing, D = 0.8 m (in orange) results in higher shading losses.

Figure 6 .
Figure 6.Monthly shading losses percentage for 2.12 H  m,


Criterion for Inter-Row Spacing on day December 21 The results for the s c of the sun elevation angle for the criteria c s    on winter solstice day needed for Equation (

Figure 7 .
Figure 7. Sun elevation angle on winter solstice for criteria c blue-Equation (6)) and for the spacing criterion c s    (in orange-Equation (13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 8 reveals that the inter-row spacing may be equally determined by Equation (6) or by Equation (13) up to the collector azimuths 20 o c   .The inter-row spacing based on Equation (13) increases rapidly for 40 o c   , resulting in a large spacing.Figure 9 describes the shading losses percentage as a function of the collector azimuth for a constant spacing 1.05 D  m (in blue) and for a spacing criterion c s    (in orange).The figure clearly emphasizes the reduction in the shading losses by applying the criteria c s    .

Figure 6 .
Figure 6.Monthly shading losses percentage for 2.12 H  m,


Criterion for Inter-Row Spacing on day December 21 The results for the s c of the sun elevation angle for the criteria c s    on winter solstice day needed for Equation (

Figure 7 .
Figure 7. Sun elevation angle on winter solstice for criteria c -Equation (13)).The vertical black line indicates the sunset azimuth angle on day December 21.Figure8reveals that the inter-row spacing may be equally determined by Equation (6) or by Equation (13) up to the collector azimuths 20 o c   .The inter-row spacing based on Equation (13) increases rapidly for large spacing.Figure 9 describes the shading losses percentage as a function of the collector azimuth for a constant spacing ).The figure clearly emphasizes the reduction in the shading losses by applying the criteria c s
) up to the collector azimuths γ c = 20 • .The inter-row spacing based on Equation (13) increases rapidly for γ c > 40 • , resulting in a large spacing.Figure 9 describes the shading losses percentage as a function of the collector azimuth for a constant spacing D = 1.05 m (in blue) and for a spacing criterion γ c = γ s (in orange).The figure clearly emphasizes the reduction in the shading losses by applying the criteria γ c = γ s .

Figure 8 .
Figure 8. Variation of the inter-row spacing D with collector azimuth c  for 2.12 H  m,

Figure 9 .
Figure 9. Variation of shading losses percentage with collector azimuth c  for 2.12 H  m, the inter-row spacing with the collector azimuth angle is depicted in Figure10for a constant spacing 2.836 D  m (in blue-Equation (6)) and for the spacing criterion c s    (in orange-Equation (13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 11 describes the shading losses percentage as a function of the collector azimuth for a constant spacing

Energies 2023 , 17 Figure 8 .
Figure 8. Variation of the inter-row spacing D with collector azimuth c  for 2.12 H  m,

Figure 9 .
Figure 9. Variation of shading losses percentage with collector azimuth c  for 2.12 H  m, the inter-row spacing with the collector azimuth angle is depicted in Figure10for a constant spacing 2.836 D  m (in blue-Equation (6)) and for the spacing criterion c s    (in orange-Equation (13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 11 describes the shading losses percentage as a function of the collector azimuth for a constant spacing

Energies 2023 , 17 Figure 8 .
Figure 8. Variation of the inter-row spacing D with collector azimuth c  for 2.12 H  m,

Figure 9 .
Figure 9. Variation of shading losses percentage with collector azimuth c  for 2.12 H  m, the inter-row spacing with the collector azimuth angle is depicted in Figure10for a constant spacing 2.836 D  m (in blue-Equation (6)) and for the spacing criterion c s    (in orange-Equation (13)).The vertical black line indicates the sunset azimuth angle on day December 21. Figure 11 describes the shading losses percentage as a function of the collector azimuth for a constant spacing 2.386 D  m (in blue) and for a spacing criterion c s    (in orange).Less shading losses may be obtained by applying the criteria c s    .

Figure 10 .
Figure 10.Variation of the inter-row spacing D with collector azimuth c  for 2.12 H  m,

Figure 11 .
Figure 11.Variation of shading losses percentage with collector azimuth c  for row spacing D is bounded by a lower limit of 0.8 m for maintenance purposes.The results for maximum yearly incident energy, the optimal number of the collector rows K , and the optimal collector inclination angle  , for different azimuth angles of the collector c  are depicted in Figure12for the latitude optimization process tends to a minimum inter-row spacing of 0.8 m for all azimuth angles of the collectors in order to increase the number of collector rows in the given field and, hence, to obtain high incident energies.The relation between the optimal inclination angle  (in orange) and the azimuth angle c  , as depicted in Figure12, is of a stepwise shape.Therefore, the relationship between shading losses (in blue) and the parameters  and c  is complicated.For example, shading losses of 1.91% correspond to

Figure 11 .
Figure 11.Variation of shading losses percentage with collector azimuth c  for

N
The inter-row spacing D is bounded by a lower limit of 0.8 m for maintenance purposes.The results for maximum yearly incident energy, the optimal number of the collector rows K , and the optimal collector inclination angle  , for different azimuth angles of the collector c  are depicted in Figure12for the latitude optimization process tends to a minimum inter-row spacing of 0.8 m for all azimuth angles of the collectors in order to increase the number of collector rows in the given field and, hence, to obtain high incident energies.The relation between the optimal inclination angle  (in orange) and the azimuth angle c  , as depicted in Figure12, is of a stepwise shape.Therefore, the relationship between shading losses (in blue) and the parameters  and c  is complicated.For example, shading losses of 1.91% correspond to

Figure 12 .
Figure 12.Optimal design parameters for yearly maximum incident energy (in blue) on a PV field as a function of the collector azimuth angle γ c : number of collector rows K, optimal inclination angle β (in orange), φ = 32 • 6 N (Equations (14) and (15)).

Figure 13 .
Figure 13.Optimal design parameters for yearly maximum incident energy (in blue) on a PV field as a function of the collector azimuth angle c  : number of collector rows K , optimal inclination Facing the South 3.5.1.Constant Inter-Row Spacing-Maintenance Spacing, shows the variation of the percentage of the annual shading losses with the collector azimuth c figure shows that a larger inclination angle  results in less shading.

Figure 13 .
Figure 13.Optimal design parameters for yearly maximum incident energy (in blue) on a PV field as a function of the collector azimuth angle γ c : number of collector rows K, optimal inclination angle β (in orange), φ = 52 • 12 N (Equations (14) and (15)).

17 Figure 14 .
Figure 14.Variation of the annual shading losses percentage with collector azimuth c  for the inter-row spacing with the collector azimuth angle is depicted in Figure15for a constant spacing vertical black line indicates the sunset azimuth angle on day December 21.Figure15reveals that the inter-row spacing may be reduced for sloped planes.Figure16describes the annual shading losses percentage as a function of the collector azimuth for a constant spacing 1.05 D  m (Equation (6)), for spacing criterion c s respectively.The figure clearly emphasizes the reduction in the shading losses by applying the criteria c s    , and larger slopes  result in less shading.

Figure 14 .
Figure 14.Variation of the annual shading losses percentage with collector azimuth γ c for H = 2.12 m, β = 20 • , L = 40 m, φ = 32 • 6 N, and for a constant inter-row spacing D = 0.8 m, (Equation (7)).3.5.2.γ s = γ c Criterion for Inter-Row Spacing on day December 21, φ = 32 • 6 N The variation of the inter-row spacing with the collector azimuth angle is depicted in Figure 15 for a constant spacing D = 1.05 m, for spacing criterion γ c = γ s (for ε = 0), and for ε = 5 • and ε = 10 • , respectively for given parameters H = 2.12 m, β = 20 • , L = 40 m, φ = 32 • 6 N. The vertical black line indicates the sunset azimuth angle on day December 21.Figure 15 reveals that the inter-row spacing may be reduced for sloped planes.Figure 16 describes the annual shading losses percentage as a function of the collector azimuth for a constant spacing D = 1.05 m (Equation (6)), for spacing criterion γ c = γ s (Equation (13)), and for ε = 5 • and ε = 10 • , respectively.The figure clearly emphasizes the reduction in the shading losses by applying the criteria γ c = γ s , and larger slopes ε result in less shading.

Figure 15
reveals that the inter-row spacing may be reduced for sloped planes.
Figure 16 describes the annual shading losses percentage as a function of the collector azimuth for a constant spacing D = 1.05 m (Equation (6)), for spacing criterion γ c = γ s (Equation (13)), and for ε = 5 • and ε = 10 • , respectively.The figure clearly emphasizes the reduction in the shading losses by applying the criteria γ c = γ s , and larger slopes ε result in less shading.
vertical black line indicates the sunset azimuth angle on day December 21.Figure15reveals that the inter-row spacing may be reduced for sloped planes.Figure16describes the annual shading losses percentage as a function of the collector azimuth for a constant spacing The figure clearly emphasizes the reduction in the shading losses by applying the criteria c s    , and larger slopes  result in less shading.

Figure 15 .
Figure 15.Variation of the inter-row spacing south D 

Figure 15 .
Figure 15.Variation of the inter-row spacing D south ε

Energies 2023 , 17 Figure 16 .
Figure 16.Variation of the annual shading losses percentage with collector azimuth c  for

3. 6 . 1 .
No-Shading Criteria for Inter-Row Spacing and Shading on Winter Solstice at Solar Noon,

.
17 shows the variation of the percentage of the annual shading losses with the collector azimuth c The figure shows that a larger inclination angle  results in less shading.

Figure 17 .
Figure 17.Variation of the annual shading losses percentage with collector azimuth c  for

3. 6 .
Collectors Deployed on Sloped Planes Facing the North 3.6.1.No-Shading Criteria for Inter-Row Spacing and Shading on Winter Solstice at Solar Noon, φ = 32 • 6 N The inter-row spacing D north ε

17 Figure 16 .
Figure 16.Variation of the annual shading losses percentage with collector azimuth c  for

3. 6 . 1 .
No-Shading Criteria for Inter-Row Spacing and Shading on Winter Solstice at Solar Noon,

.
17 shows the variation of the percentage of the annual shading losses with the collector azimuth c The figure shows that a larger inclination angle  results in less shading.

Figure 17 .
Figure 17.Variation of the annual shading losses percentage with collector azimuth c  for

Energies 2023 ,
16,  x FOR PEER REVIEW 15 of 17 the required spacing becomes impractical even for smaller azimuths.Smaller inter-row spacing is required for collectors deployed on southern slopes and larger for northern slopes, with respect to horizontal planes.Figure19describes the annual shading losses percentage as a function of the collector azimuth c  for a constant spacing , respectively.The figure clearly emphasizes the reduction in the shading losses by applying the criteria c s    , and larger southern slopes   result in less shading.

Figure 18 .
Figure 18.Variation of the inter-row spacing south D 

Figure 19 .
Figure 19.Variation of the annual shading losses percentage with collector azimuth c  for

Figure 18 .
Figure 18.Variation of the inter-row spacing D south ε

Figure 18 .
Figure 18.Variation of the inter-row spacing south D 

Figure 19 .
Figure 19.Variation of the annual shading losses percentage with collector azimuth c  for