Next Article in Journal
Recent Advances in Carbon Dioxide Conversion: A Circular Bioeconomy Perspective
Previous Article in Journal
A Critical Survey of Environmental Content in United States Undergraduate Mechanical Engineering Curricula
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Extraction of Three Diode Solar Photovoltaic Model Using Improved Grey Wolf Optimizer

1
Department of Electrical Engineering, Faculty of Engineering, Aswan University, Aswan 81542, Egypt
2
Department of Electrical Engineering, Yeungnam University, Gyeongsan 38541, Korea
3
Korea Electric Power Company (KEPCO), Deajon 24056, Korea
*
Authors to whom correspondence should be addressed.
Sustainability 2021, 13(12), 6963; https://0-doi-org.brum.beds.ac.uk/10.3390/su13126963
Submission received: 22 April 2021 / Revised: 14 June 2021 / Accepted: 16 June 2021 / Published: 21 June 2021
(This article belongs to the Special Issue Sustainability of Photovoltaics)

Abstract

:
The enhancement of photovoltaic (PV) energy systems relies on an accurate PV model. Researchers have made significant efforts to extract PV parameters due to their nonlinear characteristics of the PV system, and the lake information from the manufactures’ PV system datasheet. PV parameters estimation using optimization algorithms is a challenging problem in which a wide range of research has been conducted. The idea behind this challenge is the selection of a proper PV model and algorithm to estimate the accurate parameters of this model. In this paper, a new application of the improved gray wolf optimizer (I-GWO) is proposed to estimate the parameters’ values that achieve an accurate PV three diode model (TDM) in a perfect and robust manner. The PV TDM is developed to represent the effect of grain boundaries and large leakage current in the PV system. I-GWO is developed with the aim of improving population, exploration and exploitation balance and convergence of the original GWO. The performance of I-GWO is compared with other well-known optimization algorithms. I-GWO is evaluated through two different applications. In the first application, the real data from RTC furnace is applied and in the second one, the real data of PTW polycrystalline PV panel is applied. The results are compared with different evaluation factors (root mean square error (RMSE), current absolute error and statistical analysis for multiple independent runs). I-GWO achieved the lowest RMSE values in comparison with other algorithms. The RMSE values for the two applications are 0.00098331 and 0.0024276, respectively. Based on quantitative and qualitative performance evaluation, it can be concluded that the estimated parameters of TDM by I-GWO are more accurate than those obtained by other studied optimization algorithms.

1. Introduction

