October 2011, Volume 2, No.

5 International Journal of Chemical and Environmental Engineering

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

S. Sheik Mohammed-a, a Faculty of Engineering, Dhofar University, PB. No. 2509, Salalah, Sultanate of Oman, PC-211. a Corresponding Author E-mail: [email protected]
Abstract This paper presents modeling of Photovoltaic (PV) module using MATLAB/Simulink. The model is developed based on the mathematical model of the PV module. Two particular PV modules are selected for the analysis of developed model. The essential parameters required for modeling the system are taken from datasheets. I-V and P-V characteristics curves are obtained for the selected modules with the output power of 60W and 64W from simulation and compared with the curves provided by the datasheet. The results obtained from the simulation model are well matched with the datasheet information.

Keywords: Modeling, PV module, Simulation, MATLAB

1. Introduction
The concentration on the use of fossil fuels for energy supply is the main threat for the stability of the global climate system and our natural living conditions. To conserve our globe, the scientific community gave evidence that mankind has to decrease the green house gases emissions, mainly CO2 and methane, by 60 - 70% as a minimum until the year 2050 [1]. In order not to harm our natural living spaces and threaten their resilience, a renewed compatibility would require a suitable form of energy alternatives sources that should be independent, easily accessible, and low in cost and should be environmentally clean. Renewable energy, and in particular power generation from solar energy using Photovoltaic (PV) has emerged in last decades since it has the aforesaid advantages and less maintenance, no wear and tear. The main applications of PV systems are in either stand-alone systems such as water pumping, domestic and street lighting, electric vehicles, military and space applications [2-3] or grid-connected configurations like hybrid systems and power plants [4]. The main aim of this paper is to provide a reader with the fundamental knowledge on design and building the blocks of PV module based on the mathematical equations using MATLAB/Simulink. The principle and operation of the PV cell and the fundamental characteristics of PV cell are discussed in chapter II. In chapter III the mathematical model of the ideal PV cell and also the practical PV cell are described. The

simulation model developed using MATLAB/Simulink and the results obtained are presented and discussed in chapter IV.

2. Operation and Characteristics of PV or Solar Cells

2.1Principle of Operation of Solar Cell

An array of solar cells converts solar energy into a usable amount of direct current (DC) electricity. The solar panels are mainly made out of semiconductor material and silicon being the most abundantly used semiconductor.Solar cells are connected in series to increase the output voltage. Similarly, the cells in parallel will yield a higher current. Series connected cells are called as PV modules and the interconnection series and parallel combination of solar cells is an array. The operation of solar cells may be described from a PN junction where there are diffusion currents and drift currents for the direct and reverse polarization, respectively. Usually, the cells operate in reverse direction so that the current drift is desirable. When the PN junction is exposed to light, photons with energy greater than the gap of energy are absorbed, causing the emergence of electron-hole pairs. These carriers are separated under the influence of electric fields within the junction, creating a current that is proportional to the incidence of solar irradiation [5].
2.2 Characteristics of Solar Cell:

Solar cells naturally exhibit a nonlinear I-V and P-V

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

characteristics which vary with the solar irradiation and cell temperature. The typical I-V and P-V characteristics of solar cell are shown in figure 1.

saturation or leakage current of the diode, q is the electron charge [1.60 x 10-19 oC], k is the Boltzmann constant [1.38x10-23J/K], Tc is the cell operating temperature (oK) and A is the ideality constant of diode. Both k and Tc should have the same temperature unit, either Kelvin or Celsius. The ideality constant varies depends on PV technology [6]. Ideality constant of different PV technology is presented in Table I.
Table 1: Ideality constant (A) for PV technology Technology Si-mono Ideality Factor 1.2 1.3 1.8 3.3 5 1.5 1.5 1.3

Figure 1. Characteristics of solar cell

Si-poly a-Si:H a-Si:H tandem a-Si:H triple CdTe

The fundamental parameters related to solar cell are short circuit current (Isc), open circuit voltage (Voc), maximum power point (MPP), efficiency of solar cell and fill factor. Short Circuit Current is the current corresponds to the short circuit condition when the impedance is low and it is calculated when the voltage equals to zero. I (at V=0) = Isc (1)

CIS AsGa

Maximum Power Pointis the operating point at which the power is maximum across the load. Pm = VmIm (4)

ISC occurs at the beginning of the forward-bias sweep and is the maximum current value in the power quadrant. For an ideal cell, this maximum current value is the total current produced in the solar cell by photon excitation. Open Circuit Voltage is the voltage when the open circuit occurs and there is no current passing through the cell. The Open circuit voltage is calculated when the voltage equals to zero. V (at I=0) = Voc (2)

where, Vm is the maximum voltage and Im is the maximum current. Efficiency of solar cell is the ratio between the maximum power and the incident light power.

(5) Pin is taken as the product of the solar irradiation of the incident light (G=/1000), measured in W/m2, with the surface area (Ac) of the solar cell in m2. Pin = G. Ac (6)

VOC is also the maximum voltage difference across the cell for a forward-bias sweep in the power quadrant. Voc= Vm for forward-bias power quadrant. The open circuit voltage can be expressed as:

Fill Factor (FF) is essentially a measure of quality of the solar cell. It is calculated by comparing the maximum power to the theoretical power (Pt) that would be output at both the open circuit voltage and short circuit current together (3) where,Ipv is the light generated current (it is directly proportional to the solar irradiation), Is is the (7)

351

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

Typical fill factors range from 0.5 to 0.82. The fill factor diminishes as the cell temperature is increased. (12)

3. Mathematical model of PV cell

The ideal photovoltaic cell is shown in figure 2. The mathematical equation for I-V characteristics of the ideal cell is given by equations (8), (9) [7]. Isc is the short-circuit current of cell at 25C and 1000W/m2, KI is the short-circuit current temperature co- efficient of cell, Tc and Tr are the working temperature of cell and reference temperature respectively in oK. The diode saturation current of the cell varies with the cell temperature, which is expressed in [11] as, where,

(13)
Figure 2. Equivalent circuit of Ideal PV Cell

where, Irs is the reverse saturation current of a cell at a reference temperature and a solar irradiation, Eg is the band gap energy of the semiconductor used in the cell (Eg 1.12 eV for the polycrystalline Si) at 25 oC (8) [11].

where, Id is the Shockley equation expressed as

and it can be

(9) A general mathematical description of I-V output characteristics for a PV cell has been studied for over the pass four decades [8]-[10]. The equivalent circuit of solar cell comprised of a current source connected in antiparallel with the diode, series resistance of the cell Rs and the shunt resistance of the cell Rsh is shown in figure 3. Based on figure 3, the output current of the solar cell can be calculated as
Figure.3 Equivalent circuit of solar cell

The reverse saturation current of a cell Irs is

(14) The work proposed in [12] claims that the only two unknowns in equation (11) are Rs and Rsh, and it may be found by making the maximum power calculated from I-V curve of the model (Pm,m) equal to the maximum power obtained from the experiment i.e., the power given in the datasheet (Pm,e) and solving the resulting equations for Rs, using iteration method by slowly increasing the value of Rs from zero. (i.e.,) the iteration starts with Rs=0. The equation which satisfies Pm,m = Pm,e is given as,

(10)

(11) where, Ns is number of cells in series for a PV module. The light generated current of the solar cell is mainly depends on the solar irradiation level and its working temperature, which is expressed as

352

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

Table 2. Electrical Characteristics of MSX 60/MSX 64 PV Module

Description Maximum Power (Pm)

MSX 64 64 W 17.5 V 3.66A 62W 4.0a 21.3A

MSX 64 60W 17.1A 3.5A 58W 3.8A 21.1A

o

(15)

Voltage @ Pmax (Vm) Current at @ Pmax (Im) Guaranteed Minimum Pm Short Circuit current (Isc)

(16)

Open Circuit voltage (Voc) Temperature co-eff of Voc Temperature co-eff of Isc Temperature co-eff of power NOCT2

-(8010) mV/ C (0.0650.15) %/oC -(0.50.05) %/oC 472 oC

4. Modeling and Simulation of PV Module

The blocks of the model are developed using MATLAB /Simulink based on Equation (11).The Solarex MSX60/MSX64 PV modules are chosen for modeling. The typical electrical characteristics of both MSX60 and MSX 64 are given in table 2. These modules consist of 36 polycrystalline silicon solar cells electrically configured as two series strings of 18 cells each [13]. The blocks developed using MATLAB/simulink for the PV module is shown in figures 4, 5, 6 and 7. The values of Rs and Rsh are calculated using equations (15), (16) and the optimum values which satisfies Pm,m=Pm,e are Rs = 0.1 and Rsh =161.344 The I-V characteristics and P-V characteristics curves obtained from the simulation for MSX 60 and MSX 64 PV modules at Tc=25C (298.15 K) and G=1 are presented in figures 8, 9, 10 and 11. The guaranteed minimum output power of MSX 60 module is 58W and MSX 64 module is 62W [13]. The output obtained from the model exactly matches with the data provided by the datasheet for both I-V characteristics and P-V characteristics of MSX 60 and MSX 64.

5. Conclusion
A PV module model based on the mathematical model of solar cell is developed using MATLAB/Simulink blocks. The essential input parameters such as Vm, Im, Voc, Isc, Ns, KI, Tc and G are taken from the manufacturers datasheet for the typical 60W and 64W modules selected for analysis. The I-V and P-V characteristics outputs are generated using the developed model for the selected modules and the obtained results are well matched with the datasheet information.

Figure 4. Simulation model for calculation of Ipv

Figure 5. Simulation model for calculation of Irs

353

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

Figure 6. Simulation model for calculation of Is

Figure 7. Simulation model for calculation of I

Figure 8. Simulation results of I-V characteristics of MSX60 at o Tc=25 C and G=1

Figure 10. Simulation results of I-V characteristics of MSX64 PV at o Tc=25 C and G=1

Figure 9. Simulation results of P-V Characteristics of MSX60 at o Tc=25 C and G=1

Figure 11. Simulation results of P-V characteristics of MSX64 at Tc=25 C and G=1

354

Modeling and Simulation of Photovoltaic module using MATLAB/Simulink

REFERENCES
[1] K. Tpfer Decarbonisation of energy supply central challenge for mankind Clean Power from Deserts, The DESERTEC Concept for Energy, Water and Climate Security, WhiteBook 4th Edition, Feb 2009. S. Leva, D. Zaninelli, Technical and Financial Analysis for Hybrid Photovoltaic Power Generation Systems, WSEAS Transactions on Power Systems, vol.5, no.1, May 2006, pp.831-838 S. Leva, D. Zaninelli, R. Contino, Integrated renewable sources for supplying remote power systems, WSEAS Transactions on Power Systems, vol.2, no.2, February 2007, pp.41-48 J.Schaefer, Review of Photovoltaic Power Plant Performance and Economics, IEEE Trans. Energy Convers., vol. EC-5,pp. 232-238, June, 1990. Markvart T.: Solar Electricity, John Wiley & Sons, 1994. C. C. Hua and C. M. Shen, Study of maximum power tracking techniques and control of dc-dc converters for photovoltaic power system, Proceedings of 29th annual IEEE Power Electronics Specialists Conference, vol. 1, 1998, pp. 86-93. [12] [7] [8] H. S. Rauschenbach. Solar cell array design handbook. Van Nostrand Reinhold, 1980. S. W. Angrist, , Direct Energy Conversion, Allyn and th Bacon, Inc., 4 edition, 1982, pp. 177-227. O. Wasynczuk, Dynamic behavior of a class of photovoltaic power systems, IEEE Transactions on Power Apparatus and Systems, vol. PAS-102, no. 9, 1983, pp. 3031-3037. J. C. H. Phang, D. S. H. Chan, and J. R. Philips, Accurate analytical method for the extraction of solar cell model parameters, Electronics Letters, vol. 20, no. 10, 1984, pp.406408. W. De Soto, S. A. Klein, and W. A. Beckman. Improvement and validation of a model for photovoltaic array performance. Solar Energy, 80(1):7888, January 2006. Marcelo Gradella Villalva, Jonas Rafael Gazoli, Ernesto Ruppert FilhoModeling and Circuit-based Simulation of photovoltaic arrays 10th Brazilian Power Electronics Conference (COBEP), 2009. Solarex MSX60/64 datasheet. Available online http://www.californiasolarcenter.org/newssh/pdfs/SolarexMSX64.pdf at:

[2]

[9]

[3]

[10]

[4]

[11]

[5] [6]

[13]

355