Abstract

This paper presents a new current-voltage (I-V) model for solar cells. It has been proved that series resistance of a solar cell is related to temperature. However, the existing five-parameter model ignores the temperature dependence of series resistance and then only accurately predicts the performance of monocrystalline silicon solar cells. Therefore, this paper uses Chebyshev polynomials to describe the relationship between series resistance and temperature. This makes a new parameter called temperature coefficient for series resistance introduced into the single-diode model. Then, a new six-parameter model for solar cells is established in this paper. This new model can improve the accuracy of the traditional single-diode model and reflect the temperature dependence of series resistance. To validate the accuracy of the six-parameter model in this paper, five kinds of silicon solar cells with different technology types, that is, monocrystalline silicon, polycrystalline silicon, thin film silicon, and tripe-junction amorphous silicon, are tested at different irradiance and temperature conditions. Experiment results show that the six-parameter model proposed in this paper is an I-V model with moderate computational complexity and high precision.

1. Introduction

Photovoltaic (PV) power generation system directly converts solar energy into electrical energy by using PV arrays. To obtain higher energy efficiency, PV power generation systems need to establish their simulation models to get optimized parameters. PV power generation system is mainly composed of PV arrays, controller [13], and inverter [46]. PV arrays, the core devices of PV power generation system, usually consist of solar cells in series and/or in parallel. There are two kinds of popular simulation models for solar cells, that is, the single-diode model and the double-diode model [710]. Because of having fewer parameters, the single-diode model is simpler than the double-diode model. This makes the single-diode model widely used. Some efforts have been made to improve the accuracy of the single-diode model. According to the number of parameters, there are mainly three kinds of models, that is, the four-parameter model, the five-parameter model, and the seven-parameter model.

The four-parameter model in [11] includes ideality factor, diode reverse saturation current, light-generated current, series resistance. The four-parameter model has the fewest cell parameters, and then its expression is the simplest. However, the predicted accuracy of the four-parameter model is very limited. And the four-parameter model is validated only for monocrystalline silicon modules. On the basis of the four-parameter model, the five-parameter models are obtained by introducing shunt resistance [1121]. The five-parameter models are more precise than the four-parameter model in [11]. Compromising the number of parameters and the approximate accuracy, the five-parameter models are the most commonly used. The five-parameter models proposed in [11, 1921] utilize the reciprocals of the slopes at the open-circuit point and short-circuit point to calculate the cell parameters. The two slope values are not usually provided by the manufacturers’ datasheet. And they can be obtained by using enough data pairs of experimental characteristic under the certain condition. This makes it very complicated to obtain the slope values. Therefore, these five-parameter models actually have seven parameters rather than five parameters. The five-parameter model in [12] does not utilize the above-mentioned two slope parameters and can accurately predict the performance of monocrystalline and polycrystalline silicon modules. However, the five-parameter model in [12] exists large modeling errors for amorphous silicon solar cells. This may be due to the facts that series resistance in [12] is assumed to be a constant and the temperature dependence of series resistance is ignored. Here, temperature refers to the temperature of solar cells. The seven-parameter models proposed in [13, 22] add two additional parameters on the basis of the five-parameter model. For example, in [13], the two additional parameters are temperature coefficient for series resistance and the diode reverse saturation current radiation dependence, respectively. Temperature coefficient for series resistance describes the temperature dependence of series resistance. The value of series resistance changes exponentially with the cell temperature. The diode reverse saturation current radiation dependence describes the influence of irradiance on diode reverse saturation current, which is obtained by using the current and the voltage at maximum power point under the irradiance of and the temperature of . The extraction process for this parameter is very complicated. Although the seven-parameter model in [13] provides a more accurate characteristic model, it is not widely used because of the complexity of parameter extraction. Therefore, we propose a new six-parameter model.

In this paper, Chebyshev polynomials are developed to describe the relationship between series resistance and temperature. This makes a new parameter called temperature coefficient for series resistance introduced into the single-diode model. Therefore, the single-diode model has six cell parameters in this paper, that is, ideality factor, light-generated current, diode reverse saturation current, series resistance, shunt resistance, and temperature coefficient for series resistance. The new six-parameter model provides a simpler and more reasonable expression for series resistance and accurately predicts characteristic for silicon solar cells. Five kinds of silicon solar cells with four different technology types, that is, monocrystalline silicon, polycrystalline silicon, thin film silicon and triple-junction amorphous silicon, are tested to validate the six-parameter model in this paper.

2. The Single-Diode Model of Solar Cells

A solar cell is essentially a very large area p-n junction diode. Under illumination, a single-diode model of a solar cell can be described by an equivalent circuit, shown in Figure 1. The equation of a single-diode model is shown as follows:where () is light-generated current, () is diode reverse saturation current, () is shunt resistance, and () is series resistance, is ideality factor which is defined as follows:where is the number of solar cells in series in a PV module, is diode ideality factor, and is thermal voltage which is defined as follows:where is Boltzmann’s constant ( J/K), () is cell temperature, and is electron charge ( C).

3. Chebyshev Polynomials Are Used to Model Series Resistance in the Six-Parameter Model

Temperature dependence of series resistance has been quantitatively analyzed and various expressions are obtained [18, 21, 2327]. Most of the expressions are nonlinear and implicit formulations about temperature. In [18, 21, 24], the relationship between series resistance and temperature is expressed as the exponential function of temperature. In [25], series resistance is literally expressed as the linear function of temperature. In fact, the expression also contains other electrical parameters at the same time, such as the short-circuit current and the voltage and the current at maximum power point, which are closely related to temperature. Therefore, the expression of series resistance in [25] is actually complex nonlinear. In addition, [13, 27] have given explicit representation of series resistance about temperature. Reference [27] has proved that series resistance of a solar cell belongs to positive temperature coefficient type. However, some of the series resistance of solar cells monotonically increases with the increase of temperature [23], but some of the series resistance of solar cells monotonically decreases with the increase of temperature [28, 29]. Therefore, for different kinds of solar cells, the expression in [27] for temperature dependence of series resistance is limited. As is known, the expression of characteristic is a transcendental equation and series resistance is an exponential function of temperature. When series resistance with the exponential function is substituted into characteristic model, the processes of characteristic modeling and parameter extraction will become complex and have computational burden. Therefore, in this paper, we follow the relationship between series resistance and temperature in [13] to develop a new expression of series resistance about temperature. Chebyshev polynomials are used to describe temperature dependence of series resistance. It is worth mentioning that the new expression of series resistance about temperature is an explicit form in this paper.

3.1. Chebyshev Polynomials Is Used to Model Series Resistance

In this paper, we assume the range of temperature is from 233 K to 353 K and the magnitude order of is . Here, denotes temperature coefficient for series resistance. To satisfy the need of accuracy and simplification, we select cubic Chebyshev polynomials [30, 31] to establish the expression for temperature dependence of series resistance. The remainder of this cubic Chevbyshev polynomials is shown as follows in this section.

Since the magnitude order of is , we have

Letwhere denotes the temperature of at Standard Reference Condition (the irradiance is 1000 W/m2, the temperature is 298 K, and air mass is 1.5, SRC). Then, we have . Since , the following formula holds:

According to (6), we have

Substituting (4) into (7), we can obtain .

Now, we use cubic Chebyshev polynomials [30, 31] to describe the expression for temperature dependence of series resistance as follows:where is the series resistance, is the series resistance at SRC, , , , , , and is defined as follows:where is defined as follows:

3.2. Six-Parameter Model

Manufacturers usually provide some data at SRC, such as short-circuit current , open-circuit voltage , current and voltage at maximum power point, , , and their temperature coefficient , , , and .

To obtain current corresponding to different voltage at various conditions, , , , , and must be known. Firstly, we calculate these parameters at SRC, that is, . Here, the parameters with subscript “ref” denote the values at SRC. According to (8), we know that series resistance is related to the temperature coefficient for series resistance . Therefore, also needs to be calculated. Thus, we need six pieces of independent information to calculate the six parameters. In this paper, we make full use of the definitions of , , , , , , , , and manufacturers’ data at SRC to obtain the above-mentioned six parameters.

According to the definition of short-circuit current at SRC, we have when . Equation (1) is written as follows:

According to the definition of open-circuit voltage at SRC, we have when . Equation (1) is written as follows:

According to the definition of maximum-power point at SRC, we have when . Equation (1) is written as follows:

At maximum-power point, the derivation of power is zero, we can obtain

For now, we have four pieces of independent information with six unknown parameters. We need to add two pieces of information. and (the temperature coefficient for maximum power) are used as the added information. It is worth mentioning that the units of and are V/K and 1/K, respectively.

The temperature coefficient for open-circuit voltage () can be defined as follows:

The temperature coefficient for maximum-power point () can be defined as follows:

In (15) and (16), temperature can change from 288 K to 308 K. can be given by manufacturers, is also given by manufacturers or calculated by using the values of and . But and are unknown in (15) and (16). To compute and , we have to obtain the values of , , , and at any temperature and irradiance conditions. Therefore, (17) show the expressions of , , , and at any temperature and irradiance conditions.

At any conditions, , , , and can be expressed as follows [7]:where is Boltzmann’s constant and is the material band gap. is the material band gap at SRC shown in Table 1, which depends on different technology types of solar cells. In addition, the value of at any conditions is shown as follows:

Form (5) and (8), the expression of series resistance can be rewritten as follows:

According to the data provided by manufacturers at SRC, we can obtain , , , , , and by using (11)–(16). Then, (17)–(19) are used to determine , , , , and of characteristics for different solar cells at any temperature and irradiance conditions. Once the five parameters are obtained at certain temperature and irradiance condition, we can easily obtain the current corresponding to different voltage according to (1). Therefore, (11)–(19) are called the six-parameter model. , , , and at any conditions are significant and have practical application. However, (11)–(14) are special representation of (1) at SRC for three important points: (), (), and (). That is to say, (11)–(14) can directly obtain , , , and at any conditions when subscript “ref” is deleted.

4. Model Validation and Discussion

In this section, we select five kinds of PV modules with four technology types to validate the six-parameter model proposed in this paper. They are Siemens SP75 and Siemens SM-55 (monocrystalline silicon), Solarex MSX-64 (polycrystalline silicon), Anstropower APX-90 (thin film silicon), and USSC US-21 (tripe-junction amorphous silicon), respectively. These five PV modules are composed of different number of solar cells in series, shown in Table 1. The reference values of SP75, MSX-64, APX-90 are provided by NIST (Nation Institute of Standards and Technology). The reference values of US21, SM-55 are provided by SNL (Sandia National Laboratory), shown in Table 1. To analyse the accuracy of the six-parameter model proposed in this paper, the absolute error (AE), the relative error (RE), and the root mean squared error (RMSE) are used which are defined as follows:where is the number of experimental values. is the current of tested models. is the current of experimental values.

To validate the effectiveness of the six-parameter model proposed in this paper, we compare the new six-parameter model with the five-parameter model in [14] and the seven-parameter model in [13]. We compute the parameter values by Engineering Equation Solver (EES) software. The parameter values of the six-parameter model in this paper, the five-parameter model in [14], and the seven-parameter model [13] at SRC are shown in Table 2.

We extract the experimental data in [12] by the software called Get Data. GetData software can make every data point amplified 15 times the original size. Thus, we can control that the abscissa of data point (voltage) has an absolute error of less than 0.04 V and the ordinate of data point (current) has an absolute error of less than 0.006 A. The experimental data obtained by GetData software has enough number of significant digit to achieve enough precision. Therefore, the experimental data obtained by GetData software can be taken as real experimental data provided in [14]. Then, we use the software Matlab to plot the figures of characteristic curves and the absolute errors between compared models and experimental data. The characteristic curves of SP-75 module, MSX-64 module, and APX-90 module at different temperature and irradiance conditions are shown in Figures 24, Figures 57, and Figures 810, respectively. The corresponding absolute errors of characteristic curves for these three modules are shown in Figures 1119. AEs of the six-parameter model in this paper is very close to the seven-parameter model in [13] at every point. For higher irradiance, the absolute errors of the six-parameter model in this paper are less than that of the five-parameter model in [14] within the range of maximum-power point to open-circuit voltage point, shown in Figures 11-12, 14-15, and 17-18. For lower irradiance, AEs of the six-parameter model in this paper and the five-parameter model in [14] are almost the same. In addition, We use the maximum absolute error and RMSE to further analyze the effectiveness of the six-parameter model in this paper. RMSE shows overall error. We select 97 data sets for SP-75 module, 101 data sets for MSX-64 module and 108 data sets for APX-90 module to compute RMSEs, see Table 3. From Table 3, we can find that that RMSEs of the six-parameter model in this paper are almost the same as that of the seven-parameter model in [13] and are far less than that of the five-parameter model in [14]. Therefore, in estimating of characteristic, the six-parameter model in this paper is obviously superior to the five-parameter model in [14], but slightly less than seven-parameter model in [13].

characteristic of PV modules are widely used to estimate the maximum power point. Next, we take the data of SM-55 module and US-21 module given in [32] as experimental data to validate the accuracy of the new six-parameter model in estimating the maximum power point. Tables 4 and 5 show and for the two PV modules at the temperature of and the irradiance of 1000  to 100 . For SM-55 module and US-21 module, the relative errors of the six-parameter model in this paper are almost the same as that of the seven-parameter model in [13]. For SM-55 module, the relative errors of the six-parameter model in this paper are far less than that of five-parameter model in [14] except for irradiance of 1000  and 100 . For US-21 module, the relative errors of the six-parameter model in this paper are far less than that of five-parameter model in [14] except for irradiance of 1000 . This shows that the six-parameter model in this paper can give a better estimation of maximum power point at lower irradiance for tripe-junction amorphous silicon module.

Therefore, the six-parameter model in this paper, used to estimate characteristic and the maximum power point, is obviously superior to the five-parameter model in [14] under most weather conditions, and is close to seven-parameter model in [13]. In fact, the six-parameter model in this paper makes a good compromise between the five-parameter model in [14] and the seven-parameter model in [13] at accuracy and computational complexity. Therefore, the six-parameter model is an model with moderate computational complexity and high precision.

5. Conclusions

We propose a new six-parameter model for solar cells in this paper. At calculation accuracy, the six-parameter model is better than the five-parameter model. The prediction of characteristic of six-parameter model is very close to the seven-parameter model, but less than seven-parameter model. However, at the computational complexity, the six-parameter model is superior to the seven-parameter model. Experiment results show that the new six-parameter model can provide a good prediction for silicon solar cells with different technology types, that is, monocrystalline silicon, polycrystalline silicon, thin film silicon, and tripe-junction amorphous silicon. That is to say, the model in this paper is more universal. It is worth mentioning that the relationship between the series resistance and temperature is described by using Chebyshev polynomials in this paper. Therefore, the six-parameter model proposed in this paper is an model with moderate computational complexity and high precision.

Nomenclature

:Ideality factor parameter which is defined as
:Ideality factor parameter at SRC
:Energy bandgap (eV)
:Energy bandgap at SRC (eV)
:Solar irradiance ()
:Solar irradiance at SRC ()
:Current (A)
:Light-generated current (A)
:Light-generated current at SRC (A)
:The current at maximum power point (A)
:The current at maximum power point at SRC (A)
:Diode reverse saturation current (A)
:Diode reverse saturation current at SRC (A)
:Short-circuit current (A)
:Short-circuit current at SRC (A)
:Boltzmann’s constant ( J/K)
:Diode ideality factor
:Number of solar cells in series in one module
:Power at maximum power point (W)
:Power at maximum power point at SRC (W)
:Electron charge (C)
:Temperature coefficient for series resistance
:Series resistance (Ω)
:Series resistance at SRC (Ω)
:Shunt resistance (Ω)
:Shunt resistance at SRC (Ω)
:Cell temperature (K)
:Cell temperature at SRC (K)
:Voltage (V)
:The voltage at maximum power point (V)
:The voltage at maximum power point at SRC (V)
:Open-circuit voltage (V)
:Open-circuit voltage at SRC (V)
:Thermal voltage (V)
:Temperature coefficient for maximum power current (A/K)
:Temperature coefficient for short-circuit current (A/K)
:Temperature coefficient for maximum power voltage (V/K)
:Temperature coefficient for open-circuit voltage (V/K)
:Temperature coefficient for maximum power (1/K)
:Temperature changing factor.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Program for New Century Excellent Talents in University under Grant NCET-11-1005, the Nature Science Foundation of Liaoning Province under Grants 201102005 and 201402014, the First Batch of Science and Technology Projects of Liaoning Province under Grant 2011402001, and Liaoning BaiQianWan Talents Program under Grant 2012921061.