Solar energy plays a large role in saving the environment, helping people socially and economically, and creating jobs and research. Solar energy has many advantages and it is a clean energy source that has many different uses [1,2,3]. In contrast, fossil fuel resources are decreased and have a negative environmental impact; hence, solar energy is considered one of the most promising alternative energy sources [4,5].
Nowadays, solar energy has many applications. One of these applications is installing solar cells in residential units and connecting them to the main electrical network. Connecting residential solar energy systems to the main network requires an accurate model to analyze and describe the voltage and current changes to these systems [6]. There are also many other applications of solar energy, which also require for its development to an accurate mathematical model for the solar system, such as solar water heating, solar heating of buildings, solar-distillation, solar-pumping, solar drying of agricultural and animal products, and solar furnaces [7,8,9].
One of the biggest challenges limiting the usage of solar energy is the huge cost of manufacturing and installation. PV cell has nonlinear current-voltage characteristics and non-linear power-voltage characteristics. Due to these nonlinear properties, the actual operation of the solar cells varies with different operating conditions such as irradiance/temperature conditions, which change over the day and over the seasons. These changes make it difficult to get the best efficiency of the solar cell by relying on the operating conditions of the factory [10].
Modeling of PV systems is considered one of the important stages in the design and manufacturing of PV systems, because of the important information collected from the PV model that helps to predict the PV characteristics under different operating conditions, and thus increase the operating efficiency of these cells [11].
Different models have been proposed in the literature for the PV systems. Single diode model (SDM) and double diode model (DDM) are considered the most popular models in literature. PV cell is a semiconductor PN junction that converts sunlight to an electrical current. The ideal model of the solar cell can be represented by an ideal current source, but the real model of a solar cell should take into consideration the losses in light and current that can be represented electrically by diodes. The model that has a greater number of diodes represents more losses. SDM has one diode that represents the loss in the quasi-neutral zone. A SDM has five estimated parameters [11,12,13,14]. DDM is developed to represent the loss in recombination at lower irradiance by adding a second diode to the model. The total estimated parameters in DDM are seven parameters [15,16,17].
The triple diode model (TDM) is developed to represent the leakage in grain boundaries in PV system by adding a third diode to the model. The total estimated parameters in TDM are nine parameters [18,19,20,21,22,23,24,25,26]. Although the accuracy of the model increased by increasing the number of diodes, the complexity of the model also increased.
The challenge of parameters estimation of these models using optimization algorithms has been discussed in many previous works. A review about the application of evolutionary algorithms (EA) on PV parameters estimation has been presented in a previous study [27]. Different types of evolutionary algorithms (EA) have been discussed in [27], such as EA based on swarm intelligence, EA based on physical theory and EA based on hybridization between different types. They concluded that hybridized techniques have a better accuracy for PV parameters estimation, but with less convergence speed. A new method that combines harries hawk optimization (HHO) with computation method to extract the parameters of the PV TDM has been presented in [28]. Four parameters are extracted by equations and the other five parameters are estimated though HHO. An improved equilibrium optimizer was presented in [29], to estimate the PV parameters of SDM and DDM. The enhancement is achieved using linear reduction diversity technique to improve the diversity of the population until getting the best estimated model parameters that achieve the least RMSE between real measure data and the model which is considered as the best solution. Another improvement for moth-flam optimizer was presented in [30], to estimate the PV parameters of SDM and DDM. The improvement is based on double flame generation. Grasshopper Optimization Algorithm was presented to estimate the PV parameters of TDM in [31]. An intelligent grey wolf optimizer was presented in [32], to estimate the PV panel parameters of SDM and DDM. The idea behind intelligent gray wolf optimizer is the incorporation of opposition based learning to the conventional gray wolf optimizer (GWO) to enhance the exploration and exploitation phases.
This paper is considered a completion in this track, but here the study considers an application of another improvement for GWO which was presented in [33] and discusses the effect of this improvement on PV parameters estimation of TDM. The improved GWO algorithm is called I-GWO. I-GWO is developed with the aim of improving population, exploration and exploitation balance and convergence of the original GWO. The I-GWO is applied to estimate the parameters of TDM for 57 mm diameter commercial silicon R.T.C France solar cell (1000 W/m2 at 33 °C). For more comprehensive results, the I-GWO is applied to estimate the parameters of TDM for Photowatt-PWP201 PV module. This module contains 36 polycrystalline silicon cells connected in series and operating at an irradiance of 1000 W/m2 and temperature of 45 °C. The obtained results in the two applications from I-GWO are compared with the original GWO and other recent optimization algorithms. The obtained results in the two applications are evaluated using different evaluation methods. The accuracy of the studied algorithms are evaluated using root mean square error (RMSE) and absolute error between the real output current and the estimated current by TDM. The speed of the algorithms is compared through the convergence curves. The robustness of the algorithms is compared through the statistical analysis for different independent runs. Minimum, maximum, average and standard deviation statistical values are calculated for all runs. The behaviors of the current-voltage characteristics and power-voltage characteristics at different operating temperatures are discussed for the optimized TDM in the two applications.

2. Mathematical Mode and Optimization Problem

This section discusses the main components of TDM and PV panel as well as the main objective functions of the studied optimization problem.

2.1. Model of TDM and PV Panel

The TDM consists of an ideal current source to represent the ideal solar cell, three diodes to represent different current loses, series and shunt resistance [34,35]. The TDM equivalent circuit is shown in Figure 1. Each component is discussed as:
Current source represents the photo-generated current ( I p h ) .
  • First diode represents the effect of diffusion current   ( I s 1 ) .
  • First diode represents the effect of diffusion current   ( I s 1 ) .
  • Second diode represents the effect of recombination current   ( I s 2 ) .
  • Third diode represents the effect of grain boundaries and large leakage current   ( I s 3 ) .
  • Series resistance represents the semiconductor material resistance at neutral regions   ( R s ) .
  • Shunt resistance represents the current leakage resistance across the P-N junction of PV system   ( R s h ) .
Equations (1) and (2) describe the total output current ( I ) of the mathematical model of TDM.
I = I p h I s 1 I s 2 I s 3 I s h
I = I p h I s 1 [ exp ( q ( V + R s I ) η 1 K T ) 1 ] I s 2 [ exp ( q ( V + R s I ) η 2 K T ) 1 ] I s 3 [ exp ( q ( V + R s I ) η 3 K T ) 1 ] ( V + R s I ) R s h
A combination between series solar cells represents a PV module and a combination between parallel modules represents a PV array as shown in Figure 2. The total output current (I) of the mathematical model of TDM for PV module is described by Equation (3).
I / N p = I p h I s 1 [ exp ( q ( V / N s + R s I / N p ) η 1 K T ) 1 ] I s 2 [ exp ( q ( V / N s + R s I / N s ) η 2 K T ) 1 ] I s 3 [ exp ( q ( V / N s + R s I / N p ) η 3 K T ) 1 ] ( V / N s + R s I / N p ) R s h

2.2. Objective Function

