This function corrects the surface temperature climatic constants by applying n-factors with the methodology prescribed in the District Cooling Guide (ASHRAE, 2013).
n-factors corresponds to the ratio of freezing (/thawing)
degree days or index between the air and the ground surface and can
be used to couple the air and ground surface temperatures.
For example, a freezing n-factor of 1.35 means that during
the freezing season, the daily average ground surface temperature
is on average 1.35 times colder than the air (using the freezing
temperature of water as a reference).
In its guide, ASHRAE suggests to first apply the n-factors
to the uncorrected ground temperature at zero burial depth, and to
fit a new sinusoidal curve to the result.
In this function, the freezing n-factor is applied to days
where the daily mean air temperature is below the freezing
temperature (0degC), whereas the thawing n-factor is applied
to the remaining days.
Tground -
Tfreezing = ηfreezing * (Tair -
Tfreezing) if Tair <=
Tfreezing
Tground - Tfreezing = ηthawing *
(Tair - Tfreezing) if Tair
> Tfreezing
The sinusoidal curve is then fitted using the analytical solution proposed in Appendix B of the District Heating Guide.
n-factors are specific to a climate and surface cover, and should be extrapolated from other sites with caution. As a first approximation, tabulated values of n-factors are available in Lunardini (1981) and Freitag and McFadden (1997).
ASHRAE (2013). District Cooling Guide. ASHRAE, Atlanta,
GA.
ASHRAE (2013). District Heating Guide. ASHRAE, Atlanta,
GA.
D.W. Riseborough (2003). Thawing and freezing indices in the
active layer. Proceedings of the 8th International Conference
on Permafrost, Swets & Zeitlinger.
V.J. Lunardini (1981). Heat Transfer in Cold Climates. Van
Nostrand Reinhold Company.
D.R. Freitag and T. McFadden (1997). Introduction to Cold
Regions Engineering. American Society of Civil Engineers.
function surfaceTemperature extends Modelica.Icons.Function; input ClimaticConstants.Generic cliCon "Surface temperature climatic conditions"; input Real nFacTha "Thawing n-factor (TAir > 0degC)"; input Real nFacFre "Freezing n-factor (TAir <= 0degC)"; output ClimaticConstants.Generic corCliCon "Corrected surface temperature climatic conditions"; end surfaceTemperature;
parameter declaration on a protected
variable which causes an error in Dymola 2025x beta1.