Estimation of Single-Diode and Two-Diode Solar Cell Parameters by Using a Chaotic Optimization Approach

Estimation of single-diode and two-diode solar cell parameters by using chaotic optimization approach (COA) is addressed. The proposed approach is based on the use of experimentally determined current-voltage (I-V) characteristics. It outperforms a large number of other techniques in terms of average error between the measured and the estimated I-V values, as well as of time complexity. Implementation of the proposed approach on the I-V curves measured in laboratory environment for different values of solar irradiation and temperature prove its applicability in terms of accuracy, effectiveness and the ease of implementation for a wide range of practical environment conditions. The COA-based parameter estimation is, therefore, useful for PV power converter designers who require fast and accurate model for PV cell/module.


Introduction
The contribution of solar energy in total electric energy production is growing constantly. As the price of solar inverters and solar panels constantly decreases, most countries are basing their energy policy on higher use of solar energy. Studies on energy networks, and especially testing of the integration of solar energy sources into power networks, requires accurate calculation of the solar output power, as well as accurate modeling of solar cells. For that reason, modeling of solar cells (corresponding equivalent circuit and accurate parameters value) represents a very popular research field.
In the literature, two basic models of the equivalent circuits of solar cell can be found, namely the single-diode model (SDM) [1] and the double-diode model (DDM) [2]. DDM considers the composite effect of the neutral region of the junction, and, therefore, models the solar cells more accurately [3]. However, it is characterized by seven unknown parameters. Because of the complexity of DDM, some authors reduce the number of unknown parameters [3,4], which can greatly affect the model accuracy [5]. In this work, we focus on both SDM and DDM, without neglecting any of the model parameters of the solar cell.
For the estimation of solar cell parameters, two main sets of "input data" and corresponding estimations can be found, namely (a) estimation based on datasheet information [6,7] and (b) estimation based on experimental data [8].
The former uses the datasheet information (open circuit voltage, short circuit current, voltage and current value at maximum power point characteristics) provided by photovoltaic (PV) manufacturers under standard test conditions. However, recent research [7] on the usage of datasheet values for solar cell parameter estimation shows that current-voltage characteristic is not unique when designers focus on three datasheet points (open circuit, short circuit and maximum power). It is shown that by observing only three points, we can have multiple I-V characteristics, although in reality a solar cell has a single defined I-V characteristic corresponding to a specific set of cell parameters. To obtain unique and accurate I-V characteristic of PV cell, experimental data on more than three major points are necessary [7]. Research [7] has also implied that only approaches based on experimental data generate accurate models.
Evaluating the performance of solar cells (or PV panels) requires as accurate an estimation of the equivalent circuit solar cell parameters as possible. The approaches used for this purpose can be categorized as follows: (a) analytical techniques [9][10][11][12][13], (b) numerical extraction [13] and (c) meta-heuristic techniques .
Analytical techniques provide mathematical expressions for solving equivalent circuit parameters based on some input data (manufacturer data or/and data obtained from measurements). A review and comparative assessment of non-iterative methods for the extraction of the single-diode model parameters of photovoltaic modules is given in [11]. In general, analytical techniques provide rapid solution. On the other hand, these techniques give erroneous results when the estimated and measure solar cell output characteristics are compared [12].
Numerical techniques are based on curve fitting, usually via iterative methods. However, the application of curve fitting to nonlinear diode equations is quite complex, making numerical determination of solar cell parameters unpopular [14].
Recently, meta-heuristic algorithms for solar cell parameter estimation have been proposed . They impose no restrictions on the problem formulation, they are excellent in dealing with nonlinear equations, and they can be applied for different numbers of unknown parameters.
Of all the proposed techniques, none excels in terms of accuracy and efficiency with respect to others. This was our main incentive for doing research in this field. We propose both accurate and efficient parameter optimization of solar cell SDM and DDM through chaotic optimization approach (COA).
Recently, COA has been used in solving various optimization problems: parameter identification of Jiles-Atherton hysteresis model [59], single-phase transformer parameter estimation [60], design of PID parameters for automatic voltage regulation of synchronous machine [61], antenna array radiation pattern synthesis [62,63]. The main advantages of COA over other optimization techniques are easy implementation and short execution time [64]. It should be noted that different versions of chaotic algorithm have also been used in solar cell parameters estimation, namely chaotic heterogeneous comprehensive learning particle swarm optimizer variants [16], chaotic asexual reproduction optimization [33], mutative-scale parallel chaos optimization algorithm [41], chaos-embedded gravitational search algorithm [65], chaotic improved artificial bee colony algorithm [66], improved chaotic whale optimization algorithm [67], etc. Unlike methods proposed in [16,33,41,[65][66][67], this paper will use COA based on Logistic map for solar cell parameter estimation. Logistic map has a simple form with one variable and one control parameter, and it can produce chaotic behavior similar to more complex chaotic systems. [68]. The power of this optimization approach is demonstrated in [60], where it was used for the estimations of the transformer's parameters.
The effectiveness of the proposed approach will be evaluated on different solar cells (different with respect to solar cell voltage and current level) found in the literature and the laboratory environment. Furthermore, COA-based parameter estimation will be compared with 50 various literature techniques for SDM, and with 12 various techniques for DDM. Also, we will compare parameters obtained by using the proposed method on the measured data with analytically and numerically obtained parameters. Finally, we will apply COA on the measured I-V characteristics using the laboratory environment.
The paper is organized as follows. SDM and DDM for solar cells are described in Section 2. In Section 3, COA and its implementation for solar cell parameters estimation are described. The comparison of solar cell parameters estimation accuracy obtained by using the COA-based and other methods, for one solar cell and one solar module, is presented in Section 4. The experimental setup for measuring I-V curves is presented in Section 5, along with the COA-based parameter estimation results. The concluding remarks are given in Section 6.

Mathematical Modeling of Single and Double Solar Cells
SDM is commonly used model for solar cell representation [1], and the equivalent circuit is shown in Figure 1a. The I-V relationship for this model can be described by the following equation: where I pv is the photo-generated current, R s the series parasitic resistance, R p the parallel parasitic resistance, I 0 the saturation current, n is ideality factor and V th = k B T/q is the thermal voltage (k B is Boltzmann constant equal to 1.38 × 10 −23 J/K, T the temperature and q the electron charge equal to 1.602 × 10 −19 C).
Energies 2019, 12, x FOR PEER REVIEW 3 of 15 numerically obtained parameters. Finally, we will apply COA on the measured I-V characteristics using the laboratory environment. The paper is organized as follows. SDM and DDM for solar cells are described in Section 2. In Section 3, COA and its implementation for solar cell parameters estimation are described. The comparison of solar cell parameters estimation accuracy obtained by using the COA-based and other methods, for one solar cell and one solar module, is presented in Section 4. The experimental setup for measuring I-V curves is presented in Section 5, along with the COA-based parameter estimation results. The concluding remarks are given in Section 6.

Mathematical Modeling of Single and Double Solar Cells
SDM is commonly used model for solar cell representation [1], and the equivalent circuit is shown in Figure 1a. The I-V relationship for this model can be described by the following equation: where Ipv is the photo-generated current, Rs the series parasitic resistance, Rp the parallel parasitic resistance, I0 the saturation current, n is ideality factor and Vth = kBT/q is the thermal voltage (kB is Boltzmann constant equal to 1.38•× 10 −23 J/K, T the temperature and q the electron charge equal to 1.602•× 10 −19 C). The equivalent circuit with DDM for the solar cell is shown in Figure 1b. Therefore, unlike the SDM model, the DDM model of the solar cell, in addition to the rectifying diode, includes one more diode to consider the space charge recombination current [48]. The I-V characteristic of DDM is given as where Io1 and Io2 are the diffusion and saturation currents, whereas n1 and n2 are the diffusion and recombination diode ideality factors [48]. The ideality factor is discussed in [69,70], whereas [71] presents a method for ideality factor calculation.

COA and Objective Function
COA is a very powerful optimization technique that has found numerous scientific applications [59][60][61][62][63]. This approach is based on the theory of chaos, which is, in a mathematical sense, described by ordinary differential equations or by an iterative map [64].
Different chaotic systems, including the logistic map, lozi map, tent map and Lorenz system, can be found in the literature. In this paper, we will base COA on the logistic map [59][60][61][62][63][64].
The task of COA is to estimate a set of unknown parameters X which minimizes the objective function (OF). In our case, for SDM, X = [Rs, Rp, Ipv, Io, n], and for DDM, X = [Rs, Rp, Ipv, Io1, Io2, n1, n2]. The equivalent circuit with DDM for the solar cell is shown in Figure 1b. Therefore, unlike the SDM model, the DDM model of the solar cell, in addition to the rectifying diode, includes one more diode to consider the space charge recombination current [48]. The I-V characteristic of DDM is given as where I o1 and I o2 are the diffusion and saturation currents, whereas n 1 and n 2 are the diffusion and recombination diode ideality factors [48]. The ideality factor is discussed in [69,70], whereas [71] presents a method for ideality factor calculation.

COA and Objective Function
COA is a very powerful optimization technique that has found numerous scientific applications [59][60][61][62][63]. This approach is based on the theory of chaos, which is, in a mathematical sense, described by ordinary differential equations or by an iterative map [64].
Different chaotic systems, including the logistic map, lozi map, tent map and Lorenz system, can be found in the literature. In this paper, we will base COA on the logistic map [59][60][61][62][63][64].
The task of COA is to estimate a set of unknown parameters X which minimizes the objective function (OF). In our case, for SDM, X = [R s , R p , I pv , I o , n], and for DDM, X = [R s , R p , I pv , I o1 , I o2 , n 1 , n 2 ].
whereas for DDM it reads where P is the number of measured I-V pairs from the I-V characteristics, and V t and I t represent the voltage and current value of pair t. Figure 2 presents the search procedure, i.e., the COA flowchart. The detailed description of COA flowchart can be found in [59].
whereas for DDM it reads where P is the number of measured I-V pairs from the I-V characteristics, and Vt and It represent the voltage and current value of pair t. Figure 2 presents the search procedure, i.e., the COA flowchart. The detailed description of COA flowchart can be found in [59].
In this paper, the following COA parameters were used: M = 1000, N = 50000. The COA-based estimation is compared with other approaches through the root mean square error (RMSE), defined as follows: where Iest,k and Imeas,k represent the estimated and the measured values of solar output current in point k, respectively.  In this paper, the following COA parameters were used: M = 1000, N = 50,000. The COA-based estimation is compared with other approaches through the root mean square error (RMSE), defined as follows: where I est,k and I meas,k represent the estimated and the measured values of solar output current in point k, respectively.

Simulation Results
To evaluate COA for solar cell parameters estimation, we first applied the proposed method to an experimental current-voltage characteristic extracted from the manufacturer's datasheets of a well-known R.T.C. France solar cell operating under standard test conditions.
The values of parameters obtained by using COA for the R.T.C. France solar cell are summarized, by year of publication, in Table 1, for SDM and Table 2 for DDM. These values are compared with the values of parameters published in recent papers (column Reference) for the same experimental data. During the estimation process, the parameter ranges for SDM estimation were R s (Ω) ∈ [0.02, 0.05], [50,55] and n ∈ [1.35, 1.6], whereas for DDM,  Tables 1 and 2 report parameters as they appear in the cited papers with no modification. However, in some papers in the Energy Conversion and Management journal (in Table 1 marked by *), inaccuracies occurred in parameter estimation of the PV cell using metaheuristic techniques. Namely, the results proposed in [15][16][17][18]21] do not correspond to the objective function [56].
The presented results, especially the value of RMSE, show that COA offers solar characteristics closer to the measured characteristics than the other existing methods, i.e., it outperforms other methods in terms of accuracy. In addition, by observing Tables 1 and 2, it is also evident that DDM characterizes solar cells more accurately than SDM, which supports the conclusion regarding DDM accuracy noted in [3].
It can be seen that COA outperforms several other techniques, such as evaporation rate-based water cycle algorithm (ER-WPA) [19] and cat swarm optimization (CSO) [24] for SDM, and with the generalized opposition-flower pollination algorithm-nelder-mead simplex method (GOFPANM) [53], by a small margin. However, the implementation of COA is simpler than implementation of ER-WPA, CSO and GOFPANM. Furthermore, COA is computationally less demanding than CSO since CSO requires changing the operation mode during the estimation process [24]. On the other side, GOFPANM is a hybrid algorithm which combines local and global search as well as different algorithms during estimation [53]. In general, most evolutionary algorithms have the complexity of O((np + C of p)N i ), where O is the big O notation, n is the dimension of the parameter space, p is the population size, N i is the number of iterations and C of is the complexity of the OF. The complexity of COA is O(QC of ), where Q is the number of points in the parameter space in which the OF is calculated. Therefore, the proposed COA-based estimation has significantly lower computational complexity than evolutionary algorithms.
To show the additional advantage of COA over other techniques, we conducted a comparison in terms of required time for one iteration. In that sense, in MATLAB 2015 (MathWorks, Natick, MA, USA) we have implemented the following algorithms for solar cell parameter estimation: evaporation rate-based water cycle algorithm (ER-WCA) [19], cuckoo search (CS) [46] and harmony search (HS) [48]. ER-WCA algorithm has a very good accuracy, very close to that obtained by the proposed method (see Table 1). On the other hand, HS and CS also have a good accuracy (~10 −4 ). All computer simulations were carried out on a PC with Intel(R) Core (TM) i3-7020U CPU @ 2.30 GHz and 4 GB RAM. The obtained results, i.e., the mean, maximal and minimal required time per one iteration, obtained over 20 runs, are presented in Table 3. Clearly, the COA-based algorithm is the most efficient method, as it is characterized by the lowest value of required time per iteration. Note, in order to draw a fair comparison between the considered algorithms, MATLAB implementation follows the same rules for each algorithm (e.g., avoiding loops and using array operations such as dot product and matrix product whenever possible).  [19] 0.021063 0.024145 0.019492 CS [46] 0.029179 0.037177 0.027130 HS [48] 0.021103 0.023264 0.020393 The measured I-V and P-V characteristics and the corresponding simulated characteristics, for parameters obtained by using COA, are shown in Figure 3. Very good agreement can be seen between the measured and estimated curves. Also, the difference between the DDM and SDM simulated curves is small but consistent and always in favor of DDM. In addition, in Table 4, we presented the estimated value of the unknown DDM parameters of the BPSolar MSX-60 module. These parameters are obtained by using COA as well as by using analytical, numerical, iteration and Newton methods presented in the literature. In the COA-based estimation, the ranges of parameters were R s (Ω) ∈ [0.  The measured I-V and P-V characteristics and the corresponding simulated characteristics, for parameters obtained by using COA, are shown in Figure 3. Very good agreement can be seen between the measured and estimated curves. Also, the difference between the DDM and SDM simulated curves is small but consistent and always in favor of DDM. In addition, in Table 4, we presented the estimated value of the unknown DDM parameters of the BPSolar MSX-60 module. These parameters are obtained by using COA as well as by using analytical, numerical, iteration and Newton methods presented in the literature. In the COA-based estimation, the ranges of parameters were The measured and corresponding simulated I-V and P-V characteristics, for parameters obtained by using COA and other methods, are shown in Figures 4 and 5, respectively. It is evident that COA outperforms the other methods in terms of approaching the measured characteristics.
Based on all of the presented results, it can be concluded that COA can precisely estimate the solar cell/module circuit parameters, outperforming the other metaheuristics as well as analytical or numerical methods in terms of estimation accuracy.   The measured and corresponding simulated I-V and P-V characteristics, for parameters obtained by using COA and other methods, are shown in Figures 4 and 5, respectively. It is evident that COA outperforms the other methods in terms of approaching the measured characteristics.         Based on all of the presented results, it can be concluded that COA can precisely estimate the solar cell/module circuit parameters, outperforming the other metaheuristics as well as analytical or numerical methods in terms of estimation accuracy.

Experimental Results and Analysis
To check the applicability and efficiency of COA for solar cell parameter estimation, we also observed solar cells from the Clean Energy Trainer setup. The main motivation to use these solar modules is that this setup enables adjustable solar insolation, USB data monitoring for PC-supported data acquisition and analysis, as well as highly advanced didactic software for system control and real-time data plotting.
The observed system contains of: software-designed to facilitate system control, parameter monitoring, data acquisition and graphical representation of the collected data.
The experimental setup, installed in Laboratory for Automatics, at the Faculty of Electrical engineering, University of Montenegro, is presented in Figure 6.
3. lamp-special double spotlight lamp that simulates sunlight. It provides the optimal light spectrum for the solar module.
4. USB Data Monitor-used for data acquisition. Also, it is connected to the computer and software through the USB port.
5. load-simulates electric consumer load. 6. software-designed to facilitate system control, parameter monitoring, data acquisition and graphical representation of the collected data.
The experimental setup, installed in Laboratory for Automatics, at the Faculty of Electrical engineering, University of Montenegro, is presented in Figure 6. Firstly, we measured the I-V characteristics for insolation of 1285 W/m 2 and temperature of 42 °C. For the measured I-V pairs, we determined single and double diode solar cell parameters (see Table 5). The parameter ranges for SDM estimation were Firstly, we measured the I-V characteristics for insolation of 1285 W/m 2 and temperature of 42 • C. For the measured I-V pairs, we determined single and double diode solar cell parameters (see Table 5). The parameter ranges for SDM estimation were and I o2 (A) ∈ 5 × 10 −8 , 15×10 −8 . Then we measured the I-V and P-V characteristics for different values of insolation and temperature. The corresponding simulated characteristics were determined by taking into account the change of parameters with insolation and temperature (see [13]). The measured and estimated I-V and P-V characteristics for different values of insolation and temperature are presented in Figures 7-10. The agreement between the measured and estimated characteristics is evident (see zoomed parts in these figures). Finally, we repeated the estimation procedure on all measured I-V characteristics. The estimated values of parameters were in range of ±4% of the initially observed, which confirms that we can use any of the measured characteristics for parameter estimation. On the other hand, by observing the data provided in Table 5, even for this module, it is evident that DDM is more accurate than SDM. other hand, by observing the data provided in Table 5, even for this module, it is evident that DDM is more accurate than SDM.    other hand, by observing the data provided in Table 5, even for this module, it is evident that DDM is more accurate than SDM.    other hand, by observing the data provided in Table 5, even for this module, it is evident that DDM is more accurate than SDM.

Conclusion
Modeling of solar cells is a very popular research direction, which is supported by numerous recent contributions in the literature. This paper proposes COA as a very successful approach for this purpose.
The proposed method is verified using practical data from various manufacturers. Its accuracy is confirmed by comparing its RMSE with numerous metaheuristics and non-metaheuristics methods for different solar cells. Experimental testing of COA applicability for parameter estimation is also implemented in laboratory environment. In all considered scenarios, a high level of accuracy is demonstrated. Apart from this, excellent matching of the simulated I-V and P-V curves with the measured characteristics additionally confirms the COA accuracy and its applicability for parameter estimation.
In future work, our attention will be focused on the usage of COA for estimation of solar cell parameters when solar cell output current is represented through the Lambert W function.

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

Conclusions
Modeling of solar cells is a very popular research direction, which is supported by numerous recent contributions in the literature. This paper proposes COA as a very successful approach for this purpose.
The proposed method is verified using practical data from various manufacturers. Its accuracy is confirmed by comparing its RMSE with numerous metaheuristics and non-metaheuristics methods for different solar cells. Experimental testing of COA applicability for parameter estimation is also implemented in laboratory environment. In all considered scenarios, a high level of accuracy is demonstrated. Apart from this, excellent matching of the simulated I-V and P-V curves with the measured characteristics additionally confirms the COA accuracy and its applicability for parameter estimation.
In future work, our attention will be focused on the usage of COA for estimation of solar cell parameters when solar cell output current is represented through the Lambert W function.

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