In TDM, nine estimated parameters are considered [ R s ,   R sh ,   I ph , I s 1 , I s 2 , I s 3 , η 1 , η 2 ,   η 3 ] and they can be represented in one vector x = [ x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 , x 9 ] . The optimization problem in this work is to estimate the best values of TDM parameters to reduce the error between experimental and estimated output current; hence the objective function of the current optimization problem is described by Equation (4). The best estimated parameters are those that achieve the lowest RMSE for the objective function, as shown by Equation (5) [36,37]. The objective function of the TDM for PV module is described by Equation (6) [10]. The main steps of the optimization process are summarized in Figure 3.
f T D M ( V , I , X ) = I X 3 + X 4 [ exp ( q ( V + R s I ) X 7 K T ) 1 ] + X 5 [ exp ( q ( V + R s I ) X 8 K T ) 1 ] + X 6 [ exp ( q ( V + R s I ) X 9 K T ) 1 ] + ( V + X 1 I ) X 2
R M S E = 1 N K = 1 N f 2 ( V t m , I t m , X )
where, Vtm and Itm are the measured voltage and current, X is the estimated parameters, N is the maximum number of the measured data.
f T D M ( V , I , X ) = I X 3 + X 4 [ exp ( q ( V / N s + R s I / N p ) X 7 K T ) 1 ] + X 5 [ exp ( q ( V / N s + R s I / N p ) X 8 K T ) 1 ] + X 6 [ exp ( q ( V / N s + R s I / N p ) X 9 K T ) 1 ] + ( V / N s + X 1 I / N p ) X 2

3. Grey Wolf Optimizers

In this section, both original GWO and I-GWO are discussed.

3.1. Grey Wolf Optimizer (GWO)

GWO is inspired by the grey wolf’s manner in attacking their prey [38,39,40,41,42]. GWO determines the way to the best solution based on the first three best solutions inspired from three wolf leaders. The main steps of GWO are described as follows:
The wolf X P ( t ) detects the prey position X P ( t ) to surround it, as shown by Equation (7).
D = | C × X P ( t ) X ( t ) |
Movement towards the prey Equation (8).
X ( t + 1 ) = X P ( t ) A × D
where, A and C are given by Equations (9) and (10), respectively.
A = 2 × A r 1 a ( t )
C = 2 × r 2
where, a ( t ) is a factor decreased depending on the maximum number of iterations as described in (11).
a ( t ) = 2 ( 2 × t ) / M a x I t e r
Tracking the prey considering the best three positions for the three wolf leaders α , β and δ , with positions   X α , X β   and X δ . Each wolf updates its position depending on the three best leaders, as given in Equation (12).
X i 1 ( t ) = X α ( t ) A i 1 × D α ( t )
X i 2 ( t ) = X β ( t ) A i 2 × D β ( t )
X i 3 ( t ) = X δ ( t ) A i 3 × D δ ( t )
where, D α ,   D β and D δ   are calculated from Equation (13).
D α = | C 1 × X α X t |
D β = | C 2 × X β X t |
D δ = | C 3 × X δ X t |
The final best position is updated depending on the three best leaders’ positions as given in Equation (14).
X ( t + 1 ) = X i 1 ( t ) + X i 2 ( t ) + X i 3 ( t ) 3  

3.2. Improved Grey Wolf Optimizer (I-GWO)

