A Study of Developing a Prediction Equation of Electricity Energy Output via Photovoltaic Modules

: Various equations are being developed and applied to predict photovoltaic (PV) module generation. Currently, quite diverse methods for predicting module generation are available, with most equations showing accuracy with ≤ 5% error. However, the accuracy can be determined only when the module temperature and the value of irradiation that reaches the module surface are precisely known. The prediction accuracy of outdoor generation is actually extremely low, as the method for predicting outdoor module temperature has extremely low accuracy. The change in module temperature cannot be predicted accurately because of the real-time change of irradiation and air temperature outdoors. Calculations using conventional equations from other studies show a mean error of temperature difference of 4.23 ◦ C. In this study, an equation was developed and veriﬁed that can predict the precise module temperature up to 1.64 ◦ C, based on the experimental data obtained after installing an actual outdoor module.


Introduction
As the market for photovoltaic modules has increased, studies are being conducted continually to increase the PV module power generation. At the same time, research on PV module power generation prediction has become as important as the research on increasing the efficiency of the PV modules. Most PV modules have a spec mark of PV module power generation in standard test conditions (STC) (25 • C and 1000 W/m 2 ), but actual outdoor installations do not reach this spec. One reason for the failure to reach the spec is irradiation, the factor that has the most influence on PV module power generation. The period during the day when irradiation is maintained at 1000 W/m 2 is extremely short. In addition, module generation does not satisfy the STC standard of actual air temperature of 25 • C, even if the irradiation was 1000 W/m 2 , because of the internal temperature of the module. Even if the air temperature was 25 • C, the irradiation energy that remains after the PV module power generation increases the module temperature [1,2].
Equations that predict PV module power generation are pretty accurate as long as the values of irradiation and module temperature are accurate. An error range of approximately 1-2% with actual power generation is obtained when based on the equation developed by the Sandia National Laboratory of the United States Department of Energy [3]. When this equation is used, the prediction of PV module power generation can be predicted when the irradiation and the temperature of the PV module are known. However, in reality, research on predicting the temperature of modules is lacking. The performance of the module generation system depends on irradiation, outside temperature, module temperature, and wind speed. It is reported that it is affected by changes in meteorological and environmental conditions. The deterioration of the system performance due to the increase in the PV module temperature should be considered for the accurate performance evaluation of the system as well as the process of calculating the capacity and generation amount of PV module generation system. It has been confirmed that research is needed to evaluate the accuracy and validity of predicting PV module temperature using the prediction model in meteorological conditions [4].
Conventional methods for predicting the temperature of actual modules have the following form [5].
where T m : module temperature (K) k: coefficient where In other papers in the same form as Equation (1), equations for predicting temperature have been proposed as Equation (2) [3] and Equation (3) [4].
Regarding Equation (1), coefficients have been suggested by considering variables such as wind speed, outdoor temperature, efficiency, and generation for many predictive models' module temperature. PV module efficiency is influenced by the insolation, air temperature, wind speed, and relative humidity [6]. The prediction of module temperature with respect to irradiation is convenient when the change of module temperature attains a complete thermal equilibrium, when the air temperature and irradiation are maintained for a certain period. However, this is a continuous phenomenon in outdoor conditions. Irradiation and temperature continuously change in real time in outdoor conditions and the sampling of data has to be conducted in a very short time interval (approximately 30 s to 1 min). The existing equations show a large error in predicting real time changes in the module temperature. When the temperature of the module is predicted through the temperature prediction equation in the same form as in actual Equation (1), the form shown in Figure 1 is shown. The paper shown in Figure 1 is a simulation value using the temperature prediction equation of another paper. The maximum difference in measurements and predictions is up to 9.7 • C. Therefore, in this paper, an outdoor testbed was established and an equation was developed that can predict module temperature changes in real time.

Outdoor Testbed Setup
As shown in Figure 2, a testbed with fixed modules was established. Up to eight modules can be installed simultaneously in this testbed. The wind direction, irradiation, air temperature, and module temperature were measured, and the data sampling interval was set to 30 s. The testbed was installed in Gyeongsan-si, South Korea (north latitude 36 • ), with the installation angle being 30 • . An electronic load was used in this study and the maximum power point tracking (MPPT) was continuously traced with an MPPT program. The datamonitoring facility was installed directly adjacent to the test bed and air-conditioned because of the inverter heat. All the data were gathered in the monitoring facility and saved to one server.
In Figure 3, the x-axis is shown in the order of year, month, day, and hour. The y-axis is two, meaning W/m 2 and • C. This is a graph obtained by extracting GNI and ambient temperature data for a day. In this experiment, 60 Cell 265 W multi-crystal modules were used. The multi-crystal module size is 1650 mm × 980 mm × 32 mm (including frame). This PV module type is commonly used in solar energy generation at present.
The temperature sensor was attached to the center part of the cell, assuming that the temperature distribution on the total module area was actually not large. The actual temperature difference between the center and the corner part of the module was shown to be approximately 1-2 • C [2]. The temperature sensor used a T type sensor and checked the temperature. The model name of the T type is TT-T-24-SLE-1000 and the insulation type has a range of −267 • C to 260 • C for the PFA type. Special Limits of Error is ±0.5 • C or 0.4% and it provides high temperature measurement accuracy. Hence, the testbed was used. In order to increase the accuracy, a weather station was also installed at the MegaWatt (MW) verification development center near Yeungnam University and was used to measure irradiation, ambient temperature, wind speed, wind direction, rainfall, humidity etc. Since both sites showed similar results, the temperature and wind speed data were simulated at the test site for PV modules.

Outdoor PV Module Temperature Measurement and Modeling
The equation relevant to the module temperature is based on the following equation: [7]. Input = Output + Generation + Accumulation (4) where Input: irradiation; Output: thermal energy emitted from the PV module to the atmosphere; Generation: energy generated by the PV module through irradiation; Accumulation: input − output − remaining energy after generation.
Input is irradiation; Output is the energy going out to the atmosphere; Generation is the energy that is transmitted after generating electricity; and Accumulation is the remaining energy after the PV module is generated, and it is the remaining energy related to temperature rise. Unlike the previous Equation (1), the decision equation is created. In Equation (4), Output and Generation are Output Energy, and Accumulation is the remaining energy after the module works, so Equation (5) is as follows [7].
Input Energy = Output Energy + Remaining Energy (5) where Input Energy: irradiation; Output Energy: output, going out to the atmosphere + generation; Remaining Energy: energy accumulation, after the module works.
It can be deduced from Equation (5) that the energy increasing the module temperature is the remaining energy, which is converted to heat and which subsequently increases the module temperature [5]. Obviously, the heat cannot be measured precisely because the module is not composed of only one material. However, it is assumed that the module is composed of one new integrated material; therefore, the module temperature can be predicted by multiplying the converted value from the remaining energy. The only energy that could enter from the outdoor condition becomes the irradiation. Therefore, the input energy is the irradiation. The output energy is composed of electric energy generated by the module and thermal energy emitted by the difference between the module and the air temperature. The output energy is based on the convective heat transfer Equation (6) by convective energy heat exchange [8].
The convective heat transfer coefficient and temperature difference are moved to the left side and the remaining right side becomes Equation (6) [8].
Remaining energy is the Input − Output energy balance. It is assumed that the system (target volume) is PV module based on the energy balance equation of chemical engineering thermodynamics [7][8][9]. If Equation (8) is expressed as remaining energy, which is equal to Input − Output and is expressed as Equation (9) The Equation (7) left term is the amount of accumulation, which is a term based on how much the system changes over time.
W is the term for energy access and [(H + 1 2 u 2 + zg ) . m] is the term for material access [10]. That is, the right term has the meaning of Input energy − Output energy as a whole. Since PV module doesn't have material entry, the term for material access is "0". Therefore, W is Input energy − Output energy. It becomes Equation (10). Because the only Input energy entering is insolation. So, Input energy is irradiance. Output energy is P m − loss energy as above. In the loss energy, there is heat energy by convective energy heat exchange.

