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

Hydraulic conductivities of sands with different gradation and grain shape were estimated experimentally at a relative density (Dr) 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 (d10, d20, d30, d50, d60, cu, cc, n, Io). 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.

using quantitative relations have been developed in the literature. A commonly accepted equation was proposed by Hazen (1911) and given k = cd 2 10 for predicting the hydraulic conductivity of saturated sands. Where k is hydraulic conductivity, c is constant, and d 10 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 2 )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 50 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 andD' 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).
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 grainsize 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- Muszynski and Stanley, 2012) 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 2 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    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 10 , d 20 , d 30 , d 50 , d 60 , 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 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 10 ) 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,    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  (USBR 1990) estimates k values using the effective grain size

Researcher/organization Equation Limitations
Hazen k = 6 × 10 −4 × g v × [1 + 10(n − 0.26)] × (d 10 ) 2 C u < 5 0.1 < d 10 < 3.0 (1−n) 2 × (d 10 ) 2 0.5 < d 10 < 4.0 Terzaghi Alyamani and Sen k = 1300 × [I 0 + 0.025(d 50 − d 10 )] 2 -Breyer k = 6 × 10 −4 × g v × log 500 Cu × (d 10 ) 2 0.06 < d 10 < 0.6 1 < C u < 20 (d 20 ), 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 10 , d 50 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 10 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. 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  (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  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  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).

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.  Fig. 9 Effects of grain shape on hydraulic conductivity values of tested samples 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 10 , d 20 , d 30 , d 50 , d 60 , 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.