# Solar Collectors[LINK]

Solar collectors are devices that convert solar energy into thermal energy by raising the temperature of a circulating heat transfer fluid. The fluid can then be used to heat water for domestic hot water usage or space heating. Flat-plate solar collectors using water as the heat transfer fluid, Integral-Collector Storage solar collectors using water and unglazed transpired solar collectors using air are currently the only types of collector available in EnergyPlus.

## Flat-Plate Solar Collectors[LINK]

The input object SolarCollector:FlatPlate:Water provides a model for flat-plate solar collectors that are the most common type of collector. Standards have been established by ASHRAE for the performance testing of these collectors (ASHRAE 1989; 1991) and the Solar Rating and Certification Corporation (SRCC) publishes a directory of commercially available collectors in North America (SRCC 2003).

The EnergyPlus model is based on the equations found in the ASHRAE standards and Duffie and Beckman (1991). This model applies to glazed and unglazed flat-plate collectors, as well as banks of tubular, i.e. evacuated tube, collectors.

### Solar and Shading Calculations[LINK]

The solar collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account. Likewise, the collector surface can shade other surfaces, for example, reducing the incident radiation on the roof beneath it.

The thermal efficiency of a collector is defined as the ratio of the useful heat gain of the collector fluid versus the total incident solar radiation on the gross surface area of the collector.

- η=(q/A)Isolar

where

q = useful heat gain

A = gross area of the collector

I_{solar} = total incident solar radiation

Notice that the efficiency is only defined for I_{solar} > 0.

An energy balance on a solar collector with double glazing shows relationships between the glazing properties, absorber plate properties, and environmental conditions.

- qA=Isolarτg1τg2αabs−T4abs−T4g2Rrad−Tabs−Tg2Rconv−Tabs−TairRcond

where

_{g1} = transmittance of the first glazing layer

_{g2} = transmittance of the second glazing layer

_{abs} = absorptance of the absorber plate

R_{rad} = radiative resistance from absorber to inside glazing

R_{conv} = convective resistance from absorber to inside glazing

R_{cond} = conductive resistance from absorber to outdoor air through the insulation

T_{abs} = temperature of the absorber plate

T_{g2} = temperature of the inside glazing

T_{air} = temperature of the outdoor air

The equation above can be approximated with a simpler formulation as:

- qA=FR[Isolar(τα)−UL(Tin−Tair)]

where

F_{R} = an empirically determined correction factor

(τα) = the product of all transmittance and absorptance terms

U_{L} = overall heat loss coefficient combining radiation, convection, and conduction terms

T_{in} = inlet temperature of the working fluid

Substituting this into Equation (1),

- η=FR(τα)−FRUL(Tin−Tair)Isolar

A linear correlation can be constructed by treating F_{R}(τα) and F_{R}U_{L} as characteristic constants of the solar collector:

- η=c0+c1(Tin−Tair)Isolar

Similarly, a quadratic correlation can be constructed using the form:

- η=c0+c1(Tin−Tair)Isolar+c2(Tin−Tair)2Isolar

Both first- and second-order efficiency equation coefficients are listed in the *Directory of SRCC Certified Solar Collector Ratings*.

### Incident Angle Modifiers[LINK]

As with regular windows the transmittance of the collector glazing varies with the incidence angle of radiation. Usually the transmittance is highest when the incident radiation is normal to the glazing surface. Test conditions determine the efficiency coefficients for normal incidence. For off-normal angles, the transmittance of the glazing is modified by an *incident angle modifier* coefficient.

- Kτα=(τα)(τα)n

Additional testing determines the incident angle modifier as a function of incident angle. This relationship can be fit to a first-order, linear correlation:

- Kτα=1+b0(1cosθ−1)

or a second-order, quadratic correlation:

- Kτα=1+b0(1cosθ−1)+b1(1cosθ−1)2

The incident angle modifier coefficients *b*_{0} and *b*_{1} are usually negative, although some collectors have a positive value for *b*_{0}. Both first- and second-order incident angle modifier equation coefficients are listed in the *Directory of SRCC Certified Solar Collector Ratings*.

The SRCC incident angle modifier equation coefficients are only valid for incident angles of 60 degrees or less. Because these curves can be valid yet behave poorly for angles greater than 60 degree, the EnergyPlus model cuts off collector gains for incident angles greater than 60 degrees.

For flat-plate collectors, the incident angle modifier is generally symmetrical. However, for tubular collectors the incident angle modifier is different depending on whether the incident angle is parallel or perpendicular to the tubes. These are called bi-axial modifiers. Some special flat-plate collectors may also exhibit this asymmetry. The current model cannot yet handle two sets of incident angle modifiers. In the meantime it is recommended that tubular collectors be approximated with caution using either the parallel or perpendicular correlation.

Incident angle modifiers are calculated separately for sun, sky, and ground radiation. The net incident angle modifier for all incident radiation is calculated by weighting each component by the corresponding modifier.

- Kτα,net=IbeamKτα,beam+IskyKτα,sky+IgndKτα,gndIbeam+Isky+Ignd

For sky and ground radiation the incident angle is approximated using Brandemuehl and Beckman’s equations:

θsky=59.68−0.1388ϕ+0.001497ϕ2

θground=90.0−0.5788ϕ+0.002693ϕ2

where ϕ is the surface tilt in degrees.

The net incident angle modifier is then inserted into the useful heat gain equation :

- qA=FR[IsolarKτα,net(τα)n−UL(Tin−Tair)]

Equation (4) is also modified accordingly.

- η=FRKτα,net(τα)n−FRUL(Tin−Tair)Isolar

### Outlet Temperature[LINK]

Outlet temperature is calculated using the useful heat gain q as determined by Equation , the inlet fluid temperature T_{in}, and the mass flow rate available from the plant simulation:

- qA=˙mcp(Tout−Tin)

where

˙m = fluid mass flow rate through the collector

c_{p} = specific heat of the working fluid

Solving for T_{out},

- Tout=Tin+q˙mcpA

If there is no flow through the collector, T_{out} is the stagnation temperature of the fluid. This is calculated by setting the left side of Equation (13) to zero and solving for T_{in} (which also equals T_{out} for the no flow case).

ASHRAE. 1989. ASHRAE Standard 96-1980 (RA 89): Methods of Testing to Determine the Thermal Performance of Unglazed Flat-Plate Liquid-Type Solar Collectors. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

ASHRAE. 1991. ASHRAE Standard 93-1986 (RA 91): Methods of Testing to Determine the Thermal Performance of Solar Collectors. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

Duffie, J. A., and Beckman, W. A. 1991. Solar Engineering of Thermal Processes, Second Edition. New York: Wiley-Interscience.

Solar Rating and Certification Corporation. 2004. Directory of SRCC Certified Solar Collector Ratings, OG 100. Cocoa, Florida: Solar Rating and Certification Corporation.

## Integral-collector-storage (ICS) Solar Collector[LINK]

Solar collectors with integral storage unit models use SolarCollector:IntegralCollectorStorage object, and the characteristics parameter inputs of this collector are provided by the SolarCollectorPerformance:IntegralCollectorStorage object. This model is based on detailed Energy Balance equations of solar collectors that integrates storage in it. This model has two options to represent the collector bottom outside boundary conditions: AmbientAir, and OtherSideConditionsModel. AmbientAir simply applies outside air temperature using combined convection and radiation conductance, and the OtherSideConditionsModel applies combined radiation and convection models that exiats in a naturally ventilated cavity to represent the collector bottom outside boundary condition. The later boundary condition accounts for the shading of the collector on the underlying surface, hence, the ICS collector can be assumed as an integral part of the building envelope. Schamtic diagram of a rectangular ICS solar collector is shown in Figure 295 below:

### Solar and Shading Calculations[LINK]

The solar collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account. Likewise, the collector surface shades the roof surface beneath it, hence no direct solar radiation incident on the roof surface. The collector and the roof outside boundary conditions should be specified as OtherSideConditionModel to account for solar collector shading impact on the roof surface.

### Mathematical Model[LINK]

The integral-collector-storage (ICS) solar collector is represented using two transient energy balance equations shown below. These equations represent the energy balance equation for the absorber plate, and the water in the collector.

mpCpdTpdt=A⋅[(τα)e⋅It−hpw(Tp−Tw)−Ut(Tp−Ta)]

mwCwdTwdt=A⋅[hpw(Tp−Tw)−Ub(Tw−TOSC)−Us(Tw−Ta)]−˙mwCw(Tw−Twi)

Where,

*m*_{p}C_{p} = thermal capacity of the absorber surface, J/°C

A = collector gross area, m^{2}

(τα)e = transmittance-absorptance product of the absorber plate and cover system

*I*_{t} = total solar irradiation, (W/m^{2})

*h*_{pw} = convective heat transfer coefficient from absorber plate to water, (W/m2°K)

*U*_{t} = overall heat loss coefficient from absorber to the ambient air, (W/m2 °K)

*T*_{p} = absorber plate average temperature, (°C)

*T*_{w} = collector water average temperature, (°C)

*T*_{a} = ambient air temperature, (°C)

*m*_{w}C_{pw} = thermal capacity of the water mass in the collector, (J/°C)

*U*_{s} = area-weighted conductance of the collector side insulation, (W/m^{2}°K)

*U*_{b} = conductance of the collector bottom insulation, (W/m^{2}°K)

*T*_{osc} = outside temperature of bottom insulation determined from the other side condition model, (°C)

*T*_{wi} = Entering makeup or mains water temperature, (°C)

˙mwCw = water capacity flow through the collector, (W/°C)

The other side condition model boundary condition represented by the *T*_{osc}, allows us to apply a realistic outside boundary condition for a collector mounted on a building roof. This also accounts for the shading impact of the collector on the under-laying surface (roof). On the other hand if ambient air boundary condition is specified, then the collector does not shade the underlying surface it is mounted on.

The two energy balance equation can be written as non-homogeneous first order DE with constant coefficients. The initial conditions for these equations are the absorber plate average temperature and the collector water average temperature at previous time steps.

dTpdt=a1Tp+a2Tw+a3

dTwdt=b1Tp+b2Tw+b3

a1=−(A⋅hpw+A⋅Ut)/(mp⋅Cp)

a2=A⋅hpwTw/(mp⋅Cp)

a3=A⋅(τα)e⋅It+A⋅UtTa

b1=A⋅hpwTp/(mw⋅Cw)

b2=−(A⋅hpw+A⋅Ub+A⋅Us+˙mwCw)

b3=(A⋅UbTosc+A⋅UsTa+˙mwCwTwi)

The two coupled first order differential equation are solved analytically. Auxiliary equation of the the coupled homogeneous differential equation is given by:

λ2−(a1+b2)+(a1b2−a2b1)=0

This auxiliary quadratic equation has always two distinct real roots (_{1} and _{2}) hence the solution of the homogeneous equation is exponential, and the general solutions of the differential equations are given by:

Tp=c1el1t+c2el2t+A

Tw=r1c1el1t+r2c2el2t+B

The constant terms *A* and *B* are the particular solution of the non-homogeneous differential equations, the coefficients of the exponential terms (*c*_{1}, *c*_{2}, *r*_{1}, and *r*_{2}) are determined from the initial conditions of the absorber and collector water temperatures (*T*_{p0}, *T*_{w0}) and are given by:

r1=(l1−a1)/a2;r2=(l2−a1)/a2

A=(−a3b2+b3a2)/(a1b2−b1a2);B=(−a1b3+b1a3)/(a1b2−b1a2)

c1=(r2Tp0−Tw0−r2A+B)/(r2−r1);c2=(Tw0−r1Tp0+r1A−B)/(r2−r1)

### Thermal Network Model:[LINK]

The thermal network model requires energy balance for each of the collector covers as well. The heat balance equation of the collector covers is assumed to obey steady state formulation by ignoring their thermal mass. The thermal-network representation of the ICS collector is shown in Figure 296. Also, the heat balance at each cover surface requires knowledge of the amount of solar fraction absorbed, which is determined from the ray tracing analysis. For the thermal network model shown above the overall top heat loss coefficient is determined from combination of the resistances in series as follows:

Ut=[R1+R2+R3]−1

Or

Ut=[1hc,c1−a+hr,c1−a+1hc,c2−c1+hr,c2−c1+1hc,p−c2+hr,p−c2]−1

The convection and radiation heat transfer coefficients in equation above are calculated based on temperatures at the previous time step and determined as described in the *Heat Transfer Coefficients* section.

**Collector Cover Heat Balance**

Ignoring the thermal mass of the collector cover, a steady state heat balance equations are formulated for each cover that allow us to determine cover temperatures. The cover surface heat balance representation is shown in Figure 297 below.

The steady state cover heat balance equation is given by:

qLWR,1+qCONV,1+qsolar,abs+qLWR,2+qCONV,2=0

Linearizing the longwave radiation exchange and representing the convection terms using the classical equation for Newton’s law of cooling, the equations for the temperatures of covers 1 and 2 are given by:

Tc1=αc1It+hr,c1−aTa+hc,c1−aTa+hr,c2−c1Tc2+hc,c2−c1Tc2hr,c1−a+hc,c1−a+hr,c2−c1+hc,c2−c1

Tc2=αc2It+hr,c2−c1Tc1+hc,c2−c1Tc1+hr,p−c2Tp+hc,p−c2Tphr,c2−c1+hc,c2−c1+hr,p−c2+hc,p−c2

Where,

_{c} = the weighted average solar absorptance of covers *1* and *2*, (-)

*h*_{r,c1-a} = adjusted radiation heat transfer coefficient between cover *1* and the ambient air, (W/m^{2}K)

*h*_{c,c1-a} = convection heat transfer coefficient between cover *1* and the ambient, (W/m^{2}K)

*h*_{r,c2-c1} = radiation heat transfer coefficient between covers *1* and *2*, (W/m^{2}K)

*h*_{c,c2-c1} = convection heat transfer coefficient between covers *1* and *2*, (W/m^{2}K)

*h*_{r,p-c2} = radiation heat transfer coefficient between covers *2* and the absorber plate, (W/m^{2}K)

*h*_{c,p-c2} = convection heat transfer coefficient between covers *2* and the absorber plate, (W/m^{2}K)

*q*_{LWR,1} = longwave radiation exchange flux on side *1* of the collector cover, (W/m^{2})

*q*_{CONV,1} = convection heat flux on side *1* of the collector cover, (W/m^{2})

*q*_{LWR,2} = longwave radiation exchange flux on side *2* of the collector cover, (W/m^{2})

*q*_{CONV,2} = convection heat flux on side *2* of the collector cover, (W/m^{2})

*q*_{solar,abs} = net solar radiation absorbed by the collector cover, (W/m^{2})

*R* = thermal resistance for each section along the heat flow path, (m^{2}K/W)

### Other Side Condition Model[LINK]

ICS Solar Collectors are commonly mounted on building heat transfer surfaces hence the collectors shade the underlying heat transfer surface and require a unique boundary condition that reflects the air cavity environment created between the bottom of the collector surface and the underlying surface. The other side condition model that allows us to estimate the other side temperature, *T*_{osc}, can be determined based on steady state heat balance using the known collector water temperature at the previous time step.

Ignoring thermal mass of the collector bottom insulation, steady state surface heat balance can be formulated on the outer plane of the collector bottom surface facing the cavity as shown in Figure 4. The heat balance equation on the outer plane of the collector bottom surface is given by:

qcond+qconv,cav+qrad,cav=0

Substituting the equations for each term in the above equation yields:

UL(Tw−Tosc)+hc,cav(Ta,cav−Tosc)+hr,cav(Tso−Tosc)=0

Simplifying yields the bottom insulation other side condition temperature:

Tosc=ULTw+hc,cavTa,cav+hr,cavTsoUL+hc,cav+hr,cav

The cavity air temperature is determined from cavity air heat balance as follows:

Ta,cav=hc,cavATosc+˙mventCpTa+hc,cavATsohc,cavA+˙mventCp+hc,cavA

Where

*h*_{r,cav} = linearized radiation coefficient for underlying surface in the cavity, (W/m^{2}K)

*h*_{c,cav} = convection coefficient for underlying surface in the cavity, (W/m^{2}K)

*T*_{so} = the outside face temperature of the underlying heat transfer surface, (ºC)

˙mvent = air mass flow rate due to natural ventilation, (kg/s)

*q*_{cond} = conduction heat flux though the insulation bottom and, (W/m^{2})

*q*_{conv,cav} = convection heat flux between the collector bottom outside surface and the cavity air, (W/m^{2})

*q*_{rad,cav} = longwave radiation exchange flux between the collector bottom outside surface and the outside surface of the underlying surface, (W/m^{2})

The cavity air temperature is determined from the cavity air energy balance. The air heat balance requires the ventilated cavity air natural ventilation rates. The calculation of the ventilation rate is described else where in this document. The SurfaceProperty:ExteriorNaturalVentedCavity, object is required to describe the surface properties, the characteristics of the cavity and opening for natural ventilation.

### Heat Transfer Coefficients[LINK]

The equations used to determine for the various heat transfer coefficients in the absorber and water heat balance equations are given below. The absorbed solar energy is transferred to the water by convection. Assuming natural convection dominated heat transfer for a hot surface facing down and a clod surface facing down the following correlation for Nusselt number by Fujii and Imura (1972). The Nusselt number for hot surface facing down ward is given by:

Nu=0.56(Gr⋅Pr⋅cosθ)1/5105<Gr⋅Pr<1011

The Nusselt number for hot surface upward and cold surface facing down is given by:

Nu=0.13(Gr⋅Pr)1/3Gr⋅Pr<5.0×108

Nu=0.16(Gr⋅Pr)1/3Gr⋅Pr>5.0×108

Gr=gβv(Tp−Tw)L3c/ν2

Pr=ν/α

Tr=Tp−0.25(Tp−Tw)

hw=Nu⋅k/Lc

Where,

θ = angle of inclination of the collector to the vertical, radians

*g* = gravitation force constant, 9.806 (m/s^{2})

*T*_{r} = reference properties where the thermo-physical properties are calculated, (°C)

*L*_{c} = characteristic length for the absorber plate, (m)

*k* = thermal conductivity of water at reference temperature, (W/mK)

ν = kinematic viscosity of water at reference temperature, (m*^{2}/s)

α = thermal diffusivity of water at reference temperature, (m*^{2}/s)

*β*_{v} = volumetric expansion coefficient evaluated at Tv, Tv =Tw+0.25(Tp-Tw), (K-1)

*Nu* = Nusselt number calculated for water properties at the reference temperature, (-)

*Gr* = Grashof number calculated for water properties at the reference temperature, (-)

*Pr* = Prandtle number calculated for water properties at the reference temperature, (-)

The various radiation and convection heat transfer coefficients are given by the following equations. The convection heat transfer coefficients between the covers and the absorber plate are estimated from the empirical correlation for the Nusselt number for air gap between two parallel plates developed by Hollands et al. (1976) is:

Nua=1+1.44{1−1708(sin1.8β)1.6Racosβ}{1−1708Racosβ}++{(Racosβ5830)1/3−1}+

hc=Nu⋅k/L

hrp−c2=σ(Tp+Tc2)(T2p+T2c2)1/εp+1/εc2−1

hrc1−c2=σ(Tc1+Tc2)(T2c1+T2c2)1/εc1+1/εc2−1

The long wave radiation exchange coefficients between the outer collector cover and the sky and ground referencing the ambient air temperature for mathematical simplification are given.

hrc1−s=Fsεc1σ(Tc1+Ts)(T2c1+T2s)(Tc1−Ts)(Tc1−Ta)

hrc1−g=Fgεc1σ(Tc1+Tg)(T2c1+T2g)(Tc1−Tg)(Tc1−Ta)

hcc1−a=hcc1−s+hcc1−g

The convection heat transfer coefficient from the outer cover to the surrounding air is given by:

hcc1−a=2.8+3.0Vw

When the bottom surface boundary condition is AmbientAir, the combined conductance from the outer cover to the surrounding is calculated from the equation below (Duffie and Beckman, 1991).

hcomb=5.7+3.8Vw

The overall loss coefficient through the bottom and side of the collector-storage is estimated as follows:

Ub=ULb(Ab/A)

Us=[1ULs(As/A)+1hcomb]−1

Where,

_{c1} = thermal emissivity of collector cover *1*, (-)

_{c2} = thermal emissivity of collector cover *2*, (-)

*F*_{s} = view factor from the collector to the sky, (-)

*F*_{g} = view factor from the collector to the ground, (-)

*T*_{c1} = temperature of collector cover *1*, (K)

*T*_{c2} = temperature of collector cover *2*, (K)

*T*_{s} = sky temperature, (K)

*T*_{g}= ground temperature, (K)

*k* = thermal conductivity of air, (W/m K)

*L* = air gap between the covers, (m)

*β* = inclination of the plates or covers to the horizontal, (radian)

*V*_{w} = wind speed, (m/s)

*U*_{Lb} = user specified bottom heat loss conductance, W/m^{2}K

*U*_{Ls} = user specified side heat loss conductance, W/m^{2}K

*A*_{b} = collector bottom heat transfer area, m^{2}

*A*_{s} = collector side area, m^{2}

*h*_{comb} = combined conductance from the outer cover to the ambient air, W/m^{2}K

**Transmittance-Absorptance Product**

The transmittance-absorptance product of solar collector is determined using ray tracing method for any incident angle (Duffie and Beckman, 1991). This requires optical properties of the cover and absorber materials and the the transmittance-absorptance product for any incident angle is given by:

(τα)θ=τα1−(1−α)ρd

The transmittance of the cover system for single and two cover are given by:

τ=12[(τ1⋅τ21−ρ1ρ2)⊥+(τ1⋅τ21−ρ1ρ2)∥]

ρ=12[(ρ1+τ⋅ρ2⋅τ1τ2)⊥+(ρ1+τ⋅ρ2⋅τ1τ2)∥]

The effective transmittance, reflectance and absorptance of a single cover are given by:

τ=τa2⎧⎨⎩1−r⊥1+r⊥[1−r2⊥1−(r⊥τa)2]+1−r∥1+r∥⎡⎣1−r2∥1−(r∥τa)2⎤⎦⎫⎬⎭

ρ=12⎧⎨⎩[r⊥+(1−r⊥)2τ2ar⊥1−(r⊥τa)2]+⎡⎣r∥+(1−r∥)2τ2ar∥1−(r∥τa)2⎤⎦⎫⎬⎭

α=(1−τa)2{(1−r⊥1−r⊥τa)+(1−r∥1−r∥τa)}

The transmittance of the cover system with absorption only considered _{a}, is defined as:

τa=exp(−KLcosθ2)

θ2=sin−1(sinθ1⋅n1n2)

The reflectance of un-polarized radiation on passing from medium 1 with reflective index *n*_{1} to medium 2 with reflective index *n*_{2} is given by:

r⊥=sin2(θ2−θ1)sin2(θ2+θ1)

r∥=tan2(θ2−θ1)tan2(θ2+θ1)

The sky and ground reflected diffuse radiations average equivalent incident angles are approximated by Brandemuehl and Beckman correlation (Duffie and Beckman, 1991) as follows:

θsd=59.68−0.1388β+0.001497β2

θgd=90−0.5788β+0.002693β2

where,

τ = transmittance of the cover system, (-)

τ1 = transmittance of the cover 1, (-)

τ2 = transmittance of the cover 2, (-)

α = absorptance of the absorber plate, (-)

ρd = diffuse reflectance of the inner cover, (-)

*L* = thickness of a cover material, (m)

*K* = extinction coefficient of a cover material, (m^{-1})

θ1 = angle of incidence, degree

θ2 = angle of refraction, degree

r⊥ = parallel component of reflected un-polarized radiation, (-)

r∥ = perpendicular component of reflected un-polarized radiation, (-)

β = slope of the collector, degree

θsd = equivalent incident angle for sky diffuse solar radiation, degree

θgd = equivalent incident angle for ground diffuse solar radiation, degree

The integral collector storage unit thermal performance parameters are calculated as follows:

Qdelivered=˙mwCw(Tw−Twi)

QStored=mwCwdTwdt

QSkinLoss=A⋅[Ut(Tp−Tw)+Ub(Tw−TOSC)+Us(Tw−Ta)]

ηthermal=mwCwdTwdt+˙mwCw(Tw−Twi)A⋅It

Duffie, J.A., and W.A. Beckman. 1991. *Solar Engineering of Thermal Processes,* 2d ed. New York: John Wiley & Sons.

Kumar, R. and M.A. Rosen. Thermal performance of integrated collector storage solar water heater with corrugated absorber surface. Applied Thermal Engineering: 30 (2010) 1764-1768.

Fujii, T., and H. Imura. Natural convection heat transfer from aplate with arbitrary inclination. International Journal of Heat and Mass Transfer: 15(4), (1972), 755-764.

## Photovoltaic Thermal Flat-Plate Solar Collectors[LINK]

Photovoltaic-Thermal solar collectors (PVT) combine solar electric cells and thermal working fluid to collect both electricity and heat. Athough there are currently comparatively few commercial products, PVT research has been conducted for the past 30 years and many different types of collectors have been studied. Zondag (2008) and Charalambous et. al (2007) provide reviews of the PVT literature. Because PVT is much less commercially-mature, there are no standards or rating systems such as for thermal-only, hot-water collectors. EnergyPlus currently has one simple model based on user-defined efficiencies but a more detailed model based on first-principles and a detailed, layer-by-layer description is under development.

The PVT models reuse the PV models for electrical production. These are described elsewhere in this document in the section on Photovoltaic Arrays-Simple Model

### Simple PVT Thermal Model[LINK]

The input object SolarCollector:FlatPlate:PhotovoltaicThermal provides a simple PVT model that is provided for quick use during design or policy studies. The user simply provides values for a thermal efficiency and the incident solar heats the working fuild. The model also includes a cooling mode for air-based systems where a user-provided surface emmittance is used to model cooling of the working fluid to the night sky (water-based cooling will be made available once a chilled water storage tank is available). No other details of the PVT collector’s construction are required as input data.

The simple model can heat either air or liquid. If it heats air, then the PVT is part of HVAC air system loop with air nodes connected to an air system. If it heats liquid, then the PVT is part of plant loop with nodes connected to a plant loop and the plant operating scheme determines flows.

Air-system-based PVT modeling include a modulating bypass damper arrangement. Control logic decides if the air should bypass the collector to better meet setpoint. The model requires a drybulb temperature setpoint be placed on the outlet node. The model assume the collector is intended and available for heating when the incident solar is greater than 0.3 W/m^{2} and otherwise it is intended for cooling. The inlet temperature is compare to the setpoint on the outlet node to determine if cooling or heating are beneficial. If they are, then the PVT thermal models are applied to condition the air stream. If they are not beneficial, then the PVT is completely bypassed and the inlet node is passed directly to the outlet node to model a completely bypassed damper arrangement. A report variable is available for bypass damper status.

Plant-based PVT do not include a bypass (although one could be used in the plant loop). The collector requests its design flow rate but it otherwise relies on the larger plant loop for control.

When the PVT themal collector is controlled to be “on,” in heating mode, and working fluid is flowing, the model calculates the outlet temperature based on the inlet temperature and the collected heat using the following equations.

Qtherm=Asurf⋅factiv⋅GT⋅ηthermal

where,

Qtherm is the thermal energy collected [W]

Asurf is the net area of the surface [m^{2}]

factiv is the fraction of surface aire with active PV/T collector, and

ηthermal is the thermal conversion efficiency.

Tout=Tin+Qtherm˙mcp

where,

Tout is the temperature of the working fluid leaving the PV/T

Tin is the temperature of the working fluid entering the PV/T

˙m is the entire mass flow of the working fluid through the PV/T

cp is the specific heat of the working fluid.

For air-based systems, the value of Tout is then compared to the temperature setpoint on the outlet node. If Tout exceeds the desired outlet temperature, Tset,out, then a bypass fraction is calculated to model a modulating bypass damper using:

fbypass=(Tset,out−Tout)(Tin−Tout)

When the PVT themal collector is controlled to be “on,” in cooling mode, and working fluid is flowing, the model calculates the outlet temperature based on the inlet temperature and the heat radiated and convected to the ambient using a heat balance on the outside face of the collector:

˙mcp(Tin−Tout)=˙QLWR+˙Qconv

Where,

˙QLWR is the net rate of long wavelength (thermal) radiation exchange with the air, night sky, and ground. See the section “External Longwave Radiation” in the Outside Surface Heat Balance, for full discussion of how this modeled in EnergyPlus using linearized radiation coefficients.

˙Qconv is the net rate of convective flux exchange with outdoor air. See the section “Exterior/External Convection” in the Outside Surface Heat Balance, for full discussion of how this modeled in EnergyPlus. The surface roughness is assumed to be “very smooth.”

The simple model assumes that the effective collector temperature, Tcol, is the average of the working fluid inlet and outlet temperatures so that we can make the following substitution:

Tout=2Tcol−Tin

Substituting and solving for Tcol we obtain the following model for collector temperatures during a (possible) cooling process:

Tcol=2˙mcpTin+Asurffactiv(hr,gndTgnd+hr,skyTsky+hr,airTair+hc,extTair)2˙mcp+Asurffactiv(hr,gnd+hr,sky+hr,air+hc,ext)

Then the outlet temperature can be calculated and heat losses determined. However, the model allows only sensible cooling of the air stream and limits the outlet temperature to not go below the dewpoint temperature of the inlet.

PVT collectors have a design volume flow rate for the working fluid that is autosizable. For air-based systems used as pre-conditioners, the volume flow rate is sized to meet the maximum outdoor air flow rate. For water-based systems on the supply side of a plant loop, each of the PVT collectors are sized to the overall loop flow rate. For water-based systems on the demand side of a plant loop, the collectors are sized using a rule-of-thumb for typical flow rates per unit of collector area. This rule-of-thumb is based on a constant factor of 1.905x10^{-5} m^{3}/s-m^{2}^^that was developed by analyzing SRCC data set for conventional solar collectors (see data set SolarCollectors.idf) and averaging the ratio for all 171 different collectors.

Charalambous, P.G., Maidment, G.G., Kalagirou, S.A., and Yiakoumetti, K., Photovoltaic thermal (PV/T) collectors: A review. Applied Thermal Engineering 27 (2007) 275-286.

Zondag, H.A. 2008. Flat-plate PV-Thermal collectors and systems: A review. Renewable and Sustainable Energy Reviews 12 (2008) 891-959.

## Unglazed Transpired Solar Collectors[LINK]

The input object SolarCollector:UnglazedTranspired provides a model for transpired collectors that are perhaps one of the most efficient ways to collect solar energy with demonstrated instantaneous efficiencies of over 90% and average efficiencies of over 70%. They are used for preheating outdoor air needed for ventilation and processes such as crop drying.

In EnergyPlus, an unglazed transpired solar collector (UTSC) is modeled as a special component attached to the outside face of a heat transfer surface that is also connected to the outdoor air path. A UTSC affects both the thermal envelope and the HVAC air system. From the air system’s point of view, a UTSC is heat exchanger and the modeling needs to determine how much the device raises the temperature of the outdoor air. From the thermal envelope’s point of view, the presence of the collector on the outside of the surface modifies the conditions experienced by the underlying heat transfer surfaces. EnergyPlus models building performance throughout the year and the UTSC will often be “off” in terms of forced airflow, but the collector is still present. When the UTSC is “on” there is suction airflow that is assumed to be uniform across the face. When the UTSC is “off” the collector acts as a radiation and convection baffle situated between the exterior environment and the outside face of the underlying heat transfer surface. We distinguish these two modes of operation as *active* or *passive* and model the UTSC component differently depending on which of these modes it is in.

### Heat Exchanger Effectiveness[LINK]

The perforated absorber plate is treated as a heat exchanger and modeled using a traditional effectiveness formulation. The heat exchanger effectiveness, εHX , is determined from correlations derived from small-scale experiments. Two correlations available in the literature are implemented in EnergyPlus. The first is based on research by Kutscher at the National Renewable Energy Laboratory. The second is based on the research by Van Decker, Hollands, and Brunger at the University of Waterloo. Because both correlations are considered valid, the choice of which correlation to use is left to the user.

### Kutscher Correlation[LINK]

Kutscher’s (1994) correlation encompasses surface convection between the collector and the incoming outdoor air stream that occurs on the front face, in the holes, and along the back face of the collector. The correlation uses a Reynolds number based on the hole diameter as a length scale and the mean velocity of air as it passes through the holes as the velocity scale:

ReD=VhDν

where,

Vh is the velocity through the holes [m/s]

D is the hole diameter [m]

ν is the kinematic viscosity of air [m^{2}/s]

The correlation is a function of Reynolds number, hole geometry, the free stream air velocity, and velocity through the holes:

NuD=2.75[(PD)−1.2Re0.43D+0.011σReD(U∞Vh)0.48]

where,

P is the pitch, or distance between holes, [m],

D is the diameter of the holes, [m],

σ is the porosity, or area fraction of the holes, [dimensionless],

Vh is the mean velocity of air passing through the holes, [m/s],

U∞ is the free stream velocity (local wind speed) [m/s].

The Nusselt number is formulated as:

NuD=UDk

where,

U is the overall heat transfer coefficient based on log mean temperature difference, [W/m^{2}·K], and

k is the thermal conductivity of air [W/m·K].

The heat exchanger effectiveness is:

εHX=1−e[−UA˙mcp]

Kutscher’s relation was formulated for triangular hole layout, but based on Van Decker et al. (2001) we allow using the correlation for square hole layout and scale P by a factor of 1.6.

### Van Decker, Hollands, and Brunger Correlation[LINK]

Van Decker et. al. extended Kutscher’s measurements to include a wider range of collector parameters including plate thickness, pitch, suction velocities, and square hole patterns. Their model formulation differs from Kutscher’s in that the model was built up from separate effectiveness models for the front, back, and holes of the collector. Their published correlation is:

εHX=[1−(1+ResMax(1.733Re−1/2w,0.02136)−1)]×[1−(1+0.2273Re1/2b)−1]×e(−0.01895PD−20.62ReDtD)

where,

Res=VsPv

Rew=U∞Pv

Reb=VhPv

Vs is the average suction velocity across the front face of the collector [m/s]

t is the collector plate thickness

### Heat Exchanger Leaving Temperature[LINK]

Using either of the correlations above allows determining the heat exchanger effectiveness from known values. By definition the heat exchanger effectiveness is also:

εHX=Ta,HX−TambTs,coll−Tamb

where,

Ta,HX is the temperature of the air leaving the collector and entering the plenum [ºC]

Ts,coll is the temperature of the collector’s absorber plate, [ºC], and

Tamb is the temperature of the ambient outdoor air [ºC].

By rewriting equation to solve for Ta,HX we see that the temperature of the heated outdoor air entering the plenum can be determined once the collector surface temperature is known,

Ta,HX=εHXTs,coll+(1−εHX)Tamb

### Collector Heat Balance[LINK]

The collector is assumed to be sufficiently thin and high-conductivity so that it can be modeled using a single temperature (for both sides and along its area). This temperature Ts,coll is determined by formulating a heat balance on a control volume that just encapsulates the collector surface. The heat balances are formulated separately for active and passive modes and are diagrammed in the following figure.

Observe that for the passive case, we do not use the heat exchanger relations to directly model the interaction of ventilating air with the collector. This is because these relations are considered to not apply when the UTSC is in passive mode. They were developed for uni-directional flow (rather than the balanced-in-and-out flow expected from natural forces) and for specific ranges of suction face velocity. Therefore, this heat transfer mechanism is handled using classical surface convection models (as if the collector was not perforated). (Air exchanges are modeled as ventilation in the plenum air heat balance but do not interact with the hole edges in the collector surface.)

When the UTSC is active, the heat balance on the collector surface control volume is:

q′′αsol+q′′LWR,Env+q′′conv,wind−q′′HX+q′′LWR,plen+q′′source=0

where:

q′′αsol is absorbed direct and diffuse solar (short wavelength) radiation heat flux.

q′′LWR,Env is net long wavelength (thermal) radiation flux exchange with the air and surroundings.

q′′conv,wind is surface convection flux exchange with outdoor air under high wind and rain conditions. Note that this term is usually assumed to be zero in UTSC model development but we add the term to allow for deteriorated performance of the UTSC under poor conditions.

q′′HX is heat exchanger flux from collector to incoming outdoor air.

q′′LWR,plen is net long wavelength (thermal) radiation flux exchange with the outside face of the underlying surface(s).

q′′source is a source/sink term that accounts for energy exported out of the control volume when the collecter’s absorber plate is a hybrid device such as a photovoltaic panel.

While the heat balance on the passive collector surface control volume is:

q′′αsol+q′′LWR,Env+q′′conv,Env+q′′LWR,plen+q′′conv,plen+q′′source=0

where:

q′′conv,Env = surface convection flux exchange with outdoor air.

q′′conv,plen = surface convection flux exchange with plenum air.

All terms are positive for net flux to the collector except the heat exchanger term, which is taken to be positive in the direction from collector to incoming air stream. Each of these heat balance components is introduced briefly below.

### External SW Radiation[LINK]

q′′αsol is calculated using procedures presented elsewhere in this manual and includes both direct and diffuse incident solar radiation absorbed by the surface face. This is influenced by location, surface facing angle and tilt, shading surfaces, surface face material properties, weather conditions, etc.

### External LW Radiation[LINK]

q′′LWR,Env is a standard radiation exchange formulation between the surface, the sky, the ground, and the atmosphere. The radiation heat flux is calculated from the surface absorptivity, surface temperature, sky, air, and ground temperatures, and sky and ground view factors. Radiation is modeled using linearized coefficients.

### External Convection[LINK]

q′′conv,Env is modeled using the classical formulation: q′′conv = h_{co} (T_{air} - T_{o}) where h_{co}, is the convection coefficient. This coefficient will differ depending on whether or not the UTSC is active or passive. When the UTSC is passive, h_{co} is treated in the same way as an outside face with ExteriorEnvironment conditions. When the UTSC is active, the special suction airflow situation of a transpired collector during operation means that h_{co}is often zero because the suction situation can eliminate mass transport away from the collector. However when the winds are high, the strong turbulence and highly varying pressures can cause the suction flow situation to breakdown. Therefore, we include the q′′conv,wind term in the heat balance and use a special coefficient hc,wind to model this lost heat transfer. In addition, when it is raining outside, we assume the collector gets wet and model the enhanced surface heat transfer using a large value for hc,wind.

### Heat Exchanger[LINK]

q′′HX is modeled using the classical formulation q′′HX=˙mcp(Ta,HX−Tamb)A where Ta,HX is determined using correlations described above. When the UTSC is active, the air mass flow is determined from the operation of the outdoor air mixer component. When the UTSC is off, this term is zero.

### Plenum LW Radation[LINK]

q′′LWR,plen is a standard radiation exchange formulation between the collector surface and the underlying heat transfer surface located across the plenum. Radiation is modeled using linearized coefficients.

### Plenum Convection[LINK]

q′′conv,plen is modeled using the classical formulation: q′′conv = h_{cp} (T_{air} - T_{o}) where h_{cp}, is the convection coefficient. This coefficient is taken as zero when the UTSC is operating because of the suction airflow situation. When the UTSC is off, the value for h_{cp} is obtained from correlations used for window gaps from ISO (2003) standard 15099.

Substituting models into and solving for Ts,coll yields the following equation when the UTSC is active (“on”):

Ts,coll=(Isα+hr,atmTamb+hr,skyTsky+hr,gndTamb+hr,plenTso+hc,windTamb+˙mcpATamb−˙mcpA(1−εHX)Tamb+q′′source)(hr,atm+hr,sky+hr,gnd+hr,plen+hc,wind+˙mcpAεHX)

and substituting into yields the following equation when the UTSC is passive (“off”):

Ts,coll=(Isα+hcoTamb+hr,atmTamb+hr,skyTsky+hr,gndTamb+hr,plenTso+hc,plenTa,plen+q′′source)(hco+hr,air+hr,sky+hr,gnd+hr,plen+hc,plen)

where,

Is is the incident solar radiation of all types [W/m^{2}],

α is the solar absorptivity of the collector [dimensionless],

hr,atm is the linearized radiation coefficient for the surrounding atmosphere [W/m^{2}·K],

Tamb is the outdoor drybulb from the weather data, also assumed for ground surface [ºC],

hr,sky is the linearized radiation coefficient for the sky [W/m^{2}·K],

Tsky is the effective sky temperature [ºC],

hr,gnd is the linearized radiation coefficient for the ground [W/m^{2}·K],

hr,plen is the linearized radiation coefficient for the underlying surface [W/m^{2}·K],

hc,wind is the convection coefficient for the outdoor environment when the UTSC is active and winds are high or it is raining [W/m^{2}·K],

Tso is the temperature of the outside face of the underlying heat transfer surface [ºC],

˙m is the air mass flow rate when in active mode [kg/s],

cp is the specific heat of air at constant pressure [J/kg·K],

A is the area of the collector [m^{2}],

hco is the convection coefficient for the outdoor environment [W/m^{2}·K],

hc,plen is the convection coefficient for the surfaces facing the plenum [W/m^{2}·K], and

Ta,plen is the air drybulb for air in the plenum and entering the outdoor air system [ºC].

### Plenum Heat Balance[LINK]

The *plenum* is the volume of air located between the collector and the underlying heat transfer surface. The plenum air is modeled as well-mixed. The uniform temperature of the plenum air, Ta,plen, is determined by formulating a heat balance on a control volume of air as diagrammed below.

Note that we have formulated the control volumes with slight differences for the active and passive cases. For the active case, the suction air situation and heat exchanger effectiveness formulations dictate that the collector surface control volume encompass part of the air adjacent to both the front and back surfaces of the collector. However for the passive case, the collector surface control volume has no air in it and the plenum air control volume extends all the way to the surface of the collector.

When the UTSC is active, the heat balance on the plenum air control volume is:

˙Qair+˙Qco=0

where,

˙Qair is the net rate of energy added by suction air convecting through the control volume.

˙Qco is the net rate of energy added by surface convection heat transfer with the underlying surface.

When the UTSC is passive, the heat balance on the plenum air control volume is:

˙Qvent+˙Qco+˙Qc,coll=0

where,

˙Qvent is the net rate of energy added from infiltration - where outdoor ambient air exchanges with the plenum air.

˙Qc,coll is the net rate of energy added by surface convection heat transfer with the collector.

Substituting into and solving for Ta,plen yields the following equation for when the UTSC is active:

Ta,plen=(˙mcpTa,HX+hc,plenATso)(˙mcp+hc,plenA)

And substituting into yields the following equation when the UTSC is passive:

Ta,plen=(hc,plenATso+˙mventcpTamb+hc,plenATs,coll)(hc,plenA+˙mventcp+hc,plenA)

where,

˙mvent is the air mass flow from natural forces [kg/s]

The literature on UTSC does not appear to address the passive mode of operation and no models for ˙mvent have been identified. Nevertheless, natural buoyancy and wind forces are expected to drive air exchange between the plenum and ambient and some method of modeling ˙mvent is needed. Reasoning that the configuration is similar to single-side natural ventilation, we elect to use correlations for natural ventilation presented in Chapter 26. of ASHRAE HOF (2001).

˙mvent=ρ−˙Vtot

where,

ρ is the density of air [kg/m^{3}], and

−˙Vtot=−˙Vwind+−˙Vthermal is the total volumetric flow rate of air ventilating in and out of the plenum.

−˙Vwind=CvAinU∞

−˙Vthermal=CDAin√2gΔHNPL(Ta,plen−Tamb)/Ta,plen (if Ta,plen>Tamb)

−˙Vthermal=CDAin√2gΔHNPL(Tamb−Ta,plen)/Tamb (if Tamb>Ta,plen and UTSC vertical)

Cv is the effectiveness of the openings that depends on opening geometry and the orientation with respect to the wind. ASHRAE HoF (2001) indicates values ranging from 0.25 to 0.6. In the UTSC model, this value is available for user input and defaulted to 0.25.

CD is the discharge coefficient for the opening and depends on opening geometry. In the UTSC model, this value is available for user input and defaulted to 0.65.

Mass continuity arguments lead to modeling the area of the openings as one half of the total area of the holes, so we have:

Ain=Aσ2

g is the gravitational constant taken as 9.81 [m/s^{2}].

ΔHNPL is the height from midpoint of lower opening to the Neutral Pressure Level. This is taken as one-fourth the overall height of the UTSC if it is mounted vertically. For tilted collectors, the nominal height is modified by the sine of the tilt. If the UTSC is mounted horizontally (e.g. on the roof) then the ΔHNPL is taken as the gap thickness of the plenum.

If the UTSC is horizontal and Tamb>Ta,plen then −˙Vthermal=0 because this is a stable situation.

### Underlying Heat Transfer Surface[LINK]

The UTSC is applied to the outside of a heat transfer surface. This surface is modeled using the usual EnergyPlus methods for handling heat capacity and transients - typically the CTF method. These native EnergyPlus Heat Balance routines are used to calculate Tso. The UTSC model is coupled to the underlying surface using the OtherSideConditionsModel mechanism. The UTSC model provides values for hr,plen, Ts,coll, hc,plen, and Ta,plen for use with the Heat Balance Model calculations for the outside face of the underlying surface (described elsewhere in this manual).

### Solar and Shading Calculations[LINK]

The transpired collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account.

### Local Wind Speed Calculations[LINK]

The outdoor wind speed affects terms used in modeling UTSC components. The wind speed in the weather file is assumed to be measured at a meteorological station located in an open field at a height of 10 m. To adjust for different terrain at the building site and differences in the height of building surfaces, the local wind speed is calculated for each surface.

The wind speed is modified from the measured meteorological wind speed by the equation (ASHRAE 2001):

U∞=Vmet(δmetzmet)amet(zδ)a

where z is the height of the centroid of the UTSC, z_{met} is the height of the standard metereological wind speed measurement, and a and δ are terrain-dependent coefficients. δ is the boundary layer thickness for the given terrain type. The values of a and δ are shown in the following tables:

Terrain-Dependent Coefficients (ASHRAE 2001).
1 |
Flat, open country |
0.14 |
270 |

2 |
Rough, wooded country |
0.22 |
370 |

3 |
Towns and cities |
0.33 |
460 |

4 |
Ocean |
0.10 |
210 |

5 |
Urban, industrial, forest |
0.22 |
370 |

The UTSC can be defined such that it has multiple underlying heat transfer surfaces. The centroid heights for each surface are area-weighted to determine the average height for use in the local wind calculation.

### Convection Coefficients[LINK]

UTSC modeling requires calculating up to three different coefficients for surface convection heat transfer. These coefficients are defined in the classic way by:

hc=Tair−Tsurfq′′conv

First, hco is the convection coefficient for the collector surface facing the outdoors when the UTSC is passive. It is modeled in exactly the same way as elsewhere in EnergyPlus and will depend on the user setting for Outside Convection Algorithm - Outside Surface Heat Balance entry elsewhere in this document.

Second, hc,plen is the convection coefficient for surfaces facing the plenum. This coefficient is applied to just the underlying surface’s convection when the UTSC is active and to both the collector and the underlying surface when the UTSC is passive. When the UTSC is active, we use the convection correlation for forced air developed by McAdams (1954) as published by ASHRAE HoF (2001):

hc,plen=5.62+3.9Vp

where,

Vp is the mean velocity in the plenum determined from Vp=˙m2ρAp where Ap is the effective cross section area of the plenum perpendicular to the primary flow direction. When the UTSC is passive, we model the convection in the same way used in EnergyPlus to model air gaps in windows. These correlations vary by Rayleigh number and surface tilt and are based on the work of various research including Hollands et. al., Elsherbiny et. al., Wright, and Arnold. The formulations are documented in ISO (2003) standard 15099. For the UTSC implementation, the routines were adapted from Subroutine NusseltNumber in WindowManager.f90 (by F. Winkelmann), which itself was derived from Window5 subroutine “nusselt”.

Third, hc,wind is the convection coefficient used to degrade the UTSC performance under environmental conditions with high wind or rain. If the weather file indicates it is raining, then we set hc,wind = 1000.0 which has the effect of making the collector the same temperature as ambient air. The heat exchanger effectiveness correlations described above account for a moderate amount of wind, but the correlations appear confined to the range 0 to 5.0 m/s. Therefore we set hc,wind equal to zero if U∞ is <= 5.0 m/s. If U∞ is > 5.0 m/s then we use the McAdams correlation but with a reduced velocity magnitude:

$$ hc,wind=5.62+3.9(U∞−5.0)

### Radiation Coefficients[LINK]

UTSC modeling requires calculating up to four different linearized coefficients for radiation heat transfer. Whereas radiation calculations usually use temperature raised to the fourth power, this greatly complicates solving heat balance equations for a single temperature. Linearized radiation coefficients have the same units and are used in the same manner as surface convection coefficients and introduce very little error for the temperature levels involved.

The radiation coefficient, hr,plen , is used to model thermal radiation between the collector surface and the outside face of the underlying heat transfer surface. We assume a view factor of unity. It is calculated using:

hr,plen=σSBecolleso(T4s,coll−T4so)(Ts,coll−Tso)

where,

all temperatures are converted to Kelvin,

σSB is the Stefan-Boltzmann constant,

ecoll is the longwave thermal emmittance of the collector, and

eso is the longwave the

## Solar Collectors[LINK]

Solar collectors are devices that convert solar energy into thermal energy by raising the temperature of a circulating heat transfer fluid. The fluid can then be used to heat water for domestic hot water usage or space heating. Flat-plate solar collectors using water as the heat transfer fluid, Integral-Collector Storage solar collectors using water and unglazed transpired solar collectors using air are currently the only types of collector available in EnergyPlus.

## Flat-Plate Solar Collectors[LINK]

The input object SolarCollector:FlatPlate:Water provides a model for flat-plate solar collectors that are the most common type of collector. Standards have been established by ASHRAE for the performance testing of these collectors (ASHRAE 1989; 1991) and the Solar Rating and Certification Corporation (SRCC) publishes a directory of commercially available collectors in North America (SRCC 2003).

The EnergyPlus model is based on the equations found in the ASHRAE standards and Duffie and Beckman (1991). This model applies to glazed and unglazed flat-plate collectors, as well as banks of tubular, i.e. evacuated tube, collectors.

## Solar and Shading Calculations[LINK]

The solar collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account. Likewise, the collector surface can shade other surfaces, for example, reducing the incident radiation on the roof beneath it.

## Thermal Performance[LINK]

The thermal efficiency of a collector is defined as the ratio of the useful heat gain of the collector fluid versus the total incident solar radiation on the gross surface area of the collector.

where

q = useful heat gain

A = gross area of the collector

I

_{solar}= total incident solar radiationNotice that the efficiency is only defined for I

_{solar}> 0.An energy balance on a solar collector with double glazing shows relationships between the glazing properties, absorber plate properties, and environmental conditions.

where

_{g1}= transmittance of the first glazing layer_{g2}= transmittance of the second glazing layer_{abs}= absorptance of the absorber plateR

_{rad}= radiative resistance from absorber to inside glazingR

_{conv}= convective resistance from absorber to inside glazingR

_{cond}= conductive resistance from absorber to outdoor air through the insulationT

_{abs}= temperature of the absorber plateT

_{g2}= temperature of the inside glazingT

_{air}= temperature of the outdoor airThe equation above can be approximated with a simpler formulation as:

where

F

_{R}= an empirically determined correction factor(τα) = the product of all transmittance and absorptance terms

U

_{L}= overall heat loss coefficient combining radiation, convection, and conduction termsT

_{in}= inlet temperature of the working fluidSubstituting this into Equation (1),

A linear correlation can be constructed by treating F

_{R}(τα) and F_{R}U_{L}as characteristic constants of the solar collector:Similarly, a quadratic correlation can be constructed using the form:

Both first- and second-order efficiency equation coefficients are listed in the

Directory of SRCC Certified Solar Collector Ratings.## Incident Angle Modifiers[LINK]

As with regular windows the transmittance of the collector glazing varies with the incidence angle of radiation. Usually the transmittance is highest when the incident radiation is normal to the glazing surface. Test conditions determine the efficiency coefficients for normal incidence. For off-normal angles, the transmittance of the glazing is modified by an

incident angle modifiercoefficient.Additional testing determines the incident angle modifier as a function of incident angle. This relationship can be fit to a first-order, linear correlation:

or a second-order, quadratic correlation:

The incident angle modifier coefficients

band_{0}bare usually negative, although some collectors have a positive value for_{1}b. Both first- and second-order incident angle modifier equation coefficients are listed in the_{0}Directory of SRCC Certified Solar Collector Ratings.The SRCC incident angle modifier equation coefficients are only valid for incident angles of 60 degrees or less. Because these curves can be valid yet behave poorly for angles greater than 60 degree, the EnergyPlus model cuts off collector gains for incident angles greater than 60 degrees.

For flat-plate collectors, the incident angle modifier is generally symmetrical. However, for tubular collectors the incident angle modifier is different depending on whether the incident angle is parallel or perpendicular to the tubes. These are called bi-axial modifiers. Some special flat-plate collectors may also exhibit this asymmetry. The current model cannot yet handle two sets of incident angle modifiers. In the meantime it is recommended that tubular collectors be approximated with caution using either the parallel or perpendicular correlation.

Incident angle modifiers are calculated separately for sun, sky, and ground radiation. The net incident angle modifier for all incident radiation is calculated by weighting each component by the corresponding modifier.

For sky and ground radiation the incident angle is approximated using Brandemuehl and Beckman’s equations:

θsky=59.68−0.1388ϕ+0.001497ϕ2

θground=90.0−0.5788ϕ+0.002693ϕ2

where ϕ is the surface tilt in degrees.

The net incident angle modifier is then inserted into the useful heat gain equation :

Equation (4) is also modified accordingly.

## Outlet Temperature[LINK]

Outlet temperature is calculated using the useful heat gain q as determined by Equation , the inlet fluid temperature T

_{in}, and the mass flow rate available from the plant simulation:where

˙m = fluid mass flow rate through the collector

c

_{p}= specific heat of the working fluidSolving for T

_{out},If there is no flow through the collector, T

_{out}is the stagnation temperature of the fluid. This is calculated by setting the left side of Equation (13) to zero and solving for T_{in}(which also equals T_{out}for the no flow case).## References[LINK]

ASHRAE. 1989. ASHRAE Standard 96-1980 (RA 89): Methods of Testing to Determine the Thermal Performance of Unglazed Flat-Plate Liquid-Type Solar Collectors. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

ASHRAE. 1991. ASHRAE Standard 93-1986 (RA 91): Methods of Testing to Determine the Thermal Performance of Solar Collectors. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

Duffie, J. A., and Beckman, W. A. 1991. Solar Engineering of Thermal Processes, Second Edition. New York: Wiley-Interscience.

Solar Rating and Certification Corporation. 2004. Directory of SRCC Certified Solar Collector Ratings, OG 100. Cocoa, Florida: Solar Rating and Certification Corporation.

## Integral-collector-storage (ICS) Solar Collector[LINK]

Solar collectors with integral storage unit models use SolarCollector:IntegralCollectorStorage object, and the characteristics parameter inputs of this collector are provided by the SolarCollectorPerformance:IntegralCollectorStorage object. This model is based on detailed Energy Balance equations of solar collectors that integrates storage in it. This model has two options to represent the collector bottom outside boundary conditions: AmbientAir, and OtherSideConditionsModel. AmbientAir simply applies outside air temperature using combined convection and radiation conductance, and the OtherSideConditionsModel applies combined radiation and convection models that exiats in a naturally ventilated cavity to represent the collector bottom outside boundary condition. The later boundary condition accounts for the shading of the collector on the underlying surface, hence, the ICS collector can be assumed as an integral part of the building envelope. Schamtic diagram of a rectangular ICS solar collector is shown in Figure 295 below:

Schematic diagram of rectangular Integrated Collector Storage unit

## Solar and Shading Calculations[LINK]

The solar collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account. Likewise, the collector surface shades the roof surface beneath it, hence no direct solar radiation incident on the roof surface. The collector and the roof outside boundary conditions should be specified as OtherSideConditionModel to account for solar collector shading impact on the roof surface.

## Mathematical Model[LINK]

The integral-collector-storage (ICS) solar collector is represented using two transient energy balance equations shown below. These equations represent the energy balance equation for the absorber plate, and the water in the collector.

mpCpdTpdt=A⋅[(τα)e⋅It−hpw(Tp−Tw)−Ut(Tp−Ta)]

mwCwdTwdt=A⋅[hpw(Tp−Tw)−Ub(Tw−TOSC)−Us(Tw−Ta)]−˙mwCw(Tw−Twi)

Where,

m= thermal capacity of the absorber surface, J/°C_{p}C_{p}A = collector gross area, m

^{2}(τα)e = transmittance-absorptance product of the absorber plate and cover system

I= total solar irradiation, (W/m_{t}^{2})h= convective heat transfer coefficient from absorber plate to water, (W/m2°K)_{pw}U= overall heat loss coefficient from absorber to the ambient air, (W/m2 °K)_{t}T= absorber plate average temperature, (°C)_{p}T= collector water average temperature, (°C)_{w}T= ambient air temperature, (°C)_{a}m= thermal capacity of the water mass in the collector, (J/°C)_{w}C_{pw}U= area-weighted conductance of the collector side insulation, (W/m_{s}^{2}°K)U= conductance of the collector bottom insulation, (W/m_{b}^{2}°K)T= outside temperature of bottom insulation determined from the other side condition model, (°C)_{osc}T= Entering makeup or mains water temperature, (°C)_{wi}˙mwCw = water capacity flow through the collector, (W/°C)

The other side condition model boundary condition represented by the

T, allows us to apply a realistic outside boundary condition for a collector mounted on a building roof. This also accounts for the shading impact of the collector on the under-laying surface (roof). On the other hand if ambient air boundary condition is specified, then the collector does not shade the underlying surface it is mounted on._{osc}The two energy balance equation can be written as non-homogeneous first order DE with constant coefficients. The initial conditions for these equations are the absorber plate average temperature and the collector water average temperature at previous time steps.

dTpdt=a1Tp+a2Tw+a3

dTwdt=b1Tp+b2Tw+b3

a1=−(A⋅hpw+A⋅Ut)/(mp⋅Cp)

a2=A⋅hpwTw/(mp⋅Cp)

a3=A⋅(τα)e⋅It+A⋅UtTa

b1=A⋅hpwTp/(mw⋅Cw)

b2=−(A⋅hpw+A⋅Ub+A⋅Us+˙mwCw)

b3=(A⋅UbTosc+A⋅UsTa+˙mwCwTwi)

The two coupled first order differential equation are solved analytically. Auxiliary equation of the the coupled homogeneous differential equation is given by:

λ2−(a1+b2)+(a1b2−a2b1)=0

This auxiliary quadratic equation has always two distinct real roots (

_{1}and_{2}) hence the solution of the homogeneous equation is exponential, and the general solutions of the differential equations are given by:Tp=c1el1t+c2el2t+A

Tw=r1c1el1t+r2c2el2t+B

The constant terms

AandBare the particular solution of the non-homogeneous differential equations, the coefficients of the exponential terms (c,_{1}c,_{2}r, and_{1}r) are determined from the initial conditions of the absorber and collector water temperatures (_{2}T,_{p0}T) and are given by:_{w0}r1=(l1−a1)/a2;r2=(l2−a1)/a2

A=(−a3b2+b3a2)/(a1b2−b1a2);B=(−a1b3+b1a3)/(a1b2−b1a2)

c1=(r2Tp0−Tw0−r2A+B)/(r2−r1);c2=(Tw0−r1Tp0+r1A−B)/(r2−r1)

## Thermal Network Model:[LINK]

The thermal network model requires energy balance for each of the collector covers as well. The heat balance equation of the collector covers is assumed to obey steady state formulation by ignoring their thermal mass. The thermal-network representation of the ICS collector is shown in Figure 296. Also, the heat balance at each cover surface requires knowledge of the amount of solar fraction absorbed, which is determined from the ray tracing analysis. For the thermal network model shown above the overall top heat loss coefficient is determined from combination of the resistances in series as follows:

Ut=[R1+R2+R3]−1

Or

Ut=[1hc,c1−a+hr,c1−a+1hc,c2−c1+hr,c2−c1+1hc,p−c2+hr,p−c2]−1

The convection and radiation heat transfer coefficients in equation above are calculated based on temperatures at the previous time step and determined as described in the

Heat Transfer Coefficientssection.Thermal network diagram for ICS Solar Collector

Collector Cover Heat BalanceIgnoring the thermal mass of the collector cover, a steady state heat balance equations are formulated for each cover that allow us to determine cover temperatures. The cover surface heat balance representation is shown in Figure 297 below.

Collector Cover Surface Heat Balance

The steady state cover heat balance equation is given by:

qLWR,1+qCONV,1+qsolar,abs+qLWR,2+qCONV,2=0

Linearizing the longwave radiation exchange and representing the convection terms using the classical equation for Newton’s law of cooling, the equations for the temperatures of covers 1 and 2 are given by:

Tc1=αc1It+hr,c1−aTa+hc,c1−aTa+hr,c2−c1Tc2+hc,c2−c1Tc2hr,c1−a+hc,c1−a+hr,c2−c1+hc,c2−c1

Tc2=αc2It+hr,c2−c1Tc1+hc,c2−c1Tc1+hr,p−c2Tp+hc,p−c2Tphr,c2−c1+hc,c2−c1+hr,p−c2+hc,p−c2

Where,

= the weighted average solar absorptance of covers_{c}1and2, (-)h= adjusted radiation heat transfer coefficient between cover_{r,c1-a}1and the ambient air, (W/m^{2}K)h= convection heat transfer coefficient between cover_{c,c1-a}1and the ambient, (W/m^{2}K)h= radiation heat transfer coefficient between covers_{r,c2-c1}1and2, (W/m^{2}K)h= convection heat transfer coefficient between covers_{c,c2-c1}1and2, (W/m^{2}K)h= radiation heat transfer coefficient between covers_{r,p-c2}2and the absorber plate, (W/m^{2}K)h= convection heat transfer coefficient between covers_{c,p-c2}2and the absorber plate, (W/m^{2}K)q= longwave radiation exchange flux on side_{LWR,1}1of the collector cover, (W/m^{2})q= convection heat flux on side_{CONV,1}1of the collector cover, (W/m^{2})q= longwave radiation exchange flux on side_{LWR,2}2of the collector cover, (W/m^{2})q= convection heat flux on side_{CONV,2}2of the collector cover, (W/m^{2})q= net solar radiation absorbed by the collector cover, (W/m_{solar,abs}^{2})R= thermal resistance for each section along the heat flow path, (m^{2}K/W)## Other Side Condition Model[LINK]

ICS Solar Collectors are commonly mounted on building heat transfer surfaces hence the collectors shade the underlying heat transfer surface and require a unique boundary condition that reflects the air cavity environment created between the bottom of the collector surface and the underlying surface. The other side condition model that allows us to estimate the other side temperature,

T, can be determined based on steady state heat balance using the known collector water temperature at the previous time step._{osc}Illustration for Other Side Condition Model

Ignoring thermal mass of the collector bottom insulation, steady state surface heat balance can be formulated on the outer plane of the collector bottom surface facing the cavity as shown in Figure 4. The heat balance equation on the outer plane of the collector bottom surface is given by:

qcond+qconv,cav+qrad,cav=0

Substituting the equations for each term in the above equation yields:

UL(Tw−Tosc)+hc,cav(Ta,cav−Tosc)+hr,cav(Tso−Tosc)=0

Simplifying yields the bottom insulation other side condition temperature:

Tosc=ULTw+hc,cavTa,cav+hr,cavTsoUL+hc,cav+hr,cav

The cavity air temperature is determined from cavity air heat balance as follows:

Ta,cav=hc,cavATosc+˙mventCpTa+hc,cavATsohc,cavA+˙mventCp+hc,cavA

Where

h= linearized radiation coefficient for underlying surface in the cavity, (W/m_{r,cav}^{2}K)h= convection coefficient for underlying surface in the cavity, (W/m_{c,cav}^{2}K)T= the outside face temperature of the underlying heat transfer surface, (ºC)_{so}˙mvent = air mass flow rate due to natural ventilation, (kg/s)

q= conduction heat flux though the insulation bottom and, (W/m_{cond}^{2})q= convection heat flux between the collector bottom outside surface and the cavity air, (W/m_{conv,cav}^{2})q= longwave radiation exchange flux between the collector bottom outside surface and the outside surface of the underlying surface, (W/m_{rad,cav}^{2})The cavity air temperature is determined from the cavity air energy balance. The air heat balance requires the ventilated cavity air natural ventilation rates. The calculation of the ventilation rate is described else where in this document. The SurfaceProperty:ExteriorNaturalVentedCavity, object is required to describe the surface properties, the characteristics of the cavity and opening for natural ventilation.

## Heat Transfer Coefficients[LINK]

The equations used to determine for the various heat transfer coefficients in the absorber and water heat balance equations are given below. The absorbed solar energy is transferred to the water by convection. Assuming natural convection dominated heat transfer for a hot surface facing down and a clod surface facing down the following correlation for Nusselt number by Fujii and Imura (1972). The Nusselt number for hot surface facing down ward is given by:

Nu=0.56(Gr⋅Pr⋅cosθ)1/5105<Gr⋅Pr<1011

The Nusselt number for hot surface upward and cold surface facing down is given by:

Nu=0.13(Gr⋅Pr)1/3Gr⋅Pr<5.0×108

Nu=0.16(Gr⋅Pr)1/3Gr⋅Pr>5.0×108

Gr=gβv(Tp−Tw)L3c/ν2

Pr=ν/α

Tr=Tp−0.25(Tp−Tw)

hw=Nu⋅k/Lc

Where,

θ = angle of inclination of the collector to the vertical, radians

g= gravitation force constant, 9.806 (m/s^{2})T= reference properties where the thermo-physical properties are calculated, (°C)_{r}L= characteristic length for the absorber plate, (m)_{c}k= thermal conductivity of water at reference temperature, (W/mK)ν = kinematic viscosity of water at reference temperature, (m*

^{2}/s)α = thermal diffusivity of water at reference temperature, (m*

^{2}/s)β= volumetric expansion coefficient evaluated at Tv, Tv =Tw+0.25(Tp-Tw), (K-1)_{v}Nu= Nusselt number calculated for water properties at the reference temperature, (-)Gr= Grashof number calculated for water properties at the reference temperature, (-)Pr= Prandtle number calculated for water properties at the reference temperature, (-)The various radiation and convection heat transfer coefficients are given by the following equations. The convection heat transfer coefficients between the covers and the absorber plate are estimated from the empirical correlation for the Nusselt number for air gap between two parallel plates developed by Hollands et al. (1976) is:

Nua=1+1.44{1−1708(sin1.8β)1.6Racosβ}{1−1708Racosβ}++{(Racosβ5830)1/3−1}+

hc=Nu⋅k/L

hrp−c2=σ(Tp+Tc2)(T2p+T2c2)1/εp+1/εc2−1

hrc1−c2=σ(Tc1+Tc2)(T2c1+T2c2)1/εc1+1/εc2−1

The long wave radiation exchange coefficients between the outer collector cover and the sky and ground referencing the ambient air temperature for mathematical simplification are given.

hrc1−s=Fsεc1σ(Tc1+Ts)(T2c1+T2s)(Tc1−Ts)(Tc1−Ta)

hrc1−g=Fgεc1σ(Tc1+Tg)(T2c1+T2g)(Tc1−Tg)(Tc1−Ta)

hcc1−a=hcc1−s+hcc1−g

The convection heat transfer coefficient from the outer cover to the surrounding air is given by:

hcc1−a=2.8+3.0Vw

When the bottom surface boundary condition is AmbientAir, the combined conductance from the outer cover to the surrounding is calculated from the equation below (Duffie and Beckman, 1991).

hcomb=5.7+3.8Vw

The overall loss coefficient through the bottom and side of the collector-storage is estimated as follows:

Ub=ULb(Ab/A)

Us=[1ULs(As/A)+1hcomb]−1

Where,

= thermal emissivity of collector cover_{c1}1, (-)= thermal emissivity of collector cover_{c2}2, (-)F= view factor from the collector to the sky, (-)_{s}F= view factor from the collector to the ground, (-)_{g}T= temperature of collector cover_{c1}1, (K)T= temperature of collector cover_{c2}2, (K)T= sky temperature, (K)_{s}T= ground temperature, (K)_{g}k= thermal conductivity of air, (W/m K)L= air gap between the covers, (m)β= inclination of the plates or covers to the horizontal, (radian)V= wind speed, (m/s)_{w}U= user specified bottom heat loss conductance, W/m_{Lb}^{2}KU= user specified side heat loss conductance, W/m_{Ls}^{2}KA= collector bottom heat transfer area, m_{b}^{2}A= collector side area, m_{s}^{2}h= combined conductance from the outer cover to the ambient air, W/m_{comb}^{2}KTransmittance-Absorptance ProductThe transmittance-absorptance product of solar collector is determined using ray tracing method for any incident angle (Duffie and Beckman, 1991). This requires optical properties of the cover and absorber materials and the the transmittance-absorptance product for any incident angle is given by:

(τα)θ=τα1−(1−α)ρd

The transmittance of the cover system for single and two cover are given by:

τ=12[(τ1⋅τ21−ρ1ρ2)⊥+(τ1⋅τ21−ρ1ρ2)∥]

ρ=12[(ρ1+τ⋅ρ2⋅τ1τ2)⊥+(ρ1+τ⋅ρ2⋅τ1τ2)∥]

The effective transmittance, reflectance and absorptance of a single cover are given by:

τ=τa2⎧⎨⎩1−r⊥1+r⊥[1−r2⊥1−(r⊥τa)2]+1−r∥1+r∥⎡⎣1−r2∥1−(r∥τa)2⎤⎦⎫⎬⎭

ρ=12⎧⎨⎩[r⊥+(1−r⊥)2τ2ar⊥1−(r⊥τa)2]+⎡⎣r∥+(1−r∥)2τ2ar∥1−(r∥τa)2⎤⎦⎫⎬⎭

α=(1−τa)2{(1−r⊥1−r⊥τa)+(1−r∥1−r∥τa)}

The transmittance of the cover system with absorption only considered

, is defined as:_{a}τa=exp(−KLcosθ2)

θ2=sin−1(sinθ1⋅n1n2)

The reflectance of un-polarized radiation on passing from medium 1 with reflective index

n_{1}to medium 2 with reflective indexn_{2}is given by:r⊥=sin2(θ2−θ1)sin2(θ2+θ1)

r∥=tan2(θ2−θ1)tan2(θ2+θ1)

The sky and ground reflected diffuse radiations average equivalent incident angles are approximated by Brandemuehl and Beckman correlation (Duffie and Beckman, 1991) as follows:

θsd=59.68−0.1388β+0.001497β2

θgd=90−0.5788β+0.002693β2

where,

τ = transmittance of the cover system, (-)

τ1 = transmittance of the cover 1, (-)

τ2 = transmittance of the cover 2, (-)

α = absorptance of the absorber plate, (-)

ρd = diffuse reflectance of the inner cover, (-)

L= thickness of a cover material, (m)K= extinction coefficient of a cover material, (m^{-1})θ1 = angle of incidence, degree

θ2 = angle of refraction, degree

r⊥ = parallel component of reflected un-polarized radiation, (-)

r∥ = perpendicular component of reflected un-polarized radiation, (-)

β = slope of the collector, degree

θsd = equivalent incident angle for sky diffuse solar radiation, degree

θgd = equivalent incident angle for ground diffuse solar radiation, degree

The integral collector storage unit thermal performance parameters are calculated as follows:

Qdelivered=˙mwCw(Tw−Twi)

QStored=mwCwdTwdt

QSkinLoss=A⋅[Ut(Tp−Tw)+Ub(Tw−TOSC)+Us(Tw−Ta)]

ηthermal=mwCwdTwdt+˙mwCw(Tw−Twi)A⋅It

## References:[LINK]

Duffie, J.A., and W.A. Beckman. 1991.

Solar Engineering of Thermal Processes,2d ed. New York: John Wiley & Sons.Kumar, R. and M.A. Rosen. Thermal performance of integrated collector storage solar water heater with corrugated absorber surface. Applied Thermal Engineering: 30 (2010) 1764-1768.

Fujii, T., and H. Imura. Natural convection heat transfer from aplate with arbitrary inclination. International Journal of Heat and Mass Transfer: 15(4), (1972), 755-764.

## Photovoltaic Thermal Flat-Plate Solar Collectors[LINK]

Photovoltaic-Thermal solar collectors (PVT) combine solar electric cells and thermal working fluid to collect both electricity and heat. Athough there are currently comparatively few commercial products, PVT research has been conducted for the past 30 years and many different types of collectors have been studied. Zondag (2008) and Charalambous et. al (2007) provide reviews of the PVT literature. Because PVT is much less commercially-mature, there are no standards or rating systems such as for thermal-only, hot-water collectors. EnergyPlus currently has one simple model based on user-defined efficiencies but a more detailed model based on first-principles and a detailed, layer-by-layer description is under development.

The PVT models reuse the PV models for electrical production. These are described elsewhere in this document in the section on Photovoltaic Arrays-Simple Model

## Simple PVT Thermal Model[LINK]

The input object SolarCollector:FlatPlate:PhotovoltaicThermal provides a simple PVT model that is provided for quick use during design or policy studies. The user simply provides values for a thermal efficiency and the incident solar heats the working fuild. The model also includes a cooling mode for air-based systems where a user-provided surface emmittance is used to model cooling of the working fluid to the night sky (water-based cooling will be made available once a chilled water storage tank is available). No other details of the PVT collector’s construction are required as input data.

The simple model can heat either air or liquid. If it heats air, then the PVT is part of HVAC air system loop with air nodes connected to an air system. If it heats liquid, then the PVT is part of plant loop with nodes connected to a plant loop and the plant operating scheme determines flows.

Air-system-based PVT modeling include a modulating bypass damper arrangement. Control logic decides if the air should bypass the collector to better meet setpoint. The model requires a drybulb temperature setpoint be placed on the outlet node. The model assume the collector is intended and available for heating when the incident solar is greater than 0.3 W/m

^{2}and otherwise it is intended for cooling. The inlet temperature is compare to the setpoint on the outlet node to determine if cooling or heating are beneficial. If they are, then the PVT thermal models are applied to condition the air stream. If they are not beneficial, then the PVT is completely bypassed and the inlet node is passed directly to the outlet node to model a completely bypassed damper arrangement. A report variable is available for bypass damper status.Plant-based PVT do not include a bypass (although one could be used in the plant loop). The collector requests its design flow rate but it otherwise relies on the larger plant loop for control.

When the PVT themal collector is controlled to be “on,” in heating mode, and working fluid is flowing, the model calculates the outlet temperature based on the inlet temperature and the collected heat using the following equations.

Qtherm=Asurf⋅factiv⋅GT⋅ηthermal

where,

Qtherm is the thermal energy collected [W]

Asurf is the net area of the surface [m

^{2}]factiv is the fraction of surface aire with active PV/T collector, and

ηthermal is the thermal conversion efficiency.

Tout=Tin+Qtherm˙mcp

where,

Tout is the temperature of the working fluid leaving the PV/T

Tin is the temperature of the working fluid entering the PV/T

˙m is the entire mass flow of the working fluid through the PV/T

cp is the specific heat of the working fluid.

For air-based systems, the value of Tout is then compared to the temperature setpoint on the outlet node. If Tout exceeds the desired outlet temperature, Tset,out, then a bypass fraction is calculated to model a modulating bypass damper using:

fbypass=(Tset,out−Tout)(Tin−Tout)

When the PVT themal collector is controlled to be “on,” in cooling mode, and working fluid is flowing, the model calculates the outlet temperature based on the inlet temperature and the heat radiated and convected to the ambient using a heat balance on the outside face of the collector:

˙mcp(Tin−Tout)=˙QLWR+˙Qconv

Where,

˙QLWR is the net rate of long wavelength (thermal) radiation exchange with the air, night sky, and ground. See the section “External Longwave Radiation” in the Outside Surface Heat Balance, for full discussion of how this modeled in EnergyPlus using linearized radiation coefficients.

˙Qconv is the net rate of convective flux exchange with outdoor air. See the section “Exterior/External Convection” in the Outside Surface Heat Balance, for full discussion of how this modeled in EnergyPlus. The surface roughness is assumed to be “very smooth.”

The simple model assumes that the effective collector temperature, Tcol, is the average of the working fluid inlet and outlet temperatures so that we can make the following substitution:

Tout=2Tcol−Tin

Substituting and solving for Tcol we obtain the following model for collector temperatures during a (possible) cooling process:

Tcol=2˙mcpTin+Asurffactiv(hr,gndTgnd+hr,skyTsky+hr,airTair+hc,extTair)2˙mcp+Asurffactiv(hr,gnd+hr,sky+hr,air+hc,ext)

Then the outlet temperature can be calculated and heat losses determined. However, the model allows only sensible cooling of the air stream and limits the outlet temperature to not go below the dewpoint temperature of the inlet.

PVT collectors have a design volume flow rate for the working fluid that is autosizable. For air-based systems used as pre-conditioners, the volume flow rate is sized to meet the maximum outdoor air flow rate. For water-based systems on the supply side of a plant loop, each of the PVT collectors are sized to the overall loop flow rate. For water-based systems on the demand side of a plant loop, the collectors are sized using a rule-of-thumb for typical flow rates per unit of collector area. This rule-of-thumb is based on a constant factor of 1.905x10

^{-5}m^{3}/s-m^{2}^^that was developed by analyzing SRCC data set for conventional solar collectors (see data set SolarCollectors.idf) and averaging the ratio for all 171 different collectors.## References[LINK]

Charalambous, P.G., Maidment, G.G., Kalagirou, S.A., and Yiakoumetti, K., Photovoltaic thermal (PV/T) collectors: A review. Applied Thermal Engineering 27 (2007) 275-286.

Zondag, H.A. 2008. Flat-plate PV-Thermal collectors and systems: A review. Renewable and Sustainable Energy Reviews 12 (2008) 891-959.

## Unglazed Transpired Solar Collectors[LINK]

The input object SolarCollector:UnglazedTranspired provides a model for transpired collectors that are perhaps one of the most efficient ways to collect solar energy with demonstrated instantaneous efficiencies of over 90% and average efficiencies of over 70%. They are used for preheating outdoor air needed for ventilation and processes such as crop drying.

In EnergyPlus, an unglazed transpired solar collector (UTSC) is modeled as a special component attached to the outside face of a heat transfer surface that is also connected to the outdoor air path. A UTSC affects both the thermal envelope and the HVAC air system. From the air system’s point of view, a UTSC is heat exchanger and the modeling needs to determine how much the device raises the temperature of the outdoor air. From the thermal envelope’s point of view, the presence of the collector on the outside of the surface modifies the conditions experienced by the underlying heat transfer surfaces. EnergyPlus models building performance throughout the year and the UTSC will often be “off” in terms of forced airflow, but the collector is still present. When the UTSC is “on” there is suction airflow that is assumed to be uniform across the face. When the UTSC is “off” the collector acts as a radiation and convection baffle situated between the exterior environment and the outside face of the underlying heat transfer surface. We distinguish these two modes of operation as

activeorpassiveand model the UTSC component differently depending on which of these modes it is in.## Heat Exchanger Effectiveness[LINK]

The perforated absorber plate is treated as a heat exchanger and modeled using a traditional effectiveness formulation. The heat exchanger effectiveness, εHX , is determined from correlations derived from small-scale experiments. Two correlations available in the literature are implemented in EnergyPlus. The first is based on research by Kutscher at the National Renewable Energy Laboratory. The second is based on the research by Van Decker, Hollands, and Brunger at the University of Waterloo. Because both correlations are considered valid, the choice of which correlation to use is left to the user.

## Kutscher Correlation[LINK]

Kutscher’s (1994) correlation encompasses surface convection between the collector and the incoming outdoor air stream that occurs on the front face, in the holes, and along the back face of the collector. The correlation uses a Reynolds number based on the hole diameter as a length scale and the mean velocity of air as it passes through the holes as the velocity scale:

ReD=VhDν

where,

Vh is the velocity through the holes [m/s]

D is the hole diameter [m]

ν is the kinematic viscosity of air [m

^{2}/s]The correlation is a function of Reynolds number, hole geometry, the free stream air velocity, and velocity through the holes:

NuD=2.75[(PD)−1.2Re0.43D+0.011σReD(U∞Vh)0.48]

where,

P is the pitch, or distance between holes, [m],

D is the diameter of the holes, [m],

σ is the porosity, or area fraction of the holes, [dimensionless],

Vh is the mean velocity of air passing through the holes, [m/s],

U∞ is the free stream velocity (local wind speed) [m/s].

The Nusselt number is formulated as:

NuD=UDk

where,

U is the overall heat transfer coefficient based on log mean temperature difference, [W/m

^{2}·K], andk is the thermal conductivity of air [W/m·K].

The heat exchanger effectiveness is:

εHX=1−e[−UA˙mcp]

Kutscher’s relation was formulated for triangular hole layout, but based on Van Decker et al. (2001) we allow using the correlation for square hole layout and scale P by a factor of 1.6.

## Van Decker, Hollands, and Brunger Correlation[LINK]

Van Decker et. al. extended Kutscher’s measurements to include a wider range of collector parameters including plate thickness, pitch, suction velocities, and square hole patterns. Their model formulation differs from Kutscher’s in that the model was built up from separate effectiveness models for the front, back, and holes of the collector. Their published correlation is:

εHX=[1−(1+ResMax(1.733Re−1/2w,0.02136)−1)]×[1−(1+0.2273Re1/2b)−1]×e(−0.01895PD−20.62ReDtD)

where,

Res=VsPv

Rew=U∞Pv

Reb=VhPv

Vs is the average suction velocity across the front face of the collector [m/s]

t is the collector plate thickness

## Heat Exchanger Leaving Temperature[LINK]

Using either of the correlations above allows determining the heat exchanger effectiveness from known values. By definition the heat exchanger effectiveness is also:

εHX=Ta,HX−TambTs,coll−Tamb

where,

Ta,HX is the temperature of the air leaving the collector and entering the plenum [ºC]

Ts,coll is the temperature of the collector’s absorber plate, [ºC], and

Tamb is the temperature of the ambient outdoor air [ºC].

By rewriting equation to solve for Ta,HX we see that the temperature of the heated outdoor air entering the plenum can be determined once the collector surface temperature is known,

Ta,HX=εHXTs,coll+(1−εHX)Tamb

## Collector Heat Balance[LINK]

The collector is assumed to be sufficiently thin and high-conductivity so that it can be modeled using a single temperature (for both sides and along its area). This temperature Ts,coll is determined by formulating a heat balance on a control volume that just encapsulates the collector surface. The heat balances are formulated separately for active and passive modes and are diagrammed in the following figure.

Observe that for the passive case, we do not use the heat exchanger relations to directly model the interaction of ventilating air with the collector. This is because these relations are considered to not apply when the UTSC is in passive mode. They were developed for uni-directional flow (rather than the balanced-in-and-out flow expected from natural forces) and for specific ranges of suction face velocity. Therefore, this heat transfer mechanism is handled using classical surface convection models (as if the collector was not perforated). (Air exchanges are modeled as ventilation in the plenum air heat balance but do not interact with the hole edges in the collector surface.)

Transpired Collector Heat Balance

When the UTSC is active, the heat balance on the collector surface control volume is:

q′′αsol+q′′LWR,Env+q′′conv,wind−q′′HX+q′′LWR,plen+q′′source=0

where:

q′′αsol is absorbed direct and diffuse solar (short wavelength) radiation heat flux.

q′′LWR,Env is net long wavelength (thermal) radiation flux exchange with the air and surroundings.

q′′conv,wind is surface convection flux exchange with outdoor air under high wind and rain conditions. Note that this term is usually assumed to be zero in UTSC model development but we add the term to allow for deteriorated performance of the UTSC under poor conditions.

q′′HX is heat exchanger flux from collector to incoming outdoor air.

q′′LWR,plen is net long wavelength (thermal) radiation flux exchange with the outside face of the underlying surface(s).

q′′source is a source/sink term that accounts for energy exported out of the control volume when the collecter’s absorber plate is a hybrid device such as a photovoltaic panel.

While the heat balance on the passive collector surface control volume is:

q′′αsol+q′′LWR,Env+q′′conv,Env+q′′LWR,plen+q′′conv,plen+q′′source=0

where:

q′′conv,Env = surface convection flux exchange with outdoor air.

q′′conv,plen = surface convection flux exchange with plenum air.

All terms are positive for net flux to the collector except the heat exchanger term, which is taken to be positive in the direction from collector to incoming air stream. Each of these heat balance components is introduced briefly below.

## External SW Radiation[LINK]

q′′αsol is calculated using procedures presented elsewhere in this manual and includes both direct and diffuse incident solar radiation absorbed by the surface face. This is influenced by location, surface facing angle and tilt, shading surfaces, surface face material properties, weather conditions, etc.

## External LW Radiation[LINK]

q′′LWR,Env is a standard radiation exchange formulation between the surface, the sky, the ground, and the atmosphere. The radiation heat flux is calculated from the surface absorptivity, surface temperature, sky, air, and ground temperatures, and sky and ground view factors. Radiation is modeled using linearized coefficients.

## External Convection[LINK]

q′′conv,Env is modeled using the classical formulation: q′′conv = h

_{co}(T_{air}- T_{o}) where h_{co}, is the convection coefficient. This coefficient will differ depending on whether or not the UTSC is active or passive. When the UTSC is passive, h_{co}is treated in the same way as an outside face with ExteriorEnvironment conditions. When the UTSC is active, the special suction airflow situation of a transpired collector during operation means that h_{co}is often zero because the suction situation can eliminate mass transport away from the collector. However when the winds are high, the strong turbulence and highly varying pressures can cause the suction flow situation to breakdown. Therefore, we include the q′′conv,wind term in the heat balance and use a special coefficient hc,wind to model this lost heat transfer. In addition, when it is raining outside, we assume the collector gets wet and model the enhanced surface heat transfer using a large value for hc,wind.## Heat Exchanger[LINK]

q′′HX is modeled using the classical formulation q′′HX=˙mcp(Ta,HX−Tamb)A where Ta,HX is determined using correlations described above. When the UTSC is active, the air mass flow is determined from the operation of the outdoor air mixer component. When the UTSC is off, this term is zero.

## Plenum LW Radation[LINK]

q′′LWR,plen is a standard radiation exchange formulation between the collector surface and the underlying heat transfer surface located across the plenum. Radiation is modeled using linearized coefficients.

## Plenum Convection[LINK]

q′′conv,plen is modeled using the classical formulation: q′′conv = h

_{cp}(T_{air}- T_{o}) where h_{cp}, is the convection coefficient. This coefficient is taken as zero when the UTSC is operating because of the suction airflow situation. When the UTSC is off, the value for h_{cp}is obtained from correlations used for window gaps from ISO (2003) standard 15099.Substituting models into and solving for Ts,coll yields the following equation when the UTSC is active (“on”):

Ts,coll=(Isα+hr,atmTamb+hr,skyTsky+hr,gndTamb+hr,plenTso+hc,windTamb+˙mcpATamb−˙mcpA(1−εHX)Tamb+q′′source)(hr,atm+hr,sky+hr,gnd+hr,plen+hc,wind+˙mcpAεHX)

and substituting into yields the following equation when the UTSC is passive (“off”):

Ts,coll=(Isα+hcoTamb+hr,atmTamb+hr,skyTsky+hr,gndTamb+hr,plenTso+hc,plenTa,plen+q′′source)(hco+hr,air+hr,sky+hr,gnd+hr,plen+hc,plen)

where,

Is is the incident solar radiation of all types [W/m

^{2}],α is the solar absorptivity of the collector [dimensionless],

hr,atm is the linearized radiation coefficient for the surrounding atmosphere [W/m

^{2}·K],Tamb is the outdoor drybulb from the weather data, also assumed for ground surface [ºC],

hr,sky is the linearized radiation coefficient for the sky [W/m

^{2}·K],Tsky is the effective sky temperature [ºC],

hr,gnd is the linearized radiation coefficient for the ground [W/m

^{2}·K],hr,plen is the linearized radiation coefficient for the underlying surface [W/m

^{2}·K],hc,wind is the convection coefficient for the outdoor environment when the UTSC is active and winds are high or it is raining [W/m

^{2}·K],Tso is the temperature of the outside face of the underlying heat transfer surface [ºC],

˙m is the air mass flow rate when in active mode [kg/s],

cp is the specific heat of air at constant pressure [J/kg·K],

A is the area of the collector [m

^{2}],hco is the convection coefficient for the outdoor environment [W/m

^{2}·K],hc,plen is the convection coefficient for the surfaces facing the plenum [W/m

^{2}·K], andTa,plen is the air drybulb for air in the plenum and entering the outdoor air system [ºC].

## Plenum Heat Balance[LINK]

The

plenumis the volume of air located between the collector and the underlying heat transfer surface. The plenum air is modeled as well-mixed. The uniform temperature of the plenum air, Ta,plen, is determined by formulating a heat balance on a control volume of air as diagrammed below.Note that we have formulated the control volumes with slight differences for the active and passive cases. For the active case, the suction air situation and heat exchanger effectiveness formulations dictate that the collector surface control volume encompass part of the air adjacent to both the front and back surfaces of the collector. However for the passive case, the collector surface control volume has no air in it and the plenum air control volume extends all the way to the surface of the collector.

Transpired Collector Plenum Air Heat Balance

When the UTSC is active, the heat balance on the plenum air control volume is:

˙Qair+˙Qco=0

where,

˙Qair is the net rate of energy added by suction air convecting through the control volume.

˙Qco is the net rate of energy added by surface convection heat transfer with the underlying surface.

When the UTSC is passive, the heat balance on the plenum air control volume is:

˙Qvent+˙Qco+˙Qc,coll=0

where,

˙Qvent is the net rate of energy added from infiltration - where outdoor ambient air exchanges with the plenum air.

˙Qc,coll is the net rate of energy added by surface convection heat transfer with the collector.

Substituting into and solving for Ta,plen yields the following equation for when the UTSC is active:

Ta,plen=(˙mcpTa,HX+hc,plenATso)(˙mcp+hc,plenA)

And substituting into yields the following equation when the UTSC is passive:

Ta,plen=(hc,plenATso+˙mventcpTamb+hc,plenATs,coll)(hc,plenA+˙mventcp+hc,plenA)

where,

˙mvent is the air mass flow from natural forces [kg/s]

The literature on UTSC does not appear to address the passive mode of operation and no models for ˙mvent have been identified. Nevertheless, natural buoyancy and wind forces are expected to drive air exchange between the plenum and ambient and some method of modeling ˙mvent is needed. Reasoning that the configuration is similar to single-side natural ventilation, we elect to use correlations for natural ventilation presented in Chapter 26. of ASHRAE HOF (2001).

˙mvent=ρ−˙Vtot

where,

ρ is the density of air [kg/m

^{3}], and−˙Vtot=−˙Vwind+−˙Vthermal is the total volumetric flow rate of air ventilating in and out of the plenum.

−˙Vwind=CvAinU∞

−˙Vthermal=CDAin√2gΔHNPL(Ta,plen−Tamb)/Ta,plen (if Ta,plen>Tamb)

−˙Vthermal=CDAin√2gΔHNPL(Tamb−Ta,plen)/Tamb (if Tamb>Ta,plen and UTSC vertical)

Cv is the effectiveness of the openings that depends on opening geometry and the orientation with respect to the wind. ASHRAE HoF (2001) indicates values ranging from 0.25 to 0.6. In the UTSC model, this value is available for user input and defaulted to 0.25.

CD is the discharge coefficient for the opening and depends on opening geometry. In the UTSC model, this value is available for user input and defaulted to 0.65.

Mass continuity arguments lead to modeling the area of the openings as one half of the total area of the holes, so we have:

Ain=Aσ2

g is the gravitational constant taken as 9.81 [m/s

^{2}].ΔHNPL is the height from midpoint of lower opening to the Neutral Pressure Level. This is taken as one-fourth the overall height of the UTSC if it is mounted vertically. For tilted collectors, the nominal height is modified by the sine of the tilt. If the UTSC is mounted horizontally (e.g. on the roof) then the ΔHNPL is taken as the gap thickness of the plenum.

If the UTSC is horizontal and Tamb>Ta,plen then −˙Vthermal=0 because this is a stable situation.

## Underlying Heat Transfer Surface[LINK]

The UTSC is applied to the outside of a heat transfer surface. This surface is modeled using the usual EnergyPlus methods for handling heat capacity and transients - typically the CTF method. These native EnergyPlus Heat Balance routines are used to calculate Tso. The UTSC model is coupled to the underlying surface using the OtherSideConditionsModel mechanism. The UTSC model provides values for hr,plen, Ts,coll, hc,plen, and Ta,plen for use with the Heat Balance Model calculations for the outside face of the underlying surface (described elsewhere in this manual).

## Solar and Shading Calculations[LINK]

The transpired collector object uses a standard EnergyPlus surface in order to take advantage of the detailed solar and shading calculations. Solar radiation incident on the surface includes beam and diffuse radiation, as well as radiation reflected from the ground and adjacent surfaces. Shading of the collector by other surfaces, such as nearby buildings or trees, is also taken into account.

## Local Wind Speed Calculations[LINK]

The outdoor wind speed affects terms used in modeling UTSC components. The wind speed in the weather file is assumed to be measured at a meteorological station located in an open field at a height of 10 m. To adjust for different terrain at the building site and differences in the height of building surfaces, the local wind speed is calculated for each surface.

The wind speed is modified from the measured meteorological wind speed by the equation (ASHRAE 2001):

U∞=Vmet(δmetzmet)amet(zδ)a

where z is the height of the centroid of the UTSC, z

_{met}is the height of the standard metereological wind speed measurement, and a and δ are terrain-dependent coefficients. δ is the boundary layer thickness for the given terrain type. The values of a and δ are shown in the following tables:TerrainDescriptionExponent, aLayer Thickness, δ (m)The UTSC can be defined such that it has multiple underlying heat transfer surfaces. The centroid heights for each surface are area-weighted to determine the average height for use in the local wind calculation.

## Convection Coefficients[LINK]

UTSC modeling requires calculating up to three different coefficients for surface convection heat transfer. These coefficients are defined in the classic way by:

hc=Tair−Tsurfq′′conv

First, hco is the convection coefficient for the collector surface facing the outdoors when the UTSC is passive. It is modeled in exactly the same way as elsewhere in EnergyPlus and will depend on the user setting for Outside Convection Algorithm - Outside Surface Heat Balance entry elsewhere in this document.

Second, hc,plen is the convection coefficient for surfaces facing the plenum. This coefficient is applied to just the underlying surface’s convection when the UTSC is active and to both the collector and the underlying surface when the UTSC is passive. When the UTSC is active, we use the convection correlation for forced air developed by McAdams (1954) as published by ASHRAE HoF (2001):

hc,plen=5.62+3.9Vp

where,

Vp is the mean velocity in the plenum determined from Vp=˙m2ρAp where Ap is the effective cross section area of the plenum perpendicular to the primary flow direction. When the UTSC is passive, we model the convection in the same way used in EnergyPlus to model air gaps in windows. These correlations vary by Rayleigh number and surface tilt and are based on the work of various research including Hollands et. al., Elsherbiny et. al., Wright, and Arnold. The formulations are documented in ISO (2003) standard 15099. For the UTSC implementation, the routines were adapted from Subroutine NusseltNumber in WindowManager.f90 (by F. Winkelmann), which itself was derived from Window5 subroutine “nusselt”.

Third, hc,wind is the convection coefficient used to degrade the UTSC performance under environmental conditions with high wind or rain. If the weather file indicates it is raining, then we set hc,wind = 1000.0 which has the effect of making the collector the same temperature as ambient air. The heat exchanger effectiveness correlations described above account for a moderate amount of wind, but the correlations appear confined to the range 0 to 5.0 m/s. Therefore we set hc,wind equal to zero if U∞ is <= 5.0 m/s. If U∞ is > 5.0 m/s then we use the McAdams correlation but with a reduced velocity magnitude:

$$ hc,wind=5.62+3.9(U∞−5.0)

## Radiation Coefficients[LINK]

UTSC modeling requires calculating up to four different linearized coefficients for radiation heat transfer. Whereas radiation calculations usually use temperature raised to the fourth power, this greatly complicates solving heat balance equations for a single temperature. Linearized radiation coefficients have the same units and are used in the same manner as surface convection coefficients and introduce very little error for the temperature levels involved.

The radiation coefficient, hr,plen , is used to model thermal radiation between the collector surface and the outside face of the underlying heat transfer surface. We assume a view factor of unity. It is calculated using:

hr,plen=σSBecolleso(T4s,coll−T4so)(Ts,coll−Tso)

where,

all temperatures are converted to Kelvin,

σSB is the Stefan-Boltzmann constant,

ecoll is the longwave thermal emmittance of the collector, and

eso is the longwave the