Photovoltaic Panel Parameters Estimation Using Grey Wolf Optimization Technique "2279

In different photovoltaic PV applications, it is very important to model the PV cell. However, the model parameters are usually unavailable in the datasheet provided by the manufacturers and they change due to degradation. This paper presents a method for identifying the optimal parameters of a PV cell. This method is based on the one diode model using the grey wolf algorithm as well as datasheets. An algorithm is implemented in a SIMULINK simulator for making the I-V and P-V characteristics. This approach is found to be useful for designers due to its simplicity, fastness, and accuracy. The final results are compared to demonstrate the efficiency and accuracy of the proposed method.


Introduction
Among all renewable energy sources, solar energy has acquired the highest growth rate worldwide in recent years. The major application of solar energy is photovoltaic (PV) power generation, which saw growth of more than 22% in 2019 and is predicted to output over 720 TWh [1]. PV systems are easy to install, noise-free, and can directly convert solar energy to electrical energy. For an accurate study in different PV applications, it is very important to model the basic device of the PV cell. However, the model parameters are usually unavailable in the datasheet provided by the manufacturers and their values change over time due to the PV degradation [2]. Thus, how to estimate appropriate parameters is of high importance and has attracted immense interest among researchers. The one diode model (ODM) is considered the most suitable model used to characterize the solar cells/modules [3][4][5][6][7] in comparison to the double diode model (DDM) and the three-diode model (TDM) as it has a minimum number of parameters and a good level of accuracy. The five electrical parameters of the ODM are: photocurrent (I ph ), diode ideality factor (n), reverse saturation current (I 0 ), shunt resistance (R sh ), and series resistance (R S ).
Several methods have been developed to extract the ODM parameters which are classified into three main categories [5,7]: -Analytical methods [7] -Numerical methods [3-6] -Artificial intelligence or optimization methods [8][9][10] Numerical methods are widely used in the literature since they provide a good compromise between speed of calculation, simplicity, and accuracy. These numerical methods are utilized to solve a system of a few non-linear equations related with the PV cell.
The objective of this work is to simulate the I-V and P-V characteristics using a one diode model (ODM) associated with a developed algorithm that permits finding the appropriate value of the ideality factor and hence extracting the other needed parameters simply from datasheet information provided by the manufacturer.
The objective of this work is to simulate the I-V and P-V characteristics diode model (ODM) associated with a developed algorithm that permits find propriate value of the ideality factor and hence extracting the other needed simply from datasheet information provided by the manufacturer.

Photovoltaic Cell
A PV cell is basically a semiconductor diode, as shown in Figure 1. Acco principle theory of photoelectric effect, for any material exposed to light genera carriers, when sunlight that is basically photons of different frequencies and e the solar cell surface and is absorbed by a semiconducting material, such as s trons will be excited from atoms, making them free and ready to move. Solar cells produce electricity with very small voltage. For the common tion silicon solar cell, it is approximately 0.5-0.6 volts. They are available in modules or panels to provide sufficient voltage and current for real life applic

