modelSolarCollector_L1_constProp

Solar flat plate collector model (EN 12975) with effective heat capacity for transient behavior, constant properties

Extends from TransiEnt.Basics.Icons.SolarThermalCollector.

Information

1. Purpose of model

A simple solar collector providing useful energy gain as recommended by EN 12975 steady state thermal performance equation [1]. Effective heat capacity was added for transient behavior according to [2].

2. Level of detail, physical effects considered, and physical insight

No physical but parameter based model. Model considers heat loss to ambience but no convection, and pressure losses (linear or quadratic).

No TILMedia models used because they produce warnings at temperatures below 0°C even if no heat flow is generated.

3. Limits of validity

Without effective heat capacity only valid for steady state performance.

Model ignores wind speed.

4. Interfaces

waterIn: water inlet

waterOut: water outlet

T_in: output of temperature at inlet

T_out: output of temperature at outlet

G: output of total irradiance

irradiance_direct_measured_input: input of direct irradiance

irradiance_diffuse_horizontal_input: input of diffuse horizontal irradiance

5. Nomenclature

(no remarks)

6. Governing Equations



  Q= area*((G_total*eta_0)-a1*(T_m-T_amb)-a2*(T_m-T_amb)^2)

G_total = iam*direct_irradiance.irradiance+IAM_diffuse*(diffuse_irradiance.irradiance+ground_reflected_irradiance.irradiance)

    iam= if biaxial then iam_obj.value_longitudinal*iam_obj.value_transversal
  else
   iam_obj.value

 T_m=0.5*(T_in+T_out)



  waterIn.m_flow + waterOut.m_flow = 0

  waterOut.p =  if noFriction then waterIn.p+gravAcc*fluidIn.d*(z1-z2) else waterIn.p+gravAcc*fluidIn.d*(z1-z2)-(a*waterOut.m_flow+b*waterOut.m_flow^2)


h=c*T+h_0

7. Remarks for Usage

This solar collector can be used for temperatures below 0°C as well.

8. Validation

not validated yet but tested in the check model "TransiEnt.Producer.Heat.SolarThermal.Check.TestCollectorFluidCycle_constProp"

9. References

[1] P. Kovacs, “Quality Assurance in solar thermal heating and cooling technology – keeping track with recent and upcoming developments. A guide to the standard EN 12975,” Brussels, 2012.

[2] T. Osório and M. J. Carvalho, “Testing of solar thermal collectors under transient conditions,” Sol. Energy, vol. 104, pp. 71–81, 2014.

10. Version History

Model created by Tobias Toerber (tobias.toerber@tuhh.de), Jul 2015

Model modified by Sascha Guddusch (sascha.guddusch@tuhh.de), May 2016

Model modified by Anne Senkel (anne.senkel@tuhh.de), Mar 2017

Model modified by Lisa Andresen (andresen@tuhh.de), Apr. 2017

Model modified by Oliver Schülting (oliver.schuelting@tuhh.de), May 2018: added possibility to use inputs for solar irradiation

Model modified by Oliver Schülting (oliver.schuelting@tuhh.de), Aug 2018: added boolean 'UseStationaryCalculationMethod' to deactive instationary part of heat flow calculation '-c_eff*der_T'. Results are nearly the same but unrealistic peaks can be avoided if calculation is stationary

Model modified by Carsten Bode (c.bode@tuhh.de), Nov 2018: replaced fluid objects by equations, eliminated m_flowInv and x

Model modified by Anne Hagemeier (anne.hagemeier@umsicht.fraunhofer.de), Aug 2021: enabled the model to be useable without fluid ports

Parameters

TypeNameDefaultDescription
SI.Densityd985.269Density of medium
SI.SpecificHeatCapacityc4190.96Specific heat capacity of medium
SI.SpecificEnthalpyh_059.6504Specific enthalpy of medium at 0°C
General
TILMedia.VLEFluidTypes.BaseVLEFluidmediumsimCenter.fluid1
BooleanuseHomotopysimCenter.useHomotopytrue = homotopy method is used during initialisation
General › General
SI.HeatFlowRateQ_flow_nNominal heat flow rate (for cost calculation)
SI.AreaareaAperture area
Realc_effEffective thermal capacity of the collector
SI.IrradianceG_min0Minimum Irradiance before collector is working
BooleanUseStationaryCalculationMethodtrueChoose if calculation uses only stationary equations or not
BooleanuseFluidPortstrueTrue if fluid ports shall be used
Coefficients for thermal performance
Realeta_0Zero-loss collector efficiency
Reala1Heat loss coefficient at (T_m - T_amb) = 0
Reala2Temperature dependent heat loss coefficient
Pressure drop
BooleannoFrictiontruetrue = assume no pressure loss due to friction
Integern_serial1Number of collectors in series (max. 12)
Reala0Linear pressure drop coefficient
Realb0Quadratic pressure drop coefficient
SI.Heightz10Height inlet
SI.Heightz20Height outlet
Irradiance › Losses
RealSoiling0Average annual losses of radiation in % due to soiling
Irradiance › Solartime
SI.Anglelongitude_localModelica.Units.Conversions.from_deg(10)longitude of the local position, east positive, 10 East for Hamburg
SI.Anglelongitude_standardModelica.Units.Conversions.from_deg(15)needed for calculation of coordinated universal time (utc), 15 for central european time, 30 for central european summer time
Irradiance › Extraterrestrial Irradiance
SI.AnglelatitudeModelica.Units.Conversions.from_deg(53.55)latitude of the local position, north posiive, 53,55 North for Hamburg
SI.AngleslopeModelica.Units.Conversions.from_deg(53.55)slope of the tilted surface, assumption
SI.AnglesurfaceAzimuthAngle0surface azimuth angle
Irradiance › Skymodel
Realreflectance_ground0.2reflectance of the ground
Booleandirect_normaltrueIs the direct irradiance measured on a surface normal to irradiance?
Booleanuse_input_datafalsechoose if input data is given by inputs - if not, simCenter data is used
IAM › General
Integerkind1IAM for direct Irradiance
Realconstant_iam_dir1constant IAM for direct irradiation
Realconstant_iam_diff1constant IAM for diffuse irradiation
Realconstant_iam_ground1constant IAM for ground-reflected irradiation
Realb01assumption: constant b0-value for IAM=1-b0*(1/cos(theta)-1)
Realiam_SRCC{1, 1, 1, 1, 1, 1, 1, 1}IAM for theta = 0, 10, 20, ..., 70
Modelica.Units.NonSI.Angle_degtheta{0, 10, 20, 30, 40, 50, 60, 70}

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.Thermal.FluidPortOutwaterPortOut
TransiEnt.Basics.Interfaces.Thermal.FluidPortInwaterPortIn
Basics.Interfaces.General.TemperatureOutT_in_int
TransiEnt.Basics.Interfaces.General.TemperatureOutT_out
Modelica.Blocks.Interfaces.RealOutputGG_total
Modelica.Blocks.Interfaces.RealOutputQ_flow_out
Modelica.Blocks.Interfaces.RealOutputm_flow_int
TransiEnt.Basics.Interfaces.Ambient.IrradianceInirradiance_direct_measured_input
TransiEnt.Basics.Interfaces.Ambient.IrradianceInirradiance_diffuse_horizontal_input
Modelica.Blocks.Interfaces.RealInputm_flow
Basics.Interfaces.General.TemperatureInT_inflow
Basics.Interfaces.General.TemperatureOutT_in

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter
TransiEnt.ModelStatisticsmodelStatistics
TransiEnt.Producer.Heat.SolarThermal.Base.IrradianceOnATiltedSurfaceirradiance
Modelica.Units.NonSI.Time_daytotaldays365total days of the year, standard=365, leap year=366
SI.TemperatureT_ambsimCenter.T_amb_var + 273.15Ambient temperature in K
SI.TemperatureT_mAverage temperature of the medium
SI.TemperatureSlopeder_TDerivative of T_m
SI.IrradianceG_totalTotal irradiance
SI.HeatFlowRateQ_flow_collectorHeat flow provided by solar collector
RealetaEfficiency factor
Modelica.Units.SI.SpecificEnthalpyh_outEnthalpy of the medium flowing out of the collector
Modelica.Units.SI.SpecificEnthalpyh_inEnthalpy of the medium flowing into the collector
TransiEnt.Producer.Heat.SolarThermal.Base.IAMIAM
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectHeatingPowercollectHeatingPower
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectCostsGeneralcollectCosts
TransiEnt.Components.Statistics.Collectors.LocalCollectors.CollectGwpEmissionsHeatcollectGwpEmissions
Base.HeatFlow_SolarThermalheatFlow2fluid
ClaRa.Components.Sensors.SensorVLE_L1_m_flowmassFlowSensorVLE

Contents

NameDescription
Skymodelchoose between HDKR and isotropic sky model
CostRecordSolarThermal