d(mU) dt
indicates that PV module generates electricity and the remaining energy is converted into heat, which increases temperature of the PV module. Based on this energy balance, the following Equation (11) was inferred [10].
where t: time step; ∆T: temperature rise of module at t + 1 time ( • C); n: module temperature conversion coefficient, with respect to energy (( • C m 2 )/W); I r : irradiation (W/m 2 ); P m : module power generation at t time (W/m 2 ); c 0 : thermal conductivity transferred from module to air (W/(m 2 • C)); T m : module temperature ( • C); T A : air temperature ( • C).
The left term of Equation (11) is the change of the system due to the right term, i.e., the work after the energy entry occurs and the right term can be seen as energy entry in and out from the surface. Equations (10) and (11)  is integrated for t from 0 to 1 (0 is the equilibrium state without energy transfer, 1 is when energy transfer occurs), it becomes mU, where U is the specific internal energy. Since mU = mc∆T (∆T = T 1 − T 0 ), when the area is divided, it becomes mc∆T a∆t [10].
c 0 is based on the convective heat transfer equation. n is based on the equation for the accumulation of the energy balance equation. While n looks similar to specific heat, it is viewed as a unit area rather than J/mol·K per unit mass. PV module internal energy is the place where the remaining total energy exchanged through convective energy will go to power generation by efficiency prior to the increase in temperature. In other words, the temperature increases by providing sensible heat to PV module. Since n and c 0 are complex, they can't be viewed individually. So, their individual experimental values cannot be obtained. As a result, the following Equation (15) is derived.
The interval of time t varies depending on the measurement time. In this study, the time interval was maintained at 30 s. The initial module temperature was assumed equivalent to the air temperature, as there was no generation without any irradiation. Starting from the following irradiation, the module power generation energy and the energy transferred to the air were calculated, with respect to actual irradiation, in order to predict the module temperature. The open-circuit voltage (Voc), voltage at maximum-power point (Vmp), short-circuit current (Isc), and current at the maximum-power point (Imp) equation presented by the Sandia National Laboratory was used to calculate the module generation with respect to irradiation [3]. 265 W PV module Isc = 9.0922 A; Voc = 38.1612 V; Imp = 8.5144 A; Vmp = 30.2515 V; and comparison target 260 W PV module Isc = 8.96 A; Voc = 39.213 V; Imp = 8.5 A; Vmp = 31.361 V. The elements n and c 0 in Equation (15) can be characterized from the collected data. As shown in Figure 4, when the module starts to generate electricity, the PV module temperature rises to generate electricity. Assuming 1000 W/m 2 of solar radiation is generated, electricity is converted to 15%-20%. The remaining energy is energy related to heat, which is consumed outside the module rather than PV module inside like convective energy. However, PV module temperature rises due to the presence of remaining thermal energy generated and accumulated. There is a large difference between ambient temperature and PV module temperature. PV module temperature is inevitably higher than the ambient temperature. The higher the PV temperature, the lower the generation due to the characteristics of the PV module temperature change. In addition, the reason that voltage is affected by temperature and current is affected by insolation is also due to characteristics of the PV module temperature change. The conversion coefficient n and c 0 is calculated as follows. First, a region, such as part B in Figure 4, where the module temperature does not rise or fall, is identified with clear irradiation and low wind speed (2 m/s or below). Temperature that does not change (rise or fall) means that ∆T = 0. Therefore, the following equation becomes valid: After identifying c 0 here, the n value can be determined by using the data in part A of Figure 4. c 0 can be expanded as below to include the wind effect [3][4][5].
The initial c 0 becomes c 1 , as it is the conversion coefficient identified when the wind speed was low. After identifying c 1 and n, c 2 can be identified, based on the data of the date with wind speed that is higher than a certain value. In addition, it is necessary to verify the value of n, c 0 by looking at the entire point through the graph in Figure 4. In the case of 200~900 W/m 2 or more based on insolation, the average values n and c 0 are obtained by dividing a total of 8 sections and after substituting into Equation (15) and simulating. However, due to KS C IEC/TS 61724-2:2016, the minimum POA insolation intensity during the four seasons is more than 450 W/m 2 . Then, 6 equations with different n, c 0 are expressed. Among the regressing analysis methods, if the RMS (Root Mean Square) error is used and a combination of coefficients when the RMS error of the module temperature values is minimized is used, the result of conversion coefficient can be more accurately expressed [11,12].

Results
As a result of applying Equation (15) to the measured 265 W PV module in this study, the extracted coefficients are n = 2.33 ( • C m 2 )/W, c 1 = 33.26 W/(m 2 • C), and c 2 = 5.31 W/(m 3 • C)). The measured and calculated daily change in module temperature for one day is presented in Figure 5. The mean difference was 1.64 • C and the maximum difference was 5.63 • C. In addition, the difference between the actual energy generation and the predicted generation was approximately 5.2%, as shown in Figure 6. The results that were better than the mean error of 4.23 • C, calculated with the Equation (1), are shown. Figure 7 shows the results of the module temperature predicted by using the Sandia National Laboratory equations [12,13].  Instances of substantial error are found in previous research results. The existing mass-produced Si-type modules have similar values. Figure 8 indicates the results of a comparison between the predicted values and the measurement values of another 260 W mass-produced module.
The mean error of the measurements and predictions is approximately 7.53%. This shows that the Si-type modules have a similar trend.

Conclusions
With regard to the temperature rise of the PV module, the prediction of the module temperature showed slight differences in a shorter time than when thermal equilibrium occurs and as compared with the existing research results. In other words, the temperature or generation is usually predicted at intervals of more than one hour, but this experiment can be predicted in real-time in seconds or minutes. Therefore, the shorter the interval, the smaller is the error. For this, a new prediction equation was suggested and was implemented in this experiment. The results indicate that prediction was possible within an error range of approximately 5% in comparison with the actual generation. Obviously, not all the modules have the same conversion coefficient; however, the Si-type modules that use similar materials show similar results.
The daily irradiation, wind direction, and wind speed data could vary; however, the annual average does not show substantial changes. If there were an annual or longer period data of irradiation, air temperature, and wind speed, the annual generation in the actual outdoor condition could be predicted by weeklong tests or 10-day tests on the new module. In addition, the results of the current study can be comparatively analyzed with several modules. Such comparison is needed to develop modules with improved performance.
In addition, studies on the influence of humidity on the module temperature should be performed. Clearly, the temperature changes in the module would be different in dry desert regions and in humid tropical regions. The n, c 1 , and c 2 values that were calculated in this study could vary depending on regions, and new variables could be added. Therefore, to promote research in this wide field, the environmental data of various climates should be accumulated and shared. This experiment was supported by the equipment and technology of Yeungnam University Daekyung solar cell regional innovation center.

Conflicts of Interest:
The authors declare no conflict of interest.