Next Article in Journal
Joint Successful Transmission Probability, Delay, and Energy Efficiency Caching Optimization in Fog Radio Access Network
Previous Article in Journal
Fast Flux and Torque Control of a Double Inverter-Fed Wound Machine Considering All Coupling Interferences
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recent Meta-Heuristic Algorithms with a Novel Premature Covergence Method for Determining the Parameters of PV Cells and Modules

1
Department of Computer Science, Faculty of Computers and Informatics, Zagazig University, Zagazig 44519, Egypt
2
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
3
Department of Computational Mathematics, Science, and Engineering (CMSE), College of Engineering, Michigan State University, East Lansing, MI 48824, USA
4
Department of Computer Science and Engineering, Soonchunhyang University, Asan 31538, Korea
5
Electrical Power and Machines Department, Faculty of Engineering, Zagazig University, Zagazig 44519, Egypt
*
Authors to whom correspondence should be addressed.
Submission received: 22 June 2021 / Revised: 14 July 2021 / Accepted: 28 July 2021 / Published: 31 July 2021
(This article belongs to the Section Power Electronics)

Abstract

:
Currently, the incorporation of solar panels in many applications is a booming trend, which necessitates accurate simulations and analysis of their performance under different operating conditions for further decision making. In this paper, various optimization algorithms are addressed comprehensively through a comparative study and further discussions for extracting the unknown parameters. Efficient use of the iterations within the optimization process may help meta-heuristic algorithms in accelerating convergence plus attaining better accuracy for the final outcome. In this paper, a method, namely, the premature convergence method (PCM), is proposed to boost the convergence of meta-heuristic algorithms with significant improvement in their accuracies. PCM is based on updating the current position around the best-so-far solution with two-step sizes: the first is based on the distance between two individuals selected randomly from the population to encourage the exploration capability, and the second is based on the distance between the current position and the best-so-far solution to promote exploitation. In addition, PCM uses a weight variable, known also as a controlling factor, as a trade-off between the two-step sizes. The proposed method is integrated with three well-known meta-heuristic algorithms to observe its efficacy for estimating efficiently and effectively the unknown parameters of the single diode model (SDM). In addition, an RTC France Si solar cell, and three PV modules, namely, Photowatt-PWP201, Ultra 85-P, and STM6-40/36, are investigated with the improved algorithms and selected standard approaches to compare their performances in estimating the unknown parameters for those different types of PV cells and modules. The experimental results point out the efficacy of the PCM in accelerating the convergence speed with improved final outcomes.

1. Introduction

Solar energy converted to electric power using a photovoltaic (PV) system offers considerable opportunities to overcome the drawbacks of the traditional energy sources in terms of unavailability, environmental pollution, and global warming [1,2,3,4,5]. Some of theadvantages which show the importance of solar cells as mentioned in [6] are stated as follows:
  • Less operational cost [7];
  • Low maintenance [8];
  • Reducing air population [9].
Solar energy also varies daily and is affected by transient obstructions due to cloud cover; thus, there is an essential need to optimize the performance of the PV system to reach the best output, especially under irradiance and temperature variation. Simulation, evaluation, and control of PV systems require an accurate model that accurately represents the nonlinear current-voltage (I-V) characteristics curve of the PV cells. Consequently, over the last decades, several different models have been proposed for simulating this curve [10].
Based on the PV model, analytical, deterministic, and meta-heuristic methods have been proposed to estimate the unidentified parameters of the three principal types of PV models, namely, single diode model (SDM), double diode model (DDM), and triple diode model (TDM), to improve the performance of the PV system. Analytical methods are based on the solution of a series of mathematical equations and have the advantages of ease of implementation and speed. However, there is not a sufficient match between the simulated and measured I-V curves [11,12]. Deterministic methods have constraints on the model, such as convexity and differentiability, depend heavily on the initial guess, and are easily trapped in local minima [13]. PV models are often implicit, multimodal, and nonlinear, which means that they are poorly solved using deterministic methods, such as Lambert W-functions [14], iterative approach [15], and the Newton–Raphson method [12,16].
Meta-heuristic algorithms have been shown to overcome several real-world problems in a reasonable time with high accuracy [17,18,19,20,21,22] and therefore have been widely applied for tackling the parameter identification problem (PIP) of PV models. This section reviews the major applications and their shortfalls.
Xiong et al. [23] proposed the competitive swarm optimization (CSO) approach based on an advanced variant of particle swarm optimization (PSO) for identifying the unidentified parameters of the PV models. However, the authors found that this advanced variant of PSO still suffers from falling into local minima for complex multimodal optimization problems such as the PIP of PV models because of its weak exploration capability. Therefore, CSO was improved by two strategies to improve its exploration capability. The first was the winner-leading search strategy proposed to make the losers explores more regions within the search space. The second strategy is the Gaussian mutation operator that was proposed to improve the exploration operator of CSO to escape the local minima problem. The resultant version of CSO integrated with the winner-leading search and the Gaussian mutation was called WLCSODGM and was extensively investigated on four PV models and compared with 12 optimization algorithms to check its superiority.
Diab et al. [24] proposed the coyote optimization algorithms (COA) for tackling the PIP of SDM, DDM, and TDM and was observed on multi-crystalline, mono-crystalline, and thin-film PV modules under various irradiance and temperature levels. Long et al. [10] integrated both the grey wolf optimizer and cuckoo search (CS) algorithm (GWOCS) for identifying the parameters of DDM, SDM, and PV modules based on the test points measured under various operation conditions. For increasing the diversity among the members of the GWO, the opposition learning strategy was integrated to balance the exploration and exploitation operators of GWOCS while avoiding local minima and moving accurately and quickly toward the optimal solution. To validate the performance of GWOCS, ten complex mathematical functions were investigated, in addition to estimating the parameters of the different PV models as a harder problem with several local minima.
Ridha et al. [25] proposed the boosted Harris hawks optimization algorithm (BHHO) for estimating the parameters of the solar cell SDM by taking into consideration the sensitivity under various sunlight and temperature conditions. BHHO was enhanced by renting the exploration operator of the flower pollination algorithm and the vigorous mutation scheme of the DE to move quickly the individuals in a direction of the optimal solution and explore extensively the search space of the problem to find the most promising regions. The experimental outcomes show that BHHO is superior to some of the well-known meta-heuristic optimization algorithms. Ultimately, Table 1 provides a brief review of some of the optimization algorithms proposed in the last two years for the SDM, DDM, and TDM PV models.
From the preceding review, it is clear that most of the algorithms suffer from low convergence speed, consuming a large number of iterations without any great benefit. Therefore, in this paper, a new method called the premature convergence method (PCM) has been proposed to help the optimization algorithms in utilizing the iterations as much as possible to accelerate convergence and achieve better accuracy. PCM is based on a controlling factor used to determine if the current particle will be updated around the best-so-far solution in the direction of one of the following:
(a)
Two solutions selected randomly from the population;
(b)
The current position and the best-so-far solution;
(c)
Balancing between the previous two steps.
PCM is integrated with three well-known meta-heuristic algorithms: Harris hawks optimization [38], moth–flame optimizer [39], and equilibrium optimizer (EO) [40] to identify its efficacy for estimating the parameter of the SDM. The three algorithms were here used because they have a high exploration rate at the beginning of the optimization process that may lead to the consumption of a large number of iterations without reaching better outcomes to accelerate convergence in a direction of the optimal solution. Moreover, the search direction of these algorithms is not based on the direction in which the current population move, and that may make the algorithms search randomly for better outcomes in the exploration case or move in the same direction of the best-so-far solution that may promote falling into local minima in the exploitation case.
A PV cell based on the RTC France and three PV modules (Photowatt-PWP201, STM6-40/36, and Ulta 85-P) were used to investigate the performance of the PCM-improved versions of each algorithm compared with the standard version. The experimental results show the effectiveness of the integration of PCM with each of the three optimization algorithms, particularly with EO.
The remainder of this paper is organized as follows: Section 2 describes the SDM and PV module model. Furthermore, Section 3 describes the meta-heuristic algorithms and the premature convergence method. Section 4 gives the experimental results of the studied algorithms in addition to some discussions on those results. Finally, Section 5 provides conclusions and some discussion on future work.

2. Mathematical Descriptions of the Problem

This section described in detail the mathematical model of the single diode model and the PV module model, in addition to the objective function used in this work.

2.1. Single Diode Model (SDM)

As illustrated in Figure 1, the SDM is simple, comprising the photo-generated current I p h [41], the diode current as estimated using Equation (1), the shunt resistor current as calculated according to Equation (3), and I that indicates the output current of the SDM.
I D = I s d ( e x p ( V + I × R s n × V t ) 1 )
where I s d is the diode current, V indicates the output voltage, R s is the series resistance, n is the ideality factor of the diode, and V t is the junction thermal voltage and mathematically formulated as
V t = k × T q
T is the temperature of the junction in kelvin, k is the Boltzmann constant (1.3806503 × 10−23 J/K), and q refers to the electron charge (1.60217646 × 10−19 C).
I s h = V + I × R s R s h
Rsh is the shunt resistance. By replacing ID and Ish in Equation (4) with their formula defined in Equations (1) and (3), I is as follows:
I = I p h I D I s h
I = I p h I s d ( e x p ( q ( V + I × R s ) n × k × T ) 1 ) V + I × R s R s h
For the SDM, five known parameters (Iph, Isd, n, Rs, Rsh) therefore need to be estimated efficiently to optimize the performance of the PV system.

2.2. Photovoltaic (PV) Module Model

The PV cell produces extremely small amounts of useful electricity on its own, and therefore, Ns PV cells are connected in series to form a practical solar cell unit known as a PV module to supply the required output current and power. The PV module can be formulated as Equation (5) with the difference that V t is computed as follows: V t = (Ns kT)/q [42].
For the nonstandard conditions unlike STC, the above-stated mathematical model should be altered to show the performances under varied cell temperatures and changed radiation levels. The equation’s model should be adapted to accommodate such changes due to G and T variations as follows [6,35,36,37,43,44]:
I p h = G G S T C [ I p h S T C + K i ( T 25 ) ]
V o c = V o c S T C + K v ( T 25 )
E g = E g S T C [ 1 2.677 × 10 4 ( T 25 ) ]
I s d i = I s d i S T C ( T 25 ) 3 · e x p ( q · E g ( T 25 ) 25 a i · K · T ) ,     i 1 : 3
R s h = R s h S T C · ( G S T C G )
where K i and K v define the thermal coefficient of current and voltage, respectively, E g denotes the semiconductor band-energy, V o c denotes open-circuited voltage, and others are the normalized parameters at STC.

3. Meta-Heuristic Algorithms and the Premature Convergence Method

The meta-heuristic algorithms have sought to find the near-optimal solutions of the optimization problem based on two operators—exploration and exploitation. In the exploration operator, the algorithm explores most regions of the search space searching extensively for the promising region that may significantly contain the optimal solution. Then, within the first half of the optimization process, the exploitation operator may be applied to search around the current promising region for better solutions. In the second half of the optimization process, the exploitation is mandatorily applied as an attempt to exploit the last current promising region that may contain the global solution.
In this section, three optimization algorithms are described, in addition to the PCM. First, before describing each algorithm separately, it is useful to explain the stage shared between them—the initialization step. In the initialization step, a population with N individuals of d dimensions, where d is the number of the unknown parameters in the problem, are distributed within the search space of the problem using Equation (11) as the start point upon which the optimization process relies.
X i = L m i n + ( U m a x L m i n ) × r   ,     i N
where X i is a vector to load the position of the ith individual, r is a vector generated randomly at the range of 0 and 1, and L m i n and U m a x are two vectors, including the lower bound and upper bound of each unknown parameter in the solved problem.

3.1. Objective Function

The root mean squared error (RMSE) between the measured and estimated current that is computed under the estimated parameters and the Newton–Raphson method [43] will be used as an objective function to evaluate the quality of each solution in the population. This objective function defined according to RMSE is formulated as
RMSE = f ( X i ) = 1 M × k = 1 M ( I m I e ) 2  
where I m indicates the measured current, M is the number of the measured test points, and I e is the simulated current and defined as
I k + 1 = I k F ( I ) F ( I )
F(I) is computed according to Equation (5) and F’(I) is the first derivative of the same equation with respect to I. In common practice, five iterations are sufficient to solve the above, as stated in (13) with a very acceptable tolerance using the NR method.

3.2. Premature Convergence Method (PCM)

In this section, a novel method known as a premature convergence method is proposed to help the optimization algorithm in accelerating the convergence speed with avoiding becoming trapped in local minima. This method also includes a control factor r generated randomly between 0 and 1. According to the control factor, the exploitation capability will be significantly encouraged when r > 0.5, while the exploration operator is applied to the current position if r < 0.5, and balancing between the exploration and exploitation operators is achieved in a case of r = 0.5. Finally, each solution under the PCM will be updated according to the following equation:
n X i ( i t + 1 ) = X * + ( 1 r ) ( X a ( i t ) X b ( t ) ) + ( r ) ( X * X i ( i t ) )
where a and b are two integers selected randomly between 1, and N and represent the indices of two individuals within the population. n X i is a vector used to store the next position of the ith individual, and X * is the best-so-far solution obtained by any optimization algorithm. Finally, the steps of this method are presented in Algorithm 1.
Algorithm 1 The steps of PCM
Input: current population X i ,   i = 1 , 2 , 3 , 4 , , N
1.   for each i individual
2.       Update X i using Equation (14)
3.       If f ( X i ) > f ( n X i )
4.            X i = n X i ;
5.       end
6.   end for
Output: return X
After describing the main steps of PCM, in the next subsections, this method is integrated with three well-known optimization algorithms to identify its influence on the performance of those algorithms.

3.3. Meta-Heuristic Algorithms

3.3.1. Equilibrium Optimizer

