DeSoto PV

Download as pdf or txt
Download as pdf or txt
You are on page 1of 7

IMPROVEMENT AND VALIDATION OF A MODEL FOR

PHOTOVOLTAIC ARRAY PERFORMANCE



W. De Soto S.A. Klein
Solar Energy Laboratory Solar Energy Laboratory
University of Wisconsin-Madison University of Wisconsin-Madison
Madison, WI 53706 Madison, WI 53706

W.A. Beckman
Solar Energy Laboratory
University of Wisconsin-Madison
Madison, WI 53706


ABSTRACT

The 5-Parameter model described in this paper uses data
provided by the manufacturer, together with semi-empirical
correlation equations and weather data, to predict annual
energy production. Annual data obtained from a building
integrated photovoltaic facility at the National Institute of
Standards and Technology (NIST) were provided for four
different cell technologies. The 5-Parameter model and the
King [1] model are compared with the measured data
provided by NIST over a wide range of operating
conditions. Both models show excellent results.
Fig. 1: Equivalent circuit for the 5-Parameter model
Five parameters must be known in order to determine the
current and voltage, and thus the power delivered to the
load. These are: the light current I
L
, the diode reverse
saturation current I
o
, the series resistance R
s
, the shunt
resistance R
sh
, and a modified ideality factor a N
s
n
I
k T
c
/q.
The electron charge q, and Boltzmann's constant k are
known, n
I
is the usual ideality factor, N
s
is the number of
cells in series and T
c
is the cell temperature. The current-
voltage relationship at a fixed cell temperature and solar
radiation for the circuit in Fig. 1 is expressed in Eqn. (1).

1. INTRODUCTION
The electrical power output from a photovoltaic panel
depends primarily on the incident solar radiation, the cell
temperature, and solar incidence angle. Manufacturers
typically provide only limited operational data for
photovoltaic panels that include: open circuit voltage (V
oc
),
the short circuit current (I
sc
), the maximum power current
(I
mp
) and voltage (V
mp
), the temperature coefficients at open
circuit voltage and short circuit current (
Voc
and
Isc
,
respectively), and the nominal operating cell temperature
(NOCT). These data are available only at Standard Rating
Conditions (SRC), for which the irradiance is 1,000 W/m
2

and the cell temperature (T
c
) is 25C (except for the NOCT
which is determined at 800 W/m
2
, 20C). The purpose of
this paper is to demonstrate how the performance
information provided at SRC can be used to determine the
photovoltaic system performance at any specified condition.
1
s
V IR
s
a
L o
sh
V IR
I I I e
R
+

+
=


(1)
3. CALCULATING THE REFERENCE PARAMETERS

To evaluate the five parameters in Eqn. (1), five
independent pieces of information are needed. In general,
these five parameters are functions of the solar radiation
incident on the cell and cell temperature. Reference values
of these parameters are determined for a specified operating
condition such as SRC. Three current-voltage pairs are
normally available from the manufacturer at SRC: the short
circuit current, the open circuit voltage and the current and

2. CURRENT-VOLTAGE RELATIONSHIP
The electrical power available from a photovoltaic (PV)
device can be modeled with the circuit shown in Fig. 1 [2,
3].
voltage at the maximum power point. A fourth piece of
information results from recognizing that the derivative of
the power at the maximum power point is zero. Although
both the temperature coefficient of the open circuit voltage
(
Voc
) and the temperature coefficient of the short circuit
current (
Isc
) are known, only
Voc
is used to find the five
reference parameters.
Isc
is used when the cell is operating
at conditions other than reference conditions.

Substituting the known information into Eqn, (1) results into
Eqns. (2-7) which determine the circuit parameters (a
ref
,
I
o,ref
, I
L,ref
, R
s,ref
, and R
sh,ref)
at reference conditions.

For short circuit current: I=I
sc,ref
, V=0
, ,
, ,
, , ,
,
1
sc ref s ref
ref
I R
a sc ref s ref
sc ref L ref o ref
sh ref
I R
I I I e
R


=



(2)
For open circuit voltage: I=0, V=V
oc,ref


,
,
, ,
,
0
oc ref
ref
V
a oc ref
L ref o ref
1
sh ref
V
I I e
R


=



(3)
At the maximum power point: I=I
mp,ref
, V=V
mp,ref


, , ,
, , ,
, , ,
,
1
mp ref mp ref s ref
ref
V I R
a
mp ref L ref o ref
mp ref mp ref s ref
sh ref
I I I e
V I R
R
+


=

(4)
At the maximum power point,

( )
0
mp mp
mp mp
d IV dI
I V
dV dV
= = . (5)
where
mp
/ dI dV is given by:


1
1
mp mp s
mp mp s
V I R
o
a
sh
V I R
mp
o s s
a
sh
I
e
a R dI
dV
I R R
e
a R
+
+


=
+ +
(6)
The temperature coefficient of open circuit voltage is given
by:
, ,
0
c
oc ref oc T
Voc
ref c
I
V V V
T T T

=

=

(7)
To evaluate
Voc
numerically, it is necessary to know
, the open circuit voltage at some cell temperature
near the reference temperature. The cell temperature used
for this purpose is not critical since values of T
c
ranging
from 1 to 10 K above or below T
ref
provide essentially the
same result. In order to apply Eqn. (7), it is necessary to
obtain expressions for the temperature dependence of the
three parameters a, I
o
and, I
L
. The temperature dependence
of all of the parameters in the model is considered in the
following section.
,
c
oc T
V
4. DEPENDENCE OF OPERATING CONDITIONS
From the definition of a, the modified ideality factor is a
linear function of cell temperature (assuming n
I
is
independent of temperature) so that:

,
c
ref c ref
T a
a T
= (8)
where T
c,ref
and a
ref
are the cell temperature and modified
ideality factor for reference conditions, while T
c
and a are
the cell temperature and modified ideality factor parameter
for the new operating conditions.
Messenger and Ventre [4] present an equation from diode
theory for the diode reverse saturation current, I
o
. The ratio
of their equation at the new operating temperature to that at
the reference temperature yields:

3
, ,
1
exp
ref c
g g
T T
o c
o ref c ref
E E
k T T
I T
I T











= (9)
where k is Boltzmann's constant and E
g
is the material band
gap. The values of the material band gap energies at 25C
for the four cell types investigated in this study can be found
in Table A1. E
g
exhibits a small temperature dependence
[5] which, for silicon, can be represented as indicated in
Eqn. (10) where
,
ref
g T
=1.121 eV for silicon cells. Equation
(10) was used for all of the cells considered in this study
although the value of
,
ref
E
g T
used for the triple junction
amorphous cell type was 1.6 eV.
E
(
,
1 0.0002677
ref
g
ref
g T
E
T T
E
=
)
(10)
The light current, (I
L
), is nearly a linear function of incident
solar radiation. Some pyranometers in fact use the short
circuit current of a solar cell as a measure of the incident
solar radiation. The light current (I
L
) is observed to depend
on the absorbed solar irradiance (S), the cell temperature
(T
c
), the short circuit current temperature coefficient (
Isc
),
and the air mass modifier (M). The light current I
L
for any
operating conditions is related to the light current at
reference conditions by:
,
[ (
L L ref Isc c c ref
ref ref
S M
I I T
S M
= +
,
)] T (11)
where S
ref
, M
ref,
I
L,ref
, T
c,ref
are the parameters at reference
conditions, while S, M, I
L
, and T
c
are the values for specified
operating conditions. When using Eqn. (11) to find the
reference parameters, S = S
ref
and M = M
ref
. The air mass
modifier is assumed to be a function of the local zenith
angle and is discussed below.
The information needed to determine the reference
parameters is now complete. Eqns. (3-7) relate the model to
the known reference conditions. To evaluate Eqn. (7) it is
necessary to use the temperature dependence of a, I
o
and I
L

as given by Eqns. (8-11).
The final task to complete the model is to investigate the
operating condition dependence of the series resistance R
s
,
and the shunt resistance, R
sh
. The series resistance impacts
the shape of current and voltage curve near the maximum
power point. The effect on the I-V curve is small and,
although methods of adjusting R
s
as a function of operating
conditions have been investigated [6], R
s
is assumed
constant at its reference value, R
s,ref
in this study.
The shunt resistance (R
sh
) controls the slope of the I-V curve
at the short circuit condition; large shunt resistances result
in a horizontal slope. The shunt resistance appears to
change with absorbed solar radiation for all of the cells
although the effect is most noticeable for cell types that
have a relatively small shunt resistance at SRC. If
experimental data were generally available at more than one
solar radiation value, it would be possible to develop a
relation between the shunt resistance and absorbed
radiation. However, this information is not normally
available. An observation apparent from an examination of
the slopes of the I-V curves at short circuit conditions based
on the experimental data is that the effective shunt
resistance increases (and the slope thus decreases) as
absorbed radiation is reduced. Eqn (12) is proposed to
describe this effect.

,
ref sh
sh ref
S R
R S
= (12)

5. THE INCIDENCE ANGLE MODIFIER, K


The incidence angle is directly involved in the
determination of the radiation incident on the surface of the
PV device. In addition, the incidence angle affects the
amount of solar radiation transmitted through the protective
cover and absorbed by the cell. As the incidence angle
increases, the amount of radiation reflected from the cover
increases with the most noticeable effects occurring at
incidence angles greater than 65.
The effects of reflection and absorption as a function of
incidence angle are expressed in terms of the incidence
angle modifier, K

(), defined as the ratio of the radiation


absorbed by the cell at some incidence angle divided by
radiation absorbed by the cell at normal incidence. The
incidence angle can be calculated as a function of the panel
slope and location and time. Panels that are mounted on a
vertical surface exacerbate the incidence angle effects
because much of the available beam solar radiation strikes
the panel surface at angles greater than 65. Nevertheless,
vertically mounted panels are of interest because of the
applicability of this orientation to installation on building
faades. The experimental data that are available to validate
the model presented in this paper were taken on a vertical
surface.
The incidence angle modifier for a PV panel differs
somewhat from that of a flat-plate solar collector in that the
glazing may be bonded to the cell surface, thereby
eliminating one air-glazing interface and the glazing surface
may be treated so as to reduce reflection losses. Equations
(13) and (14) based on Snells and Boughers laws, as
reported in Duffie and Beckman [2], are used to calculate
the incidence angle modifier. The angle of refraction (
r
) is
determined from Snells law
( arcsin sin
r
n ) = (13)
where is the incidence angle and n is an effective index of
refraction of the cell cover. A good approximation of the
transmittance of the cover system is:
( )
( )
( )
( )
( )
( )
2 2
cos
2 2
sin tan 1
1
2
sin tan
r
r r KL
r r
e





= +

+ +

(14)
where K is the glazing extinction coefficient and L is the
glazing thickness. In this study the value of K is assumed to
be 4 m
-1
, the value for water white glass and the glazing
thickness is assumed to be 2 mm, a reasonable value for
most PV cell panels. The refractive index is set to 1.526,
the value for glass.

To obtain the incidence angle modifier (K

), Equation (15)
needs be evaluated for incidence angles of 0 and . The
ratio of these two transmittances yields K

:
( )
( )
( ) 0
K

= (15)
Separate incidence angle modifiers are needed for beam,
diffuse, and ground-reflected radiation, but each can be
calculated in the same way. Average angles for isotropic
diffuse and ground-reflected radiation are provided as a
function of the slope of the panel by Fig. 5.4.1 of Duffie and
Beckman [2]. Although these average angles for diffuse
radiation were obtained for thermal collectors, they were
found to yield reasonable results for PV systems.

6. THE AIR MASS MODIFIER, M
Air mass is the ratio of the mass of air that the beam
radiation has to traverse at any given time and location to
the mass of air that the beam radiation would traverse if the
sun were directly overhead. Selective absorption by species
in the atmosphere causes the spectral content of irradiance
to change, altering the spectral distribution of the radiation
incident on the PV panel. King et al. [7] developed an
empirical relation to account for changes in the spectral
distribution resulting from changes in air mass:

2 3
0 1 2 3 4
ref
M
a a AM a AM a AM a AM
M
= + + + +
4
(16)
In Equation (16), AM is the air mass and a
o
, a
1
, a
2
, a
3
, and a
4

are constants for different PV materials. These constants
are listed in Table A1, along with other values used for this
study. However, if one set of constants is chosen and used
for all cell types, the difference in the results obtained is
insignificant so that the values of a
i
for the poly-crystalline
cell in Table A1 can be used for all cases.

7. CALCULATING THE ABSORBED RADIATION, S
The major factor affecting the power output from a PV
device is the solar radiation absorbed on the cell surface, S,
which is a function of the amount of radiation and the
incidence angle. Radiation data are not normally known on
the plane of the PV panel, so it is necessary to estimate the
absorbed solar radiation using horizontal data and incidence
angle information. The total absorbed irradiance S consists
of beam, diffuse, and ground reflected components. The
beam radiation G
b
is the solar radiation received from the
sun that has not been scattered by the atmosphere while the
diffuse radiation G
d
is the total horizontal radiation G minus
the beam radiation. The required beam-diffuse breakup was
done using the Erbs correlation [2]. Equation (11) provides
a means of estimating the absorbed radiation, S, assuming
that both diffuse and ground-reflected radiation are isotropic
[2].
, , ,
(1 ) (1 )
2 2
b d
beam b d g
ref ref ref ref
G G S Cos G Cos
R K K K
S G G G

+
= + +
(11)

In Equation (11), is the ground reflectance, is the slope
of the panel, K
,b
is the incidence angle modifier at the
beam incidence angle, K
,d
and K
,g
are the incidence angle
modifiers at effective incidence angles for isotropic diffuse
and ground-reflected radiation, and R
beam
is the ratio of
beam radiation on a tilted surface to that on a horizontal.

8. VALIDATION OF THE 5-PARAMETER MODEL
The data used for this study were provided by Fanney et al.
[8] from a building integrated photovoltaic facility at NIST
in Gaithersburg, Maryland. The experimental data provide
at five-minute intervals one year (Jan 1-Dec. 31, 2000) of
meteorological data and measured cell temperature along
with current and voltage values for four different
photovoltaic cell technology types (single crystalline, poly-
crystalline, silicon thin film, and 3-junction amorphous)
installed on a vertical surface.
King [1] has developed a model to predict the energy
production from a PV array. This model determines the
current at five strategically located points on the current-
voltage (I-V) curve and uses this information to identify the
entire I-V curve. The accuracy of the King model is
excellent. However, its disadvantage is that it requires a
large amount of experimental data.
Figures 3-6 show results obtained with the 5-Parameter and
King models. Note that the reference parameters were
determined at the SRC condition, 1,000 W/m
2
and 25C.

Fig. 3: Current vs voltage for the single-crystalline cell type.


Fig. 4: Current vs voltage for the poly-crystalline cell type.
Fig. 5: Current vs voltage for the silicon thin film cell type.


Fig. 6: Current vs voltage for the amorphous triple junction
cell type.

Figures 3-6 show very good agreement for both the 5-
Parameter model and the King model with the experimental
data for all four cell types over a range of operating
conditions.

9. CONCLUSIONS
The 5-Parameter model presented in this paper uses only
data provided by the manufacturer with semi-empirical
equations to predict the cell I-V curve (and thus the power)
for any operating condition. The model requires a one-time
calculation of the five parameters (a
ref,
I
o,ref
, I
L,ref
,
Rs,ref
, and
R
sh,ref
) at reference conditions. These values are then used
with in the model to calculate the parameters at other
operating conditions, making it possible to predict the
power output at any operating conditions. The shunt
resistance was observed to be a function of absorbed
radiation and an empirical equation was proposed to
describe this effect.
Comparisons with experimental data provided by NIST [8]
have shown that the 5-Parameter model can be an accurate
tool for the prediction of energy production for single-
junction cell types. The 5-Parameter model uses only data
provided by the manufacturer, and in contrast to the King
model, does not require parameters that need to be
predetermined by additional experiments. The predictions
from the 5-Parameter model are shown to agree well with
both the King model results and the NIST [8] measurements
for all four cell types over a range of operating conditions.
The differences between the experimental data and the 5-
Parameter model could be reduced if additional
experimental data, e.g., I-V curves at two radiation levels,
were used to determine the reference parameters.

10. ACKNOWLEDGMENTS
The authors thank the Graduate Engineering Research
Students (GERS) for their financial support and Hunter
Fanney and Mark Davis from NIST for providing the data
used to validate the models. We especially wish to thank
Mark Davis for his help and insight. We would also like to
thank Michal Kummert for his help with data
manipulation.

11. REFERENCES
1. King, D.L., (2000), Sandias PV Module Electrical
Performance Model (Version, 2000), Sandia National
Laboratories, Albuquerque, New Mexico, September 5.
2. Duffie, J.A., and Beckman, W.A. (1991), Solar
Engineering of Thermal Processes 2nd Edition, John
Wiley & Sons, Inc., New York.,
3. Nelson, Jenny, (2003), The Physics of Solar Cells,
Imperial College Press, London.
4. Messenger, R. A., and Ventre, J. (2004), Photovoltaic
Systems Engineering 2nd Edition, CRC Press LLC,
Boca Raton, Florida.
5. Van Zeghbroeck, B., Principles of Semiconductor
Devices, (2004) http://ece-www.colorado.edu/~bart/
book/ book/chapter2/ch2_3.htm.
6. DeSoto, Widalys, (2004), M.S. Thesis, Mechanical
Engineering, U. of Wisconsin-Madison, Improvement
and Validation of a Model for Photovoltaic Array
Performance.
7. King, D.L., Kratochvil, J.A., Boyson, W.E., and Bower,
W.I., (1998) Field Experience with a New
Performance Characterization Procedure for
Photovoltaic Arrays, 2
nd
World Conference and
Exhibition on Photovoltaic Solar energy Conversion,
Vienna, Austria, July 6-10.
8. Fanney, A. H., Dougherty, B.P. and Davis, M.W.,
(2002), Evaluating Building Integrated Photovoltaic
Performance Models, Proceedings of the 29th IEEE
Photovoltaic Specialists Conference (PVSC), May 20-
24th, New Orleans, Louisiana.
9. Fanney, A.H, Davis, M.W., and Dougherty, B.P.,
Short-Term Characterization of Building-Integrated
Photovoltaic Panels,, Proceedings of the Solar Forum,
Sunrise on the Reliable Energy Economy, ASES, Reno,
NV, June 15-19, 2002.


Table A1: Values provided by NIST for the different cell types
Type of Cell
Silicon
thin film
Single-
crystalline
Poly-
crystalline
3-Junction
amorphous
At Reference Conditions
P
mp,ref
(W) 103.96 133.40 125.78 57.04
I
sc,ref
(A) 5.11 4.37 4.25 4.44
V
oc,ref
(V) 29.61 42.93 41.50 23.16
I
mp,ref
(A) 4.49 3.96 3.82 3.61
V
mp,ref
(V) 23.17 33.68 32.94 16.04
NOCT (K) 316.15 316.85 316.45 311.05
Temperature Coefficients

Isc
(A/K) 0.00468 0.00175 0.00238 0.00561

Isc
(1/K) 0.000916 0.000401 0.000560 0.001263

Imp
(A/K) 0.00160 -0.00154 0.00018 0.00735

Imp
(1/K) 0.000358 -0.000390 0.000047 0.002034

Voc
(V/K) -0.12995 -0.15237 -0.15280 -0.09310

Voc
(1/K) -0.004388 -0.003549 -0.003682 -0.004021

Vmp
(V/K) -0.13039 -0.15358 -0.15912 -0.04773

Vm
p (1/K) -0.005629 -0.004560 -0.004830 -0.002976
King Model Parameters Determined by NIST [9]. c
4
through c
7
were obtained from Sandia
http://www.sandia.gov/pv/docs/Database.htm
a
0
0.938110 0.935823 0.918093 1.10044085
a
1
0.062191 0.054289 0.086257 -0.06142323
a
2
-0.015021 -0.008677 -0.024459 -0.00442732
a
3
0.001217 0.000527 0.002816 0.000631504
a
4
-0.000034 -0.000011 -0.000126 -1.9184E-05
b
0
0.998980 1.000341 0.998515 1.001845
b
1
-0.006098 -0.005557 -0.012122 -0.005648
b
2
8.117E-04 6.553E-04 1.440E-03 7.25E-04
b
3
-3.376E-05 -2.730E-05 -5.576E-05 -2.916E-05
b
4
5.647E-07 4.641E-07 8.779E-07 4.696E-07
b
5
-3.371E-09 -2.806E-09 -4.919E-09 -2.739E-09
c
0
0.9615 0.9995 1.0144 1.072
c
1
0.0368 0.0026 -0.0055 -0.098
c
2
0.2322 -0.5385 -0.3211 -1.8457
c
3
-9.4295 -21.4078 -30.2010 -5.1762
c
4
0.967 0.9980 0.9931 1.059
c
5
0.033 0.0020 0.0069 -0.059
c
6
1.12 1.159 1.104 1.188
c
7
-0.120 -0.159 -0.104 -0.188
n
D
1.357 1.026 1.025 3.09
Other Parameters
N
s
40 72 72 22
E
g
(eV) at 25C 1.12 1.12 1.14 1.6

You might also like