Next Article in Journal
Manufacturing Challenges of a Modular Transverse Flux Alternator for Aerospace
Next Article in Special Issue
The Estimation of the Influence of Household Appliances on the Power Quality in a Microgrid System
Previous Article in Journal
Kinematics Measurement and Power Requirements of Fruitflies at Various Flight Speeds
Previous Article in Special Issue
Analysis of Primary Energy Factors from Photovoltaic Systems for a Nearly Zero Energy Building (NZEB): A Case Study in Lithuania
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature Analysis of the Stand-Alone and Building Integrated Photovoltaic Systems Based on Simulation and Measurement Data

by
Adam Idzkowski
*,
Karolina Karasowska
and
Wojciech Walendziuk
Faculty of Electrical Engineering, Bialystok University of Technology, Wiejska St. 45D, PL-15351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Energies 2020, 13(16), 4274; https://doi.org/10.3390/en13164274
Submission received: 12 July 2020 / Revised: 11 August 2020 / Accepted: 17 August 2020 / Published: 18 August 2020
(This article belongs to the Special Issue Innovations-Sustainability-Modernity-Openness in Energy Research 2020)

Abstract

:
Sunlight is converted into electrical energy due to the photovoltaic effect in photovoltaic cells. Energy yield of photovoltaic systems depends on the solar array location, orientation, tilt, tracking and local weather conditions. In order to determine the amount of energy produced in a photovoltaic system, it is important to analyze the operation of the photovoltaic (PV) arrays in real operating conditions and take into account the impact of external factors such as irradiance, ambient temperature or the speed of blowing wind, which is the natural coolant of PV panels. The analysis was carried out based on mathematical models and actual measurement data, regarding the dependence of the average temperature of PV arrays on variable and difficult to predict ambient conditions. The analysis used standard (nominal operating cell temperature (NOCT)), King, Skoplaki, Faiman and Mattei thermal models and the standard model for flat-plate photovoltaic arrays. Photovoltaic installations PV1, PV2a and PV2b, being part of the hybrid power plant of the Bialystok University of Technology, Poland, were the objects of the research. In the case of a free-standing solar system, the Skoplaki model proved to be the best method for determining the average temperatures of the PV arrays. For building-integrated PV systems, a corrected value of the mounting coefficient in the Skoplaki model was proposed, and the original results were compared. The comparison of the accuracy measures of the average operating temperatures for three micro-power plants, differently mounted and located, is presented.

Graphical Abstract

1. Introduction

In the last few years, the development of smart grids, which contain distributed photovoltaic plants and energy storage, has been observed. Power dissipation is intelligently managed [1] in them.
The aim of the paper is to describe the thermal analysis of three small photovoltaic arrays being part of a hybrid power plant of the Bialystok University of Technology, Poland. The motivation for writing this paper is to present to solar panel systems users that important issues such as mounting and ventilation are essential to minimize panel temperature increases.
The research is based on the measurement data obtained from the power plant data backup system collected during 2016–2019, and the simulations results. Moreover, information from scientific articles and books describing the influences of the weather conditions changes (irradiance, temperature or wind speed) on the parameters of the photovoltaic panels performances was considered. Basic notions concerning the functioning of photovoltaic arrays are described in books [2,3,4]. Energy yields of photovoltaic installations dependent on external conditions were analyzed and presented in papers [5,6,7,8,9]. Mathematical models and the quantities necessary for conducting the analysis are presented in papers [10,11]. The statistical analysis was worked out according to the dependencies described in article [12].
Irradiance is the most important and, usually, the only external condition influencing the work of modules that is taken into consideration. This parameter is crucial, because it is related to the number of photons that can be absorbed by the cells and participate in the photovoltaic energy conversion, while determining the potential of the solar radiation energy is essential for predicting the energy gain from the installation. The solar energy resources characteristic for a certain area are connected with the area location (latitude) and the local weather conditions, which are dependent on the landform. Individual characteristics of the array, such as orientation towards the azimuth or the module slope angles, are also of great significance. They determine the energy gain and the system efficiency [13].
The energetic efficiency of the installation is conditioned by a number of factors, such as the installation location, the way it was mounted and the local weather conditions; the amount of energy produced by a photovoltaic array depends on insolation, ambient temperature and wind speed [14]. The greater the irradiance, the higher the ambient temperature value, which causes the temperature increase of the photovoltaic cells and affects their performance. Additionally, the temperature increase changes the physical properties of the semi-conductive material building the cells, which results in a change of their current-voltage characteristics curve and the values of the obtained electric parameters [15,16].
The changing external conditions have a significant influence on the efficiency and durability of the photovoltaic cells; however, because of their stochastic character, the process is very complex. In the discussion over the quality of photovoltaic arrays, the functioning, lighting and temperature conditions of their work are usually considered. However, the influence of the speed of the wind, which is their natural cooler, plays a significant role, too [10,11].
Precisely estimating the average temperature of the modules (Tmod) is a key factor of a better evaluation of the photovoltaic (PV) array efficiency. The most frequently used mathematical models that estimate the temperature of modules include only the irradiance and ambient temperature, but the latest research proves that, also, the parameters of the blowing wind contribute to the increase of the PV panels’ efficiency. Increasing the heat dissipation intensity of the PV module influences its performance parameters [17]. Numerous authors, e.g., Skoplaki [18], Mattei [19], Faiman [20] and King [21], have undertaken intense work on new mathematical models, determining the average temperature of photovoltaic modules of micro-power plants, which include the effect of the blowing wind. The models were tested and validated by other authors, e.g., Koehl [22] and Kurtz [23].
The aim of the research presented in this article was to determine the most suitable model based on the lowest differences between the estimated and observed values of the average temperature of the photovoltaic modules. Moreover, the analysis of the accuracy of this temperature estimation for the whole of 2019, as well as for selected months (April, June, August and September, when the average temperature differences between modules and ambient were highest) of 2019, by introducing such measures as NRMSE (the normalized mean square error), NMBE (the normalized mean bias error) and the Pearson’s correlation coefficient, was done.
The contributions of the present work are an original comparison of the above-mentioned measures for three micro-power plants that are differently mounted and located. As a result of this analysis, the values of the parameters appearing in mathematical models from the literature [10,18,19] could be verified or corrected. These values are different for free-standing and building-integrated photovoltaic installations. The authors used a different value of the mounting coefficient in the Skoplaki model for building-integrated PV systems than the one appearing in the literature [18]. As it turned out, this had a great impact on the accuracy measures of the average estimated daily and monthly temperatures of the photovoltaic modules.

2. Mathematical Models Determining the Temperature of a Photovoltaic Module

The real temperature of a module depends on the material properties of the used semiconductor, the quality and the hermetic configuration of the module and on the surrounding environment and the weather conditions. This is a result of physical phenomena occurring in an illuminated PV cell, associated with the process of the light absorption, generation and transport of charge carriers; their recombination and many processes occurring simultaneously in the environment of the module [17].

2.1. Standard Model (Nominal Operating Cell Temperature (NOCT))

The nominal operating cell temperature (NOCT) describes the influence of all factors on the temperatures of cells in the modular constructions and provides necessary information concerning their thermal properties. The values of particular parameters in NOCT conditions are frequently used while determining the operating temperatures of PV modules. In order to estimate the average temperature value in real operational conditions, the following equation is applied [19]:
T m o d = T a + G 0 G 0 , N O C T ( T m o d , N O C T T a , N O C T ) ,
where:
  • Ta—ambient temperature (°C),
  • G0—the in-plane irradiance (W/m2),
  • G0,NOCT = 800 W/m2—the irradiance at which the NOCT is defined,
  • Ta,NOCT = 20 °C—the ambient temperature at which the NOCT is defined and
  • Tmod,NOCT = 44 °C–46 °C—the technology-dependent nominal operating cell temperature [10].

2.2. King Model

Another equation, called a two-component model, worked out in Sandia National Laboratory by King [21,24], is also used for this purpose:
T m o d = T a + G 0   e ( a + b v f ) ,
where:
  • a—the coefficient describing the effect of the radiation on the module temperature in the King model,
  • b—the coefficient describing the effect of cooling by the wind in the King model and
  • vf—the wind speed at a height of 10 meters (m/s).
The formula was applied by Kurtz to get the model coefficients a and b [10,23].

2.3. Skoplaki Models

Skoplaki is the author of two other formulas of estimating the temperature of modules that include the wind speed. One of the models (called further as Skoplaki), apart from the weather conditions, includes also various types of mounting the PV arrays and is expressed as follows [11]:
T m o d = T a + ω   0.32 G 0 8.91 + 2.0 v w   ,
where:
  • ω—mounting coefficient defined as the ratio of the Ross parameter for the mounting situation at-hand to the Ross parameter for a well-ventilated free-standing case. It takes on the values of 1, 1.2, 1.8 or 2.4, respectively, for free-standing installations, flat roofs, sloping roofs and facade-integrated photovoltaics, and
  • vw—the local wind speed around the module (m/s).
Whereas another model worked out by the same author includes certain properties of solar cells, such as the efficiency, the temperature power coefficient and the transmission and absorption coefficients [10,18]:
T m o d = T a + G 0 G 0 , N O C T ( T m o d , N O C T T a , N O C T ) h w , N O C T h w ( v ) [ 1 η S T C τ α ( 1 β S T C T a , S T C ) ] ,
where:
  • hw—wind convection coefficient (W·m−2 °C−1);
  • hw,NOCT—wind convection coefficient (W·m−2·°C−1) for wind speed at NOCT conditions, i.e., vw = 1 m/s;
  • v—wind speed (m/s);
  • τ∙α—the effective transmittance-absorbance product of the module;
  • ηSTC—efficiency coefficient of maximal power under standard test conditions (STC): irradiance G0,STC = 1000 W/m2, ambient temperature Ta,STC = 25 °C and air mass AM = 1.5 and
  • βSTC—temperature coefficient of maximal power (Pmax) under STC (%/°C); it is technology-dependent [25].
Additionally, in this model, there are three possible ways of parameterization of the wind convection heat transfer coefficient hw. Depending on the way of parameterization, the models are called: Skoplaki 1, Skoplaki 2 and Skoplaki 3.
For particular variants of the above model, the wind convection heat transfer coefficient hw is described with the following formulas [10,18]:
  • Skoplaki 1:
    h w = 8.91 + 2.00 v f   ,
  • Skoplaki 2:
    h w = 5.7 + 2.8 v w   .
For transformation between two different values of the wind speed, in Formulas (5) and (6), the following equation was applied [10,18]:
v w = 0.68 v f 0.5 ,
  • Skoplaki 3:
    h w = 8.3 + 2.2 v w 1 ,
    where:
  • vw1—wind directions perpendicular (±45°) to the module’s surface, and
    h w = 5.7 + 2.8 v w 2 ,
    where:
  • vw2—wind directions parallel (±45°) to the module’s surface.

2.4. Faiman Model

Another method of calculating the temperature of a PV module, including the wind conditions, was proposed by Faiman and is expressed as follows [10,22]:
T m o d = T a +   G 0 U 0 + U 1 v w ,
where:
  • U0—coefficient describing the effect of the radiation on the module temperature (W·m−2·°C−1), and
  • U1—coefficient describing the cooling by the wind (W·s·m−3·°C−1).
These coefficients were specified by Koehl in article [22] for selected technologies of PV modules.

2.5. Mattei Models

The mathematical model created by Mattei [19] is yet another way of determining the average temperature of photovoltaic modules:
T m o d =   U P V ( v ) T a + G 0 [ τ α η S T C ( 1 β S T C T c , S T C ) ] U P V ( v ) + β S T C η S T C G 0 ,
where:
  • UPV(v)—thermal losses coefficient from module to the surroundings (W·m−2·°C−1).
This model may take the name Mattei 1 or Mattei 2, depending on the method of determining the UPV(v) heat transfer coefficient, which is described by the equations presented below [10,19]:
  • Mattei 1:
    U P V ( v w ) = 26.6 + 2.3 v w ,   and
  • Mattei 2:
    U P V ( v w ) = 24.1 + 2.9 v w .

3. The Subject of the Research Description

Three photovoltaic installations being part of a hybrid power plant at the Bialystok University of Technology, Poland, initiated in 2015, are the subject of this research (Figure 1) [26]. The average annual irradiance value for the area of Bialystok in 2007–2016 was 1055 kWh/m2, and the average annual insolation was 1758.7 h [27].
The effective work of the discussed photovoltaic power plant was only possible during sunny periods, i.e., when the sky was cloudless and in the period between dawn and dusk. The shortest operating period in Bialystok lasts only 7 h and 37 min, which is only 30% of a day. Whereas the longest operating period in the summer is 16 h and 58 min, which stands for 70% of a day. Moreover, in winter, many days in Bialystok are cloudy, which is the reason why the value of insolation in this town within the period from October to March is very low. In the summer, the number of sunny days is about 50% to 70%, while, in winter, it is only 20% to 25%. As a result, in winter, the period of effective work of a PV power plant in Bialystok is very short, and the values of generated energy are rather low [6]. For the purpose of this work, the following installations were analyzed.
PV1—a photovoltaic installation of the rated peak power of 3 kWp, located on the roof of the building of Bialystok University of Technology Academic Business Incubator. Figure 2 presents a photograph of the PV generator composed of 12 polycrystalline ESP 250 6P modules set in two rows and mounted optimally, which means towards the azimuth facing the south and with the slope angle of 38°. The modules are connected with an inverter produced by SMA (Niestetal, Germany), type SUNNY BOY 3000 HF, of the maximum power of 3.15 kWp DC.
PV2a—a photovoltaic installation of the rated peak power of 1.5 kWp (Figure 3). It is located on the south-east wall of the building of Bialystok University of Technology Academic Business Incubator. Its generator includes 6 ESP 250 6P polycrystalline modules mounted at the angle of 90° towards the ground level and with the azimuth angle of 160°. The electric energy processing (DC/AC) was executed with the use of a GoodWe (Suzhou, China) NS-1500 inverter of the maximum power of 1.8 kWp DC.
PV2b—a photovoltaic installation of the rated peak power of 1.5 kWp, which, similarly to installation PV2a, is composed of 6 ESP 250 6P polycrystalline modules connected to a GoodWe NS-1500 inverter of the maximum power of 1.8 kWp DC. Its location, however, is different, as the PV generator is mounted on the south-west wall of the same building at the angle of 90° to the ground level and the azimuth angle of 250°. The view of this installation is presented in Figure 4.
All the analyzed photovoltaic installations exploit identical polycrystalline ESP 250 6P modules produced by Europe Solar Production (Munich, Germany). Table 1 and Table 2 contain the technical parameters of the modules and inverters installed in the described arrays.
The hybrid power plant of the Bialystok Technical University is also equipped with a specially created measurement system that records all important parameters of the photovoltaic panels and inverter performances, as well as the weather conditions.
Temperature and solar radiation measurements were conducted with the use of a distributed measurement system consisting of a central unit and measurement modules. The NI (National Instruments, Austin, TX, USA) PXIe-8108 measurement unit was equipped with an Intel® Core™ 2 Duo T9400 2.53 GHz dual core processor and an 80 GB HDD. Additionally, extension modules enabling the analysis of current parameters of the hybrid power station were placed in it. However, the main aim of the central unit was acquiring data from the measurement modules attached to the photovoltaic panels PV1, PV2a and PV2b. The NI-9148 external measurement modules (8-Slot, Spartan-3 FPGA and Ethernet CompactRIO Chassis), which were situated at the photovoltaic panels, were attached to the central unit via a LAN network. The modules conducted measurements of the temperature parameters of the photovoltaic modules (NI-9213+ K-type thermocouples, Table 3), as well as the global solar radiation (NI-9219 + LP PYRA 03 sensors, Table 4).
The temperature measurements of the photovoltaic modules were done by drilling a hole in the back surface of a module (Figure 5) and placing there a K-type thermocouple. This solution was considered to be sufficient, as the accuracy of the applied sensor was ±1.5 °C, and a small hole in the panel did not affect its mechanical structure.
Moreover, environmental parameters were also monitored with the use of a WS501-UMB (Lufft, Fellbach, Germany) meteorological station, which was connected to an NI PXIe-8108 computer through an RS485 interface. It should be mentioned that this compact weather station of a relatively good accuracy exploits the following technologies for the measurements of weather conditions: wind speed and direction—ultrasonic sensor, air temperature—negative temperature coefficient (NTC) thermistor, relative humidity (RH)—capacitive sensor, air pressure—microelectromechanical (MEMS) capacitive sensor and solar radiation—Kipp & Zonen (Delft, The Netherlands) pyranometer sensor [31].
The system of acquiring, storing, processing and visualizing the measurement data records and processes all information about the weather conditions, including the information concerning solar radiation, wind direction and speed, ambient temperature and the temperature of particular PV modules building the analyzed installation. The measured electrical and nonelectrical quantities are stored in a server with a database, and the current power plant performances are presented on a special website [26]. Selected environmental parameters and those related with the hybrid system performance are stored and shared for conducting calculations and scientific analysis. It is also possible to export the data as spreadsheet files, which enables precise data analysis in external software.

4. The Analysis of Average Monthly and Daily Estimated Temperatures of the PV Modules

The changing weather conditions have a significant impact on the performance of photovoltaic installations. The radiation intensity and the temperature of the PV modules are two main parameters that directly influence the amount of energy produced in a PV system. The biggest efficiency may be obtained through minimizing the temperature of a photovoltaic module and maximizing the insolation on its surface.

4.1. The Influence of Weather Conditions on the Operating Temperatures of PV Modules

The increase of irradiance is the main factor of the module power increase, but it also causes the increase of the operating temperature. This temperature value is also correlated with the ambient temperature, which can be seen in Figure 6. It presents the characteristics of the average observed monthly temperature of the PV1 modules, the ambient temperature and the wind speed during 2019. In order to determine the sheets presented below, the average monthly values of the parameters mentioned above, obtained from the data storage system of the Bialystok Technical University power plant, were used. The average measured monthly irradiance values in 2019 are presented in Table 5. They are the highest in April, June and August. The analysis takes into account the working time of the installation only from 8 a.m. to 4 p.m., which enables avoiding a calculation error resulting from significant changes of the irradiance and the temperature values before sunrise and after sunset.
It can be observed in Table 5 that, in April and in September, the monthly average operating temperatures of the PV modules had similar values (about 28 °C); however, the weather parameters were slightly different. Comparing April to June, the monthly average values of the wind speed and the monthly average values of irradiance were approximately the same, but there was a significant difference in the average ambient temperature (about 12 °C). In turn, comparing April to August, the monthly average values of irradiance were similar, but there were significant differences in the wind speed (about 1 m/s lower in August) and in the ambient temperature (about 10 °C higher in August).
The average temperature of the PV modules rises with the increase of the average ambient temperature, whereas its dependence from the average wind speed is not noticeable. Describing the correlation between these two quantities is quite complicated because of the fact that the wind conditions change constantly.

4.2. The Parameters of Mathematical Models Found in the Literature

In the literature, we can find a few ways of parametrizing the connection between the temperature of the modules, irradiance, ambient temperature and the wind speed, and the most frequently used for this purpose formulas are presented in Section 2 of this paper. Table 6 collects the values of the parameters occurring in the above-mentioned mathematical models.

4.3. The Graphical Comparison of the Selected Models for Estimating the Temperatures of the PV Modules

Charts representing the average observed and estimated monthly and daily temperature values of the PV modules for installations PV1, PV2a and PV2b were determined. Figure 7, Figure 8 and Figure 9 present the effects of the analysis for the whole of 2019, while Figure 10, Figure 11 and Figure 12 show the results obtained for April and September of 2019. Lines of different colors represent the average values of the module temperatures determined with the use of mathematical models, whereas the black dotted line stands for the average values of the measured temperatures of the photovoltaic modules.
Analyzing the characteristics presented in Figure 10, we can notice that the average daily values of the module temperatures for installation PV1 closest to the measurement values were obtained with the use of the Skoplaki model (the measured values are about 2–5% higher than the estimated values). When the Faiman and the standard models were applied, the estimated values were similar to the observed ones, but the values obtained on the basis of the Mattei 1, Skoplaki 1 and Skoplaki 2 models were significantly different from them. It is visible in Figure 11 that, for installation PV2a, determining the values that are similar to the measured ones is possible due to conducting calculations using the Skoplaki 1 and the standard methods. The effects of the calculations done with the use of the Faiman and Mattei 1 models have too-low values in relation to the temperature values determined on the basis of the measurement data. According to the characteristics obtained from installation PV2b (Figure 12), the results obtained with the use of the standard method best reflect the values determined in the data storage system (the differences between them are 0.1–3.5 °C). Good results are also noted for the Skoplaki 1 method. The results of calculations conducted with the use of other mathematical models differ significantly from the observed temperature values.

4.4. Methodology and Data Analysis

The methodology is based on monthly and daily data averaging and using the normalized error measures. The normalized mean square error (NRMSE), normalized mean bias error (NMBE) and Pearson’s correlation coefficient k values were carried out in accordance with the following formulas [12]:
R M S E = i = 1 n ( T m o d C i T m o d M i ) 2 n ,
N R M S E = R M S E T ¯ m o d M 100 % ,
M B E = i = 1 n ( T m o d C i T m o d M i ) n ,
N M B E = M B E T ¯ m o d M 100 % ,
k = i = 1 n ( T m o d C i T ¯ m o d C ) ( T m o d M i T ¯ m o d M ) [ i = 1 n ( T m o d C i T ¯ m o d C ) 2 ] [ i = 1 n ( T m o d M i T ¯ m o d M ) 2 ] ,
where:
  • T m o d C i —the estimated value of the temperature,
  • T ¯ m o d C —the average value of the estimated temperature,
  • T m o d M i —the observed (measured) value of the temperature,
  • T ¯ m o d M —the average value of the observed temperature,
  • i—the summing index and
  • n—the number of data used.
In the first part of the data analysis (Table 7), the average monthly values of the module temperatures in 2019, obtained with the use of the standard method, the Skoplaki, Skoplaki 1, Skoplaki 2, Faiman, Mattei 1, Mattei 2 and King models, were taken as estimated values. The measured average monthly values of the module temperatures of particular installations, determined on the basis on the data storage system, were taken as observed (measured) values. Due to the annual period of the analysis, the number of data used (estimated and corresponding observed) n was 12, whereas the values of the summing index i were consistent with the number of investigated months. The number of data points in a 30-day month was 122,000. The data was asynchronously acquired, the sampling rate was about eight per minute. The RAW data without filtration was used to create models. There were not unwanted spikes or missing data points in the datasets.
In the second part of the data analysis (Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13), the average daily temperature values in April, June, August and September of 2019 were taken as calculated values. The number of data used (estimated and corresponding observed) n was 30, and the values of the summing index i were consistent with the number of investigated days in a month. The number of data points in a day was 4066.
Summing up, the estimated average values of the module temperatures closest to the observed (measured) values in the case of installation PV1 were obtained with the use of the Skoplaki model (3); for installation PV2a, it was the Skoplaki (3) and Skoplaki 1 (5) models, whereas, for installation PV2b—the standard (1) and Skoplaki 1 (5) models. Applying other mathematical methods in the calculations leads to underestimating or overestimating the obtained values in relation to the values from the data storage system.
After many simulations, the authors observed that the value of the mounting coefficient ω representing the ventilation effect in the Skoplaki model, which in the literature is advised as 2.4 for facade-integrated installations, could be lower. The most optimal value for both the PV2a and PV2b installations is 2. This confirms the NMRSE and NMBE values presented in Table 12 and Table 13, which are lower for ω = 2 in the datasets from April, June, August and September.

4.5. The Impact of Wind Speed

An additional verification of the mathematical models was performed in Table 14. It was based on two datasets that were acquired on two sunny days of August 2019. Two days when the values of the average irradiance G were similar but the wind speeds vm were significantly different (1 m/s vs. 2.7 m/s) were chosen. The smallest differences between the average observed and estimated module temperatures ( T ¯ d i f f = T ¯ m o d M T ¯ m o d C ) were obtained for the Skoplaki* model with a modified mounting coefficient (PV1: ω = 1.2 and PV2a and PV2b: ω = 2). It proved to be a good model for all installations. The Skoplaki 1 model is equally well-fitting in the case of the PV2a and PV2b installations. Tolerable results are also reported for the Faiman model with the parameters proposed by Koehl. However, the Mattei 1, 2 and the King models are poorly fitting for the parameters from Table 6.

5. Conclusions

Photovoltaic installations PV1, PV2a and PV2b, belonging to a hybrid power plant at the Bialystok University of Technology, were subject to the research. A precise analysis of the real influence of weather conditions, obtained during long-term measurements, enabled to learn in detail the behavior of the arrays in changeable and unpredictable environmental conditions that have occurred lately in Poland. Due to that, for the purpose of the presented research, tests concerning the dependence of the PV system operating temperatures on the weather conditions, with a special concern for the wind speed, were conducted on the basis of the collected measurement data and mathematical models. For this purpose, different mathematical models: King, Skoplaki, Skoplaki 1, Skoplaki 2, Faiman, Mattei 1 and Mattei 2 and the standard model were analyzed. They enabled to determine the operating temperatures of photovoltaic modules with weather conditions, including wind, taken into consideration. Average values of the ambient temperature, wind speed and the irradiance for the tested installations in 2019 obtained from the data storage system of the university power plant were used to conduct the calculations. The graphs created as a result of the tests showed the characteristics of the temperature mean value changes for the PV1, PV2a and PV2b installations.
On the basis of the conducted analysis, it can be stated that it is hard to choose a general theoretical model for determining the operating temperature of the PV modules that would be relevant for each installation. In the case of a detached array (PV1), the Skoplaki model appeared to be the best one (NRMSE was 6% maximum and NMBE—5% maximum; see Table 7, Table 8, Table 9, Table 10 and Table 11). Whereas, for the installation PV2b mounted on the south-west wall, the values obtained from the data storage system were better expressed by the standard model (NOCT) and the Skoplaki 1 model. In the case of the NOCT model, the NRMSE was 10% maximum and the NMBE—8% maximum. In turn, the Skoplaki and Skoplaki 1 models were more suitable for installation PV2a (on the south-east wall). The error values were slightly higher—the maximum NRMSE was 11% (Skoplaki 1) and 16% (Skoplaki). This may result from variations in the ventilation of the modules for various types of mountings of the installations. In the detached systems, the blowing wind has a positive influence on the modules’ work, because it receives the heat excess from their back surfaces. However, in the installations mounted on a facade, it is not possible, because the air flow is limited by the wall of the building.
The results of the conducted analysis (Table 14) confirm a significant influence of the cooling effect, which naturally decreases the operating temperature of the PV modules and improves the negative influence of the temperature increase on their operating parameters. A good cohesion of the results of the calculations based on the Skoplaki model with the observed measurement data proves the need to consider also the influence of the wind speed while determining the operating temperature of photovoltaic detached installations. The mounting issue is essential as well. Based on earlier observations of the accuracy measures, the authors of this paper proposed a lower value of the mounting coefficient ω than in the literature that was more relevant for the PV2a and PV2b facade-based installations. After this modification, there was a significant decrease of the NRMSE and the NMBE values (Table 12 and Table 13).
The average daily values of the errors were larger than the values given by the authors of publication [10]. In this publication, 15 min, hourly and daily RMSE values were presented for different panel technologies. Moreover, the RMSE measure was used, and the presented results varied from 2 °C to 7 °C for the polycrystalline (p-Si) panels, i.e., from 5.7% to 20% at 35 °C. We decided to use the normalized NRMSE measure, which is related to the average operating temperature of PV panels. The achieved results were up to 62% (Table 10). For the Skoplaki model, it was 6% (PV1), 10% (PV2a) and 11% (PV2b) maximums. This model was the most suitable for the studied PV system configurations.

Author Contributions

Conceptualization, A.I. and K.K.; methodology, A.I. and K.K.; software, K.K.; validation, A.I., K.K. and W.W.; formal analysis, W.W.; investigation, K.K.; resources, K.K.; data curation, K.K.; writing—original draft preparation, K.K., A.I. and W.W.; writing—review and editing, W.W. and A.I.; visualization, W.W. and K.K.; supervision, A.I.; project administration, A.I. and funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by project WZ/WE-IA/2/2020 of the Bialystok University of Technology and financed from a subsidy provided by the Ministry of Science and Higher Education.

Acknowledgments

The research could take place thanks to the earlier realization of two projects: a hybrid system of small wind and photovoltaic power designed to supply electricity to the Research and Teaching Center of the Electrical Faculty at Zwierzyniecka Street 10, the project No. WND-RPPD. 05. 02. 00-20-034/12, entitled “Improving the energy efficiency of the infrastructure of the Bialystok University of Technology using renewable energy sources”, Priority Axis V. Development of environmental infrastructure, Action 5.2 Development of local environmental infrastructure, co-financed by the European Regional Development Fund under the Regional Operational Programme of the Podlasie Region. Additionally, the operation of the production system in urban (urbanized) conditions is examined in the framework of project No. WND-RPPD. 01. 01. 00-20-015/ 12, entitled “Study of the effectiveness of active and passive methods of improving the energy efficiency of infrastructure using renewable energy sources” priority axis I. Increase innovation and support entrepreneurship in the region, Actions 1.1. Creating conditions for the development of innovation, co-financed by the European Regional Development Fund and the state budget under the Regional Operational Programme of the Podlasie Region.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alazab, M.; Khan, S.; Krishnan, S.S.R.; Pham, Q.V.; Reddy, M.P.K.; Gadekallu, T.R. A multidirectional LSTM model for predicting the stability of a smart grid. IEEE Access 2020, 8, 85454–85463. [Google Scholar] [CrossRef]
  2. Markvart, T. (Ed.) Solar Electricity, 2nd ed.; Wiley: Hoboken, NJ, USA, 2000; ISBN 978-0-471-98852-6. [Google Scholar]
  3. Petrone, G.; Ramos-Paja, C.A.; Spagnuolo, G.; Xiao, W. Photovoltaic Sources Modeling; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015; ISBN 9781118679036. [Google Scholar]
  4. Xiao, W. Photovoltaic Power System: Modeling, Design, and Control; Wiley: Hoboken, NJ, USA, 2017; ISBN 9781119280347. [Google Scholar]
  5. Gökmen, N.; Hu, W.; Hou, P.; Chen, Z.; Sera, D.; Spataru, S. Investigation of wind speed cooling effect on PV panels in windy locations. Renew. Energy 2016, 90, 283–290. [Google Scholar] [CrossRef]
  6. Kusznier, J.; Wojtkowski, W. Impact of climatic conditions on PV panels operation in a photovoltaic power plant. In Proceedings of the 15th Selected Issues of Electrical Engineering and Electronics, WZEE 2019, Zakopane, Poland, 8–10 December 2019. [Google Scholar]
  7. Nkurikiyimfura, I.; Safari, B.; Nshingabigwi, E. A simulink model of photovoltaic modules under varying environmental conditions. IOP Conf. Ser. Earth Environ. Sci. 2018, 159, 012024. [Google Scholar] [CrossRef] [Green Version]
  8. Perovic, B.; Klimenta, D.; Jevtic, M.; Milovanovic, M. A transient thermal model for flat-plate photovoltaic systems and its experimental validation. Elektronika Ir Elektrotechnika 2019, 25, 40–46. [Google Scholar] [CrossRef] [Green Version]
  9. Mora Segado, P.; Carretero, J.; Sidrach-de-Cardona, M. Models to predict the operating temperature of different photovoltaic modules in outdoor conditions. Prog. Photovolt. Res. Appl. 2015, 23, 1267–1282. [Google Scholar] [CrossRef]
  10. Schwingshackl, C.; Petitta, M.; Wagner, J.E.; Belluardo, G.; Moser, D.; Castelli, M.; Zebisch, M.; Tetzlaff, A. Wind effect on PV module temperature: Analysis of different techniques for an accurate estimation. Energy Procedia 2013, 40, 77–86. [Google Scholar] [CrossRef] [Green Version]
  11. Veldhuis, A.J.; Nobre, A.; Reindl, T.; Ruther, R.; Reinders, A.H.M.E. The influence of wind on the temperature of PV modules in tropical environments, evaluated on an hourly basis. In Proceedings of the IEEE Photovoltaic Specialists Conference, Tampa, FL, USA, 16–21 June 2013. [Google Scholar]
  12. Frydrychowicz-Jastrzebska, G.; Bugała, A. Modeling the distribution of solar radiation on a two-axis tracking plane for photovoltaic conversion. Energies 2015, 8, 1025–1041. [Google Scholar] [CrossRef]
  13. Armstrong, S.; Hurley, W.G. A thermal model for photovoltaic panels under varying atmospheric conditions. Appl. Eng. 2010, 30, 1488–1495. [Google Scholar] [CrossRef]
  14. Chen, J.L.; He, L.; Yang, H.; Ma, M.; Chen, Q.; Wu, S.J.; Xiao, Z. Empirical models for estimating monthly global solar radiation: A most comprehensive review and comparative case study in China. Renew. Sustain. Energy Rev. 2019, 108, 91–111. [Google Scholar] [CrossRef]
  15. Idzkowski, A.; Walendziuk, W.; Borawski, W. Analysis of the temperature impact on the performance of photovoltaic panel. In Proceedings of the SPIE—The International Society for Optical Engineering; Romaniuk, R.S., Ed.; International Society for Optics and Photonics: Bellingham, WA, USA, 2015; Volume 9662, p. 96620. [Google Scholar]
  16. Kandemir, E.; Borekci, S.; Cetin, N.S. Conventional and soft-computing based MPPT methods comparisons in direct and indirect modes for single stage PV systems. Elektronika Ir Elektrotechnika 2018, 24, 45–52. [Google Scholar] [CrossRef] [Green Version]
  17. Jaszczur, M.; Teneta, J.; Hassan, Q.; Majewska, E.; Hanus, R. An experimental and numerical investigation of photovoltaic module temperature under varying environmental conditions. Heat Transf. Eng. 2019. [Google Scholar] [CrossRef]
  18. Skoplaki, E.; Boudouvis, A.G.; Palyvos, J.A. A simple correlation for the operating temperature of photovoltaic modules of arbitrary mounting. Sol. Energy Mater. Sol. Cells 2008, 92, 1393–1402. [Google Scholar] [CrossRef]
  19. Mattei, M.; Notton, G.; Cristofari, C.; Muselli, M.; Poggi, P. Calculation of the polycrystalline PV module temperature using a simple method of energy balance. Renew. Energy 2006, 31, 553–567. [Google Scholar] [CrossRef]
  20. Faiman, D. Assessing the outdoor operating temperature of photovoltaic modules. Prog. Photovolt. Res. Appl. 2008, 16, 307–315. [Google Scholar] [CrossRef]
  21. King, D.L.; Boyson, W.E.; Kratochvil, J.A. Photovoltaic array performance model. Sandia Rep. 2004, 8, 1–41. [Google Scholar] [CrossRef] [Green Version]
  22. Koehl, M.; Heck, M.; Wiesmeier, S.; Wirth, J. Modeling of the nominal operating cell temperature based on outdoor weathering. Sol. Energy Mater. Sol. Cells 2011, 95, 1638–1646. [Google Scholar] [CrossRef]
  23. Kurtz, S.; Whitfield, K.; Miller, D.; Joyce, J.; Wohlgemuth, J.; Kempe, M.; Dhere, N.; Bosco, N.; Zgonena, T. Evaluation of high-temperature exposure of rack-mounted photovoltaic modules. In Proceedings of the 34th IEEE Photovoltaic Specialists Conference, Philadelphia, PA, USA, 7–12 June 2009. [Google Scholar]
  24. Mavromatakis, F.; Kavoussanaki, E.; Vignola, F.; Franghiadakis, Y. Measuring and estimating the temperature of photovoltaic modules. Sol. Energy 2014, 110, 656–666. [Google Scholar] [CrossRef]
  25. HOMER Energy. How HOMER Calculates the PV Array Power Output. Available online: https://www.homerenergy.com/products/pro/docs/latest/how_homer_calculates_the_pv_array_power_output.html (accessed on 13 June 2020).
  26. Bialystok University of Technology. Hybrid Power Plant of the Bialystok University of Technology. Available online: http://elektrownia.pb.edu.pl/ (accessed on 13 June 2020).
  27. European Commission. JRC Photovoltaic Geographical Information System (PVGIS)—European Commission. Available online: https://re.jrc.ec.europa.eu/pvg_tools/en/tools.html (accessed on 13 June 2020).
  28. Europe Solar Production. Polycrystalline Photovoltaic Module Europe Solar Production Premium Quality Solar Module Data Sheet ESP 6P 240-255 Wp. Available online: http://europe-solarproduction.com/media/690/ESP-Polycrystalline-Solar-Module-Datasheets-ESP-6P-series.pdf (accessed on 17 August 2020).
  29. SMA Solar Technology AG. SUNNY BOY 2000HF/2500HF/3000HF—Installation Guide. Available online: https://www.solartradesales.co.uk/Cache/Downloads/SunnyBoy-HF-Installation-guide-3.pdf (accessed on 17 August 2020).
  30. Solahart Goodwe Single Phase Small Domestic Inverter. GW1500-NS & GW3000-NS. Available online: https://www.solahart.com.au/media/5305/ih0113_gw1500-ns-gw3000-ns_single-phase-inverters_june-2019_web.pdf (accessed on 13 June 2020).
  31. Compact Weather Sensors—WS501-UMB Smart Weather Sensor. Available online: https://www.lufft.com/products/compact-weather-sensors-293/ws501-umb-smart-weather-sensor-1839/ (accessed on 24 July 2020).
  32. National Instruments. NI 9213 Datasheet. Available online: https://www.ni.com/pdf/manuals/378021a_02.pdf (accessed on 24 July 2020).
  33. National Instruments. NI 9219 Datasheet. Available online: https://www.ni.com/pdf/manuals/377223a_02.pdf (accessed on 24 July 2020).
  34. Delta OHM. LP Pyra 02—LP Pyra 03—LP Pyra 12. Available online: https://www.deltaohm.com/ver2012/download/LP_PYRA02_03_12_uk.pdf (accessed on 24 July 2020).
Figure 1. The general scheme of the hybrid power plant at the Bialystok University of Technology, where: WT1 and WT2—wind turbines, PV1—optimal located photovoltaic installation, PV2a and PV2b—wall-mounted photovoltaic installations, PV3—tracker photovoltaic installation (not analyzed in this paper), WS—weather station and EC&DC—electrical converters and the main data analysis center.
Figure 1. The general scheme of the hybrid power plant at the Bialystok University of Technology, where: WT1 and WT2—wind turbines, PV1—optimal located photovoltaic installation, PV2a and PV2b—wall-mounted photovoltaic installations, PV3—tracker photovoltaic installation (not analyzed in this paper), WS—weather station and EC&DC—electrical converters and the main data analysis center.
Energies 13 04274 g001
Figure 2. A view of a photovoltaic installation with panels set optimally.
Figure 2. A view of a photovoltaic installation with panels set optimally.
Energies 13 04274 g002
Figure 3. A view of the photovoltaic installation on the south-east wall.
Figure 3. A view of the photovoltaic installation on the south-east wall.
Energies 13 04274 g003
Figure 4. A view of the photovoltaic installation on the south-west wall.
Figure 4. A view of the photovoltaic installation on the south-west wall.
Energies 13 04274 g004
Figure 5. A view of the thermocouple sensor location in photovoltaic modules.
Figure 5. A view of the thermocouple sensor location in photovoltaic modules.
Energies 13 04274 g005
Figure 6. The characteristics of the average observed monthly temperature of the PV1 installation modules, the ambient temperature and the wind speed in 2019.
Figure 6. The characteristics of the average observed monthly temperature of the PV1 installation modules, the ambient temperature and the wind speed in 2019.
Energies 13 04274 g006
Figure 7. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV1 installation in 2019.
Figure 7. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV1 installation in 2019.
Energies 13 04274 g007
Figure 8. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV2a installation in 2019.
Figure 8. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV2a installation in 2019.
Energies 13 04274 g008
Figure 9. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV2b installation in 2019.
Figure 9. The characteristics of the average observed and estimated monthly values of the module temperatures of the PV2b installation in 2019.
Energies 13 04274 g009
Figure 10. The characteristics of the average observed and estimated daily values of the module temperatures of the PV1 installation in September of 2019.
Figure 10. The characteristics of the average observed and estimated daily values of the module temperatures of the PV1 installation in September of 2019.
Energies 13 04274 g010
Figure 11. The characteristics of the average observed and estimated daily values of the module temperatures of the PV2a installation in September of 2019.
Figure 11. The characteristics of the average observed and estimated daily values of the module temperatures of the PV2a installation in September of 2019.
Energies 13 04274 g011
Figure 12. The characteristics of the average observed and estimated daily values of the module temperatures of the PV2b installation in September of 2019.
Figure 12. The characteristics of the average observed and estimated daily values of the module temperatures of the PV2b installation in September of 2019.
Energies 13 04274 g012
Table 1. A sheet of technical parameters of the photovoltaic module used in the analyzed arrays [28]. STC: standard test conditions and NOCT: nominal operating cell temperature.
Table 1. A sheet of technical parameters of the photovoltaic module used in the analyzed arrays [28]. STC: standard test conditions and NOCT: nominal operating cell temperature.
CompanyEurope Solar Production
ModelESP 250 6P
Dimensions1640 × 990 × 40
Electrical Data at STC
Module Efficiency ηSTC15.3%
Peak Power Watts Pmax,STC250Wp
Maximum Power Voltage Vmpp,STC30.93V
Maximum Power Current Impp,STC8.08A
Electrical Data at NOCT
Peak Power Watts Pmax,NOCT182Wp
Maximum Power Voltage Vmpp,NOCT28.04V
Maximum Power Current Impp,NOCT6.49A
Temperatures Ratings
Temperature Coefficient of Pmax (βSTC)−0.46%/°C
Temperature Coefficient of Open-Circuit Voltage VOC−0.34%/°C
Temperature Coefficient of Short-Circuit Voltage Isc0.07%/°C
Table 2. A datasheet of the technical parameters of the inverters used in installations PV1 [29], PV2a and PV2b [30].
Table 2. A datasheet of the technical parameters of the inverters used in installations PV1 [29], PV2a and PV2b [30].
CompanySMA Solar Technology AGGoodWe
ModelSunny Boy 3000 HFGW 1500-NS
DC Input
Maximum Input Voltage700V500V
Minimum Input Voltage175V80V
Rated Input Voltage530V360V
MPP Voltage Range210–560V80–450V
Maximum Input Current15A10A
AC Output
Maximum Output Current15A7.5A
Rated Power at 230 V, 50 Hz3000W1500W
Maximum Apparent AC Power3000VA1500VA
Efficiency
Maximum efficiency ηmax96.3%97.0A
European weighted efficiency ηEU95.4%96.0W
Table 3. The selected parameters of the measurement units used for temperature measurements of photovoltaic modules and for measuring the output voltage of the pyranometer sensor [32,33]. PV: photovoltaic.
Table 3. The selected parameters of the measurement units used for temperature measurements of photovoltaic modules and for measuring the output voltage of the pyranometer sensor [32,33]. PV: photovoltaic.
ParameterMeasurement of PV Module Temperature
NI-9213—Thermocouple Module, Voltage Measurement with Conversion to Temperature
Pyranometer Data Acquisition
NI-9219—Universal Multipurpose Module, Voltage
Measurement
Number of channels16 thermocouple channels4 analog input channels
Analog-to-digital converter (ADC) resolution24 bits24 bits
Type of ADCDelta-sigma
(with analog prefiltering)
Delta-sigma
(with analog prefiltering)
Sampling rate1 S/s2 S/s
Gain error in high-resolution mode at −40 °C to 70 °C0.07% typical,
0.15% maximum
0.3% typical,
0.4% maximum
Offset error in high-resolution mode at −40 °C to 70 °C4 μV typical,
6 μV maximum
6 μV typical,
180 μV maximum
OtherCold-junction compensation
−40 °C to 70 °C
1.1 °C typical, 2.1 °C maximum
Stability of the gain drift ± 20 (ppm of reading/°C)
Table 4. The selected parameters of the LP PYRA 03 pyranometer (Delta OHM, Caselle di Selvazzano, Italy) used in the measurement system [34].
Table 4. The selected parameters of the LP PYRA 03 pyranometer (Delta OHM, Caselle di Selvazzano, Italy) used in the measurement system [34].
ParameterValue
Classsecond (ISO 9060)
Typical sensitivity 10 μV/(W/m2)
Impedance33–45 Ω
Measuring range0 ÷ 2000 W/m2
Viewing Field2π sr
Spectral Field305–2800 nm
Operating temperature−40 °C–80 °C
Table 5. Average measured monthly values of the weather parameters and operating temperatures of the PV1 modules.
Table 5. Average measured monthly values of the weather parameters and operating temperatures of the PV1 modules.
MonthAverage Wind Speed vw (m/s)Average Ambient Temperature Ta (°C)Average Irradiance G0 (W/m2)Average Operating Temperature Tmod (°C)
I2.84−3.27120.201.24
II3.202.71208.068.88
III3.905.74339.0515.08
IV2.9511.95556.2027.81
V2.1610.68310.3420.10
VI2.8423.48592.9439.75
VII3.2018.95445.2831.28
VIII2.2220.60533.6036.37
IX3.2015.73458.1428.17
X2.4010.76280.3918.98
XI2.845.2966.388.54
XII2.442.6781.775.99
Table 6. The juxtaposition of the parameters used in the calculations of the temperatures of the PV modules [10,11,18,19,21,22].
Table 6. The juxtaposition of the parameters used in the calculations of the temperatures of the PV modules [10,11,18,19,21,22].
ModelParameterValueQuantities
FaimanU030.02W·m−2·°C−1
U16.28W·s·m−3·°C−1
KingaStand-alone power system PV1−3.56-
Building-integrated power systems PV2a and PV2b−2.81
KingbStand-alone power system PV1−0.0750s·m−1
Building-integrated power systems PV2a and PV2b−0.0455
Skoplaki ω Stand-alone power system PV1 (flat-roof)1.2-
Building-integrated power systems PV2a and PV2b2.4 (2.2-2.6)
Skoplaki 1 and 2 τ·α0.90-
Mattei 1 and 20.81
Table 7. The dataset of the normalized root mean square error (NRMSE), normalized mean bias error (NMBE) values and the Pearson’s correlation coefficient for the average monthly values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models during 2019. K: Pearson’s correlation coefficient.
Table 7. The dataset of the normalized root mean square error (NRMSE), normalized mean bias error (NMBE) values and the Pearson’s correlation coefficient for the average monthly values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models during 2019. K: Pearson’s correlation coefficient.
PV1PV2aPV2b
NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)
Standard8.275.311.0020.16−19.531.0010.24−8.051.00
Skoplaki5.29−4.871.007.951.921.0016.9310.490.99
Skoplaki 124.2219.500.8310.98−10.390.869.00−0.180.90
Skoplaki 253.3944.960.9912.076.060.9920.6014.030.99
Faiman/Koehl14.53−13.821.0033.65−31.931.0019.35−18.831.00
Mattei 122.72−7.920.9842.98−33.490.9840.48−30.070.95
Mattei 222.60−7.350.9842.83−33.190.9840.56−29.920.95
King52.99−47.890.9858.09−53.960.9739.22−37.920.99
Table 8. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in April, 2019.
Table 8. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in April, 2019.
PV1PV2aPV2b
NRMS E (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)
Standard12.757.990.9620.44−16.630.909.43−6.810.97
Skoplaki5.93−4.920.9914.4910.300.9720.3417.880.96
Skoplaki 126.4924.430.488.63−5.440.528.733.760.56
Skoplaki 259.6955.580.9420.3416.000.9727.1223.320.91
Faiman/Koehl16.20−15.670.9934.29−33.090.9522.08−21.250.97
Mattei 122.21−15.480.9841.39−37.800.8939.05−34.850.88
Table 9. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in June, 2019.
Table 9. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in June, 2019.
PV1PV2aPV2b
NRMS E (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)
Standard8.286.430.9914.81−13.930.988.27−4.870.99
Skoplaki3.01−0.871.009.034.980.9616.3613.400.99
Skoplaki 124.2022.480.987.77−5.030.986.164.370.97
Skoplaki 252.0148.560.9613.929.790.9424.2218.720.98
Faiman/Koehl9.64−9.321.0023.91−23.310.9815.47−12.650.98
Mattei 129.87−27.840.9948.43−47.280.7747.33−44.930.93
Table 10. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in August, 2019.
Table 10. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in August, 2019.
PV1PV2aPV2b
NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)
Standard7.163.670.9817.95−17.030.985.04−3.760.99
Skoplaki2.90−0.461.0012.738.340.9718.2815.990.98
Skoplaki 128.0825.000.947.43−3.020.929.237.060.93
Skoplaki 262.0455.840.9622.9417.180.9626.3222.750.96
Faiman/Koehl9.74−9.501.0026.63−25.570.9711.41−10.560.98
Mattei 127.44−24.500.9846.13−43.520.9343.80−40.680.82
Table 11. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in September, 2019.
Table 11. The dataset of the NRMSE, NMBE values and the Pearson’s correlation coefficient for the average daily values of the module temperatures for the PV1, PV2a and PV2b installations calculated with the use of mathematical models in September, 2019.
PV1PV2aPV2b
NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)NRMSE (%)NMBE (%)k (-)
Standard11.258.780.9814.56−12.700.976.44−1.830.97
Skoplaki3.75−1.981.0016.349.280.9619.4915.770.97
Skoplaki 126.1421.690.569.96−3.650.599.235.230.66
Skoplaki 254.7846.420.9622.2313.590.9323.9219.150.96
Faiman/Koehl11.32−10.681.0028.34−26.210.9414.62−13.080.97
Mattei 121.84−14.640.9739.54−34.610.8937.03−31.680.72
Table 12. The dataset of the NRMSE values for the average daily and monthly values of the module temperatures for the PV2a and PV2b installations calculated with the use of the Skoplaki model for mounting coefficient ω = 2.4 from the literature and for ω = 2 proposed by the authors.
Table 12. The dataset of the NRMSE values for the average daily and monthly values of the module temperatures for the PV2a and PV2b installations calculated with the use of the Skoplaki model for mounting coefficient ω = 2.4 from the literature and for ω = 2 proposed by the authors.
PV2aPV2b
NRMSE (%)NRMSE (%)NRMSE (%)NRMSE (%)
ω ω = 2 ω = 2.4 ω = 2
April14.497.8120.3410.54
June9.036.6716.367.94
August12.736.8618.2810.64
September16.349.7219.4911.36
whole year 20197.958.4116.9310.24
Table 13. The dataset of the NMBE values for the average daily and monthly values of the module temperatures for the PV2a and PV2b installations calculated with the use of the Skoplaki model for mounting coefficient ω = 2.4 from the literature and for ω = 2 proposed by the authors.
Table 13. The dataset of the NMBE values for the average daily and monthly values of the module temperatures for the PV2a and PV2b installations calculated with the use of the Skoplaki model for mounting coefficient ω = 2.4 from the literature and for ω = 2 proposed by the authors.
PV2aPV2b
NMBE (%)NMBE (%)NMBE (%)NMBE (%)
ω = 2.4 ω = 2 ω = 2.4 ω = 2
April10.30−1.6917.887.07
June4.98−2.8213.406.22
August8.34−0.9915.998.68
September9.28−0.5315.777.80
whole year 20191.92−7.4310.492.38
Table 14. The differences between the average observed and estimated module temperatures at two average wind speeds on 24 August 2019 and 28 August 2019 (from 8 a.m. to 4 p.m.).
Table 14. The differences between the average observed and estimated module temperatures at two average wind speeds on 24 August 2019 and 28 August 2019 (from 8 a.m. to 4 p.m.).
Installation Name and DateFaimanSkoplakiSkoplaki*Skoplaki 1Mattei 1Wind SpeedModule and Ambient Temp.Irradiance
T ¯ d i f f
(°C)
T ¯ d i f f
(°C)
T ¯ d i f f
(°C)
T ¯ d i f f
(°C)
T ¯ d i f f
(°C)
v ¯ w
(m/s)
T ¯ m o d M ; T ¯ a
(°C)
G ¯ 0
(W/m2)
PV1: 24 Aug.4.3−1.4−1.4−18.417.61.047.9; 22.8755.5
PV1: 28 Aug.3.1−1.1−1.1−12.613.62.744.5; 25.4750.5
PV2a: 24 Aug.13.1−9.0−3.0 ( ω   = 2)−2.526.81.050.1; 22.8518.6
PV2a: 28 Aug.9.9−7.5−2.7 ( ω   = 2)−1.221.92.746.8; 25.4534.7
PV2b: 24 Aug.3.1−10.2−6.6 ( ω   = 2)−6.317.21.034.4; 22.8311.8
PV2b: 28 Aug.5.0−5.0−2.2 ( ω   = 2)−1.418.62.737.0; 25.4315.7

Share and Cite

MDPI and ACS Style

Idzkowski, A.; Karasowska, K.; Walendziuk, W. Temperature Analysis of the Stand-Alone and Building Integrated Photovoltaic Systems Based on Simulation and Measurement Data. Energies 2020, 13, 4274. https://doi.org/10.3390/en13164274

AMA Style

Idzkowski A, Karasowska K, Walendziuk W. Temperature Analysis of the Stand-Alone and Building Integrated Photovoltaic Systems Based on Simulation and Measurement Data. Energies. 2020; 13(16):4274. https://doi.org/10.3390/en13164274

Chicago/Turabian Style

Idzkowski, Adam, Karolina Karasowska, and Wojciech Walendziuk. 2020. "Temperature Analysis of the Stand-Alone and Building Integrated Photovoltaic Systems Based on Simulation and Measurement Data" Energies 13, no. 16: 4274. https://doi.org/10.3390/en13164274

APA Style

Idzkowski, A., Karasowska, K., & Walendziuk, W. (2020). Temperature Analysis of the Stand-Alone and Building Integrated Photovoltaic Systems Based on Simulation and Measurement Data. Energies, 13(16), 4274. https://doi.org/10.3390/en13164274

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop