Effect of various model parameters on solar photovoltaic cell simulation: a SPICE analysis

In this paper, all the models of PV cell, namely ideal single-diode model, single-diode Rs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{\rm s}$$\end{document} model, single-diode Rp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{\rm p}$$\end{document} model, the two-diode model, and the three-diode model, have been discussed. SPICE simulation is done to evaluate the impact of model parameters on the operation of PV cell. The effects of the parameters are discussed. The photocurrent, IL,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_{\rm L},$$\end{document} is proportional to irradiance, and the series resistance, Rs,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{\rm s},$$\end{document} reduces the short-circuit current and fill factor. The parallel resistance, Rp,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{\rm p},$$\end{document} reduces the open-circuit voltage, and both the diffusion diode and recombination diode reduce the open-circuit voltage value and fill factor. Finally, it is shown that an increase in cell operating temperature reduces the open-circuit voltage and fill factor and thus degrades the performance significantly.


Background
A photovoltaic (PV) cell generates electricity when it is illuminated by the sun or some other light sources. Small PV cells can be grouped to form panels, and panels can be grouped to form arrays. The tremendous growth of PV industry and the increased number of installed PV systems all over the world raised the need for supervision and simulation tools for PV systems. To understand the PV system in a better way, modeling the PV system in various operating and weather conditions is necessary (Soto et al. 2006;Carrero et al. 2007). The simulation is also useful for various purposes like: (1) to analyze and study the behavior of power converters when they are attached with the PV system (Camps et al. 2015;Eccher et al. 2015), (2) to simulate the behavior of maximum power point tracker (MPPT) (Bendib et al. 2015;Ishaque et al. 2014), and (3) to estimate the efficiency of the PV system (Khazaei et al. 2015;Sivakumar et al. 2015). In addition to these, the simulation can also be used to analyze various operational conditions like: partial shading, change in irradiance, and fault conditions (Brano et al. 2010(Brano et al. , 2012. On the other hand, manufacturers of the PV modules provide electrical parameters only at standard test conditions (STC) which are irradiance = 1000 W/m 2 , 1.5 air mass (AM), and cell temperature of T cell = 25 °C (Soto et al. 2006). Manufacturers datasheet only provides the short-circuit current I SC , open-circuit voltage V OC , the voltage at maximum power point V MP , the current at maximum power point I MP , and the temperature coefficients at open-circuit voltage and short-circuit current (Sera et al. 2007). However, the operating range of the PV panel is large and operating and weather conditions also vary. Therefore, performance models are built to predict the performance of the PV system at any operating and weather condition. The model predicts the I-V characteristic of a PV system as a function of irradiance, angle of incidence of solar radiation, the spectrum of sunlight, and temperature (Chatterjee et al. 2011). For the last few years, substantial amount of work has been done to develop simulation models and extraction of model parameters of photovoltaic systems (Chouder et al. 2012;Ma et al. 2014;Villalva et al. 2009;Castaner and Silvestre 2002;Bal et al. 2012;Xie et al. 2014). Specially, Ma et al. (2014b) have developed a model to simulate the performance characteristics of crystalline silicon photovoltaic modules/strings/arrays. Zeroual et al. (1998) designed and constructed a closed loop sun-tracker with microprocessor management. Chin et al. (2011) have done a design, modeling, and testing of standalone single axis active solar tracker using MATLAB/Simulink. Chin (2012) has also done model-based simulation of an intelligent microprocessor-based standalone solar tracking system. Models for numerical device simulations of crystalline silicon solar cells were also reviewed by Altermatt (2011). De Blas et al. (2002 found suitable models for characterizing photovoltaic devices. Celik and Acikgoz (2007) have done the modeling and experimental verification of the operating current of monocrystalline photovoltaic modules using four-and five-parameter models. Chatterjee et al. (2011) have identified suitable photovoltaic source models. However, the effects of individual model parameters were not clearly reviewed in the present literature.
The objective of this work is to analyze the effects of model parameters on the simulation of PV cell. PSPICE is used to analyze and simulate the effects of parameters on photovoltaic cell performance.

How a photovoltaic cell works?
A photovoltaic cell is usually a semiconductor device that converts sunlight into electricity by the means of photovoltaic effect (Archer and Hill 2001). When light falls on a solar cell, the incoming photons can be absorbed, reflected, or passed through it. For a photon to be absorbed by the solar cell, the energy of the photon must be greater than the band gap energy of the cell (Luque and Hegedus 2011). The photon is then absorbed to generate pairs of mobile charge carriers (for example, electron and hole) which are then separated by the structure of the device (such as a p-n junction). This produces a potential difference and thus produces electrical current. The photovoltaic effect is shown by various materials. In most of the cases, semiconductor materials (like silicon) in the form of p-n junction are commercially used to produce solar cells (Green et al. 2015). The I-V characteristic of a p-n junction diode is given by the Shockley equation as (Luque and Hegedus 2011): where I D is the current generated by the diode, V D is the voltage across the diode, I 0 is the reverse saturation current of the diode, q = 1.602 × 10 19 Coulombs is the electron charge, T cell is the cell temperature in Kelvin, n is the diode ideality factor (dimensionless), and (1) I D = I 0 exp qV D nkT cell − 1 k = 1.38065 × 10 23 J/K is the Boltzmann constant. The reciprocal term of (q/nkT cell ) is called the thermal voltage of the diode. Therefore, thermal voltage of the diode is: Now Eq.
(1) becomes: The fill factor of a solar cell is given as: A semiconductor p-n junction can be made to operate as a solar cell. Figure 1 shows the basic structure of a PV cell. When light is incident on the cell, the photons of light generate free electron-hole pairs which are then attracted toward the junction. Figure 2 illustrates the operation of the solar cell. In the cell, electron flows from n-type silicon to p-type silicon. The generated current is proportional to the light intensity Archer and Hill (2001).

Different parameter models of solar cell
There are various electrical circuit models of photovoltaic cells which have been widely described in the literature (Saloux et al. 2011;Ishaque et al. 2011a, b;Ishaque and Salam 2011;Ishaque et al. 2011). Usually, a PV cell is represented by a current source (I ph ), one or two diodes (D) in parallel with the current source, a series resistance (R S ) or a parallel resistance (R p ) or both connected. The datasheet of a PV panel provided by the manufacturer usually includes few electrical and thermal parameters as given in Table 1. All these values are measured at standard test condition which is at 1.5 air mass (AM), 1000 W/m 2 , and 25 °C temperature. Usually, the models are known by the unknown number of values of the parameters which cannot be deduced directly from datasheet values. In the preceding section, the models and their associated parameters are described. (2)

Ideal cell model
In ideal condition, the solar cell is electrically equivalent to a current source in parallel with a diode as shown in Fig. 3a. The light-generated current, also known as photocurrent, is represented as I L , the diode current as I D , and the net current and terminal voltage of solar cell as I cell and V cell , respectively. The net current I cell available from the PV cell is given as: Substituting the value of I D from Eq. (3) into Eq. (4), we get: The photocurrent I L increases linearly with solar insolation. If diode current I D is small, more current is delivered by the PV cell. The diode ideality factor n indicates how closely the diode follows the ideal diode equation (Archer and Hill 2001). The value of n greater than 1 represents non-ideal condition, whereas n = 1 represents ideal behavior of the diode. As we can see from Eq. (3) that the ideal cell model has three parameters to find which are photocurrent I L , dark current I 0 , and diode ideality factor A. Therefore, this ideal model is also called the 3-p (three-parameter) model as shown in Table 2. This ideal cell model can be used to demonstrate the basic concept of PV cell, but is never used to simulate the actual operating conditions. This is because there are several limitations present in a practical PV cell (Castaner and Silvestre 2002). The series resistive losses, the shunt resistive losses, the recombination at the space-charge region of solar cells which is non-ohmic, and the non-ideality of the diffusion diode all indicate that ideal single-diode model is insufficient to accurately represent the maximum power delivered by the PV cell (Luque and Hegedus 2011).

Single-diode R s model
In fact, the current generated in the photovoltaic cell travels through semiconductor material which are not heavily doped and thus show resistivity (Duffie and Beckman 2006;Ma et al. 2014a). Besides this, the resistance of the metal grid, contacts, and current-collecting wires also contribute to the total series resistive losses. Usually, a lumped resistor, R s , is added in series with the ideal circuit model to represent these series losses as shown in Fig. 3b (Ma et al. 2014a). This model is called singlediode R s model.
Because of the addition of the lumped series resistor, R s , the net current flow through the cell, I cell , now becomes: From Eq. (4), it can be seen that there are now four unknown parameters: photocurrent I L , reverse saturation current I 0 , diode ideality factor A, and the newly added lumped series resistor R s . This model is therefore named as 4-p model. However, in recent research it is shown that the 4-p model which ignores the shunt resistance effect does not perfectly fit the experimental I-V and P-V data (Dongue et al. 2012).

Single-diode R p model
As the solar cells are made out of large area wafers, or from large area thin film material, a number of shunt resistive losses occur at n layer of the p-n junction PV cell. Localized short circuit or short circuiting of the cell

Electrical parameters
Maximum power rating

Thermal parameters
Normal operating cell temperature NOCT (°C) Temperature coefficient: short-circuit current (A/°C) Temperature coefficient: open-circuit voltage V (°C)

Standard test conditions (STC)
Air mass AM = 1.5 border are the most common form of shunt losses (Castaner and Silvestre 2002). This is represented generally by a lumped resistor, R sh , in parallel with the ideal device. The single-diode R p model which has both the series resistance R s and shunt resistance R p is shown in Fig. 3c.
As there are five unknown parameters in this model, it is called the 5-p model. Table 2 illustrates the models along with their parameters. The current through the parallel resistor, I P , can be given as:  Single-diode Rp model 5 I L , I 01 , n 1 , R s , R p I 02 , the reverse saturation current of the second diode Two-diode model 7 I L , I 01 , I 02 , n 1 , n 2 , R s , R p I 03 , the reverse saturation current of the third diode Three-diode model 9 I L , I 01 , I 02 , I 03 , n 1 , n 2 , n 3 , R s , R p n 1 , the first diode ideality factor n 2 , the second diode ideality factor n 3 , the third diode ideality factor R s , lumped series resistance R p , shunt resistance Sarkar Renewables (2016) 3:13 The governing equation of this 5-p model now becomes: This single-diode R p model or 5-p model is the most widely used and accepted model found in the literature (Chatterjee et al. 2011;Chenni et al. 2007;Chin et al. 2015;Ciulla et al. 2014;Mahmoud et al. 2012;Patel and Agarwal 2008;Brano et al. 2010Brano et al. , 2012. Most of the research has been done on the improved 5-p model and its parameter (I L , I 0 , V T , R s , R p ) extraction procedure (Soto et al. 2006;Brano et al. 2010Brano et al. , 2012Sera et al. 2007;Chatterjee et al. 2011;Villalva et al. 2009;Ishaque et al. 2011).

Two-diode model
Even though the commonly used 5-p model can achieve acceptable level of accuracy, the saturation current of the photovoltaic cell is the linear superposition of charge diffusion and recombination in the space-charge layer (Luque and Hegedus 2011) and therefore the more accurate electrical representation of the PV cell can be done by two Shockley diodes in parallel with a current source and associated series and shunt resistances (Castaner and Silvestre 2002). This model is more relevant at low-voltage bias, i.e., at low irradiance level operation of the PV cell (Luque and Hegedus 2011). Figure 3d illustrates the circuit of this model.
The equivalent circuit equation of the two-diode model is given as: We can see form Table 2 that the number of parameters is increased to seven in two-diode model. To reduce the number of unknowns, the ideality factor of the first diode is usually taken as 1 and the second diode ideality factor is taken as 2. This assumption is done based on the approximation of Shockley-Read-Hall recombination in the space-charge layer of the photodiode (Chih-Tang et al. 1957

Three-diode model
To take the influence of grain boundaries and leakage current through the peripheries into consideration, another diode can be added in parallel with the two other diodes of the two-diode model (Nishioka et al. 2007). This model has nine parameters (see Table 2) and a really cumbersome solution process. Most of the time, ignoring few parameters is done to reduce the number of equations of this model. Figure 3e shows the circuit diagram of this model. The equation of the current of this model is: This model is usually used to extract the electrical properties of multicrystalline Si solar cells. Because of the complexity of the parameter extraction process of threediode model, the use of this model is quite limited (Ma et al. 2014a).

Experimental setup
Siemens SM50 module is connected in series with variable resistors, and current and voltages are measured for varying parameters. An Eppley PSP pyranometer is used to take the radiation measurement. The pyranometer gives output in a range of 0-10 mV having sensitivity of approximately 8 µV/Wm −2 . The SM50 module has 36 monocrystalline silicon solar cells. The electrical and thermal parameters of this module are enlisted in Table 3. Two digital multimeters are used to take current and voltage measurements. Figure 4 shows the experimental setup.

Results and discussion
Simulation Program with Integrated Circuit Emphasis or SPICE is an open-source circuit simulator developed by Nagel andPederson of University of California, Berkeley, in 1973 (Nagel andPederson 1973). With time, SPICE became popular and adopted by various institutions. Several proprietary software tool gradually evolved from the basic SPICE. Cadence OrCAD PSPICE is a similar proprietary software tool used to simulate circuit using SPICE. In this study, we have simulated the parameters shown in Table 2 to check their effects on the PV system (14) performance. For this purpose, we have taken a hypothetical small PV cell with an area of 5 sq. inch, short-circuit current density 34.3mA/cm 2 , and dark current I 0 as 10 −11 A/cm 2 . Using a PSPICE function named PROBE, we got the simulation curves. The schematic diagram of the simulation is given in Fig. 5.

Effect of photocurrent, I L
The photocurrent I L is directly proportional to the irradiance received by the PV cell. The value of I L can be calculated as: where J sc is the short-circuit current density (mA/cm 2 ) , A is the area of the PV cell (cm 2 ), and G is the irradiance (W/m 2 ). By changing the value of G, we got several curves as shown in Fig. 6a. It is evident from the curves that the short-circuit current reduces significantly with reducing G. The fill factor (Eq. 4) also reduces to a great extent. Figure 6b shows the simulated power curves of the PV cell with varying irradiance. Power also reduces proportionally with irradiance.

Effects of the series resistance, R s
From Fig. 7, it can be seen that with the increase in series resistance R s the open-circuit voltage remains the same, but both the short-circuit current and fill factor reduce significantly. In practice, the value of the lumped series resistance R s is kept quite low to prevent this drastic degradation in performance.

Effects of the shunt resistance, R p
By doing a dc sweep using PSPICE, the effect of the shunt resistance can be found. Figure 8 shows the I-V curves with varying shunt resistance R s . From the figure, it can be seen that the value of open-circuit current remains same, but short-circuit current and fill factor reduce with the increase in shunt resistance.

Effects of the diffusion and recombination diode, D 1 and D 2
Usually the diffusion diode, D 1 , has an ideality factor close to unity. But in practice, because of the fabrication process, some non-ideality may remain. Similarly, the recombination diode should have the ideality factor as 2.
But it also exhibits some extents of non-ideality because  Fig. 5 Schematic diagram of PV cell simulation in PSPICE of the manufacturing process. The impact of both D 1 and D 2 is shown in Fig. 9a, b. It is evident that in both cases the short-circuit current remains same, but the open-circuit voltage and fill factor reduce to great extent.

Temperature effects, T cell
Because of the semiconductive property of the solar cell, temperature has a strong effect on its performance. In practice, solar cell can easily be warmed up to 60-65 °C. Using the built-in feature temperature sweep of PSPICE, we can check the effect of temperature on solar cell operation. Figure 10 shows the effect of cell temperature on its performance. It is seen that increased cell temperature reduces the open-circuit voltage and fill factor and thus reduces its performance significantly.

Experimental results
The experimental I-V curves provided by the manufacturer at different illumination and ambient temperature are shown in Fig. 11a. The measured I-V curve at 45 °C and approximately 800 W/m 2 is shown at Fig. 11b. From the figures, it is clear that with the increase in insolation short-circuit current increases drastically, but opencircuit voltage does not increase that much. With the increase in temperature, short-circuit current does not decrease, but open-circuit voltage decreases a lot. The experimental results match with the simulation perfectly.

Conclusion
The development of PV system raised the need of simulation of PV system. In this paper, all the modeling methods have been discussed and SPICE simulation is done to evaluate the impact of model parameters on the operation of PV cell. The models discussed here are ideal single-diode model, single-diode R s model, single-diode R p model, the two-diode model, and the three-diode model. All the model parameters are enlisted in Table 2. The effect of the parameters are discussed with PSPICE simulation. The photocurrent, I L , is proportional to irradiance, and the series resistance, R s , reduces the short-circuit

Competing interests
The author declare that he has no competing interests.