Effect of temperature on the sorption curves of earthen materials

Unfired earth is a sustainable construction material with low embodied energy, but its development requires a better evaluation of its moisture–thermal buffering abilities and its mechanical behavior. Both of them are known to strongly depend on the amount of water contained in its porous network and its evolution with external conditions (temperature, humidity), which can be assessed through several sorption–desorption curves at different temperature. However, the direct measurement of these curves is particularly time consuming (up to 2 month per curve) and thus, indirect means of their determination appear of main importance for evident time saving and economical reasons. In this context, this paper focuses on the prediction of the evolution of sorption curves with temperature on earth plasters and compacted earth samples. For that purpose, two methods are proposed. The first one is an adaptation of the isosteric method, which gives the variation of relative humidity with temperature at constant water content. The second one, based on the liquid–gas interface equilibrium, gives the variation of water content with temperature at constant relative humidity. These two methods lead to quite consistent and complementary results. It underlines their capability to predict the sorption curves of the tested materials at several temperatures from the sole knowledge of one sorption curve at a given temperature. Finally, these predictions are used to scan the range of temperature variation within which the evolution of water content with temperature at constant humidity could be neglected or should be taken into account.

controlling properties [3][4][5]. At building scale, it was shown that the use of hygroscopic materials leads to a significant reduction of moisture variation amplitudes, which thus induces energy savings on ventilation and heating [6]. This is due to the microstructure of the earth, which enables hydric exchanges between the environment and water molecules on the pore surfaces through condensation/evaporation and sorption/desorption phenomena. This affinity with water molecules also significantly impacts mechanical behavior of earth materials. For example, the decrease in strength with moisture, which is well known in soil mechanics, has been demonstrated for rammed earth [7][8][9]. Consequently, it appears that the liquid water content of a rammed earth wall is a key parameter in order to understand the behavior and the strength of this material.
At the macroscopic scale, sorption isotherm characterizes the water uptake with increasing ambient humidity at a constant temperature, while the desorption isotherms characterizes water expulsion with decreasing ambient humidity at a constant temperature. These curves are strongly non-linear. Furthermore, a hysteresis can be observed between sorption and desorption. This phenomenon is quite common and has been widely studied by many authors for a large variety of materials like soil, wood and hemp concrete [10][11][12][13]. In addition, sorption-desorption curves are known to vary with temperature. The general tendency is a small decrease of the water content at constant relative humidity as temperature rise.
It follows that a proper determination of the evolution of the water content of an hygroscopic material with the external conditions (temperature, humidity) would require to perform several sorptiondesorption curves at different temperature. Knowing that the measurement of a single sorption-desorption loop on earthen material lasts from two weeks to more than two months in function on the experimental protocol, this option appears to be not acceptable for most of earthen construction projects for quite evident duration and economical reasons. Thus, any procedure which would allow to accurately assess the sorption and desorption curves from indirect measurement is of main importance.
In this context, this paper aims at proposing a method to predict the evolution of sorption curves with temperature. This effect has already been studied for a quite large range of porous materials including concrete [14], bio-based [15] and even clayey [16] materials. These studies are mainly based on the chemical equilibrium between adsorbed water and water vapor, which allows to link this temperature dependence of the sorption curves to the isosteric heat. However, since no clear data exists on the evolution of isosteric heat with temperature and relative humidity (or water content), these methods do not manage to predict the evolution of sorption-desorption curves with temperature. To overcome this problem, this paper proposes two different methods. The first one is based on the analysis of the heat of phase change in order to predict the isosteric heat evolutions. The second one is based on the link between morphology of the porous network and the water content. These two methods were applied on two kind of earthen materials: compacted earth blocks and samples of earth plasters and the obtained results were compared to experimental data which were obtained at 10 and 40°C. At end, this approach was used in order to discuss on the impact of temperature dependency of sorption curves on hygrothermal modeling.
Finally, it is important to note that a good understanding of the earthen constructions requires taking into account their large variability, which is due to the inner nature of the material (depending on the soil location and history) and on the great numbers of construction technics (compacted earth, adobe, cob, rammed earth, plaster, extruded bricks [17]…). It is not possible, however, to explore all the possibilities at the same time. In this paper the choice of compacted earth and earth plaster has been made since it allows to encompass the behavior of a quite important number of systems. It follows that the results presented in this paper cannot be directly generalized to all earthen construction techniques.

Material
Two kind of earth samples are tested in this paper: compacted earths and earth plasters. The raw materials used to realize the compacted earth samples come from existing centenarian rammed-earth constructions located in the Rhône-Alpes region in southeastern France: either STR from Saint-Trivier de Courtes or CRA from Cras-sur-Reyssouze. This choice ensures that the earths are suitable to construct earthen buildings. The complete characterization of these two earths is given in [8]. To sum up, the plasticity index and the Methylene Blue Value of CRA are almost twice as high as those for STR, while their grading curves are similar. These two values allow us to have a first insight on the nature of the clays within the material: the clayey proportion of STR should mainly be composed of Kaolinite and Illite, while much more active clays are expected in CRA. Having said that, the X-ray diffraction analysis of these two earths shows that their clayey proportions are both mainly composed by Illite.
The first step of the compacted samples fabrication protocol consists in crushing the rammed earth pieces of wall (which were taken directly from the construction sites) into small particles sieved through a 2-mm sieve to obtain a homogenous material. An homemade double-compaction process is then used to realize cylindrical samples of diameter 3.5 and 7 cm in height with a loading charge of 4 MPa. This charge is in the range of the one applied for the realization of compacted earth blocks with a manual press. The water content of fabrication was fixed in order to reach the higher dry density. This fabrication procedure ensures a global homogeneity of the resulting samples and a good repeatability of their main characteristics (density, compressive strength, deformability, vapor resistance factor, sorption-desorption curves). The compressed samples made with the earth CRA were referenced as C1, while the compressed samples made with the earth STR were referenced as C2.
Two earth plaster formulations, denoted by P1 and P2, were studied. The plasters P1 were composed by 25 wt% of Kaolinite (clay with a low specific surface compared to other clays; yet with a high sorption capacity compared to most minerals, whose particle size distribution is 43%\2 lm and 95%\80 lm), 73 wt% of sand (particle size lower than 2 mm), and 2 wt% of straw particles (size between 30 and 50 mm). The plasters P2 were composed by 31 wt% of Ascal10 (commercial fine calcareous-clay material whose particle size distribution is: 14%\2 lm and 97%\80 lm), 67 wt% of the same sand (particle size lower than 2 mm), and 2 wt% of straw particles (size between 10 and 30 mm). Each plaster was mixed with water (30 wt% for P1 and 24 wt% for P2) before being cast in a specific formwork (50 Â 50 cm Â 2 cm).
All the samples (either compacted earths or earth plasters) were dried in a controlled environment in a conditioning room at 20 ±2°C and 55 ± 5% of relative humidity. Periodical measurements were realized to follow the drying stage of the samples, and they were not tested before a stable mass was reached. A summary of samples composition is given in the Table 1.

Determination of the sorption isotherms
Several methods exist to estimate the isothermal sorption and desorption curves, but the two most widely used are the desiccator and dynamic gravimetric vapor sorption methods.
The dynamic gravimetric sorption method, commonly called the DVS (dynamic vapor sorption) method, consists in measuring uptake and loss of moisture by flowing a carrier gas at a specified relative humidity (or partial pressure) over a small sample (from several milligrams to several grams depending on the device used) suspended from the weighing mechanism of an ultrasensitive recording microbalance. Variations in the gas's relative humidity are automatically calculated by the device when the target condition in mass stability is reached. This method was used to draw the sorption-desorption curves at 23 and 30°C for all the formulations with the DVS Intrinsic Ó from Surface Measurement Systems. The results are reported in the Fig. 1. Only the sorption curve is drawn to ease the reading. Whatever the material, a small but visible decrease in water content with temperature was observed. It is however important to keep in mind that the difference of temperature between the curves were quite limited (only 7°C).
The other method used in this paper to assess the evolution of water content with relative humidity and temperature is the desiccator method, which is precisely described in the international standard ISO 12751. The sorption stage consists in successively putting a previously dried sample in several environments of increasing relative humidity and constant temperature. The sample is periodically weighed and it stays within a given environment until constant mass. The desorption stage consists in successively putting a sample previously equilibrated at 95% (at least) in several environments of decreasing relative humidity until mass is constant and at constant temperature.
The relative humidity within the desiccators is fixed by equilibrium with saturated saline solutions through the activity of its solvent (that is liquid water). This latter tends to decrease with salt concentration, but it is also impacted by the nature of the salt, the temperature, etc. The value of the equilibrium relative humidity at saturation for a dozen salts of different solubility levels and at several temperatures is given in the standard. However, this value remains quite theoretical, and it can be modified by a great number of factors like impurities in salts or water, nonsaturation of the saline solution at the liquid-gas interface… To avoid this problem, the effective relative humidity in the desiccators near the samples have been measured with Rotronic HygroLog HL-NT sensors. The salts which were used in this study as well as the target and measured relative humidity within the desiccators are summarized in Table 2.
Using this protocol, a sorption-desorption loop can be made in approximatively 2 months for earthen materials, while only a period of 2 weeks is necessary if the DVS method is used. On the other hand, the desiccator method can test several specimens at the same time, and it is the only way to test materials with high levels of heterogeneity because of the small dimensions of DVS samples.
Whatever method used, the isothermal sorption curves should be intrinsic to the tested material. However, direct comparisons between the curves obtained by the DVS or the desiccator methods is not possible since, as it is already underlined by Labat et al. [18], the reference dry mass used to calculated the water content is intrinsically different for the two   methods: the dry mass considered by the DVS is obtained by drying the sample through a flow of dry air, while the dry mass for the desiccator method is commonly estimated from the drying at 105°C. One option to compare the two methods should be to dry the DVS sample in an oven at 105°C after (or before) the sorption-desorption test. If this dry mass, denoted by m 105 C d , is taken as the reference, the water content given by the DVS, denoted by w DVS , must be corrected as follows: where m 23 C d is the dry mass obtained at 23°C under the dry atmosphere of the DVS. However, the oven drying stage is particularly complicated to perform. Indeed, since the sample mass for the DVS is extremely low (around 1 g), any dust supply or release may significantly change the results. To avoid this complexity, another option is to modify the desiccator experimental protocol by adopting a drying process similar to the DVS's one. For that purpose, a vacuum drying at 23°C within a dessicator with silica gel, insuring a relative humidity lower than 5% during at least the whole drying period (approx. 2 weeks), was used in this study.
The isothermal sorption curves of compacted earth C1 and earth plaster P1 at 23°C for the two methods (either DVS or dessicator) and the two kinds of reference dry mass (either m 105 C d or m 23 C d ) are showed in the Fig. 2. This comparison gives some confidence on the consistency between DVS and dessicator methods if the same protocol is used to estimate the dry mass. In the following of this study, for convenient purpose, no correction was applied on the DVS results and the dry mass of the desiccator method, which is also used to calculate the dry density, was determined using the vacuum drying procedure at 23°C.
3 Sorption curve and enthalpy of phase change

Integral and differential enthalpy of phase change
The natural definition of the specific enthalpy of desorption is the difference between the specific enthalpies of the vapor and of the adsorbed water, denoted by Dh v . If thermodynamical equilibrium is assumed between adsorbed water and water vapor, it is linked to the specific entropy of these two phases (denoted by s L for adsorbed water and by s v for water vapor) through the equation (proof in ''Appendix''): where T is the temperature at which the phase change occurs. Using the expression of s L and s v reported in the ''Appendix'', the Eq. (2) becomes: where L e ¼ 2:26 MJ=kg is the specific latent heat of evaporation of free water at boiling temperature, denoted by T e ¼ 373 K, R ¼ 8:31 J=mol=K is the perfect gas constant, M H 2 O ¼ 18 g=mol is the molar mass of water, while C v ¼ 1616 J=kg=K and C L ¼ 4185 J=kg=K are respectively the specific heat of vapor and adsorbed water, which are assumed to be constant. Finally, u is the relative humidity, which is defined as the ratio between the partial pressure of vapor at which the phase change occurs (p v ) and the partial pressure of vapor at saturation (p sat v ). The last term of (3), namely ðÀRT=M H 2 O Þ ln u leads to an increase in Dh v when the equilibrium relative humidity decreases. As is depicted by the Eq. (2) this stems from the quite logical increase of vapor entropy when the concentration of vapor molecules within the air phase decreases. Note that u in this expression is the relative humidity at the local sorption/desorption front (that is the interface between adsorbed water and water vapor). For free pure water, the relative humidity at the evaporation/condensation front is always equal to 1, and thus this last term vanishes. However, it is almost always neglected even for confined adsorbed in-pore water (for which u\1 at the liquid-gas interface), while as already discussed by Soudani et al. [5], the variations of Dh v with T and u appear to be on the same order of magnitude. Let us eventually notice that the expression (3) is based on the assumption that the entropy of adsorbed water only depends on temperature. This assumption will surely be wrong at very low saturation where strong interactions occur between adsorbed water and pore walls. In such case, an additional corrective term should be added, but its quantification is out of the scope of this paper.
Another way to assess the enthalpy of phase change is to use the differential specific internal energy and entropy of the adsorbed phase, respectively denoted by _ u L and _ s L , and defined as where w ¼ m L =m d , with m L the mass of adsorbed and/ or liquid water and m d the dry mass, is the gravimetric water content, while U L and S L ¼ q d ws L are the internal energy and the entropy of the adsorbed phase per unit of initial total volume of the material and q d is the dry density. The combined use of these notations with the chemical equilibrium condition provides: By analogy with the Eq. (2), D _ h ¼ h v À _ u L is defined as the differential enthalpy of desorption. It is also commonly called isosteric heat. This form of the enthalpy of phase change is particularly interesting because it satisfies the following Clapeyron-like equation (proof in ''Appendix''): It follows that D _ h can be experimentally determined from the analysis of the sorption-desorption curve variations with temperature. If the two temperatures remain sufficiently close, it is possible to ignore the variation of D _ h with temperature, and the integration of the Eq. (31), which is reported in the ''Appendix'', provides: where x is the water content, which is held constant, T 1 and T 2 are the two temperatures at which the sorption-desorption curves are drawn, T m ¼ ðT 1 þ T 2 Þ=2 is the average temperature, while p v ðT 1 ; wÞ and p v ðT 2 ; wÞ are the vapor pressure at equilibrium with the water content w for the temperatures T 1 and T 2 , respectively. The combination of equations (2) and (5) leads to the following relation between Dh v and D _ h: Under the assumption that s L only depends on the temperature, the term os L =ow becomes null and the isosteric heat (D _ h) can henceforth be assumed to be equal to the integral enthalpy of the phase change (Dh v ). A way to check this assumption is to compare the values of Dh v given by the Eq. (2) and the differential enthalpy estimated from the isosteric method on sorption curves at 23 and 30°C for each tested materials. The results are reported in the Fig. 3. The experimental points appear to be quite scattered, in particular at low humidity. This may be due to measurement uncertainties, since an error of 3% in the measurement of the relative humidity leads to a theoretical relative error of approximatively 20% in D _ h at 10%rh. This error becomes lower than 1% when the relative humidity is higher than 30%. Anyway, this graph underlines that the global tendency of calculated values of Dh v and the estimated ones for D _ h is the same, that is an increase of the enthalpy of desorption when the relative humidity decreases. Thus, even if they can not be considered as hard proofs, these results give some confidence on the validity of this approach for the materials considered in this study.

Isosteric method: variation of relative humidity with temperature at constant water content
Alternatively, the combined use of the Eq. (7) and the assumption Dh v ¼ D _ h can be used to estimate the evolution of sorption curve with temperature: Throughout the rest of the paper, this way to estimate the variation of sorption curve with temperature is called the ''isosteric method''. The difference between the sorption curves at 23 and 30°C is not sufficient in order to scan the accuracy of this method. In consequence, the comparison is made between sorption curves of the same samples at 10 and 40°C. Due to technical limitation of the DVS which is used in this study, the sorption curves are obtained with the desiccator method with the salts reported in the Table 2. In addition, only the compacted earth C1 and the earth plaster P1 have been tested.
Calculations are made using the sorption curve at 23°C obtained with the DVS as a reference. In other words, the Eq. 9 is used with T 1 ¼ 296 K while T 2 ¼ 283 K to build the sorption curve at 10°C, and T 2 ¼ 313 K for the sorption at 40°C. The results, are shown in Fig. 4. In order to ease the reading of the results, the water content was also drawn in function of vapor pressure instead of the relative humidity in Fig. 5.
Whatever the tested material, results underlines a good correlation between the predicted curves and the experimental data.

Theoretical background
One other way to assess the evolution of sorption curves is the analysis of the mechanical equilibrium of the liquid-gas interface, which can be depicted by the Young-Laplace law: where P atm is the atmospheric pressure, P L is the liquid pressure, j L;G is the average curvature of the liquidgas interface and c L;G is the apparent surface tension between the liquid water and the in-pore gas. It is quite difficult to assess this latter since its value may be modified by the interactions between the water molecules and the pore walls. In this study, the assumption was made that the tabulated expression determined for pure unconfined water, that is c L;G ¼ 72:74 À 0:1513ðT À 293:13Þ in [mN/m] [19], can be used.
If the deformation of in-pore space (due to mechanical loading or thermal dilation for example) and hysteresis phenomena are ignored, j L;G can directly be linked to the morphology of the porous network. This link has been already widely used in order to estimate the pore size distribution of porous material from their sorption/desorption curves [20,21]. Under this strong assumption, which should however be checked, in particular for earthen materials with a high content in vegetal fibres (like straw or hemp), j L;G can be written in the form: where F À1 is a monotonous bijective decreasing function of the water content, denoted by w, which satisfies F À1 ðw sat Þ ¼ 0, with w sat being the water content at saturation. Alternatively, the chemical equilibrium between the adsorbed water and its vapor leads to the celebrated Kelvin's law: where q L is the density of water, R is the perfect gas constant, M H 2 O is the molar mass of water, while T and u are the temperature and the relative humidity at the local sorption/desorption front. The combined use of (12), (10), and (11) allows to write x in the form: where F is the inverse function of F À1 , and The existence of this monotonous bijective function F between the water content and X have been checked for the earth plaster P1 and the compacted earth C1 using the DVS's sorption curves at 23 and 30°C as well as the desiccator's sorption curves at 10, 23 and 40°C. The results are reported in the Fig. 6. They lead to quite consistent results, which tends to validate the existence of this function, at least for the tested materials.
Using this formalism, and assuming that the surface tension only depends on temperature, the variation of the water content with temperature at constant hygrometry can be linked to the slope of the isothermal sorption curve through the equation (proof in ''Appendix''): is the slope of the sorption curve at the temperature T for the relative humidity u.

Interface method: Variation of water content with temperature at constant relative humidity
Similarly to what it is done in the previous section, the Eq. (15) is used to estimate the variation of water content at constant humidity between 10 and 40°C for the compacted earth C1 and the plaster P1: Throughout the rest of the paper, this way to estimate the variation of sorption curve with temperature is called the ''interface method''. The results, using the sorption curve at 23°C obtained with the DVS as the reference curve, are showed in Fig. 5. They underline a good correlation between the experimental data and calculated values for earth plaster P1 and compacted earth C1 for the sorption curves at 10 and 40°C.

Discussion
Both methods (either isosteric or interface) seem to quite well predict the variation of the sorption curve with temperature of the earthen materials tested in this study. However, information which is given by these two methods is intrinsically different. Isosteric method provides the variation of equilibrium relative humidity with temperature at constant water content while interface method gives the variation of water content with temperature at constant relative humidity.
Hygroscopic and hygrothermal models are based on the mass conservation of the in-pore water (both adsorbed, liquid and vapor). Assuming the equilibrium between all the in-pore water phases, and if the mass variation of vapor is strongly negligible towards the mass variation of liquid and adsorbed water, it leads to: where g L is the mass flow vector of liquid water, which depends on liquid pressure gradient and gravity through the Darcy's law and g v is the mass flow vector of vapor, which is composed by a diffusive term and an advective term, and which depend on gradients of vapor and total gas pressures. The developed equation of this relation can be found in a quite large amount of references, see for example [5,22,23]. Since the water content is a function of temperature and relative humidity, the first term of the Eq. (18), can be written as: It follows that, to solve the Eq. (18), the information provided by the interface method is more convenient since it directly gives the term ow=oT, and the combination between (15) and (19) provides: Fig. 6 Relation between the water content and X for the earth plaster P1 and the compacted earth C1 from the sorption curves obtained at 10, 23 and 40°C What's more, the knowledge of the function F ðXÞ, which can be experimentally obtained by plotting x as a function of X ¼ ðÀq L RTÞ=ðc L;G M H 2 O uÞ, would allow to assess the expression of n for each relative humidity and temperature through the Eq. (32) in ''Appendix''. In consequence, under the assumptions and conditions of this study, only a single sorption curve is required to assess the evolution of water content with both temperature and relative humidity. However, the interface method is based on the empirical knowledge of the c L;G ðTÞ function. Thus, for each new material, the existence of a function c L;G ðTÞ which allow to reach a single w ¼ FðXÞ curve must be verified. Along the same line, the assumption that h used in the isosteric method must be verified for each new tested material, especially when strong interactions may occurs between the adsorbed layer and the solid particles that forms the material.
Having said that, if the interface method is validated, the Eq. (20) can be used in order to estimate the impact of the temperature variation of the sorption curve on hygroscopic (and thus hygrothermal) calculations. For that purpose, the data from the Handbook of Chemistry and Physics [19], d ln c L;G =dT % À2:6 Â 10 À3 for temperature ranging from 0 and 50°C is considered. This gives values of u ln u 1 T À d ln c L;G dT which remain at most on the order of 2 Â 10 À3 K À1 (cf. Fig. 7). As a consequence, a 5°C variation in temperature appears to have the same impact on the water content as a 1% variation in relative humidity, which is negligible. On the other side, for temperature variations higher than 25°C, the variation of water content with temperature at constant humidity can become on the same order than a 5% variation in relative humidity, which may become significant.

Conclusion
This paper presents a study on the impact of temperature on the sorption curves of earth plasters and compacted earth samples. For that purpose, two means of measurement are used, either the DVS (dynamic gravimetric sorption method) or the desiccator and saline solutions method (defined in the international standard ISO 12571 [24]). A good consistency was found between these two methods providing that the sample dry masses were determined using a similar protocol.
The prediction of sorption curves at several temperatures was made following two approaches. The first one, denoted by the ''Isosteric method'' was based on the well known Clapeyron relation. For that purpose, assuming that entropy of adsorbed water only depends on temperature, a simple expression of the latent heat of desorption was proposed. This method provides the variation of the equilibrium relative humidity with temperature at constant water content.
The second one, denoted by the ''interface method'', gives variation of water content with temperature at constant relative humidity. It relies on the assumption of existence of a bijective function between the water content and the curvature of liquidgas interface, which is in turn bijectively linked to the porous network geometry.
The sorption curves estimated with these two methods were quite similar and were consistent with the experimental data in the studied range of temperature (between 10 and 40°C). This first result is quite interesting since it underlines that it is possible to predict the evolution of water content curve with temperature from the sole knowledge of one isothermal sorption curve.
Finally, using the interface method, the comparative impact assessment of temperature and relative humidity on water content have been realized. It follows that if temperature variations are below 25°C, the induced variation in water content would be, at maximum, in the same range than the one induced by a 5% variation in relative humidity. Thus even if the choice of taking into account the temperature humidity and temperature dependency of sorption curves will depend on the required precision, the assumption that w ¼ wðuÞ seems quite acceptable for moderate variations of temperature. Therefore, this study provides a simple tool which can assist in making that decision. To conclude, even if the developments of this paper focuses on earthen materials submitted to quite moderate variations of temperature (realistic for the indoor atmosphere), the methodology which is developed may be applied to other materials and/or conditions. In that view, it may be interesting to analyze materials for which the sorption curve exhibits stronger dependency in temperature and/or for larger ranges of temperatures.
where s L and s v are the specific entropies of adsorbed water and vapor, respectively, expressed as: The superscript 0 is used to denote the reference state, defined by P 0 G ¼ P 0 L ¼ P atm and T 0 ¼ T e ¼ 373:15 K, where T e is the boiling temperature at the atmospheric pressure. In this condition, vapor pressure at saturation is also equal to atmospheric pressure: while s 0 v À s 0 L , which is the difference in entropy, at P atm and T e between pure vapor and liquid water, is directly linked to the latent heat of evaporation at P atm and T e of water, denoted by L e , though the equation: The natural definition of the specific enthalpy of desorption is the difference between the specific enthalpies of the vapor and of the adsorbed water, denoted by Dh v . Since at equilibrium the specific free enthalpies of liquid water and its vapor are equal (g L ¼ g v ), the combined use of (21) and (22a-22b) makes it possible to write Dh v ¼ h v À h L as: At first order, the dependency partial vapor pressure at saturation can be depicted by the relation (see for example [5,25]): The combined use of this expression with the equations (26), (23), and (24) provides: where u ¼ p v =p sat v is the relative humidity.

Differential enthalpy
The use of the equation of state g L ¼ oF L =oðq d wÞ, where F L is the free energy of the adsorbed water (per unit of initial total volume), while accounting for the relation F L ¼ U L À TS L , where U L is the internal energy of the adsorbed water (per unit of initial total volume), leads to: which, when combined with the chemical equilibrium condition (g L ¼ g v ), finally provides: Clapeyron's equations An extensive use of the equations that define D _ h provides: The derivation of (29) with respect of the temperature reciprocal, at constant water content, while using the Eq. (30) eventually provides the following Clausius-Clapeyron-like equation: Water content and pore size distribution The partial derivation of the Eq. (13) with respect to the relative humidity (at constant temperature) provides: where Along the same line, the partial derivation of (13) with respect to the temperature (at constant relative humidity) provides: The combination of (32) and (34) eventually provides the equation: