Enhanced Vibrating Particles System Algorithm for Parameters Estimation of Photovoltaic System ()
1. Introduction
Many disadvantages affect the availability of fossil fuels: the fluctuating prices, the environmental pollution and the fact that they are not abundant [1]. The Energy produced by the sun is the most widespread, free and clean among all renewable energy resources. In recent years, interest to use PV as power generation has increased because of its many advantages [2]. Photovoltaic solar installations around the world are down from 89.5 GW in 2012 to just over 303 GW in 2016 [3]. The PV designers need reliable and accurate tools to predict the power produced by a PV [4]. A PV array comprises several photovoltaic cells connected in series and parallel according to the output power desired. Cells are made from semiconductor materials that produce an electric current when illuminated; the intensity of the current depends on the quantity of solar irradiance [5] [6]. Many factors like solar radiation [7], location latitude influence the output power of the solar system [8].
To design and assess the operation of a PV system, a PV model should be implemented with appropriate accurateness that one can employ to predict the reliable I-V and P-V output characteristics under normal operation [9]. To do this, many models have been developed in the literature. Several parameters need to be accurately extracted with good precision for the purpose to evaluate the performance of a PV system. These intrinsic parameters are: saturation current, series resistance, diode ideality factor saturation current, generated photocurrent and shunt resistance. As output power is proportional to solar irradiance, an estimate of the intrinsic parameters of the PV is necessary in order to evaluate its performance [10]. To extract these intrinsic parameters, we can use either the manufacturer’s datasheet or experimentally measure the voltage and current from the PV [11].
Many methods in the literature have been developing to extract PV parameters. These methods can be classified into three categories: analytical methods, numerical methods and evolutionary methods. In the analytical method, a set of transcendental equations is solved to extract parameters from solar cell [12]. The main advantage of the analytical method is the speed of calculation and reasonably accurate results. Analytical methods are simple. They have a reduced calculation time. Sometimes, just one iteration is necessary to reach the result [10].
Explicit modeling from current and voltage characteristic is used by [13] where from a single diode model, a Pade’s approximate method is used to extract the parameters. In [14], Lambert W-function is used to extract parameters. In [15], the author used Taylor’s series expansion to extract the five parameters in the single diode model. [16] developed a sample model to extract just four parameters without shunt resistance. In [17], analytical methods are compared with a curve-fitting tool method, and the result shows that the analytic method is more accurate. Analytical methods work properly under standard weather conditions; but when weather conditions change, analytical methods become ineffective [18].
Numerical extraction techniques based on some algorithm fit the points on the PV characteristic curve. Compared to the analytical method, an accurate result can be attained since the algorithm tries to consider all points on the characteristic curve [19]. In the literature, the Newton-Raphson method is the most used [19] [20] [21] ; In [22], a numerical method is proposed for modeling and the simulation of PV. The method finds the five parameters from the current-voltage characteristic by using three points of the curve (maximum power, open circuit and short circuit). In [23], the Levenberg Marquardt algorithm is implemented to estimate the parameters from the I-V characteristic of a single diode PV model.
The main drawbacks of numerical techniques such as Newton Raphson are the need for extensive computations for convergence and fail to result in accurate results when the number of parameters to be estimated increases and a close approximation of initial conditions [24].
Despite the efficiency of the numerical methods, their slow convergence does not always guarantee the best result because they can converge through a local minimum and the choice of the initial condition is not often easy [11].
To overcome the drawback of analytical and iterative methods, metaheuristic algorithms have been developed. They are nature-inspired algorithms using probabilities to find the best result. They have shown their effectiveness in solving difficult problems. Their main advantage is that they do not need continuity and differentiability of the objective function In the last decade, metaheuristics have been frequently applied for parameter estimation of circuit model parameters of solar PV cells. The main develops in recent research are: genetic algorithm (GA) [25], grey wolf optimization (GWO) [26], particles swarm optimization (PSO) [27], moth-flame optimization algorithm (MFOA) [28], harmony search (HS) [29], artificial neural network (ANN) [30], multi-verse optimizer (MVO) [31], bond-graph based modelling [32], cuckoo search (CS) [33], bacterial foraging optimization [34], multiple learning backtracking search algorithm (MLBSA) [35], whale optimization algorithm (WAO) [36], salp swarm-inspired algorithm (SSA) [37] … New metaheuristic algorithms have been also recently developed to solve mathematic and engineering problems. [38] used World Cup Optimization (WCO) algorithm to find the optimal parameters of PID controller; in [39] a new algorithm based on Variance Reduction of Guassian Distribution is proposed; a new algorithm based on the invasive weed by the quantum computing is proposed by [40] ; [41] combined Gravitational Search Algorithm (GSA) and Particle Swarm Optimization (PSO) to train wavelet neural networks.
Until today in the research and industry domain, there is no method of extracting PV parameter that has been introduced in the manufacturing of PV. By another, “The no-free-lunch theorem” remarked that: there is no algorithm able to solve all optimization problems, where it is important to propose new algorithms for solving engineering optimization problems [42]. In this paper, a new algorithm based on EVPS is used to extract parameters from a PV system. The inspiration of vibrating particles system comes from the free vibration of one degree of freedom systems with viscous damping. The algorithm has been applied to extract the best parameters of a PV cell and module under uniform and partial shading conditions. Five recent algorithms (SSA [37], GWO, MFOA [28], WAO [36], MVO [31] ) are also implemented on the same computer with the parameters gave by authors. The result obtained from the EVPS is compared with other recent methods in the literature and different results obtain to demonstrate the high quality of the algorithm.
The rest of this paper is presented as follows: In Section 2, PV cell modeling is presented; Section 3 presents the problem formulation for extracting parameters from single and double diode model; the inspiration and the mathematical model of vibrating particles system are proposed in Section 4; Section 5 presents the different case study with different results and Section 6 is the conclusion.
2. Photovoltaic Cell Modeling
Many models of PV cell have been developed in the literature; but there are two models mostly used: single and double diode model.
2.1. Single Diode Model
Most of the literature uses this model. The main reason for this widely used is their simplicity and the least number of parameters. In the single diode model, there are five parameters to be extracted. Figure 1 shows the electric diagram of the model.
The current I at the output of a PV module can be expressed using Kirchhoff’s theorem in Equation (1).
(1)
The diode current can be express as:
(2)
where
(3)
The current through parallel resistance is
Figure 1. Electrical diagram of one diode PV cell.
(4)
By replacing Equation (2) & Equation (3) into Equation (1), we have the output current at the output of a PV module.
(5)
The parameters which characterize this equation are:
. These five parameters can be determined by all the method described in Section 1.
2.2. Double Diode Model
This double-diode model (Figure 2) has a better accuracy than one diode model, but also more complex because of the numbers of parameters [11]. The model has been used by many authors [11] [43] [44] [45].
The current I, at the output of a PV module can be expressed using Kirchhoff’s theorem in Equation (6).
(6)
The diodes currents can be express as:
(7)
(8)
where
(9)
The current through parallel resistance is
(10)
By replacing Equations (7)-(10) into Equation (6), we have the output current at the output of a PV module:
Figure 2. Electrical diagram of double diode PV cell.
(11)
The parameters which characterize this equation are:
.
These seven parameters can be determined by all the methods described in Section 1.
3. Problem Formulation
To extract the PV parameters using EVPS, the problem can be set as an optimization problem with the objective to minimize the difference between measured and estimated current. The objective function (OF) is defined as the root mean square error (RMSE) where the error function is defined as the difference between estimated and experimental currents. It’s expressed as follows:
(12)
where:
: is the objective function to minimize;
N: is the number of points (Vi, Ii) measured;
: is the measured current;
: is the estimated current;
: Parameters to estimate.
For a single diode model, the objective function is expressed as:
(13)
with
the parameters to estimate.
For double diode model, the objective function is:
(14)
with
the parameters to estimate
In this paper, EVPS algorithm is used to minimize Equation (13) and Equation (14).
4. Optimization Algorithm
4.1. Inspiration
The inspiration of vibrating particles system comes from the free vibration of a single degree of freedom system with viscous damping. The VPS contains a number of population solutions that represent the particle system. The particles are randomly initialized in an n-dimensional search space and step-by-step, they approach their equilibrium positions [46].
4.2. The Vibrating Particles System Algorithm
As other meta-heuristic algorithms, VPS has population particles which are considered as the parameters of the problem. The initial positions of particles are firstly generated randomly.
(15)
represents the jth position of the ith particle;
and
are respectively the initial and the final position
is a random number between [0, 1].
Three equilibrium positions affected by different weights are defined for each particle. During each generation, the particle positions are updated by learning from them. The equilibrium positions are:
- HB: Historically best location;
- GP: Good particle;
- BP: Bad particle.
To include the effect of the damping level in the vibration, a descending function is introduced:
(16)
and
: represent respectively the current iteration, the maximum iteration and a constant.
The next position is updated by the following equations:
(17)
: represent the jth position of the ith particle;
: parameters to measure the best value of HB, GP, BP;
: random numbers between [0, 1].
4.3. Enhanced Vibrating Particles System Algorithm
4.3.1. Description
EVPS algorithm has been initially developed by [47] ; the main advantages of EVPS are to avoid slow convergence, local minimum and increase the number of space search.
In EVPS, we introduced two new parameters: “memory” and “OHB (one of the best historically locations in the whole population)” [47]. HB in VPS is replaced by the memory. The memory now saves the best historically positions of the whole population. OHB is one row of memory whose selection is random. The next changing is the replacement of Equation (17) by Equation (18).
(18)
(±1) are applied randomly. BP, GP and OHB are independently determined for each particle [47]. All details concerning VPS algorithm can be found in [47].
4.3.2. Details Algorithms
The steps to compute EVPS algorithms are described as follows:
Step 1: Initializing of EVPS’s parameters.
- Initialize VPS parameters (size of the population, number of optimization variables, memory size, maximum number of iterations, lower and upper bound of the variables, parameters for handling the side constraints, w1 and w2).
- Initializing particles positions using Equation (19).
(19)
Step 2: Search.
- Evaluate the objective function for each particle.
- For each particle, select “memory” and “OHB (one of the best historically locations in the whole population).
- Update particle’s position.
(20)
Updating the next position by
(21)
A parameter like k in a range of [0, 1] must be defined to specify if BP must be considered in the new position. For each population, k is compared with a random number (rand) uniformly distributed in the range of [0, 1]; if k < rand, then w3 = 0 and w1 = 1 − w2.
Step 3: Handling the Side Constraints.
If a particle went out of the boundary, it must be updated by harmony search-based side constraints handling approach. The method consists to determine if the violating particle should be updated either by the best historically particle or randomly in the search space.
Step 4: Out memory and best positions.
4.4. Application of Enhanced Vibrating Particles System to Photovoltaic Parameters Extraction
4.4.1. Problem Statement
The goal is to determine the global optimum, which is the best value of the OF (RMSE). If we replace memory (OHB) by the global optimum and the positions (
) by the estimated parameters, therefore (
), the best position automatically moves towards it. However, the problem is that the global solution of the optimization problems is unknown. In this case, the optimal solution obtained is the global optimum and presumed as the best selection of the memory.
4.4.2. Problem Formulation
The objective function of Equation (22) is used to find the best parameters of single and double diode model; the difference between single and double diodes model is the numbers of parameters (5 to single diode model and 7 to double diodes model). The formulation is:
(22)
where
is the best parameters; ng is the number of parameters;
the RMSE; N the number of points
measured;
the measured current;
the estimated current;
is the lower bound, and
the upper bounds;
the design constraints and Nc the number of constraints.
The objective function of each model is formulated as follows:
- Single diode model
For a single diode model, the objective function is expressed as:
(23)
with
the five estimated parameters which correspond respectively to
.
- Double diode model
For single diode model, the objective function is expressed as:
(24)
with
the seven estimated parameters which correspond respectively to
.
The flowchart algorithm is represented in Figure 3.
Figure 3. Flowchart of the EVPS algorithm.
5. Experiment and Results
This section presents different results and implementation of the algorithm. Four case studies have been implemented. The first two cases of study have been implemented in Matlab 2017a. In the first case study; Photowatt-PWP201 PV which has 36 polycrystalline silicon cells, all connected in series, the irradiance is 1000 W/m2 and temperature 45˚C. The second case refers to the RTC France commercial silicon PV, irradiance is 1000 W/m2 and temperature 33˚C. These two cases were for the first time initiated by [48] and it’s largely used today in research as a test system. To show the performance of the algorithm, a Matlab-Simulink model has been implemented in different irradiance conditions. The third case refers to polycrystalline SW255. The fourth case implements real experimental data from the Sharp ND-R250A5 PV module provided to us by [11].
5.1. Case Study 1
A single diode model has been implemented in this case study to extract the five parameters of the Photowatt-PWP201 PV which is a 11.5 W PV module. It has 36 cells connected in series. Irradiance is 1000 W/m2 and temperature 45˚C. It’s widely used in literature by many researchers. The manufacturer’s characteristics of the PV module at STC are listed in Table 1 and the lower and upper bound are expressed in Table 2. The 26 I-V measured data has been collected from [35] [11]. Five recent algorithms (SSA, GWO, MFOA, WAO and MVO) are also implemented in this case study in order to make the comparison. The initial parameters of each algorithm are given by the authors in each paper.
The initial parameters of EVPS are: number of search agents = 50; maximum number of iterations = 1000; alpha = 0.05; w1 = 0.3; w2 = 0.3.
5.2. Result of Case Study 1
Table 3 presents the result of the five estimated parameters of the Photowatt-PWP201 PV module and the best OF (RMSE) after 20 independent tests is 2.4267 × 10−3. In Table 4 the results of the 5 other algorithms are presented. The Average value of RMSE shows the constant of the algorithm after many tests. The comparison with other methods in the literature is presented in Table 5 to
Table 1. PWP201 PV module parameter.
Table 3. Best result for case study 1.
Table 4. Comparison results of five algorithms for 20 independent tests.
Table 5. Comparison of the results obtained from the Photowatt-PWP201 PV with other methods in the literature.
show the superiority of the algorithm. Figure 4 shows the I-V characteristics of the measured and estimated curve of the Photowatt-PWP201 PV under a 1000 W/m2 irradiance and 45˚C temperature. In Figure 5, the P-V characteristic of the measured and estimated curve is presented. Figure 6 presents the convergence of each algorithm.
5.3. Case Study 2
In this case, seven parameters of RTC France PV have been extracted. The irradiance of the RTC France PV is 1000 W/m2 and temperature 33˚C. The typical electrical characteristics of the PV cell at STC are listed in Table 6; the lower and upper bound are expressed in Table 7. The 26 I-V measured data have been collected from [11]. The initial parameters of EVPS are the same as in case study 1.
5.4. Result of Case Study 2
Table 8 presents the result of the seven estimated parameters of the PV cell; the best OF (RMSE) after 20 tests is 9.8510e−4. In Table 9, the results of the 5 other algorithms are presented. In Table 10, the comparison with other methods in the literature is presented to show the superiority of the algorithm. Figure 7 shows the I-V characteristic of the measured and estimated curve of the PV cell under a 1000 W/m2 and 33˚C. In Figure 8, the P-V characteristic of the measured and estimated curve is presented.
5.5. Case Study 3
The case study 3 consists of the implementation of a Matlab/Simulink model at
Figure 4. I-V characteristic of the measured and estimated curve.
Figure 5. P-V characteristic of the measured and estimated curve.
Table 6. RTC France cell parameter.
Figure 6. Comparaison of convergence curve.
Figure 7. I-V characteristic of the measured and estimated curve.
Figure 8. P-V characteristic of the measured and estimated curve.
Table 9. Comparison results of five algorithms from the RTC France for 20 independent tests.
Table 10. Comparison of the RTC France PV obtained results with the ones in the literature.
different irradiance conditions. The case study 3 refers to the polycrystalline SW255. The using manufacturer data at STC is reported in Table 11 [54]. The Simulink model is presented in Figure 9. The block solar cell is configured with the STC condition of Table 11. The experiment data at each irradiance condition (E(W/m2) = 1000, 800, 600, 400, 200) at 25˚C temperature are exported to Matlab via the blocks Workspace.
5.6. Result of Case Study 3
Table 12 and Table 13 present respectively the result of the five and seven estimated parameters of the polycrystalline SW255 module at different irradiance conditions. In Table 14, the results of the 5 other algorithms are presented. Figure 10 shows the I-V characteristic of the measured and estimated curve under different irradiance. In Figure 11, the P-V characteristic of the measured and estimated curve is presented.
5.7. Case Study 4
One diode model has been implemented in this last case to extract the five parameters based on experimental data of the Sharp ND-R250A5 PV module. The PV has 60 cells in series. Irradiance is 1040 W/m2 and temperature 59˚C. The typical electrical characteristics of the Sharp ND-R250A5 PV module at STC are listed in Table 15 and the lower and upper bound are expressed in Table 16. The 36 I-V measured data has been provided to us by [11], where irradiance and temperature have been measured by the sensor Ingenieurbüro Si-13TC-T (Figure 12).
5.8. Result of Case Study 4
Table 17 presents the result of the five estimated parameters of the Sharp ND-R250A5 PV module. The best OF (RMSE) is 11.252719 × 10−3. Figure 13 shows the I-V characteristics of the measured and estimated curve of the ND-R250A5 PV under a 1040 W/m2 irradiance and 59˚C temperature. In Figure 14, the P-V characteristic of the measured and estimated curve is presented. Figure 15 presents the convergence curve.
Table 11. Datasheet parameters of SW255 at STC.
Table 12. Five estimated parameters of the polycrystalline SW255 module at different irradiance conditions.
Table 13. Seven estimated parameters of the polycrystalline SW255 module at different irradiance conditions.
Table 14. Comparison results of five algorithms for 20 independent tests.
Table 15. Sharp ND-R250A5 module parameter.
Figure 9. Simulink model of polycrystalline SW255.
Figure 10. I-V characteristic of the measured and estimated curve.
Figure 11. P-V characteristic of the measured and estimated curve.
Figure 13. I-V characteristic of the measured and estimated curve.
Figure 14. P-V characteristic of the measured and estimated curve.
Figure 15. Convergence curve during the parameter extraction for the one diode model ND_R250A5 solar cell.
6. Conclusions
In this paper, we have presented a novel bio-inspired optimizer of a very recent heuristic-based on technique, namely enhanced vibrating particles system to extract the best values of parameters of a photovoltaic cell. The particles are randomly initialized in an n-dimensional search space and Step-by-Step, they approach their equilibrium positions.
To show the performance of the algorithm, many cases have been implemented from one and two diodes model. The current-voltage and power-voltage characteristic of measured and estimated data show the best accuracy of the method. The simulations result and comparisons with another method exhibit high accuracy and validity of the proposed Enhanced Vibrating particles system to extract parameters of a photovoltaic cell and module. Thus, enhanced vibrating particles system can be recommended as an efficient method not only to extract the best parameters of a PV cell and module, but also to solve optimization problems in power systems. As every algorithm, enhanced vibrating particles system has some drawback like the variability of the result at each independent test and the limit of the algorithm to solve the only mono-objective problem. In the future work, the stability of the enhanced vibrating particles system should be improved and other parameters should be added to permit the algorithm to solve multi-objective optimization in power systems.
Nomenclature
E: Solar irradiance;
: objective function to minimize;
I: cell output current [A];
: is the estimated current;
: measured current [A];
: estimated current [A];
: Diode reverse saturation currents [μA];
: diode currents [A];
current at the maximum power point [A];
photoelectric current [A];
short-circuit current [A];
k: Boltzman constant [J/K];
: temperature coefficient of Isc [A/K];
: Diode ideality factors;
N: number of the experimental I-V data pairs;
: number of cells connected in series;
OF: objective function;
q: electron charge [C];
: random numbers between [0, 1];
: parallel resistance [Ω];
: series resistance [Ω];
RMSE: root mean square error;
STC: Standard testing condition (1000 watts/m2, 25˚C);
T: Temperature [K];
V: cell output voltage [V];
: measured voltage [V];
: voltage at the maximum power point [V];
: open-circuit voltage [V];
: thermal voltage [V];
: parameters;
: position of the particle;
: Parameters to estimate.