This experimental study is an attempt to apply the constriction size concept for determination of hydraulic conductivity of filter material. Five shapes of granular filter material of different gradations were used for the experimental investigations. The controlling constriction size for selected shapes and gradations were worked out and related to the experimental values of hydraulic conductivity. The empirical model for determination of hydraulic conductivity based on controlling constriction size for each shape of granular filter material was developed. Finally, a combined empirical model for the determination of hydraulic conductivity of granular filters was developed based on shape parameters and the controlling constriction size. The model was validated with the experimental data of a past research work; the predicted values of hydraulic conductivity were found in close approximation to the values obtained from the literature. Results of the present work will help in improving the existing design criteria of protective filters, which form an integral part of the safety of hydraulic structures.

  • Inclusion of constriction size and shape parameters in the hydraulic conductivity determination of granular filters.

  • Improvement of the effectiveness of the filter design based on permeability.

  • Enhancement of water resources management in the essential hydraulic structure design where filters form a critical component.

Granular filters are important structural components of hydraulic structures provided with an aim to prevent the migration of fine particles of base soil and at the same time allow the free movement of seepage flow in order to avoid the pore pressures. Hydraulic conductivity denotes the ease with which flow takes place through any porous medium (Terzaghi 1922; Jabro 1992). The importance of hydraulic conductivity for the design of filters has been discussed by many researchers (Bertram 1940; US Army Corps 1953; Sherard et al. 1984; Kenney et al. 1985; USBR 1987; Indraratna et al. 1996; Alsakran & Zhu 2020). Boadu (2000) concluded that the measurement of hydraulic conductivity by laboratory and field methods was costly and time consuming, and insufficient field data of aquifers and hydraulic boundaries add to the limitation (Freeze & Cherry 1979; Todd & Mays 2005). These limitations have led researchers to shift to analytical approaches for the study of groundwater and the porous media (Salmasi & Azamathulla 2013; Cheng & Hsu 2021). Hazen (1893), Carman (1939), Sherard et al. (1984), Vukovic & Soro (1992), Indraratna et al. (1996) and Ishaku et al. (2011) conducted investigations to suggest empirical equations for the hydraulic conductivity determination based on particle sizes of granular soils. Kenney et al. (1984), Odong (2007) and Trani & Indraratna (2010) assessed void network of the granular medium as an important parameter for hydraulic conductivity compared to the particle size. Indraratna et al. (2007) highlighted the advantages of applying the method of constriction size for granular soils. The relationship between constriction size and hydraulic conductivity of granular soils was investigated by Indraratna et al. (2012) and Seblany (2018). Also, studies, although limited have been conducted on the influence of shape and gradation parameters on filter permeability (Rather et al. 2020). Several studies have been carried out for determination of hydraulic conductivity of granular soils based on particle size, shape parameters as well as constriction size of void network. It was concluded from the literature that there was less focus on the combined effect of constriction size and shape parameters on the hydraulic conductivity of granular filters.

The objective of the present study was to determine hydraulic conductivity of granular material using constriction size of the void network and the shape parameters of the filter material. Laboratory tests were carried out to investigate their relationship; for this purpose five shapes of granular filter material were used; further nine gradations for each shape were studied to deduce the relationship between these parameters. Experimental data was further analysed to arrive at an empirical model for prediction of hydraulic conductivity that considers constriction size of void network and shape parameters of filter material.

In the present study, hydraulic conductivity of filter material was determined using two important parameters i.e. constriction size and shape parameters. Constriction size (d*) is the smallest opening between the void spaces of filter material. This size is responsible for controlling the migration of fine particles of base material and determines the value of hydraulic conductivity. Incorporation of constriction size in the hydraulic conductivity computations refines the accuracy of results. Many theoretical concepts have been suggested to work out the constriction size from the particle gradation (Lone et al. 2005, Indraratna & Raut 2006, Indraratna et al. 2007; Raut & Indraratna 2008; Indraratna et al. 2015). In the present work, the constriction size of void network of filter material was determined by the procedure given by Lone et al. (2005). A unit assembly of three non-uniform spheres was considered and these spheres were assumed to be in contact with each other in a manner such that their centres were in a single plane forming a window whose sides are of concave arcs with radii equal to the radii of the assembly spheres. The arrangement of unit assembly is shown in Figure 1.

