This component computes the solar heat gain of the solar thermal collector. It only calculates the solar heat gain without considering the heat loss to the environment. This model uses ratings data according to EN12975. The solar heat gain is calculated using Equations 18.298 and 18.302 in the referenced EnergyPlus documentation. The calculation is modified to use coefficients from EN12975.
The equation used to calculate solar gain is a modified version of Eq 18.302 from the EnergyPlus documentation. It is
Qflow[i] = Ac/nSeg η0 (K (τα),dir Gdir (1-shaCoe)+Kdif G dif),
where Qflow[i] is the heat gained in each segment,
Ac is the area of the collector,
nSeg
is the number of segments in the collector,
η0 is the maximum efficiency of the collector,
K(τα),dir is the incidence angle modifier
for (direct) beam radiation,
Gdir is the current beam radiation on the collector,
shaCoe
is the shading coefficient,
Kdif is the incidence angle modifier
for diffuse radiation,
and Gdif is the diffuse radiation striking the surface.
The solar radiation equation indicates that the collector is divided into
multiple segments.
The number of segments used in the simulation is specified by the user via nSeg
.
The area of an individual segment is identified by dividing the collector area
by the total number of segments.
The term shaCoe
is used to define the percentage of the collector
that is shaded.
The main difference between this model and the ASHRAE model is the handling of diffuse radiation. The ASHRAE model contains calculated incidence angle modifiers for both sky and ground diffuse radiation, while this model uses a coefficient from test data for diffuse radiation.
The incidence angle modifier is calculated using Eq 18.298 from the EnergyPlus documentation. It is
K(τα),x=1+b0 (1/cos(θ)-1)+b1 (1/cos(θ)-1)2
where K(τα),dir is the incidence angle modifier for beam radiation, b0 is the first incidence angle modifier coefficient, b1 is the second incidence angle modifier coefficient, and θ is the incidence angle.
This model reduces the heat gain rate to 0 W when the fluid temperature is within 1 degree C of the maximum temperature of the medium model. The calculation is performed using the AixLib.Utilities.Math.Functions.smoothHeaviside function.
CEN 2022, European Standard 12975:2022, European Committee for Standardization
EnergyPlus 23.2.0 Engineering Reference
Name | Description |
---|---|
Medium | Medium in the system |
each
.
displayUnit
.Temperature
to ThermodynamicTemperature
to avoid an error in the pedantic model check in Dymola 2017 FD01 beta2.Medium.T_max
to make it more efficient.