Analysis of the Frost-Induced Damage of Building Enclosures on the Territory of the Czech Republic

A comprehensive analysis of environmental loads across the Czech Republic in terms of frost-induced damage is presented. Computational simulation of hygrothermal performance of eleven characteristic types of building envelopes, composed of both contemporary and historical materials, is performed at first. -e exterior boundary conditions of the computational model are defined by a set of weather data characterizing the environmental conditions in the Czech Republic, which are acquired from 64 weather stations. -e results of hygrothermal simulations are assessed using several specific damage functions. In this way, the basic datasets for the frost damage analysis are obtained. -eir application as input parameters of a specially developed correction procedure based on elevation makes them possible to obtain a continuous coverage of the geographic area of the Czech Republic. Finally, isopleths of the supposed frost damage are drawn, depending on the envelope type, and damage maps are produced which may help the engineers to enhance the building envelope design process. -e presented results indicate the necessity of paying attention to local environmental loads in the building enclosure design process and reveal both critical and favorable locations from the point of view of frost-induced damage to buildings.


Introduction
e primary objective of hygrothermal design is to ensure healthy and comfortable place for occupants, which must comply with the requirements for energy efficiency.Last but not least, a serious hygrothermal design should ensure a sufficient durability of the constructions in terms of their resistance to the environmental and climatic load [1][2][3][4].
ose loads are acting on the outer surface of the buildings, thus causing long-term damage to the construction, for example, in the form of corrosion or frost-induced damage [5,6].e frost-induced damage is one of the most typical long-term damages of building enclosures in the Middle Europe.It affects most types of building enclosures and is manifested by various kinds of failures ranging from surface cracks to structural damage [7,8].Given that building materials contain a wide range of pores, the damage is typically initialized when liquid water present in the pore system of the material is exposed to subzero temperatures.It is initiated through the nucleation, growth, and interaction of microcracks.ese processes usually occur internally and are manifested by volume expansion during cooling [9,10] and by macroscopic cracks that develop after repeated cycles of freezing and thawing.From this point of view, almost every porous building material that is exposed to the environment becomes vulnerable.
e thermal and hygric properties of construction materials are dependent on actual temperature and moisture content affecting the hygrothermal performance of the entire wall system or even the whole building [11][12][13][14].erefore, the frost-induced damage of building enclosures should be considered as a joint influence of two input factors.It is not only the material properties that need to be taken into account but also the boundary condition acting on the material's surface must be incorporated into the frost damage assessment procedures.Ignoring the characteristic climatic conditions of each location may lead to an incorrect selection of construction materials and enclosure design in terms of their durability.Moreover, building codes usually introduce constant thermal or hygric material properties, which are common for all locations and associated with standardized temperature and relative humidity conditions in the materials.e liquid moisture transport in the material is often even not considered at all.Such a definition does definitely not represent the actual operating conditions of the building enclosures.For that reason, it is quite a challenging task to properly assess the frost durability of wall assemblies in given locations.As the designers cannot rely on the traditional empirical rules and standardized methods due to the high variability of the environmental and material properties or insufficient quality of input data, more advanced techniques should be incorporated in order to assess the frost durability of building enclosures properly.
As the climate in the Middle Europe in not uniform, it is logical that even on the scale of individual countries the environmental conditions may significantly vary as well.In some countries, there can exist such localities, where local climate will not produce sufficient conditions for the frost damage to be initiated, as well as localities with so severe conditions that will induce abnormal frost load to the constructions.It should be noted that also the elevation changes over individual countries may imply various frost load levels.erefore, some serious analyses of local regions in those countries should be performed providing basic information on the climatic loads, for example, in the form of the isopleth or contour maps.
e isopleth maps providing the engineers with design thermal conductivity for the region of Northeastern Spain was presented by Pérez-Bella at al. [15].e authors analyzed weather stations across the investigated territory and calculated correction factors for the design thermal conductivity as a function of geographic coordinates.Similarly, the isopleth maps were used as an output of the investigation of water penetration risks in building facades in Brazil [16].
ose maps were also used for the territorial analysis of wind-driven rain and driving-rain wind pressure for several Spanish regions [17].Territorial investigation of air pollution impact on construction materials in the Slovak Republic was presented in another research [18].e authors of this study generated the maps of annual average surface loss of stone materials and carbon steel.In all of the abovementioned studies, the authors processed the data from sufficient number of weather stations providing reasonable data coverage allowing them to construct the isopleth maps and to draw scientific conclusions.
e main objective of this paper is to analyze the local frost loads on the territory of the Czech Republic and to provide the designers and engineers with an efficient tool displaying the critical locations in the country.To achieve that objective, 64 locations across the Czech Republic were monitored using meteorological stations and the weather data of those locations were recorded.en, the computational analysis was involved, and the effect of weather data was investigated on different types of building enclosures using several damage functions.In this way, the basic datasets for the frost damage analysis were obtained.Finally, the isopleth maps were generated using interpolation techniques, and an additional correction procedure was applied in order to correlate the obtained results in the locations insufficiently covered by the weather data.

Materials and Methods
2.1.Climatic Data.For the investigation of the frost loads in the Czech Republic, 64 locations across the country were selected and the weather data from those locations were collected.As the different parts of the Czech Republic experience both hot summers (in the region of South Moravia) and cold winters (in the mountain regions), the selection covered both lowlands and mountains across the country.All weather data were obtained from the Czech Hydrometeorological Institute, which is the official authority for meteorology, climatology, hydrology, and air quality protection in the Czech Republic.ey included hourly values of temperature, relative humidity, precipitation, wind direction, wind velocity, diffuse and direct shortwave radiation, sky longwave emission radiation, and longwave emission radiation and were provided in the form of the test reference year (TRY) [19][20][21].e list of involved locations together with their elevations is shown in Table 1.e map of weather stations is depicted in Figure 1.

Studied Building Enclosures.
e investigation of the local frost loads in the Czech Republic was carried out for various types of building enclosures covering both historical and contemporary wall assemblies.
e loadbearing materials included autoclaved aerated concrete (AAC), concrete, ceramic brick, advanced hollow bricks, and sandstone.e contemporary building envelopes were provided with different types of thermal insulation layers based on polystyrene and mineral wool.e historical masonry did not have any thermal insulation.e exterior plasters were chosen with respect to the material composition of the envelopes, such as lime-cement plaster (LC), renovation plaster for historical masonry (RPHM), or lime-pozzolan plaster that was specially developed for the advanced hollow bricks (LPC).On the interior side of all structures, 10 mm thick lime-cement plaster was assumed.e list of studied building enclosures is shown in Table 2.

Computational Simulation.
All the simulations were performed under time-dependent boundary conditions using the finite element method.Computer simulations of hygrothermal performance of building envelopes listed in Table 2 were conducted for the time period of ten years using the HEMOT simulation tool (HEat and MOisture Transport), which is a preprocessing tool for the general finite element package SIFEL (SImple Finite Elements) [22].In the simulations, a slightly modified version of Künzel's [23] mathematical model of coupled heat and moisture transport was used: Advances in Materials Science and Engineering where is the thermal conductivity, T (K) is the temperature, δ p (s) is the water vapor permeability de ned by Fick's law, p v (Pa) is the partial pressure of water vapor in the air, ρ w (kg•m −3 ) is the water density, w (m 3 •m −3 ) is the water content by volume, n (−) the porosity of the porous body, M (kg•mol −1 ) the molar mass of water vapor, R (J•K −1 •mol −1 ) is the universal gas constant, t (s) denotes time, and D g (s) is the global moisture transport function.e exact sources of hygric, thermal, and basic physical properties of building materials [24][25][26][27][28][29][30][31][32], which were used in the simulations as input parameters, are summarized in Table 3 in detail.

Applied Damage Functions.
e hygrothermal performance of studied building envelopes listed in Table 2 was evaluated using several damage functions, which were introduced recently [32].
ese damage functions allow a relative assessment of severity of environmental loads in terms of frost-induced damage.In this way, various regions of the Czech Republic and their local climatic loads can be analyzed.e applied damage functions are based on the assessment of temperature and moisture content distribution across the wall system over time with respect to pore characteristics of individual material layers.In this paper, the functions evaluate the performance in a particular point of the  Advances in Materials Science and Engineering wall assembly further referred to as "Point of investigation." is point of investigation was placed 2 mm under the exterior surface of each studied wall.ree different damage functions were applied throughout the paper, namely, time-of-frost (TOF), amount-of-solidified-water (ASW), and number of indicative freeze/thaw cycles (IFTC).However, the assessment using abovementioned damage indexes can be further modified, so that not only one point is investigated.For example, one can investigate several selected points, and then their average can be used as the function output.e definition of damage function output should be done with respect to particular construction type and application purposes.
e mathematical description of the applied damage functions is given below.However, at first, some physical simplifications need to be introduced, in order to ensure a simple applicability of those functions and to keep their fast execution rate: (i) water primarily fills the pores from smallest to larger ones; (ii) when the hygroscopic moisture content in the investigated point is exceeded, all excess moisture is considered to be in the liquid phase; (iii) the freezing point of water depresses according to the size of largest water-filled pore; (iv) part of the water in the smallest pores remains unfrozen due to freezing point depression caused by pore curvature; and (v) all liquid water which is subjected to the water/ice phase change is solidified immediately, that is, the dynamics of ice formation is neglected.
Using the assumptions (i) and (iii), the radius of largest water-filled pore can be evaluated from the pore size distribution function when the moisture content in the investigated point is known and the freezing point depression ΔT f in a cylindrical pore of radius R can be then expressed by the Gibbs-omson equation as follows: where c sl (mJ•m −2 ) is the surface-free energy (interfacial tension) of the solid/liquid interface (c sl � 31.7 ± 2.7 mJ•m −2 [33]), ] l (cm 3 •mol −1 ) is the molar volume of the liquid (] l �18.02 cm 3 •mol −1 ), and Δh b (kJ•mol −1 ) is the melting enthalpy in the unconfined (bulk) state (Δh b � 6.01 kJ•mol −1 ); all quantities were taken at the bulk coexistence temperature T 0 (T 0 � 273.15 K). e TOF damage function was inspired by the commonly used time-of-wetness damage function designed for the analysis of corrosion-related degradation [34][35][36].TOF calculates the number of hours during the year when the conditions in the investigated point of the wall cross section are favorable for ice formation, that is, the temperature is below the critical temperature while the moisture content is above the critical value.TOF ranges between 0 and 8760 and can be expressed as where T L and w L are critical values of temperature and moisture content and T i and w i are hourly values of temperature and moisture content, respectively.Note that the summation in Equation ( 4) uses the iverson bracket [37], which is a notation that denotes a number that is 1 if the condition in square brackets is satisfied and 0 otherwise.TOF is calculated only when both T i < T L and w i > w L .e critical values T L and w L are defined as where w hyg is the hygroscopic moisture content.e ASW damage function returns the annual amount of liquid water retained in the investigated point that is solidified under critical temperature.ASW works with the assumption that when the temperature drops below the critical level T L , and liquid water is present in the material (i.e., w i > w L ), and only a part of that amount in the where w i,unfrozen (m 3 •m −3 ) is the unfrozen moisture content retained in the smallest pores of the material.e unfrozen moisture content is evaluated in each calculation step in two phases.First, the critical radius R is calculated reciprocally from Equation (3) from the known temperature T i in the investigated point.en, the unfrozen moisture content w i,unfrozen is evaluated using the pore size distribution function and the known R from the previous step.
e number of indicative freeze/thaw cycles (IFTC), as the last damage function, returns a hypothetical number of freeze/thaw cycles that may occur in the investigated point in the wall assembly during the simulation year.e damage function is inspired by various frost durability tests given in the national standards [38][39][40], where the samples are cyclically loaded by freezing and thawing for certain time periods.e standardized freeze/thaw cycle is mostly dened by 2 hours of freezing and 2 hours of thawing.Such a de nition was adopted by the IFTC damage function.e critical conditions for ice formation are monitored in the point of investigation.
ey come from the same critical values of temperature and moisture content as in TOF or ASW (Equation ( 5)).In order to account the freeze/thaw cycle as valid, the freezing conditions in the point of investigation must take at least 2 hours (i.e., TOF ≥ (2) while any two consecutive freeze/thaw cycles must be separated by at least 2 hours of thawing period, that is, the conditions are not favorable for ice formation.If two or more freezing periods longer than 2 hours are separated by one or more periods of thawing shorter than 2 hours, it is still accounted as one freeze/thaw cycle.

Output Data
Postprocessing.e computational simulations were performed for the combination of 11 wall assemblies under 4 di erent orientations together with 64 di erent sets of boundary conditions, giving 2816 output les.Each output le was analyzed by 3 di erent damage functions producing 8448 damage function outputs.In order to reduce this number, the outputs were rst averaged by orientation and then grouped by the load-bearing material into 5 categories-AAC, ceramic brick, concrete, advance hollow clay brick, and sandstone (Table 2).All data were then processed by Surfer 14 which is an advanced tool for visualizing, mapping, and modelling in both 2D and 3D.For each set of damage function outputs consisting of 64 values with assigned X, Y coordinates (weather stations), Advances in Materials Science and Engineering regular mesh grids were generated using the Kriging method with a linear variogram model [41].Each grid consisted of 701 × 301 grid nodes providing reasonable resolution for generation of the isopleth maps.
As the total area of the Czech Republic is approximately 79 000 km 2 , the gridding from 64 weather stations was insu cient as each weather station covers approximately the area of 1 200 km 2 .For that reason, a correction of the generated grids was necessary to obtain satisfactory results.
e correction was based on assumption that frost damage indexes correlate with the elevation.us, for each category of load-bearing masonry, the damage function outputs were plotted against the elevation of particular weather stations and linear regression was performed.e slope of identi ed linear function was then used as a correction coe cient of generated mesh grid.e values in each grid node were then corrected according to following formula: where DF corr is the corrected value of damage function output, DF is calculated damage function output, c is the correction coe cient, A is the real elevation, and B is the interpolated elevation from the weather station altitudes.e scheme of correction procedure is shown in Figure 2. e results of linear regression for AAC are shown in Figure 3. e correction coe cients for all group categories are summarized in Table 4. e real elevation A was obtained from publicly available data.ey included all cities and villages in the Czech Republic, together with the known altitudes of all signi cant hilltops and mountain peaks.In total, 605 elevation data on cities and villages and 1127 elevation data on hilltop and mountain peaks were collected.

Results and Discussion
e isopleth maps displaying various damage function outputs were generated to characterize the frost-induced damage in the Czech Republic, using the methods described above.In total, 15 di erent maps were generated but due to similar damage function outputs of AAC-and concrete-based walls, the concrete results were excluded from the gures.e distribution of damage function outputs across the Czech Republic is shown in Figures 6-8.It is obvious that the hygrothermal performance of all studied walls correlates with each other, and the service life of surface layers is signi cantly a ected by the material of load-bearing structure.From that point of view, the sandstone masonry proved the highest resistivity to severe environmental conditions.On the other hand, the brickwork made of advanced hollow clay brick seemed to be the most susceptible to the frost-induced damage of its surface.For example, a detailed analysis of IFTC showed that the     Advances in Materials Science and Engineering building enclosures made of advanced hollow clay bricks exhibited twice higher number of freeze/thaw cycles in surface layers than enclosures made of sandstone.For that reason, a particular attention needs to be paid to the frost damage assessment of these kinds of building walls where the choice of a proper material for exterior renders is critical.e isopleth maps also indicated excessive frost loads in the mountain areas, which was expected.e critical areas were identi ed especially in Jeseníky, Šumava, and Krkonoše mountains.On the other hand, the lowest values of damage functions, that is, the mildest conditions for the frostinduced damage, were observed in the lowlands around the river Labe and in the Silesian region or South Moravia.e most severe conditions were observed in the location of Šerák in Jeseníky mountains (1328 m a.s.l.) and Filipova Hu ť in Šumava mountains (1110 m a.s.l.).e most favorable conditions were identi ed for the location of Ústí nad Labem (375 m a.s.l.) and Doksany (158 m a.s.l.), both located in the northern part of the Czech Republic in the valley of the Labe river.
As the applied damage functions were primarily designed for the relative assessment, that is, for a comparison of individual environmental loads, the obtained data should be confronted with eld or laboratory measurements.From this point of view, the IFTC damage function seems to be the most convenient as it was designed based on standardized laboratory frost damage tests [38][39][40].However, the principle of the testing methods consists in exposing fully saturated samples to alternating environment ranging from 20 to −23 °C for certain time period, which does not correspond well with the eld conditions.Even if the freezing and thawing time periods used in IFTC calculation was adopted from the standardized procedures, the materials under eld conditions become fully saturated on rare occasions only.Moreover, if the saturation occurs, it covers only minor part of the construction for a short period of time.Also, the eld temperatures drop below -20 °C only occasionally, contrary to the standardized temperature loads.For this reason, the laboratory results of frost durability provide an incomplete information on the materials' durability.e materials loaded by real weather conditions may stand signi cantly higher number of freeze/thaw cycles than indicated by laboratory tests.erefore, in future studies, it would be useful to link the damage function outputs with the results of eld measurements, providing thus the designers with even more precise datasets improving the frost damage assessment.

Conclusions
A comprehensive characterization of the geographic area of the Czech Republic in terms of the frost-induced damage to building enclosures was presented in the paper.e obtained results allowed for the generation of a series of isopleth maps  Advances in Materials Science and Engineering displaying the local frost loads across the country.e maps can be easily used by designers and engineers as a support in the building enclosure design.
e main findings can be summarized as follows: (i) e service life of exterior surface layers of building enclosures was significantly affected by the material of load-bearing structure.e sandstone masonry proved the highest resistivity to severe environmental conditions, while the brickwork made of advanced hollow clay bricks was the most susceptible to the frost-induced damage of the surface.(ii) e most severe conditions for the ice formation in building materials were observed for the mountain locations of Šerák and Filipova Hu ť, the mildest conditions were observed for the lowland location of Ústí nad Labem and Doksany.(iii) e gathered data should be supplemented with the results of field measurements rather than laboratory tests, in order to enhance the applicability of the isopleth maps.

Figure 1 :
Figure 1: Map of involved weather stations in the Czech Republic.

Figure 2 :
Figure 2: Scheme of the correction procedure.

Figure 3 :
Figure 3: Linear regression of AAC damage function outputs plotted against elevation.

e
data were combined with the weather station elevations as well, in order to obtain zero corrections for the grid nodes representing those stations (A − B 0). e map of gathered control points and generated elevation are shown in Figures 4 and 5, respectively.

Table 1 :
List of applied weather data.

Table 2 :
List of studied building enclosures.

Table 3 :
List of sources of input parameters for computational simulation.