# A New Strategy for Accurately PredictingI-V Electrical Characteristics of PV Modules Using a Nonlinear Five-Point Model

Authors: Sakaros Bogning  Dongue; Donatien Njomo; Lessly Ebengai
Journal: Journal of Energy
 (2013)

---

## Abstract

This paper presents the modelling of electricalI-V response of illuminated photovoltaic crystalline modules. As an alternative method to the linear five-parameter model, our strategy uses advantages of a nonlinear analytical five-point model to take into account the effects of nonlinear variations of current with respect to solar irradiance and of voltage with respect to cells temperature. We succeeded in this work to predict with great accuracy the I-V characteristics of monocrystalline shell SP75 and polycrystalline GESOLAR GE-P70 photovoltaic modules. The good comparison of our calculated results to experimental data provided by the modules manufacturers makes it possible to appreciate the contribution of taking into account the nonlinear effect of operating conditions data on I-V characteristics of photovoltaic modules.


---

## Body

## 1. Introduction

The performance of photovoltaic modules is highly dependent on both the availability of solar radiation in the locality and the temperature of photovoltaic cells [1]. A real knowledge and mastery of photovoltaic modules are important as are the design and technology choice for a project operating photovoltaic electricity.Important research has been carried out throughout the world in view of mastering electrical performance of photovoltaic modules [2, 3]. Numerical models are developed for this purpose [4, 5], mainly to reduce the costs of experimental studies. Unfortunately, these studies are very useful, particularly to obtain accurate results [6, 7]. It is obvious that the application of the best model will combine in a better way simplicity and accuracy [8, 9].This work consists in doing comparative study between a five-parameter three-point model hereafter referred to as five-parameter model characterized by both the linear variation of the current with respect to the irradiance and the linear variation of the voltage with respect to the temperature of the cells and an adjusted five-parameter five-point model hereafter referred to as five-point model which takes into account the nonlinear effects of both the variation of the current with respect to the irradiance and the variation of the voltage with temperature.


## 2. Single-Diode Equivalent Circuit of a PV Module

A photovoltaic module under illumination can be represented either by two-diode [10] or by one-diode equivalent circuit. In normal operating conditions, the one-diode model is the most common [11], and the photovoltaic module is equivalent to the following electrical circuit (see Figure 1).Figure 1
Equivalent electrical circuit for the one-diode model.The current-voltage (I-V) characteristic under a solar radiation and temperature set is expressed in
(1)I=Iph-Id-Ish,
where Id is the current through the diode and Ish is the current through the shunt resistor both given by (2) as follows:
(2)Id=I0[exp(V+IRSa)-1],Ish=V+IRSRsh.
Therefore,
(3)I=Iph-I0[exp(V+IRSa)-1]-V+IRSRsh,
where Iph, I0, Rs, and Rsh are the photocurrent, the reverse saturation current of the diode, the series resistance and shunt resistance, respectively. The modified ideality factor of the junction is given as
(4)a=NsA0KTcq,
where A0 is the normal ideality factor; Ns, k, Tc, and q are the number of cells in series, the Boltzmann constant, the temperature of the cells, and elementary charge, respectively.The determination of the electrical behaviour of a photovoltaic module in given operating conditions is based on the current-voltage (I-V) curve. However, (1) that gives this characteristic is implicit and cannot be easily solved. Many numerical algorithms to model the I-V curve from one-diode representation of photovoltaic module have been proposed. Some of these models are known as “resistive companion” [12], the nonlinear least squares optimization, and other iterative solutions described in [8, 9, 13]. Although many improvements have been done, these approaches still require a great effort of numerical calculations. Recently, some authors have used artificial intelligence methods (AI) such as fuzzy logic [14] and artificial neural network (ANN) [15, 16] for modelling (I-V) curves.Another interesting approach is based on the evaluation of unknown parameters of (1) based on information related to the field operation of the module in order to obtain an explicit I-V relationship. To this end, an iterative process may be necessary to evaluate the different parameters [17, 18]. One may also perform a direct analytical calculation of each parameter [19, 20]. The latter approach requires information often available in the data sheets provided by photovoltaic modules manufacturers or experimental data collected on site. Both five-parameter and five-point models described in this study use only data provided by module manufacturers, together with semiempirical correlation equations and weather data, to predict I-V characteristics of photovoltaic crystalline modules.


## 3. Five-Parameter Model

The five parameters that give name to the model are the photocurrent, the reverse saturation current, series resistance, shunt resistance, and modified ideality factor. These parameters are firstly evaluated at so-called reference operating conditions or standard test conditions (STC) for which the irradiance on the modules is 1000 W/m² and the cells temperature is 25°C with an air mass number AM = 1.5. This environmental situation is used by manufacturers of the modules to measure the data laid out on the data sheet.The five parameters defined above corresponding to STC are designated byIphref, I0ref, Rsref, Rshref, and aref (the index is relative to reference conditions). We obtain the (I-V) characteristic equation at STC as
(5)I=Iphref-I0ref[exp(V+IRsrefaref)-1]-V+IRsrefRshref.


### 3.1. Calculation of Reference Parameters

Data provided by module manufacturers for STC are used when writing the following five equations [21].At short-circuit point: V=0, I=Isc,
(6)Isc0=Iphref-I0ref[exp(Isc0Rsrefaref)-1]-Isc0RsrefRshref.
At the open-circuit point: I=0, V=Voc,
(7)0=Iphref-I0ref[exp(Vocrefaref)-1]-VocrefRshref.
At maximum power point (MPP): I=Imp, V=Vmp,
(8)Impref=Iphref-I0ref[exp(Vmpref+ImprefRsrefRshref)-1]-Vmpref+ImprefRsrefRshref.
At MPP, the derivative of the power with respect to the voltage is equal to zero, that is,
(9)[d(IV)dV]mpp,ref=0.
At the short-circuit point, the slope of the I-V curve at short circuit is given by
(10)(dIdV)Isc=-1Rshref.
A simultaneous resolution of (6)–(10) yields the five reference parameters, and then the I-V curve for STC is built using (5). As one can notice, in this strategy we made use of only three critical points of the I-V curve, namely, the short-circuit point, the maximum power point, and the open-circuit point.


### 3.2. Dependence of Operating Conditions

The determination of theI-V characteristic in all other operating conditions can be made possible by evaluating the various parameters of the I-V equation in the new weather environment using the semiempirical relationships that follow as shown in [21]. One may consider the following:
(11)aaref=TcTcref,I0I0ref=[TcTcref]3exp[EGNsaref(1-TcTcref)],Iph=GGref[Iph,ref+αcc(Tc-Tcref)],GGref=[-(I0/a)exp[(Vmp+ImpRs)/a]-1/Rsh1+(I0Rs/a)exp[(Vmp+ImpRs)/a]+Rs/Rsh]×[[1+(I0refRsrefaref)exp[(Vmpref+ImprefRsref)aref]+RsrefRshref]-1-I0refarefexp[(Vmpref+ImprefRsref)aref]-1Rshref/×(1+(I0refRsrefaref)exp[(Vmpref+ImprefRsref)aref]+RsrefRshref)-1]-1,Rsh=Rsh,ref.
Finally at the MPP, we write the following two equations:
(12)Imp=ImprefGGref,Vmp=Vmpref+βVoc(Tc-Tcref).


## 3.1. Calculation of Reference Parameters

Data provided by module manufacturers for STC are used when writing the following five equations [21].At short-circuit point: V=0, I=Isc,
(6)Isc0=Iphref-I0ref[exp(Isc0Rsrefaref)-1]-Isc0RsrefRshref.
At the open-circuit point: I=0, V=Voc,
(7)0=Iphref-I0ref[exp(Vocrefaref)-1]-VocrefRshref.
At maximum power point (MPP): I=Imp, V=Vmp,
(8)Impref=Iphref-I0ref[exp(Vmpref+ImprefRsrefRshref)-1]-Vmpref+ImprefRsrefRshref.
At MPP, the derivative of the power with respect to the voltage is equal to zero, that is,
(9)[d(IV)dV]mpp,ref=0.
At the short-circuit point, the slope of the I-V curve at short circuit is given by
(10)(dIdV)Isc=-1Rshref.
A simultaneous resolution of (6)–(10) yields the five reference parameters, and then the I-V curve for STC is built using (5). As one can notice, in this strategy we made use of only three critical points of the I-V curve, namely, the short-circuit point, the maximum power point, and the open-circuit point.


## 3.2. Dependence of Operating Conditions

The determination of theI-V characteristic in all other operating conditions can be made possible by evaluating the various parameters of the I-V equation in the new weather environment using the semiempirical relationships that follow as shown in [21]. One may consider the following:
(11)aaref=TcTcref,I0I0ref=[TcTcref]3exp[EGNsaref(1-TcTcref)],Iph=GGref[Iph,ref+αcc(Tc-Tcref)],GGref=[-(I0/a)exp[(Vmp+ImpRs)/a]-1/Rsh1+(I0Rs/a)exp[(Vmp+ImpRs)/a]+Rs/Rsh]×[[1+(I0refRsrefaref)exp[(Vmpref+ImprefRsref)aref]+RsrefRshref]-1-I0refarefexp[(Vmpref+ImprefRsref)aref]-1Rshref/×(1+(I0refRsrefaref)exp[(Vmpref+ImprefRsref)aref]+RsrefRshref)-1]-1,Rsh=Rsh,ref.
Finally at the MPP, we write the following two equations:
(12)Imp=ImprefGGref,Vmp=Vmpref+βVoc(Tc-Tcref).


## 4. Five-Point Model

The five-point model is an improved method to extract the five unknown parametersa, Iph,  I0,  Rs,  and  Rsh of (3) using five key physical quantities of the I-V characteristic: the short-circuit current, the open-circuit voltage, current and voltage at the maximum power point, and the nominal series and shunt resistances Rs0 and Rsh0, respectively. This method was first introduced in [22] using experimental values of key physical quantities. Data issued by manufacturers of modules do not contain the values of the nominal series and shunt resistances. Indeed, the exact five-point model and the approximate five-point model that use approximate values of the nominal series and shunt resistances yield comparable results [22], provided that Rs0 and Rsh0 are well defined.


### 4.1. Calculation of Key Quantities

The nominal series and shunt resistances are determined using the slope method. The quantitiesRs0 and Rsh0 are defined at STC as the slope of I-V curves at I=0 and at V=0, respectively, as it can be seen in following equations [22]:
(13)Rs0=-(dVdI)V=Voc,(14)Rsh0=-(dVdI)I=Isc.
In this paper, the concept of nonlinearity is introduced for a better prediction of the behaviour of the tested photovoltaic modules. Thus, the dependency of the others key quantities on temperature and irradiance is given by the following equations:
(15)Isc=[Iscref+μIsc(Tc-Tcref)](GGref)α,Voc=Vocref1+βln(Gref/G)(TcrefTc)γ,Imp=[Impref+μImp(Tc-Tcref)](GGref)α,Vmp=Vmpref1+βln(Gref/G)(TcrefTc)γ.
Exponents α and γ are constants describing all effects related to the nonlinear dependence of current on the intensity of solar radiation and considering all the nonlinear effects of voltage variation with the cells temperature, respectively [23]. β is a dimensionless coefficient related to the specific technology of photovoltaic module [24, 25].


### 4.2. Calculation of Parameters

After expressing the above essential quantities as functions of operating conditions, the five parameters required are evaluated using the following equations [26, 27]:
(16)Iph=Isc(1+RsRsh)+I0[exp(IscRsa)-1],I0=(Isc+VocRsh)exp(-Voca),a=(Vmp+ImpRs0-Voc)×[ln(Isc-VmpRsh-Imp)-ln(Isc-VocRsh)+ImpIsc-(Voc/Rsh)(Isc-VmpRsh-Imp)]-1,Rs=Rs0-[aI0exp(-Voca)],Rsh=Rsh0(GrefG).


### 4.3. Calculation of the Constants

The constantsα, β, and γ are related to important effects in the photovoltaic module. They are empirical, and the precision of model depends on the acuity of their calculated results. Accurate determination of these constants requires significant experimentation to find the values of current and voltage at different points for different operating conditions in view to calculate them using (17)–(19):
(17)α=ln((Isc1-μIsc(T1-T0))/Isc0)ln(G1/G0),
where Isc, G, and T are values of short-circuit current, irradiance, and temperature, respectively; and the indices 0 and 1 are related to the fact that these quantities have to be measured for two different operating conditions.The variation of the open-circuit voltage is related to the variations of the solar radiation intensity and cells temperature. To calculateβ, temperature is kept constant while irradiance varies from G0 to G1. β is obtained as
(18)β=(Voc0/Voc1)-1ln(G0/G1).
For the same value of irradiance, the temperature of the cells varies from T0 to T1, and γ is obtained as follows:
(19)γ=ln(Voc0/Voc1)ln(T1/T0).


## 4.1. Calculation of Key Quantities

The nominal series and shunt resistances are determined using the slope method. The quantitiesRs0 and Rsh0 are defined at STC as the slope of I-V curves at I=0 and at V=0, respectively, as it can be seen in following equations [22]:
(13)Rs0=-(dVdI)V=Voc,(14)Rsh0=-(dVdI)I=Isc.
In this paper, the concept of nonlinearity is introduced for a better prediction of the behaviour of the tested photovoltaic modules. Thus, the dependency of the others key quantities on temperature and irradiance is given by the following equations:
(15)Isc=[Iscref+μIsc(Tc-Tcref)](GGref)α,Voc=Vocref1+βln(Gref/G)(TcrefTc)γ,Imp=[Impref+μImp(Tc-Tcref)](GGref)α,Vmp=Vmpref1+βln(Gref/G)(TcrefTc)γ.
Exponents α and γ are constants describing all effects related to the nonlinear dependence of current on the intensity of solar radiation and considering all the nonlinear effects of voltage variation with the cells temperature, respectively [23]. β is a dimensionless coefficient related to the specific technology of photovoltaic module [24, 25].


## 4.2. Calculation of Parameters

After expressing the above essential quantities as functions of operating conditions, the five parameters required are evaluated using the following equations [26, 27]:
(16)Iph=Isc(1+RsRsh)+I0[exp(IscRsa)-1],I0=(Isc+VocRsh)exp(-Voca),a=(Vmp+ImpRs0-Voc)×[ln(Isc-VmpRsh-Imp)-ln(Isc-VocRsh)+ImpIsc-(Voc/Rsh)(Isc-VmpRsh-Imp)]-1,Rs=Rs0-[aI0exp(-Voca)],Rsh=Rsh0(GrefG).


## 4.3. Calculation of the Constants

The constantsα, β, and γ are related to important effects in the photovoltaic module. They are empirical, and the precision of model depends on the acuity of their calculated results. Accurate determination of these constants requires significant experimentation to find the values of current and voltage at different points for different operating conditions in view to calculate them using (17)–(19):
(17)α=ln((Isc1-μIsc(T1-T0))/Isc0)ln(G1/G0),
where Isc, G, and T are values of short-circuit current, irradiance, and temperature, respectively; and the indices 0 and 1 are related to the fact that these quantities have to be measured for two different operating conditions.The variation of the open-circuit voltage is related to the variations of the solar radiation intensity and cells temperature. To calculateβ, temperature is kept constant while irradiance varies from G0 to G1. β is obtained as
(18)β=(Voc0/Voc1)-1ln(G0/G1).
For the same value of irradiance, the temperature of the cells varies from T0 to T1, and γ is obtained as follows:
(19)γ=ln(Voc0/Voc1)ln(T1/T0).


## 5. Results and Discussion

The calculations were made for the Shell SP75 module from Shell Solar and the GES-P70 module from GESOLAR. Shell SP75 is a commercial module with an output of 75Wp, consisting of 36 single crystal silicon cells connected in series. GES-P70 is also a commercial module of 36 cells connected in series with an output of 70 Wp but in polycrystalline silicon material. Both panels have similar high fill factor values (≈0.75) and at first glance, their curves seem to be quite similar, but if more accurately examined, they show some differences, probably due to their different processing technology. The data of the modules are presented in Table 1.Table 1
Data of the photovoltaic modules used.
Parameters of modules at STC
Shell SP75
GES-P70
Short-circuit current
4,8 A
4,28 A
Open-circuit voltage
21,7 V
21,8 V
Maximum power point current
4,4 A
3,98 A
Maximum power point voltage
17 V
17,6 V
Temperature coefficient ofIsc
0,002 A/°C
0,0017 A/°C
Temperature coefficient ofVoc
−0,076 V/°C
−0,069 V/°C
Number of cells
36
36
Maximum power
75 W
70 WManufacturers usually provide only limited operational data for PV modules. These data are available at STC and sometimes under nominal operating cell temperature (NOCT) conditions, that is, irradiance of 800 W/m2, temperature of 45 or 47°C (depending on the manufacturer), ambient temperature of 20°C, wind speed of 1 m/s, and AM 1.5 spectrum.The constantα, β, and γ can be calculated following the constant estimation procedures given above, and the results are shown in Table 2.Table 2
Constants estimation for PV modules.
Constant
α
β
γ
Shell SP75
1,022
0,058
1,116
GES-P70
0,974
0,060
1,076Utilization of available data aforementioned in the equations above presented for the five-parameter and the five-point models yields computed parameters of theI-V characteristic at given operating conditions for both modules used. Following Figures 2 and 3 shows that I-V curves fitted from computed parameters of five-parameter and five-point models and experimental data issue by the manufacturers of both the photovoltaic modules used.CalculatedI-V curves and experimental data of Shell SP75 (a) and GES-P70 (b) at different irradiance, 25°C.
(a)
(b)CalculatedI-V curves and experimental data of Shell SP75 (a) and GES-P70 (b) at different temperatures, 1000 W/m2.
(a)
(b)Figure2 shows the I-V curves for the modules Shell SP75 Figure 2(a) and GES-P70 Figure 2(b) for different levels of irradiance. It can be seen that the five-point model strongly agrees with experimental data than the five-parameter model for both types of modules. Another significant observation is the very accurate prediction in low irradiance of the five-point model.Besides verifying the accuracy of the models to fit with the solar irradiance, it is also very important to analyze their capabilities to reproduce the way photovoltaic panel performance is affected by the silicon temperature. TheI-V curves of both modules have been realized by both mathematical models used in this work when subjected to the variation of temperature as it is shown in Figure 3. The accuracy of both models was tested for different levels of temperature. All measurements were performed at 1000 W/m².As it has been noticed previously with the variation of irradiance, Figure3 shows that five-point model better fits I-V curves of both modules than five-parameter model when subjected to the temperature variation. It is well known that irradiance and temperature of the cells strongly affect the performance of operating photovoltaic devices.However, latest observations made on modelledI-V characteristics when subjected to irradiance and temperature of the cells variations are insufficient to make a true judgement about the capability of both the five-point and the five-parameter models to reproduce the behaviour of operating photovoltaic modules. To further investigate the performance of both the models in various operating conditions of both analyzed photovoltaic modules, inaccuracies on I-V curves prediction may be quantified.Root mean square error (RMSE) is used as the statistical tool for assessing the performance of the models to predict the current. The %RMSE compares the predicted (I-V) characteristic at different operating conditions to the measured one, thus providing a clear picture of the precision of the model. Figures 4 and 5 show %RMSE for five-point and five-parameter models for each operating condition analyzed.%RMS for Shell SP75 (a) and GES-P70 (b) at different irradiations, 25°C.
(a)
(b)%RMS for Shell SP75 (a) and GES-P70 (b) at different temperatures, 1000 W/m².
(a)
(b)Figures4 and 5 show the measuring deviation of computed I-V curves from experimental I-V data issue by manufacturers of modules for the whole of operating conditions analyzed. In general, five-point model exhibits lower %RMSEs than five-parameter model for both the studied photovoltaic modules for all environmental conditions. These observations are as expected because it has been above observed that computed I-V curves using five-parameter model are more unfasten from experimental I-V data than those computed using five-point model.In this work, it has been noticed that the five-parameter model for which nonlinearities are neglected either over-estimates or underestimates the values of current while the five-point model that takes into account the effects of nonlinear shows a remarkable fitting around the key points and a sufficiently accurate prediction at the other points of theI-V characteristic for the whole set of the weather conditions analyzed. This is to be foreseeable because the formulation of the equations of the five-parameter model presents an obsolete dependence on the environmental variations. The five-parameter model like many photovoltaic models used up to now, in addition to considering that the variations of the current and voltage with irradiance and temperature of cells, respectively, at the key operating points are made in linear ways or deviate slightly from their linearities, does not well define the dependence on irradiance and temperature of the cells of certain parameters. The shunt resistance and the ideality factor are considered independent from irradiance and temperature of the cells; the variations of the parameter a (modified ideality factor) with respect to temperature are not related to ideality factor but to its definition such as the above mentioned in (4). It is obvious that the use of constant parameters determined under STC conditions must bring about inaccuracies in other operating conditions [8, 28, 29].In conventional one-diode representation of photovoltaic module, physical phenomena are considered globally, so that parameters (which are related to these phenomena) lose their meaning [30]. In the five-parameter model, variations of parameters with environmental data (irradiance and temperature) are from their physical conception. Contrary to this technique, five-point model uses only experimental data and information of key operating points at STC and at any other operating conditions. Furthermore, the effects of nonlinearities considered in key operating points are transported in any other point of the characteristic through the calculation of parameters using key physical quantities. Methods that use a minimum simplification assumption from (3) such as five-point model are expected to give more realistic values for parameters under investigation [22, 30].


## 6. Conclusion

In this work five-parameter and five-point models were used to predict electricalI-V characteristics of illuminated commercial photovoltaic modules Shell SP75 and GES-P70 of monocrystalline and polycrystalline silicon material, respectively. Validation of these models was carried out with comparison to the data provided by module manufacturers. The following conclusions were drawn from our analysis.(i)
Most parameters depend on both the cell temperature and the solar irradiance. Thus, the mastery of their behavior is crucial to correctly predict the performance of the photovoltaic cells and arrays.(ii)
For both models used in this work, the precision varies on one hand from one module to another and on the other hand from one operating condition to another. Furthermore, the five-point model is more precise than the five-parameter model as well for the two photovoltaic module types studied for the whole set of weather conditions analyzed.(iii)
Taking into account the effects of nonlinearities has made it possible to achieve better fitting of theI-V characteristics. This study therefore shows that the five-point model is an excellent tool for the prediction of the electrical response of the illuminated crystalline photovoltaic modules.



---

*Source: 321694-2013-06-06.xml*