# Modeling Of Photovoltaic Array Biology Essay

Published: Last Edited:

This essay has been submitted by a student. This is not an example of the work written by our professional essay writers.

According to the theory of semiconductors, an ideal PV cell can be considered as a current source which is in parallel with one or two diodes depending on the manufacturing technology and semi-conducting materials [11], [13]. However, in case of representing practical behavior of a PV cell, observations show that the model should include two additional elements: parallel resistance and series resistance [13]. Based on this model, the fundamental equation describing the current-voltage characteristics of a cell is given as:

where

_ Iph = Photo generated current

_ Is= Saturation current of diode

_ Rs= Cell series resistance

_ Rp= Cell shunt resistance

_ A= Diode quality factor

_ T=Cell operating temperature, Kelvin

_ q= Electronic charge, (1.6x10-19 C)

_ k= Boltzmann's constant, (1.38 x 10-23 J/K)

According to (1), the five parameters Iph, Is, Rs, Rp and A which are regarded hereafter as model parameters, should be known in order to construct the PV model. Model parameters are not usually provided by manufacturers. They can be either determined from the data provided in the datasheet or obtained through experiments. Nevertheless, model parameters do not have constant values as they all vary depending on temperature and solar irradiation levels. In other words, each set of temperature and solar irradiance corresponds to a unique set of model parameters. Therefore, in order to model PV characteristics under various environmental conditions, the effects of solar irradiance and temperature on PV parameters should be taken into account. One approach that has been widely proposed in the literature is to mathematically formulate the relationship of model parameters with temperature and solar irradiance [9], [12], [14]. Thus far, several equations have been developed. The most popular and accepted one is the equation derived for the photo -generated current, Iph, which relates it to both temperature and solar irradiance:

(2)

where, Iph(n) is the photo-generated current at the standard condition (25_C and 1000 W=m2). T and T0 are the cell actual and standard temperatures (K) respectively. Similarly, E and En are the actual and standard irradiance levels (W=m2) on the PV surface.

Unlike the photo-generated current, the diode saturation current, Is, has been reported to have a non-linear dependence on temperature. Despite its strong dependence on temperature, it is not considerably affected by irradiance changes:

(3)

where, Is(0) is the saturation current at the standard temperature and Eg is the band gap energy of semiconductor. Another electrical parameter that can be formulated is Voc:

Voc = Voc(n) + k(T - T0) (4)

In comparison to above equations, there has been less success in presenting consistent equations relating Rs and Rp to the environmental parameters. In [3] the following equations are presented:

(5)

(6)

where, k1 - k4 are constants and Rp(n) is the nominal parallel resistance at the standard condition. In many modeling techniques […], Rs and Rp are assumed to be constant under all environmental conditions, i.e. being equal to their nominal values. Such assumption simplifies the modeling technique at the cost of the accuracy. On the other hand, an accurate model is needed to analyze the PV characteristics under partially shaded condition in which PV array is modeled under both non-shaded and shaded conditions. Therefore, a comprehensive modeling of PV parameters is required to represent PV characteristics under various weather conditions.

In order to be worthwhile, the comprehensive modeling of PV should include another crucial stage which is experimental validation. Since the proposed technique in this work has been initially developed based on experimental study, it was able to generate satisfactory results in agreement with experimental results.

Modeling of PV cell based on experimental study

Formulation of model parameters based on environmental parameters

In order to model PV parameters based on environmental parameters, temperature and irradiance level, the above-mentioned equations need to be verified through experimental study. The experiments have been carried out on a single PV module subject to various temperature and uniform solar irradiation levels. By using an algorithm which is described in section III.B, the model parameters corresponding to experimental data have been determined. The model parameters together with temperature and irradiance values have been applied to the equations and ultimately the accuracy of each equation has been investigated. Table I provides eight sets of experimental data including the model parameters determined for each set. Based on the available experimental data, equations (2) and (4) have been validated and the parameters associated to each equation are determined. However, equations (3), (5) and (6) failed to satisfy the experimental conditions. Instead, the equations relating Voc, Impp and Vmpp to temperature and solar irradiance were proved to be effective. Impp and Vmpp are formulated in the same form as Isc and Voc, respectively. In other words, (2) can be applied to Iph, Isc and Impp while (4) can be applied to Voc and Vmpp.

## Table 1: Experimental data for various temperature and solar irradiance conditions

1

2

3

4

5

6

7

8

9

10

E(W/m2)

1073

1019

992

975

845

828

775

425

279

130

T (ËšC)

20

17

55

53

51.5

56

43.5

51

51

30

Isc (A)

4.3

4.01

3.98

3.781

3.44

3.52

3.34

1.87

1.35

0.69

Impp (A)

3.86

3.65

3.43

3.37

3.07

3.15

2.97

1.62

1.14

0.53

Vmpp(V )

54

54

50

50

50

48

50

52

52

54

Voc(V )

68.02

68.48

63.2

63.31

63.3

62.2

63.4

63.09

63.013

65.2

Iph (A)

4.335

4.0367

4.02

3.81

3.47

3.55

3.37

1.88

1.36

0.71

Is (μA)

268

157

283

229

188

401

96

102

733

10.3

Rs (mΩ)

17.7

20.4

15.7

16.9

18.5

21.5

19.9

18.8

25.2

71

Rp (Ω)

2.18

3.06

1.61

2.31

2.56

2.86

2.52

3.09

3.28

2.97

Equation (2) can be rewritten in the short form of:

(7)

where and

In order to find the coefficient values k0 and k1, two different sets of experimental data Iph, E and T are needed to be substituted in (7). From the available experimental data presented in Table I, all possible combinations of two sets have been attempted to ensure the effectiveness of equation (7) under various conditions. The resulted k1

has the mean value of -0.00247. Using this value for k1 in (7), k0 is determined by solving the equation for each set of data. Table II shows the results for k0. According to Table II, k0 varies in a small range and could hence be approximated by its mean value which is 21.648. The complete form of equation can be written as:

Iph = 21.648 x E(1 - 0.00247 x T) (8)

The similar procedure can be followed for Isc and Impp equations. The results for k0 corresponding to each equation are also shown in Table II. Incorporating the determined coefficients into (7) leads to:

Isc = 4.2348 x E(1 - 0.000015 x T) (9)

Impp = 7.104851 x E(1 - 0.0015 x T) (10)

The equation used for Voc and Vmpp is the simplified form of (4), expressed as:

Voc = Voc(n) + k(T - T0) = k0 + k1T (11)

Like aforementioned equations for currents, the coefficients

k0 and k1 can be determined through substitution of two sets of experimental data Voc and T into (11). This method has been attempted for all possible combinations of two out of eight sets. The resulted mean value of k1 equals -0.09019 which, upon substitution in (11), leads to a single value of k0 for each set of data. The results are shown in Table II. According to Table II, the mean value of k0 equals 93.1442 and hence (11) can be written as:

Voc = 93.1442 - 0.09019 x T (12)

The same procedure can be followed for Vmpp. Consequently, equation (13) can be written for Vmpp:

Vmpp = 118.9274 - 0.2153 x T (13)

Having developed mathematical formulas (9) - (13), the PV electrical parameters can now be calculated for any temperature and solar irradiance levels. The model parameters are now needed to be determined based on these electrical parameters. This can be easily fulfilled through developing an algorithm, capable of determining model parameters from electrical ones; voltages and currents. The algorithm has been used earlier to determine the model parameters associated with each set of experimental data.

## Table 2. Coefficient k0 based on experimental study

case

k0 - ISC

k0 - IMPP

k0 - VMPP

k0 - VOC

1

4.025

6.340344

116.437

94.6351

2

3.95281

6.41459

117.0829

94.4456

3

4.03598

6.80046

120.618

92.7823

4

3.89700

6.75799

120.187

92.7119

5

4.09448

7.07407

119.864

92.5666

6

4.27350

7.51582

118.833

91.8725

7

4.33282

7.30343

118.142

91.9451

8

4.42621

7.41130

121.757

92.3115

9

4.87674

7.92853

121.757

92.2345

10

5.33965

7.50193

119.235

92.5275

Determination of PV model parameters

In this section, the PV model parameters; Iph, Is, Rs, Rp and A are determined based on electrical parameters; Isc, Impp, Vmpp and Voc. The objective is to find the model parameters such that the resulted I-V and P-V curves accurately match the experimental curves, specifically at three key points: short-circuit, maximum power and open circuit. In order to satisfy this criterion, the model parameters are obtained by solving the fundamental equation (1), for the key points. The electrical parameters for a particular weather condition are substituted in these equations. The values for Iph, Is, Rs and Rp are then determined through an iterative procedure. As A is an empirical value which expresses the quality of diode used in PV, an initial value can be chosen in order to obtain other model parameters. Depending on PV device structure and material, A usually varies in the range of 1 to 1.5. In this paper, A is assumed to be equal 1. This value can be modified later if fine tune of the model to the experimental study is needed.

In order to determine the other four parameters, four equations are required. In addition to the equations derived for short circuit, maximum power and open circuit conditions, the forth equation governs that the models maximum power occurs at the same voltage as the measured Vmpp.

The equations for the three key points can be written as:

(14)

(15)

(16)

The forth equation can be derived using the fact that the derivative of power with respect to voltage should be zero at the maximum power point, Vmpp and Impp:

(17)

The derivative of power can be written as:

(18)

According to (18), the derivative of current with respect to the voltage should be found first:

(19)

Solving the above equation for dI/dV and substituting the value in (18) results in:

(20)

Substituting Vmpp and Impp in (20) leads to the forth equation:

(21)

The iterative procedure is illustrated in Fig.1. As shown, by using (15), Impp(cal) is calculated and compared to its experimental value. Depending on the difference, Rs is slowly incremented or decremented. Based on the new value of Rs and above equations, Rp, Iph and Is can be calculated. These values are substituted in (15) in the next iteration and the same procedure repeats. The iteration continues until the difference between Impp(cal) and Impp reaches zero and all the values become stable. Unlike the algorithm proposed in [11] in which Rs is incremented starting from zero, in this paper an initial value for Rs has been assumed based on the equation proposed in [2]:

(22)

Assuming an initial value for Rs reduces the number of iterations needed to reach the final answer. It also has a great impact on the convergence of the algorithm. The other three parameters; Rp, Iph and Is are set to initial values [2]:

(23)

(24)

(25)

Depending on the software in which the algorithm is implemented, the number of iterations may be required as well. In Matlab, the algorithm reaches the stable final answer with less than 100 iterations, whereas in PSCAD, there is no need to define a number and it automatically reaches the stable answers within less than 0.1 second.