A Thermal Model to Estimate PV Electrical Power and Temperature Proﬁle along Panel Thickness

: The production of electricity from photovoltaic panels has experienced signiﬁcant developments. To manage the energy ﬂows introduced into the electricity grid, it is necessary to estimate the productivity of PV panels under the climatic conditions. In this study, a photovoltaic panel is modelled from thermal and electrical points of view to evaluate electrical performance and identify the temperature distribution in the layers. The analysis performed is time dependent and the problem is solved using the ﬁnite difference technique. A methodology is introduced to estimate the cloudiness of the sky, which affects radiative heat exchange. The calculation method is validated using experimental data recorded in a laboratory of the University of Calabria. Temperature and electrical power are predicted with RMSE of 1.5–2.0 ◦ C and NRMSE of 1.2–2.1%, respectively. The evaluation of the temperature proﬁle inside the panel is essential to understand how heat is dissipated. The results show that the top surface (glass) is almost always colder than the back of the panel, despite being exposed to radiation. In addition, the upper surface dissipates more heat power than the lower one. Cooling systems, such as spray cooling, work better if they are installed on the back of the panel.


Introduction
Increasing the proportion of energy from renewable sources is a key element in reducing greenhouse gas emissions and keeping the average global temperature increase below 2 • C compared to pre-industrial levels.In fact, reducing the use of fossil fuels by replacing them with renewable resources contributes to the decarbonization of the energy system.
Energy communities are collective organizations that will help the transition to clean energy.They aim to use energy in the same place where it is produced, sharing it among several citizens and managing storage systems in a smart way.They make it possible to create energy-independent areas, with benefits for citizens and the health of the planet.In this context, photovoltaic (PV) fields play an important role.Predicting their electrical productivity is useful for managing their interconnection within a smart grid [1].In order to optimize the operation of energy communities, it is very useful to be able to predict the energy demand of users, the electrical supplies of the generators and the charge levels of the electrical accumulators installed in the area [2][3][4].Being able to simulate the thermal and electrical behavior of generation systems as a function of atmospheric conditions becomes important, especially if the primary source has variable availability, as in the case of solar source [5].
The estimation of producibility helps to understand whether the system solution is appropriate for the location or whether other types of installation are preferable [6], even to obtain only small advantages.Photovoltaic production can be increased in several ways, such as using panels that can rotate around an axis [7], using double-sided panels [8], or Energies 2022, 15, 7577 2 of 17 through cooling methods [9].Photovoltaic cells are the elements that convert solar radiation into electrical energy.The conversion efficiency is adversely affected by the temperature rise of the PV cells.This increase is associated with part of the absorbed solar radiation, which is converted into heat.The temperature reached by the cells is therefore an important factor to monitor.The need therefore arises to develop a mathematical model that is able to estimate the temperature field within a panel in a simple way.
Several authors proposed different methods for modelling PV modules as electrical circuits [10][11][12][13][14][15], many of them focusing on the PV cell model parameter estimation problem [11][12][13][14][15]. Saadaoui et al. [15] proposed a genetic algorithm to identify the parameters of photovoltaic panels.Abbassi et al. [16], used the Newton-Raphson method to determine the parameters of the PV panel.The model, then, was improved using a genetic algorithm.These studies aim to model the systems from an electrical point of view without assessing the thermal behavior of the panel.
The determination of the temperature of the photovoltaic cells is a different and very interesting problem, since the conversion efficiency of solar radiation into electrical energy depends on it.Sohani et al. [17] solved the thermal and electrical problem by focusing on the importance of the variation of thermophysical properties with temperature.They obtain interesting results, but the model is not quickly applicable as the temperatures are not known; therefore, a recursive method is required.Mavromatakis et al. [18] compared greatly simplified models that could be applied with a single equation.The results, however, have an excessive RMSE.They stated that a more complex method is needed to obtain more accurate results.Therefore, it can be useful to find the right compromise between implementation difficulty and accuracy.
Faiman D. [19] proposed a modified form of the Hottel-Whillier-Bliss equation, usually used for solar thermal collectors, in order to predict PV temperatures.Mattei M. et al. [20] proposed different simple models to estimate the PV temperature.The first uses the definition of NOCT and the others use the energy balance equation considering different correlations for the heat convective coefficient.They obtained a root-mean-square error (RMSE) of about 2.24 • C. Other simple methods to estimate the PV temperature were collected by Jakhrani et al. [21].They observed that most models show quite distinct values when compared among each other.The best models are those based on the energy balance equations.However, all these methods have turned out to be simple and with low accuracy.
Ceylan et al. [22] used an artificial neural network to predict module temperature.Their experimental study concerned the training of ANN according to solar radiation and air temperature detected in Turkey.For this investigation, ambient temperature was kept constant at values of 10, 20, 30 and 40 • C and the network cannot provide accurate results under conditions other than those used for training, including panel inclination and wind speed.Aly et al. [23] and Bevilacqua et al. [24] proposed a finite difference method to estimate the temperature profile.They managed to obtain very accurate results, but the method used is complex to implement because it divides the panel into numerous nodes.In the present study, a simpler finite difference approach will be used, using only the strictly necessary nodes, one for each layer.In this way, the model will be easily replicable, while maintaining good accuracy in the results.
Moreover, none of the previous studies considered the cloudiness of the sky to estimate the equivalent temperature of the celestial vault.This is an important parameter influencing radiative heat exchange that is often not measured and therefore cannot be used as input.The present work provides a method for estimating cloud cover using solar irradiance already considered as input.
The aim of this study is to develop a simple calculation method to determine the temperature profile inside the PV panel and to provide information on heat dissipation.The code will be experimentally validated with data acquired at the laboratory of the University of Calabria.The model implements the energy balance equations of the layers also considering the capacitive terms.Developing a non-stationary model allows good results to be obtained even when atmospheric conditions change suddenly, for example, when a cloud passes.Furthermore, it allows the energy balance equations to be solved in a non-recursive manner, evaluating temperatures at the next moment.The model will later be used to derive information about the role played by each layer in heat dispersion.This aspect is important because it helps to choose the best strategy in case a cooling system is to be implemented.
The rest of the paper proceeds as follows; in Section 2, the methodology to estimate the temperature profile is shown and the method to estimate the cloud cover is proposed.After providing information on the experimental setup, the model is validated in Section 3. Useful results for understanding the heat loss in the panel will then be shown.Finally, in Section 4, the results obtained are discussed.