Characteristics of the PV Cell
Electrical characteristics of PV modules are given by the producers un conditions that are known as standard test conditions (STC). Such conditions by the ambient temperature T = 25 °C, irradiation level G = 1000 W/m mass value AM = 1.5. However, in the working field, PV modules operate at peratures and somewhat lower insulation conditions. In order to determine output of the solar cell, it is important to know the expected operating tempe nominal operating cell temperature (NOCT) is defined as the temperature open circuited cells in a module under the conditions: solar irradiance G = 80 temperature T = 20 °C, and wind speed = 1 m/s. Then, the cell tempera calculated by the following equation: where G is taken in (w/m ). The typical I-V and P-V characteristics of a photovoltaic cell are shown The main three significant parameters on the photovoltaic characteristics are voltage (V ), short circuit current (I ), and maximum power point at (V , I Solar cells produce electricity with very small voltage. For the common single junction silicon solar cell, it is approximately 0.5-0.6 volts. They are available in the form of modules or panels to provide sufficient voltage and current for real life applications.

Characteristics of the PV Cell
Electrical characteristics of PV modules are given by the producers under precise conditions that are known as standard test conditions (STC). Such conditions are defined by the ambient temperature T STC = 25 • C, irradiation level G STC = 1000 W/m 2 , and the air mass value AM = 1.5. However, in the working field, PV modules operate at higher temperatures and somewhat lower insulation conditions. In order to determine the power output of the solar cell, it is important to know the expected operating temperature. The nominal operating cell temperature (NOCT) is defined as the temperature reached by open circuited cells in a module under the conditions: solar irradiance G = 800 w/m 2 , air temperature T ambiant = 20 • C, and wind speed = 1 m/s. Then, the cell temperature can be calculated by the following equation: where G is taken in (w/m 2 ). The typical I-V and P-V characteristics of a photovoltaic cell are shown in Figure 2. The main three significant parameters on the photovoltaic characteristics are open circuit voltage (V oc ), short circuit current (I sc ), and maximum power point at (V mpp , I mpp ).
simply from datasheet information provided by the manufacturer.

Photovoltaic Cell
A PV cell is basically a semiconductor diode, as shown in Figure 1. Accordin principle theory of photoelectric effect, for any material exposed to light generating carriers, when sunlight that is basically photons of different frequencies and energ the solar cell surface and is absorbed by a semiconducting material, such as silico trons will be excited from atoms, making them free and ready to move. The maximum current in photovoltaic cell is produced when there is a short circuit between its positive and negative terminals and it is denoted as I sc . As I = I sc , the voltage in the circuit is zero. The maximum voltage, denoted as V oc , occurs when there is a break in the circuit.
The maximum power achieved from a photovoltaic cell occurs at a point on the bend in the I-V curve known as the maximum power point (MPP), which can be found at voltage and current points designated as V mpp and I mpp .
Generally, these parameters are provided in the datasheet by manufacturers of a particular photovoltaic cell or module. When the PV cell is connected to an external load, the electrical characteristics of the load determine the actual point on the I-V curve at which the photovoltaic cell operates.

Single Diode Model of PV Cell
To analyze characteristics of solar cells, electrical equivalent circuits are needed and hence modeled using simulation software. Researchers have developed mathematical models to understand and predict the effect of changing conditions on photovoltaic electrical output. The lumped parameter model is one of these models classified based on the number of diodes. It is widely used and has proven to be more successful. The lumped parameter models can take the form of a one diode model, double diode model, or three diode model. Although the accuracy of the characteristics of the model improves as the number of diode increases, the required mathematical expression to obtain the output characteristics become more complex. For simplicity, in this work, the one diode model known as the five parameter model is chosen for the identification of photovoltaic cell parameters.
The one-diode model takes into account different properties of solar cell: R s is introduced as to consider the voltage drops and internal losses due to flow of current, and R sh takes into account the leakage current to the ground when the diode is reverse biased. However, this model has neglected the recombination effect of the diode. Therefore, it is still not the most accurate model.

Identification of PV Cell Parameters
The current versus voltage relationship of the single diode PV cell model ( Figure 3) is presented in Equation (2). It can be noticed that the equation of the I-V curve is nonlinear, which is difficult to solve by the analytical methods. Due to this difficulty, scientists have developed several algorithms in order to solve this equation to determine the parameters of the solar cell. In this work, metaheuristic methods, which can be adapted to solve a wide range of optimization problems, are used. These methods are designed to find a good solution among a large set of feasible solutions with less computational effort than other optimization techniques.

Grey Wolf Optimization Technique
Grey wolf optimization (GWO) is a population-based met method inspired by grey wolves (Canis lupus). The GWO algor ship hierarchy and hunting mechanism of grey wolves in natur into four groups, alpha, beta, delta, and omega, for simulating t Alpha is the first level and is the leader of the pack. Beta is the se chy of wolves, as they help alpha wolves to make a decision. D level in the pack, as members have to succumb to alpha and beta, omega. Omega wolves have the lowest position in the pack, havin dominant wolves. Figure 4 shows the grey wolf social hierarchy In addition to the social hierarchy of wolves, group huntin social behavior of grey wolves. According to Muro et al. [13], the m hunting are as follows: -Tracking, chasing, and approaching the prey.

Grey Wolf Optimization Technique
Grey wolf optimization (GWO) is a population-based meta-heuristic optimization method inspired by grey wolves (Canis lupus). The GWO algorithm mimics the leadership hierarchy and hunting mechanism of grey wolves in nature. The pack is classified into four groups, alpha, beta, delta, and omega, for simulating the leadership hierarchy. Alpha is the first level and is the leader of the pack. Beta is the second level on the hierarchy of wolves, as they help alpha wolves to make a decision. Delta represents the third level in the pack, as members have to succumb to alpha and beta, however they dominate omega. Omega wolves have the lowest position in the pack, having to succumb to all other dominant wolves. Figure 4 shows the grey wolf social hierarchy [11,12].

Grey Wolf Optimization Technique
Grey wolf optimization (GWO) is a population-based met method inspired by grey wolves (Canis lupus). The GWO algor ship hierarchy and hunting mechanism of grey wolves in natur into four groups, alpha, beta, delta, and omega, for simulating t Alpha is the first level and is the leader of the pack. Beta is the se chy of wolves, as they help alpha wolves to make a decision. D level in the pack, as members have to succumb to alpha and beta, omega. Omega wolves have the lowest position in the pack, havin dominant wolves. Figure 4 shows the grey wolf social hierarchy In addition to the social hierarchy of wolves, group huntin social behavior of grey wolves. According to Muro et al. [13], the m hunting are as follows: -Tracking, chasing, and approaching the prey. -Following, encircling, and harassing the prey until it stops m -Attacking the target.
Wolves encircle the target during the chase. The encircling tions (7) and (8)   In addition to the social hierarchy of wolves, group hunting is another interesting social behavior of grey wolves. According to Muro et al. [13], the main phases of grey wolf hunting are as follows: -Tracking, chasing, and approaching the prey. -Following, encircling, and harassing the prey until it stops moving. -Attacking the target.
Wolves encircle the target during the chase. The encircling is modeled by the Equations (7) and (8) [13]: A and → C , as computed through these two equations: Components of → a linearly reduce from 2 to 0 through the iterations and → r 1 and → r 2 as random vectors in [0, 1].
The pursuit is habitually directed by the leader alpha (α) followed by beta (β) and delta (δ) which can sometimes contribute in chasing. (δ) and (ω) look after the injured wolves in the group. Alpha (α) is considered the best result owing to the best information of the place of the target, while beta (β) and delta (δ) are the second and the third best solutions respectively in designing GWO. (ω) is the last best. Therefore, the first three best solutions obtained so far are saved and the other search agents (including ω) are obliged to update their positions according to the position of the best search agent. When the prey stops moving, the wolves terminates the chase by attacking it, as shown in Figure 5. D δ ⃗ = |C ⃗ · X δ ⃗ (t) − X ⃗ (t)| X ⃗ = X δ ⃗ (t) − A ⃗ · D δ ⃗ X ⃗ (t + 1)= X ⃗ + X ⃗ + X ⃗ 3

ODM Parameters Extraction Using a GWO Method
The proposed five-parameter estimation method is based on t characterize the I-V curve, which are the maximum power point (V point (0, I ), and open circuit point (V , 0 ). These data can be provi or extracted from the experiments.
Like any optimization algorithm, an objective function F(X) sho is based on the one diode Equation (2). For the identification of PV the number of series connected cells N is considered. Thus, the obje  In this regard, the following formulas are applied: Eng. Proc. 2022, 14, 3 6 of 10

ODM Parameters Extraction Using a GWO Method
The proposed five-parameter estimation method is based on the three points that characterize the I-V curve, which are the maximum power point (V mpp , I mpp ), short circuit point (0, I sc ), and open circuit point (V oc , 0 ). These data can be provided in the datasheet or extracted from the experiments.
Like any optimization algorithm, an objective function F(X) should first be set. F(X) is based on the one diode Equation (2). For the identification of PV module parameters, the number of series connected cells N s is considered. Thus, the objective function is: The fitness function that needs to be minimized in order to quantify the error is the root main square error (RMSE) between the approximated values in the datasheet I mpp , I sc and I = 0 at open circuit condition and the calculated ones.
Then, the NRMSE error is calculated as follows: The pseudo code of the proposed algorithm is presented in Algorithm 1.
Algorithm 1 GWO 1: Input: T, N s , V mpp , I mpp , V oc , I sc 2: Output: X α 3: Initialize the grey wolf population Xi (i = 1, 2, ... n) 4: Initialize a, A, and C 5: Calculate the fitness of each search agent by Equation (12) 6: X α = the best search agent 7: X β = the second best search agent 8: X δ = the third best search agent 9: while (t < Max number of iterations) 10: for each search agent 11: X δ = the third best search agent 12: while (t < Max number of iterations) 13: for each search agent 14: Update the position of current search agent by Equation (10) 15: end for 16: Update a, A, and C 17: Calculate the fitness of all search agents by Equation (12) 18: Update X α , X β and X δ 19: t = t + l 20: end while 21: return X α 22: Calculate NRMSE using Equation (13) 23: end procedure

Test Results and Discussion
The proposed algorithm GWO is used to extract the parameters of the ODM based on curve fitting method. The algorithm is tested on two PV modules and compared with other Eng. Proc. 2022, 14, 3 7 of 10 algorithms to prove its effectiveness. The simulator is developed using MATLAB R2016a and executed on a PC with Intel ® Core™ i5-2450M CPU processor @ 2.50 GHz, 4 GB RAM, under Windows 10 64-bit OS.
The search ranges used in the optimization of the five parameters are given Table 1. However, the GWO parameters are given in Table 2. The datasheet parameters of used PV modules are presented in Table 3.   Case Study #1: KC200GT The polycrystalline module KC200GT ODM parameters are extracted to draw the I-V and P-V characteristics using the developed PV simulator, as shown in Figure 6. The results obtained at STC are presented in Table 4. The GWO results for the KC200GT module are compared with other published optimization methods results to prove its efficiency.   As given in Table 4, the normalized mean absolute error NRMSE obtained by the proposed method (0.02072%) is the lowest, which proves the effectiveness of this technique in extracting the unknown PV parameters. Moreover, the convergence using this method is very fast, where the simulation execution time is less than 7 s. Case Study #2: STP050D-12/MEA The PV characteristics of the poly-crystalline STP050D-12/MEA module are simulated at T = 45.57 • C and G = 632 w/m 2 using the GWO based ODM PV simulator as shown in Figure 6. The obtained results are given in Table 5 and Figure 7. The computed errors are: RMSE = 1.327 × 10 −3 and NRMSE = 7.282 × 10 −2 . The simulated characteristics are compared with the module's experimental data obtained in our research laboratory [14]. The comparison results are presented in Figure 8 and Table 6.      It can be noticed that the measured data points obtained in our laboratory for the STP050D-12/MEA module highly agreed with the simulated curves. Thus, the obtained I-V and P-V characteristics for the STP050D-12/MEA module based only on the datasheet information using the GWO based ODM PV simulator gives good results.

Conclusions
PV cell characterization is a hot research topic in the field of renewable energy. Obtaining the most accurate I-V and P-V characteristics has been the main purpose of this research work. This was attempted by finding an efficient method based on the one diode model using the information provided by the manufacturers. However, not all (five) parameters are available in the datasheet. Thus, the simulation has been associated with the developed algorithm that permits finding the appropriate value of the needed parameters. Grey wolf optimization (GWO) has been implemented using SIMULINK for extracting the five parameters based only on the three critical I-V points information provided in the datasheet, namely the open circuit point (0, V oc ), the short circuit point ( I sc , 0), and the maximum power point ( I mpp , V mpp ). Furthermore, this method has been enhanced to provide the parameters under different environmental conditions. The developed ODM-GWO simulation was tested for various PV modules, under different temperatures and irradiances. The obtained I-V and P-V curves were compared to the characteristics provided in the datasheet. Moreover, the BP MSX 120 multi-crystalline module simulated curves were compared to the experimental I-V data. The measured data points conformed to the obtained curves. Thus, the good accuracy of the developed PV simulator was demonstrated.
Adding to the fact that these proposed algorithms have provided optimal results with an acceptable accuracy, the time taken by the ODM-GWO simulation execution is less than 10 s.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the data may involve confidential information of our research group.

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