Figure 1

Packing pattern of non-uniform spheres of filter assembly (Lone et al. 2005).

Figure 1

Packing pattern of non-uniform spheres of filter assembly (Lone et al. 2005).

Close modal
The packing arrangement (Lone et al. 2005) consists of three spheres of radius r1, r2, r3 and diameter D1, D2, D3, and a fourth sphere of radius r4 and diameter D4. The sphere of diameter D4 is the function of the three spheres of diameter D1, D2, D3 and other pore sizes like D5, D6, and D7 and so on. The size of fourth sphere can be calculated by Equation (1) given by Lone et al. (2005):
(1)
where, , , and .

The m, n, and curves given by Lone et al. (2005) can be used to determine the values of . The above discussed model has been used to determine the constriction size of filter mass by grouping the entire filter material into different skeleton particle sizes.

The material used for experimentation and test procedure of the experimental runs for the present study are discussed in this section.

Material

The granular filter material required for the laboratory investigations was obtained from a local river bed and grouped into different standard sizes ranging from 4.75 to 50 mm on the basis of particle diameter. Further, the material was sorted into five shapes, designated as S1, S2, S3, S4 and S5 based on their shape parameters. The shape parameters considered in the present study are sphericity, shape factor, elongation ratio, and flatness ratio; designated as Sph, SF, ER and FR respectively. Nine filter material gradations (M, M2, M3, M4, M5, M6, M7, M8 and M9) of each shape were prepared for laboratory investigations. The shape parameters of the filter material used in the present study are given in Table 1.

Table 1

Filter shape parameters

ShapesSphSFFRER
S1 0.892 0.798 1.269 0.930 
S2 0.767 0.744 1.374 0.721 
S3 0.747 0.528 1.960 0.774 
S4 0.838 0.731 1.406 0.864 
S5 0.504 0.378 2.65 0.397 
ShapesSphSFFRER
S1 0.892 0.798 1.269 0.930 
S2 0.767 0.744 1.374 0.721 
S3 0.747 0.528 1.960 0.774 
S4 0.838 0.731 1.406 0.864 
S5 0.504 0.378 2.65 0.397 

Hydraulic conductivity determination

The constant head permeability apparatus was used to determine the hydraulic conductivity of filter material (ASTM D2434-68 2006) as shown in Figure 2. The test apparatus consists of a cylinder having diameter of 250 mm, length of 600 mm with an 80 mm diameter of hopper base. The movement of filter material in the steel cylinder due to the inflowing water was prevented by providing the wire mesh of different openings ranging from 2.5 to 8 mm. At the top plate of the cylinder, a pressure gauge and two air vents were provided. A stopcock was also provided to regulate the water supply at the inlet. The outflow from the apparatus was allowed into the measuring tank with the help of a flexible pipe of reasonable height for measuring the discharge of water required for the test. Metric staffs were provided for the measurement of head.

Figure 2

Experimental set-up for determination of hydraulic conductivity.

Figure 2

Experimental set-up for determination of hydraulic conductivity.

Close modal
The apparatus was filled with the filter material up to 100–150 mm below the top of cylinder; a wire mesh of 10 mm opening was placed at the top of the cylinder. The pea gravel was placed over the filter material to prevent the direct impact of inflowing water on the filter material. The wire mesh of 10 mm opening was also placed between the pea gravel and the filter material to differentiate the two types of material. Darcy's law was used to determine the permeability of the filter material at different temperatures and then standardized to viscosity of water at 20 degree Celsius, using Equation (2).
(2)

In Equation (2), the coefficient of permeability and the viscosity of water at 20 °C and T°C are denoted by k20, kt, μ20 and μt respectively.

The prepared samples of filter material were tested and the laboratory results are given in Table 2. The experimental results clearly reveal the variation of hydraulic conductivity with the shape of filter material even for similar gradations. The relationship of hydraulic conductivity of filter material and constriction size for the five shapes are given in Figure 3(a)–3(e). The corresponding empirical models are given in Table 3.

Figure 3

(a) Variation of hydraulic conductivity with constriction sizes for Shape 1 of filter material. (b) Variation of hydraulic conductivity with constriction sizes for Shape 2 of filter material. (c) Variation of hydraulic conductivity with constriction sizes for Shape 3 of filter material. (d) Variation of hydraulic conductivity with constriction sizes for Shape 4 of filter material. (e) Variation of hydraulic conductivity with constriction sizes for Shape 5 of filter material.

Figure 3

(a) Variation of hydraulic conductivity with constriction sizes for Shape 1 of filter material. (b) Variation of hydraulic conductivity with constriction sizes for Shape 2 of filter material. (c) Variation of hydraulic conductivity with constriction sizes for Shape 3 of filter material. (d) Variation of hydraulic conductivity with constriction sizes for Shape 4 of filter material. (e) Variation of hydraulic conductivity with constriction sizes for Shape 5 of filter material.

Close modal
Table 2

Hydraulic conductivity of the filter material

Shape
S1S2S3S4S5
Gradationd*
(cm)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
M1 0.310 95,590.8 63,000.0 46,440.0 37,350.0 108,565.2 
M2 0.250 48,391.2 33,174.0 30,250.8 22,176.0 57,823.2 
M3 0.216 27,421.2 24,433.2 20,581.2 17,434.8 34,585.2 
M4 0.1575 15,156.0 14,137.2 13,921.2 13,500.0 17,604.0 
M5 0.131 13,140.0 12,708.0 12,276.0 11,952.0 14,086.8 
M6 0.130 13,957.2 11,725.2 10,062.0 9,205.2 16,225.2 
M7 0.180 20,077.2 14,832.0 13,147.2 11,674.8 22,978.8 
M8 0.217 25,657.2 19,242.0 18,039.6 16,081.2 30,600.0 
M9 0.171 15,685.2 14,209.2 10,573.2 9,838.8 18,756.0 
Shape
S1S2S3S4S5
Gradationd*
(cm)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
k
(cm/hr)
M1 0.310 95,590.8 63,000.0 46,440.0 37,350.0 108,565.2 
M2 0.250 48,391.2 33,174.0 30,250.8 22,176.0 57,823.2 
M3 0.216 27,421.2 24,433.2 20,581.2 17,434.8 34,585.2 
M4 0.1575 15,156.0 14,137.2 13,921.2 13,500.0 17,604.0 
M5 0.131 13,140.0 12,708.0 12,276.0 11,952.0 14,086.8 
M6 0.130 13,957.2 11,725.2 10,062.0 9,205.2 16,225.2 
M7 0.180 20,077.2 14,832.0 13,147.2 11,674.8 22,978.8 
M8 0.217 25,657.2 19,242.0 18,039.6 16,081.2 30,600.0 
M9 0.171 15,685.2 14,209.2 10,573.2 9,838.8 18,756.0 
Table 3

Empirical models for determination of hydraulic conductivity for different shapes of filter particles

ShapeProposed empirical modelCo-efficient of determination (R2)Remarks
S1  90% Significant 
S2  88% Significant 
S3  86% Significant 
S4  82% Significant 
S5  92% Significant 
ShapeProposed empirical modelCo-efficient of determination (R2)Remarks
S1  90% Significant 
S2  88% Significant 
S3  86% Significant 
S4  82% Significant 
S5  92% Significant 

The experimental results obtained in this study indicate that the hydraulic conductivity of filter material having similar constriction sizes differs with the filter particle shape. Hence, it is necessary to incorporate the shape parameters for prediction of hydraulic conductivity of filter material based on constriction size. The relationship of hydraulic conductivity of filter material with constriction size of void network and individual shape parameters are given in Figure 4(a)–4(e).

Figure 4

(a) Variation of hydraulic conductivity (k) with d*/SPH, (b) Variation of hydraulic conductivity (k) with d*/SF, (c) Variation of hydraulic conductivity (k) with d*/FR, (d) Variation of hydraulic conductivity (k) with d*/ER.

Figure 4

(a) Variation of hydraulic conductivity (k) with d*/SPH, (b) Variation of hydraulic conductivity (k) with d*/SF, (c) Variation of hydraulic conductivity (k) with d*/FR, (d) Variation of hydraulic conductivity (k) with d*/ER.

Close modal

The proposed empirical models for prediction of hydraulic conductivity based on controlling constriction size and individual shape parameters are given in Table 4.

Table 4

Empirical models for determination of hydraulic conductivity for different shape parameters of filter particles

Shape parameterProposed empirical modelCo-efficient of determinationRemarks
Sph  68% Significant 
SF  53% Fair 
FR  23% Non-significant 
ER  57% Fair 
Shape parameterProposed empirical modelCo-efficient of determinationRemarks
Sph  68% Significant 
SF  53% Fair 
FR  23% Non-significant 
ER  57% Fair 

The experimental results reveal that constriction size of void network and shape parameters of filter particles have significant effect on hydraulic conductivity of filter material. The relationship of hydraulic conductivity of filter material with constriction size of void network and shape parameters (sphericity, flatness ratio, shape factor and elongation ratio) is given in Figure 5.

Figure 5

Variation of hydraulic conductivity (k) with (d*/ER)*Sph*SF*FR.

Figure 5

Variation of hydraulic conductivity (k) with (d*/ER)*Sph*SF*FR.

Close modal
The model developed based on power regression is given as Equation (3). This model has R2 value of 79% which makes the model significant and hence, acceptable:
(3)
where d* is constriction size (cm) and k is hydraulic conductivity (cm/hr).

The model proposed for determining hydraulic conductivity of filter material was validated by selecting permeability data, shape parameters and constriction size of filter material used in the research work carried out by Lone et al. (2005). The computed hydraulic conductivity values, shape parameters and constriction size of filter material used in the study by Lone et al. (2005) and the predicted values of hydraulic conductivity by the proposed model Equation (3) are given in Table 5. The values of hydraulic conductivity of filter material computed by Lone et al. (2005) were in close approximation to the values of hydraulic conductivity of filter material predicted by Equation (3). The observed versus predicted hydraulic conductivity of filter material is presented in Figure 6.

Table 5

Observed and predicted values of hydraulic conductivity

Sphericity
(Sph)
Shape factor
(SF)
Flatness ratio
(FR)
Elongation ration
(ER)
Controlling constriction size
(d*)
cm
Computed hydraulic conductivity
(k)
cm/sec
Predicted hydraulic conductivity (k) cm/sec
0.602 0.592 1.845 0.557 0.108 3.7 2.52 
0.647 0.621 1.715 0.576 0.131 4.58 3.65 
0.789 0.747 1.371 0.735 0.1575 5.66 4.371 
0.802 0.714 1.44 0.765 0.216 6.6 7.42 
0.812 0.733 1.394 0.779 0.25 11.12 9.437 
Sphericity
(Sph)
Shape factor
(SF)
Flatness ratio
(FR)
Elongation ration
(ER)
Controlling constriction size
(d*)
cm
Computed hydraulic conductivity
(k)
cm/sec
Predicted hydraulic conductivity (k) cm/sec
0.602 0.592 1.845 0.557 0.108 3.7 2.52 
0.647 0.621 1.715 0.576 0.131 4.58 3.65 
0.789 0.747 1.371 0.735 0.1575 5.66 4.371 
0.802 0.714 1.44 0.765 0.216 6.6 7.42 
0.812 0.733 1.394 0.779 0.25 11.12 9.437 
Figure 6

Observed versus predicted hydraulic conductivity.

Figure 6

Observed versus predicted hydraulic conductivity.

Close modal

The present work was carried out with an objective to determine the relationship of hydraulic conductivity with constriction size and shape parameters of filter material. The main findings of the study are summarized as follows:

  • Hydraulic conductivity of filter material was found to be the function of controlling constriction size and shape parameters.

  • Empirical relations were developed for hydraulic conductivity determination of filters based on constriction size of different filter particle shapes.

  • An empirical model based on power regression of significant parameters like constriction size and the shape parameters of filter material was developed for the computation of hydraulic conductivity.

The prediction of hydraulic conductivity is a pre-requisite for the design of structures that encounter seepage flow such as hydraulic structures constructed on permeable soils. To ensure the safety of such structures, protective filters form an integral component of the set-up. The design of protective filters depends upon the hydraulic conductivity of filter material and the base material. The model proposed in the present study is more suitable for the prediction of hydraulic conductivity as it encompasses shape as well as the constriction size of the filter assembly. It is very clear from the results of this study that the shape of the filter particles and the constriction size have considerable effects on the hydraulic conductivity; therefore, their incorporation in the model for hydraulic conductivity determination will improve the accuracy of estimation.

The aim of the present work was to incorporate effect of shape into the empirical models for hydraulic conductivity determination. Because of the time limit of the study conducted and the exhaustive procedure of sorting the shapes, only five shapes of filter particles could be studied. This being one of the limitations of the study, it is highly recommended to extend this study for more shapes so as to improve the efficacy of the design criteria for filters. Also, there is scope of improving the outcome, if the number of shape parameters considered are increased and their determination is conducted in the laboratory by latest material testing techniques.

This paper presents the research work carried out by the author*, who received a PhD fellowship from the Ministry of Human Resources Development (MHRD), Government of India. The experimental work was carried out in the Fluid Mechanics Laboratory, NIT Srinagar, Jammu and Kashmir, India.

No potential conflict of interest was reported by the authors.

The data used for the development of empirical models in the present study are given in Table 2 in the Results and discussion section. The data used for the validation of proposed empirical model are given in Table 5 of the validation section.

All relevant data are included in the paper or its Supplementary Information.

Alsakran
M. A.
&
Zhu
J. G.
2020
The use of distance between soil layers in predicting the hydraulic conductivity of granular soils
.
IOP Conference Series Earth and Environmental Science
525
,
012086
.
ASTM D2434-68
2006
Standard Test Method for Permeability of Granular Soils (Constant Head)
.
ASTM International, American Society for Testing and Materials
,
West Conshohocken, PA, USA
.
Bertram
G. E.
1940
An Experimental Investigation of Protective Filters. Harvard Soil Mechanics Series No. 7, Publication No. 267
.
Graduate School of Engineering, Harvard University
,
Cambridge, MA, USA
.
Boadu
F. K.
2000
Hydraulic conductivity of soils from grain size distribution: new models
.
Journal of Geotechnical and Geo-Environmental Engineering
126
(
8
),
739
746
.
Carman
P. C.
1939
Permeability of saturated sands, soils and clays
.
Journal of Agricultural Sciences
29
(
2
),
262
273
.
Cheng
S. Y.
&
Hsu
K. C.
2021
Bayesian integration using resistivity and lithology for improving estimation of hydraulic conductivity
.
Water Resources Research
5
(
3
), e2020WR027346.
Freeze
R. A.
&
Cherry
J. A.
1979
Groundwater
.
Prentice Hall
,
Englewood Cliffs, NJ
.
Hazen
A.
1893
Some Physical Properties of Sands and Gravels with Special Reference to Their Use in Filtration. 24th Annual Report of the State Board of Health of Massachusetts
.
Wright and Potter Printing
,
Boston
.
Indraratna
B.
&
Raut
A. K.
2006
Enhanced criterion for base soil retention in embankment dam filters
.
Journal of Geotechnical and Geoenvironmental Engineering
132
(
12
),
1621
1627
.
Indraratna
B.
,
Vafai
F.
&
Dilema
E. L. G.
1996
An experimental study of filtration of a lateritic clay slurry by sand filters
.
Proceedings of Institution of Civil Engineers, Geotechnical Engineering
119
(
2
),
75
83
.
Indraratna
B.
,
Raut
A. K.
&
Khabbaz
H.
2007
Constriction-based retention criterion for granular filter design
.
Journal of Geotechnical and Geo-Environmental Engineering
133
(
3
),
266
276
.
Indraratna
B.
,
Nguyen
V. T.
&
Rujikiatkamjorn
C.
2012
Hydraulic conductivity of saturated granular soils determined using a constriction-based technique
.
Canadian Geotechnical Journal
49
(
5
),
607
613
.
Indraratna
B.
,
Israr
J.
&
Rujikiatkamjorn
C.
2015
Geometrical method for evaluating the internal instability of granular filters based on constriction size distribution
.
Journal of Geotechnical and Geoenvironmental Engineering
141
(
10
),
0405045-1
04015045-14
.
Ishaku
J. E.
,
Gadzama
E.
&
Kaigama
U.
2011
Evaluation of empirical formulae for the determination of hydraulic conductivity based on grain size analysis
.
Journal of Geology and Mining Research
3
(
4
),
105
113
.
Kenney
T. C.
,
Lau
D.
&
Ofoegbu
G. I.
1984
Permeability of compacted granular materials
.
Canadian Geotechnical Journal
21
(
4
),
726
729
.
Kenney
T. C.
,
Chahal
R.
,
Cheiu
E.
,
Ofegbu
G. I.
,
Omange
G. N.
&
Ume
C. A.
1985
Controlling constriction sizes of granular filters
.
Canadian Geotechnical Journal
22
,
32
43
.
Lone
M. A.
,
Hussain
B.
&
Asawa
G. L.
2005
Filter design criteria for graded cohesionless bases
.
Journal of Geotechnical and Geo-Environmental Engineering
131
(
2
),
251
259
.
Odong
J.
2007
Evaluation of empirical formulae for determination of hydraulic conductivity based on grain-size analysis
.
Journal of American Science
3
(
3
),
54
60
.
Rather
N. A.
,
Lone
M. A.
,
Dar
A. Q.
&
Mir
B. H.
2020
The influence of particle shape and gradation parameters on the permeability of filter media
.
International Journal of Hydrology Science and Technology
10
(
3
),
267
286
.
Raut
K.
&
Indraratna
B.
2008
Further advancement in filtration criteria through constriction based techniques
.
Journal of Geotechnical and Geo-Environmental Engineering
134
(
6
),
883
887
.
Seblany
F.
2018
Filter Criterion for Granular Soils Based on the Constriction Size Distribution
.
University de Lyon
,
Lyon, France
.
Sherard
J. L.
,
Dunnigan
L. P.
&
Talbot
J. R.
1984
Filters for silts and clays
.
Journal of Geotechnical Engineering
110
(
6
),
701
718
.
Terzaghi
K.
1922
Der Grundguch an Stauwerken and seine Verhutung [The failure of dams by piping and its prevention]
.
Die Wasserkraft
17
,
445
449
.
Todd
D. K.
&
Mays
L. W.
2005
Groundwater Hydrology
.
John Wiley and Sons
,
New York
.
US Army Corps of Engineers
1953
Filter Experiments and Design Criteria. Waterways Experiment Station, Vicksburg, MS., Technical Memorandum, No. 3, pp 360
.
USBR
1987
Design Standards-Embankment Dam, No. 13 Chapter 5, Protective Filters
.
U.S. Department of Commerce, National Technical Information Service
,
Springfield
.
Vukovic
M.
&
Soro
A.
1992
Determination of Hydraulic Conductivity of Porous Media From Grain-Size Distribution
.
Water Resources Publication, LLC Highlands Ranch
,
Colorado
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).