Thermal Model
The typical structure of a photovoltaic panel consists of five layers: glass, EVA, PV cell, EVA and insulation.Each of them behaves differently with respect to heat transfer.To simulate their behaviour when exposed to solar radiation, each layer is thermally modelled using energy balance equations.Spatial discretization is performed by dividing the panel into five nodes which are representative of the single layers.
Figure 1 shows the stratigraphy of the panel and its division into nodes.The first point of the spatial discretization is fixed on the top surface of the glass layer.The temperature of this point is called T f g .The next point is fixed in the center of the EVA (Ethylene Vinyl Acetate) layer and its temperature is T ev,sup .The photovoltaic layer is modelled through a central node and the associated temperature is T pv .The next point is in the center of the EVA layer that divides it from the thermal insulation and its temperature is called T ev,in f .Finally, the last one representing the thermal insulating tedlar (polyvinyl fluoride) layer is placed on the lower surface and its temperature is called T ted .The top and bottom nodes are placed on the edges of the layers and not in the middle to facilitate data processing.These collocations, in fact, allow the temperatures at the extremities to be directly evaluated directly by solving the balance equations.These nodes are, however, still representative of the entire layer on which they are placed, both with reference to the calculation of thermal capacitances and conductive resistances.The nodes inside the panel exchange heat only by conduction with adjacent nodes.The nodes at the ends, on the other hand, also exchange heat by convection with the outside air and by radiation with the sky and ground.Incident global irradiance G is absorbed by the front glass and PV cells, in proportion to their absorption coefficient.Accumulated power also appears in the equations because the problem is time-variant.The problem is solved using the finite difference (FD) discretization technique.The derivative terms are therefore discretized in a finite time interval, ∆t.The equations, written for the defined nodes, are: T ev,sup − T pv T pv − T ev,in f where R i is the thermal resistance of the i-th layer, defined by the ratio of the thickness s i to the thermal conductivity k i ; C i is the heat capacity per unit area of the i-th layer, defined by the product of specific heat c i density ρ i and thickness s i .The optical properties appearing in the equations are α and τ, which represent the absorption and transmission coefficients, respectively.EVA layers are considered to be perfectly transparent to solar radiation.Moreover, solar radiation absorbed by the tedlar is considered negligible.The electrical efficiency η PV depends on the incident irradiance and on the cell temperature T pv .It is obtained with the Evans equation [25]: where β R is 0.006 • C −1 and γ is 0.085 for monocrystalline silicon panels.The temperature of standard reference conditions T R for performance classification is 25 • C. The subscript R refers to quantities acquired under standard test conditions.The unknowns are the temperatures at the current instant.Temperatures with apex p are already known as they are calculated at the previous time instant.The input data required to apply the model are the ambient temperature, wind speed and global irradiance incident on the tilted plane.Convective and radiative heat transfer coefficients are determined using the following procedure.

Heat Transfer Coefficients
Convective heat transfer coefficients are determined by considering both natural and forced contributions.The correlation used to combine the two effects is provided by Churchill [26]: With reference to natural convection, the convective coefficient is determined using the correspondent formulations of the Nusselt number [27].The equation for the upper surface is: Nu n = 0.13•Ra 1 3 (8) while for the lower surface, it is: Energies 2022, 15, 7577 5 of 17 The length used for the calculation of the Rayleigh number Ra is the average length of the two main dimensions of the panel: length L 1 and width L 2 .The convective heat transfer coefficient is therefore deduced from the definition of the Nusselt number.
With reference to forced convection, the equations for determining the Nusselt number depend on the type of flow and are considered equal for the upper and lower surface.In the case of turbulent flow, the Nusselt number is a function of Reynolds Re and Prandtl Pr numbers [28]: (11) where the characteristic length L c , used to calculate the Nu and Re number, is: In conditions of forced convection with laminar outflow, the classical expression is used [29]: Finally, the convective coefficient can be determined.
In order to determine the heat transfer coefficients for radiation, the non-linearized analytical formulas are used for both the lower and upper surfaces: where the subscript i can represent the upper surface or the lower surface and the subscript j can represent the sky or the ground.The view factors are: The ground temperature is set equal to the ambient temperature: Sky temperature is an important quantity to set as it influences heat exchange by radiation.Results obtained with different sky temperature models will be compared.The first method was proposed by Grag [30].The temperature of the sky is obtained by reducing the ambient air temperature by 20 degrees.
The second method considered is provided by Swinbank [31]: Energies 2022, 15, 7577 6 of 17 where the temperatures are expressed in Kelvin.The last method is an integration of the previous one [32].The sky temperature is a function of the air temperature T a and cloud cover N expressed in octaves of the sky covered by clouds:

Proposed Model for Estimating Sky Cloudiness
To apply Equation ( 21), it is necessary to estimate the cloudiness of the sky.Sky conditions are defined by how many eighths are covered by clouds.It is difficult, with the help of low-cost sensors, to be able to determine the cloudiness of the sky, yet it has a significant influence on the radiative heat exchange.Therefore, a method is proposed to estimate the cloud cover.It consists of making a comparison with the irradiance of clear sky conditions estimated by the ASHRAE model [33].The theoretical global irradiance G t for a clear day is: G t = I bn cos i + I do F f g→sky + (I bn sin θ s + I d0 )ρ gr F f g→gr (24) in which C 1 , C 2 and C 3 are tabulated values based on the period of the year [33], ρ gr is the albedo taken equal to 0.2, θ s is the solar altitude, I bn is the direct normal irradiance, I do is the diffuse irradiance on the horizontal plane, i is the incidence angle.The method for estimating cloud cover consists in comparing the real measured irradiation G with the theoretical one for clear day G t .The proposed formulation is the following, where the numerical values are obtained during the validation phase: To avoid slight transient phenomena, the cloud cover is then adjusted by making time averages over one-hour timesteps.The advantage of this model is that it allows the sky cloudiness to be estimated with reference only to the portion of the sky seen by the panel based on its inclination.The cloudiness of the portion of the sky not seen by the panel, whose cloudiness has no influence on radiative exchange, is therefore not considered.

Experimental Setup
The data used to validate the theoretical model were taken at the Solar Energy Laboratory of DIMEG at the University of Calabria (latitude 39.37 • N, longitude 16.23 • E).The photovoltaic panel is Schüco MPE 245 PG 60 FA.The size of the panel is 1663 × 998 mm and its inclination is 30 • towards the South (Figure 2).The nominal power of the panel is 245 W and the reference yield declared by the manufacturer η PV,R , evaluated under standard conditions, is 14.5%.The NOCT is 43 • C. The thermophysical characteristics of the panel layers are listed in Table 1 [30].
The temperature is measured by means of Pt 100 1/3 sensors with an error of 0.1 • C positioned on the back of the panels.The solar irradiance is measured directly on the inclined plane by means of a SMP3 Kipp&Zonen pyranometer with a non-linearity of 0.5% and uncertainty of 2%.The wind speed is measured with uncertainty of 2.5%.The experimental data used were acquired in January 2016 and July 2016 with a timestep of 1 min.The absorption coefficients of glass α f g and photovoltaic cells α PV are 0.05 and 0.93 respectively; the glass transmission coefficient τ f g is 0.9; the emissivity values of glass ε f g and tedlar ε ted are both 0.85.

Performance Metrics
The parameters that will be used to validate the model are: the Pearson correlation coefficient r xy , the mean bias error (MBE), the root-mean-square error (RMSE), the normalized root-mean-square error (NRMSE), and the Nush-Sutcliffe Efficiency (NSE).These indices will be calculated on minute-by-minute data and are defined by the following equations, where x and y are the calculated and experimental quantities respectively.
where x ave is the average of x, y ave is the average of y, y max and y min are the maximum and minimum values recorded.

Validation of Sky Temperature Model
The equations are solved by determining the temperatures progressively at the next instant of time.The time interval considered for time discretization is 1 min in analogy with the acquisition of experimental data.The equations are structured in such a way that the system is solved in the "implicit mode" and the time interval does not generate problems of instability of the solution.The model used to determine the sky temperature has a considerable influence on the results.In order to determine the best model for estimating the sky temperature, the formulations usually used in the literature and the one proposed by this study were considered.Figure 3 shows, as an example for 26 July, the comparison between the temperature values obtained for tedlar and the experimental values measured.The first model is based only on a reduction of 20 • C in the ambient temperature and approximates the results worse than the others.The day in question was cloudy at certain times in the afternoon.The Swinbank model and the proposed model show identical solutions in the morning hours, i.e., with clear skies.In the afternoon, due to the cloudiness of the sky, the panel temperature decreases.The proposed model, under these circumstances, provides a more accurate solution than the simplified model.It is able to recognize the cloudiness and a higher temperature is assigned to the sky, which reduces the radiative heat exchange with the PV panel.Figure 4 shows the root-mean-square error on the temperature obtained with reference to the measurements for the entire months of January and July.The proposed model allows the error to be reduced compared to the other methods usually used in the literature, obtaining an RMSE of 1.84 • C.  The model that will now be used to estimate the sky temperature will be the one proposed.Validation of the model will be conducted via winter and summer operation, with reference to both clear and cloudy days, with weak and strong winds.The model that will now be used to estimate the sky temperature will be the one proposed.Validation of the model will be conducted via winter and summer operation, with reference to both clear and cloudy days, with weak and strong winds.The model that will now be used to estimate the sky temperature will be the one proposed.Validation of the model will be conducted via winter and summer operation, with reference to both clear and cloudy days, with weak and strong winds.

Validation in Winter Conditions
The results obtained for a clear winter day (31 January) and an overcast winter day (2 January) are analyzed to validate the model.The climatic conditions are shown in Figure 5.The figure shows, on the left, the data for 31 January, which was a clear day, as can be seen from the solar radiation curve.On the right, the figure shows the climate data for 2 January, which was cloudy in the middle of the day and at certain times in the afternoon.During 31 January, the thermal variation was greatest; on 2 January the air temperature fluctuated between 5 and 14 • C and there was a reduction between 11:00 and 12:00, when the sky was cloudy.At the bottom, the figure shows the trends of wind speed, which was low in the morning for both days and then rose in the afternoon.These two days have been chosen because the wind was very variable and, therefore, it is possible to monitor the results with reference to flow conditions that can be laminar or turbulent.Therefore, the sky was cloudy between 11:00 and 12:00 on 2nd January.We do not know if a small cloud covered only the solar disk or if the whole sky was cloudy.To avoid minor transient phenomena, the cloud cover is estimated by averaging it over one hour.Figure 6 shows the percentage of cloud cover estimated with the proposed model.It can be noted that high cloud cover is observed in the central hours of the day, while at 14:00, the percentage of cloud cover is shown to be very low.In fact, the sky was probably not half-covered as it would appear from the reduction in solar irradiance, but a passing cloud may have covered only the solar disk.Thanks to the temporal averages, the cloud cover is about 85% in this first case, and it is less than 5% at 14:00.Before sunrise and after sunset, when the sun is not in the sky, the values are not of interest to us.However, in these cases, cloud cover is less accurately estimated based on the relative humidity value.Therefore, the sky was cloudy between 11:00 and 12:00 on 2 January.We do not know if a small cloud covered only the solar disk or if the whole sky was cloudy.To avoid minor transient phenomena, the cloud cover is estimated by averaging it over one hour.Figure 6 shows the percentage of cloud cover estimated with the proposed model.It can be noted that high cloud cover is observed in the central hours of the day, while at 14:00, the percentage of cloud cover is shown to be very low.In fact, the sky was probably not half-covered as it would appear from the reduction in solar irradiance, but a passing cloud may have covered only the solar disk.Thanks to the temporal averages, the cloud cover is about 85% in this first case, and it is less than 5% at 14:00.Before sunrise and after sunset, when the sun is not in the sky, the values are not of interest to us.However, in these cases, cloud cover is less accurately estimated based on the relative humidity value.
be noted that high cloud cover is observed in the central hours of the day, while at 14:00, the percentage of cloud cover is shown to be very low.In fact, the sky was probably not half-covered as it would appear from the reduction in solar irradiance, but a passing cloud may have covered only the solar disk.Thanks to the temporal averages, the cloud cover is about 85% in this first case, and it is less than 5% at 14:00.Before sunrise and after sunset, when the sun is not in the sky, the values are not of interest to us.However, in these cases, cloud cover is less accurately estimated based on the relative humidity value.The values of tedlar temperature, electrical power, and electrical efficiency are shown in Figure 7.The latter is calculated using Equation ( 6).The electrical power is obtained with the following equation: The values of tedlar temperature, electrical power, and electrical efficiency are shown in Figure 7.The latter is calculated using Equation ( 6).The electrical power is obtained with the following equation: Energies 2022, 15, x FOR PEER REVIEW 11 of 18 The highest temperature value was reached on the clear day around noon and was 47.4 °C.From 12:30 onwards, there was a drastic drop in temperature, presumably due to the increase in the wind speed.In the early afternoon it was constant and around 30 °C.On the overcast day, on the other hand, the temperature started to decrease from 11:00 onwards as the radiation decreased.The values obtained from the mathematical model are very similar to those measured.The statistical parameters also confirm this qualitative consideration (Table 2).In particular, the correlation coefficient CC is very high in both cases: over 0.99.The Nash-Sutcliffe efficiency (NSE) for the clear day is 0.989, which means that the accuracy of the model is excellent.It is slightly worse on the overcast day (0.964).The Mean Bias Error (MBE) is slightly negative, and this denotes a slight underestimation of the simulated data compared to the real ones.The root-mean-square error (RSME) is approximately 1.7 °C for the clear day and 1.8 °C for the overcast day.
With reference to the electrical power, the experimental values are well matched.It is only in the early afternoon that the estimated electrical power was slightly higher than the actual power for both days considered.This is due to the slight underestimation of the temperature in the same time slot.The statistical parameters confirm the good approximation of the model.The correlation index and the NSE are very high.The MBE is about half a Watt, and it is negative for the clear day and positive for the overcast day.The root-mean-square error is less than 3 W.
The electrical efficiency varies between 12 and 15% and this is confirmed by the experimental data.The statistical indices show, again, that the model adequately approximates the actual data.At 16:00 on 2nd January, the model shows a fluctuation in electrical efficiency.The experimental data do not agree with the prediction, as the actual oscillation is greater.These errors are usually present at dawn and dusk, when, in any case, electricity production is low.The highest temperature value was reached on the clear day around noon and was 47.4 • C. From 12:30 onwards, there was a drastic drop in temperature, presumably due to the increase in the wind speed.In the early afternoon it was constant and around 30 • C. On the overcast day, on the other hand, the temperature started to decrease from 11:00 onwards as the radiation decreased.The values obtained from the mathematical model are very similar to those measured.The statistical parameters also confirm this qualitative consideration (Table 2).In particular, the correlation coefficient CC is very high in both cases: over 0.99.The Nash-Sutcliffe efficiency (NSE) for the clear day is 0.989, which means that the accuracy of the model is excellent.It is slightly worse on the overcast day (0.964).The Mean Bias Error (MBE) is slightly negative, and this denotes a slight underestimation of the simulated data compared to the real ones.The root-mean-square error (RSME) is approximately 1.7 • C for the clear day and 1.8 • C for the overcast day.With reference to the electrical power, the experimental values are well matched.It is only in the early afternoon that the estimated electrical power was slightly higher than the actual power for both days considered.This is due to the slight underestimation of the temperature in the same time slot.The statistical parameters confirm the good approximation of the model.The correlation index and the NSE are very high.The MBE is about half a Watt, and it is negative for the clear day and positive for the overcast day.The root-mean-square error is less than 3 W.
The electrical efficiency varies between 12 and 15% and this is confirmed by the experimental data.The statistical indices show, again, that the model adequately approximates the actual data.At 16:00 on 2 January, the model shows a fluctuation in electrical efficiency.The experimental data do not agree with the prediction, as the actual oscillation is greater.These errors are usually present at dawn and dusk, when, in any case, electricity production is low.

Validation in Summer Conditions
The summer validation is referred to 11th July and 23rd July.Again, these represent a clear day and an overcast day.The global irradiance, ambient temperature and wind speed are shown in Figure 8.During the clear day, the temperature varied between 20 and 35 °C, and decreased in the afternoon due to high wind speeds.On 23 July it reached a maximum of about 32.5 °C around 10:00 and decreased first due to high wind speeds and, in the afternoon, due to the absence of solar radiation.The days, again, were chosen with varying wind speeds throughout the day.The percentage cloud cover estimated for 23rd July is shown in Figure 9.The sky is estimated to have been slightly overcast in the early morning hours (07:00-09:00), and two peaks of cloudiness are identified in the afternoon equal to 49 and 77%.The percentage cloud cover estimated for 23 July is shown in Figure 9.The sky is estimated to have been slightly overcast in the early morning hours (07:00-09:00), and two peaks of cloudiness are identified in the afternoon equal to 49 and 77%.The tedlar temperatures are shown in Figure 10.On 11 July, the maximum value reached was about 65 °C, well identified by the model.After 13:00, the temperature suddenly dropped to around 50 °C due to the turbulent flow generated by the high-velocity air.On 23 July, the panel cooled down slightly in the afternoon due to the reduction in solar irradiation.The model behaves adequately but, in this case, is less accurate than on the clear day.The correlation indexes are 0.996 for clear day and 0.983 for overcast day (Table 3).The NSE also decreases in the case of overcast day, passing from 0.99 to 0.957.The MBE indicates that, in both cases, the data are slightly underestimated.The rootmean-square error is 1.4 °C for the clear day and 2.25 °C for the overcast day.
The model predicts electrical output quite well, especially when irradiance is highly variable.For the clear day, the MBE is +1.7 W and RMSE is 3.4 W. CC and NSE are, instead, excellent.The model slightly overestimates electricity production at noon.With reference to the electrical efficiency, the NSE is about 0.98.For the overcast day, on the other hand, the MBE is about the same as in the previous case, but with a higher RMSE.Again, there is an overestimation at midday, although the temperature is well identified at this time.The tedlar temperatures are shown in Figure 10.On 11 July, the maximum value reached was about 65 • C, well identified by the model.After 13:00, the temperature suddenly dropped to around 50 • C due to the turbulent flow generated by the high-velocity air.On 23 July, the panel cooled down slightly in the afternoon due to the reduction in solar irradiation.The model behaves adequately but, in this case, is less accurate than on the clear day.The correlation indexes are 0.996 for clear day and 0.983 for overcast day (Table 3).The NSE also decreases in the case of overcast day, passing from 0.99 to 0.957.The MBE indicates that, in both cases, the data are slightly underestimated.The root-mean-square error is 1.4 • C for the clear day and 2.25 • C for the overcast day.The tedlar temperatures are shown in Figure 10.On 11 July, the maximum value reached was about 65 °C, well identified by the model.After 13:00, the temperature suddenly dropped to around 50 °C due to the turbulent flow generated by the high-velocity air.On 23 July, the panel cooled down slightly in the afternoon due to the reduction in solar irradiation.The model behaves adequately but, in this case, is less accurate than on the clear day.The correlation indexes are 0.996 for clear day and 0.983 for overcast day (Table 3).The NSE also decreases in the case of overcast day, passing from 0.99 to 0.957.The MBE indicates that, in both cases, the data are slightly underestimated.The rootmean-square error is 1.4 °C for the clear day and 2.25 °C for the overcast day.
The model predicts electrical output quite well, especially when irradiance is highly variable.For the clear day, the MBE is +1.7 W and RMSE is 3.4 W. CC and NSE are, instead, excellent.The model slightly overestimates electricity production at noon.With reference to the electrical efficiency, the NSE is about 0.98.For the overcast day, on the other hand, the MBE is about the same as in the previous case, but with a higher RMSE.Again, there is an overestimation at midday, although the temperature is well identified at this time.The model predicts electrical output quite well, especially when irradiance is highly variable.For the clear day, the MBE is +1.7 W and RMSE is 3.4 W. CC and NSE are, instead, excellent.The model slightly overestimates electricity production at noon.With reference to the electrical efficiency, the NSE is about 0.98.For the overcast day, on the other hand, the MBE is about the same as in the previous case, but with a higher RMSE.Again, there is an overestimation at midday, although the temperature is well identified at this time.Thanks to the model, it is possible to calculate the temperatures of the nodes for which the balance equations were written.It is possible, moreover, to identify the temperatures at the layer interfaces using a simple electrical analogy.Figure 11 shows that temperature is almost uniform within the panel, oscillating between 61.3 and 62.1 °C.The hottest component is the PV layer, while the coldest one is the outer surface of the glass, despite being directly exposed to sunlight.The temperatures on the two sides of the panel are different.It may seem anomalous that the temperature in the tedlar surface is higher than that in the glass surface.Tedlar, in fact, is a thermal insulator and it was expected that the temperature of its outer face could be close to the temperature of the air (or at least lower than the temperature of the face of the glass, which is also directly exposed to the sun's rays).The graphs in Figure 12 correlate the temperature difference between the PV layer and both the surfaces of the front glass and tedlar with the irradiation values.The figure shows that the internal temperature gradients are linearly proportional to the irradiance  incident on the panel.The evidence found on 11 July at noon (tedlar temperature higher than glass temperature) is confirmed under all operating conditions.
In fact, blue dots in the figure are almost always placed higher than the orange dots.Only for a few points, in low irradiation conditions, is there an overlap between the two Thanks to the model, it is possible to calculate the temperatures of the nodes for which the balance equations were written.It is possible, moreover, to identify the temperatures at the layer interfaces using a simple electrical analogy.Figure 11 shows that temperature is almost uniform within the panel, oscillating between 61.3 and 62.1 • C. The hottest component is the PV layer, while the coldest one is the outer surface of the glass, despite being directly exposed to sunlight.The temperatures on the two sides of the panel are different.It may seem anomalous that the temperature in the tedlar surface is higher than that in the glass surface.Tedlar, in fact, is a thermal insulator and it was expected that the temperature of its outer face could be close to the temperature of the air (or at least lower than the temperature of the face of the glass, which is also directly exposed to the sun's rays).
The graphs in Figure 12 correlate the temperature difference between the PV layer and both the surfaces of the front glass and tedlar with the irradiation values.The figure shows that the internal temperature gradients are linearly proportional to the irradiance G incident on the panel.The evidence found on 11 July at noon (tedlar temperature higher than glass temperature) is confirmed under all operating conditions.fore, the temperature deviation with tedlar is four times lower than with the front glass surface.In summer operating conditions, the graph is similar to the previous one and the linear proportionality with the global irradiance is confirmed.The interpolating equations change slightly.With reference to the temperature difference between the PV layer and the front glass, the equation is  = 0.0009194 • , with  = 0.995; the temperature difference with the tedlar, on the other hand, has a correlation of  = 0.0003049 • , with  = 0.977.The equations are quite similar between the winter and summer case, but in the latter case it turns out that the temperature difference with the tedlar is about three times smaller than the temperature difference with the front glass.The glass surface is therefore cooler than the back of the panel.However, this does not mean that it dissipates less heat by convection and radiation.On the contrary, the graphs in Figure 13 show that it is the cooler of the two surfaces (the front glass) and dissipates the most heat, under all conditions.The figure shows the trends of  ̇ and  ̇ indicating the thermal power lost to the outside by convection and radiation from the front glass and tedlar, respectively.The negative values of  ̇ , recorded in the morning and evening, represent the moments when the air temperature was higher than the surface; here, the back of the panel is heated by the air, by radiation, and by the environment.Heat loss is greatly influenced by wind speed.In fact, the fluctuating trends are caused by the varying wind intensity. ̇ is always higher than  ̇ by about 50-100 W/m 2 .In the summer period, the behavior is similar; however, the differences between the two curves are slightly smaller.In fact, blue dots in the figure are almost always placed higher than the orange dots.Only for a few points, in low irradiation conditions, is there an overlap between the two areas.In these few cases it may result in the temperature of the glass being higher than the temperature of the tedlar.
The figure distinguishes between winter and summer conditions.With reference to winter conditions, the temperature difference T pv − T f g is linearly linked with global radiance by means of the correlation y = 0.0010004•x, with R 2 = 0.981.The temperature difference, instead, between the PV layer and the tedlar surface is linked with global irradiance by means of the equation y = 0.000261•x, with R 2 = 0.939.On average, therefore, the temperature deviation with tedlar is four times lower than with the front glass surface.In summer operating conditions, the graph is similar to the previous one and the linear proportionality with the global irradiance is confirmed.The interpolating equations change slightly.With reference to the temperature difference between the PV layer and the front glass, the equation is y = 0.0009194•x, with R 2 = 0.995; the temperature difference with the tedlar, on the other hand, has a correlation of y = 0.0003049•x, with R 2 = 0.977.The equations are quite similar between the winter and summer case, but in the latter case it turns out that the temperature difference with the tedlar is about three times smaller than the temperature difference with the front glass.
The glass surface is therefore cooler than the back of the panel.However, this does not mean that it dissipates less heat by convection and radiation.On the contrary, the graphs in Figure 13 show that it is the cooler of the two surfaces (the front glass) and dissipates the most heat, under all conditions.The figure shows the trends of .Q f g and .Q ted indicating the thermal power lost to the outside by convection and radiation from the front glass and tedlar, respectively.The negative values of .Q ted , recorded in the morning and evening, represent the moments when the air temperature was higher than the surface; here, the back of the panel is heated by the air, by radiation, and by the environment.Heat loss is greatly influenced by wind speed.In fact, the fluctuating trends are caused by the varying wind intensity.
. Q f g is always higher than .
Q ted by about 50-100 W/m 2 .In the summer period, the behavior is similar; however, the differences between the two curves are slightly smaller.

Discussion
The model allows temperature results to be estimated with simple methodology and good accuracy.The overall RMSE obtained for temperature is 1.84 °C.This result is lower than that of other works in the literature, such as Mavromatakis et al. [19] (∼2.1-2.2 °C) and Mattei et al. [21] (2.24 °C).The error is greater when compared to more complex methods.For example, Bevilacqua et al. [30] and Aly et al. [23] obtained an RMSE of 1.37 °C and between 0.94 and 1.35 °C respectively by using numerous nodes to discretize the PV panel.The methodology presented in this paper for estimating cloud cover has given very good results.It is based on the comparison between the measured irradiance and the theoretical irradiance under clear sky conditions.It has made it possible to reduce the RMSE temperature from 2.11 to 1.84 °C.
The study provided a simple method for estimating the performance of a PV panel.In addition, the proposed new sky temperature relationship can be used for research purposes in the photovoltaic and other fields.
The temperature profile shows that the highest temperature is reached by the PV cell.However, the temperature remains quite homogeneous within the entire panel: the maximum gradient is about 1 °C.This characteristic allows an easier thermal dispersion of the absorbed heat.The temperature of the glass is cooler than the temperature of the back of the panel.In addition, the upper surface dissipates more heat power than the lower surface.The study, therefore, suggests that any cooling system should be applied on the back side to increase heat exchange.

Conclusions
The paper presented a simple one-dimensional transient model for estimating the temperature profile in a PV panel and the electrical power generated.The results were experimentally validated using a PV panel and the appropriate sensors placed on the 45/C building of the DIMEG at University of Calabria.The validation was carried out under both clear and overcast conditions, for both winter and summer days.In order to simulate overcast days, a methodology is proposed and validated for recognizing cloud cover.This information is necessary as it has a significant influence on the radiative heat exchange.
The model is able to predict the temperature of the back of the panel with an RMSE of about 1.5-2.0°.The RMSE on the estimation of the delivered electrical power is 2.5-4.3W. The RMSE on the electrical conversion efficiency is 0.25-0.75%.In both summer and winter periods, the MBE on the temperature estimation is negative, with values ranging between -1.0 °C and −0.5 °C.This means that the model returns slightly lower values than the real ones.The limitation of the present study is that the problem is considered onedimensional.The temperature obtained is assumed homogeneous for the faces.This methodology does not allow for the evaluation of edge effects and does not consider the

Discussion
The model allows temperature results to be estimated with simple methodology and good accuracy.The overall RMSE obtained for temperature is 1.84 • C.This result is lower than that of other works in the literature, such as Mavromatakis et al. [19] (∼2.1-2.2 • C) and Mattei et al. [21] (2.24 • C).The error is greater when compared to more complex methods.For example, Bevilacqua et al. [30] and Aly et al. [23] obtained an RMSE of 1.37 • C and between 0.94 and 1.The study provided a simple method for estimating the performance of a PV panel.In addition, the proposed new sky temperature relationship can be used for research purposes in the photovoltaic and other fields.
The temperature profile shows that the highest temperature is reached by the PV cell.However, the temperature remains quite homogeneous within the entire panel: the maximum gradient is about 1 • C.This characteristic allows an easier thermal dispersion of the absorbed heat.The temperature of the glass is cooler than the temperature of the back of the panel.In addition, the upper surface dissipates more heat power than the lower surface.The study, therefore, suggests that any cooling system should be applied on the back side to increase heat exchange.

Conclusions
The paper presented a simple one-dimensional transient model for estimating the temperature profile in a PV panel and the electrical power generated.The results were experimentally validated using a PV panel and the appropriate sensors placed on the 45/C building of the DIMEG at University of Calabria.The validation was carried out under both clear and overcast conditions, for both winter and summer days.In order to simulate overcast days, a methodology is proposed and validated for recognizing cloud cover.This information is necessary as it has a significant influence on the radiative heat exchange.
The model is able to predict the temperature of the back of the panel with an RMSE of about 1.5-2.0• .The RMSE on the estimation of the delivered electrical power is 2.5-4.3W. The RMSE on the electrical conversion efficiency is 0.25-0.75%.In both summer and winter periods, the MBE on the temperature estimation is negative, with values ranging between −1.0 • C and −0.5 • C.This means that the model returns slightly lower values than the real ones.The limitation of the present study is that the problem is considered one-dimensional.The temperature obtained is assumed homogeneous for the faces.This methodology does not allow for the evaluation of edge effects and does not consider the nonuniformity of the photovoltaic layer.However, the errors are acceptable because they are often lower than errors obtained through other scientific studies using more complex models.
The study found that the back of the panel reaches high temperatures and dissipates less thermal power than the front face.Facilitating cooling of the back face could be a useful solution to keep cell temperatures lower and improve panel efficiency.For this role, an effective solution is spray cooling, which involves spraying water mist droplets on the surface and waiting for them to evaporate by absorbing the latent heat of evaporation.The study shows that this solution should be adopted on the back face instead of the front face, and a study of its effectiveness represents a possible future development for this work.

18 Figure 1 .
Figure 1.Stratigraphy of the PV panel and nodes.The nodes inside the panel exchange heat only by conduction with adjacent nodes.The nodes at the ends, on the other hand, also exchange heat by convection with the outside air and by radiation with the sky and ground.Incident global irradiance  is absorbed by the front glass and PV cells, in proportion to their absorption coefficient.Accumulated power also appears in the equations because the problem is time-variant.The problem is solved using the finite difference (FD) discretization technique.The derivative terms are therefore discretized in a finite time interval, ∆.The equations, written for the defined nodes, are: ℎ ( −  ) + ℎ , → ( −  ) +  , −   + 1 2  + ℎ , → ( −  ) +   =   −  ∆ (1)

Figure 1 .
Figure 1.Stratigraphy of the PV panel and nodes.
The temperature is measured by means of Pt 100 1/3 sensors with an error of 0.1 °C positioned on the back of the panels.The solar irradiance is measured directly on the inclined plane by means of a SMP3 Kipp&Zonen pyranometer with a non-linearity of 0.5% and uncertainty of 2%.The wind speed is measured with uncertainty of 2.5%.The experimental data used were acquired in January 2016 and July 2016 with a timestep of 1 min.The absorption coefficients of glass  and photovoltaic cells  are 0.05 and 0.93 respectively; the glass transmission coefficient  is 0.9; the emissivity values of glass  and tedlar  are both 0.85.

18 Figure 3 .
Figure 3. Temperature of tedlar for different sky temperature models.

Figure 4 .
Figure 4. Temperature RMSE for different sky temperature models.

Figure 4 .
Figure 4. Temperature RMSE for different sky temperature models.

Figure 4 .
Figure 4. Temperature RMSE for different sky temperature models.

18 Figure 5 .
Figure 5. Ambient temperature, solar irradiance, and wind speed on 31st January and 2nd January.

Figure 5 .
Figure 5. Ambient temperature, solar irradiance, and wind speed on 31 January and 2 January.

Figure 7 .
Figure 7. Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.

Figure 7 .
Figure 7. Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.
The summer validation is referred to 11 July and 23 July.Again, these represent a clear day and an overcast day.The global irradiance, ambient temperature and wind speed are shown in Figure 8.During the clear day, the temperature varied between 20 and 35 • C, and decreased in the afternoon due to high wind speeds.On 23 July it reached a maximum of about 32.5 • C around 10:00 and decreased first due to high wind speeds and, in the afternoon, due to the absence of solar radiation.The days, again, were chosen with varying wind speeds throughout the day.Energies 2022, 15, x FOR PEER REVIEW 12 of 18

Figure 8 .
Figure 8. Ambient temperature, solar irradiance, and wind speed on 11th July and 23rd July.

Figure 8 .
Figure 8. Ambient temperature, solar irradiance, and wind speed on 11 July and 23 July.

Figure 10 .
Figure 10.Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.

Figure 10 .
Figure 10.Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.Figure 10.Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.

Figure 10 .
Figure 10.Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.Figure 10.Tedlar Temperature, Electrical Power and Efficiency for a clear day and overcast day.

3. 4 .
Figure 11 shows the temperature profile within the entire panel assessed by the model at 12:00 on 11 July.The outside air temperature is 34.6 • C, the global irradiance measured on the panel plane is 978 W/m 2 and the wind speed is 7.1 m/s.

3. 4 .
Figure 11 shows the temperature profile within the entire panel assessed by the model at 12:00 on 11 July.The outside air temperature is 34.6 °C, the global irradiance measured on the panel plane is 978 W/m 2 and the wind speed is 7.1 m/s.Thanks to the model, it is possible to calculate the temperatures of the nodes for which the balance equations were written.It is possible, moreover, to identify the temperatures at the layer interfaces using a simple electrical analogy.Figure11shows that temperature is almost uniform within the panel, oscillating between 61.3 and 62.1 °C.The hottest component is the PV layer, while the coldest one is the outer surface of the glass, despite being directly exposed to sunlight.The temperatures on the two sides of the panel are different.It may seem anomalous that the temperature in the tedlar surface is higher than that in the glass surface.Tedlar, in fact, is a thermal insulator and it was expected that the temperature of its outer face could be close to the temperature of the air (or at least lower than the temperature of the face of the glass, which is also directly exposed to the sun's rays).

Figure 11 .
Figure 11.Temperature of the panel at noon of July 11th.

Figure 11 .
Figure 11.Temperature of the panel at noon of 11 July.

Figure 12 .
Figure 12.Temperature difference between PV layer and front glass/tedlar surfaces.Winter and summer conditions.

Figure 12 .
Figure 12.Temperature difference between PV layer and front glass/tedlar surfaces.Winter and summer conditions.

Figure 13 .
Figure 13.Temperature difference between PV cells and front glass/tedlar surfaces.Winter and summer conditions.

Figure 13 .
Figure 13.Temperature difference between PV cells and front glass/tedlar surfaces.Winter and summer conditions.
35 • C respectively by using numerous nodes to discretize the PV panel.The methodology presented in this paper for estimating cloud cover has given very good results.It is based on the comparison between the measured irradiance and the theoretical irradiance under clear sky conditions.It has made it possible to reduce the RMSE temperature from 2.11 to 1.84 • C.

Table 1 .
Thermophysical characteristics of the layers.

Table 2 .
Experimental validation indices for winter conditions.

Table 2 .
Experimental validation indices for winter conditions.

Table 3 .
Experimental validation indices for summer conditions.

Table 3 .
Experimental validation indices for summer conditions.