I-GWO was developed to improve the performance of the original GWO by reducing the probability of fall in local optima [33]. The improvement is achieved by a new search strategy that includes selecting and updating steps highlighted by a dashed line in flowchart.
The detailed description of the enhancement is as follows:
-
Increasing the initialization by random distribution for N wolfs within the search range as described in Equation (15).
X i j = L B j + r a n d j [ 0 , 1 ] ( U B j L B j )
where, i [ 1 , N ] problem dimension, j [ 1 ,   D ]   problem dimension. LB and UB are the search low and upper limits, respectively.
-
The tracking behavior is enhanced using dimension learning-based hunting technique (DLH). In DLH, each wolf learns from its neighbors. The construction of the neighbors according to the calculated radius are calculated by Equations (16) and (17). The new positions are determined using Equation (18).
The best position is determined by Equation (19) considering X i G W O ( t + 1 ) that determined is by Equation (15).
R i ( t ) = X i ( t ) X i G W O ( t + 1 )
N i ( t ) = { X j ( t ) | D i ( X i ( t ) ,   X j ( t ) ) R i ( t ) ,   X j ( t ) ϵ P o p }
X i = D L H , d ( t + 1 ) = X i , d ( t ) + r a n d × ( X n , d ( d ) X r , d ( t ) )
X i ( t + 1 ) = { X i G W O ( t + 1 )   i f   f ( X i G W O < X i = D L H ) X i = D L H ( t + 1 )   o t h e r w i s e
The flowchart of I-GWO is shown in Figure 3.
The main steps of I-GWO can be described as follows:
(Initialization steps)
-
Random distribution for N wolfs within the search range, as shown by (15)
-
(Main steps)
-
For iteration = 1 to maximum iteration
-
Find three wolf leaders α , β and δ ,
-
Update three wolf leaders’ positions   X α , X β   and X δ using (12) and (13).
-
Calculate the best position X i G W O using (14)
-
(Updating improvement steps (DLH))
-
Calculate radius to construction of the neighbors using (16)
-
Determine the neighbors using Equation (17)
-
Calculate new solution X i = D L H , d ( t + 1 ) using (18)
-
Select the best position between X i = D L H , d ( t + 1 ) and X i G W O ( t + 1 ) using (19)
-
End for loop

4. Results

This section discusses the performance of I-GWO when it is applied to estimate the parameters of TDM for different real PV systems. Different evaluation methods are used for results, evaluation and analysis. Two different applications, the real data from RTC furnace and the real data of PTW polycrystalline PV panel, are considered in the analysis.

4.1. Application #1:

In this application, the real data of 57 mm diameter commercial silicon R.T.C France solar cell (under 1000 W/m2 at 33 °C) are used for parameters estimation of the TDM [43]. A comparison between the obtained root mean square error (RMSE) values of I-GWO and other recent optimization algorithms is presented in Table 1. The parameters of each algorithm are listed in Table 2. The convergence curve of all algorithms is shown in Figure 4. The I-GWO has a RMSE value better than the original GWO as well as other compared algorithms. Moreover, I-GWO has an accepted convergence speed as shown in Figure 4. The stability of I-GWO is tested from the statistical analysis of 30 independent runs as observed in Table 3 and the boxplot in Figure 5. A graphical display of the absolute error between the real measured output current values and the model output value in Equation (20) is shown in Figure 6, Figure 7 and Figure 8 to depict the characteristic curves of the current-voltage and the power curve for the real PV measured data and the calculated by different algorithms. The characteristic curve of the current-voltage and the power curve for the TDM estimated by I-GWO at different temperatures is displayed in Figure 9 and Figure 10, respectively.
C u r r e n t   A b s o l u t   e r r o r = ( I I e s t i m a t e d ) 2 2

4.2. Application #2:

In this application, the real data of Photowatt-PWP201 PV module are used [44,45,46]. This module contains 36 polycrystalline silicon cells connected in series and operating at an irradiance of 1000 W/m2 and temperature of 45 °C. A comparison between the obtained root mean square error (RMSE) value of I-GWO and other recent optimization algorithms is presented in Table 4. The convergence curves of all algorithms are shown in Figure 11. The I-GWO has a RMSE value better than the original GWO, as well as other compared algorithms, except it is slightly better than HBO. The stability of I-GWO is tested from the statistical analysis of 30 independent runs as shown in Table 5 and the boxplot in Figure 12. A graphical display of the absolute error between the real measured output current values and the model output value is shown in Figure 13, Figure 14 and Figure 15 show the characteristic curves of the current-voltage and the power curve for the real PV measured data and the calculated by different algorithms. The characteristic curves of the current-voltage and the power curves for the TDM estimated by I-GWO at different temperatures are shown in Figure 16 and Figure 17, respectively.

5. Overall Discussion

The reliability of I-GWO algorithm has been evaluated by different comparisons with the original GWO and three well known recent algorithms (HBO, MRFO and HHO). I-GWO has been tested for parameters estimation of solar cell and PV panel which considered the most popular real PV system. In the two applications, the I-GWO wins to achieve the best RMSE compared with the four studied algorithms. In Application #1, the RMSE of I-GWO is 0.00098331. This value is better than those obtained by different algorithms proposed in literature such as TLABC, GOTLBO and TLBO [46]. In addition, the accuracy of I-GWO is clear from the current absolute error figure as it has the least error over different measured data. The convergence speed of I-GWO is better than the other algorithms as observed in Figure 4 and Figure 11. The statistical analysis for the results in the two applications prove the robustness of I-GWO. In Application #1, I-GWO has the lowest minimum and standard deviation values. In Application #2, the minimum value obtained by I-GWO is very close to those obtained by MRFO and HHO and better than others, but it has the lowest standard deviation value. Testing the estimated TDM optimized by I-GWO at different operating temperatures is verified by drawing current–voltage and power–voltage characteristic curves at temperatures for the two applications as shown in Figure 9, Figure 10, Figure 16, Figure 17, respectively.

6. Conclusions

This paper discussed the parameters estimation of TDM PV system which is considered one of important topics in the real world due to its cost-impact in PV system manufacturing. This paper proposed an application for I-GWO algorithm to estimate more accurate parameters of PV model. Two real examples were discussed: 57 mm diameter commercial silicon R.T.C France solar cell and polycrystalline PV panel. Many evaluation methods were applied, namely, a comparison of the obtained results between I-GWO and others recent optimization algorithms as well as applying different evaluation factors such as RMSE, absolute error and statistical analysis. The results of the two applied applications are extended to include quantitative and qualitative results. The accuracy of algorithms is compared through RMSE and absolute error as well as statistical analysis to cover different quantitative results. The behavior of algorithms is compared though different qualitative results such as algorithm convergence speed and boxplot graph for all statistical results. This study presented a comparison and discussion of current-voltage characteristics and power-voltage characteristics between real measured data set and PV model optimized by different algorithms. The results show that I-GWO has a significant improvement in accuracy than original GWO and other algorithms and this is made clear by the RMSE for the two applications. This study is considered a starting point for research work that focuses on studying the applicability of this improved algorithm on modeling of complex and large PV systems.

Author Contributions

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

Funding

This work was supported by the “Development of Modular Green Substation and Operation Technology” of the Korea Electric Power Corporation (KEPCO).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolDescription
TDMThree Diode Model
RMSERoot Mean Square Error
PVPhoto Voltaic
Nsnumber of series solar cells
IPV module output current
NPnumber of parallel PV modules
VTerminal voltage
XP(t)Prey current position
IphPhoto generated current source
X(t)Gray wolf current position
ɳ1First Diode Ideality factor (Diffusion current components)
ɳ2Second Diode Ideality Factor (Recombination current components)
ɳ3Third diode Ideality Factor(Leakage current components)
RshEquivalent Shunt resistance for current leakage resistance across the P-N junction of solar cell
RsEquivalent Series resistance for semiconductor material at neutral regions
K=1.380 × 10−23(J/Ko) Boltzmann constant
Is1First diode current
q1.602 × 10−19 (C) Coulombs.
Is2Second diode current
T (Ko)Photo cell temperature (Kelvin)
Is3Third diode current
HBOHeap-based optimizer
X(t)Current position
X(t+1)Position in next iteration
GWOGrey Wolf Optimizer
I-GWOImproved Grey Wolf Optimizer
HHOHarries Hawk optimization
MRFOManta Ray Foraging Optimization

References

  1. Humada, A.M.; Darweesh, K.G.M.; Kamil, M.; Samen, F.M.; Naseer, K.K.; Tahseen, A.T.; Omar, I.A.; Saad, M. Modeling of PV system and parameter extraction based on experimental data: Review and investigation. Sol. Energy 2020, 199, 742–760. [Google Scholar] [CrossRef]
  2. Gielen, D.; Boshell, F.; Saygin, D.; Bazilian, M.D.; Wagner, N.; Gorini, R. The role of renewable energy in the global energy transformation. Energy Strategy Rev. 2019, 24, 38–50. [Google Scholar] [CrossRef]
  3. Araújo, K.; Boucher, J.L.; Aphale, O. A clean energy assessment of early adopters in electric vehicle and solar photovoltaic technology: Geospatial, political and socio-demographic trends in New York. J. Clean. Prod. 2019, 216, 99–116. [Google Scholar] [CrossRef]
  4. Choudhary, P.; Srivastava, R.K. Sustainability perspectives—A review for solar photovoltaic trends and growth opportunities. J. Clean. Prod. 2019, 227, 589–612. [Google Scholar] [CrossRef]
  5. Singh, S.N.; Prabhakar, T.; Sumit, T. Introduction to solar energy. In Fundamentals and Innovations in Solar Energy; Springer: Singapore, 2021; pp. 1–9. [Google Scholar]
  6. Piyush, G. Importance of Detailed Modeling of Loads/PV Systems Connected to Secondary of Distribution Transformers. Master’s Thesis, Electrical Engineering Faculty of the Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2017. [Google Scholar]
  7. Halabi, M.A.; Al-Qattana, A.; Al-Otaibi, A. Application of solar energy in the oil industry—Current status and future prospects. Renew. Sustain. Energy Rev. 2015, 43, 296–314. [Google Scholar] [CrossRef]
  8. Mekhilef, S.; Faramarzi, S.Z.; Saidur, R.; Salam, Z. The application of solar technologies for sustainable development of agricultural sector. Renew. Sustain. Energy Rev. 2013, 18, 583–594. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Sivakumar, M.; Yang, S.; Enever, K.; Ramezanianpour, M. Application of solar energy in water treatment processes: A review. Desalination 2018, 428, 116–145. [Google Scholar] [CrossRef] [Green Version]
  10. 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]
  11. Rodrigues, E.M.G.; Melicio, R.; Mendes, V.M.F.; Catalao, J.P. Simulation of a solar cell considering single-diode equivalent circuit model. RE&PQJ 2011. [Google Scholar] [CrossRef]
  12. Ma, J.; Man, K.L.; Ting, T.O.; Zhang, N.; Guan, S.U.; Wong, P.W. Wong approximate single-diode photovoltaic model for efficient I-V characteristics estimation. Sci. World J. 2013, 2013, 230471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sabadus, A.; Mihailetchi, V.; Paulescu, M. Parameters extraction for the one-diode model of a solar cell. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2017; p. 40005. [Google Scholar]
  14. Tamrakar, V.; Gupta, S.C.; Sawle, Y. Single-diode Pv cell modeling and study of characteristics of single and two-diode equivalent circuit. Electr. Electron. Eng. Int. J. 2015, 4, 12. [Google Scholar] [CrossRef]
  15. Tamrakar, V.; Gupta, S.C.; Sawle, Y. Single-diode and two-diode Pv cell modeling using matlab for studying characteristics of solar cell under varying conditions. Electr. Comput. Eng. Int. J. 2015, 4, 67–77. [Google Scholar] [CrossRef]
  16. Sulyok, G.; Summhammer, J. Extraction of a photovoltaic cell’s double-diode model parameters from data sheet values. Energy Sci. Eng. 2018, 6, 424–436. [Google Scholar] [CrossRef] [Green Version]
  17. Tanvir, A.; Sobhan, S.; Nayan, M.F. Comparative analysis between single diode and double diode model of PV cell: Concentrate different parameters effect on its efficiency. J. Power Energy Eng. 2016, 4, 31–46. [Google Scholar]
  18. Soliman, M.A.; Al-Durra, A.; Hasanien, H.M. Electrical parameters identification of three-diode photovoltaic model based on equilibrium optimizer algorithm. IEEE Access 2021, 9, 41891–41901. [Google Scholar] [CrossRef]
  19. Wang, R. Parameter identification of photovoltaic cell model based on enhanced particle swarm optimization. Sustainability 2021, 13, 840. [Google Scholar] [CrossRef]
  20. 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]
  21. El-Hameed, M.A.; Elkholy, M.M.; El-Fergany, A.A. Three-diode model for characterization of industrial solar generating units using Manta-rays foraging optimizer: Analysis and validations. Energy Convers. Manag. 2020, 219, 113048. [Google Scholar] [CrossRef]
  22. Harrag, A.; Daili, Y. Three-diode PV model parameters extraction using PSO algorithm. Rev. Energ. Renouvelables 2019, 22, 85–91. [Google Scholar]
  23. Shekoofa, O.; Wang, J. Multi-diode modeling of multi-junction solar cells In Proceedings of the 23rd Iranian Conference on Electrical Engineering, Tehran, Iran, 10–14 May 2015.
  24. Ukoima, K.N.; Ekwe, O. A three-diode model and simulation of photovoltaic (Pv) cells. J. Eng. Technol. 2019, 5, 108–116. [Google Scholar]
  25. Bayoumi, A.S.; El-Sehiemy, R.A.; Mahmoud, K.; Lehtonen, M.; Darwish, M.M. Assessment of an improved three-diode against modified two-diode patterns of MCS solar cells associated with soft parameter estimation paradigms. Appl. Sci. 2021, 11, 1055. [Google Scholar] [CrossRef]
  26. Mehdi, O.; Mohamed, S.M.; Djalel, D.O. Comprehensive three-diode model of photovoltaic array with partial shading capability. Int. J. Power Energy Convers. 2018, 9, 159–173. [Google Scholar]
  27. Saha, C.; Agbu, N.; Jinks, R.; Huda, M.N. Review article of the solar PV parameters estimation using evolutionary algorithms. MOJ Sol. Photoenergy Syst. 2018, 2, 63–75. [Google Scholar]
  28. Qais, M.H.; Hasanien, H.M.; Alghuwainem, S. Parameters extraction of three-diode photovoltaic model using computation and Harris Hawks optimization. Energy 2020, 195, 117040. [Google Scholar] [CrossRef]
  29. 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]
  30. 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]
  31. Elazab, O.S.; Hasanien, H.M.; Alsaidan, I.; Abdelaziz, A.Y.; Muyeen, S.M. Parameter estimation of three diode photovoltaic model using grasshopper optimization algorithm. Energies 2020, 13, 497. [Google Scholar] [CrossRef] [Green Version]
  32. Saxena, A.; Sharma, A.; Shekhawat, S. Parameter extraction of solar cell using intelligent grey wolf optimizer. Evol. Intell. 2020. [Google Scholar] [CrossRef]
  33. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S. An improved grey wolf optimizer for solving engineering problems. Expert Syst. Appl. 2021, 166, 113917. [Google Scholar] [CrossRef]
  34. Qais, M.H.; Hasanien, H.M.; Alghuwainem, S.; Nouh, A.S. Coyote optimization algorithm for parameters extraction of three-diode photovoltaic model of photovoltaic modules. Energy 2019, 187, 116001. [Google Scholar] [CrossRef]
  35. Qais, M.H.; Hasanien, H.M.; Alghuwainem, S. Identification of electrical parameters for three-diode photovoltaic model using analytical and sunflower optimization algorithm. Appl. Energy 2019, 250, 109–117. [Google Scholar] [CrossRef]
  36. Mohamed, A.; Reda, M.; Attia, E.; Mohamed, A.; Askar, S.S. Parameters identification of PV triple-diode model using improved generalized normal distribution algorithm. Mathematics. 2021, 9, 995. [Google Scholar]
  37. 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]
  38. Mirjalili, S.; Mirjalili, S.M.; Lewisa, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  39. Roy, R.G.; Ghoshal, D. Grey wolf optimization-based second order sliding mode control for inchworm robot. Robotica 2019, 38, 1–19. [Google Scholar] [CrossRef]
  40. Luo, K.; Zhao, O. A binary grey wolf optimizer for the multidimensional knapsack problem. Appl. Soft Comput. 2019, 83, 105645. [Google Scholar] [CrossRef]
  41. Pal, A.; Bahuguna, S. Grey wolf optimizer. Future Asp. Eng. Sci. Technol. 2018, 36, 379–382. [Google Scholar]
  42. Sharma, I.; Chahar, V.; Agri, S. A comprehensive survey on grey wolf optimization. Recent Pat. Comput. Sci. 2020, 13, 1–13. [Google Scholar] [CrossRef]
  43. Abdelghany, R.Y.; Kamel, S.; Ramadan, A.; Sultan, H.; Rahmann, C. Solar cell parameter estimation using school-based optimization algorithm. In Proceedings of the IEEE International Conference on Automation/XXIV Congress of the Chilean Association of Automatic Control, Santiago, Chile, 22–26 March 2021. [Google Scholar]
  44. Ramadan, A.; Kamel, S.; Korashy, A.; Yu, J. Photovoltaic cells parameter estimation using an enhanced teaching learning based optimization algorithm. Iran. J. Sci. Technol. 2019, 44, 767–779. [Google Scholar] [CrossRef]
  45. Ramadan, A.; Kamel, S.; Hussein, N.M.; Hassan, M.H. A new application of chaos game optimization algorithm for parameters extraction of three diode photovoltaic model. IEEE Access 2021, 9, 51582–51594. [Google Scholar] [CrossRef]
  46. Liao, Z.; Chen, Z.; Li, S. Parameters extraction of photovoltaic models using triple-phase teaching-learning-based optimization. IEEE Access 2019, 7, 77629–77641. [Google Scholar] [CrossRef]
Figure 1. PV TDM.
Figure 1. PV TDM.
Sustainability 13 06963 g001
Figure 2. PV module.
Figure 2. PV module.
Sustainability 13 06963 g002
Figure 3. I-GWO flowchart.
Figure 3. I-GWO flowchart.
Sustainability 13 06963 g003
Figure 4. The convergence curve of all studied algorithms.
Figure 4. The convergence curve of all studied algorithms.
Sustainability 13 06963 g004
Figure 5. Boxplot of all algorithms for 30 independent runs.
Figure 5. Boxplot of all algorithms for 30 independent runs.
Sustainability 13 06963 g005
Figure 6. Current absolute error of all studied algorithms.
Figure 6. Current absolute error of all studied algorithms.
Sustainability 13 06963 g006
Figure 7. Current-voltage curve for experimental current and TDM current for different algorithms.
Figure 7. Current-voltage curve for experimental current and TDM current for different algorithms.
Sustainability 13 06963 g007
Figure 8. Power-voltage curve for experimental power and TDM power for different algorithms.
Figure 8. Power-voltage curve for experimental power and TDM power for different algorithms.
Sustainability 13 06963 g008
Figure 9. Current-voltage curve for experimental current and TDM current for I-GWO at different temperatures.
Figure 9. Current-voltage curve for experimental current and TDM current for I-GWO at different temperatures.
Sustainability 13 06963 g009
Figure 10. Power-voltage curve for experimental power and TDM power for I-GWO at different temperatures.
Figure 10. Power-voltage curve for experimental power and TDM power for I-GWO at different temperatures.
Sustainability 13 06963 g010
Figure 11. The convergence curve of all algorithms.
Figure 11. The convergence curve of all algorithms.
Sustainability 13 06963 g011
Figure 12. Boxplot for all algorithms for 30 independent runs.
Figure 12. Boxplot for all algorithms for 30 independent runs.
Sustainability 13 06963 g012
Figure 13. Current absolute error of all algorithms.
Figure 13. Current absolute error of all algorithms.
Sustainability 13 06963 g013
Figure 14. Current-voltage curve for experimental current and TDM current for different algorithms.
Figure 14. Current-voltage curve for experimental current and TDM current for different algorithms.
Sustainability 13 06963 g014
Figure 15. Power-voltage curve for experimental power and TDM power for different algorithms.
Figure 15. Power-voltage curve for experimental power and TDM power for different algorithms.
Sustainability 13 06963 g015
Figure 16. Current-voltage curve for experimental current and TDM current for I-GWO at different temperatures.
Figure 16. Current-voltage curve for experimental current and TDM current for I-GWO at different temperatures.
Sustainability 13 06963 g016
Figure 17. Power-voltage curve for experimental power and TDM power for I-GWO at different temperatures.
Figure 17. Power-voltage curve for experimental power and TDM power for I-GWO at different temperatures.
Sustainability 13 06963 g017
Table 1. Estimated parameters and RMSE for I-GWO and other optimization algorithms.
Table 1. Estimated parameters and RMSE for I-GWO and other optimization algorithms.
IGWOGWOHBOMRFOHHO
Rs(Ω)0.03670.0438210.0400.036340.018333
Rsh(Ω)54.88859.6459.99753.924692.25011
Iph(A)0.76070.76190.76080.760780.764769
Isd1(A)2.27 × 10−71.38 × 10−106.97 × 10−72.67 × 10−83.82 × 10−6
Isd2(A)3.14 × 10−72.51 × 10−101.00 × 10−101.54 × 10−82.71 × 10−6
Isd3(A)2.34 × 10−74.07 × 10−61.594723.17 × 10−71.28 × 10−6
N11.92561.66771.0091.90761.882327
N21.96001.01811.00821.86741.891646
N31.45001.91821.30831.4751.891677
RMSE0.000983310.001910.0011200.0009860020.007721
Table 2. Parameters setting for each studied algorithm.
Table 2. Parameters setting for each studied algorithm.
AlgorithmParameter Setting
I-GWOr1 = rand()r2 = rand()
GWOr1 = rand()r2 = rand()
HBOdegree = 3
MRFONP = 1000S = 2
Table 3. The statistical results of all algorithms.
Table 3. The statistical results of all algorithms.
MinimumAverageMaximumSTD
I-GWO0.0009830.0009840.0009856.60404 × 10−7
GWO0.0012980.007510.0193190.010231343
HBO0.001120.0016066670.00240.000692917
MRFO0.0009860.0009870.0009891.25983 × 10−6
HHO0.00772050.01174750.019170.006435824
Table 4. Estimated parameters and RMSE for I-GWO and other optimization algorithms.
Table 4. Estimated parameters and RMSE for I-GWO and other optimization algorithms.
I-GWOGWOHBOMRFOHHO
Rs(Ω)1.1986837731.6661381.1995821.2106093731.804672596
Rsh(Ω)986.336588660983.629799.9841045358.9956246
Iph(A)1.0305088461.1533991.0304471.0320379751.025407497
Isd1(A)1.25 × 10−61.48 × 10−103.48 × 10−76.20 × 10−86.07 × 10−10
Isd2(A)7.77 × 10−71.83 × 10−101.00 × 10−101.97 × 10−61.24 × 10−9
Isd3(A)1.54 × 10−62.35 × 10−103.19 × 10−61.09 × 10−64.89 × 10−10
N149.266722628.0999748.6026748.226630729.89953654
N248.3842266128.1307549.5294748.0882814333.55517196
N348.2559337546.0348948.5655848.0944238229.53540241
RMSE0.00242762910.033230.00242810.00246090.024273
Table 5. The statistical results of all algorithms.
Table 5. The statistical results of all algorithms.
MinimumAverageMaximumSTD
I-GWO0.0024276290.0024320.0024385.26003 × 10−6
GWO0.033230.0405630.045230.006429101
HBO0.00242810.0024650.0025285.50757 × 10−5
MRFO0.00246090.0025540.0026419.0185 × 10−5
HHO0.0242730.0248060.0252730.000503322
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ramadan, A.-E.; Kamel, S.; Khurshaid, T.; Oh, S.-R.; Rhee, S.-B. Parameter Extraction of Three Diode Solar Photovoltaic Model Using Improved Grey Wolf Optimizer. Sustainability 2021, 13, 6963. https://0-doi-org.brum.beds.ac.uk/10.3390/su13126963

AMA Style

Ramadan A-E, Kamel S, Khurshaid T, Oh S-R, Rhee S-B. Parameter Extraction of Three Diode Solar Photovoltaic Model Using Improved Grey Wolf Optimizer. Sustainability. 2021; 13(12):6963. https://0-doi-org.brum.beds.ac.uk/10.3390/su13126963

Chicago/Turabian Style

Ramadan, Abd-ElHady, Salah Kamel, Tahir Khurshaid, Seung-Ryle Oh, and Sang-Bong Rhee. 2021. "Parameter Extraction of Three Diode Solar Photovoltaic Model Using Improved Grey Wolf Optimizer" Sustainability 13, no. 12: 6963. https://0-doi-org.brum.beds.ac.uk/10.3390/su13126963

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