Valuation and Determination of Seven and Five Parameters of Photovoltaic Generator by Iterative Method

Author(s)
Ahmed Yahfdhou^{1,2}^{*},
M. M. Menou^{2},
A. M. Yahya^{2},
Ne. D. Eida^{3},
A. K. Mahmoud^{2},
Issakha Youm^{1,4}

Affiliation(s)

^{1}
Laboratory of Semiconductors and Solar Energy (LASES), Faculty of Sciences and Technology, University of Cheikh Anta Diop, Dakar, Senegal.

^{2}
Applied Research Laboratory of Renewable Energy (LRAER), Faculty of Sciences and Technology, University of Sciences Technology and Medicine, Nouakchott, Mauritania.

^{3}
Unit of Electromechanical Research, Institute Superior of Technological Education, Rosso, Mauritania.

^{4}
Center of the Study and Research of the Renewable Energy (CERER), Dakar, Senegal.

ABSTRACT

The mathematical modeling of solar cells is essential for any optimization operation of the efficiency or the diagnosis of photovoltaic generator. The photovoltaic module is generally represented by an equivalent circuit whose parameters are experimentally calculated by using the characteristic current-tension, I-V. The precise determination of these parameters stays a challenge for the researchers, making to a big difference in the models and the digital methods dedicated to their characterizations. In the present paper, We are interested to characterize the parameters of single diode and two diodes models, in order to plan the behavior of the photovoltaic generator under real functioning conditions. We developed an identification method of the parameters using Newton Raphson method by using the software Matlab/Simulink. This method is the faster technique which allows the identification of several parameters and can be used in real time applications. The results of the proposed method show an accordance with the experimental and simulated characteristics of photovoltaic generator.

The mathematical modeling of solar cells is essential for any optimization operation of the efficiency or the diagnosis of photovoltaic generator. The photovoltaic module is generally represented by an equivalent circuit whose parameters are experimentally calculated by using the characteristic current-tension, I-V. The precise determination of these parameters stays a challenge for the researchers, making to a big difference in the models and the digital methods dedicated to their characterizations. In the present paper, We are interested to characterize the parameters of single diode and two diodes models, in order to plan the behavior of the photovoltaic generator under real functioning conditions. We developed an identification method of the parameters using Newton Raphson method by using the software Matlab/Simulink. This method is the faster technique which allows the identification of several parameters and can be used in real time applications. The results of the proposed method show an accordance with the experimental and simulated characteristics of photovoltaic generator.

1. Introduction

The needs of electrical energy are increasing during these last years but production constraints like pollution [1] and global warming [2] lead to the development of renewable energy sources, particularly photovoltaic energy [3] .

To surmount the problem of modeling solar panels, we need to precise the knowledge of parameters of PV cell, essentially for conception, quality control and for evaluation of their implementation. These parameters are often determined from experimental data for best weather conditions [4] [5] . The behavior of module is usually described by current-tension characteristic (I-V), the look of which depends considerably on values of the parameters such as the current of saturation, the photocurrent, the serial resistance, shunt and of the dark current [6] [7] .

The equivalent system whose parameters are calculated to achieve the characteristic I = f (V) is represented in this work. Approximatively, we offered the real model simulations, validated from data recorder in LRAER of USTM in Nouakchott, Mauritania. Subsequently, we proposed the characteristics of sunshine, in real time, to go up, the effectiveness of our method proposed for the analysis of system behavior in front of the local meteorological disturbances. We develop, at first, the model for a diode, to take account of physical phenomena with the aim of better representation of the electrical behavior of GPV from the ideal model. For precisely exam, we study the two-diode model. However, the inclusion of the additional diode increases the parameters to look and makes it difficult for extraction process of the model parameter values. This reasoning is supported by using mathematical models of digital resolution, such as Newton Raphson, expected behavior of photovoltaic generator in real time conditions in Mauritania. We are proposed methods of parameter identification to use the iterative method using the Matlab software. This method proposed in this work is a technique that can be interesting for the development of knowledge on the identification of several parameters. We arrive to show also that it provides results that can be used in applications in real time.

2. Models of the Photovoltaicarray

The solar panel is the main source of energy of the whole PV system. It is made of pv cells associated in series and parallels in addition to protection ports. To proceed to its analysis, the basic idea will be based on the study of the photovoltaic cell, the hard core of the converter of light into electricity. Several models of photovoltaic cells are proposed in the literature reviews [8] - [10] , among these mathematical models the one with single diode is generally considered as the most adapted to model a photovoltaic cell in normal functioning.

2.1. Single Diode Model

This photovoltaic cell is characterized by its equivalent diagram (Figure 1) which consists of a source of electric current which models the conversion of the luminous flow to electrical energy; a diode models the junction of the cell. To take into account physical phenomena at the level of the cell, the model is completed by two resistances series Rs and shunt R_{sh} as seen in the equivalent diagram. The serial resistance represents the accumulated resistance of, the front junction, the basic resistances and the fore and back faces. The resistance shunt result from losses by recombination’s owed essentially to the thickness; it is to reduce following penetration of metallic impurities in the junction. The equivalent circuit of a photovoltaic cell is represented by Figure 1.

This model can be expressed by the following equation, in which I and V are respectively the current and the tension of a PV cell [11] .

. (1)

It is nonlinear equation consisting of two variables (I and V) and five parameters to be determined. These parameters are:

・ I_{ph}: photo-current, equivalent current proportional to period of sunlight received by the cell.

・ I_{s}: is the reverse saturation current.

・ a: is the diode ideality factor.

・ R_{s}: is the equivalent series Resistance.

・ R_{sh}: is the parallel Resistance

2.2. Two Diodes Model

We have, this time two diodes to the present the phenomena of polarization of the junction p-n. These diodes symbolize the recombination of the minority carriers, both on the surface of the material and on its volume. The diagram of the photovoltaic generator becomes in this case that of Figure 2.

Characteristic I-V of the model with two diodes is described by this equation [12] :

(2)

where I_{s}_{1} and I_{s}_{2} are the reverse saturation currents of the first and second diode

Figure 1. Conventional single diode model.

Figure 2. Equivalent electric model of solar cells, two diodes.

respectively, V_{t}_{1} and V_{t}_{2} are the thermal voltages of the respective diodes. a_{1} and a_{2} represent the diode ideality constants.

2.3. Determination of Single Diode Model Parameters

The Equation (1) is an implicit equation there I and V which can be solved thanks to the method of Newton-Raphson or analytically with approaches, this equation has five parameters I_{ph}, I_{s}, h_{as}, R_{s} and R_{sh} are determined.

The current I_{PH} of photovoltaic cells varies according to their temperature, the period of sunlight which receives and the coefficient of temperature at the short current circuit. By making informed reference I_{PH}_{,REF} measured on the standard test condition (G_{REF} = 1000 W/m^{2} and T_{REF} = 25˚C).The current I_{PH} for a period of sunshine and a temperature given can be calculated by the following expression:

(3)

where:

G: Solar irradiance received by the module area [W/m^{2}].

α: Temperature coefficient of the short circuit current [A/K].

G_{REF}: Solar irradiance at the Standard Test Condition [1000W/m^{2}].

T_{REF}: temperature at the Standard Test Condition [25˚C].

I_{PH}_{,REF}: Light-generated at the reference condition, is practically equal to short circuit current at the STC, ICC, REF.

The expression of the diode reverse saturation current is given by:

. (4)

where V_{oc}_{,ref} open circuit voltage at the reference condition, b is the temperature coefficient of open circuit voltage and I_{cc}_{,ref} is the short circuit current at standard test condition. Constants a and b are supplied by the manufacturer of the photovoltaic cell. The serial resistances and parallel are initialized by the following equations

(5)

. (6)

2.4. Determination of Two Diodes Model Parameters

This model is model mathematically described by the Equation (2) which it expresses according to six unknown parameters. In this paper we consider that both current of saturation are equal to facilitate calculates him of running, we can write:

. (7)

This approach is to introduce by [13] by resting. The resistance shunt can have expressed by:

. (8)

V_{m} and I_{m} are the voltage and the current at maximum power point, P_{max} is power at maximum power point delivered by the module, and these data generally are supplied by the manufacturer datasheets (Table 1). The initial value of the resistance shunt is expressed by this formula:

. (9)

In this case we can choose the serial resistance of the initial value equal zero.

3. Iterative Method of Extraction of the PV Parameters

In this paper we use Newton Raphson’s method to determine the parameters of two models but also to solve Equations ((1) and (2)). This method is described by this expression [14] :

. (10)

Such as is by-product of the function f(x) = 0, x_{n} is the present value and x_{n}_{+1} is the next value. The algorithm based on the methods of Newton Raphson for the determinations of the photovoltaic parameters is presented in Figure 3.

4. Experimental System Description

The experimental device consists of a hybrid system (Figure 4) of electricity production (photovoltaic, wind, diesel and storage system) of power 5.7 kW coupled with a desalination plant of waters (reverse osmosis) and other equipment’s. The installation of the LRAER consists of the following elements:

・ A photovoltaic generator compound of 16 panels ATERSA (AP-7105/AP-75), delivers a power of 1.2 KWc,

Table 1. Electrical data of the ATERSA solar module at 25˚C, 1000W/m^{2} by the manufacturer.

Figure 3. Matching algorithm (Newton Raphson Method).

・ Two wind generators producing a power of 3 KW,

・ A system of electrochemical storage compound of 24 batteries,

・ Diesel Group,

・ A unity of data acquisitions such as the climatologically data and the parameters of the system.

To measure the radiation and the temperature we will use available sensors at the LRAER, including the following measuring devices: one thermometer: probe PT-100 a 4 wire for measuring room temperature (Figure 4(b)), a pyranometer of CMP 3 for measuring radiation in a horizontal plane (Figure 4(c)).

5. Results and Discussions

The improved iterative method technique is proposed to identify the parameters of PV module in this section. The efficiency of the improved NRM-based parameters identification method is verified by identifying the experimental data of PV module under different irradiance and temperature conditions. Comparisons with the optimization algorithm for identification are also presented for the experimental data, which are generated using the PV module models. This algorithm is programmed and implemented in MATLAB environment to identify the PV module parameters using the two models of the photovoltaic module, the parameters identified by this implementation of the algorithm are presented to Table 2.

(a)(b)(c)(d)(e)

Figure 4. Platform experimental laboratory of applied research to renewable energy in Mauritania. (a) Slides from tests and diagnostic of the PV module, (b) thermometer has probe PT-100, (c) pyrometer, (d) acquisitions of data and (e) PV system generator installed at the laboratory.

Table 2. Parameters for simulated by iterative method.

Figure 5 and Figure 6 show the I-V and P-V curves of the ATERSA photovoltaic array adjusted with the NRM. The curves exactly match with the experimental data at the three remarkable points provided by the datasheet: short current circuit, maximum power, and open circuit voltage. Table 1 shows the experimental parameters of the array obtained from the datasheet.

This method for adjusting R_{s} and R_{sh} based on the fact that there is an only pair (R_{s}, R_{sh}) that warranties that P_{max}, m = P_{max}, e = V_{m}*I_{m} at the (V_{m}, I_{m}) point of the

Figure 5. Adjusting I-V (a) and P-V (b) curves for different variations of R_{s} and R_{sh} (single diode model) at the standard test conditions (STC).

Figure 6. Adjusting I-V (a) and P-V (b) curves for different variations of R_{s} and R_{sh} (two diodes model) at the STC.

I?V curve, the maximum power calculated by the model (P_{max}, m) is equal to the maximum experimental power from the datasheet (P_{max}, e) at the maximum power point (MPP), This point on a power (I-V) curve that has the highest value of the product of its corresponding voltage and current, or the highest power output.

Figure 7 and Figure 8 represent the V?I and V-P curves for single diode model under

Figure 7. A comparison of V-I curve using the identified parameters from the NRM (solid line) and the experimental data (dots) for ATERSA (single diode).

Figure 8. A comparison of V-P curve using the identified parameters from the NRM (solid line) and the experimental data (dots) for ATERSA (single diode).

the standard test conditions of G = 1000 W/m^{2}, T = 298.2 K. Based on these curves, the necessary specifications were calculated (Table 2).

Figure 9 and Figure 10 have V-I and V-P curves for the dual diode model according to the standard test conditions, we find that there is a very good correlation between the model and experimental datasheet.

Figure 9. A comparison of V-I curve using the identified parameters from the NRM (solid line) and the experimental data (dots) for ATERSA (two diodes model).

Figure 10. A comparison of V-P curve using the identified parameters from the NRM (solid line) and the experimental data (dots) for ATERSA (two diodes model).

Figure 11 represents the profile of sunlight and temperature noted at applied research laboratory for renewable energy for three days of the month of August. We use the weather data for the validation of models based on the extraction of five and seven parameters of the PV generator, then we simulate the power delivered by the PV system during these test days.

Meteorological data show that there was remarkable perturbation on the third day caused by the transit of clouds. The maximum irradiation saved to quarter past eleven does not exceed 800 W/m^{2} while on the first day it was of the order of eight hundred eighty and for the second day it slightly exceeds nine hundred.

The variation of temperature during the three days shows that the maximum temperature is observed in the third day around 11 and a half at about 39˚C although the 1st and the second day do not exceed respectively 32 and 35 degrees Celsius.

The results in Figure 12 show a comparison between the power estimated by single diode model and the experience of our generator that consists of 16 identical panels, as well as the relative error between the simulated and the real power during three days of tests on our site.

These results obtained show a very good correlation between simulated values and experimental data at real records using the data acquisition system during the period of testing of the photovoltaic system. The relative error calculation shows that simulated model power is very close to the actual power delivered by the PV system despite the remarkable changes in weather conditions on the site.

Figure 13 shows a comparison of the power obtained by the model of two diodes and the power of real or experimental PVG during three days of the test, and also the variation of relative error between these powers.

Figure 11. The sunshine (a) and temperature (b) profiles recorded for three days at experimental site research center applied to renewable energy.

Figure 12. Comparison between the simulated power and experimental raised power (a) for three days in the case of single diode model, (b) relative error.

Figure 13. Comparison between the simulated power and experimental raised power (a) for three days in the case of single diode model, (b) relative error.

The relative error of power generated by the two models (see Table 3) shows that the model with two diodes is more adapted to the behavior forecast of the photovoltaic system.

Table 3. Relative power of error for both models.

6. Conclusion

In this paper, the seven and five parameters of a solar PV module are extracted using an iterative method based on Newton Raphson Method (NRM). The feasibility of the NR method has been validated by synthetic and experimental I?V and P-V data set of solar module (mono-crystalline). The proposed method is very accurate and converges to the rapidly solution. The results also presented good agreement between modeling data and outdoor measurements of photovoltaic generator for single and two-diode models, even in low irradiance levels. The mean absolute relative error between estimated power and measurements is less than 1.1% for single diode model parameters and less than 1% for two-diode model.

Cite this paper

Yahfdhou, A. , Menou, M. , Yahya, A. , Eida, N. , Mahmoud, A. and Youm, I. (2016) Valuation and Determination of Seven and Five Parameters of Photovoltaic Generator by Iterative Method.*Smart Grid and Renewable Energy*, **7**, 247-260. doi: 10.4236/sgre.2016.79019.

Yahfdhou, A. , Menou, M. , Yahya, A. , Eida, N. , Mahmoud, A. and Youm, I. (2016) Valuation and Determination of Seven and Five Parameters of Photovoltaic Generator by Iterative Method.

References

[1] Askarzadeh, A. and Razazadeh, A. (2013) Extraction of Maximum Power Point in Solar Cells Using Bird Mating Optimizer-Based Parameters Identification Approach. Solar Energy, 90, 123-133.

http://dx.doi.org/10.1016/j.solener.2013.01.010

[2] Alrashidi, M.R., Alhajri, M.F., El-Naggar, K.M. and Al-Othman, A.K. (2011) A New Estimation Approach for Determining the I-V Characteristics of Solar Cells. Solar Energy, 85, 1543-1550.

http://dx.doi.org/10.1016/j.solener.2011.04.013

[3] Seifi, M., Chesoh, A.B., Abdwahab, N.I. and Hasan, N.I. (2013) A Comparative Study of PV Models in Matlab/Simulink. Word Academy of Science, Engineering and Technology, 74, 108-113.

[4] Celik, A.N. and Acikgoz., N. (2007) Modeling and Experimental Verification of the Operating Current of Mono-Crystalline Photovoltaic Modules Using Four and Five Parameter Models. Applied Energy, 84, 1-5.

http://dx.doi.org/10.1016/j.apenergy.2006.04.007

[5] Zadeh, A.A. and Rezazadeh, A. (2013) Artificial Bee Swarm Optimization Algorithm for Parameters Identifications of Solar Cell Modules. Applied Energy, 102, 943-949.

http://dx.doi.org/10.1016/j.apenergy.2012.09.052

[6] Sandrolini, L., Artioli, M. and Reggiani, U. (2010) Numerical Method for Extraction of Photovoltaic Modules Double-Diode Model Parameters through Cluster Analysis. Applied Energy, 87, 442-451.

http://dx.doi.org/10.1016/j.apenergy.2009.07.022

[7] Dakichi, F., Ouakarfi, B., Fakkar, A. and Belbounaguia, N. (2014) Parameters Identification of Solar Cell Model Using Levenberg-Marquardt Algorithm Combined with Simulated Annealing. Solar Energy, 110, 781-788.

http://dx.doi.org/10.1016/j.solener.2014.09.033

[8] Bonkoungou, D., Koalaga, Z. and Njomo, D. (2013) Modeling and Simulation of Photovoltaic Module Considering Single-Diode Equivalent Circuit Model in Matlab. International Journal of Emerging Technology and Advanced Engineering, 3, 493-502.

[9] Yahya, M., Youm, I. and Kader, A. (2011) Behavior and Performance of a Photovoltaic Generator in Real Time. International Journal of the Physical Science, 6, 4361-4367.

[10] Ishaque, K., Salam, Z. and Syafaruddin (2011) A Comprehensive MATLAB Simulink PV System Simulator with Partial Shading Capability Based on Two-Diode Model. Solar Energy, 85, 2217-2227.

http://dx.doi.org/10.1016/j.solener.2011.06.008

[11] Yahfdhou, A., Mahmoud, A. and Youm, I. (2013) Modeling and Optimization of Photovoltaic Generator with Matlab/Simulink. International Journal of I Tech and E Engineering, 3, 108-111.

[12] Wolf, P. and Benda, V. (2013) Identification of PV Solar Cells and Modules Parameters by Combining Statistical and Analytical Methods. Solar Energy, 93, 151-157.

http://dx.doi.org/10.1016/j.solener.2013.03.018

[13] Ishaque, K., Salam, Z. and Thateri, H. (2011) Simple Fast and Accurate Two Diode Model for Photovoltaic Modules. Solar Energy Materials & Solar Cells, 95, 586-594.

http://dx.doi.org/10.1016/j.solmat.2010.09.023

[14] Bonkoungou, D., Koalaga, Z. and Njomo, D. (2013) Modeling and Simulation of Photovoltaic Module Considering Single-Diode Equivalent Circuit Model in Matlab. International Journal of Emerging Technology and Advanced Engineering, 3, 493-502.

[1] Askarzadeh, A. and Razazadeh, A. (2013) Extraction of Maximum Power Point in Solar Cells Using Bird Mating Optimizer-Based Parameters Identification Approach. Solar Energy, 90, 123-133.

http://dx.doi.org/10.1016/j.solener.2013.01.010

[2] Alrashidi, M.R., Alhajri, M.F., El-Naggar, K.M. and Al-Othman, A.K. (2011) A New Estimation Approach for Determining the I-V Characteristics of Solar Cells. Solar Energy, 85, 1543-1550.

http://dx.doi.org/10.1016/j.solener.2011.04.013

[3] Seifi, M., Chesoh, A.B., Abdwahab, N.I. and Hasan, N.I. (2013) A Comparative Study of PV Models in Matlab/Simulink. Word Academy of Science, Engineering and Technology, 74, 108-113.

[4] Celik, A.N. and Acikgoz., N. (2007) Modeling and Experimental Verification of the Operating Current of Mono-Crystalline Photovoltaic Modules Using Four and Five Parameter Models. Applied Energy, 84, 1-5.

http://dx.doi.org/10.1016/j.apenergy.2006.04.007

[5] Zadeh, A.A. and Rezazadeh, A. (2013) Artificial Bee Swarm Optimization Algorithm for Parameters Identifications of Solar Cell Modules. Applied Energy, 102, 943-949.

http://dx.doi.org/10.1016/j.apenergy.2012.09.052

[6] Sandrolini, L., Artioli, M. and Reggiani, U. (2010) Numerical Method for Extraction of Photovoltaic Modules Double-Diode Model Parameters through Cluster Analysis. Applied Energy, 87, 442-451.

http://dx.doi.org/10.1016/j.apenergy.2009.07.022

[7] Dakichi, F., Ouakarfi, B., Fakkar, A. and Belbounaguia, N. (2014) Parameters Identification of Solar Cell Model Using Levenberg-Marquardt Algorithm Combined with Simulated Annealing. Solar Energy, 110, 781-788.

http://dx.doi.org/10.1016/j.solener.2014.09.033

[8] Bonkoungou, D., Koalaga, Z. and Njomo, D. (2013) Modeling and Simulation of Photovoltaic Module Considering Single-Diode Equivalent Circuit Model in Matlab. International Journal of Emerging Technology and Advanced Engineering, 3, 493-502.

[9] Yahya, M., Youm, I. and Kader, A. (2011) Behavior and Performance of a Photovoltaic Generator in Real Time. International Journal of the Physical Science, 6, 4361-4367.

[10] Ishaque, K., Salam, Z. and Syafaruddin (2011) A Comprehensive MATLAB Simulink PV System Simulator with Partial Shading Capability Based on Two-Diode Model. Solar Energy, 85, 2217-2227.

http://dx.doi.org/10.1016/j.solener.2011.06.008

[11] Yahfdhou, A., Mahmoud, A. and Youm, I. (2013) Modeling and Optimization of Photovoltaic Generator with Matlab/Simulink. International Journal of I Tech and E Engineering, 3, 108-111.

[12] Wolf, P. and Benda, V. (2013) Identification of PV Solar Cells and Modules Parameters by Combining Statistical and Analytical Methods. Solar Energy, 93, 151-157.

http://dx.doi.org/10.1016/j.solener.2013.03.018

[13] Ishaque, K., Salam, Z. and Thateri, H. (2011) Simple Fast and Accurate Two Diode Model for Photovoltaic Modules. Solar Energy Materials & Solar Cells, 95, 586-594.

http://dx.doi.org/10.1016/j.solmat.2010.09.023

[14] Bonkoungou, D., Koalaga, Z. and Njomo, D. (2013) Modeling and Simulation of Photovoltaic Module Considering Single-Diode Equivalent Circuit Model in Matlab. International Journal of Emerging Technology and Advanced Engineering, 3, 493-502.