Photovoltaic Arrays[LINK]
The Photovoltaics.f90 module includes three different models referred to as “Simple”, “Equivalent OneDiode” and “Sandia” and the choice will determine the mathematical models (and input data) used to determine the energy produced by solar/electric conversion panels. The EnergyPlus photovoltaic array models are called one at a time at the HVAC system timestep along with other electrical generation components such as gas turbines and diesel engines.
All of the photovoltaic models share the same models for predicting incident solar radiation that are also used for the solar thermal calculations and are described in the section Climate, Sky and Solar/Shading Calculations.
Note that some of the terminology used to discussed photovoltaics overlaps with terminology used to discuss Fortran programs. The word module may refer to a PV panel or to a fortran90 programming entity. Model may refer to a manufacturers production model for a specific type of PV module or to a mathematical model used for engineering analysis. Array may refer to a collection of PV modules wired together or to a mathematical variable with multiple elements.
The PV modules are assumed to always run when the total incident solar is greater than 0.3 Watts. If the incident solar is less than 0.3, then the modules produce no power.
PV arrays are managed by an electric load center. The load center is a “load” with respect to generating equipment but is a “supply center” for the rest of the building. PV arrays need to be connected to ElectricLoadCenter:Distribution objects that have a DC buss type.
The Generator:PV:Simple object describes about the simplest model for predicting photovoltaic energy production. In this model the user specifies the efficiency with which surfaces convert incident solar radiation to electricity. (In the other models this efficiency is determined as part of the model.) The full geometric model for solar radiation is used, including sky models, shading, and reflections, to determine the incident solar resource. The model accepts arbitrary conversion efficiencies and does not require actual production units be tested to obtain empirical performance coefficients. (The Energy+.idd sets the range of conversion efficiencies to be on [0..1], but the user could alter the Energy+.idd to extend this range if desired.)
Mathematical Description[LINK]
Table 89. Nomenclature for Simple Photovoltaic model
Mathematical variable

Description

P

Electrical power produced by photovoltaics [W]

A_{surf}

Net area of surface [m^{2}]

f_{activ}

Fraction of surface area with active solar cells

G_{T}

Total solar radiation incident on PV array [W/m^{2}]

h_{cell}

Module conversion efficiency

h_{invert}

DC to AC conversion efficiency

The usable electrical power produced by a PV surface are calculated using:
P=Asurf⋅factiv⋅GT⋅ηcell⋅ηinvert
On the right hand side of this equation, only G_{T} is calculated by EnergyPlus and the rest are user inputs. Power levels are assumed constant over the timestep to arrive at energy production.
There are two modes that can be selected by the user that govern how the PV system is coupled to the building surfaces. If the integration mode is selected as ‘DECOUPLED’ then no adjustments are made to account for energy extracted in the form of electricity. If the integration mode is selected as ‘INTEGRATED’ then the energy extracted in the form of electricity is removed from surface heat transfer calculations using a sink term. This sink term is lagged from the previous timestep.
Equivalent OneDiode Model[LINK]
This model predicts the electrical performance of a photovoltaic (PV) array. This model is also known as the “TRNSYS PV” model.
Mathematically speaking, the EnergyPlus PV module employs equations for an empirical equivalent circuit model to predict the currentvoltage characteristics of a single module. This circuit consists of a DC current source, diode, and either one or two resistors. The strength of the current source is dependent on solar radiation and the IV characteristics of the diode are temperaturedependent. The results for a single module equivalent circuit are extrapolated to predict the performance of a multimodule array.
The module employs a “fourparameter” equivalent circuit to model crystalline (both mono and poly) PV modules developed at the University of Wisconsin – Madison [2]. The values of these parameters cannot normally be obtained directly from manufacturers’ catalogs. However, the PV module will automatically calculate them from commonly available data. The PV module also includes an optional incidence angle modifier correlation to calculate how the reflectance of the PV module surface varies with the angle of incidence of solar radiation.
The module determines PV current as a function of load voltage. Other OUTPUTS include current and voltage at the maximum power point along the IV curve, opencircuit voltage, short circuit current as well as electrical load met and unmet.
Table 90. General Nomenclature for the PV model
Mathematical variable

Description

b

Slope of PV array [degrees]

g

Empirical PV curvefitting parameter

e_{g}

Semiconductor bandgap [eV]

h_{c}

Module conversion efficiency

m_{Isc}

Temperature coefficient of shortcircuit current [A/K]

m_{Voc}

Temperature coefficient of opencircuit voltage [V/K]

q

Angle of incidence for solar radiation [degrees]

ta

Module transmittanceabsorptance product

ta_{normal}

Module transmittanceabsorptance product at normal incidence

G_{T}

Total radiation incident on PV array

G_{T,beam}

Beam component of incident radiation

G_{T,diff}

Diffuse component of incident radiation

G_{T,gnd}

Groundreflected component of incident radiation

G_{T,NOCT}

Incident radiation at NOCT conditions

G_{T,ref}

Incident radiation at reference conditions

I

Current

I_{L}

Module photocurrent

I_{L,ref}

Module photocurrent at reference conditions

I_{o}

Diode reverse saturation current

I_{o,ref}

Diode reverse saturation current at reference conditions

I_{sc}

Shortcircuit current

I_{sc,ref}

Shortcircuit current at reference conditions

I_{mp}

Current at maximum power point along IV curve

I_{mp,ref}

Current at maximum power point along IV curve, reference conditions

IAM

Dimensionless incidence angle modifier

K

Boltzmann constant [J/K]

NP

Number of modules in parallel in array

NS

Number of modules in series in array

N_{s}

Number of individual cells in module

P

PV output power

P_{max}

PV output power at maximum power point along IV curve

Q

Electron charge constant

R_{s}

Module series resistance [W]

R_{sh}

Module shunt resistance [W]

T_{c}

Module temperature [K]

T_{c,NOCT}

Module temperature at NOCT conditions [K]

T_{c,ref}

Module temperature at reference conditions [K]

U_{L}

Array thermal loss coefficient

V

Voltage

V_{mp}

Voltage at maximum power point along IV curve

V_{mp,ref}

Voltage at maximum power point along IV curve, reference conditions

V_{oc}

Opencircuit voltage

V_{oc,ref}

Opencircuit voltage at reference conditions [V]

Mathematical Description[LINK]
PV Section 1: FourParameter Model[LINK]
The fourparameter equivalent circuit model was developed largely by Townsend [1989] and is detailed by Duffie and Beckman [1991]. The model was first incorporated into a component for the TRNSYS simulation package by Eckstein [1990]. The EnergyPlus module employs the Eckstein model for crystalline PV modules, using it whenever the shortcircuit IV slope is set to zero or a positive value as modified by Ulleberg [2000]. The four parameter model assumes that the slope of the IV curve is zero at the shortcircuit condition:
(dIdV)v=0=0
This is a reasonable approximation for crystalline modules. The “four parameters” in the model are I_{L,ref}, I_{o,ref,} g, and R_{s}. These are empirical values that cannot be determined directly through physical measurement. The EnergyPlus model calculates these values from manufactures’ catalog data as discussed in the following section on calculating these parameters
The fourparameter equivalent circuit is shown in the following figure:
Figure 339. Equivalent circuit in the four parameter model
V is the load voltage and I is the current flowing through the load and PV.
Determining Performance under Operating Conditions
The IV characteristics of a PV change with both insolation and temperature. The PV model employs these environmental conditions along with the four module constants I_{L,ref}, I_{o,ref,} g, and R_{s} to generate an IV curve at each timestep.
The currentvoltage equation of circuit shown in the previous figure is as follows:
I=IL−Io[exp(qγkTc(V+IRs))−1]
R_{s} and g are constants. The photocurrent I_{L} depends linearly on incident radiation:
IL=IL,refGTGT,ref
The reference insolation G_{ref} is nearly always defined as 1000 W/m^{2}. The diode reverse saturation current I_{o} is a temperature dependent quantity:
IoIo,ref=(TcTc,ref)3
Equation gives the current implicitly as a function of voltage. Once I_{o} and I_{L} are found from Eqs. 3 and 4, Newton’s method is employed to calculate the PV current. In addition, an iterative search routine finds the current (I_{mp}) and voltage (V_{mp}) at the point of maximum power along the IV curve.
Calculating I_{L,ref}, I_{o,ref,} g, and R_{s}
The Idf specification for the PV model include several values which must be read from manufacturers’ PV module catalogs. The manufactures’ values are used to determine the equivalent circuit characteristics I_{L,ref}, I_{o,ref,} g, and R_{s}. These characteristics define an equivalent circuit that is employed to find the PV performance at each timestep, as described previously. This section describes the algebra and calculation algorithms used to solve for the four equivalent circuit characteristics.
Three of these values, I_{L,ref}, I_{o,ref,} g, may be isolated algebraically. The first step is to substitute the current and voltage into Eq. at the opencircuit, short circuit, and maximum power conditions:
0=IL,ref−Io,ref[exp(qγkTc,refVoc,ref)−1]−Voc,refRsh
Isc,ref=IL,ref−Io,ref[exp(qIsc,refRsγkTc,ref)−1]−Isc,refRsRsh
Imp,ref=IL,ref−Io,ref[exp(qγkTc,ref(Vmp,ref+Imp,refRs))−1]−Vmp,ref+Imp,refRsRsh
In each case the “1” term is may be dropped to simplify the algebra. This approximation has little influence on the right side of the equations since because the magnitude of I_{o} is very small, generally on the order of 10^{6} A. Some rearrangement then yields the following three expressions which isolate I_{L,ref}, I_{o,ref,} g:
IL,ref≈Isc,ref
γ=q(Vmp,ref−Voc,ref+Imp,refRs)kTc,refln(1−Imp,refIsc,ref)
Io,ref=Isc,refexp(qVoc,refγkTc,ref)
At this point an additional equation is needed in order to determine the last unknown parameter. Taking the analytical derivative of voltage with respect to temperature at the reference opencircuit condition derives the fourth equation. This analytical value is matched to the opencircuit temperature coefficient, a catalog specification:
∂Voc∂Tc=μvoc=γkq[ln(Isc,refIo,ref)+TcμiscIsc,ref−(3+qεAkTc,ref)]
where
A=γNs
The “TRNSYS PV model” uses an iterative search routine in these four equations to calculate the equivalent circuit characteristics. The first step is to set upper and lower bounds for the series resistance parameter R_{s}: physical constraints require the R_{s} value to lie between 0 and the value such that g = N_{s}. The initial guess for R_{s} is midway between these bounds. g and I_{o,ref} are found from Eq. and Eq. , while Eq. gives a trivial solution for I_{L,ref}. The model then employs Eq. to compare the analytical and catalog values for m_{voc}. When all other variables are held constant, the analytical value for m_{voc} increases monotonically with series resistance (Townsend 1989). If the analytical voltage coefficient is less than the catalog value, the lower bound for R_{s} is reset to the present guess value. Likewise, the upper bound is set to the current value if the calculated m_{voc} is too large. After resetting the upper or lower bound for R_{s}, a new guess value is found by averaging the bounds. This procedure repeats until R_{s} and g converge. Note that for I_{L,ref}, I_{o,ref,} g, and R_{s} are assumed to be constant and are calculated only on the first call in the simulation. Alternatively, the user may enter a known series resistance by entering a positive value in the IDF. In this case the iterative routine described above is skipped and Eqs. , , and find I_{L,ref}, I_{o,ref,} and g directly from the given value of R_{s}.
PV Section 2 : Module Operating Temperature[LINK]
The PV model uses one of five methods for determining cell temperature data. The cell temperature of a PV module is important because the hotter the temperature of the panel, the lower its electrical output. The cell temperature calculation method is chosen by the user in the EnergyPlus IDF file through a parameter choice in the IDD entry called Integration and Cell Temperature Mode.
If the value of this parameter is “Decoupled NOCT Conditions” then the cell temperature of the PV is modeled using the method from the Duffie and Beckman (1991) for estimating cell temperature. This is based upon the standard NOCT (Nominal Operating Cell Temperature) measurements to compute the module temperature Tc at each timestep. The NOCT temperature (Tc,NOCT) is the operating temperature of the module with a wind speed of 1 m/s, no electrical load, and a certain specified insolation and ambient temperature [Beckman and Duffie, 1991]. The values for insolation GT,NOCT and ambient temperature T_{a,NOCT} are usually 800 W/m^{2} and 20º C. h_{c} is the convesion efficiency of the module, which varies with ambient conditions. ta is a userdefined constant.
The equation is:
Tc=Ta+(1−{{\eta_c}}/ηcτα{\tau \alpha })({{G_T}\tau \alpha }/GTταUL{{U_L}})
If the user specifies the “Decoupled Ulleberg Dynamic” mode for calculating cell temperature, then a method developed by Ulleberg is used:
Tcellt=Tambient+(Tcellt−1−Tambient)∗e−ULCapΔt
In other words, the cell temperature is a function of the privious cell temperature and the thermal capacity of the PV module material.
If the user specifies “Integrated Surface Outside Face” for this parameter, then the temperature result from EnergyPlus’s modeling of surfaces is used for the cell temperature. Also the energy exported from the surface as electricity becomes a sink in the internal source modeling for the heat transfer surface.
If the user specifies “Integrated Transpired Collector” for this parameter, then the temperature result for the unglazed transpired collector surfaces is used for the cell temperature. Also the energy exported from the collector surface as electricity is deprecated using a source term in the collector’s temperature modeling.
If the user specifies “Integrated Exterior Vented Cavity” for this parameter, then the temperature result for the exterior cavity is used for the cell temperature. Also the energy exported from the baffle surface as electricity is deprecated using a source term in the baffle’s temperature modeling.
PV Section 3 : MultiArray Modules[LINK]
The electrical calculations discussed in the sections above deal only with a single module. The EnergyPlus PV component may be used to simulate arrays with any number of modules. The IDF defines the number of modules in series (NS) and modules in parallel (NP) for the entire array. The total number of modules in the array is the product of NS and NP. When simulating a single module only, both NS and NP are set to 1. The singlemodule values for all currents and voltages discussed in PV Section 1 are multiplied by NP or NS to find values for the entire array. This approach neglects module mismatch losses.
With the above equations, and the assumption that the panels operate at the maximum power point, it is a direct calculation to determine DC power production. The performance of an array of identical modules is assumed to be linear with the number of modules in series and parallel. The inverter efficiency is applied linearly to derate the energy production. The inverter capacity forms a limit for power production from a PV generator. A ‘load’ is passed the PV array acting as a generator and various trivial calculations compare PV production to this load. If the PV array is associated with a surface that is associated with a zone, then if the zone has any multipliers associated with it, electricity production will be multiplied accordingly.
Duffie, John A. and William A. Beckman. 1991. Solar Engineering of Thermal Processes. New York: John Wiley & Sons, Inc..
Eckstein, Jürgen Helmut. 1990. Detailed Modeling of Photovoltaic Components. M. S. Thesis – Solar Energy Laboratory, University of Wisconsin, Madison: 1990.
Ulleberg, Øystein. HYDROGEMS Component Library for TRNSYS 15 User Manual, Institute for Energy Technology, Kjeller, Norway
The third model available in EnergyPlus for predicting the electricity generated by photovoltaics is referred to as the Sandia model. This model is based on work done at Sandia National Lab, Albuquerque, NM by David King – with the help of many others. The model consists of a series of empirical relationships with coefficients that are derived from actual testing. Once the coefficients for a particular module are available, it is straightforward matter to use the model equations to calculate five select points on the currentvoltage curve.
The implementation in EnergyPlus is also based on work done by Greg Barker (2003) for the National Renewable Energy Lab who implemented the Sandia model in FORTRAN77 as a custom type (Type101) for the TRNSYS computer program.
There are several climate and solar orientation inputs to the model that are managed elsewhere in EnergyPlus including: incident beam solar, incident diffuse solar, incidence angle of beam solar, solar zenith Angle, outdoor drybulb, wind speed, and elevation.
Mathematical Description[LINK]
This section presents the mathematical description of the Sandia model from a draft report by King et, al. (2003). The core of the model predicts the performance of a single PV module. The following nomenclature and equations summarize the Sandia model.
Table 91. Nomenclature for Sandia PV model
Mathematical variable

Description

I_{sc}

Shortcircuit current (A)

I_{mp}

Current at the maximumpower point (A)

I_{x}

Current at module V = 0.5 V_{oc}, defines 4th point on IV curve

I_{xx}

Current at module V = 0.5 (V_{oc} + V_{mp}), defines a 5th point on the IV curve

V_{oc}

Opencircuit voltage (V)

V_{mp}

Voltage at maximumpower point (V)

P_{mp}

Power at maximumpower point (W)

fd

Fraction of diffuse irradiance used by module

N_{s}

Number of cells in series in a module’s cellstring

N_{p}

Number of cellstrings in parallel in module

k

Boltzmann’s constant, 1.38066E23 (J/k)

q

Elementary charge, 1.60218E19 (coulomb)

T_{c}

Cell temperature inside module (°C)

δ(T_{c})

‘Thermal voltage’ per cell at temperature T_{c}, approximately 1 volt for a typical 26cell crystalline silicon module

E_{e}

‘Effective’ solar irradiance

E_{b}

Beam solar irradiance

E_{diff}

Diffuse solar irradiance

C_{0} , C_{1}

Empirical coefficients relating I_{mp} to E_{e} , C_{0} + C_{1} = 1 (both dimensionless)

C_{2} , C_{3}

Empirical coefficients relating V_{mp} to E_{e}(C_{2} dimensionless, C_{3} is 1/V)

C_{4} , C_{5}

Empirical coefficients relating I_{x} to E_{e}, C_{4} + C_{5} = 1 (both dimensionless)

C_{6} , C_{7}

Empirical coefficients relating I_{xx} to E_{e,} C_{6} + C_{7} = 1 (both dimensionless)

n

Empirically determined ‘diode factor’ for individual cells

AMa

Absolute Air Mas

AOI

Solar angleofincidence (degrees) from normal

f_{1}(AM_{a})

Empirical polynomial function used to relate shortcircuit current to the solar spectrum via air mass

f_{2}(AOI)

Empirical polynomial function used to relate shortcircuit current to the solar angleofincidence

a_{0}, a_{1}, a_{2}, a_{3}, a_{4}

Empirical coefficients for f_{1}(AM_{a}) polynomial

b_{0}, b_{1}, b_{2}, b_{3}, b_{4},b_{5},b_{6}

Empirical coefficients for f_{1}(AOI) polynomial

T_{o}

Reference cell temperature for rating, typically fixed at 25°C

I_{sco}

Short circuit current at reference conditions

I_{mpo}

Max power point current at reference conditions

V_{mpo}

Voltage at max power at reference conditions

V_{oco}

Open circuit voltage at reference conditions

I_{xo}

Current at V = 0.5 Voc and at reference conditions

I_{xxo}

Current at V = 0.5 (Vmp + Voc) and at reference conditions

α_{Isc}

Normalized temperature coefficient for I_{sc} (1/°C)

α_{Imp}

Normalized temperature coefficient for I_{mp} (1/°C)

β_{Voc}(E_{e})

Temperature coefficient for module opencircuitvoltage as function of E_{e}

β_{Voco}

Temperature coefficient for module opencircuitvoltage at reference conditions

m_{βVoco}

Coefficient for irradiance dependence of opencircuitvoltagetemperature coefficient, often zero (V/°C)

β_{Vmp}(E_{e})

Temperature coefficient for module maximumpowervoltage as a function of E_{e}

β_{Vmpo}

Temperature coefficient for module maximumpowervoltage at reference conditions

m_{βVmpo}

Cofficient for irradiance dependence of maximumpowervoltagetemperature coefficient, often zero (V/°C)

T_{m}

PV module temperature at back suface (°C)

T_{a}

Ambient outdoor drybulb temperature (°C)

E

Solar irradiance incident on module surface (W/m^{2})

WS

Wind speed at standard 10m height (m/s)

a

Empirical coefficient relating module temperature at low wind and high solar irradiance

b

Empirical coefficient relating module temperature decrease with increasing wind speed

T_{c}

Temperature of solar cell inside module (°C)

E_{o}

Reference solar irradiance (1000 W/m^{2})

ΔT

Temperature difference between T_{c} and T_{m}at E_{o}(°C), (This is d(Tc) in Sandia database)

The current implementation in EnergyPlus focuses on determining performance at the maximum powerpoint but also calculates, and reports, four other points on the IV curve so that the data are available for analyses outside of EnergyPlus. The equations below use the module performance parameters that are available in a database provided by Sandia National Laboratory (see www.sandia.gov/pv). The following equations form the basis of the Sandia model implemented in EnergyPlus:
Isc=Isco⋅f1(AMa)⋅{(Eb⋅f2(AOI)+fd⋅Ediff)/Eo}⋅{1+αIsc⋅(Tc−To)}
Imp=Impo⋅{Co⋅Ee+C1⋅E2e}⋅{1+αImp⋅(Tc−To}
Voc=Voco+Ns⋅δ(Tc)⋅ln(Ee)+βVoc(Ee)⋅(Tc−To)
Vmp=Vmpo+C2⋅Ns⋅δ(Tc)⋅ln(Ee)+C3⋅N3⋅{δ(Tc)⋅ln(Ee)}2+βVmp(Ee)⋅(Tc−To)
Pmp=Imp⋅Vmp
Ix=Ixo⋅{C4⋅Ee+C5⋅E2e}⋅{1+(αIsc)⋅(Tc−To)}
Ixx=Ixxo⋅{C6⋅Ee+C7⋅E2e}⋅{1+(αImp)⋅(Tc−To)}
where,
Ee=Isc/[Isco⋅{1+αIsc⋅(Tc−To)}]
δ(Tc)=n⋅k⋅(Tc+273.15)/q
f1(AMa)=a0+a1AMa+a2(AMa)2+a3(AMa)3+a4(AMa)4
f2(AOI)=b0+b1⋅AOI+b2(AOI)2+b3(AOI)3+b4(AOI)4+b5(AOI)5
βVoc(Ee)=βVoco+mβVoc⋅(1−Ee)
βVmp(Ee)=βVmpo+mβVmp⋅(1−Ee)
Tm=E⋅{ea+b⋅WS}+Ta
Tc=Tm+EEo⋅ΔT
With the above equations, and the assumption that the panels operate at the maximum power point, it is a direct calculation to determine DC power production. The performance of an array of identical modules is assumed to be linear with the number of modules in series and parallel. The inverter efficiency is applied linearly to derate the energy production. The inverter capacity forms a limit for power production from a PV generator. A ‘load’ is passed the PV array acting as a generator and various trivial calculations compare PV production to this load. If the PV array is associated with a surface that is associated with a zone, then if the zone has any multipliers associated with it, electricity production will be multiplied accordingly.
The equation above for T_{m} is used to predict backofmodule temperature when the mode ‘SANDIA RACK’ is selected. This would be appropriate for most rack mounted PV installations. If the user selects ‘EPLUS INTEGRATED’ then the backofmodule temperature is obtained from the outside face surface temperature calculated by the full complement of Heat Balance models using Conduction Transfer Functions that is native to EnergyPlus. And energy exported from the surface is accounted for using a source/sink term that is lagged from the previous timestep (pingpong).
King, D.L., Boyson, W.E., Kratochvil J.A. 2003. Photovoltaic Array Performance Model. Sandia National Laboratories, Albuquerque, NM 87185, November 2003 currently in DRAFT
Barker, G. 2003 (unknown). Predicting LongTerm Performance of Photovoltaic Arrays. Submitted as deliverable for NREL subcontract LAX13048002. Currently in DRAFT
Davis, M.W., Fanney, A.H., and Dougherty B.P. 2002. Measured Versus Predicted Performance of Building Integrated Photovoltaics. from the conference Solar 2002, Sunrise on the Reliable Energy Economy, June 1519, 2002, Reno, NV. Available from NIST website.
King, D.L. 1996. Photovoltaic Module and Array Performance Characterization Methods for All System Operating Conditions. Sandia National Laboratory. Albuquerque, NM 87185
Photovoltaic Arrays[LINK]
The Photovoltaics.f90 module includes three different models referred to as “Simple”, “Equivalent OneDiode” and “Sandia” and the choice will determine the mathematical models (and input data) used to determine the energy produced by solar/electric conversion panels. The EnergyPlus photovoltaic array models are called one at a time at the HVAC system timestep along with other electrical generation components such as gas turbines and diesel engines.
All of the photovoltaic models share the same models for predicting incident solar radiation that are also used for the solar thermal calculations and are described in the section Climate, Sky and Solar/Shading Calculations.
Note that some of the terminology used to discussed photovoltaics overlaps with terminology used to discuss Fortran programs. The word module may refer to a PV panel or to a fortran90 programming entity. Model may refer to a manufacturers production model for a specific type of PV module or to a mathematical model used for engineering analysis. Array may refer to a collection of PV modules wired together or to a mathematical variable with multiple elements.
The PV modules are assumed to always run when the total incident solar is greater than 0.3 Watts. If the incident solar is less than 0.3, then the modules produce no power.
PV arrays are managed by an electric load center. The load center is a “load” with respect to generating equipment but is a “supply center” for the rest of the building. PV arrays need to be connected to ElectricLoadCenter:Distribution objects that have a DC buss type.
Simple Model[LINK]
The Generator:PV:Simple object describes about the simplest model for predicting photovoltaic energy production. In this model the user specifies the efficiency with which surfaces convert incident solar radiation to electricity. (In the other models this efficiency is determined as part of the model.) The full geometric model for solar radiation is used, including sky models, shading, and reflections, to determine the incident solar resource. The model accepts arbitrary conversion efficiencies and does not require actual production units be tested to obtain empirical performance coefficients. (The Energy+.idd sets the range of conversion efficiencies to be on [0..1], but the user could alter the Energy+.idd to extend this range if desired.)
Mathematical Description[LINK]
Table 89. Nomenclature for Simple Photovoltaic modelThe usable electrical power produced by a PV surface are calculated using:
P=Asurf⋅factiv⋅GT⋅ηcell⋅ηinvert
On the right hand side of this equation, only G_{T} is calculated by EnergyPlus and the rest are user inputs. Power levels are assumed constant over the timestep to arrive at energy production.
There are two modes that can be selected by the user that govern how the PV system is coupled to the building surfaces. If the integration mode is selected as ‘DECOUPLED’ then no adjustments are made to account for energy extracted in the form of electricity. If the integration mode is selected as ‘INTEGRATED’ then the energy extracted in the form of electricity is removed from surface heat transfer calculations using a sink term. This sink term is lagged from the previous timestep.
Equivalent OneDiode Model[LINK]
This model predicts the electrical performance of a photovoltaic (PV) array. This model is also known as the “TRNSYS PV” model.
Mathematically speaking, the EnergyPlus PV module employs equations for an empirical equivalent circuit model to predict the currentvoltage characteristics of a single module. This circuit consists of a DC current source, diode, and either one or two resistors. The strength of the current source is dependent on solar radiation and the IV characteristics of the diode are temperaturedependent. The results for a single module equivalent circuit are extrapolated to predict the performance of a multimodule array.
The module employs a “fourparameter” equivalent circuit to model crystalline (both mono and poly) PV modules developed at the University of Wisconsin – Madison [2]. The values of these parameters cannot normally be obtained directly from manufacturers’ catalogs. However, the PV module will automatically calculate them from commonly available data. The PV module also includes an optional incidence angle modifier correlation to calculate how the reflectance of the PV module surface varies with the angle of incidence of solar radiation.
The module determines PV current as a function of load voltage. Other OUTPUTS include current and voltage at the maximum power point along the IV curve, opencircuit voltage, short circuit current as well as electrical load met and unmet.
Table 90. General Nomenclature for the PV modelMathematical Description[LINK]
PV Section 1: FourParameter Model[LINK]
The fourparameter equivalent circuit model was developed largely by Townsend [1989] and is detailed by Duffie and Beckman [1991]. The model was first incorporated into a component for the TRNSYS simulation package by Eckstein [1990]. The EnergyPlus module employs the Eckstein model for crystalline PV modules, using it whenever the shortcircuit IV slope is set to zero or a positive value as modified by Ulleberg [2000]. The four parameter model assumes that the slope of the IV curve is zero at the shortcircuit condition:
(dIdV)v=0=0
This is a reasonable approximation for crystalline modules. The “four parameters” in the model are I_{L,ref}, I_{o,ref,} g, and R_{s}. These are empirical values that cannot be determined directly through physical measurement. The EnergyPlus model calculates these values from manufactures’ catalog data as discussed in the following section on calculating these parameters
The fourparameter equivalent circuit is shown in the following figure:
Figure 339. Equivalent circuit in the four parameter model
V is the load voltage and I is the current flowing through the load and PV.
Determining Performance under Operating Conditions
The IV characteristics of a PV change with both insolation and temperature. The PV model employs these environmental conditions along with the four module constants I_{L,ref}, I_{o,ref,} g, and R_{s} to generate an IV curve at each timestep.
The currentvoltage equation of circuit shown in the previous figure is as follows:
I=IL−Io[exp(qγkTc(V+IRs))−1]
R_{s} and g are constants. The photocurrent I_{L} depends linearly on incident radiation:
IL=IL,refGTGT,ref
The reference insolation G_{ref} is nearly always defined as 1000 W/m^{2}. The diode reverse saturation current I_{o} is a temperature dependent quantity:
IoIo,ref=(TcTc,ref)3
Equation gives the current implicitly as a function of voltage. Once I_{o} and I_{L} are found from Eqs. 3 and 4, Newton’s method is employed to calculate the PV current. In addition, an iterative search routine finds the current (I_{mp}) and voltage (V_{mp}) at the point of maximum power along the IV curve.
Calculating I_{L,ref}, I_{o,ref,} g, and R_{s}
The Idf specification for the PV model include several values which must be read from manufacturers’ PV module catalogs. The manufactures’ values are used to determine the equivalent circuit characteristics I_{L,ref}, I_{o,ref,} g, and R_{s}. These characteristics define an equivalent circuit that is employed to find the PV performance at each timestep, as described previously. This section describes the algebra and calculation algorithms used to solve for the four equivalent circuit characteristics.
Three of these values, I_{L,ref}, I_{o,ref,} g, may be isolated algebraically. The first step is to substitute the current and voltage into Eq. at the opencircuit, short circuit, and maximum power conditions:
0=IL,ref−Io,ref[exp(qγkTc,refVoc,ref)−1]−Voc,refRsh
Isc,ref=IL,ref−Io,ref[exp(qIsc,refRsγkTc,ref)−1]−Isc,refRsRsh
Imp,ref=IL,ref−Io,ref[exp(qγkTc,ref(Vmp,ref+Imp,refRs))−1]−Vmp,ref+Imp,refRsRsh
In each case the “1” term is may be dropped to simplify the algebra. This approximation has little influence on the right side of the equations since because the magnitude of I_{o} is very small, generally on the order of 10^{6} A. Some rearrangement then yields the following three expressions which isolate I_{L,ref}, I_{o,ref,} g:
IL,ref≈Isc,ref
γ=q(Vmp,ref−Voc,ref+Imp,refRs)kTc,refln(1−Imp,refIsc,ref)
Io,ref=Isc,refexp(qVoc,refγkTc,ref)
At this point an additional equation is needed in order to determine the last unknown parameter. Taking the analytical derivative of voltage with respect to temperature at the reference opencircuit condition derives the fourth equation. This analytical value is matched to the opencircuit temperature coefficient, a catalog specification:
∂Voc∂Tc=μvoc=γkq[ln(Isc,refIo,ref)+TcμiscIsc,ref−(3+qεAkTc,ref)]
where
A=γNs
The “TRNSYS PV model” uses an iterative search routine in these four equations to calculate the equivalent circuit characteristics. The first step is to set upper and lower bounds for the series resistance parameter R_{s}: physical constraints require the R_{s} value to lie between 0 and the value such that g = N_{s}. The initial guess for R_{s} is midway between these bounds. g and I_{o,ref} are found from Eq. and Eq. , while Eq. gives a trivial solution for I_{L,ref}. The model then employs Eq. to compare the analytical and catalog values for m_{voc}. When all other variables are held constant, the analytical value for m_{voc} increases monotonically with series resistance (Townsend 1989). If the analytical voltage coefficient is less than the catalog value, the lower bound for R_{s} is reset to the present guess value. Likewise, the upper bound is set to the current value if the calculated m_{voc} is too large. After resetting the upper or lower bound for R_{s}, a new guess value is found by averaging the bounds. This procedure repeats until R_{s} and g converge. Note that for I_{L,ref}, I_{o,ref,} g, and R_{s} are assumed to be constant and are calculated only on the first call in the simulation. Alternatively, the user may enter a known series resistance by entering a positive value in the IDF. In this case the iterative routine described above is skipped and Eqs. , , and find I_{L,ref}, I_{o,ref,} and g directly from the given value of R_{s}.
PV Section 2 : Module Operating Temperature[LINK]
The PV model uses one of five methods for determining cell temperature data. The cell temperature of a PV module is important because the hotter the temperature of the panel, the lower its electrical output. The cell temperature calculation method is chosen by the user in the EnergyPlus IDF file through a parameter choice in the IDD entry called Integration and Cell Temperature Mode.
If the value of this parameter is “Decoupled NOCT Conditions” then the cell temperature of the PV is modeled using the method from the Duffie and Beckman (1991) for estimating cell temperature. This is based upon the standard NOCT (Nominal Operating Cell Temperature) measurements to compute the module temperature Tc at each timestep. The NOCT temperature (Tc,NOCT) is the operating temperature of the module with a wind speed of 1 m/s, no electrical load, and a certain specified insolation and ambient temperature [Beckman and Duffie, 1991]. The values for insolation GT,NOCT and ambient temperature T_{a,NOCT} are usually 800 W/m^{2} and 20º C. h_{c} is the convesion efficiency of the module, which varies with ambient conditions. ta is a userdefined constant.
The equation is:
Tc=Ta+(1−{{\eta_c}}/ηcτα{\tau \alpha })({{G_T}\tau \alpha }/GTταUL{{U_L}})
If the user specifies the “Decoupled Ulleberg Dynamic” mode for calculating cell temperature, then a method developed by Ulleberg is used:
Tcellt=Tambient+(Tcellt−1−Tambient)∗e−ULCapΔt
In other words, the cell temperature is a function of the privious cell temperature and the thermal capacity of the PV module material.
If the user specifies “Integrated Surface Outside Face” for this parameter, then the temperature result from EnergyPlus’s modeling of surfaces is used for the cell temperature. Also the energy exported from the surface as electricity becomes a sink in the internal source modeling for the heat transfer surface.
If the user specifies “Integrated Transpired Collector” for this parameter, then the temperature result for the unglazed transpired collector surfaces is used for the cell temperature. Also the energy exported from the collector surface as electricity is deprecated using a source term in the collector’s temperature modeling.
If the user specifies “Integrated Exterior Vented Cavity” for this parameter, then the temperature result for the exterior cavity is used for the cell temperature. Also the energy exported from the baffle surface as electricity is deprecated using a source term in the baffle’s temperature modeling.
PV Section 3 : MultiArray Modules[LINK]
The electrical calculations discussed in the sections above deal only with a single module. The EnergyPlus PV component may be used to simulate arrays with any number of modules. The IDF defines the number of modules in series (NS) and modules in parallel (NP) for the entire array. The total number of modules in the array is the product of NS and NP. When simulating a single module only, both NS and NP are set to 1. The singlemodule values for all currents and voltages discussed in PV Section 1 are multiplied by NP or NS to find values for the entire array. This approach neglects module mismatch losses.
With the above equations, and the assumption that the panels operate at the maximum power point, it is a direct calculation to determine DC power production. The performance of an array of identical modules is assumed to be linear with the number of modules in series and parallel. The inverter efficiency is applied linearly to derate the energy production. The inverter capacity forms a limit for power production from a PV generator. A ‘load’ is passed the PV array acting as a generator and various trivial calculations compare PV production to this load. If the PV array is associated with a surface that is associated with a zone, then if the zone has any multipliers associated with it, electricity production will be multiplied accordingly.
References[LINK]
Duffie, John A. and William A. Beckman. 1991. Solar Engineering of Thermal Processes. New York: John Wiley & Sons, Inc..
Eckstein, Jürgen Helmut. 1990. Detailed Modeling of Photovoltaic Components. M. S. Thesis – Solar Energy Laboratory, University of Wisconsin, Madison: 1990.
Ulleberg, Øystein. HYDROGEMS Component Library for TRNSYS 15 User Manual, Institute for Energy Technology, Kjeller, Norway
Sandia Photovoltaic Performance Model[LINK]
The third model available in EnergyPlus for predicting the electricity generated by photovoltaics is referred to as the Sandia model. This model is based on work done at Sandia National Lab, Albuquerque, NM by David King – with the help of many others. The model consists of a series of empirical relationships with coefficients that are derived from actual testing. Once the coefficients for a particular module are available, it is straightforward matter to use the model equations to calculate five select points on the currentvoltage curve.
The implementation in EnergyPlus is also based on work done by Greg Barker (2003) for the National Renewable Energy Lab who implemented the Sandia model in FORTRAN77 as a custom type (Type101) for the TRNSYS computer program.
There are several climate and solar orientation inputs to the model that are managed elsewhere in EnergyPlus including: incident beam solar, incident diffuse solar, incidence angle of beam solar, solar zenith Angle, outdoor drybulb, wind speed, and elevation.
Mathematical Description[LINK]
This section presents the mathematical description of the Sandia model from a draft report by King et, al. (2003). The core of the model predicts the performance of a single PV module. The following nomenclature and equations summarize the Sandia model.
Table 91. Nomenclature for Sandia PV modelThe current implementation in EnergyPlus focuses on determining performance at the maximum powerpoint but also calculates, and reports, four other points on the IV curve so that the data are available for analyses outside of EnergyPlus. The equations below use the module performance parameters that are available in a database provided by Sandia National Laboratory (see www.sandia.gov/pv). The following equations form the basis of the Sandia model implemented in EnergyPlus:
Isc=Isco⋅f1(AMa)⋅{(Eb⋅f2(AOI)+fd⋅Ediff)/Eo}⋅{1+αIsc⋅(Tc−To)}
Imp=Impo⋅{Co⋅Ee+C1⋅E2e}⋅{1+αImp⋅(Tc−To}
Voc=Voco+Ns⋅δ(Tc)⋅ln(Ee)+βVoc(Ee)⋅(Tc−To)
Vmp=Vmpo+C2⋅Ns⋅δ(Tc)⋅ln(Ee)+C3⋅N3⋅{δ(Tc)⋅ln(Ee)}2+βVmp(Ee)⋅(Tc−To)
Pmp=Imp⋅Vmp
Ix=Ixo⋅{C4⋅Ee+C5⋅E2e}⋅{1+(αIsc)⋅(Tc−To)}
Ixx=Ixxo⋅{C6⋅Ee+C7⋅E2e}⋅{1+(αImp)⋅(Tc−To)}
where,
Ee=Isc/[Isco⋅{1+αIsc⋅(Tc−To)}]
δ(Tc)=n⋅k⋅(Tc+273.15)/q
f1(AMa)=a0+a1AMa+a2(AMa)2+a3(AMa)3+a4(AMa)4
f2(AOI)=b0+b1⋅AOI+b2(AOI)2+b3(AOI)3+b4(AOI)4+b5(AOI)5
βVoc(Ee)=βVoco+mβVoc⋅(1−Ee)
βVmp(Ee)=βVmpo+mβVmp⋅(1−Ee)
Tm=E⋅{ea+b⋅WS}+Ta
Tc=Tm+EEo⋅ΔT
With the above equations, and the assumption that the panels operate at the maximum power point, it is a direct calculation to determine DC power production. The performance of an array of identical modules is assumed to be linear with the number of modules in series and parallel. The inverter efficiency is applied linearly to derate the energy production. The inverter capacity forms a limit for power production from a PV generator. A ‘load’ is passed the PV array acting as a generator and various trivial calculations compare PV production to this load. If the PV array is associated with a surface that is associated with a zone, then if the zone has any multipliers associated with it, electricity production will be multiplied accordingly.
The equation above for T_{m} is used to predict backofmodule temperature when the mode ‘SANDIA RACK’ is selected. This would be appropriate for most rack mounted PV installations. If the user selects ‘EPLUS INTEGRATED’ then the backofmodule temperature is obtained from the outside face surface temperature calculated by the full complement of Heat Balance models using Conduction Transfer Functions that is native to EnergyPlus. And energy exported from the surface is accounted for using a source/sink term that is lagged from the previous timestep (pingpong).
References[LINK]
King, D.L., Boyson, W.E., Kratochvil J.A. 2003. Photovoltaic Array Performance Model. Sandia National Laboratories, Albuquerque, NM 87185, November 2003 currently in DRAFT
Barker, G. 2003 (unknown). Predicting LongTerm Performance of Photovoltaic Arrays. Submitted as deliverable for NREL subcontract LAX13048002. Currently in DRAFT
Davis, M.W., Fanney, A.H., and Dougherty B.P. 2002. Measured Versus Predicted Performance of Building Integrated Photovoltaics. from the conference Solar 2002, Sunrise on the Reliable Energy Economy, June 1519, 2002, Reno, NV. Available from NIST website.
King, D.L. 1996. Photovoltaic Module and Array Performance Characterization Methods for All System Operating Conditions. Sandia National Laboratory. Albuquerque, NM 87185
Documentation content copyright © 19962015 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.