# Effect of Suction Cycles and Suction Gradients on the Water Retention Properties of a Hard Clay

Authors: Liufeng Chen; Hua Peng
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

The effect of suction cycles and suction gradients on a hard clay is investigated. The cylindrical samples of the hard clay are prepared to carry out the hydration and dehydration tests with different suction gradient and suction cycles. The results show that the suction gradient has little effect on the suction-water content relation, while the suction cycle has great effect on it, particularly the first cycle of hydration and dehydration. The apparent moisture diffusion coefficient of the hard clay has been identified by the use of a two-dimensional diffusion model. The moisture diffusion coefficient varies between 4.10−11 m2/s and 2.10−10 m2/s and it decreases during dehydration while the relative humidity is less than 85%. The results also show that the suction cycles play little effect on the moisture diffusion coefficient.


---

## Body

## 1. Introduction

In Belgium, a hard clay has been considered as an ideal geological formation for high-level radioactive nuclear waste disposal in Belgium [1, 2], because it has several favourable factors such as low permeability, plasticity, and high retention capacity for radionuclide. The hydromechanical behaviour of the host rock is one of the most important issues in the feasibility studies of the deep geological repository and many studies concerning coupled hydromechanical behaviour and unsaturated behaviour have been carried out [3–5]. Earlier studies show that the mechanical behaviour of clay strongly depends on water content [6, 7] and dehydration and hydration can induce a chemical-physical reaction and change the microstructure of clay and degrade the properties of clay [8–10]. The excavation stage and the open-drift stage of a nuclear waste repository could last a few decades during which the ventilation and temperature conditions will change [11]. The hydromechanical properties of clay could be changed due to cyclic moisture variation. The studies of Wang et al. in 2013 [10] show that different suction gradients play a notable effect on the microstructure of the COx argillaceous rock and larger suction gradients create more microcracks in the clayey matrix than smaller suction gradient during the same suction range. The microstructure change could affect the water retention properties of the material. The purpose of this study was to investigate the effect of suction cycles and suction gradients on the water retention properties of Boom clay.


## 2. Material and Methods

### 2.1. Material and Sample Preparation

In this study, the hydration and dehydration tests have been performed on Boom clay that is a plastic clay. Boom clay contains an average of 55% clay minerals including illite and smectite, 20%–25% quartz, 5%–10% feldspar, and the small content of pyrite (<5%) and calcite (<5%) [12]. The dry density of the material varies between 1.6 and 1.7 g/cm3 [3]. The total porosity of Boom clay is between 36% and 40% [1, 12] and the hydraulic conductivity is around 10−12 m/s. Young’s modulus of the Boom clay is around 300 MPa and the uniaxial compression strength is about 2 MPa.The cylindrical sample with a diameter of 36 mm and a length of 10 mm approximately was machined from the cores extracted from the Praclay gallery in theMol underground research laboratory at a depth of 223 m in Belgium. Special care was taken to assure that the sample was in saturated state. The measurements show that the tested samples had an average water content of 25% at the initial state and dry density of 1.63 g/cm3.


### 2.2. Experimental Method

In order to determine the water retention properties of Boom clay, the suction was imposed by controlling the relative humidity in an airtight container. The relative humidity can be imposed by the use of supersaturated saline solution, which creates a fixed vapour pressure in a sealed container. As shown by Delage et al. in 1998 [13], different values of relative humidity corresponding to the saline solution are listed and they firmly depend on the temperature. The tested samples were placed into the desiccators containing saturated saline solutions at their bottom. Different supersaturated saline solutions (LiCl, CaCl2·6H2O, K2CO3·2H2O, NaBr, NaCl, (NH4)2SO4, KCl, ZnSO4·7H2O, and CuSO4·5H2O) were chosen to cover a large range of relative humidity from 98% to 11%. The corresponding suction at 25°C varies from 2.76 to 301.51 MPa, following the Kelvin law [13]:(1)Pc=-RTMwVwln⁡PP0=-RTMwVwln⁡HR100.Pc is the suction, equal to the absolute difference between the air and water pressures; R is the constant of perfect gases (R=8.3143 J·mol·K−1); T is the absolute temperature (T=298 K at 25°C); Mw is the molar masse of water (Mw=18.016 g·mol−1); Vw is specific volume of water (Vw=0.99565 g·cm−3); P is the partial water vapour pressure; and P0 is the fact that, at state saturated, their ration is the relative humidity. Note that all the tests are carried out at a constant environment temperature.The saturated saline solutions used and the corresponding suctions/relative humidity at 20°C have been presented by Delage et al. in [13]. Note that all the tests are carried out at a constant environment temperature (25°C).


### 2.3. Test Program

To study the effects of suction gradients and suction cycles on water retention properties of Boom clay, four different dehydration and hydration paths were programmed as shown in Table1. In order to obtain the representative results, a series of samples are submitted to a same hydration or dehydration path. Path A was performed on the six samples (A1–A6) with different heights and it consisted of one cycle of dehydration and hydration with small suction gradient. Path B was performed on five samples (B1–B5) and it consists of two drying/wetting cycles with small suction gradient. Path C was performed on eight samples (C1–C8) and it is composed of two and a half cycles of hydration and dehydration with large suction gradient. Path D was performed on eight samples (D1–D8) and it consists of three drying/wetting cycles with larger suction gradient. Because of the low permeability of Boom clay, each moisture stage lasted more than one week to reach a stable state. Note that more than five samples were simultaneously tested following the same hydration and dehydration paths to get better confidence in the results.Table 1
Four different dehydration and hydration paths.
Type
Dehydration and hydration path
A
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
90
%
→
85
%
→
81
%
→
75
%
→
59
%
→
43
%
→
32.5
%
→
11
%
→
32.5
%
→
43
%
→
59
%
→
75
%
→
81
%
→
85
%
→
90
%
→
98
%
B
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
C
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
81
%
→
59
%
→
11
%
→
98
%
→
81
%
→
59
%
→
11
%
→
59
%
→
81
%
→
98
%
→
81
%
→
59
%
→
11
%.
D
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%.When the samples reached mass stabilization, the sample was weighted rapidly and then placed in the desiccator to subject to the next hydric loading. Once the whole hydration and dehydration cycles finished, the tested samples were dried in the oven at 105°C to determine the dry weight. The volume of the sample was measured by immersing the waxed samples into water.


## 2.1. Material and Sample Preparation

In this study, the hydration and dehydration tests have been performed on Boom clay that is a plastic clay. Boom clay contains an average of 55% clay minerals including illite and smectite, 20%–25% quartz, 5%–10% feldspar, and the small content of pyrite (<5%) and calcite (<5%) [12]. The dry density of the material varies between 1.6 and 1.7 g/cm3 [3]. The total porosity of Boom clay is between 36% and 40% [1, 12] and the hydraulic conductivity is around 10−12 m/s. Young’s modulus of the Boom clay is around 300 MPa and the uniaxial compression strength is about 2 MPa.The cylindrical sample with a diameter of 36 mm and a length of 10 mm approximately was machined from the cores extracted from the Praclay gallery in theMol underground research laboratory at a depth of 223 m in Belgium. Special care was taken to assure that the sample was in saturated state. The measurements show that the tested samples had an average water content of 25% at the initial state and dry density of 1.63 g/cm3.


## 2.2. Experimental Method

In order to determine the water retention properties of Boom clay, the suction was imposed by controlling the relative humidity in an airtight container. The relative humidity can be imposed by the use of supersaturated saline solution, which creates a fixed vapour pressure in a sealed container. As shown by Delage et al. in 1998 [13], different values of relative humidity corresponding to the saline solution are listed and they firmly depend on the temperature. The tested samples were placed into the desiccators containing saturated saline solutions at their bottom. Different supersaturated saline solutions (LiCl, CaCl2·6H2O, K2CO3·2H2O, NaBr, NaCl, (NH4)2SO4, KCl, ZnSO4·7H2O, and CuSO4·5H2O) were chosen to cover a large range of relative humidity from 98% to 11%. The corresponding suction at 25°C varies from 2.76 to 301.51 MPa, following the Kelvin law [13]:(1)Pc=-RTMwVwln⁡PP0=-RTMwVwln⁡HR100.Pc is the suction, equal to the absolute difference between the air and water pressures; R is the constant of perfect gases (R=8.3143 J·mol·K−1); T is the absolute temperature (T=298 K at 25°C); Mw is the molar masse of water (Mw=18.016 g·mol−1); Vw is specific volume of water (Vw=0.99565 g·cm−3); P is the partial water vapour pressure; and P0 is the fact that, at state saturated, their ration is the relative humidity. Note that all the tests are carried out at a constant environment temperature.The saturated saline solutions used and the corresponding suctions/relative humidity at 20°C have been presented by Delage et al. in [13]. Note that all the tests are carried out at a constant environment temperature (25°C).


## 2.3. Test Program

To study the effects of suction gradients and suction cycles on water retention properties of Boom clay, four different dehydration and hydration paths were programmed as shown in Table1. In order to obtain the representative results, a series of samples are submitted to a same hydration or dehydration path. Path A was performed on the six samples (A1–A6) with different heights and it consisted of one cycle of dehydration and hydration with small suction gradient. Path B was performed on five samples (B1–B5) and it consists of two drying/wetting cycles with small suction gradient. Path C was performed on eight samples (C1–C8) and it is composed of two and a half cycles of hydration and dehydration with large suction gradient. Path D was performed on eight samples (D1–D8) and it consists of three drying/wetting cycles with larger suction gradient. Because of the low permeability of Boom clay, each moisture stage lasted more than one week to reach a stable state. Note that more than five samples were simultaneously tested following the same hydration and dehydration paths to get better confidence in the results.Table 1
Four different dehydration and hydration paths.
Type
Dehydration and hydration path
A
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
90
%
→
85
%
→
81
%
→
75
%
→
59
%
→
43
%
→
32.5
%
→
11
%
→
32.5
%
→
43
%
→
59
%
→
75
%
→
81
%
→
85
%
→
90
%
→
98
%
B
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
C
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
81
%
→
59
%
→
11
%
→
98
%
→
81
%
→
59
%
→
11
%
→
59
%
→
81
%
→
98
%
→
81
%
→
59
%
→
11
%.
D
S
a
t
u
r
a
t
e
d
s
t
a
t
e
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%
→
75
%
→
32.5
%
→
75
%
→
98
%.When the samples reached mass stabilization, the sample was weighted rapidly and then placed in the desiccator to subject to the next hydric loading. Once the whole hydration and dehydration cycles finished, the tested samples were dried in the oven at 105°C to determine the dry weight. The volume of the sample was measured by immersing the waxed samples into water.


## 3. Experimental Results

### 3.1. Water Retention Curves

The successive samples weights(M) at the moisture equilibrium were measured and the corresponding water contents (w) were thus determined by comparing the measured weight (M) with the dried weight (Md) obtained at the end of the test: w=(M-Md)/Md. All the results concerning the water content show that the samples submitted to the same hydration and dehydration path have close water content change during each moisture phase as shown in Figure 1 which illustrates the water content evolution of the three samples (A1, A2, and A3) following path A. Therefore, the water content evolution observed of the Boom clay in centimetre should be representative to study the water retention curves of the material. We will present only the results of the samples A1, B1, C1, and D1 in order to focus on the effect of suction cycles and suction gradients on the water content properties.Figure 1
Hydration and dehydration path and the corresponding water content evolution of the samples A1, A2, and A3.The gravitational water content versus suction of the samples A1, B1, C1, and D1 is plotted in Figures2(a), 2(b), 2(c), and 2(d), respectively. Similar to the results of Bernier et al. in 1997 [3] and Li et al. in 2006 [1], the water retention curves show that the water content decreases nonlinearly with the increase of suction. The water content variation is relatively larger than that obtained by CIEMAT [3] and close to the results obtained by the researchers of Universitat Politècnica de Catalunya (UPC) (Figure 3). The water content decreases greatly during the first desaturation and it losses more than 10% when the suction increases to 29 MPa corresponding to a relative humidity of 81%. The large variation of the water content should be related to both the initial state and the pore size distribution of Boom clay. During the first dehydration, the initial state of the sample is fully saturated and the continuous water will be vapoured and it will induce the shrinkage of Boom clay and lead to a large variation of water content. The first cycle and the second cycle of dehydration and hydration of the four samples (A1, B1, C1, and D1) are, respectively, regrouped in Figures 4(a) and 4(b), which show that the water content of the samples at the same suction level is very close and there is little effect of the suction gradient on the water retention curves. Figures 2(b), 2(c), and 2(d) show also that there is large difference between the first hydric cycle and the second one. However, there is small difference between the second hydric cycle and the third hydric cycle. Moreover, the water content is nearly reversible after the second and third hydric cycle (Figure 2(d)). Therefore, suction cycles play an important role on the water retention curves, mainly the first dehydration and hydration cycle. The little effect of suction gradient on the water retention properties should be related to small microstructure change during hydration and dehydration which is different from the COx argillaceous rock. The effect of suction cycle on the water retention properties should be mainly related to the pore size distribution as Bernier et al. in 1997 [3] and Coussy in 2010 [15] indicated.Figure 2
The water retention curves of the samples A1, B1, C1, and D1.
(a)
(b)
(c)
(d)Figure 3
The comparison of the water retention curves of different studies.Figure 4
(a) Effect of suction gradient on the water retention curves during the first hydric cycle and (b) the second hydric cycle.
(a)
(b)


### 3.2. Moisture Diffusion Coefficient

The dehydration and hydration tests with different paths lasted more than six months. The long test duration should be firmly related to the transfer moisture properties. In fact, the moisture transfer process is complex and it depends on the relative humidity level. At high relative humidity, the capillary phase remains continuous and the moisture transfer is governed by the capillary forces; when relative humidity decreases, the liquid water and the gas (vapour and air) coexist and the moisture transfer is controlled by the exchange between evaporation and condensation; at low relative humidity, the pores are mainly occupied by gas and the moisture transfer is controlled by vapour diffusion [14]. In this study, the apparent moisture diffusion coefficient is used to characterize the transfer properties of Boom clay. Assuming the validity of the mass conservation law and Fick’s law [15], the moisture transfer of the cylindrical sample of Boom clay can be expressed using a macroscopic model in two dimensions to simulate the moisture diffusion in the radial direction and the vertical direction: (2)∂me∂t=De∂2me∂r2+1r∂me∂r+∂2me∂z2.With 0≤r<a, 0<z<l  for t>0, where me is the volumetric mass water content (kg·m−3) defined as me=ρeϕSr (kg·m−3); ρe is the water density; ϕ is porosity; Sr is the saturation degree; a is the radius of sample (m) and l is the length (m); and De is the apparent moisture diffusion coefficient (m2·s−1).While the initial conditions (porosity, geometric parameters) and boundary conditions (initial and final saturation) are given, the numerical solution of (2) can be obtained by means of the finite difference time domain method. The moisture equilibrium time only depends on De. The apparent moisture diffusion coefficient is determined by searching an optimal parameter De to minimize the difference between the experimental masse evolution and the numerical solution of (2). Figure 5(a) illustrates the experimental weight evolution of the sample D1 when the relative humidity decreased from 75% to 32% to compare with its numerical solution with a coefficient of diffusion of 1.4 10−10m2·s−1. Figure 5(b) illustrates the case of hydration. The numerical prediction data fits well with the experimental results. Most of the hydric phases are analysed and the moisture diffusion coefficient of the sample D1 varies from 4 10−11 m2·s−1 to 2 10−10m2·s−1. It is much smaller at the high relative humidity of 98% than at the lower humidity level as shown in Figures 6(a) and 6(b) which gives evidence also of the little effect of hydration and dehydration cycle on moisture diffusion coefficient. In fact, the moisture diffusion coefficient is firmly related to the water permeability, which decreases with the diminution of the relative humidity. The results of this study confirm it as shown in Figure 6, even if the apparent diffusion coefficient varies very little during hydration and dehydration when the relative humidity is less than 85%.Figure 5
(a) Experimental results and the numerical prediction of the weight evolution of the sample D1 during dehydration while the relative humidity (RH) decreases from 75% to 32% (b) and during hydration with the increase of RH from 32% to 75%.
(a)
(b)Figure 6
(a) The apparent moisture diffusion coefficient of the sample B1 (b) and the sample D1.
(a)
(b)


## 3.1. Water Retention Curves

The successive samples weights(M) at the moisture equilibrium were measured and the corresponding water contents (w) were thus determined by comparing the measured weight (M) with the dried weight (Md) obtained at the end of the test: w=(M-Md)/Md. All the results concerning the water content show that the samples submitted to the same hydration and dehydration path have close water content change during each moisture phase as shown in Figure 1 which illustrates the water content evolution of the three samples (A1, A2, and A3) following path A. Therefore, the water content evolution observed of the Boom clay in centimetre should be representative to study the water retention curves of the material. We will present only the results of the samples A1, B1, C1, and D1 in order to focus on the effect of suction cycles and suction gradients on the water content properties.Figure 1
Hydration and dehydration path and the corresponding water content evolution of the samples A1, A2, and A3.The gravitational water content versus suction of the samples A1, B1, C1, and D1 is plotted in Figures2(a), 2(b), 2(c), and 2(d), respectively. Similar to the results of Bernier et al. in 1997 [3] and Li et al. in 2006 [1], the water retention curves show that the water content decreases nonlinearly with the increase of suction. The water content variation is relatively larger than that obtained by CIEMAT [3] and close to the results obtained by the researchers of Universitat Politècnica de Catalunya (UPC) (Figure 3). The water content decreases greatly during the first desaturation and it losses more than 10% when the suction increases to 29 MPa corresponding to a relative humidity of 81%. The large variation of the water content should be related to both the initial state and the pore size distribution of Boom clay. During the first dehydration, the initial state of the sample is fully saturated and the continuous water will be vapoured and it will induce the shrinkage of Boom clay and lead to a large variation of water content. The first cycle and the second cycle of dehydration and hydration of the four samples (A1, B1, C1, and D1) are, respectively, regrouped in Figures 4(a) and 4(b), which show that the water content of the samples at the same suction level is very close and there is little effect of the suction gradient on the water retention curves. Figures 2(b), 2(c), and 2(d) show also that there is large difference between the first hydric cycle and the second one. However, there is small difference between the second hydric cycle and the third hydric cycle. Moreover, the water content is nearly reversible after the second and third hydric cycle (Figure 2(d)). Therefore, suction cycles play an important role on the water retention curves, mainly the first dehydration and hydration cycle. The little effect of suction gradient on the water retention properties should be related to small microstructure change during hydration and dehydration which is different from the COx argillaceous rock. The effect of suction cycle on the water retention properties should be mainly related to the pore size distribution as Bernier et al. in 1997 [3] and Coussy in 2010 [15] indicated.Figure 2
The water retention curves of the samples A1, B1, C1, and D1.
(a)
(b)
(c)
(d)Figure 3
The comparison of the water retention curves of different studies.Figure 4
(a) Effect of suction gradient on the water retention curves during the first hydric cycle and (b) the second hydric cycle.
(a)
(b)


## 3.2. Moisture Diffusion Coefficient

The dehydration and hydration tests with different paths lasted more than six months. The long test duration should be firmly related to the transfer moisture properties. In fact, the moisture transfer process is complex and it depends on the relative humidity level. At high relative humidity, the capillary phase remains continuous and the moisture transfer is governed by the capillary forces; when relative humidity decreases, the liquid water and the gas (vapour and air) coexist and the moisture transfer is controlled by the exchange between evaporation and condensation; at low relative humidity, the pores are mainly occupied by gas and the moisture transfer is controlled by vapour diffusion [14]. In this study, the apparent moisture diffusion coefficient is used to characterize the transfer properties of Boom clay. Assuming the validity of the mass conservation law and Fick’s law [15], the moisture transfer of the cylindrical sample of Boom clay can be expressed using a macroscopic model in two dimensions to simulate the moisture diffusion in the radial direction and the vertical direction: (2)∂me∂t=De∂2me∂r2+1r∂me∂r+∂2me∂z2.With 0≤r<a, 0<z<l  for t>0, where me is the volumetric mass water content (kg·m−3) defined as me=ρeϕSr (kg·m−3); ρe is the water density; ϕ is porosity; Sr is the saturation degree; a is the radius of sample (m) and l is the length (m); and De is the apparent moisture diffusion coefficient (m2·s−1).While the initial conditions (porosity, geometric parameters) and boundary conditions (initial and final saturation) are given, the numerical solution of (2) can be obtained by means of the finite difference time domain method. The moisture equilibrium time only depends on De. The apparent moisture diffusion coefficient is determined by searching an optimal parameter De to minimize the difference between the experimental masse evolution and the numerical solution of (2). Figure 5(a) illustrates the experimental weight evolution of the sample D1 when the relative humidity decreased from 75% to 32% to compare with its numerical solution with a coefficient of diffusion of 1.4 10−10m2·s−1. Figure 5(b) illustrates the case of hydration. The numerical prediction data fits well with the experimental results. Most of the hydric phases are analysed and the moisture diffusion coefficient of the sample D1 varies from 4 10−11 m2·s−1 to 2 10−10m2·s−1. It is much smaller at the high relative humidity of 98% than at the lower humidity level as shown in Figures 6(a) and 6(b) which gives evidence also of the little effect of hydration and dehydration cycle on moisture diffusion coefficient. In fact, the moisture diffusion coefficient is firmly related to the water permeability, which decreases with the diminution of the relative humidity. The results of this study confirm it as shown in Figure 6, even if the apparent diffusion coefficient varies very little during hydration and dehydration when the relative humidity is less than 85%.Figure 5
(a) Experimental results and the numerical prediction of the weight evolution of the sample D1 during dehydration while the relative humidity (RH) decreases from 75% to 32% (b) and during hydration with the increase of RH from 32% to 75%.
(a)
(b)Figure 6
(a) The apparent moisture diffusion coefficient of the sample B1 (b) and the sample D1.
(a)
(b)


## 4. Conclusions

In order to study the effect of suction cycles and suction gradients on the water retention properties of the hard clay, a series of samples are prepared to carry out dehydration and hydration tests by means of suction controlled technique. Four hydric paths with different suction cycles and suction gradients are programmed and each hydric test has lasted more than six months. The results concerning water content show that the samples following the same hydric path have a very close water content variation during the same suction change and suction gradient has little effect on the suction water content relation, while the suction cycle has great effect on it, particularly the first cycle of hydration and dehydration will change much the suction-water content relation.Because of the very low permeability, each hydric phase during dehydration and hydration lasted more than one week. The moisture transfer properties of Boom clay were quantified by the use of a two-dimensional diffusion model. The apparent moisture diffusion coefficient of Boom clay varies between 4 10−11m2/s and 2 10−10m2/s. It decreases during dehydration while the relative humidity is less than 85%. The moisture diffusion coefficient is relatively small when the relative humidity is close to 98% and it is about 4 10−11m2/s. The results show also that the suction cycles play little effect on the moisture diffusion coefficient.



---

*Source: 313682-2015-12-08.xml*