It is prerequisite to predict the behaviour of photovoltaic (PV) modules in a particular geographical area where the system is to be installed for their better performance and increasing lifetime. For that, models are the easiest and acceptable tools to characterise the behaviour of PV modules in any location. The purpose of this study was to develop an empirical model to predict the influence of temperature on the performance of four different PV module technologies, namely, polycrystalline, monocrystalline, amorphous, and thin film in an outdoor environment. The model has been developed by fitting of one year experimental data using the least squares method. The estimated results of the developed model were validated with real-time data (winter and summer season) and a comparison of other existing model estimates using error analysis with 95% confidence interval. The proposed model estimations confirm that the monocrystalline module performs better in winter and polycrystalline in summer as compared to amorphous and thin film in the study area. During analysis, it is revealed that developed model results are more precise and appropriate among other existing model estimations. It is concluded that the proposed model estimations could be used for the prediction of PV module temperature in similar environmental conditions as that of the study area with more accuracy and confidence. It ultimately helps to develop cost-effective and efficient PV systems.
Quaid-e-Awam University of Engineering, Science and Technology1. Introduction
The intensities of solar radiation, ambient temperature, wind speed, relative humidity, configuration, and method of mounting are considered to be responsible for variations in the power output of photovoltaic (PV) modules [1–6]. PV module temperature is one of the key parameters which affect the performance of photovoltaic (PV) modules after solar radiations [2, 4, 7, 8]. Photovoltaic power output is proportional to the PV module operating temperature [9, 10]. Since the change of PV module temperature depends on the variation of ambient temperature, as ambient temperature increases, the module temperature increases and vice versa [4, 11]. It is because, the increase of temperature reduces the band gap of a PV module and increases the energy of the electrons in the material, which ultimately increases the recombination rate of internal carriers caused by the increasing amount of carrier concentrations [9, 12–14]. Consequently, it slightly increases the short-circuit current and considerably decreases open-circuit voltage [2, 9, 15]. Weather conditions affect the PV module temperature; therefore, its influence is necessary to be quantified. This can be done with the help of modeling, which eventually helps to design better systems for proper functioning. Several attempts have been made by different authors from different countries to exemplify the behaviour of PV modules. Some models are intuitive, and others are analytical, numerical, or empirical. Nevertheless, the majority of models are validated in indoor environments of developed countries with the exception of a few in outdoor conditions. It is a very challenging task to develop a model which represents the behaviour of various module technologies simultaneously in outdoor environments. An exact module temperature estimation model is indispensable to achieve reliable data of PV module power output [5, 11, 16–22]. The models used for the prediction of module temperature can be categorized in different ways: steady-state or dynamic, explicit or implicit, etc. [2, 9, 11]. In steady-state modeling, all parameters are assumed to be independent of time (with small time interval, i.e., an hour). However, such models are useful for specific locations and module technologies, while in the dynamic models, some parameters are considered to be varied with respect to time. Dynamic models are preferable for high-resolution input data. Explicit models predict the value of photovoltaic module temperature directly, whereas the implicit correlations involve variables that themselves depend on module temperature. In implicit models, an iteration procedure is compulsory to get the outputs [2, 5, 9, 23–30]. Nevertheless, the selection of an appropriate model is crucial for the design and sizing of photovoltaic systems. The use of an inappropriate model gives faulty predictions thus making the systems over- or undersized. The oversized system becomes costly alternative, whereas undersizing causes malfunctioning of the system. This problem can be controlled through proper sizing and designing of system components with the help of precise modeling and using of long-term reliable data [9, 20, 31-34]. Unfortunately, long-term data are not available in developing countries [31] including Pakistan [35], and the reliability of data is also questionable. Actually, photovoltaic module temperature models are submodels of power output models, as these models predict the effect of temperature on the performance of photovoltaic modules. Most of such models estimate the temperature of photovoltaic modules in indoor conditions but not in outdoor environments [36–38]. The main objective of this study was to develop a simple empirical model for the estimation of the temperature effect on four different PV module technologies, namely, polycrystalline, monocrystalline, amorphous, and thin film in an outdoor environment.
2. Existing Photovoltaic Module Temperature Models
In [39], the researchers consider only one basic climatic variable such as the ambient temperature (Ta) in their study. It is clear that one input variable does not reflect the whole behave of the environment. The developed model is given in equation (1) and also used by [40].
(1)Tm=1.411×Ta−6.414.
Muzathik [38] suggested three variable models with ambient temperature Ta (°C), global solar radiation Gsr (W/m2), and wind speed Wv (m/s). The model and coefficients of each variable are provided as given in Equation 2.
(2)Tm=0.943×Ta+0.0195×Gsr−1.528×Wv+0.3529.
In addition, [2] proposed a simple and semiempirical model for the calculation of module temperature as given in equation (3). The author considered Ta in (°C), Gsr in (W/m2), and Wv in (m/s). The same model is reported by [41].
(3)Tm=Ta+0.255.7+3.8×Wv×Gsr.
Duffie and Beckman [42] proposed a novel mathematical approach for the calculation of photovoltaic module temperature in controlled nominal operating cell temperature (NOCT) conditions: 0.8 kW/m2 solar radiation, 20°C ambient temperature, and 1 m/s wind speed. The model depends on the input of Ta (°C), Gsr (W/m2), Wv (m/s), and NOCT conditions as given in equation (4). Furthermore, the model is adopted by [9].
(4)Tm=Ta+9.55.7+3.8×WvGsrGsr−NOCTTm−NOCT−Ta−NOCT1−ηmτα.
Risser and Fuentes [43] also proposed three variable models with the same variables as that of Muzathik [38] as given in equation 5. The author considered Ta in (°C), Gsr in (W/m2), and Wv in (m/s). The same model is tested by [19].
(5)Tm=1.31×Ta+0.0282×Gsr−1.65×Wv+3.81.
The authors [2, 38, 43] proposed new temperature models which were based on three basic input variables (solar radiations “Gsr,” ambient temperature “Ta,” and wind speed “Wv”). The researchers proposed linear models in their studies, but the behaviour of climatic data is parabolic with respect to time. In the morning hours, the intensities of Gsr and Ta are directly proportional, but in the evening, these are less related due to the slight decreasing trend of temperature as compared to the sharp decrease of solar radiations. The authors [9, 42] proposed a mathematical approach for the calculation of photovoltaic module temperature based on NOCT conditions. Such conditions could not be familiarized with a real outdoor condition.
Almaktar et al. [40] proposed a temperature model which depends on four climatic variables, namely, solar radiations “Gsr,” ambient temperature “Ta,” wind speed “Wv,” and relative humidity “Rh” as given in Equation 6.
(6)Tm=0.77×Ta+0.023×Gsr−0.137×Wv−0.206×Rh+26.97.
It is already mentioned that the behaviour of climatic data is parabolic in nature with respect to the time of the day. Therefore, in this study, an empirical, nonlinear, multivariate, and least squares model was developed and proposed to calculate the PV module temperature in an outdoor environment. Table 1 shows the well-known PV module temperature models.
The study was conducted in Nawabshah city, Shaheed Benazirabad District, Sindh, Pakistan, as shown in Figure 1. It is one of the hottest places and located at 26.14°N, 68.23°E [44] and mean 37 m above sea level [45].
Geographical location of the study area (Google Maps).
Map of Pakistan
Study area
3.2. Experimental Setup
An experimental setup was installed at the Energy and Environment Engineering Department, QUEST, Nawabshah. Four generic photovoltaic modules (polycrystalline, monocrystalline, amorphous, and thin film) were used in this study, and their specifications are given in Table 2. The photovoltaic modules were fixed on an iron structure, facing true south at an inclination of 12° to the horizontal plane. Figure 2 shows the schematic diagram and experimental setup. The data of each PV module was recorded for a whole year from November, 2015, to October, 2016, for the development of a suitable model for the prediction of module temperature. The data was measured at an interval of 1 hour from 07 to 18 hours daily. Global solar radiation Gsr (W/m2), ambient temperature Ta (°C), wind speed Wv (m/s), and relative humidity Rh (%) were measured with HP-2000. Photovoltaic module temperature was recorded with the help of Prova-830 (8 channel thermocouple data logger). A total of eight numbers (two on each PV module) of k-type thermocouples were pasted on the surface and backside of the photovoltaic modules as shown in Figure 3 [14, 46]. The accuracy of equipment used for data measurement is given in Table 3. The temperature sensors were pasted on the surface and backside of photovoltaic module as that of [8], and then, the average of temperature was taken as the module operating temperature [14, 46, 47].
Electrical characteristics of examined photovoltaic modules [14, 46].
Module parameters
Unit
Module technologies
(p-Si)
(m-Si)
(a-Si)
Thin film
SUN-40P
SUN-40M
TPS-40
GS-50
Voc
V
43
21.5
29
62
Isc
A
1.29
2.55
2.3
1.42
Vmax
V
35
17.5
18
43
Imax
A
1.14
2.29
2.2
1.17
Pmax
W
40
40
40
50
Area of module
m2
0.27
0.24
0.76
0.74
Experimental setup.
Schematic flow diagram of experimental setup
Experimental setup for data logging
Pasted thermocouple sensors on the surface and backside of each photovoltaic module.
Pasted thermocouples on surface and backside of each PV modules
Positions of thermocouple sensors (surface and backside of PV modules)
Accuracy of equipment used in this study.
Parameters
Unit
Weather station (HP-2000)
Module temperature recorder (Prova-830)
Accuracy
Accuracy
Gsr
W/m2
±15%
—
Ta
°C
±1.0°C
—
Wv
m/s
±1 m/s (windspeed<5) ±10% (windspeed>5)
—
Rh
%
±5%
—
Tm
°C
—
±0.1% or 1.0°C
4. Proposed Empirical Model
In this section, we develop a model for the estimation of the temperature effect on different photovoltaic (PV) module technologies, namely, polycrystalline, monocrystalline, amorphous, and thin film in the outdoor environment. In model development, one dependent variable (module temperature) and four basic independent climatic variables (global solar radiation, ambient temperature, wind speed, and relative humidity) were adopted. Furthermore, the correlation of the dependent variable with each independent variable was analyzed. The correlation of module temperature (Tm) with the global solar radiation (Gsr) was found to be 0.89217, ambient temperature (Ta) 0.73765, wind speed (Wv) 0.075766, and relative humidity (Rh) -0.55918. The relationship between climatic parameters and module temperature was found to be nonlinear because of the parabolic curve. Thus, it was deduced from the curve fitting that polynomial models might be suitable models, as these cover the maximum number of measured data points. Further scrutiny of models was made by fitting the data with different degrees of polynomials (1-9 degrees). It was found that the 2nd degree polynomial model covers the maximum number of data points of the measured data. Thus, an empirical second degree multivariate nonlinear model was proposed with fitting of data with the least squares method. It was assumed that photovoltaic module temperature (Tm) is the function of four variables, namely, Gsr, Ta, Wv, and Rh. Thus, the basic function of PV module temperature (Tm) is given in equation (7).
(7)Tm=fGsr,Ta,Wv,Rh.
The general form of the model would be given in Equation 8.
(8)Tm=a1Gsr+a2Ta+a3Wv+a4Rh+a5×b1Gsr+b2Ta+b3Wv+b4Rh+b5.
By expanding equation (8) with the combination of all four independent variables, equation (9) is developed, which demonstrates the output and input parameters and all involved coefficients.
(9)Tm=a1Gsrb2Ta+a1Gsrb3Wv+a1Gsrb4Rh+a2Tab1Gsr+a2Tab3Wv+a2Tab4Rh+a3Wvb1Gsr+a3Wvb2Ta+a3Wvb4Rh+a4Rhb1Gsr+a4Rhb2Ta+a4Rhb3Wv+a3Wv2b3+a5b4Rh+a3Wvb5+a2Tab5+a5b3Wv+a4Rh2b4+a1Gsrb5+a5b1Gsr+a4Rhb5+a5b2Ta+a2Ta2b2+a1Gsr2b1+a5b5.
Let Tm_meas be the measured module temperature and Tm_est be the estimated module temperature. The least squares method assumes that the sum of the squares of the residuals (error) is less. Therefore, it can be estimated using Equation 10.
(10)Ei=min∑i=1nTm_measi−Tm_estiGsri,Tai,Wvi,Rhi;β2.where i=1,2,⋯,n, as n=4392 and β is the set of the coefficients of the model. The minimum value of E occurs when the gradient is zero. The model contains m=25 parameters; therefore, the gradient equation is 25. Furthermore, the minimum values of E and ri are calculated through equations (11) and (12).
(11)∂Ei∂βj=2∑i=1nri∂ri∂βj=0.where j=1,2,⋯,m=25.
(12)ri=Tm_measi−Tm_estiGsri,Tai,Wvi,Rhi;β.
The model equation (12) is complex and time-consuming. Thus, it requires to be simplified for easy computation and application. For that, symbolic derivatives of equation 12 were put in MAPLE software, by producing a system of equations with the coefficients βj. Then, the obtained system of equations from MAPLE software was solved iteratively in MATLAB software. The coefficients of the developed model were approximated with an error tolerance of 0.0001. The general form of the developed model for the estimation of all four types of photovoltaic (PV) module temperatures is shown in equation (13), and model coefficients are given in Table 4.
(13)Tm=α0+β1Gsr+β2Gsr2+γ1Ta+γ2Ta2+δGsrTa+λWv+ζRh.where α0, β1, β2, γ1, γ2, δ, λ, and ζ are least squares coefficients of the proposed model.
Proposed model coefficients.
Model coefficients
Photovoltaic module technologies
Polycrystalline
Monocrystalline
Amorphous
Thin film
α
22.5505
31.3750
33.9800
32.4500
β1
0.03753
0.03858
0.03622
0.03340
β2
−5.71×10−7
−1.91×10−6
0.00000
−1.974×10−6
γ1
0.005892
0.6672
0.1191
0.2982
γ2
0.01179
0.0000
0.01078
0.007552
δ
-0.0002703
-0.0002805
-0.000245
-0.0001666
λ
-0.6070
-6.4460
-5.0350
-4.9540
ζ
-0.0960
-0.2100
-0.1691
-0.1935
5. Statistical Analysis
Statistical analysis was conducted to see the variation between models’ estimated and measured results. The coefficient of determination (R2) [48, 49], root mean square error (RMSE), and mean absolute error (MAE) [40, 48–50] were used as statistical indicators as given in equation (14), respectively. The root mean square error (RMSE) and mean absolute error (MAE) are considered in °C. The statistical analysis was done at 95% confidence level.
(14)R2=∑i=1nTm_esti−T¯m_esti2∑i=1nTm_measi−T¯m_measi2,RMSE=∑i=1nTm_esti−Tm_measi2n,MAE=1n∑i=1nTm_esti−Tm_measi.where T¯m_est is the average estimated module temperature and T¯m_meas is the average measured module temperature.
6. Results and Discussion6.1. Weather Conditions
The average hourly global solar radiation (Gsr), maximum and minimum ambient temperature (Ta), wind speed (Wv), and relative humidity (Rh) of a whole year from November, 2015, to October, 2016, are shown in Figures 4–7. The yearly average total global solar radiations were found to be 6224.35 kWh/m2/day with a maximum average of 835.25 W/m2 at 12 hours and a minimum average of 86.02 W/m2 at 07 hours. The values of global solar radiations are given in Figure 4. The maximum Ta was noted as 34.67°C at 15 hours and the minimum as 21.35°C at 07 hours with a yearly average of 30.11°C during the study period. The ambient temperature values are shown in Figure 5. Similarly, Figure 6 displays the wind speed. The maximum yearly average Wv was recorded as 2.60 m/s at 16 hours and the minimum as 1.30 m/s at 07 hours with a yearly average of 2.14 m/s. Likewise, the maximum yearly average Rh was noted as 76.90% at 07 hours and the minimum as 26.25% with a yearly average of 42.66%. The Rh is given in Figure 7. The yearly average values of climatic conditions like Gsr, Ta, Wv, and Rh are given in Table 5. Relative humidity was found inversely proportional to the intensity of global solar radiation and ambient temperature.
Yearly hourly average values of global solar radiation.
Yearly hourly average values of ambient temperature.
Yearly hourly average values of wind speed.
Yearly hourly average values of relative humidity.
Yearly average values of climatic conditions.
Gsr (W/m2)
Ta (°C)
Wv (m/s)
Rh (%)
Maximum values
835.25
34.67
2.60
76.90
Minimum values
86.02
21.35
1.30
26.25
Average values
518.69
30.11
2.14
42.66
6.2. Validation of Proposed Model Results
The proposed model results were validated by comparing its estimations with measured data of winter season (for the months of December and January 2017) and summer season (for the months of May and June 2017) and with estimations of other existing models.
6.3. Proposed and Existing Model Estimations versus Measured Data of Winter Season
The comparison of the proposed model estimation versus measured and other existing model results of polycrystalline, monocrystalline, amorphous, and thin film modules are shown in Figures 8–11 for winter season, respectively. The average proposed model has estimated 0.19°C (0.61%), 0.48°C (1.61%), 0.06°C (0.19%), and 0.07°C (0.24%) low module temperatures for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively, than measured ones. It was found that monocrystalline has a least percentage of module temperature than other module technologies. It was revealed that Rahman et al. [39], Muzathik [38], Skoplaki et al. [2] and Duffie and Beckman [42] models’ predicted results are lower than those of the measured and proposed models’ estimated module temperature because of the lower number of input parameters. Risser and Fuentes [43] and Almaktar et al. [40] models gave higher average module temperature of 3.40°C (10.91%) and 6.32°C (16.83%) for polycrystalline, 4.61°C (15.38%) and 7.53°C (25.09%) for monocrystalline, 3.34°C (10.68%) and 6.25°C (20.00%) for amorphous, and 4.62°C (15.43%) and 7.54°C (25.14%) for thin film modules, respectively, than the proposed model estimations.
Proposed model estimation versus measured data and other existing model module temperature values of polycrystalline module during winter.
Proposed model estimation versus measured data and other existing model module temperature values of monocrystalline module during winter.
Proposed model estimation versus measured data and other existing model module temperature values of amorphous module during winter.
Proposed model estimation versus measured data and other existing model module temperature values of thin film module during winter.
6.4. Proposed and Existing Model Estimations versus Measured Data of Summer Season
The comparison of the proposed model estimation versus measured and other existing model results of polycrystalline, monocrystalline, amorphous, and thin film modules are shown in Figures 12–15 for summer season, respectively. The proposed model gave 0.43°C (0.84%) higher module temperature for polycrystalline and 1.31°C (2.65%), 0.90°C (1.76%), and 1.07°C (2.15%) lower module temperature for monocrystalline, amorphous, and thin film modules than the measured module temperature. It was found that monocrystalline estimates a least percentage of module temperature than amorphous and thin film modules. It was found that Rahman et al. [39], Muzathik [38], Skoplaki et al. [2], and Duffie and Beckman [42] models’ predicted results are lower than those of the measured and proposed models’ estimated values. Risser and Fuentes [43] and Almaktar et al. [40] models gave higher module temperatures of 15.79°C (30.42%) and 10.21°C (19.68%) for polycrystalline, 18.41°C (37.35%) and 12.83°C (26.03%) for monocrystalline, 16.64°C (32.60%) and 11.07°C (21.68%) for amorphous, and 18.04°C (36.32%) and 12.46°C (25.09%) for thin film modules, respectively, than the proposed model’s estimated values. It was found that the proposed model estimates a low temperature with 1.61% in winter and 2.65% in summer from monocrystalline than other measured modules.
Proposed model estimation versus measured data and other existing model module temperature values of polycrystalline module during summer.
Proposed model estimation versus measured data and other existing model module temperature values of monocrystalline module during summer.
Proposed model estimation versus measured data and other existing model module temperature values of amorphous module during summer.
Proposed model estimation versus measured data and other existing model module temperature values of thin film module during summer.
7. Error Analysis
The error analysis of the proposed model estimations was checked with measured data of winter season and summer season and with estimations of other existing models. The coefficient of determination (R2), root mean square error (RMSE) (°C), and mean absolute error (MAE) (°C) of each PV module of winter season (months of December and January) are summarized in Tables 6–8 and of the season of summer (months of May and June) in Tables 9–11, respectively.
Coefficient of determination (R2) of each PV module in winter.
S. no.
Name of model
Coefficient of determination (R2)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.996
0.998
0.992
0.994
(2)
Rahman et al. [39]
0.646
0.707
0.650
0.725
(3)
Muzathik [38]
0.970
0.982
0.968
0.986
(4)
Skoplaki et al. [2]
0.946
0.965
0.942
0.969
(5)
Duffie and Beckman [42]
0.955
0.970
0.949
0.974
(6)
Risser and Fuentes [43]
0.972
0.982
0.969
0.986
(7)
Almaktar et al. [40]
0.924
0.953
0.925
0.962
Root mean square error (RMSE) of each PV module in winter.
S. no.
Name of model
Root mean square error (RMSE) (°C)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.955
0.673
0.898
0.763
(2)
Rahman et al. [39]
12.357
11.009
12.177
10.809
(3)
Muzathik [38]
9.935
8.678
9.854
8.560
(4)
Skoplaki et al. [2]
6.992
5.713
6.900
5.560
(5)
Duffie and Beckman [42]
6.468
5.194
6.407
5.057
(6)
Risser and Fuentes [43]
4.104
5.093
7.566
5.243
(7)
Almaktar et al. [40]
7.335
8.243
4.481
8.344
Mean absolute error (MAE) of each PV module in winter.
S. no.
Name of model
Mean absolute error (MAE) (°C)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.782
0.636
0.721
0.666
(2)
Rahman et al. [39]
11.165
9.884
11.069
9.813
(3)
Muzathik [38]
9.636
8.426
9.701
8.413
(4)
Skoplaki et al. [2]
6.512
5.302
6.576
5.289
(5)
Duffie and Beckman [42]
6.047
4.838
6.112
4.825
(6)
Risser and Fuentes [43]
3.407
4.616
3.987
4.685
(7)
Almaktar et al. [40]
6.322
7.532
6.543
7.544
Coefficient of determination (R2) of each PV module in summer.
S. no.
Name of model
Coefficient of determination (R2)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.996
0.995
0.993
0.992
(2)
Rahman et al. [39]
0.496
0.402
0.477
0.476
(3)
Muzathik [38]
0.991
0.980
0.990
0.984
(4)
Skoplaki et al. [2]
0.976
0.943
0.974
0.971
(5)
Duffie and Beckman [42]
0.983
0.955
0.983
0.979
(6)
Almaktar et al. [40]
0.982
0.946
0.973
0.974
(7)
Risser and Fuentes [43]
0.981
0.980
0.984
0.984
Root mean square error (RMSE) of each PV module in summer.
S. no.
Name of model
Root mean square error (RMSE) (°C)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.996
1.330
1.262
1.502
(2)
Rahman et al. [39]
8.458
7.179
8.137
7.292
(3)
Muzathik [38]
7.444
4.790
6.665
5.234
(4)
Skoplaki et al. [2]
5.150
3.172
4.590
3.530
(5)
Duffie and Beckman [42]
4.363
2.655
3.853
2.952
(6)
Almaktar et al. [40]
10.321
13.114
11.209
12.596
(7)
Risser and Fuentes [43]
15.833
18.564
16.694
18.107
Mean absolute error (MAE) of each PV module in summer.
S. no.
Name of model
Mean absolute error (MAE) (°C)
(p-Si)
(m-Si)
(a-Si)
Thin film
(1)
Proposed model
0.832
1.176
1.078
1.180
(2)
Rahman et al. [39]
7.603
6.363
7.286
6.434
(3)
Muzathik [38]
7.053
4.436
6.200
4.808
(4)
Skoplaki et al. [2]
4.447
2.813
3.954
3.147
(5)
Duffie and Beckman [42]
3.799
2.239
3.334
2.646
(6)
Almaktar et al. [40]
10.218
12.835
11.071
12.463
(7)
Risser and Fuentes [43]
15.796
18.413
16.649
18.041
In winter season, the maximum R2 was given by the proposed model with 0.996, 0.998, 0.992, and 0.994 and the minimum by Rahman et al. [39] model with 0.646, 0.707, 0.650, and 0.725 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively. Similarly, the minimum RMSE was noted by the proposed model with 0.955, 0.673, 0.898, and 0.763 and the maximum by Rahman et al. [39] model with 12.357, 11.009, 12.177, and 10.809 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively, than other existing models. Likewise, the minimum MAE was noted by the proposed model with 0.782, 0.636, 0.721, and 0.666 and the maximum by Rahman et al. [39] model with 11.165, 9.884, 11.069, and 9.813 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively, than other existing models.
In summer season, the maximum R2 was given by the proposed model with 0.996, 0.995, 0.993, and 0.992 and the minimum by Rahman et al. [39] model with 0.496, 0.402, 0.477, and 0.456 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively. Similarly, the minimum RMSE was noted by the proposed model with 0.996, 1.330, 1.262, and 1.502 and the maximum by Risser and Fuentes [43] model with 15.833, 18.564, 16.694, and 18.107 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively, than other existing models. Likewise, the minimum MAE was noted by the proposed model with 0.832, 1.176, 1.078, and 1.180 and the maximum by Risser and Fuentes [43] model with 15.796, 18.413, 16.649, and 18.041 for polycrystalline, monocrystalline, amorphous, and thin film modules, respectively, than other existing models.
The proposed model gave the maximum R2 and minimum RMSE and MAE than other existing model estimations. Thus, the proposed model results are more appropriate than other existing model results. It was observed that the models of one input variable show the minimum coefficient of determination and maximum root mean square error as well as the mean absolute error and vice versa for more input variable models.
8. Conclusions
Photovoltaic (PV) operating temperature plays an important role in the PV conversion process after solar radiation. It is a very challenging task to develop a model which represents the behaviour of various module technologies in outdoor environments. An empirical second degree polynomial multivariate model was developed using the least squares data fitting method to estimate the module temperature in outdoor conditions. It was validated by comparing the proposed model estimations with real-time measured data of winter and summer season and other existing model estimations through error analysis. It was revealed that the proposed model estimated the least temperature for monocrystalline module with 0.48°C (1.61%) in winter season and 1.31°C (2.66%) in summer season than other examined module technologies. Risser and Fuentes [43] and Almaktar et al. [40] models gave a higher average of module temperature with 4.62°C (15.43%) and 7.54°C (25.14%) for thin film in winter and 18.41°C (37.35%) and 12.83°C (26.03%) for monocrystalline modules in summer than the proposed model estimated values. It was found that Rahman et al. [39] model shows the least behavior of module temperature than the measured, proposed model estimations and other existing model estimation values in both seasons. The proposed model gave around 0.998 coefficient of determination for monocrystalline and low root mean square error and mean absolute error in both seasons. It is concluded that the proposed model is more appropriate for the estimation of photovoltaic module temperature in outdoor conditions because the proposed model gave a maximum coefficient of determination and minimum root mean square error and mean absolute error in both seasons. It is recommended that the time interval of data recording may be reduced from 1 hour to minutes and PV module technologies with the same ratings may be used for a comparison purpose. The performance and effect of temperature on both free standing and building integrated systems may be checked and verified in outdoor environments.
NomenclatureGsr:
Global solar radiation (W/m2)
Ta:
Ambient temperature (°C)
Wv:
Wind speed (m/s)
Rh:
Relative humidity (%)
Tm:
Module temperature (°C)
°C:
Degree centigrade
m/s:
Meter per second
%:
Percentage
PV:
Photovoltaic
p-Si:
Polycrystalline
m-Si:
Monocrystalline
a-Si:
Amorphous
Voc:
Open-circuit voltage (V)
Isc:
Short-circuit current (A)
Vmax:
Maximum voltage (V)
Imax:
Maximum current (A)
Pmax:
Maximum power (W)
kW/m²:
Kilowatt per square meter
kWh/m²/d:
Kilowatt hour per square meter per day
W/m²:
Watt per square meter
mW/cm2:
Milliwatts per square centimeter
Tm:
Module temperature (°C)
Tm_meas:
Measured module temperature (°C)
Tm_est:
Estimated module temperature (°C)
Tmp_meas:
Polycrystalline measured module temperature (°C)
Tmp_est:
Polycrystalline estimated module temperature (°C)
Tmm_meas:
Monocrystalline measured module temperature (°C)
Tmm_est:
Monocrystalline estimated module temperature (°C)
Tma_meas:
Amorphous measured module temperature (°C)
Tma_est:
Amorphous estimated module temperature (°C)
Tmt_meas:
Thin film measured module temperature (°C)
Tmt_est:
Thin film estimated module temperature (°C)
m:
Meter
m²:
Meter square
ηm:
Efficiency of module
τα:
Transmittance of glass
kr:
Ross coefficient
R2:
Coefficient of determination
RMSE:
Root mean square error (°C)
MAE:
Mean absolute error (°C)
SUN:
Sanak SK Union
TPS:
Topray Solar, Shenzhen
NOCT:
Nominal operating cell temperature
STC:
Standard test conditions
QUEST:
Quaid-e-Awam University of Engineering, Science and Technology
Proposed model coefficientα0:
Alpha0 (general coefficients)
β1:
Global solar radiation (Gsr)
β2:
Global solar radiation (Gsr2)
γ1:
Ambient temperature (Ta)
γ2:
Ambient temperature (Ta2)
δ:
Global solar radiation and ambient temperature (Gsr∗Ta)
λ:
Wind speed (Wv)
ξ:
Relative humidity (Rh).
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors have no conflict of interest.
Authors’ Contributions
The authors have worked and contributed equally to this paper. The research article is submitted with the approval of all authors.
Acknowledgments
The authors acknowledge the financial support provided by the Quaid-e-Awam University of Engineering, Science and Technology (QUEST), Nawabshah (QUEST NO. 15).
RouholaminiA.PourgharibshahiH.FadaeinedjadR.AbdolzadehM.Temperature of a photovoltaic module under the influence of different environmental conditions – experimental investigation201637326627210.1080/01430750.2014.9528422-s2.0-84961242076SkoplakiE.BoudouvisA. G.PalyvosJ. A.A simple correlation for the operating temperature of photovoltaic modules of arbitrary mounting200892111393140210.1016/j.solmat.2008.05.0162-s2.0-50249096178SharmaV.ChandelS. S.Performance and degradation analysis for long term reliability of solar photovoltaic systems: a review20132775376710.1016/j.rser.2013.07.0462-s2.0-84882260617JakhraniA. Q.OthmanA. K.RigitA. R. H.SamoS. R.Comparison of solar photovoltaic module temperature models201114318BarykinaE.HammerA.Modeling of photovoltaic module temperature using Faiman model: sensitivity analysis for different climates201714640141610.1016/j.solener.2017.03.0022-s2.0-85014806752MeralM. E.DincerF.A review of the factors affecting operation and efficiency of photovoltaic based electricity generation systems20111552176218410.1016/j.rser.2011.01.0102-s2.0-79952613274Luketa-HanlinA.SteinJ.2012SAND2012-4307CSandia National LaboratoriesMarinić-KragićI.NižetićS.Grubišić-ČaboF.PapadopoulosA. M.Analysis of flow separation effect in the case of the free-standing photovoltaic panel exposed to various operating conditions2018174536410.1016/j.jclepro.2017.10.3102-s2.0-85038872536JakhraniA. Q.OthmanA. K.RigitA. R. H.SamoS. R.Determination and comparison of different photovoltaic module temperature models for Kuching, Sarawak2011 IEEE Conference on Clean Energy and Technology (CET)June 2011Kuala Lumpur, Malaysia231236LeeY.TayA. A. O.Finite element thermal analysis of a solar photovoltaic module20121541342010.1016/j.egypro.2012.02.0502-s2.0-84860525274SkoplakiE.PalyvosJ. A.On the temperature dependence of photovoltaic module electrical performance: a review of efficiency/power correlations200983561462410.1016/j.solener.2008.10.0082-s2.0-63449141127YanH.WangX.YaoM.YaoX.Band structure design of semiconductors for enhanced photocatalytic activity: the case of TiO2201323440240710.1016/j.pnsc.2013.06.0022-s2.0-84881465957KalogirouS. A.20142nd EditionUSAAcademic Press ElsevierJatoiA. R.SamoS. R.JakhraniA. Q.Influence of temperature on electrical characteristics of different photovoltaic module technologies201872859110.14710/ijred.7.2.85-91JakhraniA. Q.OthmanA. K.RigitA. R. H.BainiR.SamoS. R.LingL. P.2012NED University Journal of Research, Thematic Issue on EnergySharanS. N.MathurS. S.KandpalT. C.Analytical performance evaluation of combined photovoltaic-thermal concentrator-receiver systems with linear absorbers198727436136510.1016/0196-8904(87)90113-02-s2.0-0023597023SandnesB.RekstadJ.A photovoltaic/thermal (PV/T) collector with a polymer absorber plate: experimental study and analytical model2002721637310.1016/S0038-092X(01)00091-32-s2.0-0037122465ChowT. T.Performance analysis of photovoltaic-thermal collector by explicit dynamic model200375214315210.1016/j.solener.2003.07.0012-s2.0-0142026175SkoplakiE.PalyvosJ. A.Operating temperature of photovoltaic modules: a survey of pertinent correlations2009341232910.1016/j.renene.2008.04.0092-s2.0-51649111368JakhraniA. Q.SamoS. R.KambohS. A.LabadinJ.RigitA. R. H.An improved mathematical model for computing power output of solar photovoltaic modules20142014934670410.1155/2014/3467042-s2.0-84897541330CancroC.FerlitoS.GraditiG.Forecasting the working temperature of a concentrator photovoltaic module by using artificial neural network-based model1766AIP Conference Proceedings2016Freiburg, Germany1090004090007PindadoS.CubasJ.Roibás-MillánE.Bugallo-SiegelF.Sorribes-PalmerF.Assessment of explicit models for different photovoltaic technologies2018116, article 135310.3390/en110613532-s2.0-85049656346MatteiM.NottonG.CristofariC.MuselliM.PoggiP.Calculation of the polycrystalline PV module temperature using a simple method of energy balance200631455356710.1016/j.renene.2005.03.0102-s2.0-28844507884FaimanD.Assessing the outdoor operating temperature of photovoltaic modules200816430731510.1002/pip.8132-s2.0-46749108912SantiagoI.Trillo-MonteroD.Moreno-GarciaI. M.Pallarés-LópezV.Luna-RodríguezJ. J.Modeling of photovoltaic cell temperature losses: a review and a practice case in South Spain201890708910.1016/j.rser.2018.03.0542-s2.0-85044618308JonesA. D.UnderwoodC. P.A thermal model for photovoltaic systems200170434935910.1016/S0038-092X(00)00149-32-s2.0-0035825922ArmstrongS.HurleyW. G.A thermal model for photovoltaic panels under varying atmospheric conditions20103011-121488149510.1016/j.applthermaleng.2010.03.0122-s2.0-77951096926Torres LoberaD.ValkealahtiS.Dynamic thermal model of solar PV systems under varying climatic conditions20139318319410.1016/j.solener.2013.03.0282-s2.0-84877811850JakhraniA. Q.OthmanA. K.RigitA. R. H.SamoS. R.KambohS. A.A novel analytical model for optimal sizing of standalone photovoltaic systems201246167568210.1016/j.energy.2012.05.0202-s2.0-84867229381AlmaktarM.Abdul RahmanH.HassanM. Y.SaehI.Artificial neural network-based photovoltaic module temperature estimation for tropical climate of Malaysia and its impact on photovoltaic system energy yield201523330231810.1002/pip.24242-s2.0-84924075377TekeA.YıldırımH. B.ÇelikÖ.Evaluation and performance comparison of different models for the estimation of solar radiation2015501097110710.1016/j.rser.2015.05.0492-s2.0-84931267297ZhouW.LouC.LiZ.LuL.YangH.Current status of research on optimum sizing of stand-alone hybrid solar-wind power generation systems2010872380389Al-NajideenM. I.AlrwashdehS. S.Design of a solar photovoltaic system to cover the electricity demand for the faculty of Engineering-Mu'tah University in Jordan201734440445AlsadiS.KhatibT.Photovoltaic power systems optimization research status: a review of criteria, constrains, models, techniques, and software tools2018810, article 176110.3390/app81017612-s2.0-85054049446AhmedM. A.AhmedF.AkhtarA. W.Distribution of total and diffuse solar radiation at Lahore, Pakistan20104013743SteinJ. S.CameronC. P.BourneB.KimberA.PosbicJ.JesterT.A standardized approach to PV system performance model validation2010 35th IEEE Photovoltaic Specialists ConferenceJune 2010Honolulu, HI, USA00107900108410.1109/pvsc.2010.56146962-s2.0-78650154788AmmarM. B.ChaabeneM.ChtourouZ.Artificial neural network based control for PV/T panel to track optimum thermal and electrical power20136537238010.1016/j.enconman.2012.08.0032-s2.0-84867429207MuzathikM. A.Photovoltaic modules operating temperature estimation using a simple correlation201444151158RahmanH. A.NorK. M.HassanM. Y.MajidM. S.Empirical models for the correlation of global solar radiation under Malaysia environment20114418641870AlmaktarM.RahmanH. A.HassanM. Y.RahmanS.Climate-based empirical model for PV module temperature estimation in tropical environment201349419220110.3103/S0003701X130400262-s2.0-84892571042CoskunC.ToygarU.SarpdagO.OktayZ.Sensitivity analysis of implicit correlations for photovoltaic module temperature: a review20171641474148510.1016/j.jclepro.2017.07.0802-s2.0-85027517413DuffieJ. A.BeckmanW. A.20063rd EditionHoboken, New JerseyJohn Wiley & Sons, Inc.RisserV. V.FuentesM. K.Linear regression analysis of flat-plate photovoltaic system performance data5th Photovoltaic Solar Energy Conference1984Athens, Greece623627SoomroA.SamoS. R.JatoiA. R.JakhraniA. Q.Energy conservation in office building by utilization of daylighting201610116AhmadW.FatimaA.AwanU. K.AnwarA.Analysis of long term meteorological trends in the middle and lower Indus Basin of Pakistan-a non-parametric statistical approach201412228229110.1016/j.gloplacha.2014.09.0072-s2.0-84908565882JatoiA. R.SamoS. R.JakhraniA. Q.Influence of ambient temperature and solar radiations on photovoltaic module’s temperature and power output20161024347HuangB. J.YangP. E.LinY. P.LinB. Y.ChenH. J.LaiR. C.ChengJ. S.Solar cell junction temperature measurement of PV module201185238839210.1016/j.solener.2010.11.0062-s2.0-78651353394Massi PavanA.MellitA.LughiV.Explicit empirical model for general photovoltaic devices: experimental validation at maximum power point201410110511610.1016/j.solener.2013.12.0242-s2.0-84892449890DeoR. C.ŞahinM.Application of the artificial neural network model for prediction of monthly standardized precipitation and evapotranspiration index using hydrometeorological parameters and climate indices in eastern Australia2015161-162658110.1016/j.atmosres.2015.03.0182-s2.0-84928136722SavvakisN.TsoutsosT.Performance assessment of a thin film photovoltaic system under actual Mediterranean climate conditions in the island of Crete2015901435145510.1016/j.energy.2015.06.0982-s2.0-84959330504