Next Article in Journal
The Communication System and its Impacts on Line Current Differential Protection in Distributed Feeder Automation
Next Article in Special Issue
On-Site Identification of the Material Composition of PV Modules with Mobile Spectroscopic Devices
Previous Article in Journal
Cost Optimization of a Stand-Alone Hybrid Energy System with Fuel Cell and PV
Previous Article in Special Issue
Theoretical and Numerical Study of a Photovoltaic System with Active Fluid Cooling by a Fully-Coupled 3D Thermal and Electric Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Diode Models of PV Modules: A Comparison of Conventional Approaches and Proposal of a Novel Model

1
Power System Department, Hanoi University of Science and Technology, Hanoi 11615, Vietnam
2
Department of Electrical, Electronic and Computer Engineering, Gifu University, Gifu 501-1193, Japan
*
Author to whom correspondence should be addressed.
Energies 2020, 13(6), 1296; https://doi.org/10.3390/en13061296
Submission received: 22 December 2019 / Revised: 4 March 2020 / Accepted: 5 March 2020 / Published: 11 March 2020
(This article belongs to the Special Issue Photovoltaic Modules)

Abstract

:
In this paper, the seven traditional models of photovoltaic (PV) modules are reviewed comprehensively to find out the appropriate model for reliability. All the models are validated using the Matlab code and graphical comparisons between models are made. The accuracy and convergence of each model is evaluated using the data of manufactured PV panels. Then, a novel model is proposed showing its consistent performance. The three most key parameters of the single-diode model are self-revised to adapt to various types of PV modules. This new method is verified in three types of PV panels’ data measured by the National Renewable Energy Laboratory (NREL), USA. The validated data show promising results when the error RMSEs’ range of the proposed model is under 0.36.

1. Introduction

The rapid exhaustion of conventional energy resources, such as coal, crude oil, and natural gas, has been threatening to energy security in the world. Since renewable resources are clean and inexhaustible, the penetration of renewable energy has been arising from time to time in academia, industry, business, and government. Concerning a 100% renewable goal in the next few decades worldwide, renewable energy has been becoming a hot topic in the research community. In this context, Photovoltaic (PV) energy becomes one of the prominent renewable energy resources. With an increase of more than 30% in 2018, the power generation from PV is estimated to be over 580 TWh [1]. Because of the technical improvements in PV and the dramatic drop in the price of solar PV panels, solar PV systems continue developing. In terms of generation, by 2050, solar PV could become the second-largest power generation source, behind the wind power. This growth would meet 25% of the total energy demand globally [2]. To estimate the economic feasibility of a solar PV system, the evaluation of the output power of a PV plant needs to be addressed. Identifying the standard I-V characteristics of specific PV panels is a key to estimating accurately the output power of PV panels. As a result, researchers have proposed various methods to predict the performance of a PV module based on modeling it, all of which can be divided into analytical methods and metaheuristic algorithms [3]. Some optimization techniques are employed in metaheuristic algorithms, and have been applied to obtain the model parameters as follows: Simulated annealing [4], bacterial foraging algorithm [5], genetic algorithm [6], differential evolution [7], partial algorithm [8], artificial bee colony [9], simplified swarm optimization [10], etc. While these models impose no restrictions on the problem formulation [11], metaheuristic algorithms might take long computational time and have convergence difficulty. Meanwhile, in [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33], analytics methods utilize a set of mathematical equations and assumptions from the datasheet provided by manufacturers. The trade-off between model accuracy and computational time makes the analytics methods much more popular in engineering applications than metaheuristic methods.
The best way to analyze the behavior of the PV generator is to adopt an equivalent circuit and to analyze relevant equations describing it [14]. A PV cell is described by an equivalent circuit consisting of a current source, at least one diode, and one resistor [15] based on the intrinsic characteristics. The single-diode model (SDM) [13,16,17] and the double-diode model [11,18,19] have gained the most researchers’ attentions, specifically the SDM. While the double-diode model has a certain advantage, that is representing the recombination loss in the depletion region [20], the double-diode model requires solving implicit nonlinear seven-parameters [11]. This has led to more difficulties in the calculation of the initial values, long computational time, and algorithm complexity. Meanwhile, the SDM gets the desired compromise between accuracy and simplicity [20].
Different methods have been proposed to extract the five parameters of the SDM. In [12,13,14,15,16,17,18,19,20,21], a set of implicit equations accompanied by assumptions and simplifications were used to formulate the equations obtaining the model parameters.
An iterative and an analytical method are conducted by [13,17,22], who proposed a method to estimate parameters through a trial-and-error approach.
Meanwhile, model parameters were identified by using mathematical manipulations in [14,15,23,24,25,26]. Specifically, the Lambert-W function at five experimental points was applied in [14,15,24,25,26].
References [3,34,35] reviewed the remarkable models from 2002 to 2017. The comparisons among models were discussed to choose the model that fits the applications in reality. However, all models are only investigated in their performances on one type of PV module. Consequently, the selection of a suitable model for an application is decided only based on computational time and model accuracy, which is insufficient to assess properties of various types of PV panels. In this article, we conduct experiments on three main types of PV modules under varied solar irradiance and cell temperature, such a performance of a PV system largely depends on cell temperature (T)–solar irradiance (G) conditions (T-G) [17]. Since the new method has proved its advantages, it is verified in three types of PV panels, the data of which are tested by the National Renewable Energy Laboratory (NREL), USA. Since the results have shown the promise of this method, it could be used in predicting the performance of a PV panel.
The main contributions of this paper are four folded:
  • Review previous models to build I-V curves of PV modules
  • Compare the accuracy of these models
  • Propose a higher performance model
  • Validate the proposed model by the real PV module’s data
The rest of this paper is organized as follows. Section 2 briefly introduces the fundamental model of the SDM. Section 3 describes reviewed models. Section 4 discusses the reviewed model and then Section 5 describes the proposed model. Section 6 shows the numerical results validating the reviewed and proposed models. Finally, the conclusion is drawn in Section 7.

2. Equivalent Circuit of the Single-Diode Model

Figure 1 illustrates the well-known equivalent circuit of the SDM, which consists of one current source, one diode, one series resistor, and one shunt resistor. The series resistor is added to take into account the voltage which drops across the transport resistances of the solar cell. The shunt resistor represents the effect of the leakage current in the p-n interface of the diode and the edges, as well as the shunt resistance occurring across the solar cell surface [27]. Normally, the estimated value of shunt resistance is very high compared to the series resistance.
The relationship between the output current-voltage is expressed by the following equation:
I = I p v I d V + I R s R s h
I p v is the PV current, generated by electron-hole pairs within a solar cell. I d is the Shockley diode current, which is described by Equation (2), where I s a t is the reverse saturation current measuring the leakage or recombination of minority carriers across the p n junction in reverse bias [25], which is calculated by means of Equation (3). The SDM assumes that the Shockley current can be described by a single exponential dependence modified by the diode ideality factor a [25]. E g a p is the bandgap energy of the semiconductor material ( E g a p = 1.12 eV for the crystalline Si at 25 ° C [30,31]).
I d = I s a t e q V a k T 1
I s a t = C T 3 e E g a p k T
In fact, PV generators are made of a number of PV cells connected in series and parallel, not referred to as a single PV cell. The operating voltage of the PV cell is a few hundreds of millivolts, while the current generating at high irradiance levels is of some amperes. As a result, to reach the desired voltage, a number of connected-series cells are arranged into PV modules. So some researchers have used a modified expression of Equation (1) as follows:
I = I p v I s a t e V a N s V t 1 V + I R s R s h
V t is defined as the thermal voltage of the cell:
V t = k T q

3. Methodology Extracting Model Parameters from Datasheet Values

All reviewed models use information from the datasheets provided by the manufacturers. Usually, manufacturers provide the values of significant points, i.e., its short-circuit current ( I s c ), its open-circuit voltage ( V o c ), its current and voltage at the MPP ( I m p p , V m p p ), the temperature coefficient of the short-circuit current ( K i ), and the temperature coefficient of the open-circuit voltage ( K v ). Unfortunately, those values are provided only at the standard test condition (STC). This condition is defined as in the temperature of 25 ° C, the irradiance of 1000 W / m 2 , and the air mass of 1.5.
In order to estimate effectively the performance of the PV module, the resolution is divided into two steps: Estimating five parameters of the SDM model (the ideality factor of diode a, the PV current I p v , the reverse saturation current I s a t , the series resistance R s , the shunt resistance R s h ) at the STC, then scaling this model at other T-G conditions of cell temperature and solar irradiance. In these reviewed articles, all methods provide three equations by applying Equation (4) in the open-circuit, short-circuit, and MPP, giving the results:
  • At short-circuit point ( I = I s c , V = 0 ):
    I s c = I p v I s a t exp R s I s c a N s V t 1 I s c R s R s h
  • At open-circuit point ( I = 0 , V = V o c ):
    0 = I p v I s a t exp V o c a N s V t 1 V o c R s h
  • At the MPP ( I = I m p p , V = V m p p ):
    I m p p = I p v I s a t exp V m p p + R s I m p p a N s V t 1 V m p p + R s I m p p R s h
Three Equations (6)–(8) are utilized in all reviewed methods. To obtain five model parameters, two additional equations or conditions are required. Since solving five transcendental equations contains the exponential functions and five variables, it cannot be solved straight away. Different approaches identifying parameters are presented in seven methods as follows.

3.1. Celik and Acikgoz Method-2007

Celik and Acikgoz [16] introduced an analytical method to extract five parameters of the SDM. The authors use two following definitions to estimate the values of the series and shunt resistances:
R s o = d V d I I = 0 ; V = V o c
1 R s h o = d I d V I = I s c , V = 0 = I s a t a N s T exp V o c a N s T + 1 R s h 1 + R s I s a t a N s T exp V o c a N s T + 1 R s h 1 R s h o = d I d V I = I s c ; V = 0
where R s o and R s h o are the reciprocals of the slopes at the open-circuit point and short-circuit point, respectively. Evaluation of the series and shunt resistances is not provided by manufacturers. Since manufacturers do not provide the information to evaluate the series and shunt resistances, these data have to be graphically extracted from the I-V curve in the datasheet, which may lead to measurement error affecting the results. After that, the following five equations, Equations (11)–(15), are proposed to calculate model parameters but do not refer to any explanations about these equations.
R s h = R s h o
a = V m p p + I m p p R s o V o c V t N s ln I s c V m p p R s h I m p p ln I s c V o c R s h + I m p p I s c V o c R s h
I s a t = I s c V o c R s h e V o c a N s V t
R s = R s o a N s V t I s a t e V o c a N s V t
I p v = I s c 1 + R s R s h + I s a t e I s c R s a N s V t 1
To scale the model at other temperature-irradiance (T-G) conditions, the values of I s c and V o c are calculated as follows:
I s c G , T = I s c G G r e f + K i T T r e f
V o c G , T = V o c + a N s V t ln G G r e f + K v T T r e f
As a result of Equations (12)–(15), four other parameters, except for R s h , can be calculated at the new T-G condition.

3.2. Villalva et al., 2009

An iterative method was used by Villalva et al. to extract five parameters. The process is started with the initial values of series and shunt resistances as follows:
R s = 0
R s h = R s h , min = V m p p I s c I m p p V o c V m p p I m p p
This process is stopped when it reaches the pre-defined tolerance of maximum power ( t o l p ). Substituting the values of the current, voltage, and power at the MPP ( I m p p , V m p p , P m a x , e ) on Equation (8), the shunt resistance can be rewritten by the relationship with series resistance as Equation (20). This equation is used to update the value of the shunt resistance based on the value of the series resistance for each loop.
R s h = V m p p V m p p + I m p p R s V m p p I p v V m p p I s a t e V m p p + I m p p R s N s a V t + V m p p I s a t P max , e
In Equation (6), the second term is assumed to be able to be neglected since it is too low compared to the PV current [30,31]. Therefore, the PV current and the reverse saturation current are, respectively, calculated as the following equations.
I p v = R s + R s h R s h I s c
I s a t = I p v V o c R s h e V o c a N s V t 1
In the paper published in 2009 [22], the author chose the initial value of the ideality diode factor equal to 1.3. However, in the updated version in 2014 on the website of the author, a is calculated by the following equation:
a = K v V o c T N s V t K i I p v 3 T E g a p k . T 2

3.3. Femia 1 et al.-2012

The value of the shunt resistance is considered high. By assuming R s h infinity, the last term in Equation (4) can be neglected. Equation (4) now becomes:
I = I p v I s a t e V + I . R s N s a V t 1
At the short-circuit condition, Equation (24) becomes:
I s c = I p v I s a t e I s c R s a N s V t 1
Since I p v , r e f I s a t , r e f e I s c R s a N s V t 1 , the PV current could be assumed equal to the short-circuit current in the STC:
I p v , r e f = I s c , r e f
The diode ideality factor is calculated by means of the following equation, which is explained in detail in [14].
a = K v V o c T N s V t K i I p v 3 T E g a p k T 2
The reverse saturation current is calculated by rewriting of Equation (24) under open-circuit condition with the assumption that e I s c R s a N s V t 1 .
I s a t = I p v e V o c a N s V t
so that the temperature coefficient C in Equation (3) can be evaluated as follows:
C = I s a t T 3 e E g a p k T
The series resistance can be determined by using the MPP data:
I m p p = I p v I s a t e V m p p + I m p p R s a N s V t 1 I p v I s a t e V m p p + I m p p R s a N s V t
Substituting Equation (28) into Equation (30) yields:
I m p p = I p v I p v e V o c + V m p p + R s I m p p a N s V t
Thus the series resistance can be calculated by means of the following equation:
R s = N s a V t ln 1 I m p p I p v + V o c V m p p I m p p
Equations (26)–(28) and (32) allow calculation of the four values of the unknown parameters ( I p v , I s a t , a, R s ) in Equation (4).

3.4. Femia 2 et al.-2012

In this method, based on the change of the variable x = V m p p + R s I m p p a V t , the series and the shunt resistance can be rewritten as functions of x as follows:
R s = x a V t V m p p I m p p ; R s h = x a V t I p v I m p p I s a t ( e x 1 )
The slope of the P-V curve at the MPP is equal to zero.
P V m p p = 0 V I V m p p = 0 I m p p + V m p p I V m p p = 0
In Equation (4) derivating the output current respect to the output voltage under the MPP gives the following equation:
I m p p + V m p p I V m p p = I m p p V m p p 1 R s h + I s a t a V t e V m p p + I m p p R s a V t 1 + R s R s h + R s I s a t a V t e V m p p + I m p p R s a V t = 0
By neglecting the small quantity terms, substituting Equation (33) into Equation (35) gives:
2 V m p p I m p p I p v I s a t + I p v + I s a t a V t x + I s a t e x a V t x + V m p p 2 V m p p a V t = 0
Simplifying the first two terms in Equation (36) gives the following equation, which is a function of x:
V m p p 2 I m p p I p v I s a t + I s a t e x a V t x + V m p p 2 V m p p a V t = 0
Based on the means of the Lambert W function, which is the solution of the equation f x = x e x :
x = l a m b e r t W V m p p ( 2 I m p p I p v I s a t ) e V m p p ( V m p p 2 a V t ) a 2 V t 2 a I s a t V t + 2 V m p p a V t V m p p 2 a 2 V t 2
The value obtained by Equation (38) is substituted into Equation (33) so that the values of the shunt and series resistances result. Three other parameters can be easily extracted using Equations (6)–(8).

3.5. Brano et al.-2010

Equation (4) is written in a different form, which expresses the five parameters’ dependency on the solar irradiance and the cell temperature.
I α G , T = α G I p v T I s a t α G , T e α G V + K T T r e f + I R s α G a N s T 1 α G V + K T T r e f + I R s R s h
where α G = G G r e f . When G = G r e f and T = T r e f , Equation (39) will become Equation (4). In order to estimate five parameters at the STC, two Equations (40) and (41) are used to find R s h o and R s o . With five non-linear equations, an iterative method was used to find five parameters. This procedure is started with two approximations: I p v , r e f = I s c , r e f and R s h = R s h o
After initializing a and R s , three other parameters I o , I p v , and R s h are estimated based on Equations (6), (8) and (40). After that, a is recalculated by Equation (7) until satisfying the pre-defined tolerance. Then, R s is evaluated by means of Equation (41). This double process is repeated until both R s and a reach the desired accuracies.
d I d V I = 0 , V = V o c = 1 R s o = I s a t a N s T exp V o c a N s T + 1 R s h 1 + R s I s a t a N s T exp V o c a N s T + 1 R s h
d I d V I = I s c , V = 0 = 1 R s h o = I s a t a N s T exp V o c a N s T + 1 R s h 1 + R s I s a t a N s T exp V o c a N s T + 1 R s h
In other T-G conditions, the reverse saturation current I s a t can be calculated by using the relationships:
I s a t α G , T = α G I p v , r e f + K i ( T T r e f ) V o c α G , T R s h e V o c α G , T a N s T 1
I s a t α G , T can be approximated by means of:
I s a t α G , T = e α G 0 . 2 1 0 . 2 ln I s a t 1 , T I s a t 0 . 2 , T + ln I s a t 0 . 2 , T

3.6. Cubas et al.-2014

Accompanied by three remarkable equations, Equations (6)–(8), the diode ideality factor is estimated to be the same as the one of the Villalva et al. method (2009). The last equation to obtain the final parameter is the differentiation of power with respect to the voltage at the MPP equal to zero. Based on some simplifications and assumptions, five parameters can be found by the following equations:
a V t V m p p 2 I m p p I s c V m p p I s c + V o c I m p p I s c V m p p I m p p R s a V t V m p p I s c V o c I m p p = exp V m p p + I m p p R s V o c a V t
R s h = V m p p I m p p R s V m p p R s I s c I m p p a V t V m p p I m p p R s I s c I m p p a V t I m p p
I s a t = R s h + R s I s c V o c R s h exp V o c a V t
I p v = R s h + R s R s h I s c
Equation (44) is an implicit expression of the series resistance. It can be solved by a solved function in Matlab. The way to scale SDM’s parameters at other T-G conditions is not mentioned in Cubas’s paper.

3.7. Laudani et al., 2014

In Laudani et al., 2014, five parameters can be classified into two groups: Three dependent variables G s h , I s a t , and I p v and two independent variables R s and a. G s h , I s a t , and I p v are rewritten as functions of R s , a .
G s h = E x p o c ( I m p p I s c ) + E x p m p p I s c E x p s c I m p p A 1 E x p s c + A 2 E x p m p p + A 3 E x p o c
I s a t = V o c ( I s c I m p p ) V m p p I s c A 1 E x p s c + A 2 E x p m p p + A 3 E x p o c
I p v = I s c V o c ( E x p m p p 1 ) + I s c V m p p ( 1 E x p o c ) + I m p p V o c ( 1 E x p s c ) A 1 E x p s c + A 2 E x p m p p + A 3 E x p o c
where V t = k T q ; E x p s c = exp I s c R s N s a V t ; E x p o c = exp V o c N s a V t ; E x p m p p = exp V m p p + I m p p R s N s a V t ; P 1 = V m p p I m p p ; P 2 = ( V o c V m p p ) I m p p ; P 3 = ( V o c V m p p ) ( I s c I m p p ) ; P 4 = V m p p ( I s c I m p p ) ; A 1 = V m p p + R s I m p p V o c ; A 2 = V o c R s I s c ; A 3 = R s I s c R s I m p p V m p p ; G s h = R 1 s h ;
The differentiation of power with respect to the voltage at the MPP at the STC is used.
I m p p V m p p = I s a t N s a V t E x p m p p + G s h 1 + R s I s a t N s a V t E x p m p p + G s h R s
By making some simplifications and substituting Equations (48)–(50) into Equation (51), the following condition is formulated:
f 1 ( R s , a ) = ( P 2 P 1 ) E x p s c + ( P 1 P 4 ) E x p o c + ( ( P 1 P 3 ) I m p p R s V m p p a N s V t + ( P 4 P 2 ) ) E x p m p p = 0
By using the Lambert fuction, R s max is expressed as a function of a:
R s max a = V m p p I m p p + N s a V t I m p p 1 + W 1 exp V o c a N s V t 2 V m p p N s a V t
There is a value of a max , which is defined by the intersection between the R s curve (by solving f 1 ) and the R s max curve (according to Equation (53)). The authors suggested using the heuristic rule: a = 0 . 9 × a max . After finding a, the R s is computed by numerically solving Equation (53). As a result, the values of I s a t , I p v , R s h are calculated by means of Equations (6)–(8) starting from the R s and a obtained previously.

4. Discussion on Reviewed Modules

4.1. Categorize Methods

Based on the approaches to obtaining the model parameters, seven methods can be categorized into three types, the analytical method, the iterative method, and the numerical method. These methods are easy to program in Matlab, using only available data in datasheets from manufacturers. However, each method has advantages and disadvantages.
In Villalva et al. [22] and Brano et al. [17], iterative methods were employed. To start the loop, these methods require to calculate some initial parameters. These values usually are approximately estimated by making assumptions. All assumptions, coupled with initial values for each model, are summarized in Table 1. In the Brano method, finding two initial calculations ( R s o , R s h o ) might be difficult since the number of points in the I-V curve chosen to calculate two values has a significant impact on the model’s result and model divergence. In the Villalva method, the diode ideality factor in Equation (23) is estimated based on the relationship with bandgap energy, which is unavailable for some PV panels. This equation, therefore, is not applicable to other PV panels, such as CIGS or poly-junction ones. Generally speaking, the iterative method gets the compromise between the algorithm complexity and the model accuracy since it is simpler than the numerical method and more accurate compared to the analytical one.
On the other hand, Celik et al. [16], Femia 1 [14] and Cubas [24] used the analytical method, posing a set of assumptions to formulate equations to calculate model parameters. These equations usually can be solved straight away, reducing the calculation time. However, since they come with assumptions to simplify equations, the model accuracy also reduces. Assumptions are listed in Table 1.
When numerical method comes into play, mathematical functions are introduced to solve the dependent relationship between current output and voltage output, specifically the Lambert-W function in [14,32]. Because the SDM applying the numerical method is solved without any assumptions thus they are more accurate at the expense of calculation time.

4.2. When Changing from One to Another Type of PV Panel

One thing that should be considered when applying these methods in applications is how effective they are when changing from one to another type of PV panel. Theoretically speaking, since each characteristic of the PV cell is different, the model parameters consequently change. The series, shunt resistances, and diode ideality factors are considered to manipulate the model [33]. Two other parameters, the reverse saturation current and the PV current, are assumed to be dependent on the three aforementioned parameters. Because the approach of each method is different, the model performance is also impacted. In Celik et al. [16] and Brano et al. [17], the diode ideality factor is incorporated in each loop, which fine-tunes the model precision. Two initial calculations, Equations (11) and (14), take into account the resistances of the model for other kinds of PV cells. In this way, the model performance is more effective for each type of PV panel. However, in calculating Equations (11) and (14), as mentioned above, the number of points of the I-V characteristic has to be chosen. Therefore, these calculations might lead to algorithm divergence when applying to some types of PV panel.
In the Cubas et al. [24] and Laudani et al. [32] method, the diode ideality factor is guessed at first. The series and shunt resistances are calculated by specific equations, which cannot self-revise when changing from one to another PV panel type. Femia et al. 1 and Femia et al. 2 [14] introduce an equation estimating the diode ideality factor based on the bandgap energy of PV panel material, which is not available for some types of PV panels. As a result, these values have to be predicted in programming.
In Villalva et al. [22], by updating the value of series and shunt resistances after a loop, these values are utilized to fit model performance compared to the provided one. The diode ideality factor is suggested as measuring the same as the Femia et al. 1 and Femia et al. 2 [14] method, which is not possible when the material of the PV panel is not silicon crystalline.

5. Proposed Method

Herein, we propose an iterative method, which takes into account the change of model parameters when changing PV panel type and does not require any initial calculations. Since the second term in Equation (6) can be neglected [30,31], the PV current is rewritten as:
I p v = I s c R s + R s h R s h
By substituting Equation (54) into Equation (8), the shunt resistance can be expressed as a function of the series resistance, as in:
R s h = I s c R s V m p p I m p p R s I m p p + I s a t exp V m p p + I m p p R s a N s V t 1 I s c
From Equation (55), the maximum value of the series resistance is calculated by assuming the denominator of the right side of Equation (55) is zero. So, it is expressed by the following equation:
R s , max = a V t ln I s c I m p p I s a t 1 V m p p I m p p
P max , e = V m p p I m p p = V m p p I p v I s a t exp V m p p + I m p p R s a N s V t 1 V m p p + R s I m p p R s h
While Equation (55), can be used as a way to update the value of the shunt resistance for each loop, it contains the simplified Equation (54), contributing the model error. By rearranging Equation (57), the shunt resistance can be calculated in an alternative way as in Equation (58):
R s h = f R s = V m p p V m p p + I m p p R s V m p p I p v V m p p I s a t exp V m p p + I m p p R s a N s V t + V m p p I s a t P m a x , e
The iterative process starts by R s = 0 ; R s h = f R s . After that, the PV current can be calculated from Equation (54) by rearranging Equation (6).
I p v = R s + R s h R s h I s c + I s a t exp R s I a N s V t 1
From Equation (8), the diode ideality factor is calculated by means of Equation (60). Finally, the reverse saturation current is calculated by Equation (61).
a n = V m p p + I m p p R s V t log I p v I m p p V m p p + I m p p R s R s h I s a t + 1
I s a t = I p v V o c R s h e V o c a N s V t 1
This process is continued with the value of the series resistance ranging from 0 ; R s , max . For each of the values of series resistance, four other parameters will be calculated, respectively, by Equations (55), (54), (60) and (61).
To improve the certainty of the algorithm, it is necessary to determine a stopping condition, which does not require any initial guess. For two points on the P-V curve, one on the left side and the other on the right side, there is a five-parameter. When substituting this five-parameter into Equation (62), which expresses the slope of the P-V curve, this value with respect to the point at the left side is greater than zero and for the point at the right side would be smaller than zero.
d P d V = d ( V I ) d V = I + V d I d V = I V I s a t a N s V t exp V + I R s a N s V t + 1 R s h 1 + R s I s a t a N s V t exp V + I R s a N s V t + R s R s h
As can be seen in Figure 2, when the series resistance increases, the peak of the P-V curve will shift from the right side of this fixed point to the left side. Consequently, the derivative of the power with respect to the voltage, as expressed in Equation (62), will monotonically decrease from positive to negative. As a result, the process stops when the condition Equation (63) is satisfied.
d P d V i d P d V i 1 < 0
Using the interpolation method, the value of R s is calculated as Equation (64).
R s = R s , i + R s , i R s , i 1 d P d V i 1 d P d V i 1 d P d V i
where i is the number of point, in which the derivative of the power with respect to the voltage is greater than zero. From Equations (58)–(61) and (64) five model parameters are identified.
Finally, the procedure is described by a flowchart in Figure 3.
Since PV systems operate at other T-G conditions further from the STC in reality, it is necessary to scale the SDM at other weather conditions. However, only some reviewed methods proposed the way to scale the model at operation conditions. In the original papers, only Brano et al. in 2010 [17] and Celik et al. in 2007 [16] proposed the way to the scale model at other T-G conditions other than the STC. However, in the Celik model [16], the authors only estimated two values of the open-circuit voltage and short-circuit current, using Equations (15) and (16). Since some datasheets do not provide the coefficient of the MPP, the current and voltage at the MPP cannot be scaled, which also means the diode ideality factor cannot be estimated. In this case, five Equations (11)–(15) are unable to scale the model. So the method to scale the model in [28] is defaulted for all models to verify model parameters at other weather conditions. The series and shunt resistance are considered to vary in inverse linear mode with solar irradiance, as expressed in Equations (65) and (66). The diode ideality factor is assumed unchanged from its value at the STC. Two remaining parameters, I s a t and I p v , are estimated by Equations (61) and (68). Since all methods use one way to scale the models at other levels of the solar irradiance and the cell temperature, the model parameters at the STC have a significant impact on the model performance at these conditions.
R s G , T = R s , r e f G r e f G
R s h G , T = R s h , r e f G r e f G
a = a r e f
I p v = I p v , r e f + K i T T r e f G G r e f

6. Numerical Results

6.1. Investigated Models

In this article, seven reviewed models and one proposed model are investigated exhaustively on monocrystalline (Shell SQ150-PC), multicrystalline (Kyocera 175GHT-2) and heterojunction (Sanyo HIT 240HDE4) PV modules. For each PV cell type, seven methods are implemented in the conditions, the cell temperature and solar irradiance of which vary.
Table 2 presents data for evaluating the model parameters provided by manufacturers only at the STC. These data are slightly different from the tabular data provided by manufacturers in the datasheets because they are graphically extracted from the characteristic voltage-current (I-V)curves in the datasheets.

6.2. Accuracy Validation

All models are verified in their performances by both metric forms and graphs. A set of graphs are depicted to validate the models’ precisions. Besides, for each graph, all results regarding the performance metric of the model are reported. The performance metric, namely the root mean squared error (RMSE), is defined as follows:
R M S E = 1 i 1 i I m , i I e , i 2
where I m , i and I e , i are the measured current values of the model, and the actual current values, respectively, and i is the total number of I-V pairs in the experimental and measured I-V curves.

6.3. Models Performances

Figure 4 depicts the I-V characteristic curves, including the curves provided by manufacturers and eight estimated curves of monocrystalline Shell SQ150-PC PV panel at various levels of solar irradiance and cell temperatures, respectively. The data show that eight estimated curves highly agree with the original curves. In Figure 4a, all curves almost overlap; meanwhile, in Figure 4c there is still a little inaccuracy after the bends of the curves.
For detailed information, Figure 4b,d reports the RMSE between individual methods in two cases, ranging solar irradiance and ranging cell temperature. In Figure 4b, from 200 W / m 2 to 800 W / m 2 , Laudani’s method reaches the highest accuracy, with the RMSEs ranging from 0.011 to 0.041. Meanwhile, Celik’s method has the lowest accuracy, with the RMSEs ranging from 0.029 to 0.080, while Laudani’s method has the highest accuracy when scaling at other conditions; at the STC, its RMSE is greater than Cubas’s and the proposal’s, which indicates its drawback when applying to the STC. One other remarkable point that should be taken into account is that the accuracies of all methods increase when applying at other levels of cell temperature and solar irradiance. As a consequence of the dependence of the explicit expressions of Ipv, I sat, a, Rs on the value of Isc and Voc, identifying these two values is vitally essential to scale model at other T-G conditions. However, two equations Equations (16) and (17), do not have the desired accuracies [17]. So the graphical values extracted from the manufactured datasheet are employed. Accordingly, the accuracy of the model increases at other T-G conditions. Figure 4d shows the exactnesses of eight methods at the solar irradiance of G = 1000 W / m 2 and cell temperature ranging from 20 ° C to 60 ° C . As can be seen in Figure 4d, Cubas’s method shows the best performance compared to other methods, with RMSEs ranging from 0.016 to 0.053. From 20 ° C to 40 ° C , Cubas’s method is the most accurate method and Celik’s method and Brano’s method are the least accurate ones, while at 50 ° C and 60 ° C the RMSEs of Cubas’s method are higher than the RMSEs of the proposed method and Villalva’s, Femia 2’s, and Laudani’s methods; the differences are not large. The effectiveness of the model increases from 20 ° C to 60 ° C .
For the Kyocera KC175GHT-2 PV panel, Figure 5a illustrates the precisions of methods at 25 ° C and various levels of solar irradiance. In Figure 5a, all of the predicted curves highly agree with the issued ones. Meanwhile, the estimated curves and issued curves at 1000 W / m 2 and various levels of cell temperature are plotted in Figure 5c. In Figure 5c, some disagreements still happen after the MPP when applying some methods to obtain parameters of the SDM.
As can be seen in Figure 5b, at high levels of solar irradiance (1000 W / m 2 , 800 W / m 2 and 600 W / m 2 ) the proposed model has the lowest RMSEs. However, at lower levels of solar irradiance (200 W / m 2 and 400 W / m 2 ) Cubas’s model and Laudani’s model have the best accuracy. In Figure 5c, generally speaking, after the MPP, eight estimated curves have agreements from the maker curves at the MPP area. As observed in Figure 5d, although Villalva’s model has the lowest errors at the STC, its errors are high at 50 ° C and 75 ° C . The proposed model shows high accuracy when RMSEs range from 0.06 to 0.1. Since the proposal’s model has low and stable RMSEs at various cell temperatures, it shows the effectiveness when scaling in other levels of cell temperature.
Figure 6 represents the performances of the SDM applied eight methods on Sanyo HIT-240HDE4 PV panel. In Figure 6a, all reviewed models show a large disagreement compared to original curve at the MPP area. This might be explained by the difference in property of this type of PV panel. In order to improve the exactness of the model, in the proposed method, the value of the diode ideality factor is modified. The data in Figure 6a shows the effectiveness of the proposed model compared to others.
As in single-crystalline and multi-crystalline PV cells, when applying eight methods to obtain I-V curves of Sanyo HIT-240HDE4 at G = 1000 W / m 2 and cell temperature ranging from 0 ° C to 75 ° C , some methods show inaccuracies after the MPP, as reported in Figure 6c.
In Figure 6b, the proposal shows the best performance since the RMSEs range from 0.021 to 0.176. While at the STC, the RMSE of proposed method is a bit higher than the RMSE of Laudani’s method, its RMSEs are lower compared to RMSEs of other methods at other levels of solar irradiance.
As can be seen from Figure 6d, Laudani’s method achieves the best performance at 25 ° C , 50 ° C , 75 ° C since the RMSEs range from 0.053 to 0.144. At 0 ° C , the Brano’s method reaches highest accuracy. This can be explained by applying Equations (40) and (41), return values of R s o and R s h o , which lead to divergence of algorithm. So we suggest a couple value of ( R s o , R s h o ) to execute the algorithm of this method.
Shortly, when applying eight methods to identify parameters of the SDM on three types of PV panels, there is not the best model that can achieve the highest accuracy in all conditions. Since the RMSEs of the reviewed models range drastically when the levels of solar irradiance and cell temperature change, they have some drawbacks. While the proposed model does not reach the highest accuracy in all conditions, the difference between its RMSE and the best method’s RMSE in each condition is not large. In addition, its RMSEs are more stable, which indicates its high reliability. Therefore, an investigation is conducted to validate the performance of the proposed model with the measured data.

6.4. Experimental Validation of Proposal Model

In order to verify the performance of the proposed model at operation conditions, an experiment is conducted on three types of PV panels, that are single-crystalline silicon Cocoa xSi12922, multi-crystalline silicon Cocoa mSi0166, and amorphous crystalline (HIT) Cocoa HIT05667. This experiment is conducted at Cocoa, Floria, USA by NREL. The data is available at the NREL website [36]. Since the cell temperature does not change much, the solar irradiance has the most significant impact on the model performance. The estimated curves are performed in the following cases:
  • Single-crystalline silicon PV panel
    As can be seen in Figure 7, for the single-crystalline silicon Cocoa xSi12922 PV panel, the predicted curves have high agreements with the actual curves. At high levels of solar irradiance and cell temperature (G = 603.8 W / m 2 ; T = 24.8 ° C to G = 1030.2 W / m 2 ; T = 33.3 ° C ), the proposal curves overestimate the output currents. On the other hand, at lower levels of solar irradiance and cell temperature (G = 35.3 W / m 2 ; T = 18.9 ° C to G = 459.9 W / m 2 ; T = 29.5 ° C ), the proposal curves underestimate the output currents.
  • Multi-crystalline silicon PV panel
    For the multi-crystalline silicon Cocoa mSi0166 PV panel in Figure 8, the model underestimates the current output. At low levels of solar irradiance and cell temperature, the disagreements between predicted curves and issued curves are larger, which can be explained by the uncertainties of experimental data tending to be bigger at low levels of T-G conditions. The predicted curves at high solar irradiance (867.2 W / m 2 and 1030.8 W / m 2 ) have inaccuracies with actual ones.
  • HIT
    For the amorphous silicon (HIT) Cocoa HIT05667 PV panel in Figure 9, there are disagreements between estimated curves and issued curves; in particular, it tends to extend after the MPP.
    In Figure 7 and Figure 8, the RMSEs of the model are lowest at irradiances (400 W / m 2 to 600 W / m 2 ) and highest at 800 W / m 2 to 1000 W / m 2 . There are two reasons for this tendency. The effectiveness of the SDM applying the proposed method is significantly affected by applying Equations (16) and (17). Consequently, in Figure 7, Figure 8 and Figure 9, the SDM shows more inaccuracies in open-circuit voltage. The RMSEs of the SDM at the solar irradiance from 30 W / m 2 to 270 W / m 2 is higher than the RMSEs of the SDM at the solar irradiance from 400 W / m 2 to 600 W / m 2 because the uncertainty of measuring instruments is higher when measuring at low levels of solar irradiance.
    As can be seen in Figure 9, the RMSEs are smaller at low levels of solar irradiance and cell temperature. The RMSEs tend to be higher at high levels of solar irradiance (635.1 W / m 2 to 1031.3 W / m 2 ) when they range from 0.2 to 0.3. At low levels of solar irradiance (35.5 W / m 2 to 1687 W / m 2 ) the RMSEs are lower with the RMSEs ranging from 0.05 to 0.1.
    The tendency of RMSEs in Figure 7 and Figure 8 is not as same as in Figure 9 because in the HIT PV panel, the fill factor is different from the two aforementioned PV types. The error of Equations (16) and (17) contributing to the SDM error surpasses the error of the uncertainty of measurement at low levels of solar irradiance.
In the above investigation, the model is tested by fixed levels of solar irradiance and cell temperature and various values of the external load, R L . For given values of solar irradiance, cell temperature, and external load, the I-V characteristics of PV panels present the function I = f ( V ) and the external load characteristic I = V / R L [17]. To make a more comprehensive validation, the proposed method is implemented for particular T-G conditions and constant values of the external load, which are characterized by 3000 pairs of ( T , G ). Figure 10 performs the predicted points and actual points for the Cocoa xSi12922 PV panel, the Cocoa mSi0166 PV panel, and the Cocoa HIT05667 PV panel, separately. The data show that the proposed method underestimates the output current for three types of panel. In addition, the simulation results tend to be more precise at lower values of T and G. The RMSE is 0.3544, 0.25, and 0.2888 for Cocoa xSi12922 PV panel, Cocoa mSi0166 PV panel, and Cocoa HIT05667 PV panel, respectively.
Four remarkable values, the open-circuit point voltage, the short-circuit current, the voltage, and current at the MPP, play the main role in obtaining the SDM as they provide the input for solving three main equations, Equations (6)–(8). In fact, expressions like Equations (16) and (17) cannot achieve the desired model precision. In addition, the coefficient of maximum power is not provided in the datasheets of some PV types; thus, the voltage and current at the MPP cannot be scaled. So, the SDM cannot be solved at other T-G conditions by the algorithm applied for solving the SDM at the STC; instead, a set of equations, Equations (61) and (65)–(68) are used to scale for each model parameter. This way might lead to the uncertainty of methods and reduce model accuracy.

7. Conclusions

In this article, seven methods of identifying parameters using the SDM have been reviewed and discussed comprehensively. The advantages and disadvantages of each method are analyzed. While some methods present good results in estimating the I-V curves in some cases, their reliability is not high when changing to other types of PV panels. A new method is proposed, which shows consistent performance for all investigated PV panels under various conditions of cell temperature and solar irradiance. This new method is verified in three types of PV panels, which are measured by NREL, USA. For three types of PV panels, single-crystalline silicon, multi-crystalline silicon, and amorphous silicon, the experimental results show RMSEs below 0.36. This error is an acceptable tolerance in practical applications requiring high accuracies such as PV power forecast and MPPT. Accurate I-V characteristics of PV panels contribute to the success in forecasting the maximum power point in the MPPT algorithm and the output power in PV power forecast.

Author Contributions

T.N.-D., H.N.-D., T.L.-V. and H.T. conceived of the methodology. T.N.-D. and T.L.-V. developed the theory and performed the computations. T.N.-D. led the research efforts and led the preparation of this paper. All authors provided critical feedback and helped shape the research, analysis and manuscript. Besides, all authors discussed the results and contributed to the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Hanoi University of Science and Technology grant number T2018-PC-058.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

adiode ideality factor [-]
a r e f diode ideality factor at the Standard Test Condition (STC) [-]
a m a x maximum value of diode ideality factor [-]
E g a p bandgap energy of the semiconductor material [J]
Gsolar irradiance [ W / m 2 ]
G r e f solar irradiance at the STC: 1000 [ W / m 2 ]
Icurrent generated by the PV modules [A]
I d Shockley diode current [A]
I m p p current at the maximum power point (MPP) [A]
I p v photovoltaic current [A]
I p v , r e f photovoltaic current at the STC [A]
I s a t reverse saturation current [A]
I s c short-circuit current [A]
I s c , r e f short-circuit current at the STC [A]
I s c ( G , T ) short-circuit current at other cell temperature (T)- solar irradiance (G) conditions (T-G conditions) [A]
kBoltzmann constant: 1 . 381 × 10 23 [ J / K ]
K i thermal coefficient of the short-circuit current [ A / o C ]
K v thermal coefficient of the open-circuit voltage [ V / o C ]
N s number of series-connected cells [-]
Ppower of the PV module [W]
P m a x , e experimental maximum power of the panel [W]
qelectron charge: 1 . 60217646 × 10 19 [ C ]
R s h shunt resistance [ Ω ]
R s h , m i n minimum shunt resistance [ Ω ]
R s h ( G , T ) shunt resistance at other levels of the cell temperature and solar irradiance [ Ω ]
R s h o reciprocal of the slope of the current-voltage (I-V) characteristic of the panel for V = 0 and I = I s c [ Ω ]
R s series resistance [ Ω ]
R s ( G , T ) series resistance at other levels of the cell temperature and solar irradiance [ Ω ]
R s o reciprocal of the slope of the I-V characteristic of the panel for V = V o c and I = 0 [ Ω ]
t o l p the pre-defined tolerance of maximum power at the STC [-]
T r e f temperature at the STC: 298 . 15 [K]
Tcell temperature [K]
Vvoltage generated by the PV modules [V]
V m p p voltage at the MPP [V]
V o c open-circuit voltage [V]
V o c ( G , T ) open-circuit voltage at other T-G conditions [V]
V t thermal voltage of the diode [V]

References

  1. IEA. Tracking Power; IEA: Paris, France, 2019; Available online: https://www.iea.org/reports/trackingpower-2019 (accessed on 20 January 2020).
  2. IRENA. Future of Solar Photovoltaic: Deployment, Investment, Technology, Grid Integration and Socio-Economic Aspects (A Global Energy Transformation: Paper); International Renewable Energy Agency: Abu Dhabi, UAE, 2019. [Google Scholar]
  3. Bader, S.; Ma, X.; Oelmann, B. One-diode photovoltaic model parameters at indoor illumination levels—A comparison. Sol. Energy 2019, 180, 707–716. [Google Scholar] [CrossRef]
  4. El-Naggar, K.M.; AlRashidi, M.R.; AlHajri, M.F.; Al-Othman, A.K. Simulated Annealing algorithm for photovoltaic parameters identification. Sol. Energy 2012, 86, 266–274. [Google Scholar] [CrossRef]
  5. Krishnakumar, N.; Venugopalan, R.; Rajasekar, N. Bacterial foraging algorithm based parameter estimation of solar PV model. In Proceedings of the 2013 Annual International Conference on Emerging Research Areas and 2013 International Conference on Microelectronics, Communications and Renewable Energy, Kanjirapally, India, 4–6 June 2013. [Google Scholar]
  6. Zagrouba, M.; Sellami, A.; Bouaïcha, M.; Ksouri, M. Identification of PV solar cells and modules parameters using the genetic algorithms: Application to maximum power extraction. Sol. Energy 2010, 84, 860–866. [Google Scholar] [CrossRef]
  7. Ishaque, K.; Salam, Z. An improved modeling method to determine the model parameters of photovoltaic (PV) modules using differential evolution (DE). Sol. Energy 2011, 85, 2349–2359. [Google Scholar] [CrossRef]
  8. Khanna, V.; Das, B.K.; Bisht, D.; Singh, P.K. A three diode model for industrial solar cells and estimation of solar cell parameters using PSO algorithm. Renew. Energy 2015, 78, 105–113. [Google Scholar] [CrossRef]
  9. Oliva, D.; Cuevas, E.; Pajares, G. Parameter identification of solar cells using artificial bee colony optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  10. Lin, P.; Cheng, S.; Yeh, W.; Chen, Z.; Wu, L. Parameters extraction of solar cell models using a modified simplified swarm optimization algorithm. Sol. Energy 2017, 144, 594–603. [Google Scholar] [CrossRef]
  11. Chen, Y.; Sun, Y.; Meng, Z. An improved explicit double-diode model of solar cells: Fitness verification and parameter extraction. Energy Convers. Manag. 2018, 169, 345–358. [Google Scholar] [CrossRef]
  12. De Blas, M.A.; Torres, J.L.; Prieto, E.; García, A. Selecting a suitable model for characterizing photovoltaic devices. Renew. Energy 2002, 25, 371–380. [Google Scholar] [CrossRef]
  13. Wang, G.; Zhao, K.; Shi, J.; Chen, W.; Zhang, H.; Yang, X.; Zhao, Y. An iterative approach for modeling photovoltaic modules without implicit equations. Appl. Energy 2017, 202, 189–198. [Google Scholar] [CrossRef]
  14. Femia, N.; Petrone, G.; Spagnuolo, G.; Vitelli, M. PV Modeling Power Electronics and Control Techniques for Maximum Energy Harvesting in Photovoltaic Systems; CRC Press: Boca Raton, FL, USA, 2017; pp. 1–33. [Google Scholar]
  15. Hussein, A. A simple approach to extract the unknown parameters of PV modules. Turkish J. Electr. Eng. Comput. Sci. 2017, 25, 4431–4444. [Google Scholar] [CrossRef]
  16. Celik, A.N.; Acikgoz, N. Modelling and experimental verification of the operating current of mono-crystalline photovoltaic modules using four- and five-parameter models. Appl. Energy 2007, 84, 1–15. [Google Scholar] [CrossRef]
  17. Lo Brano, V.; Orioli, A.; Ciulla, G.; Di Gangi, A. An improved five-parameter model for photovoltaic modules. Sol. Energy Mater. Sol. Cells 2010, 94, 1358–1370. [Google Scholar] [CrossRef]
  18. Gow, J.A.; Manning, C.D. Development of a photovoltaic array model for use in power-electronics simulation studies. IEE Proc. Electr. Power Appl. 1999, 146, 193–200. [Google Scholar] [CrossRef]
  19. Dehghanzadeh, A.; Farahani, G.; Maboodi, M. A novel approximate explicit double-diode model of solar cells for use in simulation studies. Renew. Energy 2017, 103, 468–477. [Google Scholar] [CrossRef]
  20. Wang, G.; Zhao, K.; Qiu, T.; Yang, X.; Zhang, Y.; Zhao, Y. The error analysis of the reverse saturation current of the diode in the modeling of photovoltaic modules. Energy 2016, 115, 478–485. [Google Scholar] [CrossRef]
  21. Arab, A.H.; Chenlo, F.; Benghanem, M. Loss-of-load probability of photovoltaic water pumping systems. Sol. Energy 2004, 76, 713–723. [Google Scholar] [CrossRef]
  22. Villalva, M.G.; Gazoli, J.R.; Filho, E.R. Comprehensive approach to modeling and simulation of photovoltaic arrays. IEEE Trans. Power Electron. 2009, 24, 1198–1208. [Google Scholar] [CrossRef]
  23. Cubas, J.; Pindado, S.; Victoria, M. On the analytical approach for modeling photovoltaic systems behavior. J. Power Sources 2014, 247, 467–474. [Google Scholar] [CrossRef] [Green Version]
  24. Ghani, F.; Duke, M. Numerical determination of parasitic resistances of a solar cell using the Lambert W-function. Sol. Energy 2011, 85, 2386–2394. [Google Scholar] [CrossRef]
  25. Jain, A.; Kapoor, A. Exact analytical solutions of the parameters of real solar cells using Lambert W-function. Sol. Energy Mater. Sol. Cells 2004, 81, 269–277. [Google Scholar] [CrossRef]
  26. Ghani, F.; Duke, M.; Carson, J. Numerical calculation of series and shunt resistances and diode quality factor of a photovoltaic cell using the Lambert W-function. Sol. Energy 2013, 91, 422–431. [Google Scholar] [CrossRef]
  27. Rusirawan, D.; Farkas, I. Identification of model parameters of the photovoltaic solar cells. Energy Procedia 2014, 57, 39–46. [Google Scholar] [CrossRef] [Green Version]
  28. De Soto, W.; Klein, S.A.; Beckman, W.A. Improvement and validation of a model for photovoltaic array performance. Sol. Energy 2006, 80, 78–88. [Google Scholar] [CrossRef]
  29. Walker, G. Evaluating MPPT converter topologies using a matlab PV model. J. Electr. Electron. Eng. 2001, 21, 49–55. [Google Scholar]
  30. Kennerud, K.L. Analysis of Performance Degradation in CDS Solar Cells. IEEE Trans. Aerosp. Electron. Syst. AES-5 1969, 6, 912–917. [Google Scholar] [CrossRef]
  31. Rodriguez, P.; Sera, D.; Teodorescu, R.; Rodriguez, P. PV Panel Model Based on Datasheet Values PV panel model based on datasheet values. In Proceedings of the IEEE International Symposium on Industrial Electronics, Vigo, Spain, 4–7 June 2007; pp. 2392–2396. [Google Scholar]
  32. Laudani, A.; Fulginei, F.R.; Salvini, A. Identification of the one-diode model for photovoltaic modules from datasheet values. Sol. Energy 2014, 108, 432–446. [Google Scholar] [CrossRef]
  33. Bana, S.; Saini, R.P. A mathematical modeling framework to evaluate the performance of single diode and double diode based SPV systems. Energy Rep. 2016, 2, 171–187. [Google Scholar] [CrossRef] [Green Version]
  34. Ciulla, G.; Lo Brano, V.; Di Dio, V.; Cipriani, G. A comparison of different one-diode models for the representation of I-V characteristic of a PV cell. Renew. Sustain. Energy Rev. 2014, 32, 684–696. [Google Scholar] [CrossRef]
  35. Humada, A.M.; Hojabri, M.; Mekhilef, S.; Hamada, H.M. Solar cell parameters extraction based on single and double-diode models: A review. Renew. Sustain. Energy Rev. 2016, 56, 494–509. [Google Scholar] [CrossRef] [Green Version]
  36. Marion, W.; Anderberg, A.; Deline, C.; Glick, S.; Muller, M.; Perrin, G.; Rodriguez, J.; Rummel, S.; Terwilliger, K.; Nrel, T.J.S. User’s Manual for Data for Validating Models for PV Module Performance. 2014. Available online: https://www.nrel.gov/docs/fy14osti/61610.pdf (accessed on 6 October 2019).
Figure 1. Equivalent circuit of single-diode model (SDM) for photovoltaic (PV) cells.
Figure 1. Equivalent circuit of single-diode model (SDM) for photovoltaic (PV) cells.
Energies 13 01296 g001
Figure 2. P-V curves plotted for different values of R s .
Figure 2. P-V curves plotted for different values of R s .
Energies 13 01296 g002
Figure 3. The flowchart to obtain model parameters.
Figure 3. The flowchart to obtain model parameters.
Energies 13 01296 g003
Figure 4. Comparison between the calculated (eight models) and maker I-V curves of Shell SQ150-PC: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Figure 4. Comparison between the calculated (eight models) and maker I-V curves of Shell SQ150-PC: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Energies 13 01296 g004
Figure 5. Comparison between calculated (eight models) and maker I-V curves of Kyocera KC175GHT-2: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Figure 5. Comparison between calculated (eight models) and maker I-V curves of Kyocera KC175GHT-2: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Energies 13 01296 g005
Figure 6. Comparison between calculated (eight models) and issued I-V curves of Sanyo HIT-240HDE4: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Figure 6. Comparison between calculated (eight models) and issued I-V curves of Sanyo HIT-240HDE4: (a) at T = 25 ° C ; (b) RMSE at T = 25 ° C ; (c) at G = 1000 [ W / m 2 ]; (d) RMSE at G = 1000 [ W / m 2 ].
Energies 13 01296 g006
Figure 7. Comparison between the calculated (proposed model) and experimental curves of operation conditions of single-crystalline Cocoa xSi12922 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Figure 7. Comparison between the calculated (proposed model) and experimental curves of operation conditions of single-crystalline Cocoa xSi12922 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Energies 13 01296 g007
Figure 8. Comparison between the calculated (proposed model) and experimental curves of operation conditions of the multi-crystalline Cocoa mSi0166 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Figure 8. Comparison between the calculated (proposed model) and experimental curves of operation conditions of the multi-crystalline Cocoa mSi0166 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Energies 13 01296 g008
Figure 9. Comparison between the calculated (proposed model) and experimental curves of operation conditions of the amorphous silicon Cocoa HIT05667 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Figure 9. Comparison between the calculated (proposed model) and experimental curves of operation conditions of the amorphous silicon Cocoa HIT05667 PV panel: (a) Graphical comparison; (b) RMSE of proposed model.
Energies 13 01296 g009
Figure 10. The actual and simulated points of: (a) Cocoa xSi12922 PV panel, (b) Cocoa mSi0166 PV panel, (c) Cocoa HIT05667 PV panel.
Figure 10. The actual and simulated points of: (a) Cocoa xSi12922 PV panel, (b) Cocoa mSi0166 PV panel, (c) Cocoa HIT05667 PV panel.
Energies 13 01296 g010
Table 1. Comparison among reviewed models.
Table 1. Comparison among reviewed models.
MethodAssumptions (A)
Intial Guesses (IG), Initial Caculations (IC)
a R s R sh
CelikAnalyticalIC: R s o , R s h o
A: R s h = R s h o
f R s o , R s h o
(Equation (12)
Self-revised
(Equation (14)
Self-revised
(Equation (11)
VillalvaIterativeIG: t o l P , R s = 0
R p = R p , min = V m p I s c I m p V o c V m p I m p
A: I p v , r e f I s a t , r e f exp I R s a N s V t 1
Constant
(Equation (23))
Self-revised
(Iterative)
Self-revised
(Iterative)
Femia 1AnalyticalA: R s h =
I p v , r e f I s a t , r e f exp I R s a N s V t 1
Constant
(Equation (23))
Constant
(Equation (32))
Constant
( R s h = )
Femia 2NumericalA: R s I m p p , r e f 2 V m p p , r e f Constant
(Equation (23))
ConstantConstant
BranoIterativeA: e V o c / a N s V t 1 , R s R p , I p v I s c
I s a t a N s V t e I s c R s a N s V t 1 R p I s a t a N s V t e V o c a N s V t
IC: R s o , R p o
IG: R s = R s o , R p = R p o , a = 1
Self-revised
(Iterative)
Self-revised
(Iterative)
Self-revised
(Iterative)
CubasAnalyticalIG: a
A: I p v , r e f I s a t , r e f exp I R s a N s V t 1
Constant
(initial choice)
Const
(Equation (44))
Const
(Equation (45))
LaudaniNumericalIG: a = 0 . 9 × a max Constant
(initial choice)
Const
(Equation (53))
Const
(Equation (8))
Table 2. Data for evaluating the model parameters of Shell SQ150-PC, Kyocera 175GHT-2, and Sanyo HIT 240HDE4 at the standard test condition (STC).
Table 2. Data for evaluating the model parameters of Shell SQ150-PC, Kyocera 175GHT-2, and Sanyo HIT 240HDE4 at the standard test condition (STC).
Cell Type V oc [V] I sc [A] V mpp [V] I mpp [A] N s [-] K v [ V / ° C ] K i [ A / ° C ]
Shell SQ150-PC43.464.8233.734.4872−0.1610.0014
Kyocera 175GHT-228.568.097.4723.7148−0.1070.00222
Sanyo HIT-240HDE443.887.435.157.0560−0.1090.00221

Share and Cite

MDPI and ACS Style

Nguyen-Duc, T.; Nguyen-Duc, H.; Le-Viet, T.; Takano, H. Single-Diode Models of PV Modules: A Comparison of Conventional Approaches and Proposal of a Novel Model. Energies 2020, 13, 1296. https://doi.org/10.3390/en13061296

AMA Style

Nguyen-Duc T, Nguyen-Duc H, Le-Viet T, Takano H. Single-Diode Models of PV Modules: A Comparison of Conventional Approaches and Proposal of a Novel Model. Energies. 2020; 13(6):1296. https://doi.org/10.3390/en13061296

Chicago/Turabian Style

Nguyen-Duc, Tuyen, Huy Nguyen-Duc, Thinh Le-Viet, and Hirotaka Takano. 2020. "Single-Diode Models of PV Modules: A Comparison of Conventional Approaches and Proposal of a Novel Model" Energies 13, no. 6: 1296. https://doi.org/10.3390/en13061296

APA Style

Nguyen-Duc, T., Nguyen-Duc, H., Le-Viet, T., & Takano, H. (2020). Single-Diode Models of PV Modules: A Comparison of Conventional Approaches and Proposal of a Novel Model. Energies, 13(6), 1296. https://doi.org/10.3390/en13061296

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop