UTILIZATION OF MULTIPLE SHORT HEAT CONSUMPTION MEASUREMENT METHOD FOR DIAGNOSTICS OF HEAT SOURCE OF EXISTING BUILDING: PART 2: PRACTICAL USE OF THE METHOD

## Publications

/ Export Citation / / / Text size:

#### Architecture, Civil Engineering, Environment

Silesian University of Technology

Subject: Architecture , Civil Engineering , Engineering, Environmental

ISSN: 1899-0142

6
25
Visit(s)
0
Comment(s)
0
Share(s)

SEARCH WITHIN CONTENT

FIND ARTICLE

Volume / Issue / page

Archive
Volume 13 (2020)
Volume 12 (2019)
Volume 11 (2018)
Volume 10 (2017)
Volume 9 (2016)
Related articles

VOLUME 11 , ISSUE 2 (June 2018) > List of articles

### UTILIZATION OF MULTIPLE SHORT HEAT CONSUMPTION MEASUREMENT METHOD FOR DIAGNOSTICS OF HEAT SOURCE OF EXISTING BUILDING: PART 2: PRACTICAL USE OF THE METHOD

Citation Information : Architecture, Civil Engineering, Environment. Volume 11, Issue 2, Pages 99-106, DOI: https://doi.org/10.21307/ACEE-2018-027

Received Date : 07-December-2017 / Accepted: 24-January-2018 / Published Online: 04-April-2019

### ARTICLE

#### ABSTRACT

Streszczenie

W artykule zaprezentowano określanie wymaganej mocy cieplnej dla diagnostyki źródła ciepła istniejącego budynku mieszkalnego. W celu zaproponowano wyznaczanie wymaganej mocy cieplnej budynku z wykorzystaniem oryginalnej metody kilkukrotnego krótkiego pomiaru zużycia ciepła w istniejącym budynku mieszkalnym z wentylacją naturalną.

## 1. INTRODUCTION

This section presents the practical calculation of the number of ventilation air exchanges using the short measurement method for an existing multi-family building. The analyzes were carried out using archival materials and the analysis covered the days for which complete measurement data was required.

Nomenclature in the formulas is presented in the first part.

## 2. OBJECT DESCRIPTION

The object for the analysis is a thirty-year multi-family building. Building located in 3 climatic zone, with basement. The total number of dwellings in this building is 30, and the number of residents – 82 people. Exact details of the building being analyzed are provided in [1].

The building is supplied with heat from the district heating network through a double-function, exchanging, compact heat distribution node located in the cellar of the building. The heat node is a two-stage switch adapted for central heating and hot water supply. In heated rooms there are no heating cost allocators or other devices that allow you to determine the amount of heat consumed for heating the rooms.

## 3. DETERMINATION OF THE VALUES FOR THE VERIFICATION OF THE METHOD.

In order to verify the method there were determined heat loss coefficients for the building $HTVrz$, $HTrz$, $HVrz$, and the number of the ventilation air changes nV according to the equations from part 1. The average daily flow of internal heat gains ϕZ was also determined. Appropriate balance equations were obtained on the basis of the measurement data for 10 days moderately varying in daily temperature, selected from the group of 20 measurement days, mostly in February 2012. The selection of 20 days was the result of available and complete measurement data for these days. Therefore optimum days (days with a significant difference in midday temperature) are not selected for the analysis of the short-term in-situ measurement method, but the ones that may frequently occur in the practical application of the method (precise weather conditions cannot be predicted on measurement days for the method). The external climate data (air temperature, relative humidity of the air, wind velocity and direction, total and dispersed solar radiation intensity on the horizontal plane, temperature of the horizon) were received from the measurements of the meteorological station that lies 10 km away from the building. The consumption of the heat supplied to the installation, the flow and return temperatures of the mains and the installation part resulted from the registration of these quantities made by means of the heat meter mounted in the node. The building’s heat gain from solar radiation through the windows was determined using the SOLAR program [2]. After a preliminary analysis of the results of the check, it was assumed that the average daily heat gains flux from solar radiation contains 3/4 of the gains resulting from the radiation transmitted through the windows in the examined day and 1/4 of gains from radiation from the previous day. The estimation of an internal profit to determine the efficiency of use of the gains was conducted on the basis of the unit profits referred to 1m2 of the usable area of the building, accepted as 4.5W/m2 The average internal gains flux adopted to estimate the efficiency of gains use was 7kW, and the average daily exchange of air was 0.3h−1 for the same purpose. The efficiency of the node and installation (node, transmission, emission and regulation of heat supply) was determined on the basis of the heat losses of the node components and the installation, and the heat flux taken from the heat network as well as the efficiency of regulation and the adopted emission. Based on dependencies (23, 26), estimated U-values for partitions and windows and known building dimensions, the temperature $t¯e**$ was determined. The efficiency of using internal gains was estimated on the basis of the total heat gains of the building: ϕu + ϕZ + ϕZ,tr (from people, appliances and solar radiation) and heat demand for compensating the heat loss of the building (penetration through partition walls and ventilation) ϕdem, using dependency:

##### (36)
$ηu=1−e−ϕdemϕu+ϕZ+ϕZ,tr$
where:

ϕu – measured, average utility heat flux for considered day, W,

ϕZ – average internal heat gains flux for considered day, W,

ϕZ,tr – average for considered day, differential heat flux from solar radiation through transparent components,

The method check was performed by calculating the average monthly required power of the node on the basis of $HTrz$ determined by the multiple measurement method and comparing the determined volume with the average monthly heat flux delivered to the building, measured by the heat meter at the thermal node.

## 4. DETERMINATION OF LOSS COEFFICIENTS AND THE NUMBER OF VENTILATION AIR CHANGES

The average daily values of the used measurements for the group of 20 days accepted for analysis were presented in Table 1.

##### Table 1.

The average daily values of the used measurements for the group of 20 days accepted for analysis

When omitted in generalized balance equation the arrangement (12 – part 1)

$HTVrz⋅(t¯i−t¯e**)+Hb⋅t¯e1⋅(t¯i−t¯e)=ϕu(t¯e)⋅ηtrt(t¯e)+ϕZ⋅ηehg(t¯e)+ϕZ,tr(t¯e)⋅ηehg(t¯e)$
the component containing Hb (Hb = 0) and accepting: ϕZ = 7kW, $ϕZ=ϕZ⋅ηehg(t¯e)$, and also $HTV=HTVrz$ and transformation of the equation HTVj can be determined:
##### (37)
$HTVj=ϕu(t¯e1)⋅ηtrt(t¯ej)+ϕZ⋅ηehg(t¯ej)+ϕZ,tr(t¯ej)⋅ηehg(t¯ej)(t¯i−t¯ej**)$
HTVj specified for each day of the 20-day group is given in Table 2.

##### Table 2.

The values of the coefficient HTVj for a 20 day group

The order of days accepted in Tables 1, 2, 3 is related to the increasing value of the coefficient HTVj and the creation of pairs of days which differ significantly in day temperature $t¯e$. The average value of the coefficient HTV for a group of 20 days is:

Days for which high values HTV have been obtained are also characterized by high daily gains of heat from solar radiation at relatively high average daily outside air temperatures $t¯e$. It can be assumed that such states were accompanied by simultaneous intense uncontrolled opening of windows to avoid overheating of rooms. The result of excessive ventilation of the spaces could have been the operation of the c.o. with increased efficiency due to the opening of thermostatic radiator valves.

Inserting into the equation:

##### (38)
$ϕZj=HTVrz⋅(t¯i−t¯ej**)−[ϕu(t¯ej)⋅ηtrt(t¯ej)+ϕZ,tr(t¯ej)⋅ηehg(t¯ej)]$
the values of the average daily gains of heat from internal sources ϕZj can be determined for particular days from 20-day group – Table 3.

##### Table 3.

The values of the average daily gains of heat from internal sources ϕZj for 20-day group

On the basis of the determined values ϕZj compared to the expected value (ϕZ ~ 7kW), the 10-day group presented in Table 4 was selected for further analysis.

##### Table 4.

Average daily values of measurements for a group of 10 days

The expected value of internal gains ϕZ is ϕZ ~ 7 kW, which corresponds to a unit flow of 4.5 W/m2 of the heated building’s surface. This value is relevant to

• given in [3] 3.2–6.0 W/m2 for multi-family buildings and 2.5–3.5 W/m2 for single-family houses.

• given in [4] – 7.1 W/m2 for apartments and 1.0 W/m2 for staircases in multi-family buildings and 6.8 W/m2 for single-family buildings.

Daily average heat gains significantly different from expected are due to the conditions of use of the building deviating from the average. One of the conditions of success of the method is to meet the appropriate measurement conditions. These conditions have already been mentioned. It should be emphasized that it is important to include days without rain and snow and with moderate winds. If during the minimum number of days there is rain, snow or intense wind, the measuring line should be extended by the appropriate number of days. As measurement days it should also be avoided very cloudy days and days with moderate outside temperature (not higher than 5(7)°C, as an average daily temperature) and at the same time very sunny days which would be accompanied by significant overheating of rooms and / or opening of windows by residents. Generally, however, the more intense solar radiation occurs in the days preceding the measurement day, and yet on this day – but without causing significant overheating of rooms with closed windows – the greater the difference between $(ti−t¯e**)$ and $(ti−t¯e)$, which increases the accuracy of determination of loss factors $HTrz$ and $HVrz$.

On the days of measurement and at least one day before, the internal air temperature in the heated rooms should be constant. The deviation from the stabilized temperature should not exceed 0.5°C. When measuring the windows in the building and the entrance door should be closed. All automatic air diffusers should be open. If the windows are equipped with fittings with micro ventilation, the window handles should be set to this function. In the absence of automatic diffusers it is acceptable to periodically aerate the space with a duration of up to 10 minutes and a frequency not exceeding 8 hours. Compliance with these conditions affects the accuracy of the method presented.

For proper determining of coefficient Hb it is good to include measurement of the midday temperature range: $−10°C≤t¯e≤5(7)°C$. Among the days of measurement should also be several days from $−3°C≤t¯e≤3°C$.

Estimated values of average daily heat yields diverging from expectations suggest that these conditions are exceeded. Days of exceeding can not be included in the analysis of determining the building’s heat loss coefficients and the number of air exchanges.

Accepting once more ϕZ = ϕZηehg($t¯e$) and inserting values included in Table 4 into (8 – part 1), the following sequence of equations is obtained:

1. $HTVrz⋅(t¯i+6.28)=60.50+ϕZ$

2. $HTVrz⋅(t¯i+0.69)=46.37+ϕZ$

3. $HTVrz⋅(t¯i+7.02)=63.10+ϕZ$

4. $HTVrz⋅(t¯i−6.15)=29.51+ϕZ$

5. $HTVrz⋅(t¯i+6.63)=62.87+ϕZ$

6. $HTVrz⋅(t¯i−3.83)=35.99+ϕZ$

7. $HTVrz⋅(t¯i+0.51)=47.73+ϕZ$

8. $HTVrz⋅(t¯i+7.78)=68.40+ϕZ$

9. $HTVrz⋅(t¯i+0.25)=47.34+ϕZ$

10. $HTVrz⋅(t¯i+7.38)=67.99+ϕZ$

By subtracting one another from the next pair of equations (pairing equations of similar values and differing by several degrees), to obtain a positive value of the difference, we obtain:

• 1., 2.

• 3., 4.

• 5., 6.

• 7., 8.

• 9., 10.

The solution of this redundant arrangement in regard to $HTVrz$ by means of determinants leads to:

After inserting the resulting value into (13 – part 1) using a 10-day group and taking into account ϕZ = ϕZ⋅ηehg$t¯e$ it results in:

##### (39)
$Hb⋅t¯e1⋅(t¯i−t¯e1)−ϕZ$
1. Hb ⋅ (−336.81) − ϕZ = −7.83 kW/K2

2. Hb ⋅ (−48.84) − ϕZ = −7.43 kW/K2

3. Hb ⋅ (−379.21) − ϕZ = −7.15 kW/K2

4. Hb ⋅ (61.56) − ϕZ = −6.50 kW/K2

5. Hb ⋅ (−324.36) − ϕZ = −6.34 kW/K2

6. Hb ⋅ (56.78) − ϕZ = −6.04 kW/K2

7. Hb ⋅ (−48.84) − ϕZ = −5.60 kW/K2

8. Hb ⋅ (−474.80) − ϕZ = −3.83 kW/K2

9. Hb ⋅ (−66.41) − ϕZ = −3.31 kW/K2

10. Hb ⋅ (−424.41) − ϕZ = −3.20 kW/K2

The solution of this redundant system in regard to Hb and ϕZ by means of determinants leads to:

Substituting these values into the arrangement (20 – part 1 ), taking into account ϕZ = ϕZηehg$(t¯e)$ for the following days of the considered group gives:

The solution of this arrangement are the following values of heat loss coefficients:

The value of corresponds to the number of air exchanges for $t¯e=0°C$:

For the corrected efficiencies ηtrt, to the condition shown in Table 5, and using the given procedure, it is obtained:

##### Table 5.

The average daily value of the 10-day group with corrected ηtrt

The value of ϕZ, after considering the average efficiency ηehg = 0.95, is: ϕZ = 6.86 kW ≈ 7kW.

Designated Hb have negative value.

A positive value would mean an increasing average number of air exchanges with an increase in the average temperature $t¯e$. The explanation for such a result can only be a change in the airflow resistance of the ventilation system caused by increasing unsealing, tilting, opening of the windows, with the increase of the average temperature $t¯e$. In case of constant tightness of this system and close to the constant influence of wind on the building, the change in a function of the temperature $t¯e$ the number of ventilation air changes nV ($t¯e$) should be close to the dependence:

##### (40)
$nV(t¯e1)nV(t¯e2)∼t1−t¯e1t1−t¯e2$
so it should grow with decreasing temperature $t¯e$.

The coefficient Hb is the result of the dependence:

##### (41)
$Hb⋅t¯e1⋅(t¯i−t¯e1)=[nV(t¯e=0)−nV(t¯e1)]⋅cp⋅ρp⋅(t¯i−t¯e1)$
##### (42)
$Hb⋅t¯e1=[nV(t¯e=0)]⋅nV(t¯e=0)−nV(t¯e1)nV(t¯e=0)⋅cp⋅ρp$
##### (43)
$Hb⋅t¯e1=[nV(t¯e=0)]⋅[1−ti−t¯e1ti]⋅cp⋅ρp$

The coefficient Hb is equal to :

##### (44)
$t¯e1≠ti;Hb=HVrzVB⋅t¯e1⋅[1−ti−t¯e1ti]$
and

$t¯e=0$; Hb = 0.

These dependencies refer to the turbulent air flow in the ventilation system components.

In the case of transitional air flow, the result is:

##### (45)
$t¯e1≠ti;Hb=HVrzVB⋅t¯e1⋅[1−ti−t¯e1ti3/2]$
and

$t¯e$; Hb = 0.

For example, for ti = 20°C and $t¯e1=−10°C$ coefficient Hb is:

For the considered conditions $t¯e1=−10°C$ in case of poor external windiness, the value |Hb| > 0.00178 (0.00248) kW/K2 is a sign of increased ventilation capacity due to increased opening of the air distribution slots (air flow control elements) or window opening. On the other hand, the value of |Hb| ≈ 0.00178 (0.00248) kW/K2 is the sign of absence of change in airflow resistance through the ventilation system (ventilation slots or leakiness in windows, air gaps between rooms, exhaust grates and ventilation ducts and air flow from ducts) for the considered external conditions.

In the case of value Hb > 0, the area of its validity does not for sure cover the whole range of the outside temperature of the heating period, but only a certain part thereof, for example determined by temperature . Usually at low temperatures, the tightness of the windows remains constant – the windows are closed, possibly with microventilation at constant setting. For Hb < 0 it is possible to assume the range of change of ventilation intensity according to Hb for the entire temperature range $t¯e$ of the heating season. Hb < 0 means that the change in intensity of the ventilation coincides in some way with the change in the density of the outside air, mainly caused by the change in temperature $t¯e$, with the tendency to decrease in intensity with decreasing the difference.

## 5. VERIFICATION METHOD

In order to initially verify the method the average monthly heat consumption streams for the three months (February 2012, January 2013 and February 2013) in the considered building on the basis of the results obtained for determination of heat loss coefficients. Designated monthly streams of heat consumption were compared with the results of the heat consumption measurements for central heating in the considered building. In the verification the inclusion of these months was related to the availability of measured data. Monthly average heat consumption stream for installation c.o. can be generally expressed as:

##### (46)
$ϕum=HT•rz(t¯em)⋅(ti−t¯em)+HVrz⋅(ti−t¯em)+Hb⋅t¯em⋅(ti−t¯em)−(ϕZm+ϕZ,trtm(t¯em))⋅ηehg(t¯em)ηtrt(t¯em),kW,$

The values $HT•rz(t¯em)$ for each month are derived from the dependency:

##### (47)
$HT•rz(t¯em)=HTrz−An⋅Un+An⋅Un⋅ti−t¯nmti−t¯em$
where

$t¯nm$ – the average monthly temperature in the basement.

The data for the months included in the verification are given in Tables 6 and 7.

##### Table 6.

Values of the coefficients $HT•rz(t¯em)$ used to verify the method

##### Table 7.

Selected midday values included in the method verification

Table 6 presents the values of coefficients $HT•rz(t¯em)$ determined on the basis of the ceiling area above the unheated basement An= 380 m2 and coefficient Un = 1.50 W/(m2K) and for estimated on the basis of the regulations in force during the construction period.

The monthly (for February 2012, January 2013 and February 2013) average power of the heat distribution node working in the analyzed building determined on the basis of the values determined by the method of multiple measurements for the months analyzed is: The relative divergence of the monthly (for February 2012, January 2013 and February 2013) average power of the heat distribution node working in the analyzed building is:

The relative divergence of the monthly (for February 2012, January 2013 and February 2013) average power of the heat distribution node working in the analyzed building is:

$|Δ(ϕum)o|=|(46.94−44.12)/46.94|=0.060|Δ(ϕum)o|=|(39.80−40.63)/39.80|=0.021|Δ(ϕum)o|=|(36.23−34.04)/36.23|=0.060$

Taking into account in the building’s heat balance the average monthly flow of heat gains from solar radiation through the non-transparent partitions it can be written:

The differences in calculated and measured values disclosed in the analysis are the result of the simultaneous influence of many factors. The most important are:

##### (48)
$ϕum=HT•rz(t¯em)⋅(ti−t¯em)+HVrz⋅(ti−t¯em)+Hb⋅t¯em⋅(ti−t¯em)−ϕZssm−(ϕZm+ϕZ,trm(t¯em))⋅ηehg(t¯em)ηtrt(t¯em),kW$
• the inaccuracy of thermal efficiency prediction of node and central heating installation and the efficiency of using internal heat gains,

• deviations of the actual average temperature of heated rooms from assumed in the calculation,

• variation of internal gains for individual measurement days,

• the inaccuracy of the received measurements of the amount of solar radiation due to about 10 km distance of the measuring station from the considered building,

• some inaccuracies in the measurement of the amount of heat taken from the heating network by the heating node: the standard measuring equipment of the nodes was used for the measurements,

• the difficulty in precise establishing for any day of the year the degree of utilization of heat gains from solar radiation in the heat balance of the building for the day being considered (from 0 h to 24 h) and some of gains that come from the gains from the previous day.

These mentioned and other deviations were related to the lack of control of the building’s exploitation during the measured days. The measured values of these days were archival values. The basic criterion for selecting these days was the completeness of the required measurement data. Therefore, it may be assumed that the days taken into account did not fully meet the required conditions. The major disadvantage of the data measured such way was the lack of information on the use of the building on the days of measurement. This concerns the actions of the occupants such as: time spent in building (long-term or short or full absence in the apartment during the measuring days) intensity of ventilation of the apartments (long-term windows opening or short-term half-opening), indoor temperature (constant, variable, normal or reduced), type and intensity of activities (typical or unusual). Any deviation from the average of these operating characteristics, covered by lack of knowledge, was failing to meet the required conditions for proper use of the method and was a significant impediment to proper use of this method.

## 6. CONCLUSIONS

The publication presents the original method of multiple daily measurements for determining the values allowing to designate the building’s thermal characteristics. The multiple daily measurement method allows to specify the following main values: total heat loss coefficient $HTVrz$, heat loss coefficient associated with heat transfer through the external partitions $HTrz$, heat loss coefficient for ventilation $HVrz$ and associated with it the amount of ventilation air exchange n [1/h]. The complementary value is the coefficient Hb of variation in ventilation intensity with the change in outdoor temperature.

## ACKNOWLEDGEMENTS

08/010/BK_17/0024 (BK 227/RIE1/2017) and 08/10/BKM-16/0019 supported by Polish Ministry of Science and Higher Education

## References

1. H. Foit, A. Świerc., & W. Foit, (2017) Contribution to the short measurement method for determining the termal characteristics of existing buildings: total heat transmission coefficient estimation based on double measurement, Architecture Civil Engineering Environment – ACEE, 10(1), 103–115.
[CROSSREF]
2. H. Foit, Ś. Agata & M. Szewczyk, (2013). Equivalent external temperature for energy performance of buildings, Rynek Energii, 106, 108–112.
3. Rozporządzenie Ministra Infrastruktury dnia 6 listopada 2008r. w sprawie metodologii obliczania charakterystyki energetycznej budynku i lokalu mieszkalnego lub części budynku stanowiącej samodzielną całość techniczno-użytkową oraz sposobu sporządzania i wzorów świadectw ich charakterystyki energetycznej, Dz. U. Nr 201, poz.1240. (Regulation of the Polish Minister of Infrastructure Methodology for calculating the energy performance of building).
4. Rozporządzenie Ministra Infrastruktury i Rozwoju z dnia 3 czerwca 2014 r. w sprawie metodologii obliczania charakterystyki energetycznej budynku i lokalu mieszkalnego lub części budynku stanowiącej samodzielną całość techniczno-użytkową oraz sposobu sporządzania i wzorów świadectw charakterystyki energetycznej. (Regulation of the Polish Minister of Infrastructure Methodology for calculating the energy performance of building).

### REFERENCES

1. H. Foit, A. Świerc., & W. Foit, (2017) Contribution to the short measurement method for determining the termal characteristics of existing buildings: total heat transmission coefficient estimation based on double measurement, Architecture Civil Engineering Environment – ACEE, 10(1), 103–115.
[CROSSREF]
2. H. Foit, Ś. Agata & M. Szewczyk, (2013). Equivalent external temperature for energy performance of buildings, Rynek Energii, 106, 108–112.
3. Rozporządzenie Ministra Infrastruktury dnia 6 listopada 2008r. w sprawie metodologii obliczania charakterystyki energetycznej budynku i lokalu mieszkalnego lub części budynku stanowiącej samodzielną całość techniczno-użytkową oraz sposobu sporządzania i wzorów świadectw ich charakterystyki energetycznej, Dz. U. Nr 201, poz.1240. (Regulation of the Polish Minister of Infrastructure Methodology for calculating the energy performance of building).
4. Rozporządzenie Ministra Infrastruktury i Rozwoju z dnia 3 czerwca 2014 r. w sprawie metodologii obliczania charakterystyki energetycznej budynku i lokalu mieszkalnego lub części budynku stanowiącej samodzielną całość techniczno-użytkową oraz sposobu sporządzania i wzorów świadectw charakterystyki energetycznej. (Regulation of the Polish Minister of Infrastructure Methodology for calculating the energy performance of building).