Moisture Buffering in Surface Materials Due to Simultaneous Varying Relative Humidity and Temperatures: Experimental Validation of New Analytical Formulas

: Buildings are subjected to the indoor environment, especially in non-controlled climates. Temperature and humidity variations might effect or even damage materials sensitive to moisture. For this reason, it is important to understand the response of hygroscopic materials to variable indoor environmental conditions. Existing methods looked into the dynamic sorption capacity of materials, by analysing the impact of only humidity ﬂuctuations, with temperature usually considered non-inﬂuential or non variable. However, temperature ﬂuctuations may impact the moisture capacity of the materials, as materials properties might substantially vary with temperature. Moreover, in existing protocols, the humidity variations are considered to be varying under square wave ﬂuctuations, which may not be applicable in environments, where the indoor is inﬂuenced by daily and seasonal climate variations, which presents more complex ﬂuctuation. In this study, a simulation method that can predict the impact of environmental condition on materials under simultaneous temperature and humidity ﬂuctuations was developed. Clay and gypsum plaster were analysed in the numerical model and results were then validated with experimental data. Materials were subjected to either sinusoidal and triangular temperature and RH variations and different cycle time intervals. The investigation of sinusoidal and triangular environmental variations pushed to a better understanding of materials response to different environments and to the improvement of the simpliﬁed model. The development of a simpliﬁed model can realistically predict the potential future impact of climate changes on buildings without the use of complex and memory demanding computational methods.


Introduction
To regulate temperature and humidity in buildings, it is not always possible to alter the building envelope or to install air conditioning systems, especially in historic buildings, such as churches. Consequently, indoor surfaces and hygroscopic building materials are subjected to daily and seasonal temperature and Relative Humidity (RH) fluctuations, which might degrade buildings [1,2]. For these reasons, it is important to understand the impact of temperature and humidity on buildings, as well as the response of materials to climate variations [3,4].
Indoor climate in buildings without any conditioning system is mainly determined by the outdoor climate and hygrothermal performance of the building enclosure [5]. The lack of heating and well insulated walls might generate significant daily or seasonal indoor humidity and temperature fluctuations, depending on the climate and location of the building [6,7]. The simultaneous variations of temperature and RH can produce significant variation in moisture exchange between materials and the environment [8].
Moisture transport and storage capacity of interior surface materials due to varying humidity in the indoor air is of importance both for the humidity levels of the room itself and for the moisture impact on walls surface materials [9]. There is research [10][11][12][13], standards [14] and some ongoing measurements on this topic [15] that pointed out the importance to understand the impact of moisture to materials durability, indoor thermal comfort and health. [16] took a step further, by considering periodic humidity variations and by applying moisture buffering experimental results into simulations. However, the focus of most of the studies was mainly on the sole RH variations, while temperature was always considered non-variable. The observation of daily indoor climate variations in non-conditioned buildings suggested a simultaneous opposite and quasi-sinusoidal temperature and RH fluctuations [6,7,17]. Some recent results on this area [1,18] gave new analytical expressions for the effect of combined variations of RH and temperature. Bylund et al. [18] analysed experimentally the moisture uptake of wood and introduced a calculation model that took into account simultaneous temperature variations together with the surface resistance of the material. However, the model was based on a specific material and required an elevated number of material properties and experimental analysis.
In a newly published study [19], a detailed analytical solution for temperature and RH square shaped variations was proposed. The method was developed by neglecting the surface vapour resistance and by considering a semi-infinite material [20]. The objective of this study was to further develop the analytical model from [19], by introducing triangular and sinusoidal temperature and RH variations. To validate the analytical method, results were compared with laboratory experimental tests on clay and gypsum plasters. The materials mass variations to simultaneous temperature and RH variations was observed and compared with simulations. The aim of this paper is to introduce new experimental approach to evaluate the moisture buffering capacity of materials and use experimental data to realistically simulate the materials behaviour. The effect of the surface resistance was implemented in the model and the validity of the semi-infinite approach on the moisture uptake was analysed. By improving the model and approximating the indoor daily variations to sinusoidal and triangular, it was possible to develop a reliable simulation method that is able to predict the materials response in a non-conditioned environment.

Materials
In this study, commercially available undercoat clay and gypsum plasters were analysed, due to their good moisture buffering capacity. Samples were mixed and cured, as described in [21]. Clay and gypsum were cast in 150 mm × 150 mm × 40 mm and 150 mm × 150 mm × 20 mm moulds, respectively, in order to test the sensitivity of the model to the material thickness. Specimens were air dried for 28 days, before testing. The dry density ρ dry , porosity Φ(%), the dry cup water vapour resistance factor µ(−) and moisture capacity ξ w (kg/kg) were measured, as described in [21]. Results are summarised in Table 1. Materials were stored in an environmental room at 60%RH and 22 • C.

Experimental Design
For each specimen the change in weight of the specimens was investigated, when subjected to triangular or sinusoidal simultaneous temperature and RH variations. Specimens were tested in an environmental chamber (ACS Compact Test Chambers DY110), into which a mass balance were placed to continuously measure the change in weight, as shown in Figure 1. More details of the set-up can be seen in [8] The tests followed the general guidelines for humidity variations and test set-up of the NORDTEST protocol [14]. The materials were exposed before each tests to 24 h pre-conditioning at 23 • C and 54%RH, and to six cyclic humidity and temperature variations at an air speed of 0.1 m/s. Each cycle consisted of humidity variations from 75%RH (high RH) to 33%RH (low RH), and temperature fluctuations between 18.0 • C and 28.0 • C, which is the acceptable operating temperatures in buildings in the ASHRAE Standard 55 [22]. Temperature was set to follow inverse variation than RH, as shown in Figures 2 and 3. For both triangular and sinusoidal fluctuation, three different humidification and de-humidification time intervals were defined, as shown in Table 2.  In the Triangular 12/12 and Sinusoidal 12/12 tests, considering 54%RH as starting point, the positive section of the curve represents the humidification phase, while the section below 54% is the de-humidification interval. As with the NORDTEST [14], the RH was gradually increased, until reaching its maximum 75%RH and then back to 54%RH. In the first 6 h of the de-humidification from 54%RH the humidity reached the minimum 33%RH, and then again up to 54%RH, as shown in Figures 2a and 3a. In the 12/12 test, the climatic chamber increased and decreased the RH every hour, transitioning from one step to the other by setting a "slope", which the climatic chamber automatically generated, to progressively increase/decrease the humidity until the next RH was achieved. In the Triangular 12/12, the triangular curve shape was generated by creating 24 equally distributed RH intervals, while the Sinusoidal 12/12 was achieved, by constructing a sinusoidal symmetrical equation. The temperature curve followed the same logic as the RH variations. In 12/12 test the curve started at 23 • C, reaching the minimum temperature of 18.0 • C within the first 6 hours of the cycle. Successively, the temperature jumped back to 23 • C, to increase the temperature to 28.0 • C always withing the 24 h cycle. In the Triangular 8/16 (Figure 2b) and the Sinusoidal 8/16 ( Figure 3b) the temperature and RH curves were marked by hourly interval in the humidification/cooling phase and by 2 h intervals in the de-humidification/heating phase. The Triangular 72/72 ( Figure 2c) and Sinusoidal 72/72 (Figure 3c) follow the same principle of the 12/12 h fluctuations, but 6 h intervals were applied between one step and the other, in order to understand if slower moisture variations influence the sorption capacity of the materials. A temperature and RH sensor for each specimen (Tiny Tag TV 4505) monitored the climate condition in the climatic chamber to observe the agreement between the target fluctuations and the actual measurements on both sides of the chamber. The accuracy of the sensor was ±0.5 • C for the temperature and ±3%RH for humidity.

Simulations
Steady periodic variations with the time period t p (s) in RH and temperature were investigated. The moisture transfer took place at x = 0 of a semi-infinite domain, x > 0. There was a vapour surface resistance Z (s/m) at x = 0. The moisture buffering effect, i.e., the total moisture m A (kg/m 2 ) that flows in and then out during a time period was of interest.

Equations
The moisture balance equation reads: Here, v (kg/m 3 ) is the humidity by volume and w (kg/m 3 ) is the moisture content. At the surface x = 0 there is: Here, v b (t) is the time varying boundary humidity by volume. The analysis in this paper assumed that the temperature of the surface material always followed the chamber temperature T b (t) without any delay. This was a reasonable assumption since temperature changes are much more rapid than moisture changes and that it is only the thin interior surface layer that is affected by variations in indoor cyclic moisture variations. Equation (3) states that there is a moisture tight back side of the material at the depth D (m) representing the thickness of the material. The main part of the following assumes a semi-infinite material thickness, i.e., D = ∞.

Simplified Equations
Two simplifications were introduced. The first one is that the vapor diffusion coefficient δ v (m 2 /s) is constant: The second simplification is that the slope of the sorption curve is constant. Furthermore, hysteresis is neglected, and the the average slope is used, instead. This is justified, as the investigated materials showed a very limited difference of the slope of the isotherm within the studied RH regime. The upper and lower between the adsorption and desperation limit is 0.8% and the lower one is 0.02%. A sensitivity analysis using this span would reveal an estimate of the impact. In Section 4.2.3, Equation (28) suggests a square root dependence of the slope on the total moisture uptake.
The moisture balance equation then becomes: Here, v s is the humidity by volume at saturation. Introducing the vapor moisture diffusivity a v (m 2 /s): The boundary condition becomes: The parameter d v (m) represents an equivalent thickness of the surface resistance, i.e., a material with this thickness gives the same resistance as Z. The boundary RH is denoted by ϕ b (t)(-).

Varying Temperature
Time varying temperature was considered: Here, the material temperature is equal to the boundary temperature. It is purely a function of time: A variable substitution is introduced: The moisture balance equation is transformed to: Here, there is an equation that is similar to the one-dimensional moisture balance Equation (7) with the diffusivity a v equal to 1. The equation is linear when using this transformed time variable and superposition techniques can be used. Therefore, it is necessary with a basic periodic solution to handle complex cyclic changes in the RH.
The problem can be solved in the τ-domain using superposition technique. First, the boundary values were transformed into this domain. The following general, even functions, for the boundary RH and temperature were assumed: The maximum and minimum value of the boundary RH are: The τ is defined as: Since the integrand is always positive an inverse function can be found: The boundary RH becomes: Since a periodic and even function was assume, it becomes: where a n = 2 This description suggested that the moisture problem, in the τ-domain, can be solved by summing an infinite number of solutions, one for each cosinusoidal variation using Fourier series. This solution is found in [10]:

Moisture Uptake
The moisture flow g (kg/m 2 ·s) into the material is: The integrated m(t) (kg/m 2 ) moisture uptake from time zero to time t is: The formula can be reformulated using a variable substitution, s = τ(t), and an integration between 0 and τ p .
The moisture uptake during the uptake period becomes: In a dimensional analysis of the formula the different lengths that shows up can be observed. One important is the periodic penetration depth [10]: The following non-dimensional parameters are of importance: On top of these two parameters, the shape of the curves for the RH and temperature boundary values, determining the Fourier Series coefficient a n in Equation (27), is of importance.

Experimental Results
The moisture buffering performance of clay and gypsum under triangular and sinusoidal temperature and RH fluctuations is shown in Table 3. The average peak to peak sorption capacity was calculated by considering the last three cycle of the 12 h and 8 h, as the changed in weight stabilised after the third cycle in all tests. In general, gypsum had a better moisture capacity than clay (Figures 4 and 5), by adsorbing two times more moisture regardless of the shape and time interval of the temperature and RH fluctuations. The performance of gypsum is due to its porosity and pore structure. As explained in [8], gypsum presented higher porosity than clay, but it also showed a significant presence of micro-pores that increased themoisture uptake of the material. Differences between the 12 h and 8 h curves in the triangular case either for clay and gypsum were negligible, as the same moisture uptake was measured. Between the Sinusoidal 12/12 and Sinusoidal 8/16 curves 11% and 5% variations were observed for clay and gypsum, respectively. However, differences were small and not significant, when compared to the differences between triangular and sinusoidal sorption responses. There was 18% and 20% differences between Triangular 12/12 and Sinusoidal 12/12 peak to peak sorption capacity, while between Triangular 8/16 and Sinusoidal 8/16 clay and gypsum presented 10% and 13% variations, respectively. The Triangular 72/72 and Sinusoidal 72/72 tests presented the most significant discrepancy, as clay and gypsum adsorbed 20% and 26% more in the sinusoidal case, respectively. The considerable higher sorption capacity of the 72/72 tests is related to the longer time interval, which allowed materials to adsorb more water.  As [8] also mentioned in a study on sinusoidal environmental variations, the triangular and sinusoidal curves presented a delayed response to the temperature and RH fluctuations ( Figure 6). The hygric lag varied depending on the humidification and de-humidification intervals, the length of each step, the material's characteristic, but especially on the simultaneous temperature variations. Table 4 summarises the hygric lags in all the tests, where the first value represents the time between the RH-peak and the moisture uptake peak, while the second represents the time between the lowest RH and the lowest moisture uptake value. In the Triangular and Sinusoidal 12/12 clay and gypsum showed different lags, in which clay always presented a slower response than gypsum either in the adsorption and desorption. Regardless of the materials differences, both plasters showed an asymmetry between the humidification and de-humidification. The desorption lag was 30 min for clay, and 20 and 60 min for gypsum in the Triangular 12/12 and Sinusoidal 12/12, respectively. The asymmetry indicated that in the de-humidification the response of the plasters is quicker than the adsorption, probably due to the effect of the air movement in the chamber that slowed the humidity uptake and increase the moisture release.   An opposite pattern can be seen in the Triangular 8/16 and Sinusoidal 8/16, where materials presented higher hygric lags in the desorption than the adsorption. The reason could be related to the shorter humidification interval (8 h) and longer de-humidification (16 h), as shown in Table 4. Moreover, the the longer exposure to lower humidifies and higher temperature may slow down the moisture release. The Triangular 8/16 the plasters had the same lag in the humidification, but in the desorption clay's delay increased significantly. In the Sinusoidal 8/16, clay and gypsum uptake curves had similar lags, which is in line with [8] test on the same gypsum sample.
The 72/72 test showed the biggest hygric lag, probably due to the longer steps, which allowed the materials to adsorb more moisture. In Triangular 72/72 and Sinusoidal 72/72 it was possible to see substantial difference between clay and gypsum. Clay was slower to adsorb water than gypsum, and it took more time to respond to the change of humidity in the air, due to its lower sorption capacity and water vapour permeability. Another important factor that might contribute to the slower response of clay was the different thicknesses of the materials. Clay was 20 mm thicker than gypsum, which means had more moisture storage capacity of gypsum with its 20 mm thickness. This is less noticeable for the 12 h and 8 h due to the quicker humidity variations, while in the 72/72 tests the specimens have got more time to adsorb water from the environment. As [23] investigated, the amount of moisture that can be adsorbed by materials is strongly dependent on their thickness and penetration depth. The penetration depth depends on the material characteristics and it determines the moisture buffering potential of materials. When the thickness of a material is smaller than its penetration depths, materials can buffer less moisture that what they potentially could adsorb.
To compare the experimental results to simulation, an uncertainty analysis was performed, to check the variability of the moisture uptake, due to experimental error. The size of the specimens and the variability of the moisture uptake were considered. The uncertainty of the single measurements were first analysed and than combined to calculate the overall variations of the moisture uptake. Table 5 shows the results of the analysis. Four measurements for each dimension were considered for the calculation of the surface area of the specimens, while three repeated moisture buffering tests were used for the climatic chamber and mass balance uncertainty calculation. Delays in the signal transmission between the mass balance and data logger were not detected. The temperature and RH monitoring in the climatic chamber showed a good agreement between the actual climate conditions with the targeted curves both for the triangular and sinusoidal curves. The Triangular tests are presented in Figure 7. In Figure 7a the RH presented on average 2%RH variations from target, up to 5.6%RH during the de-humidification. while temperature showed 0.43 • C variations up to 1.40 • C. In the symmetric test (Figure 7a) there was a better match between the actual and target climate. Temperature and RH presented an average variation of 0.17 • C and 1.50%RH, respectively, up to 0.79 • C and 3.5%RH always during the de-humidification. The temperature and RH fluctuations had smaller amplitudes than the target curves (10 • C and 42%RH), as shown in Table 6. In Figure 7b it was also observed a delay of 1 h of RH in the de-humidification phase that was than zeroed in the humidification. The reduced RH fluctuation may produce small variations of the moisture buffering results of the two materials. Similar consideration can be done for the sinusoidal variations, as [8] explained, and as shown in Table 6.

Triangular 12 h/12 h Variation for Gypsum
In this section, a demonstration of the application of the method is shown. As exemplification, gypsum material exposed to Triangular 12 h/12 h variation in the boundary temperature and RH-values is presented, but the analysis was performed on all cases. The following data was assumed for the vapour diffusion coefficient and the slope of the sorption isotherm: The time period, t p (s) of the process is 24 h. The RH is be highest when the temperature is lowest and vice-verse. The variations are symmetrical around the average value of both the temperature and RH. The following boundary values are assumed: Table 7 shows the amplitude parameter f A (Equations (28) and (29)). From the table it is clear that for gypsum without any coating, represented by a default Z value of 360 s/m, reduces the moisture buffering capacity with 13%. In the results presented below a value of 195 s/m is used based on a assumed air velocity of 0.1 m/s and a convection surface coefficient based on forced convection [10].
The actual moisture buffering uptake is given by Equation (28) Tables 8 and 9 show the simulated vales for the moisture buffering uptake and the time lag. The first value for the time lag represents the time between the RH-peak and the moisture uptake peak. The second one represents the time between the lowest RH and the lowest moisture uptake value. The difference in these values clearly show the asymmetry and non-linearity of the problem. As also shown in the experimental test (Table 4), the longer lag during the de-humidification in the 8/16 h tests, compared to the humidification can be also explained by analysing Equation (39). However, the time lags in Table 9 also justify the models assumption that the surface temperature follows the room temperature. The time delay for the change in temperature in the boards, due to a room temperature change, is around 1200 s for the 20 mm thick gypsum board and 1600 s for the 40 mm thick clay board. This is roughly 1/3 of an hour. For diurnal variations the ratio between this time delay and the time period is of the order 10 −2 . For the six days period the ratio is almost one order less and can further justify the approximation.  Figures 8 and 9 show the moisture uptake function m as a function of t/t p . Time zero corresponds to the time of the RH peak, and the start time of integration, i.e., we assume the uptake of moisture being zero at time zero, which is an arbitrary choice of level. Consequently, it is possible to follow the accumulation and the drying out of moisture of the specimens during the time period. The two curves with the largest amplitude represent the case with a six days variation. It is clear that the diurnal variations do not significantly differ between 12/12 h and 8/16 h cases. As also shown in the experimental test (Table 4), the longer lag during the de-humidification in the 8/16 h tests, compared to the humidification can be explained by analysing Equation (39). If a 8/8 h period was analysed, instead of the 8/16, the time t p period would correspond to 16 hours. For the 16/16 h test would be i.e., 32 h. The difference in uptake m A would then differ between the 8/8 and 16/16 h approximately of √ 32/16 = √ 2, which correspond to 40% difference. For the time delay, a pure sinusoidal variation suggests that the lag is directly proportional to the time period, which corresponds to a doubled time delay. The time delay ratio in the mixed case (8/16) in Figure 9 is 1.17-1.62. Therefore, the asymmetry is reasonable as the delay ratio is in between 1 and 2, in line with the 8/8 and 16/16 ratio.  The comparisons with the measured data shows that the simulations systematically underestimates the moisture buffering effect by 11-37%. The largest error for the case of clay is 37% (Sinusoidal 72/72) and lowest 11% (Sinusoidal 8/16). For gypsum the differences are 36% (Triangular 8/16) and 12% (Sinusoidal 8/16). Figure 10 shows a comparison between measured and simulated moisture uptake for gypsum and Triangular 12 h/12 h cycles. Figure 10. Comparison between measured (blue) and simulated (red) moisture uptake for Gypsum and triangular 12 h/12 h cycles.

Analysis of the Simulation Results
The difference between the measured and simulated results can depend on material properties, simplifications in the model, insufficient depth of samples and the assumption of the that the temperature in the material is following the boundary temperature.
Looking at the expected material dependence, from Equation (28) it was obtained: Here, b v = δ 0 v ξ is the moisture vapour effusivity. The formula shows that the slope of the sorption curve and the diffusion coefficient play an equal role for the moisture uptake. The factor f A calculated previously for all cases is not significant, as it appears to depend weakly on the RH boundary function. The same consideration applies to the square root of the integrated humidity at saturation over the time period. The previous calculation for gypsum showed that the measured value would fit if either of the parameters are increased by 66%. Considering the expression for the effusivity, it also means that a good fit happens when both material properties are increased by √ 1.66 = 1.29 i.e., 29% increase. Figure 11 shows the comparison between the experimental Triangular 12/12h case, when this increase of the material properties is applied in the simulation.
The error estimates in relation to effusivity value gives a maximum of 13.5% higher value, not fully support a 29% higher simulation value required for a good fit. By increasing the effusivity by 66% the peak value of the experimental moisture uptake was still higher than the simulated one. The measurement uncertainty is estimated to be less than 3% which would not explain the differences between the measurement and simulations. Another potential uncertainty lies in the surface resistance value. The surface resistance value is especially important for the diurnal variations since the equivalent surface resistance thickness and the periodic penetration depth is closer in magnitude. For gypsum it is round 0.5 mm. The previous calculation shows that for a zero surface resistance, the moisture buffering effect increases by 8%.

Comparison of Analytical and Numerical Results
In order to understand the impact of thicknesses that are close to the periodic penetration depth, numerical solutions for the equations (Equations (10), (11) and (14)) are presented. Only the case of gypsum was analysed. The Matlab pdepe-solver were used with a space resolution of around 0.1 mm.
From the results it can be seem that the semi-infinite analytical solution is verified by the numerical simulation and that it also works well for diurnal variations as long as the thickness is 20 mm or more (Table 10). For thinner material thicknesses the moisture buffering effect can be over or under estimated while the time lag becomes overestimated. It can also be observed that there is a maximum buffering effect for material thickness similar to the periodic penetration depth, which is result well known from the heat transfer area investigating effect of thermal mass on energy demand [24]. For 72/72 variations the semi-infinite analytical solution works well as long as the thickness is 40 mm or more (Table 11). For thinner material thicknesses the moisture buffering effect can be over or under estimated while the time lag becomes overestimated. The moisture buffering goes up with 14% for the 72/72 simulation with the thickness 20 mm instead of infinite thickness. This reduces the difference between experimental and simulation result.
The time lags is reduced between with 3.8 and 5.0 hours, when thickness is reduced down to 20 mm for infinity. The match between this simulation and the measurement becomes better for the top peak time delay but there is a larger lag for the bottom peak. Overall, it can be stated that either in the experimental and simulations temperature is a significant factor that delays the response of materials to buffer humidity. This is further justified in [8], in which the response of materials to constant and variable temperature were compared, showing that variable temperatures generate a significant lag compared to the case at constant temperature.

Conclusions
Clay and gypsum coatings were experimentally tested to investigate the moisture buffering capacity of the materials under simultaneous temperature and RH fluctuations. The dynamic sorption capacity of plasters were observed when environmental variations follow either sinusoidal and triangular variations and different interval for the humidification/de-humidification were applied (12/12 h, 8/16 h, 72/72 h). In this way, it was possible to better understand the influence of temperature and humidity on the moisture uptake of construction materials. The experimental results were compared with a novel simulation method that considers the impact of the surface resistance in the calculation of the material moisture uptake under variable environmental conditions. The model approach to consider materials as semi-infinite bodies was also verified with the analytical results.
The simulation model supported by the experimental results demonstrated clay and gypsum responded differently to the simultaneous temperature and RH fluctuations, due to the different material properties and pore structure. However, the environmental signal has a significant role in the moisture buffering capacity of materials. Both the RH and temperature signal (triangular and sinusoidal) and time intervals impact the way plasters adsorb moisture from the indoor.
Differences were found when different time frames were applied. When materials were subjected to daily cycles (12/12 h and 8/16 h), there were no significant differences between the sinusoidal and triangular test and also between the 12/12 h and 8/16 h in terms of peak to peak sorption capacity. However, when temperature and RH were slowly varied over six days, clay and gypsum stored and released more moisture, showing also differences between the sinusoidal and triangular test. The slower increase and decrease of the environmental conditions allowed the plasters to store more moisture.
The difference between the 72 h tests is linked to the different materials thickness and penetration depth. The penetration depth not only depends on the material property, but it also varies depending on the environmental conditions. For gypsum, which thickness was 20 mm, its moisture buffering performance increased in the 72 h test, due to the deeper penetration of moisture into the plaster, and due to the longer exposure to high and low environmental conditions. Another observation was the lag of both plasters when exposed to simultaneous temperature and RH variations. Simulations showed a delay of few hours in the daily variations and up to 20 h in the 72 h test. The lag is attributed to temperature variations that delays the response of materials. Moreover, the continuous variations of the environmental conditions did not allow the material to stabilise and to response quickly to the humidity variations.
In conclusion, explicit analytical formulas applied in the model were derived to calculate the moisture distribution inside a material that is exposed to cyclic variations in RH and temperature. The semi-infinite analytical solution was verified, confirming the formula can be used to estimate the uptake and release of moisture to an exposed material surface as long as the periodic penetration depth is less than half the thickness of the material.
The analytical method showed a great resemblance to the measured values, and together to the novel experimental approach of testing the sorption capacity of material will quantify more accurately the capacity of hygroscopic materials to moderate the indoor humidity and improve the hygrothermal comfort and health of people in buildings.