Evaluation of actual and estimated hydraulic conductivity of sands with different gradation and shape

Abstract

Hydraulic conductivities of sands with different gradation and grain shape were estimated experimentally at a relative density (D_r) of about 40 % and a 22 ± 2 °C of constant temperature. Narli Sand (NS) with 0.67 of sphericity (S) and 0.72 of roundness (R), and Crushed Stone Sand (CSS) with 0.55 of S and 0.15 of R values were artificially graded into sixteen different grain-size fractions (4.75–2, 2–1.18, 1.18–0.6, 0.6–0.425, 0.425–0.3, 0.3–0.075, 4.75–0.075, 2–0.075, 1.18–0.075, 0.6–0.075, 0.425–0.075, 4.75–0.6, 2–0.6, 4.75-0.425, 2–0.425, 1.18–0.425 mm). Hydraulic conductivities of the NS estimated by use of constant head test ranged from 1.61 to 0.01 cm/s, whilst those of the CSS estimated by the same test ranged from 2.45 to 0.012 cm/s. It was observed that the hydraulic conductivity values of the NS are lower than those of the CSS samples, which is likely to be the result of differences in shape, particularly in R values. The results clearly demonstrated that the hydraulic conductivity can be significantly influenced by grading characteristics (d₁₀, d₂₀, d₃₀, d₅₀, d₆₀, c_u, c_c, n, I_o). Furthermore, comparisons between results obtained in the present study and hydraulic conductivity estimated with other formulas available in the literature were made. The comparisons indicated that the best estimation of hydraulic conductivity changes based on the gradation and shape properties of the sands tested.

Background

Hydraulic conductivity, which represents the ability of a porous media to transmit water through its voids, is one of the most significant key parameters of geomaterials for many natural phenomena including the management of water resources, drinking water supply, safety of waste repositories, basin-scale hydrogeologic circulation, stability analyses, and many other problems on subsurface hydrology and geotechnical engineering (Terzaghi and Peck 1964; Moore et al. 1982; Wintsch et al. 1995; Person et al. 1996; Boadu 2000; Chapuis 2012). There have been attempts to estimate hydraulic conductivity based on grain size distribution (Mualem 1976; Freeze and Cherry 1979; Uma et al. 1989; Salarashayeri and Siosemarde 2012). Empirical (Hazen 1911; Krumbein and Monk 1942; Alyamani and Sen 1993) and predictive methods (Kozeny 1927; Carman 1937; Boadu 2000; Goktepe and Sezer 2010) of estimating the hydraulic conductivity using quantitative relations have been developed in the literature. A commonly accepted equation was proposed by Hazen (1911) and given k = cd ²₁₀ for predicting the hydraulic conductivity of saturated sands. Where k is hydraulic conductivity, c is constant, and d ₁₀ is effective diameter at which 10 % of the grains are finer. Krumbein and Monk (1942) gave an expression for the hydraulic conductivity of unconsolidated sands by an empirical equation of the form k = (760d ² _w )exp(−1.3σ _ψ ), where d _w is geometric mean diameter by weight in millimetres, σ _ψ is standard deviation of the ψ distribution function. Masch and Denny (1966) proposed the use of d ₅₀ median grain size as the representative size to correlate hydraulic conductivity with grain size. Kozeny (1927) and Carman (1937), which is widely accepted derivation for hydraulic conductivity, developed a semi-empirical formula for predicting the permeability of porous media. Koltermann and Gorelick (1995) stated that the use of geometric mean overpredicts hydraulic conductivity by several orders of magnitude for soils with significant fines content, whilst the harmonic mean grain size under predicts k by several orders of magnitude for soils with less fines content. Shepherd (1989) performed a series of statistical power regression analyses on 19 sets of published data on hydraulic conductivity of unconsolidated sediments versus grain size. Alyamani and Sen (1993) proposed an equation based on analysis of 32 samples incorporating the initial slope and the intercept of the grain-size distribution curve. Sperry and Peirce (1995) developed a model for delineating the significance of particle size/shape, and porosity in explaining the variability of hydraulic conductivity of a granular porous medium. Ishaku et al. (2011) have employed several empirical formulae to specify the hydraulic conductivity of aquifer materials in the field. Although many different techniques have been proposed to determine hydraulic conductivity value, including field methods, applications of these empirical formulae to the same porous medium material can yield different values of hydraulic conductivity because of the difficulty of including all possible variables in porous media (Vukovic and Soro 1992).

It has been long understood that grain shape characteristics have a significant effect on certain engineering properties of soils (Terzaghi 1925; Gilboy 1928; Lees 1964; Olson and Mesri 1970; Abbireddy et al. 2009; Clayton et al. 2009). Terzaghi is one of the first engineers to perform a research to understand the influences of shape characteristics by employing flat-grained constituents (Terzaghi 1925). The observations, conducted by Gilboy (1928), that any system of analysis neglecting the effect of grain shape would be incomplete. Numerous researches have been conducted due to the significance of grains’ shape and its role in the behaviour of soils for both practicing engineers and researchers. Holubec and D’Appolonia (1973) indicated that the results of dynamic penetration tests in sands depend on grains’ shape characteristics. Cornfort (1973), and Holtz and Kovacks (1981) pointed out how grain shape affects the internal friction angle (φ). Cedergen (1989) stated that grain shape affects the permeability. Grain shape also plays an important role in liquefaction potential (Kramer 1996). Wadell (1932), Krumbein (1941), Powers (1953), Holubec and D’Appolonia (1973), Youd (1973), and Cho et al. (2006) have introduced detailed explanations of grain shape. Two independent properties are basically used to describe the shape of a soil grain: (1) Roundness, a measure of the extent to which the edges and corners of a grain has been rounded (2) Sphericity (form), a measure of the extent to which a grain approaches a sphere in shape. Wadell (1932) proposed a simplified sphericity (S) parameter (D_max-insc/D_min-circ), where D_max-min is the diameter of a maximum inscribed circle and D_min-circ is the diameter of a minimum sphere circumscribing a gravel particle. Wadell (1932) defined roundness (R) as D_i-ave/D_max-insc, where D_i-ave is the average diameter of the inscribed circle for each corner of the particle. Figures 1, 2 and 3 describe R, S and a chart for comparison between them to identify grain shape (Krumbein 1941; Powers 1953).

Fig. 1

Graphical representation of roundness, R (redrawn from Muszynski and Stanley, 2012)

Fig. 2

Graphical representation of sphericity, S (redrawn from Muszynski and Stanley, 2012)

Fig. 3

Comparison chart (Santamarina and Cho 2004)

Although many field and laboratory determinations of hydraulic conductivity have been performed by engineers, geologist, hydrogeologist, and soil scientists, the fundamental relationships between the gradation and shape properties of grains and flow through them remain poorly understood and inadequately quantified. Actually, these approaches cannot yield consistent results with respect to actual hydraulic conductivity values. Therefore, this study aims to evaluate a new conceptual approach for quantifying the inherent coupling between gradation/shape of sand grains changes and hydraulic conductivity by exploiting constant head permeability tests on sixteen different grain-size fractions (4.75–2, 2–1.18, 1.18–0.6, 0.6–0.425, 0.425–0.3, 0.3–0.075, 4.75–0.075, 2–0.075, 1.18–0.075, 0.6–0.075, 0.425–0.075, 4.75–0.6, 2–0.6, 4.75–0.425, 2–0.425, 1.18–0.425 mm) of sands having two distinct shapes (rounded and angular). Furthermore, comparisons between results obtained in the present study and hydraulic conductivity estimated with other formulas available in the literature were made.

Experimental study

The materials used in the tests described in this study were Narli Sand (NS) and Crushed Stone Sand (CSS) having the distinct shapes and sizes falling between 4.75 and 2 mm, 2 and 1.18 mm, 1.18 and 0.6 mm, 0.6 and 0.425 mm, 0.425 and 0.3 mm, 0.3 and 0.075 mm, 4.75 and 0.075 mm, 2 and 0.075 mm, 1.18 and 0.075 mm, 0.6 and 0.075 mm, 0.425 and 0.075 mm, 4.75 and 0.6 mm, 2 and 0.6 mm, 4.75 and 0.425 mm, 2 and 0.425 mm, 1.18 and 0.425 mm. Narli Sand (NS) was quarried in and around Narli, Kahramanmaras in southern-central of Turkey. A commercially available Crushed Stone Sand (CSS) was supplied from the same region of Turkey, which is widely consumed in earthworks in the region. The specific gravity of the grains were found to be 2.65 for Narli Sand, and 2.68 for Crushed Stone Sand. Scanning Electron Micrograph (SEM) pictures show the physical differences/similarities among the sands used during this investigation (Fig. 4). As can be seen from the Fig. 4, Narli Sand grains have rounded, whereas the Crushed Stone Sand grains have angular shape. Figure 5 indicates the grain size distribution of the sands used during the experimental study. Roundness (R) and sphericity (S) estimations based on the study by Muszynski and Stanley (2012) were found to be 0.43, 0.67, and 0.16, 0.55 for the NS and CSS grains, respectively. The sands were tested in a constant head permeability testing apparatus at a relative density (R_d) of about 40 % and constant room temperature (20 ± 2 °C). The specimens, which were placed in a perspex cylindrical cell of about 50 cm² cross-sectional area (A), rest on a wire mesh at bottom of the cell. The volume of the water (q) flowing during a certain time (t) is measured, when a steady vertical water flow, under a constant head, is maintained through the soil specimen. Then, k values of the specimens tested were calculated using Darcy’s law (k = ql/Ah). Tables 1 and 2 present some physical characteristics of the NS and CSS samples, respectively. As can be seen from these tables the hydraulic conductivity is affected by grading characteristics d₁₀, d₂₀, d₃₀, d₅₀, d₆₀, c_u, c_c, n, and I_o.

Fig. 4

SEM pictures of the (top) CSS and (bottom) NS used during the experimental study

Fig. 5

Grain size distributions for the sands used during the experimental study

Table 1 Some physical characteristics of the NS samples

Table 2 Some physical characteristics of the CSS samples

Results and discussion

Table 3 gives a summary of the specimens used in the tests reported here. The initial relative densities of all specimens were around 40 %. The specimens were loose to medium dense. Sixteen different sizes of artificially graded NS and CSS sands, which have exactly the same gradation characteristics (d ₁₀, d ₂₀, d ₃₀, d ₅₀, d ₆₀, c _u , c _c , I _o ) (Fig. 5) within the specified ranges, have been classified as ‘poorly graded’ (SP) based on the Unified Soil Classification System (USCS9. Based on the roundness criteria and values proposed by Powers (1953), and Youd (1973), the specimens used during the experimental investigation were found to be very angular and rounded for CSS and NS grains, respectively.

Table 3 Summary of specimen data

Table 4 shows the empirical equations and their limitations for hydraulic conductivity estimates which were used to obtain the results given in Table 3. Equations developed by Hazen (1892), Kozeny-Carman (1956), Terzaghi (Odong 2007), Chapuis (2004), Slichter (1898), USBR (Vukovic and Soro 1992), NAVFAC (1974), Alyamani and Sen (1993), and Breyer (Kresic 1998) were employed in this study. Hazen (1892) proposed his formula in order to estimate the hydraulic conductivity of uniformly graded loose sand with effective grain size (d ₁₀) between 0.10 and 3.0 mm, and c _u less than 5. As can be seen from the Table 3 that hydraulic conductivity values ranged from 5.95 to 0.01 cm/s for the NS samples falling specified gradations, whilst those ranged from 8.39 to 0.02 cm/s for the CSS samples falling the same gradations. Although, presence of porosity (n) in the formula seems an advantage of the formula, this approach does not give an accurate estimates for the sands due to the limits of c _u indicated in Table 4. The authors consider that influence of the parameter c _u was neglected in his study, and thereby the grain size distribution results could yield the same c _u for various sands. Kozeny–Carman (K–C) formula, which is not applicable for neither clayey soils nor soils with effective size more than 3 mm, is one of the commonly employed approaches developed for hydraulic conductivity estimates (Carrier 2003). Actually, the Kozeny (1927) and Carman (1937) equations have been modified by certain researchers (Collins 1961; Bear 1972; de Marsily 1986), whom included the influence of both particle diameter and porosity on hydraulic conductivity. Koltermann and Gorelick (1995) compared five different approaches and found that the original Kozeny–Carman equation (Carman 1937; Bear 1972) lies approximately in the center of the possible relations. Koltermann and Gorelick (1995) used the geometric and harmonic means to calculate representative particle diameters for the high and low fraction of the coarse component, respectively. However, this approach produces a discontinuity when the fraction of the coarse component is at the intermediate level. Therefore, the authors employed the original Kozeny–Carman equation, then the Table 3 released that hydraulic conductivity values ranged from 4.93 to 0.02 cm/s for the NS samples, while those ranged from 12.36 to 0.03 cm/s for the CSS samples falling the same gradations. Estimated hydraulic conductivity values (k) by employing Terzaghi’s approach varied from 2.38 to 0.01 cm/s for the NS samples, whilst the k values varied from 4.85 to 0.01 cm/s for CSS samples. Cheng and Chen (2007) pointed out that Terzaghi’s formula is most applicable for large-grain sand. However, comparing the experimental results and the k values obtained via Terzaghi’s approach revealed that Terzaghi’s equation, which has no limitations reported (Table 4), gives more accurate results than the other equations employed for both NS and CSS samples between 1.18 and 0.075 mm, and 0.6 and 0.075 mm. Surprisingly, it gives much less accurate results for larger grains of both NS and CSS samples, including the size of 4.75–2, 2–1.18, and 4.75–0.425 mm. Therefore, the authors interpreted that grain size would not be the only parameter to make an accurate hydraulic conductivity estimate. Estimated k values via Chapuis formula gives the best correlation with measured k values for the NS samples between 0.425 and 0.075 mm. Generally speaking, estimated k values using Chapuis’s approach ranged from 6.31 to 0.01 cm/s for the NS samples, whilst those ranged from 6.33 to 0.01 cm/s for the CSS samples falling the same gradations. In the light of the Goktepe and Sezer (2010), which indicated that Chapuis method best estimates the hydraulic conductivity of fine sands, the predictions were found to be acceptable for the NS samples but not for the CSS samples. The authors considered that such difference could be because of shape properties of the sand grains. Although Goktepe and Sezer (2010) indicated that the Chapuis and Slitcher approaches are in harmony with the results, the present study shows remarkable differences between these two approaches. Considering the differences in relative density values employed in these studies, the authors’ interpretation is that such differences in the approaches could be the reason of high successes of the empirical equations. For example, the present study shows that Slitcher formula is the best fitted to the hydraulic conductivity of NS samples between 4.75 and 2 mm, 2 and 1.18 mm, 4.75 and 0.075 mm, 2 and 0.075 mm, 4.75 and 0.6 mm, 2 and 0.6 mm, 4.75 and 0.425 mm, 2 and 0.425 mm, 1.18 and 0.425 mm, and the hydraulic conductivity of CSS samples between 4.75 and 2 mm, 4.75 and 0.075 mm, 2 and 0.075 mm, 4.75 and 0.6 mm, 4.75 and 0.425 mm, 2 and 0.425 mm. However, Chapuis approach does not give similar results. The Naval Facilities Engineering Command (NAVFAC) suggested a chart to estimate the hydraulic conductivity of clean sand and gravel based on the e and d ₁₀. Predicted k values using NAVFAC varied from 8.48 to 0.01 cm/s for the NS samples, and 13.24 to 0.04 cm/s for the CSS samples. The approach proposed by the United States Bureau of Reclamation (USBR 1990) estimates k values using the effective grain size (d ₂₀), and it does not depend on the porosity (Table 4). Cheng and Chen (2007) stated that this approach is most suitable for medium-grain sand with c _u less than 5. Estimated k values using the USBR formula were found to be same for NS samples and CSS samples, which ranged from 4.46 to 0.01 cm/s, as they have the same gradations. It was observed that the USBR approach gave its best results for relatively large grain samples including those between 2 and 1.18 mm, 1.18 and 0.6 mm, 0.6 and 0.425 mm, and 1.18 and 0.425 mm. Alyamani and Sen (A–S), which is one of the widely known approaches to estimate the hydraulic conductivity, employs the grain size properties d ₁₀, d ₅₀ and I _o . Alyamani and Sen (1993) proposed their equation based on different samples that incorporates the initial slope and the intercept of the grading curve. Estimated k values using the Alyamani and Sen approach ranged from 6.16 to 0.01 cm/sec for both type of sands. As can be seen from Table 3 that the A–S approach results in same estimates for both NS and CSS samples, as they have same grading curves. Similarly, Breyer method gave the same k values for both NS and CSS samples due to the same d ₁₀ value employed in this equation. The predicted k values ranged from 7.21 to 0.01 cm/s. Plots presented in Figs. 6 and 7 indicate comparisons of measured hydraulic conductivity (k) with predictions from various models for NS samples, and CSS samples, respectively.

Table 4 Empirical equations and their limitations for permeability estimates

Fig. 6

Comparison of measured hydraulic conductivity (k) with predictions from various models for NS samples (straight line represents line of perfect equality)

Fig. 7

Comparison of measured hydraulic conductivity (k) with predictions from various models for CSS samples (straight line represents line of perfect equality)

The differences between measured and predicted hydraulic conductivity values using various equations were because of either inaccuracy in measured soil parameters or deficiency in the predictive equations. Therefore, Table 5 and 6 were complied in order to present a comparative study for the NS and CSS samples using all the formulas employed in this study, respectively. The Tables 5 and 6 show the results of calculations performed with the objective of determining hydraulic conductivity according to the nine different approaches (Hazen, Slitcher, K–C, Terzaghi, USBR, Chapuis, A–S, Breyer, NAVFAC), expressed as a relative ratio of the difference between estimated and calculated values to the estimated hydraulic value of the NS and CSS samples at sixteen different gradations (4.75–2, 2–1.18, 1.18–0.6, 0.6–0.425, 0.425–0.3, 0.3–0.075, 4.75–0.075, 2–0.075, 1.18–0.075, 0.6–0.075, 0.425–0.075, 4.75–0.6, 2–0.6, 4.75–0.425, 2–0.425, 1.18–0.425 mm). The nine approaches used for comparison were listed from the best fitting on left to the worst fitting on right. For example, the best estimation of hydraulic conductivity for the NS samples between 4.75 mm and 2 mm was found to be based on Slitcher equation, followed by Terzaghi, USBR, Kozeny–Carman, Hazen, Alyamani–Sen, Chapius, Breyer, and NAVFAC equations, respectively. The authors have observed that, as an overall view, Slitcher and Terzaghi’s approaches give the best correlation with measured k values for both NS and CSS samples, whilst Kozeny–Carman and NAVFAC approaches give the worst correlation with measured k values for both NS and CSS samples for any gradation.

Table 5 Comparisons for the NS samples

Table 6 Comparisons for the CSS samples

Nevertheless, despite the good predictions in certain grading of samples, the authors interpreted that reliability of these approaches is relatively low as that any system of analysis neglecting the effect of grain shape would be incomplete. Effect of gradation as well as grain shape on hydraulic conductivity values have been presented in Figs. 8 and 9. Effects of five different gradation including 4.75–2, 2–1.18, 1.18–0.6, 0.6–0.425, and 0.425–0.3 mm on hydraulic conductivity of NS and CSS samples were illustrated in Fig. 8. The highest value of hydraulic conductivity for the NS was found to be for the samples between 4.75 and 2 mm, and then followed by the samples between 2–1.18, 1.18–0.6, 0.6–0.425, and 0.425–0.3 mm, respectively. Effects of grain shape on hydraulic conductivity values was clearly seen in Fig. 9, which proves that samples with two different shapes could have a unique hydraulic conductivity value, likely due to the differences in shape characteristics (R, S) leading to the different void ratios (e).

Fig. 8

Effects of gradation on hydraulic conductivity values of NS and CSS samples

Fig. 9

Effects of grain shape on hydraulic conductivity values of tested samples

Conclusions

The objective of this research was to study the influences of gradation and grain shape on hydraulic conductivity of soils, which is of importance in relation to certain geotechnical problems including stability analyses, settlement and seepage computations. The samples used in the present study are composed of poorly graded Narli Sand (NS) and Crushed Stone Sand (CSS), which were found to be rounded (R = 0.72, S = 0.67) and very angular (R = 0.15, S = 0.55), respectively. Sixteen ranges of grain sizes (4.75–2, 2–1.18, 1.18–0.6, 0.6–0.425, 0.425–0.3, 0.3–0.075, 4.75–0.075, 2–0.075, 1.18–0.075, 0.6–0.075, 0.425–0.075, 4.75–0.6, 2–0.6, 4.75–0.425, 2–0.425, and 1.18–0.425 mm) of both NS and CSS samples were tested in a constant head permeability testing apparatus at a relative density (D_r) of about 40 %. Moreover, various predictive methods of estimating the hydraulic conductivity values (Hazen, Kozeny–Carman, Terzaghi, Chapuis, Slitcher, USBR, NAVFAC, Alyamani and Sen, and Breyer) have been employed to compare the measured and estimated hydraulic conductivity results. In general, the Slitcher and Terzaghi’s approaches give the best correlation with measured k values for both NS and CSS samples, whilst Kozeny–Carman and NAVFAC approaches give the worst correlation with measured k values for both NS and CSS samples for any gradation. The test results and comparative study reported here in this paper indicate following facets of behavior:

• 1. The hydraulic conductivity values of the NS samples with rounded grains were lower than those of the CSS samples with very angular grains, which is likely to be the result of shape characteristics leading different void ratios.

• 2. The hydraulic conductivity can be significantly influenced by grading characteristics including d₁₀, d₂₀, d₃₀, d₅₀, d₆₀, c_u, c_c, n, and I_o.

• 3. Gradation of the grains have a significant effect on hydraulic conductivity of both NS and CSS samples.

• 4. The comparative study on the perceptions of estimated and predicted results with other approaches available in the literature indicated that the best prediction of hydraulic conductivity changes based on the gradation and shape properties of the sands tested.