Recently, a novel physics-based optimization algorithm known as the equilibrium optimizer (EO) was proposed by Faramarzi [40] for solving global optimization problems. The mathematical model of EO is described in detail.
At the outset, the equilibrium state of the system is not known; thus, EO considers the best-so-far four particles, in addition to their average as the equilibrium candidates, and adds them in the ascending order according to their fitness values if the problem is minimized inside an equilibrium pool p e q , p o o l , as described in Equation (15); the first four particles encourage the exploration capability, while the last promotes the exploitation capability.
p e q , p o o l = [ X e q   ( 1 ) ,         X e q   ( 2 ) ,       X e q ( 3 ) ,       X e q ( 4 ) ,         X e q ( a v g ) ]
After defining the candidate solutions, the optimization process begins to update the individuals within the population by searching for other better solutions. EO seeks a reasonable balance between exploration and exploitation according to the following equation:
F = a 1 s i g n ( r 0.5 ) [ e λ ( t ) 1 ]
λ and r are two vectors containing values created randomly between 0 and 1, and t is computed as
t = ( 1 i t t m a x ) ( a 2 ( i t t m a x ) )
where it indicates the current generation, t m a x indicates the maximum of generations, a2 is a constant value to control the exploitation, and a1 is a constant value controlling the diversification (exploration) capability. Another factor G represents the generation rate and is used to enhance the intensification/exploitation operator of EO,
G = G 0 F
G 0 = G C P ( X e q λ * X )
G C P = { 0.5 r 1                     r 2 > G P   0                           o t h e r w i s e
where r1 and r2 are numbers generated randomly within 0 and 1. GP is a constant parameter used to balance the ratio between diversification and intensification operators of EO. In the end, each individual within the population is generally updated as defined in Equation (21).
X = X e q + ( X X e q ) F + R λ V ( 1 F ) ,   V = 1
Finally, the steps of EO integrated with the PCM are given in Algorithm 2 and Figure 2. Integrating PCM with EO helps to explore the promising regions obtained by it at each generation as an attempt to see if the near-optimal solution is there or not. If the optimal solution is found, that will reduce the number of function evaluations consumed by the standard algorithm even reaching or might the maximum evaluations is terminated without reaching this optimal solution. Therefore, we integrated this strategy to improve the exploitation capability of the meta-heuristic algorithms, which has a high-exploration operator at the beginning of the optimization process to avoid the initial time-consuming process by the standard algorithm in its search for the promising region; this solution already emerged in the previous generation and need only to be more focused for achieving the near-optimal solution.
Algorithm 2 The pseudocode for the PCM integrated EO algorithm (PEO)
1. Initialize a population of N individuals using Equation (6)
2.  Set a 1 = 2 ; a 2 = 1 ; G P = 0.5
3.  Initialize the equilibrium pool p e q , p o o l with a large value for fitness.
4.  it = 1;
5.  while (it < t max )
6.      for each i individual
7.      Compute the fitness value f ( x i )
8.      Updating p e q , p o o l , if x i better
9.      end for
10.      X e q ( a v g ) = (   X eq   ( 1 ) + X eq   ( 2 ) + X eq   ( 3 ) + X eq   ( 4 ) ) / 4
11.     Update
          p e q , p o o l = [   X e q   ( 1 ) , X e q   ( 2 ) , X e q (   3 ) , X e q (   4 ) , X e q ( a v g ) ]
12.     Accomplish the memory saving
13.     Compute t using Equation (17)
14.     for each i individual
15.        Select a candidate from p e q , p o o l randomly.
16.        Assign random values to two vectors, namely r , λ
17.        Compute F   based on Equation (16)
18.        Compute G C P   based on Equation (20)
19.        Calculate G 0   based on Equation (19)
20.        Calculate G based on Equation (18)
21.        Update X i using Equation (21).
22.        it++
23.     end for
24.     Update the current population using Algorithm 1.
25.     After then, applying the steps from Line 12–17 with skipping
          Line 13 in this algorithm
26.  it = it + N;
27. end while

3.3.2. Moth–Flame Optimizer

Mirjalili [39] proposed a novel meta-heuristic algorithm, called the moth–flame optimization algorithm (MFO), based on the navigation strategy of the moths in nature. According to [39], the MFO algorithm consists of three phases: initialization phase, updating phase, and stopping conditions. After distributing N moths within the search space of the problem using Equation (11), the fitness value for each is computed, and the best positions of the moths will be assigned to flames. Then, the optimization process will update the positions of moths based on the flames as defined in the following equation:
X i = F j + D i e b l cos ( 2 π l ) ,   i ,   j = 1 , 2 , 3 , ,   N
D i = | F j X i |
where l is a random number between 1 and −1, and b refers to the jth flame that is used to define the shape of the logarithmic spiral function. To avoid degrading the exploitation of the best solution, the number of flames (flame_no) must be decreased according to the adaptive mechanism strategy defined in Equation (24) to make the algorithm focus on the best-so-far solution in a hope of finding a better solution.
f l a m e _ n o = r o u n d ( N i t × N 1 t m a x )
The final steps of MFO with the PCM are presented in Algorithm 3. The standard MFO algorithm, at its beginning, will focus on searching for better solutions around the corresponding flame for each individual, then it will gradually move to focus on a number of the best-so-far solutions determined by Equation (24) even only focusing on the best-so-far one at the end. As in EO, this algorithm suffers from low exploitation capability and hence will waste several function evaluations at the start to search for the promising region, which involves the near-optimal solution for an optimization problem. A promising alternative solution to this problem is PCM, which works by exploiting the function evaluations as much as possible for reaching a better solution in less number of function evaluations.
Algorithm 3 The steps of PCM integrated MFO(PMFO)
1. Initialization step
2.  M F 1   = Calculate f ( X i ) for each moth.
3.  i t = 1
4. while ( i t < t m a x )
5.   Update f l a m e _ n o using Equation (24)
6.   if it==1
7.    F = sort X according to M F 1 ;
8.    FF = sort( M F 1 ); //sort the fitness values of moths
9.   Else
10.     M F i t   = Calculate f ( X i ) for each moth.
11.    F = sort ( X (it-1), X (it)) according to their fitness values found in
    M F i t 1 and M F i t
12.     FF = sort( M F i t 1 , M F i t );
13.   end
14.   for each X i ,   i = 1 , 2 , 3 , , N
15.       If   i f l a m e _ n o
16.         Update X i according to equation (22) with respect to its corresponding flame.
17.       else
18.         Update X i according to Equation (22) with respect to the flame   f l a m e _ n o .
19.       End if
20.        i t ++
21.   end for
22.   Update the current population using Algorithm 1.
23.    i t = i t + N
24.  end while

3.3.3. Harris Hawks Optimization Algorithm

Haidari et al. [38] proposed a meta-heuristic algorithm inspired by the chasing style and cooperative behaviors of Harris hawks, known as the HHO algorithm. The mathematical model of HHO simulates the behaviors of the hawks where several hawks cooperatively swoop on prey, often a rabbit, from different paths to surprise it. Additionally, Harris hawks have the ability to reveal different types of chasing patterns to choose the best one according to the distinct patterns of prey flight. In the exploration phase, HHO enables a trade-off between updating the Harris hawk’s perch randomly at a location near other members and wait to find prey, and perch on random tall trees with an equal probability of q as follows:
X ( i t + 1 ) = { X r ( i t ) r 1 | X r ( i t ) 2 r 2 X ( i t ) |                                                                                               q 0.5 X *   ( i t ) X m e a n ( i t ) r 3 ( L m i n + r 4 ( U m a x L m i n ) )           q < 0.5
where X (it) and X (it + 1) are two vectors including the current and the next position of the hawks. X r ( i t ) is a hawk selected randomly from the population. X*(it) is the location of the rabbit, which is also called the best-so-far solution. q, r 1 ,   r 2 ,   r 3 , and r 4 are five numerical values generated randomly. X m e a n ( i t ) indicates the mean position of solutions in the current population and computed as follows:
X m e a n ( i t ) = 1 N i = 1 N X i ( i t )
According to the escaping energy E of the rabbit, HHO can move from the exploration to the exploitation operator as defined in the following equation:
E = 2 E 0 ( 1 i t t m a x )
where E0 indicates the initial energy of the rabbit and is generated randomly between −1 and 1. If |E|≥1, hawks explore more regions to search for the rabbit position, otherwise it will exploit the current rabbit location. Based on the value of E, the hawks will create their step sizes according to a soft (|E| ≥ 0.5) or hard (|E| < 0.5) besiege. The soft besiege can be modeled as
X ( i t + 1 ) = Δ X ( i t ) E | J X *   ( i t ) X ( i t ) |
Δ X ( i t ) = X *   ( i t ) X ( i t )
J = 2 ( 1 r a n d )
J indicates the random jump strength of the rabbit; rand is a random number generated between 0 and 1. The hard besiege can be mathematically formulated as follows:
X ( i t + 1 ) = X *   ( i t ) E | Δ X ( i t ) |
where (|E| ≥ 0.5) and (p ≥ 0.5) in which the rabbit has enough energy to escape of the hawks, the soft besiege with progressive rapid dives (PRD) will be performed. Based on the actual behaviors of the hawks, they can choose the best possible dive toward the intended prey. In addition, in this phase, Lévy flight is used to mimic the escaping steps of the prey and the leapfrog of hawks. In this stage, the next position of the hawks is updated according to the following equation:
k = X *   ( i t ) E | J X *   ( i t ) X ( i t ) |
If this updated position represented in k is not better than the current position represented X(it + 1), the hawks will dive according to the Lévy flight L as follows [38]:
z = k + S × L ( d )
S is a numerical vector including random numbers. The final soft besiege can be summarized as
X ( i t + 1 ) = { k                               i f   f ( k ) < f ( X ( i t ) ) z                               i f   f ( z ) < f ( X ( i t ) )
In the hard besiege with PRD, the rabbit has too low energy to run away from hawks, when (|E| <0.5) and (p < 0.5) by using Equation (34) where Z is computed based on Equation (33) and k is updated according to the following:
X = X *   ( i t ) E | J X *   ( i t ) X m e a n ( i t ) |
Finally, the pseudo-code of the HHO hybridized with PCM is given in Algorithm 4. As in the EO and MFO, the HHO has a weak exploitation operator; therefore, due to the advantages of the PCM mentioned before, it is integrated with this algorithm to enhance its exploitation operator for attacking the promising regions obtained in each generation as an attempt to find the near-optimal solution without consuming several iterations.
Algorithm 4 The standard HHO algorithm with PCM (PHHO)
1. Initialization step.
2. Evaluate each hawk
3.  X * = best-so-far hawk
4.  i t = 1
5. while ( i t   t m a x )
6.    compute E according to Equation (27)
7.    if   ( | E | 1 )
8.         Update X(it + 1) using Equation (25)
9.   end if
10.    if   ( | E | < 1 )
11.      if ( p 0.5   & & | E | 0.5 )
12.           Reposition X(it + 1) based on soft besiege
13.      end if
14.      if ( p 0.5   & & | E | < 0.5 )
15.           Reposition X(it + 1) based on hard besiege
16.      end if
17.      if ( p < 0.5   & & | E | 0.5 )
18.           Reposition X(it + 1) based on soft besiege with PRD
19.      end if
20.      if ( p < 0.5   & & | E | 0.5 )
21.           Reposition X(it + 1) based on hard besiege with PRD
22.      end if
23.   end if
24.    if (f( X ( t + 1 )) < f( X r a b b i t ))
25.      Update X * = X(it + 1)
26.   end if
27.   i t = i t + N
28.     Update the current population using Algorithm 1.
29.      Update X * if there is better.
30.      it = it + N
31.  end while
Output: return X *

4. Results and Discussion

In this section, the influence of the PCM is observed with three well-known optimization algorithms: EO, HHO, and MFO. To validate the algorithms, the study used a PV cell based on RTC France (RTC) in addition to STM6-40/36 (STM6) module, Photowatt-PWP201 (PWP201) module, and STP6-120/36 (STP6) module as three PV modules. To illustrate the efficacy of the algorithms, several statistical metrics were used: worst, best, Avg, standard deviation (SD), and rank. In addition, Boxplots were used to depict the four quartiles of the outcomes obtained. Further, the convergence speed of each improved algorithm using PCM was compared with the standard version. For each test case, various demonstrations of experimental versus model results are presented, along with principal characteristics under different operating conditions such as varying temperatures and sun irradiances.
The experiments were conducted on a device with RAM of 32 GB, Core (TM) i7, and Windows 10. MATLAB R2019a was used to implement the algorithms. In total, 30 independent runs were carried out to determine the stability of the algorithms.

4.1. Datasets Description

The studied algorithms in this research were used to estimate the parameters of various photovoltaic (PV) models that include the SDM and PV models. For the RTC France Si cell, the measured I-V data are estimated at irradiance level (G) of 1000 W/m2 and temperature T = 33 °C [44]. On the other hand, for the PV modules, STM6, the PWP201, and Ultra 85-P modules are employed to validate the extracted parameters. The Photowatt-PWP201 module has Ns = 36 cells connected in series and the measured I-V test points are measured at G = 1000 W/m2 and T = 45 °C [44]. The mono-crystalline STM6 [45], and Ultra 85-P [46] modules consist of 36 monocrystalline PV cells connected in series and are measured under temperature levels of 51 °C, and 25 °C, respectively. The L m i n   and L m a x   of each known parameter according to the type of the PV cell and PV modules used in this study are given in Table 2, as described in the literature [13,16].

4.2. Parameter Selection

All the algorithms were executed using an equal number of function evaluations and runs t m a x = 50,000 to make a fair comparison among the algorithms. Picking the best value for N may significantly affect the performance of the algorithm; thus, different values for this parameter are observed to determine the best value with each algorithm. After observing 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, and 80 for N and depicting the outcomes in Figure 3, it is obvious that the performance of MFO is extremely poor when N is 5, but for the higher values, the performance almost is competitive; therefore, we randomly selected a population size of 80 in the next experiments since the population size higher than 5 does not significantly affect its performance. The performance of PMFO with a population size higher than 35, according to Figure 3b, is approximately the same, and therefore, 35 was selected as a population size for this algorithm within the next experiments. The best value for the population size of the HHO, PHHO, EO, and PEO according to Figure 3c–f is, respectively, 80, 50, 40, and 30 because they are competitive with the other values. From Figure 3, it is also obvious that the performance of all algorithms is significantly degraded with the small population size up to 10. Therefore, it is recommended that the population size for those investigated algorithms be assigned a value higher than 10 t to guarantee better performance.

4.3. RTC France

In this section, the standard algorithms and the PCM-improved algorithms were investigated on the SDM for the RTC France solar cell. It is worth mentioning that the improved algorithms in the following tables start with the letter P to be distinguished from the standard. Table 3 shows the optimal parameters estimated by each algorithm with the corresponding RMSE. In addition, the entire Table 4 shows the values of the statistical measures obtained by each algorithm, and Figure 4 shows the convergence speed, both of which indicate that the improved version of each algorithm outperforms the standard version for the different statistical measures as well as convergence speed: MFO has an average fitness value of 0.0012948678, while PMFO has the lower average of 0.0007731606; HHO has an average of 0.0027623201, while PHHO has a better average higher convergence speed; EO has an average fitness value of 0.0007736621, while PEO has a significant improvement with a value of 0.0007730063 and with higher convergence rate. It is worth noting that PEO is the algorithm with the lowest RMSE and better convergence speed. Based on this result, it is concluded that the PCM has a significant effect on the performance of the investigated meta-heuristic algorithms, and hence, this method considers a significant addition to the meta-heuristic algorithms for reaching better outcomes in less number function evaluations.
At this moment, various principal characteristics of this RTC Si PV cell can be generated using the cropped optimal values of the SDM parameters. Figure 5a,b illustrates the I–V and P–V plots of the PEO model against the real/experimental dataset points. On the other hand, various characteristics under varied sun irradiances (400–600–800–1000 W/m2) are revealed in Figure 5c,d, and under changeable cell, temperatures are shown in Figure 5e,f for 4 levels of temperatures at 0 °C to 75 °C in a step of 25 °C.

4.4. Photowatt-PWP201 Module

For this module, Table 5 presents the estimated parameters by each algorithm with the corresponding RMSE. Statistical measures are provided in Table 6, which shows that the performances of the algorithms integrated with PCM are significantly improved in terms of the final accuracy over the standard versions except for SCA that achieves a better outcome than the version integrated with PCM. The convergence speed of each standard algorithm against the improved version, in addition to the convergence speed among the improved versions, is given in Figure 6. From this figure, it can be inferred that PCM improved the classical algorithms for rapidly reaching better outcomes, compared to the classical ones, and consequently, this method is considered a good approach to overcome the convergence rate as the main shortcoming for most meta-heuristic algorithms. Again, this figure shows that PCM could accelerate the convergence speed of the classical EO in comparison to the other improved variants, as described also statistically in Table 5.
Once again, various principal characteristics of the PWP-201 module can be produced using the cropped optimal values of the SDM parameters. Figure 7a,b illustrates the I–V and P–V plots of the PEO model against the real/experimental dataset points. On the other hand, various characteristics under varied environmental conditions such as under varied sun irradiances (400–600–800–1000 W/m2) are revealed in Figure 7c,d, and under changeable cell, temperatures are shown in Figure 7e,f for 4 levels of temperatures at 0 °C to 75 °C in a step of 25 °C.

4.5. STM6-40/36 Module

Each algorithm was executed 30 independent times and the optimal parameter values obtained through those runs were introduced in Table 7 with the corresponding RMSE. Table 8 shows the statistical measures for the best, SD, Avg, and worst obtained in those independent runs. This table confirms the superiority of the improved algorithm over the standard version for the three algorithms investigated. Figure 8 shows the better convergence of PMFO, PHHO, and PEO over the standard versions and also clarifies the competitivity between PMFO and PEO.
Similar to the abovementioned test cases, many principal characteristics of the STP6-40 module can be produced using the cropped optimal values of the SDM parameters. Figure 9a,b illustrates the I–V and P–V plots of the PEO model against the real/experimental dataset points. On the other hand, various characteristics under varied environmental conditions such as under varied sun irradiances (400–600–800–1000 W/m2) are revealed in Figure 9c,d, and under changeable cell, temperatures are shown in Figure 9e,f for 4 levels of temperatures at 0 °C to 75 °C in a step of 25 °C.

4.6. Ultra 85-P Module

In this section, a new commercial module, called the Ultra 85-P module, was used to check the performance of the algorithms. The typical nameplate of this module is shown in [46], where it consists of 36 PV cells connected in series and could generate a maximum power of 85 W at standard conditions. This module has an efficiency of 13.4% and 70.3% as a fill factor. To check the performance of the algorithms under this module, Each algorithm was executed 30 independent times and the optimal parameter values and the corresponding RMSE obtained through those runs are presented in Table 9, which show that PEO could reach a value of 0.002551066 for RMSE as the lowest one over the others. Furthermore, Table 10 shows the statistical measures for the best, SD, Avg, and worst obtained in those independent runs. This table confirms the superiority of PEO over the standard and the other algorithms, but unfortunately, both PHHO and PMFO could not exceed the standard ones. In general, our experiments turn out that PCM could significantly improve the performance of the standard EO, and this affirms its efficiency when integrating with some algorithms as an aiding tool to explore some regions, which are intractable by those standard ones. Figure 10a shows the better convergence of PEO, compared to the others. Figure 10b,c illustrates the P–V and I–V plots of the PEO model versus the real/experimental dataset points.

4.7. Comparison of the Studied Algorithms Using Boxplot

In Figure 11, the studied algorithms are compared by drawing the boxplot of THE PWP201, and STM6 modules, and the RTC France solar cell. After running each algorithm for 30 independent runs and depicting the obtained outcomes for those PV modules in Figure 11, it is obvious that the improved algorithms, i.e., PHHO, PEO, and PMFO, outperform their standard versions in all cases. On the other hand, PEO is considered the best option in comparison to all improved variants using PCM, and this shows that hybridization between the classical EO and this method manages to build a new variant having a high ability for widely effective exploration of the search space to reach better outcomes in less number of function evaluations, compared to the classical one.

4.8. Wilcoxon Rank-Sum Test

The results obtained within 30 independent runs of each algorithm (original and improved) were compared using the Wilcoxon rank-sum test [47] at a confidence level of 5%. In Table 11, the p-value illustrates the extent of differences in the outcomes obtained by each pair of the algorithms; the h includes only two values: 0 indicates when there is no difference between the outcomes of a pair of the algorithms, and 1 when there is a difference. According to Table 11, for RTC France, the p-value obtained is less than 0.05 for all the improved algorithms against the standard versions, which confirms that there is a difference between the outcomes obtained by each pair of those algorithms. Based on this statistical test and the statistical analysis given in the previous subsections, the improved versions can be accepted as the best option in a significant number of the studied cases. In the same context, Table 12 which compares the difference between the outcomes of PEO as the best-improved variant and all the others under the Wilcoxon rank-sum test, confirms that PEO is the best since it could obtain outcomes significantly different from those of the others.

5. Conclusions

This paper proposes a new strategy known as the premature convergence method (PCM) in order to accelerate the convergence speed of meta-heuristic algorithms while improving the final accuracy of the optimization algorithms. PCM updates the current individual around the best-so-far solution based on two-step sizes: the first is based on the distance between two individuals selected randomly from the population, while the second is based on the distance between the current solution and the best-so-far solution. For the trade-off between those two steps, a weight variable is used to determine the length of each step size that was added to the best-so-far solution. This weight variable is generated randomly between 0 and 1. Additionally, when the value in this variable is greater than 0.5, emphasis is placed on the second step, and this will increase the exploitation operator of the optimization algorithm. However, if the value of this variable is small, then the exploration capability will be encouraged.
The proposed PCM was integrated with three well-known optimization algorithms—HHO, MFO, and EO—to observe its effectiveness in improving those algorithms’ ability to find the unidentified parameters of SDM. After investigating the performance of the improved optimizers on the SDM and PV module model, it is obvious that PCM has a significant effect on the performance of the optimization algorithms, especially EO, for observed cases. The numerical findings obtained by PEO for the observed cases represented in RTC France, PWP201 module, Ultra 85-P, and STM6 module, respectively, are 0.0007730063, 0.0020529606, 0.0025510660, and 0.0017219215. Future work includes improving the control factor of this method to balance exploration and exploitation capability. In addition, we will investigate the performance of PCM with other algorithms that have a balance between the exploration and exploitation at the beginning of the optimization process to investigate if it could improve their performance.

Author Contributions

Conceptualization, M.A.-B., R.M., A.E.-F. and M.A.; methodology, M.A.-B., R.M. and M.A.; software, M.A.-B., R.M., M.A. and A.E.-F.; validation, M.A., A.E.-F.; formal analysis, M.A.-B., R.M. and M.A.; investigation, Y.N., A.E.-F. and M.A.; resources, M.A.-B., M.A. and R.M.; data curation, M.A.-B., A.E.-F., R.M. and M.A.; writing—original draft preparation, M.A.-B., R.M. and M.A.; writing—review and editing, Y.N., A.E.-F. and M.A.; visualization, M.A.-B., M.A., A.E.-F. and R.M.; supervision, M.A.-B., M.A. and A.E.-F.; project administration, M.A.-B., R.M. and M.A.; funding acquisition, Y.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Soonchunhyang University Research Fund.

Institutional Review Board Statement

The study did not involve humans or animals.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ayala, H.V.H.; Coelho, L.D.S.; Mariani, V.C.; Askarzadeh, A. An improved free search differential evolution algorithm: A case study on parameters identification of one diode equivalent circuit of a solar cell module. Energy 2015, 93, 1515–1522. [Google Scholar] [CrossRef]
  2. Oliva, D.; Cuevas, E.; Pajares, G. Parameter identification of solar cells using artificial bee colony optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  3. Awadallah, M.A. Variations of the bacterial foraging algorithm for the extraction of PV module parameters from nameplate data. Energy Convers. Manag. 2016, 113, 312–320. [Google Scholar] [CrossRef]
  4. AlHajri, M.; El-Naggar, K.; AlRashidi, M.; Al-Othman, A. Optimal extraction of solar cell parameters using pattern search. Renew. Energy 2012, 44, 238–245. [Google Scholar] [CrossRef]
  5. Yeh, W.; Huang, C.; Lin, P.; Chen, Z.; Jiang, Y.; Sun, B. Simplex simplified swarm optimisation for the efficient optimisation of parameter identification for solar cell models. IET Renew. Power Gener. 2017, 12, 45–51. [Google Scholar] [CrossRef]
  6. Kiani, A.T.; Nadeem, M.F.; Ahmed, A.; Khan, I.; Elavarasan, R.M.; Das, N. Optimal PV Parameter Estimation via Double Exponential Function-Based Dynamic Inertia Weight Particle Swarm Optimization. Energies 2020, 13, 4037. [Google Scholar] [CrossRef]
  7. Li, G.; Chen, X.; Jin, Y. Analysis of the Primary Constraint Conditions of an Efficient Photovoltaic-Thermoelectric Hybrid System. Energies 2016, 10, 20. [Google Scholar] [CrossRef] [Green Version]
  8. Ullah, I.; Rasul, M.G. Recent Developments in Solar Thermal Desalination Technologies: A Review. Energies 2018, 12, 119. [Google Scholar] [CrossRef] [Green Version]
  9. Jelle, B.P. Building Integrated Photovoltaics: A Concise Description of the Current State of the Art and Possible Research Pathways. Energies 2015, 9, 21. [Google Scholar] [CrossRef] [Green Version]
  10. Long, W.; Cai, S.; Jiao, J.; Xu, M.; Wu, T. A new hybrid algorithm based on grey wolf optimizer and cuckoo search for parameter extraction of solar photovoltaic models. Energy Convers. Manag. 2020, 203, 112243. [Google Scholar] [CrossRef]
  11. Liang, J.; Ge, S.; Qu, B.; Yu, K.; Liu, F.; Yang, H.; Wei, P.; Li, Z. Classified perturbation mutation based particle swarm optimization algorithm for parameters extraction of photovoltaic models. Energy Convers. Manag. 2020, 203, 112138. [Google Scholar] [CrossRef]
  12. Li, S.; Gu, Q.; Gong, W.; Ning, B. An enhanced adaptive differential evolution algorithm for parameter extraction of photovoltaic models. Energy Convers. Manag. 2020, 205. [Google Scholar] [CrossRef]
  13. Chen, X.; Xu, B.; Mei, C.; Ding, Y.; Li, K. Teaching–learning–based artificial bee colony for solar photovoltaic parameter estimation. Appl. Energy 2018, 212, 1578–1588. [Google Scholar] [CrossRef]
  14. Ortiz-Conde, A.; Garcia-Sanchez, F.; Muci, J. New method to extract the model parameters of solar cells from the explicit analytic solutions of their illuminated I–V characteristics. Sol. Energy Mater. Sol. Cells 2006, 90, 352–361. [Google Scholar] [CrossRef]
  15. Tong, N.T.; Kamolpattana, K.; Pora, W. A deterministic method for searching the maximum power point of a PV panel. In Proceedings of the 2015 12th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology (ECTI-CON), Hua Hin, Thailand, 24–27 June 2015. [Google Scholar]
  16. Yu, K.; Liang, J.; Qu, B.; Chen, X.; Wang, H. Parameters identification of photovoltaic models using an improved JAYA optimization algorithm. Energy Convers. Manag. 2017, 150, 742–753. [Google Scholar] [CrossRef]
  17. Allaoui, M.; Ahiod, B.; El Yafrani, M. A hybrid crow search algorithm for solving the DNA fragment assembly problem. Expert Syst. Appl. 2018, 102, 44–56. [Google Scholar] [CrossRef]
  18. Abdel-Basset, M.; Mohamed, R.; Elhoseny, M.; Bashir, A.K.; Jolfaei, A.; Kumar, N. Energy-Aware Marine Predators Algorithm for Task Scheduling in IoT-based Fog Computing Applications. IEEE Trans. Ind. Inform. 2020, 17, 5068–5076. [Google Scholar] [CrossRef]
  19. Abdel-Basset, M.; Chang, V.; Mohamed, R. HSMA_WOA: A hybrid novel Slime mould algorithm with whale optimization algorithm for tackling the image segmentation problem of chest X-ray images. Appl. Soft Comput. 2020, 95, 106642. [Google Scholar] [CrossRef]
  20. Zhou, Y.; Yen, G.G.; Yi, Z. A Knee-Guided Evolutionary Algorithm for Compressing Deep Neural Networks. IEEE Trans. Cybern. 2019, 51, 1626–1638. [Google Scholar] [CrossRef]
  21. Zhang, J.; Zhang, T.; Shin, H.-S.; Wang, J.; Zhang, C. Geomagnetic Gradient-Assisted Evolutionary Algorithm for Long-Range Underwater Navigation. IEEE Trans. Instrum. Meas. 2021, 70, 2503212. [Google Scholar] [CrossRef]
  22. Xiao, D.; Prado, J.C.D.; Qiao, W. Optimal joint demand and virtual bidding for a strategic retailer in the short-term electricity market. Electr. Power Syst. Res. 2021, 190, 106855. [Google Scholar] [CrossRef]
  23. Xiao, D.; Prado, J.C.d.; Qiao, W. Winner-leading competitive swarm optimizer with dynamic Gaussian mutation for parameter extraction of solar photovoltaic models. Electr. Power Syst. Res. 2021, 190, 106855. [Google Scholar] [CrossRef]
  24. Diab, A.A.Z.; Sultan, H.M.; Do, T.D.; Kamel, O.M.; Mossa, M.A. Coyote Optimization Algorithm for Parameters Estimation of Various Models of Solar Cells and PV Modules. IEEE Access 2020, 8, 111102–111140. [Google Scholar] [CrossRef]
  25. Ridha, H.M.; Heidari, A.A.; Wang, M.; Chen, H. Boosted mutation-based Harris hawks optimizer for parameters identification of single-diode solar cell models. Energy Convers. Manag. 2020, 209, 112660. [Google Scholar] [CrossRef]
  26. Long, W.; Wu, T.; Jiao, J.; Tang, M.; Xu, M. Refraction-learning-based whale optimization algorithm for high-dimensional problems and parameter estimation of PV model. Eng. Appl. Artif. Intell. 2020, 89, 103457. [Google Scholar] [CrossRef]
  27. Abdel-Basset, M.; Mohamed, R.; Mirjalili, S.; Chakrabortty, R.K.; Ryan, M.J. Solar photovoltaic parameter estimation using an improved equilibrium optimizer. Sol. Energy 2020, 209, 694–708. [Google Scholar] [CrossRef]
  28. Ridha, H.M.; Gomes, C.; Hizam, H. Estimation of photovoltaic module model’s parameters using an improved electromagnetic-like algorithm. Neural Comput. Appl. 2020, 32, 12627–12642. [Google Scholar] [CrossRef]
  29. Ram, J.P.; Pillai, D.S.; Rajasekar, N.; Chinnaiyan, V.K. Flower Pollination Based Solar PV Parameter Extraction for Double Diode Model. In Intelligent Computing Techniques for Smart Energy Systems; Springer: Berlin/Heidelberg, Germany, 2019; pp. 303–312. [Google Scholar] [CrossRef]
  30. Hassan, K.H.; Rashid, A.T.; Jasim, B.H. Parameters estimation of solar photovoltaic module using camel behavior search algorithm. Int. J. Electr. Comput. Eng. 2021, 11, 788–793. [Google Scholar] [CrossRef]
  31. Kashefi, H.; Sadegheih, A.; Mostafaeipour, A.; Omran, M.M. Parameter identification of solar cells and fuel cell using improved social spider algorithm. COMPEL Int. J. Comput. Math. Electr. Electron. Eng. 2020, 40, 142–172. [Google Scholar] [CrossRef]
  32. Sharma, A.; Sharma, A.; Averbukh, M.; Jately, V.; Azzopardi, B. An Effective Method for Parameter Estimation of a Solar Cell. Electronics 2021, 10, 312. [Google Scholar] [CrossRef]
  33. Ye, X.; Liu, W.; Li, H.; Wang, M.; Chi, C.; Liang, G.; Chen, H.; Huang, H. Modified Whale Optimization Algorithm for Solar Cell and PV Module Parameter Identification. Complexity 2021, 2021, 1–23. [Google Scholar] [CrossRef]
  34. Yang, X.; Gong, W. Opposition-based JAYA with population reduction for parameter estimation of photovoltaic solar cells and modules. Appl. Soft Comput. 2021, 104, 107218. [Google Scholar] [CrossRef]
  35. Mokeddem, D. Parameter Extraction of Solar Photovoltaic Models Using Enhanced Levy Flight Based Grasshopper Optimization Algorithm. J. Electr. Eng. Technol. 2020, 16, 171–179. [Google Scholar] [CrossRef]
  36. Ismaeel, A.A.K.; Houssein, E.H.; Oliva, D.; Said, M. Gradient-Based Optimizer for Parameter Extraction in Photovoltaic Models. IEEE Access 2021, 9, 13403–13416. [Google Scholar] [CrossRef]
  37. Sharma, A.; Saxena, A.; Shekhawat, S.; Kumar, R.; Mathur, A. Solar Cell Parameter Extraction by Using Harris Hawks Optimization Algorithm; Springer: Singapore, 2020; pp. 349–379. [Google Scholar] [CrossRef]
  38. Heidari, A.A.; Mirjalili, S.; Faris, H.; Aljarah, I.; Mafarja, M.; Chen, H. Harris hawks optimization: Algorithm and applications. Futur. Gener. Comput. Syst. 2019, 97, 849–872. [Google Scholar] [CrossRef]
  39. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl. Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  40. Faramarzi, A.; Heidarinejad, M.; Stephens, B.; Mirjalili, S. Equilibrium optimizer: A novel optimization algorithm. Knowl. Based Syst. 2020, 191, 105190. [Google Scholar] [CrossRef]
  41. Tan, Y.T.; Kirschen, D.; Jenkins, N. A Model of PV Generation Suitable for Stability Analysis. IEEE Trans. Energy Convers. 2004, 19, 748–755. [Google Scholar] [CrossRef]
  42. Gong, W.; Cai, Z. Parameter extraction of solar cell models using repaired adaptive differential evolution. Sol. Energy 2013, 94, 209–220. [Google Scholar] [CrossRef]
  43. Nunes, H.G.G.; Pombo, J.; Mariano, S.; Calado, M.D.R.; de Souza, J.F. A new high performance method for determining the parameters of PV cells and modules based on guaranteed convergence particle swarm optimization. Appl. Energy 2018, 211, 774–791. [Google Scholar] [CrossRef]
  44. Easwarakhanthan, T.; Bottin, J.; Bouhouch, I.; Boutrit, C. Nonlinear Minimization Algorithm for Determining the Solar Cell Parameters with Microcomputers. Int. J. Sol. Energy 1986, 4, 1–12. [Google Scholar] [CrossRef]
  45. Tong, N.T.; Pora, W. A parameter extraction technique exploiting intrinsic properties of solar cells. Appl. Energy 2016, 176, 104–115. [Google Scholar] [CrossRef] [Green Version]
  46. Elkholy, M.M.; El-Hameed, M.A.; El-Fergany, A.A. Artificial ecosystem-based optimiser to electrically characterise PV generating systems under various operating conditions reinforced by experimental validations. IET Renew. Power Gener. 2021, 15, 701–715. [Google Scholar] [CrossRef]
  47. Haynes, W. Wilcoxon Rank Sum Test. In Encyclopedia of Systems Biology; Springer: New York, NY, USA, 2013; pp. 2354–2355. [Google Scholar]
Figure 1. Equivalent circuit of SDM.
Figure 1. Equivalent circuit of SDM.
Electronics 10 01846 g001
Figure 2. Flowchart of hybridization between PCM and EO.
Figure 2. Flowchart of hybridization between PCM and EO.
Electronics 10 01846 g002
Figure 3. Tuning of N parameter.
Figure 3. Tuning of N parameter.
Electronics 10 01846 g003aElectronics 10 01846 g003b
Figure 4. Convergence curves among algorithms in RTC France.
Figure 4. Convergence curves among algorithms in RTC France.
Electronics 10 01846 g004
Figure 5. Principal characteristics of RTC France Si solar cell.
Figure 5. Principal characteristics of RTC France Si solar cell.
Electronics 10 01846 g005
Figure 6. Convergence curves among algorithms in Photowatt-PWP201.
Figure 6. Convergence curves among algorithms in Photowatt-PWP201.
Electronics 10 01846 g006aElectronics 10 01846 g006b
Figure 7. Principal characteristics of PWP201 unit.
Figure 7. Principal characteristics of PWP201 unit.
Electronics 10 01846 g007aElectronics 10 01846 g007b
Figure 8. Convergence curves among algorithms in STM6.
Figure 8. Convergence curves among algorithms in STM6.
Electronics 10 01846 g008
Figure 9. Principal characteristics of STM6-40/36PV module.
Figure 9. Principal characteristics of STM6-40/36PV module.
Electronics 10 01846 g009aElectronics 10 01846 g009b
Figure 10. Graphical depictions of the performance of the algorithm in Ultra 85-P.
Figure 10. Graphical depictions of the performance of the algorithm in Ultra 85-P.
Electronics 10 01846 g010
Figure 11. Comparison among algorithms based on Boxplot.
Figure 11. Comparison among algorithms based on Boxplot.
Electronics 10 01846 g011
Table 1. Reviews of some studied algorithms for the PIP of different PV models.
Table 1. Reviews of some studied algorithms for the PIP of different PV models.
Algorithm and YearContributions and Limitations
Classified Perturbation Mutation Based PSO Algorithm (CPMPSO, 2020) [11]This algorithm divided the individuals into two categories according to the fitness values: The first one has individuals with high quality and is updated using an effective exploitation operator, while in the other, individuals had been updated using an effective exploration operator.
Its convergence speed still needs improvement. Additionally, it has difficulty in avoiding becoming trapped in local minima for DDM.
Enhanced Adaptive Differential Evolution (EJADE, 2020) [12]This algorithm used a number of improvements, namely, a crossover sorting mechanism for using the best individuals in the next generation to reach better outcomes, and a dynamic population reduction strategy to increase the convergence speed.
Whale Optimization Algorithm (WOA) based Reflecting Learning (RLWOA, 2020) [26]WOA was improved using the reflection learning strategy to reduce the probability of becoming trapped in local minima, and subsequently increasing the possibility of reaching better outcomes. However, the increased speed of this algorithm against the studied algorithms was not analyzed.
Improved Equilibrium Optimizer (IEO, 2020) [27]In this paper, EO was improved using two strategies: the first worked on accelerating the convergence, while the second was used to avoid becoming trapped into local minima. This algorithm produced good outcomes compared to four compared algorithms on three different PV models.
Improved Electromagnetism-like (IEM, 2020) [28]A nonlinear equation was used to adjust the number of individuals in each generation to increase the convergence speed. Simplifying the total force formula to increase the exploration operator to explore the most promising regions for avoiding becoming stuck in local minima problems. Only validated on SDM, its performance is not known for DDM and TDM.
Flower Pollination Algorithm (FPA, 2020) [29]The authors adapted FPA for estimating the parameters of DDM and used RTC France to validate its performance. Additionally, for verifying its performance, it was extensively compared with four studied algorithms. Unknown performance was compared with some of the recent robust algorithms published within the last two years.
Camel Behavior Search Algorithm (CBSA, 2020) [30]For estimating the parameters of SDM for the multi-crystalline KC 200GT PV module, CBSA was proposed. This algorithm was validated on the SDM of the PV solar module, but its final outcome and convergence speed still need significant improvement.
Improved Social Spider Algorithm (ISSA, 2020) [31]In this research, the social spider algorithm was proposed with an improvement in its performance to increase its exploration operator; this improvement was based on replacing the worst individuals within the populations with other solutions within the search space of the problem after a period of the iteration.
A hybrid WOA and PSO Algorithm (HWOA, 2021) [32]In this paper, a new parameter estimation algorithm based on integrating PSO with WOA and a pipeline model was proposed to accurately speed the convergence rate. The experimental outcomes affirm that this algorithm was better than all the compared for the convergence rate and accuracy.
A Modified Whale Optimization Algorithm (MWOA, 2021) [33]In this paper, MWOA was proposed to overcome stagnation into local minima, and low convergence speed by employing a mutation operator based on the levy flight, and a local search strategy to promote the exploitation capability. Thereafter, this algorithm was employed for tackling the parameter estimation of the PV models and could fulfill superior performance.
An enhanced JAYA (EJAYA, 2021) [34]Recently, a new variant of the JAYA algorithm, namely, EJAYA, has been developed to improve the standard algorithm using three effective improvements. EJAYA could be the best solution, compared to some related techniques.
Enhanced Levy Flight Based Grasshopper Optimization Algorithm (LGOA, 2021) [35]The levy flight was integrated into the grasshopper optimization algorithm to utilize its advantages in preserving the diversity among the populations in addition to enhancing the exploitation capability for proposing a new solar cell parameter estimation technique named LGOA.
Gradient-Based Optimizer (GBO, 2021) [36]To extract five, seven, and nine unknown parameters of SDM, DDM, and TDM, respectively, the gradient-based optimizer was recently proposed for tackling the global optimization problem, for which it was adapted due to having a high convergence speed with a highly local minima avoidance strategy. The experimental outcomes show the proficiency of this developed algorithm.
Harris Hawks Optimization Algorithm (HHO, 2020) [37]The application of HHO to estimate the unknown parameters of the PV models has been recently proposed to examine its efficiency in comparison to some of the other optimization algorithms. The experimental findings show the efficiency of HHO over the compared ones.
Table 2. The search boundaries of each unknown parameter.
Table 2. The search boundaries of each unknown parameter.
ParameterRTCPWP201STM6Ultra 85-P
L m i n L m a x L m i n L m a x L m i n L m a x L m i n L m a x
I p h (A)0102024.95.73
I s d (μA)010500501.0 × 10−30.5
R s (Ω)00.50200.361.0 × 10−65
R s h (Ω)010002000015000.93100
n , n 1 , n 2 121501601.02.0
Table 3. Comparison under the extracted parameters and the corresponding RMSE of RTC France.
Table 3. Comparison under the extracted parameters and the corresponding RMSE of RTC France.
Algorithms I p h ( A ) I d ( A ) R s ( Ω ) R s h ( Ω ) n RMSE
MFO [39]0.760625.03 × 10−70.0343966.637601.527290.0007758336
PMFO0.760703.84 × 10−70.0356258.120591.498930.0007730647
HHO [38]0.762357.59 × 10−70.0318847.875491.573210.0007801742
PHHO0.766437.91 × 10−70.0304024.385381.579160.0008404323
EO [40]0.759924.48 × 10−70.0351380.931811.514760.0007731125
PEO0.760793.11 × 10−70.0365552.889791.477270.0007730063
Bold results are the best.
Table 4. Comparison of statistical measures of RTC France.
Table 4. Comparison of statistical measures of RTC France.
Algorithms B e s t W o r s t A v g S D Rank
MFO [39]0.00077873320.00236312630.00129486784.2011 × 10−44
PMFO0.00077300630.00077627730.00077316065.9821 × 10−72
HHO [38]0.00080654670.00726218000.00276232011.7806 × 10−36
PHHO0.00077685500.00462309020.00144320968.0094 × 10−45
EO [40]0.00077366210.00113823740.00083801519.7150 × 10−53
PEO0.00077300630.00077300630.00077300631.2633 × 10−171
Bold results are the best option.
Table 5. Comparison under the extracted parameters per cell and the corresponding RMSE of Photowatt-PWP201.
Table 5. Comparison under the extracted parameters per cell and the corresponding RMSE of Photowatt-PWP201.
Algorithms B e s t W o r s t A v g S D Rank
MFO [39]0.00218605410.00900075760.00645579761.873147 × 10−34
PMFO0.00205296060.00220826980.00208131264.057440 × 10−52
HHO [38]0.00424718280.04178593790.01774977599.936836 × 10−36
PHHO0.00214881750.02676750290.00876835855.384805 × 10−35
EO [40]0.00211809250.00365056190.00288875864.197005 × 10−43
PEO0.00205296060.00205296060.00205296062.3811 × 10−171
Bold results are the best option.
Table 6. Comparison of statistical measures of Photowatt-PWP201.
Table 6. Comparison of statistical measures of Photowatt-PWP201.
Algorithms I p h ( A ) I d ( A ) R s ( Ω ) R s h ( Ω ) n RMSE
MFO [39]1.032961.82 × 10−61.27671656.491741.285480.0021860541
PMFO1.031432.64 × 10−61.23563821.647421.322170.0020529606
HHO [38]1.031371.09 × 10−51.050251595.179231.485180.0042471828
PHHO1.030272.29 × 10−61.25910893.380271.307850.0021488175
EO [40]1.029563.23 × 10−61.215091090.973551.342990.0021180925
PEO1.031432.64 × 10−61.23563821.641291.322170.0020529606
Bold results are the best option.
Table 7. Comparison under the extracted parameters and the corresponding RMSE of STM6.
Table 7. Comparison under the extracted parameters and the corresponding RMSE of STM6.
Algorithms I p h ( A ) I d ( A ) R s ( Ω ) R s h ( Ω ) n RMSE
MFO [39]1.662293.36 × 10−60.07400726.091281.596510.0023187518
PMFO1.663911.74 × 10−60.15382573.274391.520300.0017219251
HHO [38]1.661271.70 × 10−60.18140718.979101.517660.0026097874
PHHO1.666911.21 × 10−60.18668462.622091.481700.0021550622
EO [40]1.663072.11 × 10−60.13187624.290601.542030.0017870003
PEO1.663901.74 × 10−60.15364573.533911.520470.0017219215
Bold results are the best option.
Table 8. Comparison of statistical measures of STM6.
Table 8. Comparison of statistical measures of STM6.
Algorithms B e s t W o r s t A v g S D Rank
MFO [39]0.00231875180.02255365980.01147080646.7323 × 10−35
PMFO0.00172192510.00634364120.00214898291.1362 × 10−33
HHO [38]0.00260978740.06132816130.02651074921.6809 × 10−26
PHHO0.00215506220.04655987080.01071523181.1196 × 10−24
EO [40]0.00178700030.00308415040.00246561013.2718 × 10−42
PEO0.00172192150.00172192150.00172192155.2394 × 10−181
Bold results are the best option.
Table 9. Comparison under the extracted parameters and the corresponding RMSE of Ultra 85-P.
Table 9. Comparison under the extracted parameters and the corresponding RMSE of Ultra 85-P.
Algorithms I p h ( A ) I d ( A ) R s ( Ω ) R s h ( Ω ) n RMSE
MFO [39]5.254032.231 × 10−60.012432.472701.405010.011839298
PMFO5.722501.699 × 10−70.2474599.98871.360730.123151236
HHO [38]5.198322.763 × 10−50.010098.476941.692100.009073968
PHHO5.230422.416 × 10−50.010104.368821.674310.009076382
EO [40]5.218721.292 × 10−50.010854.417491.593880.003321018
PEO5.227071.043 × 10−50.011043.799721.568220.002551066
Bold results are the best option.
Table 10. Comparison of statistical measures of STM6.
Table 10. Comparison of statistical measures of STM6.
Algorithms B e s t W o r s t A v g S D Rank
MFO [39]0.01183929840.16327199090.12777402130.02335295375
PMFO0.12315123630.17896392160.14338827560.01385649626
HHO [38]0.00907396830.20073337080.07787402550.06809145502
PHHO0.00907638220.24531385810.09653940550.06880377073
EO [40]0.00332101800.12346499190.09680280580.04867081864
PEO0.00255106600.12314832920.07490942390.06009034561
Bold results are the best option.
Table 11. Comparison between the proposed with the others under the Mann–Whitney U test.
Table 11. Comparison between the proposed with the others under the Mann–Whitney U test.
RTCPWP201STM6
Algorithmshp-Valuehp-Valuehp-Value
MFO vs. PMFO11.114256 × 10−312.904721 × 10−114.615910 × 10−10
HHO vs. PHHO16.843226 × 10−417.012659 × 10−212.126464 × 10−4
EO vs. PEO13.017967 × 10−1113.010407 × 10−1113.006634 × 10−11
Table 12. Comparison between PEO with the others under the Mann–Whitney U test.
Table 12. Comparison between PEO with the others under the Mann–Whitney U test.
RTCPWP201STM6Ultra 85-P
Algorithmshp-Valuehp-Valuehp-Valuehp-Value
MFO13.019859 × 10−1113.019859 × 10−1113.014185 × 10−1111.76972 × 10−10
PMFO13.019859 × 10−1113.338389 × 10−1113.014185 × 10−1113.00663 × 10−11
HHO13.019859 × 10−1113.019859 × 10−1113.014185 × 10−1113.38569 × 10−2
PHHO13.019859 × 10−1113.019859 × 10−1113.014185 × 10−1117.95291 × 10−3
EO13.019859 × 10−1113.019859 × 10−1113.014185 × 10−1111.72503 × 10−6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdel-Basset, M.; Mohamed, R.; Abouhawwash, M.; Nam, Y.; El-Fergany, A. Recent Meta-Heuristic Algorithms with a Novel Premature Covergence Method for Determining the Parameters of PV Cells and Modules. Electronics 2021, 10, 1846. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics10151846

AMA Style

Abdel-Basset M, Mohamed R, Abouhawwash M, Nam Y, El-Fergany A. Recent Meta-Heuristic Algorithms with a Novel Premature Covergence Method for Determining the Parameters of PV Cells and Modules. Electronics. 2021; 10(15):1846. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics10151846

Chicago/Turabian Style

Abdel-Basset, Mohamed, Reda Mohamed, Mohamed Abouhawwash, Yunyoung Nam, and Attia El-Fergany. 2021. "Recent Meta-Heuristic Algorithms with a Novel Premature Covergence Method for Determining the Parameters of PV Cells and Modules" Electronics 10, no. 15: 1846. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics10151846

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