Next Article in Journal
Thermographic Measurement of the Temperature of Reactive Power Compensation Capacitors
Next Article in Special Issue
Power Quality Improvement in HRES Grid Connected System with FOPID Based Atom Search Optimization Technique
Previous Article in Journal
Generation Expansion Planning with Energy Storage Systems Considering Renewable Energy Generation Profiles and Full-Year Hourly Power Balance Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Extraction of Photovoltaic Cells and Modules Using Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy

by
Mehmet Yesilbudak
Department of Electrical and Electronics Engineering, Faculty of Engineering and Architecture, Nevsehir Haci Bektas Veli University, Nevsehir 50300, Turkey
Energies 2021, 14(18), 5735; https://doi.org/10.3390/en14185735
Submission received: 2 August 2021 / Revised: 2 September 2021 / Accepted: 6 September 2021 / Published: 12 September 2021

Abstract

:
With the increase in the share of solar energy in the sustainable development, accurate parameter identification plays a significant role in designing optimal solar photovoltaic systems. For this purpose, this paper extensively implements and evaluates the grey wolf optimizer with a dimension learning-based hunting search strategy, an improved version of GWO named I-GWO, in the parameter extraction of photovoltaic cells and modules. According to the experimental results, the double-diode model leads to better fitness than the other diode models in representing the physical behaviors of both photovoltaic cells and photovoltaic modules. For further performance validation, firstly, the internal parameters extracted by the I-GWO algorithm and the corresponding output current data are compared with a number of widely-used parameter extraction methods in the literature. Then, the best goodness-of-fit results achieved by the I-GWO algorithm are evaluated considering many state-of-the-art metaheuristic algorithms in the literature. The accuracy measures including root mean squared error and sum of individual absolute errors show that I-GWO is fairly promising to be the efficient and valuable parameter extraction method for both photovoltaic cells and photovoltaic modules.

1. Introduction

1.1. Overview

By the end of 2019, the overall electricity generation from bioenergy, geothermal energy, hydropower, ocean power, solar energy and wind energy systems reached about 2588 GW worldwide. Among these renewable energy sources, solar photovoltaic systems have once again emerged as the frontrunner with the additional capacity of around 115 GW in 2019 [1]. In solar photovoltaic systems, the accurate characterization of diode model parameters enables us to closely represent the nonlinear current–voltage (IV) and power–voltage (PV) curves of photovoltaic cells and modules. Therefore, the solution of parameter extraction issue is a fundamental requirement for proper design, performance evaluation, quality controls and optimal operation of photovoltaic systems [2,3]. However, this task is still complex and challenging under various operating conditions.

1.2. Literature Review

In order to tackle the parameter extraction problem of photovoltaic cells and modules, a number of techniques have been recommended in the literature. These techniques are mostly categorized into three groups: analytical methods, deterministic techniques and metaheuristic optimization algorithms [4,5].
The analytical methods are mainly based on the correctness of maximum power voltage, maximum power current, short circuit current, open circuit voltage, etc. points on the IV curve [4,6]. If these key data points are erroneously identified, the modeling accuracy drastically worsens. Instead of selecting several key data points on the IV curve, deterministic techniques and metaheuristic optimization algorithms consider all measured IV data [4,7]. Deterministic techniques, such as Levenberg–Marquardt method [8], Nelder-Mead simplex method [9], Gauss–Seidel method [10], Runge–Kutta–Merson method [11], Lambert W functions [12], least squares [13], etc., impose the limitations on the differentiability, convexity and continuity characteristics of fitness functions. This condition causes convergence to the local optima in the solution of the multimodal, multivariable and nonlinear parameter extraction issue. An analytical method based on the ideality factor variation was used to identify the electrical parameters and the absolute error of output current was found below 0.25 A for a solar module [14]. A simple iterative method was adopted to obtain the initial values of shunt resistance and its relative error did not exceed 4.38% [15]. The Newton–Raphson method was utilized to estimate the reverse saturation current and diode ideality factor, and the deviation in open circuit voltage was found 0.08 V for a solar cell [16].
Unlike the deterministic techniques, metaheuristic algorithms do not impose any limitations on the optimization problem and they have the capability of solving various complex issues. A genetic algorithm was used to improve the parameter estimation accuracy and it surpassed the quasi-Newton method [17]. However, it provided the electrical parameters with large errors. A differential evolution algorithm was utilized to enhance the parameter identification performance and it showed lower errors than the series resistance-based model [18]. Nevertheless, its improper differential scheme caused the premature convergence. Particle swarm optimization was employed to identify the uncertain parameters and it outperformed the genetic algorithm [19]. However, it required more data as an input. Artificial bee colony optimization (ABC) was applied to estimate the unknown parameters and it obtained higher precision than the particle swarm optimization [20]. Nevertheless, it suffered from the local exploitation ability. Some of other metaheuristic algorithms utilized in the literature are artificial bee swarm optimization (ABSO) [21], biogeography-based optimization with mutation strategies (BBO-M) [22], bacterial foraging algorithm (BFA) [23], bird mating optimizer (BMO) [24], grouping-based global harmony search (GGHS) [25], generalized oppositional teaching–learning-based optimization (GOTLBO) [26], improved whale optimization algorithm (IWOA) [27], pattern search (PS) [28], simulated annealing (SA) [29], self-adaptive teaching–learning-based optimization (SATLBO) [30], hybrid methods [31,32,33,34,35], etc. These metaheuristic approaches also suffer from different inherent characteristics, such as premature convergence, weak exploitation, high sensitivity to the initial population and adherence to the limited number of parameters for performance improvement [36,37,38].

1.3. Research Contributions

Grey Wolf Optimizer (GWO), which was proposed in 2014 [39], is an efficient evolutionary algorithm imitating the group hunting behavior of wolves in nature. Its success in solving the optimization problems has already been proven for several fields, such as unit commitment [40], load frequency control [41], optimal power flow [42], etc. A number of variants of the basic GWO algorithm, such as mGWO [43], EGWO [44], AgGWO [45], etc. have also been proposed in order to overcome its various weaknesses in the optimization process. Grey wolf optimizer with dimension learning-based hunting search strategy (I-GWO) [46] is one of the latest GWO variants benefit from the more population diversity, the high local optima avoidance and the good balance between exploration and exploitation. On the basis of the current knowledge in the literature, the I-GWO algorithm has not yet been extensively implemented and evaluated for the parameter extraction problem of photovoltaic cells and modules. To this end, the main goals of this paper are given below. In addition, the research framework of this study is also summarized in Figure 1.
  • I-GWO is applied to all of benchmark test models including single-diode, double-diode and three-diode models in the literature;
  • The effectiveness of I-GWO is validated in terms of multiple performance aspects including accuracy, robustness and solution quality;
  • A deeper comparison is made considering not only commonly used but also recently proposed parameter extraction models in the literature;
  • The results demonstrate that I-GWO is often superior and quite competitive for reliably estimating the internal parameters of photovoltaic cells and modules.

1.4. Structure of the Article

The rest of this paper is organized as follows: Section 2 describes the diode circuit models utilized for photovoltaic systems. Section 3 presents the mathematical background of GWO and I-GWO in detail. Section 4 and Section 5 elaborate and compare the experimental results in the parameter extraction of photovoltaic cells and modules. Section 6 discusses the overall results. Finally, Section 7 concludes this paper.

2. Diode Circuit Models of Photovoltaic Systems

In the related literature, single-diode model and double-diode model are generally used for the mathematical modeling of electrical characteristics of PV cells and modules [47]. The single-diode model has easy-to-use characteristics. However, it shows unsatisfactory performance under low irradiance and neglects the influence of recombination at junction at low voltages [27,48]. Despite that, the double-diode model is preferred at low solar radiation and it considers the influence of recombination current loss in the depletion region [49]. As an alternative to the double-diode model, the three-diode model has recently been developed for modeling the effect of grain boundaries and large leakage current [50,51]. Nevertheless, it brings large computational load owing to its high number of uncertain parameters.
In the single-diode model, the photovoltaic cell is designed as a current source connected in parallel with a rectifying diode and a shunt resistance [52,53]. Figure 2 illustrates the equivalent circuit of the single-diode model for a photovoltaic cell. The output current of the single-diode model is computed using Equations (1) and (2). In these equations, I s h , I d , I p h and I L are the shunt resistor, diode, photo-generated and output current, respectively. R s h and R s denote the shunt and series resistance, respectively. I s d represents the saturation current, V L represents the output voltage and n is the diode ideality factor. As observed from Equation (2), n , R s h , R s , I s d and I p h are the five parameters required to be extracted for the single-diode model.
In the double-diode model, the photovoltaic cell is designed as a current source connected in parallel with a diffusion diode, a recombination diode and a shunt resistance [52,53]. Figure 3 depicts the equivalent circuit of the double-diode model for a photovoltaic cell. The output current of the double-diode model is calculated utilizing Equations (3) and (4). In these equations, I d 2 and I d 1 are the second diode and first diode current, respectively. n 2 and n 1 denote the ideality factors of recombination and diffusion diodes, respectively. I s d 2 and I s d 1 represent the saturation and diffusion current, respectively. As seen from Equation (4), n 2 , n 1 , R s h , R s , I s d 2 , I s d 1 and I p h are the seven parameters needed to be identified for the double-diode model.
In the three-diode model, a third diode is added in parallel to the two diodes in the double-diode model [54,55]. Figure 4 shows the equivalent circuit of the three-diode model for a photovoltaic cell. The output current of the three-diode model is formulated as in Equations (5) and (6). In these equations, I d 3 , I d 2 and I d 1 are the third diode, second diode and first diode current, respectively. Except from n 3 , same parameters affect the value of I s d 3 . As seen from Equation (6), n 3 , n 2 , n 1 , R s h , R s , I s d 3 , I s d 2 , I s d 1 and I p h are the nine parameters to be estimated for the three-diode model.
In all of diode circuit models, T is the cell temperature ( K ), k is the Boltzmann constant with the value of 1.3806503 × 10 23   J / K and q is the electron charge with the value of 1.60217646 × 10 19   C . In addition, the material resistivity at contacts is also modeled by means of the series resistance in all of diode circuit models.
I L = I p h I d I s h  
I L = I p h I s d [ e x p ( q ( V L + R s . I L ) n . k . T ) 1 ] V L + R s . I L R s h
I L = I p h I d 1 I d 2 I s h  
I L = I p h I s d 1 [ e x p ( q ( V L + R s . I L ) n 1 . k . T ) 1 ] I s d 2 [ e x p ( q ( V L + R s . I L ) n 2 . k . T ) 1 ] V L + R s . I L R s h
I L = I p h I d 1 I d 2 I d 3 I s h  
I L = I p h I s d 1 [ e x p ( q ( V L + R s . I L ) n 1 . k . T ) 1 ] I s d 2 [ e x p ( q ( V L + R s . I L ) n 2 . k . T ) 1 ] I s d 3 [ e x p ( q ( V L + R s . I L ) n 3 . k . T ) 1 ] V L + R s . I L R s h

3. Methodology

3.1. Fundamental Concepts of Grey Wolf Optimizer

Grey wolf optimizer, which was developed by Mirjalili et al. [39], is one of the widely employed bio-inspired and population-based metaheuristic optimization techniques. It mimics the behavior of the grey wolf pack in searching, encircling and attacking the prey. There is a strict social hierarchy in the grey wolf pack with different levels of leadership and dominance. Within this hierarchy, the alpha wolves, which are the leaders, make decisions on waking up, sleeping, hunting, etc. The beta wolves, which are the subordinate wolves, help the alpha wolves in the pack activities. The delta wolves obey the alpha and beta wolves, while dominating the omega wolves. Finally, the omega wolves are the scapegoats at the bottom of the pack. According to [39], the mathematical background of GWO is summarized below.
In the process of group hunting, grey wolves encircle their prey by Equation (7), where k is the current iteration, Y ( k + 1 ) represents the next location of the grey wolf, Y ( k ) indicates the current location of the grey wolf and Y p ( k ) indicates the location of the prey.
Y ( k + 1 ) = Y p ( k ) X . | Z . Y p ( k ) Y ( k ) |
X and Z are the coefficient vectors computed with Equations (8) and (9), where r 1 and r 2 are the random vectors ranging between 0 and 1. In order to emphasize the exploration and exploitation tendencies, over the iterations, x is linearly reduced and Z takes the stochastic weights.
X = 2 . x . r 1 x
Z = 2 .   r 2
The location of the prey is unspecified in the search space. For this reason, alpha, beta and delta wolves are set as the first, second and third closest search agents to the optimal solution, respectively. Omega wolves update their locations according to these best candidate solutions. This group hunting behavior is expressed with Equations (10)–(13), where Y d e l t a , Y b e t a and Y α l p h a denote the locations of delta, beta and alpha wolves, respectively. | X | > 1 and | Z | > 1 enforce the grey wolves to go away from the current prey and to search a better one, whereas | X | < 1 and | Z | < 1 enforce the grey wolves to converge towards the prey and to attack it.
Y i 1 ( k ) = Y α l p h a ( k ) X i 1 . | Z 1 . Y α l p h a Y ( k ) |
Y i 2 ( k ) = Y b e t a ( k ) X i 2 . | Z 2 . Y b e t a Y ( k ) |
Y i 3 ( k ) = Y d e l t a ( k ) X i 3 . | Z 3 . Y d e l t a Y ( k ) |
Y i G W O ( k + 1 ) = Y i 1 ( k ) + Y i 2 ( k ) + Y i 3 ( k ) 3

3.2. Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy

GWO suffers from the premature convergence, the imbalance between exploration and exploitation phases and the lack of population diversity. In order to tackle these deficiencies, Nadimi-Shahraki et al. [46] incorporated the dimension learning-based hunting search strategy to GWO, which mainly consists of initialization, movement, and selecting and updating steps. Thus, the I-GWO algorithm benefits from not only the group hunting behavior of grey wolves but also the individual hunting behavior of each wolf. According to [46], the mathematical background of I-GWO is summarized below.
In the initialization step, N wolves are stochastically distributed in the search space in the range of [ l b j , u b j ] with Equation (14), where D represents the dimension number.
Y i j = l b j + r a n d j [ 0 , 1 ] × ( u b j l b j ) ,   i [ 1 , N ] ,   j [ 1 , D ]
In the movement step, in addition to Y i G W O ( k + 1 ) , calculated by Equations (10)–(13), Y i D L H ( k + 1 ) is generated as another candidate for the new position of Y i ( k ) . To do this, initially, a radius R i ( k ) is computed utilizing Euclidean distance among Y i ( k ) and Y i G W O ( k + 1 ) with Equation (15). Following from this, N i ( k ) , which indicates the neighbors of Y i ( k ) , is constituted with Equation (16) and multi neighbors learning is realized with Equation (17). In Equation (16), d i s t i ( Y i ( k ) , Y j ( k ) ) denotes the Euclidean distance among Y i ( k ) and Y j ( k ) . In Equation (17), Y n , d ( k ) is the d-th dimension of a stochastic neighbor chosen from N i ( k ) and Y s , d ( k ) is the d-th dimension of a stochastic wolf chosen from the entire population.
R i ( k ) = Y i ( k ) Y i G W O ( k + 1 )
N i ( k ) = { Y j ( k ) | d i s t i ( Y i ( k ) , Y j ( k ) ) R i ( k ) }
Y i D L H , d ( k + 1 ) = Y i , d ( k ) + r a n d × ( Y n , d ( k ) Y s , d ( k ) )
In the selecting and updating step, initially, the best candidate for Y i ( k + 1 ) is chosen by evaluating the fitness values of Y i G W O ( k + 1 ) and Y i D L H ( k + 1 ) with Equation (18). Following, Y i ( k ) is updated by the chosen candidate in case the fitness value of the chosen candidate is smaller than Y i ( k ) . Otherwise, Y i ( k ) remains unchanged in the population. The optimization process is completed when reaching the maximum number of iterations.
Y i ( k + 1 ) = {   Y i G W O ( k + 1 ) , Y i D L H ( k + 1 ) , f ( Y i G W O ) < f ( Y i D L H ) o t h e r w i s e

4. Experimental Results on the PV Cells

In this section, the I-GWO algorithm is implemented to the three benchmark test models including single-diode, double-diode and three-diode models in order to validate its performance on the parameter extraction problem of PV cell models. The experimental IV data are obtained from [56]. The corresponding commercial silicon solar cell (R.T.C. France) has a diameter of 57 mm and operates under standard test conditions of 1000 W/m2 at 33 °C. The search ranges for the unknown internal parameters of the mentioned benchmark test models are listed in Table 1.
In the experiments, the number of search agents and maximum iterations are set to 15 and 25,000, respectively. A number of 50 independent runs are conducted to eliminate the contingency. All experiments are executed on a 2.71 GHz Intel(R) Core(TM) personal computer with 12 GB RAM under MATLAB 2016a. In order to make a fair comparison, I-GWO is compared with some commonly used parameter extraction methods in the literature. These are ABC [20], ABSO [21], BBO-M [22], BFA [23], BMO [24], GGHS [25], GOTLBO [26], IWOA [27], PS [28], SA [29] and SATLBO [30]. It should be noted that the involved methods use the same experimental IV data along with the same PV cell models.
The main objective for the parameter extraction of PV cells is to find a set of parameter values to minimize the errors between the current I L c a l c u l a t e d (fitted from the mathematical PV cell model) and the current I L m e a s u r e d (observed from the physical PV cell) [57]. Generally, the root mean squared error R M S E in Equation (19) is used as the objective function in the literature [50,58]. The lower RMSE, the more accurate parameter values. In this equation, N is the number of experimental data and x denotes the set of the extracted parameters. x = { I p h , I s d , R s , R s h , n } for the single-diode model, x = { I p h , I s d 1 , I s d 2 , R s , R s h , n 1 , n 2 } for the double-diode model and x = { I p h , I s d 1 ,   I s d 2 , I s d 3 , R s , R s h , n 1 , n 2 , n 3 } for the three-diode model. The sum of individual absolute errors ( I A E ) , expressed in Equation (20), is also utilized for evaluating the goodness-of-fit results [57].
min F ( x ) = R M S E ( x ) = 1 N k = 1 N ( I L c a l c u l a t e d k ( x ) I L m e a s u r e d k ) 2
S u m   o f   I A E = k = 1 N | I L c a l c u l a t e d k ( x ) I L m e a s u r e d k |

4.1. Experimental Results on the Single-Diode Model of the PV Cell

The values of current and power calculated by the I-GWO algorithm for the single-diode model of the PV cell are presented in Table 2. The individual absolute errors are also included in this table. According to these results, it is clear that the current and power data generated by I-GWO are highly coinciding with the measured data. The sums of IAEs are found as 0.02152728 A and 0.00873095 W for IL-calculated and PL-calculated, respectively. In case of examining the goodness-of-fit results and the extracted parameters presented in Table 3 and Table 4, respectively, it can be easily inferred that I-GWO achieves the least value of RMSE with 9.8602 × 10−4 A. Its design coefficients are obtained as 0.76077561 A for Iph, 0.32302197 μA for Isd, 0.03637706 Ω for Rs, 53.71770917 Ω for Rsh and 1.48118398 for n. In addition, I-GWO is, respectively, followed by BMO, BBO-M, GOTLBO, SATLBO, GGHS, ABSO, IWOA, ABC, SA and PS on the basis of achieving lower RMSE. On the other hand, BFA provides the greatest value of RMSE with 2.1887 × 10−1 A.
The IV and PV characteristic curves based on the extracted parameters of I-GWO are illustrated in Figure 5. It can be concluded from this figure that the I-GWO algorithm has the high capability for representing the actual photovoltaic cell characteristics for the single-diode model of the PV cell.

4.2. Experimental Results on the Double-Diode Model of the PV Cell

The values of current and power calculated by the I-GWO algorithm for the double-diode model of the PV cell are listed in Table 5. The individual absolute errors are also given in this table. With respect to these results, it is obvious that the current and power data generated by I-GWO are in very good agreement with the measured data. The sums of IAEs are found as 0.02127500 A and 0.00877613 W for IL-calculated and PL-calculated, respectively. In the case of investigating the goodness-of-fit results and the extracted parameters listed in Table 6 and Table 7, respectively, it can be simply revealed that I-GWO accomplishes the smallest RMSE value with 9.824852 × 10−4 A. Its design coefficients are acquired as 0.76078188 A for Iph, 0.22628489 μA for Isd1, 0.74609152 μA for Isd2, 0.03673977 Ω for Rs, 55.46161769 Ω for Rsh, 1.45112760 for n1 and 1.99999856 for n2. In any case, I-GWO is sequentially followed by BMO, BBO-M, SATLBO, GOTLBO, ABSO, IWOA, GGHS, ABC, PS and SA in terms of accomplishing lower RMSE. However, BFA produces the biggest RMSE value with 2.9827 × 10−1 A.
The IV and PV characteristic curves based on the extracted parameters of I-GWO are depicted in Figure 6. It can be deduced from this figure that the I-GWO algorithm has the strong ability to represent the real photovoltaic cell characteristics for the double-diode model of the PV cell.

4.3. Experimental Results on the Three-Diode Model of the PV Cell

The values of current and power calculated by the I-GWO algorithm for the three-diode model of the PV cell are tabulated in Table 8. The individual absolute errors are also provided in this table. In light of these results, it is prominent that the current and power data generated by I-GWO shows the accurate representation of the measured data. The sums of IAEs are found as 0.02128348 A and 0.00877126 W for IL-calculated and PL-calculated, respectively. In addition, the RMSE value of I-GWO is observed as 9.8251 × 10−4 A. Its design coefficients are computed as 0.76077859 A for Iph, 0.23252760 μA for Isd1, 0.15049885 μA for Isd2, 0.54357543 μA for Isd3, 0.03670937 Ω for Rs, 55.38534211 Ω for Rsh, 1.45341362 for n1, 1.99896779 for n2 and 1.99998944 for n3. It should be noted that the three-diode model of the PV cell has recently been proposed as an alternative to single-diode and double-diode models, particularly for the modeling of solar cells in industrial applications [54]. For this reason, there have not been enough studies in the literature in order to make the benchmark test for the three-diode-based modeling of the commercial silicon solar cell of R.T.C. France, yet.
The IV and PV characteristic curves based on the extracted parameters of I-GWO are visualized in Figure 7. It can be mined from this figure that the I-GWO algorithm has the great potential in reflecting the physical photovoltaic cell behavior for the three-diode model of the PV cell.

4.4. Discussion on the Single-Diode-, Double-Diode- and Three-Diode-Based Modeling Results of PV Cells

In the case of comparing the entire goodness-of-fit results achieved by the I-GWO algorithm on the R.T.C. France solar cell, the double-diode model outperforms the single-diode and three-diode models with an RMSE value of 9.824852 × 10−4 A. Following from this, the three-diode model is found to be better than the single-diode model with the RMSE value of 9.8251 × 10−4 A. Lastly, the RMSE value of the single-diode model is obtained as 9.8602 × 10−4 A. Moreover, Figure 8 shows the convergence accuracy of I-GWO for single-diode, double-diode and three-diode models of PV cells. According to this figure, it is observed during the iteration process that the relative stability of the objective function occurs around 7500 iterations in all of the diode models. It is also obvious that I-GWO is capable of running away from the local optima and accelerating the global searching direction. Especially, I-GWO surpasses ABC, ABSO, BBO-M, BFA, BMO, GGHS, GOTLBO, IWOA, PS, SA and SATLBO by achieving the least RMSE value in both single-diode and double-diode models. Despite the parametric values identified by the mentioned models seems to be very close each other, a small difference can lead a considerable effect on the PV model. As a result, the above comparisons indicate that the I-GWO algorithm is highly competitive in the estimation of the internal parameters of photovoltaic cells.

5. Experimental Results on the PV Modules

In this section, the I-GWO algorithm is applied to the two benchmark test models involving single-diode- and double-diode-based models in order to check its performance on the parameter extraction problem of PV module models. The experimental IV data were taken from [56]. The relevant solar module (Photowatt-PWP201) consists of 36 polycrystalline silicon cells in series operating under 1000 W/m2 at 45 °C. The search ranges for the unknown internal parameters of the mentioned benchmark test models are tabulated in Table 9.
The parameter settings of I-GWO and the technical specifications of PC are the same as those used in the PV cell models. For a reasonable comparison, I-GWO is checked against some widely-reported algorithms in the literature. These are artificial bee colony-differential evolution (ABC-DE) [59], chaos particle swarm optimization (CPSO) [60], flower pollination algorithm (FPA) [61], hybrid firefly and pattern search (HFAPS) [32], improved teaching–learning-based optimization (ITLBO) [3], IWOA [27], multiple learning backtracking search algorithm (LBSA) [62], PS [28], SA [29] and teaching–learning-based optimization (TLBO) [63]. The objective function in Equation (19) and the goodness-of-fit metric in Equation (20) are employed again. Only the output current equations expressed in Equations (21) and (22) [24,64] are considered for the single-diode- and double-diode-based models of PV modules, respectively. According to these equations, the IV characteristics of single-diode- and double-diode-based models of PV modules make use of the number of parallel PV cell strings ( N p ), each having the same number of series PV cells string ( N s ).
I L = I p h . N p I s d . N p [ e x p ( q ( V L N s + R s . I L N p ) n . k . T ) 1 ] V L . N p N s + R s . I L R s h
I L = I p h . N p I s d 1 . N p [ e x p ( q ( V L N s + R s . I L N p ) n 1 . k . T ) 1 ] I s d 2 . N p [ e x p ( q ( V L N s + R s . I L N p ) n 2 . k . T ) 1 ] V L . N p N s + R s . I L R s h

5.1. Experimental Results on the Single-Diode-Based Model of the PV Module

The values of current and power calculated by the I-GWO algorithm for the single-diode-based model of the PV module are presented in Table 10. The individual absolute errors are also included in this table. According to these results, it is clear that the current and power data generated by I-GWO coincide to a great extent with the measured data. The sums of IAEs are found as 0.04892353 A and 0.51689196 W for IL-calculated and PL-calculated, respectively. In the case of examining the goodness-of-fit results and the extracted parameters presented in Table 11 and Table 12, respectively, it can be easily inferred that I-GWO achieves the least value of RMSE with 2.425075 × 10−3 A. Its design coefficients are obtained as 1.03051453 A for Iph, 3.48217802 μA for Isd, 1.20127379 Ω for Rs, 981.95296539 Ω for Rsh and 48.64274143 for n. In addition, I-GWO is, respectively, followed by HFAPS, ITLBO, IWOA, LBSA, FPA, ABC-DE, SA, CPSO, PS and TLBO on the basis of achieving lower RMSE. On the other hand, TLBO provides the most value of RMSE with 6.567087 × 10−3 A.
The IV and PV characteristic curves based on the extracted parameters of I-GWO are illustrated in Figure 9. It can be concluded from this figure that the I-GWO algorithm has the high capability of representing the actual photovoltaic module characteristics for the single-diode-based model of the PV module.

5.2. Experimental Results on the Double-Diode-Based Model of the PV Module

The values of current and power calculated by the I-GWO algorithm for the double-diode-based model of the PV module are listed in Table 13. The individual absolute errors are also given in this table. With respect to these results, it is obvious that the current and power data generated by I-GWO are in very good agreement with the measured data. The sums of IAEs are found as 0.04768683 A and 0.49347099 W for IL-calculated and PL-calculated, respectively. In addition, the RMSE value of I-GWO is observed as 2.356117 × 10−3 A. Its design coefficients are acquired as 1.03046641 A for Iph, 3.59574773 × 10−10 μA for Isd1, 3.09030003 μA for Isd2, 1.27844165 Ω for Rs, 1002.52197130 Ω for Rsh, 18.96309472 for n1 and 48.32865669 for n2. It should be noted that the double-diode-based modeling of the PV module has been conducted as an alternative to the single-diode-based modeling in recent years, and it has limited applications in the corresponding literature [65]. For this reason, there have not yet been enough studies in order to make the benchmark test for the double-diode-based modeling of the polycrystalline silicon solar module of Photowatt-PWP201.
The IV and PV characteristic curves based on the extracted parameters of I-GWO are depicted in Figure 10. It can be deduced from this figure that the I-GWO algorithm has the strong ability to represent the real photovoltaic module characteristics for the double-diode-based model of the PV module.

5.3. Discussion on the Single-Diode- and Double-Diode-Based Modeling Results of PV Modules

In the case of comparing the entire goodness-of-fit results accomplished by the I-GWO algorithm on the Photowatt-PWP201 solar module, the double-diode-based model surpasses the single-diode-based model with the RMSE value of 2.356117 × 10−3 A. The RMSE value of the single-diode-based model is computed as 2.425075 × 10−3 A. Furthermore, Figure 11 demonstrates the convergence accuracy of I-GWO for single-diode- and double-diode-based models of PV modules. With regard to this figure, it is clear that I-GWO has the capability of breaking away from the local minima and providing an encouraging searching direction towards the global optima. In consequence, the comparison results show that the I-GWO algorithm seems to be very efficient for solving the parameter estimation problem of photovoltaic modules.

6. Discussion of the Overall Results

As a deeper accuracy analysis, in Table 14, the grey wolf optimizer with dimension learning-based hunting search strategy is compared with other algorithms, which were published in the last five years.
For the single-diode model of R.T.C. France PV cell, I-GWO achieves the best RMSE (9.8602 × 10−4) similar to CSO, CWOA, DE-WOA, EHHO, ISCE, ORcr-IJADE, PGJAYA, SGDE and TLABC. It also shows better RMSE than CARO, GWO, GWOCS, IJAYA, MABC, MSSO and NM-PSO.
For the double-diode model of R.T.C. France PV cell, DE-WOA and ISCE obtain the best RMSE (9.824849 × 10−4). I-GWO provides the second best RMSE (9.824852 × 10−4), which is better than CARO, CSO, CWOA, EHHO, GWO, GWOCS, IJAYA, MABC, MSSO, NM-PSO, ORcr-IJADE, PGJAYA, SGDE and TLABC.
For the single-diode-based model of Photowatt-PWP201 PV module, I-GWO accomplishes the best RMSE (2.425075 × 10−3) similar to DE-WOA, ISCE, ORcr-IJADE, PGJAYA, SGDE and TLABC. It also demonstrates better RMSE than CARO, EHHO, GWO, GWOCS, IJAYA and NM-PSO.
According to the aforementioned comparisons, it can be observed that I-GWO can achieve the best solution performance or slightly different results from the best solution performance. In addition, it can be also revealed that I-GWO can provide better solution performance than most of the implemented algorithms.

7. Conclusions

In this paper, the grey wolf optimizer with a dimension learning-based hunting search strategy is exhaustively applied to extract the internal parameters of the diode circuit models of photovoltaic cells and modules. As a result of the entire analyses, the following reasonable and useful patterns are uncovered.
According to the experimental results on the PV cells:
  • The single-diode model- and the double-diode model-based I-GWO algorithms achieve better goodness-of-fit statistics than ABC, ABSO, BBO-M, BFA, BMO, GGHS, GOTLBO, IWOA, PS, SA and SATLBO methods. The double-diode model-based I-GWO algorithm outperforms the single-diode model-based one.
  • Although the single-diode and the double-diode circuit models are widely preferred for PV cells in the literature, the three-diode model-based I-GWO algorithm shows closer solution accuracy to the double-diode model-based one.
  • In the case of only considering the recent publications on this issue, the I-GWO algorithm ensures the best RMSE, similar to most of the implemented methods for the single-diode circuit model. However, it provides the second best RMSE for the double-diode circuit model.
  • According to the experimental results on the PV modules:
  • The single-diode model-based I-GWO algorithm accomplishes better goodness-of-fit statistics than ABC-DE, CPSO, FPA, HFAPS, ITLBO, IWOA, PS, SA and TLBO methods.
  • Despite the single-diode circuit model is widely employed for PV modules in the literature, the double-diode model-based I-GWO algorithm demonstrates better solution accuracy than the single-diode model-based one.
  • When only taking into account the recent publications on this problem, the I-GWO algorithm provides the best RMSE, similar to most of the applied methods for the single-diode circuit model.
As a result, the I-GWO algorithm shows the considerable effectiveness in terms of accuracy, robustness and solution quality aspects in solving the parameter extraction problem of PV cells and modules. Hence, it can be regarded as a competitive method in order to obtain the IV and PV curves strongly fitted to the experimental measurements. In future works, the three-diode circuit model in PV cell modeling and the double-diode circuit model in PV module modeling should get more attention in order to improve parameter extraction performance. I-GWO can be modified to accomplish the online parameter estimation for more efficient photovoltaic systems. In addition, the applications of I-GWO to other engineering problems, such as fault diagnosis, photovoltaic power prediction, the power curve modeling of wind turbines, the energy consumption modeling of large-scale power systems, etc. are also interesting for future research.

Funding

No source of funding was attained for this research activity.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

ABCArtificial bee colony optimization
ABC-DEArtificial bee colony-differential evolution
ABSOArtificial bee swarm optimization
BBO-MBiogeography-based optimization with mutation strategies
BFABacterial foraging algorithm
BMOBird mating optimizer
CPSChaos particle swarm optimization
FPAFlower pollination algorithm
GGHSGrouping-based global harmony search
GOTLBOGeneralized oppositional teaching–learning-based optimization
GWOGrey wolf optimizer
HFAPS Hybrid firefly and pattern search
IAEIndividual absolute error
I-GWOGrey wolf optimizer with dimension learning-based hunting search strategy
ITLBOImproved teaching–learning-based optimization
IWOA Improved whale optimization algorithm
LBSALearning backtracking search algorithm
PSPattern search
RMSERoot mean squared error
SASimulated annealing
SATLBOSelf-adaptive teaching–learning-based optimization
TLBOTeaching–learning-based optimization

References

  1. Renewables 2020 Global Status Report. Available online: https://www.ren21.net/gsr-2020 (accessed on 1 August 2021).
  2. Humada, A.M.; Hojabri, M.; Mekhilef, S.; Hamada, H.M. Solar cell parameters extraction based on single and double-diode models: A review. Renew. Sustain. Energy Rev. 2016, 56, 494–509. [Google Scholar] [CrossRef] [Green Version]
  3. Li, S.; Gong, W.; Yan, X.; Hu, C.; Bai, D.; Wang, L.; Gao, L. Parameter extraction of photovoltaic models using an improved teaching-learning-based optimization. Energy Convers. Manag. 2019, 186, 293–305. [Google Scholar] [CrossRef]
  4. Abbassi, R.; Abbassi, A.; Jemli, M.; Chebbi, S. Identification of unknown parameters of solar cell models: A comprehensive overview of available approaches. Renew. Sustain. Energy Rev. 2018, 90, 453–474. [Google Scholar] [CrossRef]
  5. Wang, J.; Yang, B.; Li, D.; Zeng, C.; Chen, Y.; Guo, Z.; Zhang, X.; Tan, T.; Shu, H.; Yu, T. Photovoltaic cell parameter estimation based on improved equilibrium optimizer algorithm. Energy Convers. Manag. 2021, 236, 114051. [Google Scholar] [CrossRef]
  6. Merchaoui, M.; Sakly, A.; Mimouni, M.F. Particle swarm optimisation with adaptive mutation strategy for photovoltaic solar cell/module parameter extraction. Energy Convers. Manag. 2018, 175, 151–163. [Google Scholar] [CrossRef]
  7. Fathy, A.; Rezk, H. Parameter estimation of photovoltaic system using imperialist competitive algorithm. Renew. Energy 2017, 111, 307–320. [Google Scholar] [CrossRef]
  8. Blaifi, S.; Moulahoum, S.; Taghezouit, B.; Saim, A. An enhanced dynamic modeling of PV module using Levenberg-Marquardt algorithm. Renew. Energy 2019, 135, 745–760. [Google Scholar] [CrossRef]
  9. Gao, X.; Cui, Y.; Hu, J.; Tahir, N.; Xu, G. Performance comparison of exponential, Lambert W function and Special Trans function based single diode solar cell models. Energy Convers. Manag. 2018, 171, 1822–1842. [Google Scholar] [CrossRef]
  10. Et-Torabi, K.; Nassar-Eddine, I.; Obbadi, A.; Errami, Y.; Rmaily, R.; Sahnoun, S.; El-Fajri, A.; Agunaou, M. Parameters estimation of the single and double diode photovoltaic models using a Gauss-Seidel algorithm and analytical method: A comparative study. Energy Convers. Manag. 2017, 148, 1041–1054. [Google Scholar] [CrossRef]
  11. Elbaset, A.A.; Ali, H.; El Sattar, M.A. New seven parameters model for amorphous silicon and thin film PV modules based on solar irradiance. Sol. Energy 2016, 138, 26–35. [Google Scholar] [CrossRef]
  12. Gao, X.; Cui, Y.; Hu, J.; Xu, G.; Yu, Y. Lambert W-function based exact representation for double diode model of solar cells: Comparison on fitness and parameter extraction. Energy Convers. Manag. 2016, 127, 443–460. [Google Scholar] [CrossRef]
  13. Lim, L.H.I.; Ye, Z.; Ye, J.; Yang, D.; Du, H. A linear identification of diode models from single I-V characteristics of PV panels. IEEE Trans. Ind. Electron. 2015, 62, 4181–4193. [Google Scholar] [CrossRef] [Green Version]
  14. El Achouby, H.; Zaimi, M.; Ibral, A.; Assaid, E.M. New analytical approach for modelling effects of temperature and irradiance on physical parameters of photovoltaic solar module. Energy Convers. Manag. 2018, 177, 258–271. [Google Scholar] [CrossRef]
  15. Chaibi, Y.; Allouhi, A.; Salhi, M. A simple iterative method to determine the electrical parameters of photovoltaic cell. J. Clean. Prod. 2020, 269, 122363. [Google Scholar] [CrossRef]
  16. Raj, S.; Sinha, A.K.; Panchal, A.K. Solar cell parameters estimation from illuminated I-V characteristic using linear slope equations and Newton-Raphson technique. J. Renew. Sustain. Energy 2013, 5, 033105. [Google Scholar] [CrossRef]
  17. Jervase, J.A.; Bourdoucen, H.; Al-Lawati, A. Solar cell parameter extraction using genetic algorithms. Meas. Sci. Technol. 2001, 12, 1922–1925. [Google Scholar] [CrossRef] [Green Version]
  18. Ishaque, K.; Salam, Z. An improved modeling method to determine the model parameters of photovoltaic (PV) modules using differential evolution (DE). Sol. Energy 2011, 85, 2349–2359. [Google Scholar] [CrossRef]
  19. Ye, M.; Wang, X.; Xu, Y. Parameter extraction of solar cells using particle swarm optimization. J. Appl. Phys. 2009, 105, 094502. [Google Scholar] [CrossRef]
  20. Oliva, D.; Cuevas, E.; Pajares, G. Parameter identification of solar cells using artificial bee colony optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  21. Askarzadeh, A.; Rezazadeh, A. Artificial bee swarm optimization algorithm for parameters identification of solar cell models. Appl. Energy 2013, 102, 943–949. [Google Scholar] [CrossRef]
  22. Niu, Q.; Zhang, L.; Li, K. A biogeography-based optimization algorithm with mutation strategies for model parameter estimation of solar and fuel cells. Energy Convers. Manag. 2014, 86, 1173–1185. [Google Scholar] [CrossRef]
  23. Rajasekar, N.; Kumar, N.K.; Venugopalan, R. Bacterial foraging algorithm based solar PV parameter estimation. Sol. Energy 2013, 97, 255–265. [Google Scholar] [CrossRef]
  24. Askarzadeh, A.; Rezazadeh, A. Extraction of maximum power point in solar cells using bird mating optimizer-based parameters identification approach. Sol. Energy 2013, 90, 123–133. [Google Scholar] [CrossRef]
  25. Askarzadeh, A.; Rezazadeh, A. Parameter identification for solar cell models using harmony search-based algorithms. Sol. Energy 2012, 86, 3241–3249. [Google Scholar] [CrossRef]
  26. Chen, X.; Yu, K.; Du, W.; Zhao, W.; Liu, G. Parameters identification of solar cell models using generalized oppositional teaching learning based optimization. Energy 2016, 99, 170–180. [Google Scholar] [CrossRef]
  27. Xiong, G.; Zhang, J.; Shi, D.; He, Y. Parameter extraction of solar photovoltaic models using an improved whale optimization algorithm. Energy Convers. Manag. 2018, 174, 388–405. [Google Scholar] [CrossRef]
  28. AlHajri, M.F.; El-Naggar, K.M.; AlRashidi, M.R.; Al-Othman, A.K. Optimal extraction of solar cell parameters using pattern search. Renew. Energy 2012, 44, 238–245. [Google Scholar] [CrossRef]
  29. El-Naggar, K.M.; AlRashidi, M.R.; AlHajri, M.F.; Al-Othman, A.K. Simulated annealing algorithm for photovoltaic parameters identification. Sol. Energy 2012, 86, 266–274. [Google Scholar] [CrossRef]
  30. Yu, K.; Chen, X.; Wang, X.; Wang, Z. Parameters identification of photovoltaic models using self-adaptive teaching-learning-based optimization. Energy Convers. Manag. 2017, 145, 233–246. [Google Scholar] [CrossRef]
  31. Ayang, A.; Wamkeue, R.; Ouhrouche, M.; Djongyang, N.; Salome, N.E.; Pombe, J.K.; Ekemb, G. Maximum likelihood parameters estimation of single-diode model of photovoltaic generator. Renew. Energy 2019, 130, 111–121. [Google Scholar] [CrossRef]
  32. Beigi, A.M.; Maroosi, A. Parameter identification for solar cells and module using a hybrid firefly and pattern search algorithms. Sol. Energy 2018, 171, 435–446. [Google Scholar] [CrossRef]
  33. Ram, J.P.; Babu, T.S.; Dragicevic, T.; Rajasekar, N. A new hybrid bee pollinator flower pollination algorithm for solar PV parameter estimation. Energy Convers. Manag. 2017, 135, 463–476. [Google Scholar] [CrossRef]
  34. Yu, K.; Liang, J.J.; Qu, B.Y.; 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]
  35. Yousri, D.; Allam, D.; Eteiba, M.B.; Suganthan, P.N. Static and dynamic photovoltaic models’ parameters identification using chaotic heterogeneous comprehensive learning particle swarm optimizer variants. Energy Convers. Manag. 2019, 182, 546–563. [Google Scholar] [CrossRef]
  36. 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]
  37. Cárdenas-Bravo, C.; Barraza, R.; Sánchez-Squella, A.; Valdivia-Lefort, P.; Castillo-Burns, F. Estimation of single-diode photovoltaic model using the differential evolution algorithm with adaptive boundaries. Energies 2021, 14, 3925. [Google Scholar] [CrossRef]
  38. Liu, E.J.; Hung, Y.H.; Hong, C.W. Improved metaheuristic optimization algorithm applied to hydrogen fuel cell and photovoltaic cell parameter extraction. Energies 2021, 14, 619. [Google Scholar] [CrossRef]
  39. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  40. Panwar, L.K.; Reddy, S.; Verma, A.; Panigrahi, B.K.; Kumar, R. Binary grey wolf optimizer for large scale unit commitment problem. Swarm Evol. Comput. 2018, 38, 251–266. [Google Scholar] [CrossRef]
  41. Guha, D.; Roy, P.K.; Banerjee, S. Load frequency control of interconnected power system using grey wolf optimization. Swarm Evol. Comput. 2016, 27, 97–115. [Google Scholar] [CrossRef]
  42. El-Fergany, A.A.; Hasanien, H.M. Single and multi-objective optimal power flow using grey wolf optimizer and differential evolution algorithms. Electr. Power Compon. Syst. 2015, 43, 1548–1559. [Google Scholar] [CrossRef]
  43. Mittal, N.; Singh, U.; Sohi, B.S. Modified grey wolf optimizer for global engineering optimization. Appl. Comput. Intell. Soft Comput. 2016, 2016, 7950348. [Google Scholar] [CrossRef] [Green Version]
  44. Joshi, H.; Arora, S. Enhanced grey wolf optimization algorithm for global optimization. Fundam. Inform. 2017, 153, 235–264. [Google Scholar] [CrossRef]
  45. Hu, P.; Chen, S.; Huang, H.; Zhang, G.; Liu, L. Improved alpha-guided grey wolf optimizer. IEEE Access 2019, 7, 5421–5437. [Google Scholar] [CrossRef]
  46. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S. An improved grey wolf optimizer for solving engineering problems. Expert Syst. Appl. 2020, 166, 113917. [Google Scholar] [CrossRef]
  47. Muangkote, N.; Sunat, K.; Chiewchanwattana, S.; Kaiwinit, S. An advanced onlooker-ranking-based adaptive differential evolution to extract the parameters of solar cell models. Renew. Energy 2019, 134, 1129–1147. [Google Scholar] [CrossRef]
  48. Kler, D.; Goswami, Y.; Rana, K.P.S.; Kumar, V. A novel approach to parameter estimation of photovoltaic systems using hybridized optimizer. Energy Convers. Manag. 2019, 187, 486–511. [Google Scholar] [CrossRef]
  49. Sheng, H.; Li, C.; Wang, H.; Yan, Z.; Xiong, Y.; Cao, Z.; Kuang, Q. Parameters extraction of photovoltaic models using an improved moth-flame optimization. Energies 2019, 12, 3527. [Google Scholar] [CrossRef] [Green Version]
  50. Elaziz, M.A.; Oliva, D. Parameter estimation of solar cells diode models by an improved opposition-based whale optimization algorithm. Energy Convers. Manag. 2018, 171, 1843–1859. [Google Scholar] [CrossRef]
  51. Allam, D.; Yousri, D.A.; Eteiba, M.B. Parameters extraction of the three diode model for the multi-crystalline solar cell/module using moth-flame optimization algorithm. Energy Convers. Manag. 2016, 123, 535–548. [Google Scholar] [CrossRef]
  52. Abbassi, R.; Abbassi, A.; Heidari, A.A.; Mirjalili, S. An efficient salp swarm-inspired algorithm for parameters identification of photovoltaic cell models. Energy Convers. Manag. 2019, 179, 362–372. [Google Scholar] [CrossRef]
  53. Chen, X.; Yu, K. Hybridizing cuckoo search algorithm with biogeography-based optimization for estimating photovoltaic model parameters. Sol. Energy 2019, 180, 192–206. [Google Scholar] [CrossRef]
  54. Khanna, V.; Das, B.K.; Bisht, D.; Singh, P.K. A three diode model for industrial solar cells and estimation of solar cell parameters using PSO algorithm. Renew. Energy 2015, 78, 105–113. [Google Scholar] [CrossRef]
  55. Yesilbudak, M. Accurate identification of the electrical parameters of triple-diode photovoltaic model using a metaheuristic algorithm. In Proceedings of the 9th International Conference on Smart Grid, Setubal, Portugal, 29 June–1 July 2021. [Google Scholar]
  56. 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]
  57. Gao, X.; Cui, Y.; Hu, J.; Xu, G.; Wang, Z.; Qu, J.; Wang, H. Parameter extraction of solar cell models using improved shuffled complex evolution algorithm. Energy Convers. Manag. 2018, 157, 460–479. [Google Scholar] [CrossRef]
  58. Xiong, G.; Zhang, J.; Yuan, X.; Shi, D.; He, Y.; Yao, G. Parameter extraction of solar photovoltaic models by means of a hybrid differential evolution with whale optimization algorithm. Sol. Energy 2018, 176, 742–761. [Google Scholar] [CrossRef]
  59. Hachana, O.; Hemsas, K.E.; Tina, G.M.; Ventura, C. Comparison of different metaheuristic algorithms for parameter identification of photovoltaic cell/module. J. Renew. Sustain. Energy 2013, 5, 053122. [Google Scholar] [CrossRef]
  60. Wei, H.; Cong, J.; Lingyun, X.; Deyun, S. Extracting solar cell model parameters based on chaos particle swarm algorithm. In Proceedings of the International Conference on Electric Information and Control Engineering, Wuhan, China, 15–17 April 2011; pp. 398–402. [Google Scholar] [CrossRef]
  61. Alam, D.F.; Yousri, D.A.; Eteiba, M.B. Flower pollination algorithm based solar PV parameter estimation. Energy Convers. Manag. 2015, 101, 410–422. [Google Scholar] [CrossRef]
  62. Yu, K.; Liang, J.J.; Qu, B.Y.; Cheng, Z.; Wang, H. Multiple learning backtracking search algorithm for estimating parameters of photovoltaic models. Appl. Energy 2018, 226, 408–422. [Google Scholar] [CrossRef]
  63. Patel, S.J.; Panchal, A.K.; Kheraj, V. Extraction of solar cell parameters from a single current-voltage characteristic using teaching learning based optimization algorithm. Appl. Energy 2014, 119, 384–393. [Google Scholar] [CrossRef]
  64. Jordehi, A.R. Time varying acceleration coefficients particle swarm optimisation (TVACPSO): A new optimisation algorithm for estimating parameters of PV cells and modules. Energy Convers. Manag. 2016, 129, 262–274. [Google Scholar] [CrossRef]
  65. Nunes, H.G.G.; Pombo, J.A.N.; Mariano, S.J.P.S.; Calado, M.R.A.; Felippe de Souza, J.A.M. 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]
  66. Yuan, X.; He, Y.; Liu, L. Parameter extraction of solar cell models using chaotic asexual reproduction optimization. Neural Comput. Appl. 2015, 26, 1227–1239. [Google Scholar] [CrossRef]
  67. Guo, L.; Meng, Z.; Sun, Y.; Wang, L. Parameter identification and sensitivity analysis of solar cell models with cat swarm optimization algorithm. Energy Convers. Manag. 2016, 108, 520–528. [Google Scholar] [CrossRef]
  68. Oliva, D.; Aziz, M.A.E.; Hassanien, A.E. Parameter estimation of photovoltaic cells using an improved chaotic whale optimization algorithm. Appl. Energy 2017, 200, 141–154. [Google Scholar] [CrossRef]
  69. Jiao, S.; Chong, G.; Huang, C.; Hu, H.; Wang, M.; Heidari, A.A.; Chen, H.; Zhao, X. Orthogonally adapted Harris hawks optimization for parameter estimation of photovoltaic models. Energy 2020, 203, 117804. [Google Scholar] [CrossRef]
  70. Jamadi, M.; Merrikh-Bayat, F.; Bigdeli, M. Very accurate parameter estimation of single- and double-diode solar cell models using a modified artificial bee colony algorithm. Int. J. Energy Environ. Eng. 2016, 7, 13–25. [Google Scholar] [CrossRef] [Green Version]
  71. Lin, P.; Cheng, S.; Yeh, W.; Chen, Z.; Wu, L. Parameters extraction of solar cell models using a modified simplified swarm optimization algorithm. Sol. Energy 2017, 144, 594–603. [Google Scholar] [CrossRef]
  72. Hamid, N.F.A.; Rahim, N.A.; Selvaraj, J. Solar cell parameters identification using hybrid Nelder-Mead and modified particle swarm optimization. J. Renew. Sustain. Energy 2016, 8, 015502. [Google Scholar] [CrossRef]
  73. Yu, K.; Qu, B.; Yue, C.; Ge, S.; Chen, X.; Liang, J. A performance-guided JAYA algorithm for parameters identification of photovoltaic cell and module. Appl. Energy 2019, 237, 241–257. [Google Scholar] [CrossRef]
  74. Liang, J.; Qiao, K.; Yuan, M.; Yu, K.; Qu, B.; Ge, S.; Li, Y.; Chen, G. Evolutionary multi-task optimization for parameters extraction of photovoltaic models. Energy Convers. Manag. 2020, 207, 112509. [Google Scholar] [CrossRef]
  75. Chen, X.; Xu, B.; Mei, C.; Ding, Y.; Lia, K. Teaching-learning-based artificial bee colony for solar photovoltaic parameter estimation. Appl. Energy 2018, 212, 1578–1588. [Google Scholar] [CrossRef]
Figure 1. The research framework of this study.
Figure 1. The research framework of this study.
Energies 14 05735 g001
Figure 2. Equivalent circuit of the single-diode model.
Figure 2. Equivalent circuit of the single-diode model.
Energies 14 05735 g002
Figure 3. Equivalent circuit of the double-diode model.
Figure 3. Equivalent circuit of the double-diode model.
Energies 14 05735 g003
Figure 4. Equivalent circuit of the three-diode model.
Figure 4. Equivalent circuit of the three-diode model.
Energies 14 05735 g004
Figure 5. Comparison of experimental and estimated photovoltaic cell characteristics for the single-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Figure 5. Comparison of experimental and estimated photovoltaic cell characteristics for the single-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Energies 14 05735 g005
Figure 6. Comparison of experimental and estimated photovoltaic cell characteristics for the double-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Figure 6. Comparison of experimental and estimated photovoltaic cell characteristics for the double-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Energies 14 05735 g006
Figure 7. Comparison of experimental and estimated photovoltaic cell characteristics for the three-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Figure 7. Comparison of experimental and estimated photovoltaic cell characteristics for the three-diode model of the PV cell (a) IV characteristics (b) PV characteristics.
Energies 14 05735 g007
Figure 8. The convergence curves for single-diode, double-diode and three-diode models of PV cells.
Figure 8. The convergence curves for single-diode, double-diode and three-diode models of PV cells.
Energies 14 05735 g008
Figure 9. Comparison of experimental and estimated photovoltaic module characteristics for the single-diode-based model of the PV module (a) IV characteristics (b) PV characteristics.
Figure 9. Comparison of experimental and estimated photovoltaic module characteristics for the single-diode-based model of the PV module (a) IV characteristics (b) PV characteristics.
Energies 14 05735 g009
Figure 10. Comparison of experimental and estimated photovoltaic module characteristics for the double-diode-based model of the PV module (a) IV characteristics (b) PV characteristics.
Figure 10. Comparison of experimental and estimated photovoltaic module characteristics for the double-diode-based model of the PV module (a) IV characteristics (b) PV characteristics.
Energies 14 05735 g010aEnergies 14 05735 g010b
Figure 11. The convergence curves for single-diode- and double-diode-based models of PV modules.
Figure 11. The convergence curves for single-diode- and double-diode-based models of PV modules.
Energies 14 05735 g011
Table 1. The search ranges for the parameters of photovoltaic cell models.
Table 1. The search ranges for the parameters of photovoltaic cell models.
ParameterLower BoundUpper Bound
Iph (A)01
Isd, Isd1, Isd2, Isd3 (μA)01
Rs (Ω)00.5
Rsh (Ω)0100
n, n1, n2, n312
Table 2. IAEs of I-GWO for the single-diode model of the PV cell.
Table 2. IAEs of I-GWO for the single-diode model of the PV cell.
VL-measured (V)IL-measured (A)IL-calculated (A)IAE (IL)PL-measured (W)PL-calculated (W)IAE (PL)
−0.20570.76400.764087830.00008783−0.15715480−0.157172870.00001807
−0.12910.76200.762663190.00066319−0.09837420−0.098459820.00008562
−0.05880.76050.761355390.00085539−0.04471740−0.044767700.00005030
0.00570.76050.760154060.000345940.004334850.004332880.00000197
0.06460.76000.759055260.000944740.049096000.049034970.00006103
0.11850.75900.758042380.000957620.089941500.089828020.00011348
0.16780.75700.757091680.000091680.127024600.127039980.00001538
0.21320.75700.756141370.000858630.161392400.161209340.00018306
0.25450.75550.755086870.000413130.192274750.192169610.00010514
0.29240.75400.753663860.000336140.220469600.220371310.00009829
0.32690.75050.751390940.000890940.245338450.245629700.00029125
0.35850.74650.747353820.000853820.267620250.267926340.00030609
0.38730.73850.740117180.001617180.286021050.286647380.00062633
0.41370.72800.727382180.000617820.301173600.300918010.00025559
0.43730.70650.706972600.000472600.308952450.309159120.00020667
0.45900.67550.675280110.000219890.310054500.309953570.00010093
0.47840.63200.630758250.001241750.302348800.301754740.00059406
0.49600.57300.571928350.001071650.284208000.283676460.00053154
0.51190.49900.499607040.000607040.255438100.255748850.00031075
0.52650.41300.413648850.000648850.217444500.217786120.00034162
0.53980.31650.317510200.001010200.170846700.171392000.00054530
0.55210.21200.212155050.000155050.117045200.117130800.00008560
0.56330.10350.102251430.001248570.058301550.057598230.00070332
0.5736−0.0100−0.008717430.00128257−0.00573600−0.005000320.00073568
0.5833−0.1230−0.125507320.00250732−0.07174590−0.073208420.00146252
0.5900−0.2100−0.208472260.00152774−0.12390000−0.122998630.00090137
Table 3. Comparison of the goodness-of-fit results achieved by I-GWO for the single-diode model of the PV cell.
Table 3. Comparison of the goodness-of-fit results achieved by I-GWO for the single-diode model of the PV cell.
IL-measured
(A)
IL-calculated (A)
I-GWOABCABSOBBO-MBFABMOGGHSGOTLBOIWOAPSSASATLBO
0.76400.764087830.764116490.764199500.7640060.763756380.764073010.764273950.764069470.764111040.764536860.763020.76411170
0.76200.762663190.762689880.762735990.7626040.762251770.762652440.762831780.762655300.762686770.763342810.760290.76268727
0.76050.761355390.761380270.761392500.7613170.760870460.761348370.761507880.761357110.761379310.762246510.758460.76137966
0.76050.760154060.760177260.760158420.7601350.759601290.760150460.760291770.760164580.760178290.761238730.759530.76017850
0.76000.759055260.759076930.759029800.7590530.758439420.759054790.759179460.759073790.759079770.760314050.759510.75907986
0.75900.758042380.758062600.757989890.7580560.757366720.758044730.758154180.758068110.758067150.759452770.758420.75806713
0.75700.757091680.757110460.757015170.7571200.756359390.757096540.757192030.757123650.757116650.758620510.755280.75711655
0.75700.756141370.756158510.756044390.7561820.755362900.756148370.756230900.756178280.756166470.757732280.756160.75616634
0.75550.755086870.755101750.754974900.7551380.754304790.755095410.755166000.755126520.755111900.756636000.754180.75511189
0.75400.753663860.753675030.753545090.7537230.753014500.753673060.753732400.753702690.753688440.754991750.752620.75368886
0.75050.751390940.751395530.751277660.7514530.751231460.751399410.751448230.751423460.751414240.752215490.748000.75141578
0.74650.747353820.747346590.747262870.7474140.748505430.747359610.747398660.747372480.747374250.747271720.744250.74737831
0.73850.740117180.740089650.740068390.7401680.744219510.740117890.740148540.740113290.740132000.738636380.735840.74014108
0.72800.727382180.727321290.727394110.7274160.737462590.727375260.727400230.727348100.727386800.723991060.727950.72740519
0.70650.706972600.706862550.707053190.7069850.727555330.706956410.706980350.706906380.706961080.701340900.705950.70699455
0.67550.675280110.675100320.675421970.6752690.713264370.675254260.675284030.675187710.675244290.667287450.674940.67530100
0.63200.630758250.630489740.630932430.6307280.694353440.630724650.630767990.630656130.630690240.620661050.630230.63077852
0.57300.571928350.571554760.572093700.5718870.670506820.571890470.571953510.571838410.571821560.560127690.571520.57194903
0.49900.499607040.499119130.499723170.4995630.642210850.499568920.499653070.499549630.499458520.486438730.498450.49962975
0.41300.413648850.413042190.413683200.4136120.609470830.413614350.413715900.413639480.413458720.399220350.411890.41367578
0.31650.317510200.316790100.317456730.3174850.573516990.317481470.317588540.317549820.317284100.301429830.314590.31754398
0.21200.212155050.211331290.212030150.2121420.534597200.212132570.212227260.212232640.211901860.193557740.210770.21219865
0.10350.102251430.101340780.102099240.1022450.494271110.102233410.102292670.102341410.101983770.079831750.103000.10230769
−0.0100−0.00871743−0.00969324−0.00881685−0.0087310.45354725−0.00873587−0.00874118−0.00865916−0.00898356−0.03693312−0.01094−0.00864575
−0.1230−0.12550732−0.12653296−0.12549663−0.1255370.41080836−0.12552919−0.12562452−0.12551303−0.12575903−0.16117142−0.12753−0.12541726
−0.2100−0.20847226−0.20951676−0.20829509−0.2085300.38015838−0.20850287−0.20868792−0.20857077−0.20870043−0.25151582−0.21250−0.20836716
RMSE9.8602 × 10−410.967 × 10−49.9124 × 10−49.8634 × 10−42.1887 × 10−19.8622 × 10−49.9089 × 10−49.8744 × 10−49.9487 × 10−41.4936 × 10−21.71 × 10−39.8780 × 10−4
Table 4. Comparison of the parameters extracted by I-GWO for the single-diode model of the PV cell.
Table 4. Comparison of the parameters extracted by I-GWO for the single-diode model of the PV cell.
Design CoefficientsI-GWOABCABSOBBO-MBFABMOGGHSGOTLBOIWOAPSSASATLBO
Iph (A)0.760775610.76080.760800.760780.76020.760770.760920.7607800.76080.76170.76200.7608
Isd (μA)0.323021970.32510.306230.318740.80000.324790.326200.3315520.32320.99800.47980.32315
Rs (Ω)0.036377060.03640.036590.036420.03250.036360.036310.0362650.03640.03130.03450.03638
Rsh (Ω)53.7177091753.643352.290353.3622750.869153.871653.064754.11542653.731764.102643.103553.7256
n1.481183981.48171.475831.479841.69511.481731.482171.4838201.48121.61.51721.48123
Table 5. IAEs of I-GWO for the double-diode model of the PV cell.
Table 5. IAEs of I-GWO for the double-diode model of the PV cell.
VL-measured (V)IL-measured (A)IL-calculated (A)IAE (IL)PL-measured (W)PL-calculated (W)IAE (PL)
−0.20570.76400.763985590.00001441−0.15715480−0.157151840.00000296
−0.12910.76200.762605680.00060568−0.09837420−0.098452390.00007819
−0.05880.76050.761338740.00083874−0.04471740−0.044766720.00004932
0.00570.76050.760174340.000325660.004334850.004332990.00000186
0.06460.76000.759107780.000892220.049096000.049038360.00005764
0.11850.75900.758121130.000878870.089941500.089837350.00010415
0.16780.75700.757188000.000188000.127024600.127056150.00003155
0.21320.75700.756242750.000757250.161392400.161230950.00016145
0.25450.75550.755176320.000323680.192274750.192192370.00008238
0.29240.75400.753721380.000278620.220469600.220388130.00008147
0.32690.75050.751398340.000898340.245338450.245632120.00029367
0.35850.74650.747300960.000800960.267620250.267907390.00028714
0.38730.73850.740010580.001510580.286021050.286606100.00058505
0.41370.72800.727247270.000752730.301173600.300862190.00031141
0.43730.70650.706850890.000350890.308952450.309105890.00015344
0.45900.67550.675211220.000288780.310054500.309921950.00013255
0.47840.63200.630761330.001238670.302348800.301756220.00059258
0.49600.57300.571995060.001004940.284208000.283709550.00049845
0.51190.49900.499706190.000706190.255438100.255799600.00036150
0.52650.41300.413733530.000733530.217444500.217830700.00038620
0.53980.31650.317546000.001046000.170846700.171411330.00056463
0.55210.21200.212122890.000122890.117045200.117113050.00006785
0.56330.10350.102163430.001336570.058301550.057548660.00075289
0.5736−0.0100−0.008791280.00120872−0.00573600−0.005042680.00069332
0.5833−0.1230−0.125542490.00254249−0.07174590−0.073228930.00148303
0.5900−0.2100−0.208370390.00162961−0.12390000−0.122938530.00096147
Table 6. Comparison of the goodness-of-fit results achieved by I-GWO for the double-diode model of the PV cell.
Table 6. Comparison of the goodness-of-fit results achieved by I-GWO for the double-diode model of the PV cell.
IL-measured
(A)
IL-calculated (A)
I-GWOABCABSOBBO-MBFABMOGGHSGOTLBOIWOAPSSASATLBO
0.76400.763985590.764108040.764035020.7640120.763881450.763962840.763403090.763920120.764007290.762430380.766461210.76400509
0.76200.762605680.762685060.762633930.7626220.762605950.762591490.762184250.762555310.762626060.761488910.764685660.76261657
0.76050.761338740.761378790.761347650.7613450.761435140.761332370.761065300.761302210.761357920.760624460.763055720.76134178
0.76050.760174340.760178850.760165770.7601720.760360040.760175030.760037200.760150470.760192470.759829530.761558470.76017033
0.76000.759107780.759081290.759083990.7590980.759378240.759114670.759095950.759095360.759125100.759099140.760188800.75909775
0.75900.758121130.758069480.758085010.7581060.758478870.758133080.758225560.758118930.758137990.758415560.758925930.75810646
0.75700.757188000.757119550.757143740.7571680.757652920.757203410.757401070.757194610.757205000.757745550.757740660.75717084
0.75700.756242750.756169480.756196430.7562210.756879580.756259300.756558180.756256560.756260800.757006380.756559120.75622630
0.75550.755176320.755114120.755136180.7551570.756146560.755190850.755583260.755195200.755196590.756044220.755261640.75516493
0.75400.753721380.753688200.753697150.7537080.755401130.753730180.754200830.753742870.753745040.754519740.753550640.75372013
0.75050.751398340.751408740.751398830.7513950.754566130.751398020.751912860.751419730.751425810.751845050.750900860.75141152
0.74650.747300960.747359000.747331040.7473100.753484090.747289650.747788170.747319690.747331140.746979020.746333350.74733029
0.73850.740010580.740100500.740066900.7400290.751916480.739989480.740381960.740024630.740039770.738385540.738346430.74005270
0.72800.727247270.727330020.727316990.7272700.749469830.727220990.727402920.727255240.727268870.723726660.724572010.72729318
0.70650.706850890.706868940.706912990.7068690.745774080.706826820.706718840.706852220.706856600.700990120.702834470.70688721
0.67550.675211220.675104210.675245040.6752170.740185740.675195900.674780410.675204850.675192690.666754490.669490750.67522544
0.63200.630761330.630490800.630755920.6307530.732408390.630757030.630122870.630745790.630713920.619854180.623064820.63074786
0.57300.571995060.571551990.571952910.5719760.722133150.571998800.571312290.571968690.571918510.558984570.562059740.57195709
0.49900.499706190.499110620.499641470.4996850.709440750.499711210.499162090.499668580.499604590.484961740.487254010.49965419
0.41300.413733530.413024960.413668620.4137230.694245940.413730790.413497240.413685670.413613750.397472440.398379750.41368276
0.31650.317546000.316760560.317498340.3175530.677098960.317530400.317678360.317493160.317414890.299572360.298774960.31750791
0.21200.212122890.211285240.212103200.2121510.658127340.212093490.212588270.212073640.211985950.191820790.189233440.21210409
0.10350.102163430.101274310.102168070.1022080.638137580.102126500.102795160.102130070.102023110.078507710.074388730.10216083
−0.0100−0.00879128−0.00978343−0.00879319−0.0087500.61773270−0.00881673−0.00834550−0.00879077−0.00893910−0.03747352−0.04242094−0.00879859
−0.1230−0.12554249−0.12665126−0.12556054−0.1255130.59606795−0.12554722−0.12551080−0.12549444−0.12569749−0.16060685−0.16597348−0.12556179
−0.2100−0.20837039−0.20965688−0.20845031−0.2083790.58050794−0.20833403−0.20905426−0.20827036−0.20854079−0.24982997−0.25460551−0.20843149
RMSE9.824852 × 10−411.146 × 10−49.8344 × 10−49.8272 × 10−42.9827 × 10−19.8266 × 10−410.684 × 10−49.8317 × 10−49.8580 × 10−41.5176 × 10−216.644 × 10−39.8294 × 10−4
Table 7. Comparison of the parameters extracted by I-GWO for the double-diode model of the PV cell.
Table 7. Comparison of the parameters extracted by I-GWO for the double-diode model of the PV cell.
Design CoefficientsI-GWOABCABSOBBO-MBFABMOGGHSGOTLBOIWOAPSSASATLBO
Iph (A)0.760781880.76080.760780.760830.76090.760780.760560.7607520.76080.76020.76230.76078
Isd1 (μA)0.226284890.04070.267130.591150.00940.211100.370140.8001950.67710.98890.47670.25093
Isd2 (μA)0.746091520.28740.381910.245230.04530.876880.135040.2204620.23550.00010.01000.545418
Rs (Ω)0.036739770.03640.036570.036640.03510.036820.035620.0367830.03670.03200.03450.03663
Rsh (Ω)55.4616176953.780454.621955.04946055.808162.789956.07530455.408281.300843.103555.1170
n11.451127601.44951.4651221.38091.445331.496381.99997321.61.51721.45982
n21.999998561.48851.981521.457981.52551.999971.929981.4489741.45451.192021.99941
Table 8. IAEs of I-GWO for the three-diode model of the PV cell.
Table 8. IAEs of I-GWO for the three-diode model of the PV cell.
VL-measured (V)IL-measured (A)IL-calculated (A)IAE (IL)PL-measured (W)PL-calculated (W)IAE (PL)
−0.20570.76400.763987090.00001291−0.15715480−0.157152140.00000266
−0.12910.76200.762605280.00060528−0.09837420−0.098452340.00007814
−0.05880.76050.761336620.00083662−0.04471740−0.044766590.00004919
0.00570.76050.760170670.000329330.004334850.004332970.00000188
0.06460.76000.759102810.000897190.049096000.049038040.00005796
0.11850.75900.758115200.000884800.089941500.089836650.00010485
0.16780.75700.757181640.000181640.127024600.127055080.00003048
0.21320.75700.756236780.000763220.161392400.161229680.00016272
0.25450.75550.755171860.000328140.192274750.192191240.00008351
0.29240.75400.753719740.000280260.220469600.220387650.00008195
0.32690.75050.751400630.000900630.245338450.245632870.00029442
0.35850.74650.747307630.000807630.267620250.267909790.00028954
0.38730.73850.740020790.001520790.286021050.286610050.00058900
0.41370.72800.727258750.000741250.301173600.300866950.00030665
0.43730.70650.706860300.000360300.308952450.309110010.00015756
0.45900.67550.675215570.000284430.310054500.309923940.00013056
0.47840.63200.630759470.001240530.302348800.301755330.00059347
0.49600.57300.571988160.001011840.284208000.283706130.00050187
0.51190.49900.499697170.000697170.255438100.255794980.00035688
0.52650.41300.413726460.000726460.217444500.217826980.00038248
0.53980.31650.317543650.001043650.170846700.171410060.00056336
0.55210.21200.212126570.000126570.117045200.117115080.00006988
0.56330.10350.102171680.001328320.058301550.057553300.00074825
0.5736−0.0100−0.008784860.00121514−0.00573600−0.005039000.00069700
0.5833−0.1230−0.125540520.00254052−0.07174590−0.073227780.00148188
0.5900−0.2100−0.208381150.00161885−0.12390000−0.122944880.00095512
Table 9. The search ranges for the parameters of photovoltaic module models.
Table 9. The search ranges for the parameters of photovoltaic module models.
ParameterLower BoundUpper Bound
Iph (A)02
Isd, Isd1, Isd2 (μA)050
Rs (Ω)02
Rsh (Ω)02000
n, n1, n2150
Table 10. IAEs of I-GWO for the single-diode-based model of the PV module.
Table 10. IAEs of I-GWO for the single-diode-based model of the PV module.
VL-measured (V)IL-measured (A)IL-calculated (A)IAE (IL)PL-measured (W)PL-calculated (W)IAE (PL)
0.12481.03151.029119350.002380650.128731200.128434090.00029711
1.80931.03001.027381210.002618791.863579001.858840830.00473817
3.35111.02601.025741890.000258113.438228603.437363650.00086495
4.76221.02201.024107210.002107214.866968404.877003340.01003494
6.05381.01801.022291820.004291826.162768406.188750240.02598184
7.23641.01551.019930670.004430677.348564207.380626330.03206213
8.31891.01401.016363080.002363088.435364608.455022860.01965826
9.30971.01001.010496130.000496139.402797009.407415810.00461881
10.21631.00351.000628960.0028710410.2520570510.222725640.02933141
11.04490.98800.984548390.0034516110.9123612010.874238570.03812263
11.80180.96300.959521730.0034782711.3651334011.324083520.04104988
12.49290.92550.922838900.0026611011.5621789511.528934130.03324482
13.12310.87250.872599770.0000997711.4499047511.451214050.00130930
13.69830.80750.807274380.0002256211.0613772511.058286570.00309068
14.22210.72650.728336570.0018365710.3323556510.358475560.02611991
14.69950.63450.637138060.002638069.326832759.365610890.03877814
15.13460.53450.536213070.001713078.089443708.115370390.02592669
15.53110.42750.429511290.002011296.639545256.670782800.03123755
15.89290.31850.318774410.000274415.061888655.066249830.00436118
16.22290.20850.207389410.001110593.382474653.364457730.01801692
16.52410.10100.096167080.004832921.668934101.589074400.07985970
16.7987−0.0080−0.008325440.00032544−0.13438960−0.139856630.00546703
17.0499−0.1110−0.110936490.00006351−1.89253890−1.891456000.00108290
17.2793−0.2090−0.209247200.00024720−3.61137370−3.615645060.00427136
17.4885−0.3030−0.300863420.00213658−5.29901550−5.261649860.03736564
Table 11. Comparison of the goodness-of-fit results achieved by I-GWO for the single-diode-based model of the PV module.
Table 11. Comparison of the goodness-of-fit results achieved by I-GWO for the single-diode-based model of the PV module.
IL-measured
(A)
IL-calculated (A)
I-GWOABC-DECPSOFPAHFAPSITLBOIWOALBSAPSSATLBO
1.03151.029119351.030174451.027921561.030383641.029108071.029104831.029099841.029056961.029378961.031459771.02930433
1.03001.027381211.028160191.026967971.028287401.027373981.027366741.027356021.027383161.027000301.029415041.02621400
1.02601.025741891.026270401.026014431.026324681.025738331.025727461.025711591.025802201.024778761.027493951.02334149
1.02201.024107211.024410281.024943141.024402121.024106951.024092811.024072441.024220001.022617701.025597531.02058283
1.01801.022291821.022395871.023523381.022339071.022294501.022277451.022253341.022450751.020333381.023534591.01774797
1.01551.019930671.019864761.021323831.019778511.019935861.019916281.019889481.020128401.017565421.020931361.01446920
1.01401.016363081.016160521.017603761.016071481.016370241.016348611.016320501.016590381.013663001.017116011.01009017
1.01001.010496131.010193981.011203841.010136811.010504591.010481471.010453841.010741191.007546281.010982301.00352898
1.00351.000628961.000263011.000416931.000277391.000637901.000613941.000588891.000877240.997524491.000813600.99309562
0.98800.984548390.984145020.983132100.984268930.984556840.984532750.984512400.984783230.981387610.984385290.97659274
0.96300.959521730.959084480.956857220.959344350.959528660.959505120.959491110.959726190.956388950.958956640.95129317
0.92550.922838900.922337180.919228460.922739810.922843420.922820920.922813980.922998510.919786360.921819680.91448600
0.87250.872599770.871965130.868678910.872498010.872601260.872580000.872579680.872705720.869632830.871087470.86426077
0.80750.807274380.806400800.803966900.807035590.807272600.807252420.807257110.807324780.804346010.805244220.79906808
0.72650.728336570.727108440.726460170.727829870.728331850.728312250.728319900.728338390.725361070.725788820.72036366
0.63450.637138060.635439240.637330750.636250720.637131120.637111310.637119920.637104430.634006070.634085760.62948043
0.53450.536213070.533945650.538751100.534881160.536204950.536184040.536192170.536161460.532811660.532679110.52893879
0.42750.429511290.426618230.434087670.427756220.429503230.429480310.429487770.429461120.425748680.425524130.42267978
0.31850.318774410.315218010.324865800.316651100.318767650.318741880.318749030.318744590.314572560.314360380.31243610
0.20850.207389410.203171210.214055550.205017540.207385160.207355820.207364110.207397110.202702690.202573590.20158869
0.10100.096167080.091295950.102500230.093664220.096166390.096132880.096143760.096227190.090960640.090971240.09093929
−0.0080−0.00832544−0.01375892−0.00409525−0.01070827−0.00832170−0.00835961−0.00834240−0.00820286−0.01401275−0.01388329−0.01294519
−0.1110−0.11093649−0.11690797−0.10961273−0.11309173−0.11092761−0.11097026−0.11094529−0.11074205−0.11711378−0.11684060−0.11492873
−0.2090−0.20924720−0.21570702−0.21178492−0.21104130−0.20923267−0.20928021−0.20924559−0.20897460−0.21589688−0.21548275−0.21259889
−0.3030−0.30086342−0.30773606−0.30836317−0.30213765−0.30084292−0.30089531−0.30084866−0.30050991−0.30793928−0.30741745−0.30356751
RMSE2.425075
× 10−3
3.885510
× 10−3
4.212772
× 10−3
2.742457
× 10−3
2.425088
× 10−3
2.425194
× 10−3
2.425233
× 10−3
2.430500
× 10−3
4.507511
× 10−3
4.169322
× 10−3
6.567087
× 10−3
Table 12. Comparison of the parameters extracted by I-GWO for the single-diode-based model of the PV module.
Table 12. Comparison of the parameters extracted by I-GWO for the single-diode-based model of the PV module.
Design CoefficientsI-GWOABC-DECPSOFPAHFAPSITLBOIWOALBSAPSSATLBO
Iph (A)1.030514531.03181.02861.0320911.030501.030501.030501.03041.03131.03311.031805
Isd (μA)3.482178023.27748.30103.0475383.484203.482303.471703.52333.17563.66423.280945
Rs (Ω)1.201273791.20621.07551.2175831.201301.201301.201601.20141.20531.19891.206000
Rsh (Ω)981.95296539845.24951850.10811.3721984.2813981.9823978.67711020.40714.2857833.3333548.6660
n48.6427414348.394852.243048.1312848.6449048.6428048.6313048.686648.288948.821148.44228
Table 13. IAEs of I-GWO for the double-diode-based model of the PV module.
Table 13. IAEs of I-GWO for the double-diode-based model of the PV module.
VL-measured (V)IL-measured (A)IL-calculated (A)IAE (IL)PL-measured (W)PL-calculated (W)IAE (PL)
0.12481.03151.029020430.002479570.128731200.128421750.00030945
1.80931.03001.027318560.002681441.863579001.858727470.00485153
3.35111.02601.025714130.000285873.438228603.437270630.00095797
4.76221.02201.024114570.002114574.866968404.877038390.01006999
6.05381.01801.022336840.004336846.162768406.189022770.02625437
7.23641.01551.020017930.004517937.348564207.381257730.03269353
8.31891.01401.016496880.002496888.435364608.456135940.02077134
9.30971.01001.010675490.000675499.402797009.409085570.00628857
10.21631.00351.000836620.0026633810.2520570510.224847120.02720993
11.04490.98800.984747500.0032525010.9123612010.876437610.03592359
11.80180.96300.959650980.0033490211.3651334011.325608960.03952444
12.49290.92550.922831750.0026682511.5621789511.528844790.03333416
13.12310.87250.872428740.0000712611.4499047511.448969580.00093517
13.69830.80750.806926770.0005732311.0613772511.053525010.00785224
14.22210.72650.727924720.0014247210.3323556510.352618200.02026255
14.69950.63450.636774820.002274829.326832759.360271430.03343868
15.13460.53450.535983430.001483438.089443708.111894750.02245105
15.53110.42750.429530190.002030196.639545256.671076280.03153103
15.89290.31850.318982400.000482405.061888655.069555400.00766675
16.22290.20850.207729660.000770343.382474653.369977440.01249721
16.52410.10100.096381380.004618621.668934101.592615480.07631862
16.7987−0.0080−0.007913250.00008675−0.13438960−0.132932260.00145734
17.0499−0.1110−0.110712840.00028716−1.89253890−1.887642820.00489608
17.2793−0.2090−0.209423750.00042375−3.61137370−3.618695760.00732206
17.4885−0.3030−0.301361590.00163841−5.29901550−5.270362170.02865333
Table 14. Comparison of I-GWO with other recently published algorithms (the lowest errors are highlighted in boldface).
Table 14. Comparison of I-GWO with other recently published algorithms (the lowest errors are highlighted in boldface).
ModelAlgorithmIph (A)Isd, Isd1 (μA)Isd2 (μA)Rs (Ω)Rsh (Ω)n, n1n2RMSE
Single-diode model for
R.T.C. France PV cell
CARO [66]0.760790.31724-0.0364453.08931.48168-9.8665 × 10−4
CSO [67]0.760780.323-0.0363853.71851.48118-9.8602 × 10−4
CWOA [68]0.760770.3239-0.0363653.79871.4812-9.8602 × 10−4
DE-WOA [58]0.7607760.323021-0.03637753.7185241.481184-9.8602 × 10−4
EHHO [69]0.7607750.323-0.03637553.742821.481238-9.8602 × 10−4
GWO [36]0.7699690.91215-0.0292818.1031.596658-7.5011 × 10−3
GWOCS [36]0.7607730.32192-0.0363953.6321.4808-9.8607 × 10−4
IJAYA [34]0.76080.3228-0.036453.75951.4811-9.8603 × 10−4
ISCE [57]0.760775530.32302083-0.0363770953.718527711.4811836-9.8602 × 10−4
MABC [70]0.7607790.321323-0.03638953.399991.481385-9.8610 × 10−4
MSSO [71]0.7607770.323564-0.0363753.7424651.481244-9.8607 × 10−4
NM-PSO [72]0.760770.32451-0.0363653.82581.48157-9.8605 × 10−4
ORcr-IJADE [47]0.7607760.323021-0.03637753.7185231.481184-9.8602 × 10−4
PGJAYA [73]0.76080.323-0.036453.71851.4812-9.8602 × 10−4
SGDE [74]0.760780.32302-0.03637753.718531.481184-9.8602 × 10−4
TLABC [75]0.760780.32302-0.0363853.716361.48118-9.8602 × 10−4
I-GWO0.760775610.32302197-0.0363770653.717709171.48118398-9.8602 × 10−4
Double-diode model for
R.T.C. France PV cell
CARO [66]0.760750.293150.090980.0364154.39671.473381.773219.8260 × 10−4
CSO [67]0.760780.227320.727850.03673755.38131.451511.997699.8252 × 10−4
CWOA [68]0.760770.24150.60.0366655.20161.456511.98999. 8272 × 10−4
DE-WOA [58]0.7607810.2259740.7493460.0367455.4854371.45101729.824849 × 10−4
EHHO [69]0.7607690170.5861840.2409650.03659883155.639439561.9684514491.4569104099.8360 × 10−4
GWO [36]0.7616680.403020.453380.0326572.527751.6461.55272.2124 × 10−3
GWOCS [36]0.760760.537720.248550.0366654.733121.45889.8334 × 10−4
IJAYA [34]0.76010.00504450.750940.037677.85191.21861.62479.8293 × 10−4
ISCE [57]0.760781080.225974090.749348980.0367404355.485444091.451016729.824849 × 10−4
MABC [70]0.76078210.63069220.241029920.0367121554.75500942.000005381.45685739.8276 × 10−4
MSSO [71]0.7607480.2349250.6715930.03668855.7146621.4542551.9953059.8281 × 10−4
NM-PSO [72]0.760780.238200.648100.0366855.21541.4554429.8259 × 10−4
ORcr-IJADE [47]0.7607810.2259740.7493480.0367455.4854381.45101729.824858 × 10−4
PGJAYA [73]0.76080.210310.885340.036855.81351.44529.8263 × 10−4
SGDE [74]0.760790.28070.249960.0364854.36671.4696551.932289.8441 × 10−4
TLABC [75]0.760810.423940.240110.0366754.667971.90751.456719.8414 × 10−4
I-GWO0.760781880.226284890.746091520.0367397755.461617691.45112761.999998569.824852 × 10−4
Single-diode-based model for
Photowatt-PWP201 PV module
CARO [66]1.031853.28401-1.20556841.321348.40363-2.4270 × 10−3
DE-WOA [58]1.0305143.482263-1.201271981.98214348.642835-2.425075 × 10−3
EHHO [69]1.0304986563.488188406-1.201110352984.496482448.64931708-2.425080 × 10−3
GWO [36]1.030384.9068-1.159261173.796650-2.6749 × 10−3
GWOCS [36]1.030493.465-1.2019982.756648.62367-2.4251 × 10−3
IJAYA [34]1.03053.4703-1.2016977.375248.6298-2.4251 × 10−3
ISCE [57]1.03051433.48226304-1.201271981.982280348.642835-2.425075 × 10−3
NM-PSO [72]1.03043.4888-1.2012992.941548.6498-2.4387 × 10−3
ORcr-IJADE [47]1.0305143.482263-1.201271981.98224148.642835-2.425075 × 10−3
PGJAYA [73]1.03053.4818-1.2013981.854548.6424-2.425075 × 10−3
SGDE [74]1.03053.4823-1.20127981.982248.6428-2.425075 × 10−3
TLABC [75]1.030563.4715-1.20165972.9356748.63131-2.425075 × 10−3
I-GWO1.030514533.48217802-1.20127379981.9529653948.64274143-2.425075 × 10−3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yesilbudak, M. Parameter Extraction of Photovoltaic Cells and Modules Using Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy. Energies 2021, 14, 5735. https://doi.org/10.3390/en14185735

AMA Style

Yesilbudak M. Parameter Extraction of Photovoltaic Cells and Modules Using Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy. Energies. 2021; 14(18):5735. https://doi.org/10.3390/en14185735

Chicago/Turabian Style

Yesilbudak, Mehmet. 2021. "Parameter Extraction of Photovoltaic Cells and Modules Using Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy" Energies 14, no. 18: 5735. https://doi.org/10.3390/en14185735

APA Style

Yesilbudak, M. (2021). Parameter Extraction of Photovoltaic Cells and Modules Using Grey Wolf Optimizer with Dimension Learning-Based Hunting Search Strategy. Energies, 14(18), 5735. https://doi.org/10.3390/en14185735

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