# Prediction Method for Hydraulic Conductivity considering the Effect of Sizes of Ellipsoid Soil Particles from the Microscopic Perspective

Authors: Caihua Shen; Jun Zhu; Wenbo Gu
Journal: Advances in Civil Engineering
 (2019)

---

## Abstract

In the existing research studies of hydraulic conductivity, most of them assume that soil consists of spheroidal particles and the value of hydraulic conductivity can be designated by the particle size. In the actual soil layers, the shape of soil particles is mostly ellipsoid or rod-like rather than ideal sphere. Therefore, the prediction of soil permeability using current method often deviates from the actual situation and cannot capture the anisotropy nature of soil without consideration of the effect of the axis size of soil particles in two different directions. To solve this problem, a new theoretical model with three different soil particle arrangements is introduced to derive a new hydraulic conductivity-particle size relationship considering the size difference in two directions. This model, from a microscopic perspective, divides pores into numerous pore units and obtains hydraulic conductivity of each tube unit, eventually predicting the permeability of soil layer based on an equivalence principle. The proposed equation is validated in comparison with experimental data from the existing literature and is proved to have a satisfied accuracy to predict hydraulic conductivity for a wide range of soils, from bentonite-silt mixed soils to sandy soils. The proposed model provides a new perspective for accurately predicting the hydraulic conductivity.


---

## Body

## 1. Introduction

Hydraulic conductivity is one of the most important and useful parameters in the study of the seepage process of porous media [1] and plays a vital role in the research of consolidation, settlement of soils and foundation, rheology [2–5], dewatering design [6–8], migration of gas and liquid [9–11], sandstone storage [12–15], storage of CO2 and nuclear waste [16–20], and other geological and geotechnical engineering problems. Generally, it is difficult to measure dynamic variation of the hydraulic conductivity by experimental means. For example, during the process of consolidation, the hydraulic conductivity gradually reduces with the decrease of void ratio. However, it is obviously difficult to measure the hydraulic conductivity of soil at each consolidation stage. Hence, many researchers propose theoretical equations to predict the change of the hydraulic conductivity. Unfortunately, there is lack of accuracy using current method due to the various simplifications and weaknesses. Thereby, it is important to propose a predictive model which can take into account essential factors that control the accuracy of the prediction of hydraulic conductivity.It has been well known that hydraulic conductivity is influenced by several factors, such as porosity [21–23], saturation [24], clay content [25], quantity of electric charge among soil particles [26], and soil particles geometry which includes grain size distribution [27–29], grain shape [30], and so on. Based on lots of experiments and theories, many theoretical equations have been proposed by studying the essence of percolation process [31–41]. These theoretical equations take into account some geotechnical parameters and, to some extent, can predict the value of hydraulic conductivity. For example, Ren et al. [36] introduced the concept of effective void ratio which can help to reduce the calculation deviation caused by relatively small pores and dead-end pores. Kozeny-Carman relation is the most common method used to predict hydraulic of soil, which is proposed by Kozeny and improved by Carman. Based on the Poiseuille law, the K-C equation reflects the possible relationship between hydraulic and void ratio:(1)k=CF1Ss2γwμρm2e31+e,where k is the hydraulic conductivity (m/s) and e represents void ratio which is dimensionless; γw is the unit weight of liquid (N/m3); CF is a dimensionless shape parameter which is estimated about 0.2 [42]; μ is the fluid viscosity (Pa⋅s); Ss is the specific surface area of soil (m2/kg); and ρm represents particle density of soil (kg/m3).The K-C equation was deduced from the assumption that the cross section of tubes where water flows is constant and regular and has been experimentally verified for coarse-grained soils such as sands. Because of extremely small pore size, irregular shape of pore, and electrochemical reactions between solid particles and liquid, the K-C relation is not appropriate for clay soils or silty soils which are composed of soil particles with small size. In addition, the K-C equation usesSs as the key parameter to reflect the influences of pore geometry on permeability rather than the particle size that cannot capture the effect of size effects [43–47]. According to the K-C relation, Wemaere [48] suggested that the notion of maximum soil particle size as well as minimum soil particle size should be introduced to account for the influence of particle size on the hydraulic conductivity of soil layer. Hence, equation (1) can be written as follows:(2)k=Ck−c∑fiDli0.404xDsi0.595−21SF2e31+e,where Ck−c is a test coefficient without dimension; fi is the volume components of particles with a certain size which is dimensionless as well; DliandDsi represent maximum soil particle and minimum soil particle size (mm), respectively; and SF is the shape coefficient. The effects of soil particle size and shape can be directly shown in this equation. However, the studied subjects of equation (2) are spherical soil particles that are not common in practical engineering. In fact, the natural soil layer is composed of countless soil particles; most of these are ellipsoidal or rod-like solids, rather than ideal spherical particles. The axial lengths of these nonspherical particles are not the same in different directions that would influence the shape and size of pores between particles. The differences of geometrical features of pores result in the soil anisotropy, especially anisotropic permeability. Hence, it is inappropriate to ignore the hydraulic anisotropy of soils by assuming perfect sphere particles of soil. Besides, equation (2) is not convenient to be used because several parameters cannot be easily obtained such as Ck−c. To obtain these parameters, a lot of experiments need to be conducted; this is time-consuming and may cause larger system error.The objective of this study is to present a theoretical method for predicting the hydraulic conductivity considering the influences of soil particles geometry from a microscopic view. Based on a series of assumptions, we link the hydraulic conductivity with the size of the major and minor axis and further propose a theoretical equation for the hydraulic conductivity that is suitable for a wide range of soils. The proposed model reflects the influence of soil particle size on hydraulic conductivity and embodies the characteristics of soil anisotropy.


## 2. Theoretical Method

In the existing research studies, many scholars approximately assumed that the actual shape of pore is equivalent to an ideal circular and then calculated the value of hydraulic conductivity by relative simplified method. However, from the microscopic view, the outer contour of pore tightly fits to the surfaces of soil particles and distributes periodically along the major axle of particles. Hence, the seepage flow pore actually is a channel of variable cross section, and it is impractical to simplify the pore as a cylinder.In addition, the size and shape of pore is also affected by the arrangement of soil particles. By means of scanning electron microscopy (SEM), Sun et al. [49] analyzed the geometrical characteristics of solid particles in sandy soils. It is found that the geometrical shape of sandy soils is closer to that of regular ellipsoids. The spatial distribution of soil particles, according to the approximate range of void ratio (0.6–1.5), is divided into 3 different groups: three-particle combination, four-particle combination, and five-particle combination (as shown in Figure 1). In spite of the same particles size, the size and shape of the pore could vary with different combinations. As a result, hydraulic conductivity may be completely different. Based on these assumptions, the paper proposes a theoretical equation to predict the hydraulic conductivity under 3 different arrangements of particles.Figure 1
Cross section of different particle systems. (a)–(c) show different cross sections with various particle combinations. (a) Cross section with three particles. (b) Cross section with four particles. (c) Cross section with five particles.
(a)
(b)
(c)


### 2.1. Theoretical Hypothesis

The assumptions of soil properties of the present model can be summarized as follows:(i)
Soils consist of even and equal-sized particles which orient in the major axle direction (Figure2). This assumption limits the spatial arrangement of particles and makes the model idealized. According to the assumption, the gradation of soil is neglected, and the hydraulic conductivity is directly determined by the size of soil particles.(ii)
Soil particles arrange periodically. The gap between particles, along the long axis, is small and negligible. Compared with the particle size, the gap size along the long axis between soil particles is too small, which shows little effect on hydraulic conductivity. In order to simplify the calculation and reduce the amount of unimportant parameters, the influence of gap is not considered in this paper.(iii)
Soil particles are ellipsoids with unequal axes length in two different directions. The size of soil particles along theZ and Y axis is equal and is noted by minor half-axis a, while the size of soil particles along the X is noted by semimajor axis b.(iv)
Soil particles are rigid bodies without deformation and displacement when fluid flows in the pores of soil particles. Hence, this paper does not consider the phenomenon of particles deformation and fragmentation and assumes that the deformation of soil layer is entirely caused by pore compression.(v)
Do not consider the effect of hydrated films surrounding surfaces of soil particles. This paper ignores relatively small pores and dead-end pores and assumes that all the pores are effective pores where free water can flow without obstructs.Figure 2
Overall distribution of soil particles. The figure mainly shows the soil layer which constitutes of four-particle systems.In addition, there are several assumptions of fluid properties as follows:(i)
The basic state of the fluid is continuous, full, and laminar.(ii)
Fluid properties (including the unit weightγw and volumetric flow rate Q) remain constant(iii)
Fluid only flows in the pore tubes which are surrounded by particle systems(iv)
The theory neglects surface effects and body forces such as gravity, centrifugal force, Coriolis force, and electromagnetic force


### 2.2. Horizontal Hydraulic Conductivitykx

This study establishes the theoretical model based on the relationship between the cross-sectional size of pore and hydraulic conductivity. Hence, along the surfaces of soil particles, hydraulic conductivity constantly changes with the variable cross-sectional size of tubes. The derivation of this chapter takes the four-particle arrangement as an example (Figures2–4). In order to clearly describe the process of derivation, some basic definitions of the pore have been introduced as shown in Figure 3. We defined the pore among a four-particle system as the pore period. The pore period can be divided into n (n⟶∞) pieces along the horizontal seepage direction. Each piece is represented as one unit tube, and length of that is designated dx. Countless connected pore periods constitute a complete microchannel which is used to denote the seepage channel between the upper and lower layer of soil particles. Therefore, the solution of this paper is as follows: firstly, based on the relative theories, this study proposes a predictive equation of hydraulic conductivity for each unit tube. Then, according to principle of flow equivalence, we can further obtain the predictive formula for the pore period. Finally, the theoretical model for predicting hydraulic conductivity of the microchannel will be established according to the periodicity of soil particles, which can show the permeability of whole soil layer.Figure 3
Various pore elements in the horizontal direction (unit tube, pore period, and microchannel).Figure 4
The cross section of four-particle system. (a) The projection of particle cross section inYZ plane. (b) The counterparts in XZ plane.
(a)
(b)Because the present model is derived from Darcy’s law, it is necessary to analyze the flow state in the pores. In this section, we seek to restrict the relationship between the geometric parameters and the velocity of the cross section in order to keep fluid conforming to laminar flow.


#### 2.2.1. Analysis of Flow State

Darcy [50] held that seepage velocity is directly proportional to hydraulic gradient under laminar flow conditions. When the flow rate increases to a certain value, flow state would change from laminar flow to turbulence. At this point, the regularity of percolation shows nonlinear characteristics and deviates from Darcy’s law. Most scholars use the Reynolds number as a critical condition for laminar; Bahrami et al. [51] suggested that the Reynolds number can be calculated by the following equation:(3)Re=ρvAμ,where ρ is the density of liquid (kg/m3); v is the average velocity of flow (m/s); A is the cross-sectional area of the tube (m2); and μ is the fluid viscosity (Pa⋅s). Many scholars agree that Darcy’s law is suitable when Re<ReΓ (ReΓ represents critical Reynolds number which usually is a certain value between 1 and 10) [52–54]. Equation (3) shows a strong relevance between the geometric properties of tube and flow state. Accordingly, the following parts of this section will analyze the influences of soil particle size on the variety of soil parameters, such as viscosity of the particle-liquid system, and eventually give the relationship between the particle size and velocity under Darcy’s law.(1) Viscosity of the Solid-Liquid System. Taking the four-particle system as an example, the geometric shape of pores in a cycle is shown in Figure 4. The square formed by the dotted line represents the maximum permeable area of the pore, while the black shaded part indicates the minimum permeable area. According to equation (3), μ represents the fluid viscosity which can show the frictional resistance during the flow of liquid. However, for the soil-liquid system, due to the existence of soil particles, μ actually indicates the viscosity of the solid-liquid system rather than that of the pure liquid. Many researchers proposed various theoretical equations related to the calculation of viscous coefficient of the solid-liquid system. Einstein [55] argued that there was a strong relevance between volume fraction of solids with viscosity coefficient based on the negligence of contact force among soil particles. On the basis of the study of Einstein, Batchelor [56] considered the interaction forces between particles and further obtained a more accurate calculation method. However, the abovementioned methods are mainly suitable for spherical soil particles which are uncommon in the practical soil layer. The method proposed by Brenner and Condiff [57] considering the ellipsoidal particles is adopted in this study, to obtain viscosity of the solid-liquid system. The viscosity coefficient can be expressed as(4)μ=μ01+ηφ+Oφ2,where μ0 is the viscous coefficient of the sole liquid, which is a parameter only related to the temperature and the pressure; η is a dimensionless variable, expressed as η=0.312ε/ln2ε−1.5+2−0.5/ln2ε−1.5−1.872/ε; φ is the volume fraction of solids; and Oφ2 represents the higher orders of φ. We ignore Oφ2 during calculation due to the fact that they have little effect on the results. Compared to viscosity coefficient of pure liquid, equation (4), considers the influences of particle size and aspect ratio on viscosity coefficient and more conforms to actual situation.(2) Cross-Sectional Areas of the Pore. For convenience of indication and calculation, we take the pore surrounded by four particles as a period and establish the Cartesian coordinate system in this pore period as shown in Figure 4(b), where the origin sets at the pore center. Hence, a complete sequence of periods can be seemed as a connected microchannel of arbitrary cross section. For one pore period, along the seepage direction of fluid, the cross-sectional radius of ellipsoid particles varies uniformly with the change of X coordinates. It can be expressed as follows:(5)rix=a2−axb2,where a and b are the sizes of half minor and major axis, respectively. The section area of the pore is determined entirely by the arrangement of soil particles and the size of soil particles. The formula of the section area of pore Ai can therefore be written as(6)Aix=4a2−πri2x=4a2−πa2−axb2.Furthermore, the average area of pore sectionAx∗ in a period will be obtained as(7)Ax∗=∫−bbAix2b=4−23πa2.(3) Critical Velocity. The abovementioned geometrical parameters of the pore are only related to soil particles. The average areas of pore sections as well as viscosity are determined when particles sizes are constant. At the same time, fluid velocity v is the only factor to affect the average Reynolds number Re of the whole microchannel. Accordingly, the relationship between v and Re can be described as the following equation:(8)Re=ρvAx∗μ=ρv4−2/3πa2μ01+ηφ+Oφ2.When the particle size and velocity satisfy certain relationship, the seepage state of fluid can adhere to Darcy’s law and the rest of the deduction in this paper will be effective and valid. Therefore, in this paper, 6 groups of soil particles with different sizes are selected to analyze the relationship between critical velocity and soil particle size. Figure5 shows the trial calculation with different b. It is found out that the size of the major and minor axis has different effects on the critical velocity. Only when v and Re satisfy a certain relationship (as shown in the shaded area of Figure 5), fluid state can be seemed as laminar.Figure 5
Relationship between particle size and critical velocity in laminar flow. The paper selects six groups with different values of major half-axisb, analyzing the changes between critical velocity and major half-axis a.


#### 2.2.2. Horizontal Hydraulic Conductivity of Unit Tubekxi

As mentioned earlier, the pore in one period is divided inton unit tubes. As a result, the derivation process of horizontal hydraulic conductivity of any unit tubes is as follows:Akbari et al. [58] proposed an equation for calculating the flow resistance R, which is the ratio of pressure drop Δp and volumetric flow rate Q, using the perturbation solution [59].(9)R=ΔpQ=16π2μ∫x1x2Ip∗A2dx︸frictional+ρQ1A22−1A12︸inertial,where Ip∗=Ip/A2 with Ip=∫y2+z2dA is designated as the specific polar moment of cross-sectional inertia [51] and A1 and A2 are the pore cross-sectional area at x1,x2 locations, respectively. The first item on the right side of the equation represents the flow resistance due to friction and the second represents inertia term. The inertial pressure drop in periodic channel will be negligible for low Reynolds numbers. Hence, equation (9) is reformulated as(10)R=16π2μ∫x1x2Ip∗A2dx.Equation (10) can be rewritten in the differential form. Accordingly, the flow resistance equation of the arbitrary unit tube dRi can be furthermore expressed as the following equation:(11)dRi=dΔpiQ=16π2μiIpiAi4dx.Then, the volumetric flow rate can be expressed as(12)Q=dpi16π2μiIpi/Ai4dx.Thus, according to Darcy’s law, the rate of flow in the unit tube could be indicated in the following equation, whereAti is the total area of the cross section including the tube walls:(13)Q=kxiiiAti,where along the direction of seepage flow, the hydraulic gradient can be described as(14)ii=dpiγwdx.Combining equations (12) and (14), the horizontal hydraulic conductivity of arbitrary unit tube kix can be modified as follows:(15)kxi=γw16π2μiAtiIpi/Ai4.Equation (15) shows that kxi is mainly related to three parameters, namely, Ai, Ati, and Ipi. When the soil layer is entirely composed of four-particle systems as shown in Figure 2, equation (6) can be used to calculate the area of arbitrary cross section of unit tubes Ai, while Ati can be expressed by the sum of Ai and Asi which means the cross-sectional areas of tube walls:(16)Ati=Ai+Asi=4a2+3πa2−axb2.In addition, the polar moment of inertiaIpi can be expressed as(17)Ipix=8a43−π2−329πa2−axb24−2πa2−axb22a−4a2−ax/b23π2.Thus, substituting equations (6), (16), and (17) into equation (15) yields(18)kxi=γw16π2μi4a2+4−πa2−ax/b28a4/3−π/2−32/9πa2−ax/b24−2πa2−ax/b22a−4a2−ax/b2/3π2/4a2−πri2x=4a2−πa2−ax/b24.For the sake of simplification, equation (15) is used for calculation in this study instead of equation (18). Based on the above derivation, the horizontal hydraulic conductivity of arbitrary unit tube kxi can be predicted by the major axis and minor axis size of soil particles. Combining with the equivalent flow method which will be described in the next Section 2.2.3, we can get the horizontal hydraulic conductivity of whole soil layer kx.


#### 2.2.3. Equivalent Horizontal Hydraulic Conductivitykx∗

Based on the equivalent flow method, we propose a formula to compute the whole horizontal hydraulic conductivitykx∗ which can reflect horizontal permeability of soil.Different pore elements are parallel to each other in seepage direction as shown in Figure6. In addition, the volumetric flow rate of each pore unit is equal. The velocity of flow, at xi location, can be expressed as follows:(19)vi=QAi=kxidΔpiγwdx.Figure 6
Distribution of unit tube in one period. The pore in one period is divided into countless pore units with the same widthdx.Thus,(20)dpi=Q/Aidxkxi.In a similar way, the pressure drop at both ends of the microchannel can be defined as(21)Δp=Q/Ax∗Lkx∗,where Δp= the pressure drop at both ends of the microchannel (Pa); L = channel length (m); and Ax∗= the equivalent cross-sectional area of the channel, which is represented by average cross-sectional area as equation (7). The relationship between the pressure drop of whole channel and unit tube can be expressed as follows:(22)Δp=∑dpi.Substituting equations (20) and (21) into equation (22), an equivalent formula of the horizontal hydraulic conductivity kx∗ can be obtained:(23)kx∗=LAx∗∑dx/Aixkxi=Lγw16π2Ax∗∫0LμiAtiIpi/Ai5dx.To implement the proposed equation in practice, a further simplification of equation (23) is suggested to calculate horizontal hydraulic conductivity. First, it assumes ∫0LμiIpi/Ai4dx=μ∫0LIpi/Ai4dx, where μ means the viscosity of the solid-liquid system, which can be obtained from equation (4). When the total length of microchannel is designated L, according to the theoretical hypothesis, there are m m=L/2b periods of pore. Hence, microchannel L can be expressed as L=2mb. Moreover, ∫0LμiAtiIpi/Aix5dx=∫02mbμiAtiIpi/Aix5dx. Due to the periodicity of the microchannel and the parity of μiAtiIpi/Aix5, we can obtain that ∫02mbμiAtiIpi/Aix5dx=2m∫0bμiAtiIpi/Aix5dx. Hence, equation (23) can be rewritten as(24)kx∗=bγw8π2Ax∗∫0LμiAtiIpi/Aix5dx.The equivalent hydraulic conductivity of the whole soil layer in horizontal directionkx∗ is expressed as equation (24). In practice, kx∗ can be regarded as the horizontal hydraulic conductivity of whole layer kx in order to reflect seepage law. Various microfactors such as particle shape and particle size distribution are incorporated in equation (4). Therefore, the proposed equation can well capture the microscopic mechanism of soil seepage.


#### 2.2.4. Vertical Hydraulic Conductivitykz

The derivation process of formula of the vertical hydraulic conductivity is the same as the horizontal counterpart in this paper. However, compared with horizontal pores, the shape and size of vertical microchannels are quite different as shown in Figure7. The shaded area and dotted rectangle represent the minimum seepage area and the maximum seepage area, respectively, in the vertical direction. The cross-sectional area of pore An can be expressed as the following equation:(25)Anz=4ab−πYzXz,where Yz=z, Xz=b2−b/az2.Figure 7
Projection of vertical seepage tube inXY plane.Vertical average sectional area of pore can be obtained as follows:(26)Az∗=∫−aaAnzdz2a=4−π3ab.Then, in the vertical direction, the area of the cross section of tube walls can be computed as(27)Asn=4πYzXz.UsingAsn and An, we can further obtain total area of the cross section including the tube walls which is designated Atn:(28)Atn=An+Asn=4ab+3πYzXz.Similarly, the polar moment of inertia of arbitrary vertical unit tubeIpn can be computed as(29)Ipn=43a3b+b3a−π16−49πYzX3z−π4b−4Xz3π2YzXz−π16−49πY3zXz−π4a−4Yz3π2YzXz.Based on the same equivalence principle, the equivalent hydraulic conductivity of the whole soil layer in vertical direction is given by(30)kz∗=aγw8π2Az∗∫0LμiAtnIpn/An5dx.


### 2.3. Other Particle Combinations

According to the hypothesis in Section2.1, the soil layer is composed of numerous “particle systems.” These particle systems can be divided into three-particle combination, four-particle combination, and five-particle combination according to void ratio e. Similar to the proposed method for four-particle combination, the computation equation of hydraulic conductivity for soil under two other combinations can be obtained. However, the arrangement of particles for three or five particles in the vertical plane is not as regular as that of four particles. For three-particle system as well as five-particle system, in the YZ plane, there is a certain gap between the upper and lower rows of soil particles, which means the seepage route is not connected along the Z axis as shown in Figure 8. As a result, the calculation method proposed in this paper is not applicable for calculating the vertical hydraulic conductivity of these arrangements. However, for the sake of clarity, we did not deduce new formulas for the other two cases. Equation (30) is recommended to calculate the vertical hydraulic conductivity for soil with all particle systems. This means the effect of soil particle arrangement on hydraulic conductivity is ignored during calculation; only the size of soil particles is taken into consideration. As a result, the size of pore, especially for three-particle combination, could be overestimated, and the predicted results will be inaccurate. Therefore, the vertical calculation results for three-particle and five-particle combination would merely be the reference resources. More precise calculation methods still need to be further studied.Figure 8
Projection of three-particle systems inYZ plane.


#### 2.3.1. Three-Particle Arrangement

The soil layer consisting of three-particle system usually characterizes of small pores and low hydraulic conductivity. As shown in Figure9(a), the shaded part represents the minimum seepage area and the dashed triangle indicates the maximum pore area. Under this condition, the cross-sectional area Ai and the polar moment of inertia can be expressed as(31)Aix=3a2−π2a2−axb2,Ipix=136b4πa3403b2πb−xb+xa21−x2b2+ab442+143π−33π2+42b2−2+π2x2+314−3π2x4.Figure 9
Other particle combinations. (a) The three-particle system. (b) The five-particle system.
(a)
(b)Furthermore,Ax∗,Ati can be obtained as(32)Ax∗=∫−bbAixdx2b=33−π3a2,Ati=Ai+Asi=3a2+2.5πa2−axb2.Finally, substituting formulas ofAi, Ipi, and Ax∗ into equation (24), the equivalent hydraulic conductivity in horizontal direction with three particles can be obtained.


#### 2.3.2. Five-Particle Arrangement

The size of pore formed by five particles is relatively large. The maximum seepage cross-sectional area can be simplified as regular pentagon while the minimum area is surrounded by five tangent circles, as the shaded area shown in Figure9(b). Similarly, the cross-sectional area Ai and the polar moment of inertia Ipi are given by(33)Aix=6.89a2−1.5πa2−axb2,Ipix=a3α0a21−x2/b2b4−b2x2+aα1b4+α2b2x2+α3x4b4,where α0=−5.489,α1=16.082,α2=−2.736,α3=−2.726.Thus, equation (33) yields(34)Ax∗=∫−bbAixdx2b=5.32a2,(35)Ati=Ai+Asi=6.89a2+3.5πa2−axb2.Combining equations (33)–(36), the equivalent hydraulic conductivity of whole soil layer in horizontal direction can be obtained for five-particle system.


### 2.4. Hydraulic Conductivity vs Void Ratioe

In addition to particle size, void ratio is another essential factor influencing the permeability of soil. In practice, void ratio will constantly change with the process of consolidation; as a result, the hydraulic conductivity will also vary with the change of void ratio. Based on the large number of experiments, several theoretical equations have been proposed to describe the relationship between void ratio and hydraulic conductivity. For example, Lambe and Whitman [60] argued that there is a logarithmic relationship between the hydraulic conductivity and the void ratio:(36)e=e0+Cklogkk0,where Ck denotes the coefficient of permeability change and e0 and k0 denote the initial void ratio and the initial hydraulic conductivity, respectively. This equation only considers the relationship between hydraulic conductivity and void ratio and cannot reflect the coupling effect of particle size and void ratio from the microscopic and macro perspective, respectively.This section takes the four-particle combination as an example and further derives the relation between hydraulic conductivity and pore ratio. As shown in Figure10, in the initial state, the upper and lower particles are not closely connected, and there is a certain size of gap designated as l0 which is determined by the void ratio e of soil. During the process of consolidation, the gap between particles is gradually compressed. As a result, hydraulic conductivity changed with the decline of l0 and e.Figure 10
The change of gap during the process of consolidation.Hence, the area of arbitrary cross section of unit tubesAi can be rewritten as(37)Ai=2a+l02−πa2−axb2.Furthermore, the cross-sectional areas of tube wallsAsi can be expressed as follows:(38)Asi=πa2−axb2.According to the assumption that the soil layer is composed of innumerable particle arrangement systems, void ratio in four-particle combination can represent the void ratio of the whole soil layer. Thus,(39)e=∫−bbAidx∫−bbAsidx=32a+l02−2πa22πa2.Then, the size of gapl0 can be expressed by void ratio e:(40)l0=2π1+e3−2a.Thus, the equivalent cross-sectional area of the microchannelAx∗ and the polar moment of inertia Ipi can be rewritten as(41)Ax∗=∫−bbAidx2b=4al+l2+236−πa2=4a2π1+e3−2a+2π1+e3−2a2+236−πa2,(42)Ipix=2a+l046−π2−329πa2−axb22−2πa2−axb22a+l2−4a2−ax/b23π2=2π1+e227−π2−329πa2−axb22−2πa2−axb22π1+e/32−4a2−ax/b23π2.Substituting equations (41) and (42) into equation (24), we can finally obtain the predictive equation of horizontal hydraulic conductivity considering the change of void ratio. According to the same principle, we can further derive the theoretical model for vertical hydraulic conductivity and other particle combinations.


## 2.1. Theoretical Hypothesis

The assumptions of soil properties of the present model can be summarized as follows:(i)
Soils consist of even and equal-sized particles which orient in the major axle direction (Figure2). This assumption limits the spatial arrangement of particles and makes the model idealized. According to the assumption, the gradation of soil is neglected, and the hydraulic conductivity is directly determined by the size of soil particles.(ii)
Soil particles arrange periodically. The gap between particles, along the long axis, is small and negligible. Compared with the particle size, the gap size along the long axis between soil particles is too small, which shows little effect on hydraulic conductivity. In order to simplify the calculation and reduce the amount of unimportant parameters, the influence of gap is not considered in this paper.(iii)
Soil particles are ellipsoids with unequal axes length in two different directions. The size of soil particles along theZ and Y axis is equal and is noted by minor half-axis a, while the size of soil particles along the X is noted by semimajor axis b.(iv)
Soil particles are rigid bodies without deformation and displacement when fluid flows in the pores of soil particles. Hence, this paper does not consider the phenomenon of particles deformation and fragmentation and assumes that the deformation of soil layer is entirely caused by pore compression.(v)
Do not consider the effect of hydrated films surrounding surfaces of soil particles. This paper ignores relatively small pores and dead-end pores and assumes that all the pores are effective pores where free water can flow without obstructs.Figure 2
Overall distribution of soil particles. The figure mainly shows the soil layer which constitutes of four-particle systems.In addition, there are several assumptions of fluid properties as follows:(i)
The basic state of the fluid is continuous, full, and laminar.(ii)
Fluid properties (including the unit weightγw and volumetric flow rate Q) remain constant(iii)
Fluid only flows in the pore tubes which are surrounded by particle systems(iv)
The theory neglects surface effects and body forces such as gravity, centrifugal force, Coriolis force, and electromagnetic force


## 2.2. Horizontal Hydraulic Conductivitykx

This study establishes the theoretical model based on the relationship between the cross-sectional size of pore and hydraulic conductivity. Hence, along the surfaces of soil particles, hydraulic conductivity constantly changes with the variable cross-sectional size of tubes. The derivation of this chapter takes the four-particle arrangement as an example (Figures2–4). In order to clearly describe the process of derivation, some basic definitions of the pore have been introduced as shown in Figure 3. We defined the pore among a four-particle system as the pore period. The pore period can be divided into n (n⟶∞) pieces along the horizontal seepage direction. Each piece is represented as one unit tube, and length of that is designated dx. Countless connected pore periods constitute a complete microchannel which is used to denote the seepage channel between the upper and lower layer of soil particles. Therefore, the solution of this paper is as follows: firstly, based on the relative theories, this study proposes a predictive equation of hydraulic conductivity for each unit tube. Then, according to principle of flow equivalence, we can further obtain the predictive formula for the pore period. Finally, the theoretical model for predicting hydraulic conductivity of the microchannel will be established according to the periodicity of soil particles, which can show the permeability of whole soil layer.Figure 3
Various pore elements in the horizontal direction (unit tube, pore period, and microchannel).Figure 4
The cross section of four-particle system. (a) The projection of particle cross section inYZ plane. (b) The counterparts in XZ plane.
(a)
(b)Because the present model is derived from Darcy’s law, it is necessary to analyze the flow state in the pores. In this section, we seek to restrict the relationship between the geometric parameters and the velocity of the cross section in order to keep fluid conforming to laminar flow.


### 2.2.1. Analysis of Flow State

Darcy [50] held that seepage velocity is directly proportional to hydraulic gradient under laminar flow conditions. When the flow rate increases to a certain value, flow state would change from laminar flow to turbulence. At this point, the regularity of percolation shows nonlinear characteristics and deviates from Darcy’s law. Most scholars use the Reynolds number as a critical condition for laminar; Bahrami et al. [51] suggested that the Reynolds number can be calculated by the following equation:(3)Re=ρvAμ,where ρ is the density of liquid (kg/m3); v is the average velocity of flow (m/s); A is the cross-sectional area of the tube (m2); and μ is the fluid viscosity (Pa⋅s). Many scholars agree that Darcy’s law is suitable when Re<ReΓ (ReΓ represents critical Reynolds number which usually is a certain value between 1 and 10) [52–54]. Equation (3) shows a strong relevance between the geometric properties of tube and flow state. Accordingly, the following parts of this section will analyze the influences of soil particle size on the variety of soil parameters, such as viscosity of the particle-liquid system, and eventually give the relationship between the particle size and velocity under Darcy’s law.(1) Viscosity of the Solid-Liquid System. Taking the four-particle system as an example, the geometric shape of pores in a cycle is shown in Figure 4. The square formed by the dotted line represents the maximum permeable area of the pore, while the black shaded part indicates the minimum permeable area. According to equation (3), μ represents the fluid viscosity which can show the frictional resistance during the flow of liquid. However, for the soil-liquid system, due to the existence of soil particles, μ actually indicates the viscosity of the solid-liquid system rather than that of the pure liquid. Many researchers proposed various theoretical equations related to the calculation of viscous coefficient of the solid-liquid system. Einstein [55] argued that there was a strong relevance between volume fraction of solids with viscosity coefficient based on the negligence of contact force among soil particles. On the basis of the study of Einstein, Batchelor [56] considered the interaction forces between particles and further obtained a more accurate calculation method. However, the abovementioned methods are mainly suitable for spherical soil particles which are uncommon in the practical soil layer. The method proposed by Brenner and Condiff [57] considering the ellipsoidal particles is adopted in this study, to obtain viscosity of the solid-liquid system. The viscosity coefficient can be expressed as(4)μ=μ01+ηφ+Oφ2,where μ0 is the viscous coefficient of the sole liquid, which is a parameter only related to the temperature and the pressure; η is a dimensionless variable, expressed as η=0.312ε/ln2ε−1.5+2−0.5/ln2ε−1.5−1.872/ε; φ is the volume fraction of solids; and Oφ2 represents the higher orders of φ. We ignore Oφ2 during calculation due to the fact that they have little effect on the results. Compared to viscosity coefficient of pure liquid, equation (4), considers the influences of particle size and aspect ratio on viscosity coefficient and more conforms to actual situation.(2) Cross-Sectional Areas of the Pore. For convenience of indication and calculation, we take the pore surrounded by four particles as a period and establish the Cartesian coordinate system in this pore period as shown in Figure 4(b), where the origin sets at the pore center. Hence, a complete sequence of periods can be seemed as a connected microchannel of arbitrary cross section. For one pore period, along the seepage direction of fluid, the cross-sectional radius of ellipsoid particles varies uniformly with the change of X coordinates. It can be expressed as follows:(5)rix=a2−axb2,where a and b are the sizes of half minor and major axis, respectively. The section area of the pore is determined entirely by the arrangement of soil particles and the size of soil particles. The formula of the section area of pore Ai can therefore be written as(6)Aix=4a2−πri2x=4a2−πa2−axb2.Furthermore, the average area of pore sectionAx∗ in a period will be obtained as(7)Ax∗=∫−bbAix2b=4−23πa2.(3) Critical Velocity. The abovementioned geometrical parameters of the pore are only related to soil particles. The average areas of pore sections as well as viscosity are determined when particles sizes are constant. At the same time, fluid velocity v is the only factor to affect the average Reynolds number Re of the whole microchannel. Accordingly, the relationship between v and Re can be described as the following equation:(8)Re=ρvAx∗μ=ρv4−2/3πa2μ01+ηφ+Oφ2.When the particle size and velocity satisfy certain relationship, the seepage state of fluid can adhere to Darcy’s law and the rest of the deduction in this paper will be effective and valid. Therefore, in this paper, 6 groups of soil particles with different sizes are selected to analyze the relationship between critical velocity and soil particle size. Figure5 shows the trial calculation with different b. It is found out that the size of the major and minor axis has different effects on the critical velocity. Only when v and Re satisfy a certain relationship (as shown in the shaded area of Figure 5), fluid state can be seemed as laminar.Figure 5
Relationship between particle size and critical velocity in laminar flow. The paper selects six groups with different values of major half-axisb, analyzing the changes between critical velocity and major half-axis a.


### 2.2.2. Horizontal Hydraulic Conductivity of Unit Tubekxi

As mentioned earlier, the pore in one period is divided inton unit tubes. As a result, the derivation process of horizontal hydraulic conductivity of any unit tubes is as follows:Akbari et al. [58] proposed an equation for calculating the flow resistance R, which is the ratio of pressure drop Δp and volumetric flow rate Q, using the perturbation solution [59].(9)R=ΔpQ=16π2μ∫x1x2Ip∗A2dx︸frictional+ρQ1A22−1A12︸inertial,where Ip∗=Ip/A2 with Ip=∫y2+z2dA is designated as the specific polar moment of cross-sectional inertia [51] and A1 and A2 are the pore cross-sectional area at x1,x2 locations, respectively. The first item on the right side of the equation represents the flow resistance due to friction and the second represents inertia term. The inertial pressure drop in periodic channel will be negligible for low Reynolds numbers. Hence, equation (9) is reformulated as(10)R=16π2μ∫x1x2Ip∗A2dx.Equation (10) can be rewritten in the differential form. Accordingly, the flow resistance equation of the arbitrary unit tube dRi can be furthermore expressed as the following equation:(11)dRi=dΔpiQ=16π2μiIpiAi4dx.Then, the volumetric flow rate can be expressed as(12)Q=dpi16π2μiIpi/Ai4dx.Thus, according to Darcy’s law, the rate of flow in the unit tube could be indicated in the following equation, whereAti is the total area of the cross section including the tube walls:(13)Q=kxiiiAti,where along the direction of seepage flow, the hydraulic gradient can be described as(14)ii=dpiγwdx.Combining equations (12) and (14), the horizontal hydraulic conductivity of arbitrary unit tube kix can be modified as follows:(15)kxi=γw16π2μiAtiIpi/Ai4.Equation (15) shows that kxi is mainly related to three parameters, namely, Ai, Ati, and Ipi. When the soil layer is entirely composed of four-particle systems as shown in Figure 2, equation (6) can be used to calculate the area of arbitrary cross section of unit tubes Ai, while Ati can be expressed by the sum of Ai and Asi which means the cross-sectional areas of tube walls:(16)Ati=Ai+Asi=4a2+3πa2−axb2.In addition, the polar moment of inertiaIpi can be expressed as(17)Ipix=8a43−π2−329πa2−axb24−2πa2−axb22a−4a2−ax/b23π2.Thus, substituting equations (6), (16), and (17) into equation (15) yields(18)kxi=γw16π2μi4a2+4−πa2−ax/b28a4/3−π/2−32/9πa2−ax/b24−2πa2−ax/b22a−4a2−ax/b2/3π2/4a2−πri2x=4a2−πa2−ax/b24.For the sake of simplification, equation (15) is used for calculation in this study instead of equation (18). Based on the above derivation, the horizontal hydraulic conductivity of arbitrary unit tube kxi can be predicted by the major axis and minor axis size of soil particles. Combining with the equivalent flow method which will be described in the next Section 2.2.3, we can get the horizontal hydraulic conductivity of whole soil layer kx.


### 2.2.3. Equivalent Horizontal Hydraulic Conductivitykx∗

Based on the equivalent flow method, we propose a formula to compute the whole horizontal hydraulic conductivitykx∗ which can reflect horizontal permeability of soil.Different pore elements are parallel to each other in seepage direction as shown in Figure6. In addition, the volumetric flow rate of each pore unit is equal. The velocity of flow, at xi location, can be expressed as follows:(19)vi=QAi=kxidΔpiγwdx.Figure 6
Distribution of unit tube in one period. The pore in one period is divided into countless pore units with the same widthdx.Thus,(20)dpi=Q/Aidxkxi.In a similar way, the pressure drop at both ends of the microchannel can be defined as(21)Δp=Q/Ax∗Lkx∗,where Δp= the pressure drop at both ends of the microchannel (Pa); L = channel length (m); and Ax∗= the equivalent cross-sectional area of the channel, which is represented by average cross-sectional area as equation (7). The relationship between the pressure drop of whole channel and unit tube can be expressed as follows:(22)Δp=∑dpi.Substituting equations (20) and (21) into equation (22), an equivalent formula of the horizontal hydraulic conductivity kx∗ can be obtained:(23)kx∗=LAx∗∑dx/Aixkxi=Lγw16π2Ax∗∫0LμiAtiIpi/Ai5dx.To implement the proposed equation in practice, a further simplification of equation (23) is suggested to calculate horizontal hydraulic conductivity. First, it assumes ∫0LμiIpi/Ai4dx=μ∫0LIpi/Ai4dx, where μ means the viscosity of the solid-liquid system, which can be obtained from equation (4). When the total length of microchannel is designated L, according to the theoretical hypothesis, there are m m=L/2b periods of pore. Hence, microchannel L can be expressed as L=2mb. Moreover, ∫0LμiAtiIpi/Aix5dx=∫02mbμiAtiIpi/Aix5dx. Due to the periodicity of the microchannel and the parity of μiAtiIpi/Aix5, we can obtain that ∫02mbμiAtiIpi/Aix5dx=2m∫0bμiAtiIpi/Aix5dx. Hence, equation (23) can be rewritten as(24)kx∗=bγw8π2Ax∗∫0LμiAtiIpi/Aix5dx.The equivalent hydraulic conductivity of the whole soil layer in horizontal directionkx∗ is expressed as equation (24). In practice, kx∗ can be regarded as the horizontal hydraulic conductivity of whole layer kx in order to reflect seepage law. Various microfactors such as particle shape and particle size distribution are incorporated in equation (4). Therefore, the proposed equation can well capture the microscopic mechanism of soil seepage.


### 2.2.4. Vertical Hydraulic Conductivitykz

The derivation process of formula of the vertical hydraulic conductivity is the same as the horizontal counterpart in this paper. However, compared with horizontal pores, the shape and size of vertical microchannels are quite different as shown in Figure7. The shaded area and dotted rectangle represent the minimum seepage area and the maximum seepage area, respectively, in the vertical direction. The cross-sectional area of pore An can be expressed as the following equation:(25)Anz=4ab−πYzXz,where Yz=z, Xz=b2−b/az2.Figure 7
Projection of vertical seepage tube inXY plane.Vertical average sectional area of pore can be obtained as follows:(26)Az∗=∫−aaAnzdz2a=4−π3ab.Then, in the vertical direction, the area of the cross section of tube walls can be computed as(27)Asn=4πYzXz.UsingAsn and An, we can further obtain total area of the cross section including the tube walls which is designated Atn:(28)Atn=An+Asn=4ab+3πYzXz.Similarly, the polar moment of inertia of arbitrary vertical unit tubeIpn can be computed as(29)Ipn=43a3b+b3a−π16−49πYzX3z−π4b−4Xz3π2YzXz−π16−49πY3zXz−π4a−4Yz3π2YzXz.Based on the same equivalence principle, the equivalent hydraulic conductivity of the whole soil layer in vertical direction is given by(30)kz∗=aγw8π2Az∗∫0LμiAtnIpn/An5dx.


## 2.2.1. Analysis of Flow State

Darcy [50] held that seepage velocity is directly proportional to hydraulic gradient under laminar flow conditions. When the flow rate increases to a certain value, flow state would change from laminar flow to turbulence. At this point, the regularity of percolation shows nonlinear characteristics and deviates from Darcy’s law. Most scholars use the Reynolds number as a critical condition for laminar; Bahrami et al. [51] suggested that the Reynolds number can be calculated by the following equation:(3)Re=ρvAμ,where ρ is the density of liquid (kg/m3); v is the average velocity of flow (m/s); A is the cross-sectional area of the tube (m2); and μ is the fluid viscosity (Pa⋅s). Many scholars agree that Darcy’s law is suitable when Re<ReΓ (ReΓ represents critical Reynolds number which usually is a certain value between 1 and 10) [52–54]. Equation (3) shows a strong relevance between the geometric properties of tube and flow state. Accordingly, the following parts of this section will analyze the influences of soil particle size on the variety of soil parameters, such as viscosity of the particle-liquid system, and eventually give the relationship between the particle size and velocity under Darcy’s law.(1) Viscosity of the Solid-Liquid System. Taking the four-particle system as an example, the geometric shape of pores in a cycle is shown in Figure 4. The square formed by the dotted line represents the maximum permeable area of the pore, while the black shaded part indicates the minimum permeable area. According to equation (3), μ represents the fluid viscosity which can show the frictional resistance during the flow of liquid. However, for the soil-liquid system, due to the existence of soil particles, μ actually indicates the viscosity of the solid-liquid system rather than that of the pure liquid. Many researchers proposed various theoretical equations related to the calculation of viscous coefficient of the solid-liquid system. Einstein [55] argued that there was a strong relevance between volume fraction of solids with viscosity coefficient based on the negligence of contact force among soil particles. On the basis of the study of Einstein, Batchelor [56] considered the interaction forces between particles and further obtained a more accurate calculation method. However, the abovementioned methods are mainly suitable for spherical soil particles which are uncommon in the practical soil layer. The method proposed by Brenner and Condiff [57] considering the ellipsoidal particles is adopted in this study, to obtain viscosity of the solid-liquid system. The viscosity coefficient can be expressed as(4)μ=μ01+ηφ+Oφ2,where μ0 is the viscous coefficient of the sole liquid, which is a parameter only related to the temperature and the pressure; η is a dimensionless variable, expressed as η=0.312ε/ln2ε−1.5+2−0.5/ln2ε−1.5−1.872/ε; φ is the volume fraction of solids; and Oφ2 represents the higher orders of φ. We ignore Oφ2 during calculation due to the fact that they have little effect on the results. Compared to viscosity coefficient of pure liquid, equation (4), considers the influences of particle size and aspect ratio on viscosity coefficient and more conforms to actual situation.(2) Cross-Sectional Areas of the Pore. For convenience of indication and calculation, we take the pore surrounded by four particles as a period and establish the Cartesian coordinate system in this pore period as shown in Figure 4(b), where the origin sets at the pore center. Hence, a complete sequence of periods can be seemed as a connected microchannel of arbitrary cross section. For one pore period, along the seepage direction of fluid, the cross-sectional radius of ellipsoid particles varies uniformly with the change of X coordinates. It can be expressed as follows:(5)rix=a2−axb2,where a and b are the sizes of half minor and major axis, respectively. The section area of the pore is determined entirely by the arrangement of soil particles and the size of soil particles. The formula of the section area of pore Ai can therefore be written as(6)Aix=4a2−πri2x=4a2−πa2−axb2.Furthermore, the average area of pore sectionAx∗ in a period will be obtained as(7)Ax∗=∫−bbAix2b=4−23πa2.(3) Critical Velocity. The abovementioned geometrical parameters of the pore are only related to soil particles. The average areas of pore sections as well as viscosity are determined when particles sizes are constant. At the same time, fluid velocity v is the only factor to affect the average Reynolds number Re of the whole microchannel. Accordingly, the relationship between v and Re can be described as the following equation:(8)Re=ρvAx∗μ=ρv4−2/3πa2μ01+ηφ+Oφ2.When the particle size and velocity satisfy certain relationship, the seepage state of fluid can adhere to Darcy’s law and the rest of the deduction in this paper will be effective and valid. Therefore, in this paper, 6 groups of soil particles with different sizes are selected to analyze the relationship between critical velocity and soil particle size. Figure5 shows the trial calculation with different b. It is found out that the size of the major and minor axis has different effects on the critical velocity. Only when v and Re satisfy a certain relationship (as shown in the shaded area of Figure 5), fluid state can be seemed as laminar.Figure 5
Relationship between particle size and critical velocity in laminar flow. The paper selects six groups with different values of major half-axisb, analyzing the changes between critical velocity and major half-axis a.


## 2.2.2. Horizontal Hydraulic Conductivity of Unit Tubekxi

As mentioned earlier, the pore in one period is divided inton unit tubes. As a result, the derivation process of horizontal hydraulic conductivity of any unit tubes is as follows:Akbari et al. [58] proposed an equation for calculating the flow resistance R, which is the ratio of pressure drop Δp and volumetric flow rate Q, using the perturbation solution [59].(9)R=ΔpQ=16π2μ∫x1x2Ip∗A2dx︸frictional+ρQ1A22−1A12︸inertial,where Ip∗=Ip/A2 with Ip=∫y2+z2dA is designated as the specific polar moment of cross-sectional inertia [51] and A1 and A2 are the pore cross-sectional area at x1,x2 locations, respectively. The first item on the right side of the equation represents the flow resistance due to friction and the second represents inertia term. The inertial pressure drop in periodic channel will be negligible for low Reynolds numbers. Hence, equation (9) is reformulated as(10)R=16π2μ∫x1x2Ip∗A2dx.Equation (10) can be rewritten in the differential form. Accordingly, the flow resistance equation of the arbitrary unit tube dRi can be furthermore expressed as the following equation:(11)dRi=dΔpiQ=16π2μiIpiAi4dx.Then, the volumetric flow rate can be expressed as(12)Q=dpi16π2μiIpi/Ai4dx.Thus, according to Darcy’s law, the rate of flow in the unit tube could be indicated in the following equation, whereAti is the total area of the cross section including the tube walls:(13)Q=kxiiiAti,where along the direction of seepage flow, the hydraulic gradient can be described as(14)ii=dpiγwdx.Combining equations (12) and (14), the horizontal hydraulic conductivity of arbitrary unit tube kix can be modified as follows:(15)kxi=γw16π2μiAtiIpi/Ai4.Equation (15) shows that kxi is mainly related to three parameters, namely, Ai, Ati, and Ipi. When the soil layer is entirely composed of four-particle systems as shown in Figure 2, equation (6) can be used to calculate the area of arbitrary cross section of unit tubes Ai, while Ati can be expressed by the sum of Ai and Asi which means the cross-sectional areas of tube walls:(16)Ati=Ai+Asi=4a2+3πa2−axb2.In addition, the polar moment of inertiaIpi can be expressed as(17)Ipix=8a43−π2−329πa2−axb24−2πa2−axb22a−4a2−ax/b23π2.Thus, substituting equations (6), (16), and (17) into equation (15) yields(18)kxi=γw16π2μi4a2+4−πa2−ax/b28a4/3−π/2−32/9πa2−ax/b24−2πa2−ax/b22a−4a2−ax/b2/3π2/4a2−πri2x=4a2−πa2−ax/b24.For the sake of simplification, equation (15) is used for calculation in this study instead of equation (18). Based on the above derivation, the horizontal hydraulic conductivity of arbitrary unit tube kxi can be predicted by the major axis and minor axis size of soil particles. Combining with the equivalent flow method which will be described in the next Section 2.2.3, we can get the horizontal hydraulic conductivity of whole soil layer kx.


## 2.2.3. Equivalent Horizontal Hydraulic Conductivitykx∗

Based on the equivalent flow method, we propose a formula to compute the whole horizontal hydraulic conductivitykx∗ which can reflect horizontal permeability of soil.Different pore elements are parallel to each other in seepage direction as shown in Figure6. In addition, the volumetric flow rate of each pore unit is equal. The velocity of flow, at xi location, can be expressed as follows:(19)vi=QAi=kxidΔpiγwdx.Figure 6
Distribution of unit tube in one period. The pore in one period is divided into countless pore units with the same widthdx.Thus,(20)dpi=Q/Aidxkxi.In a similar way, the pressure drop at both ends of the microchannel can be defined as(21)Δp=Q/Ax∗Lkx∗,where Δp= the pressure drop at both ends of the microchannel (Pa); L = channel length (m); and Ax∗= the equivalent cross-sectional area of the channel, which is represented by average cross-sectional area as equation (7). The relationship between the pressure drop of whole channel and unit tube can be expressed as follows:(22)Δp=∑dpi.Substituting equations (20) and (21) into equation (22), an equivalent formula of the horizontal hydraulic conductivity kx∗ can be obtained:(23)kx∗=LAx∗∑dx/Aixkxi=Lγw16π2Ax∗∫0LμiAtiIpi/Ai5dx.To implement the proposed equation in practice, a further simplification of equation (23) is suggested to calculate horizontal hydraulic conductivity. First, it assumes ∫0LμiIpi/Ai4dx=μ∫0LIpi/Ai4dx, where μ means the viscosity of the solid-liquid system, which can be obtained from equation (4). When the total length of microchannel is designated L, according to the theoretical hypothesis, there are m m=L/2b periods of pore. Hence, microchannel L can be expressed as L=2mb. Moreover, ∫0LμiAtiIpi/Aix5dx=∫02mbμiAtiIpi/Aix5dx. Due to the periodicity of the microchannel and the parity of μiAtiIpi/Aix5, we can obtain that ∫02mbμiAtiIpi/Aix5dx=2m∫0bμiAtiIpi/Aix5dx. Hence, equation (23) can be rewritten as(24)kx∗=bγw8π2Ax∗∫0LμiAtiIpi/Aix5dx.The equivalent hydraulic conductivity of the whole soil layer in horizontal directionkx∗ is expressed as equation (24). In practice, kx∗ can be regarded as the horizontal hydraulic conductivity of whole layer kx in order to reflect seepage law. Various microfactors such as particle shape and particle size distribution are incorporated in equation (4). Therefore, the proposed equation can well capture the microscopic mechanism of soil seepage.


## 2.2.4. Vertical Hydraulic Conductivitykz

The derivation process of formula of the vertical hydraulic conductivity is the same as the horizontal counterpart in this paper. However, compared with horizontal pores, the shape and size of vertical microchannels are quite different as shown in Figure7. The shaded area and dotted rectangle represent the minimum seepage area and the maximum seepage area, respectively, in the vertical direction. The cross-sectional area of pore An can be expressed as the following equation:(25)Anz=4ab−πYzXz,where Yz=z, Xz=b2−b/az2.Figure 7
Projection of vertical seepage tube inXY plane.Vertical average sectional area of pore can be obtained as follows:(26)Az∗=∫−aaAnzdz2a=4−π3ab.Then, in the vertical direction, the area of the cross section of tube walls can be computed as(27)Asn=4πYzXz.UsingAsn and An, we can further obtain total area of the cross section including the tube walls which is designated Atn:(28)Atn=An+Asn=4ab+3πYzXz.Similarly, the polar moment of inertia of arbitrary vertical unit tubeIpn can be computed as(29)Ipn=43a3b+b3a−π16−49πYzX3z−π4b−4Xz3π2YzXz−π16−49πY3zXz−π4a−4Yz3π2YzXz.Based on the same equivalence principle, the equivalent hydraulic conductivity of the whole soil layer in vertical direction is given by(30)kz∗=aγw8π2Az∗∫0LμiAtnIpn/An5dx.


## 2.3. Other Particle Combinations

According to the hypothesis in Section2.1, the soil layer is composed of numerous “particle systems.” These particle systems can be divided into three-particle combination, four-particle combination, and five-particle combination according to void ratio e. Similar to the proposed method for four-particle combination, the computation equation of hydraulic conductivity for soil under two other combinations can be obtained. However, the arrangement of particles for three or five particles in the vertical plane is not as regular as that of four particles. For three-particle system as well as five-particle system, in the YZ plane, there is a certain gap between the upper and lower rows of soil particles, which means the seepage route is not connected along the Z axis as shown in Figure 8. As a result, the calculation method proposed in this paper is not applicable for calculating the vertical hydraulic conductivity of these arrangements. However, for the sake of clarity, we did not deduce new formulas for the other two cases. Equation (30) is recommended to calculate the vertical hydraulic conductivity for soil with all particle systems. This means the effect of soil particle arrangement on hydraulic conductivity is ignored during calculation; only the size of soil particles is taken into consideration. As a result, the size of pore, especially for three-particle combination, could be overestimated, and the predicted results will be inaccurate. Therefore, the vertical calculation results for three-particle and five-particle combination would merely be the reference resources. More precise calculation methods still need to be further studied.Figure 8
Projection of three-particle systems inYZ plane.


### 2.3.1. Three-Particle Arrangement

The soil layer consisting of three-particle system usually characterizes of small pores and low hydraulic conductivity. As shown in Figure9(a), the shaded part represents the minimum seepage area and the dashed triangle indicates the maximum pore area. Under this condition, the cross-sectional area Ai and the polar moment of inertia can be expressed as(31)Aix=3a2−π2a2−axb2,Ipix=136b4πa3403b2πb−xb+xa21−x2b2+ab442+143π−33π2+42b2−2+π2x2+314−3π2x4.Figure 9
Other particle combinations. (a) The three-particle system. (b) The five-particle system.
(a)
(b)Furthermore,Ax∗,Ati can be obtained as(32)Ax∗=∫−bbAixdx2b=33−π3a2,Ati=Ai+Asi=3a2+2.5πa2−axb2.Finally, substituting formulas ofAi, Ipi, and Ax∗ into equation (24), the equivalent hydraulic conductivity in horizontal direction with three particles can be obtained.


### 2.3.2. Five-Particle Arrangement

The size of pore formed by five particles is relatively large. The maximum seepage cross-sectional area can be simplified as regular pentagon while the minimum area is surrounded by five tangent circles, as the shaded area shown in Figure9(b). Similarly, the cross-sectional area Ai and the polar moment of inertia Ipi are given by(33)Aix=6.89a2−1.5πa2−axb2,Ipix=a3α0a21−x2/b2b4−b2x2+aα1b4+α2b2x2+α3x4b4,where α0=−5.489,α1=16.082,α2=−2.736,α3=−2.726.Thus, equation (33) yields(34)Ax∗=∫−bbAixdx2b=5.32a2,(35)Ati=Ai+Asi=6.89a2+3.5πa2−axb2.Combining equations (33)–(36), the equivalent hydraulic conductivity of whole soil layer in horizontal direction can be obtained for five-particle system.


## 2.3.1. Three-Particle Arrangement

The soil layer consisting of three-particle system usually characterizes of small pores and low hydraulic conductivity. As shown in Figure9(a), the shaded part represents the minimum seepage area and the dashed triangle indicates the maximum pore area. Under this condition, the cross-sectional area Ai and the polar moment of inertia can be expressed as(31)Aix=3a2−π2a2−axb2,Ipix=136b4πa3403b2πb−xb+xa21−x2b2+ab442+143π−33π2+42b2−2+π2x2+314−3π2x4.Figure 9
Other particle combinations. (a) The three-particle system. (b) The five-particle system.
(a)
(b)Furthermore,Ax∗,Ati can be obtained as(32)Ax∗=∫−bbAixdx2b=33−π3a2,Ati=Ai+Asi=3a2+2.5πa2−axb2.Finally, substituting formulas ofAi, Ipi, and Ax∗ into equation (24), the equivalent hydraulic conductivity in horizontal direction with three particles can be obtained.


## 2.3.2. Five-Particle Arrangement

The size of pore formed by five particles is relatively large. The maximum seepage cross-sectional area can be simplified as regular pentagon while the minimum area is surrounded by five tangent circles, as the shaded area shown in Figure9(b). Similarly, the cross-sectional area Ai and the polar moment of inertia Ipi are given by(33)Aix=6.89a2−1.5πa2−axb2,Ipix=a3α0a21−x2/b2b4−b2x2+aα1b4+α2b2x2+α3x4b4,where α0=−5.489,α1=16.082,α2=−2.736,α3=−2.726.Thus, equation (33) yields(34)Ax∗=∫−bbAixdx2b=5.32a2,(35)Ati=Ai+Asi=6.89a2+3.5πa2−axb2.Combining equations (33)–(36), the equivalent hydraulic conductivity of whole soil layer in horizontal direction can be obtained for five-particle system.


## 2.4. Hydraulic Conductivity vs Void Ratioe

In addition to particle size, void ratio is another essential factor influencing the permeability of soil. In practice, void ratio will constantly change with the process of consolidation; as a result, the hydraulic conductivity will also vary with the change of void ratio. Based on the large number of experiments, several theoretical equations have been proposed to describe the relationship between void ratio and hydraulic conductivity. For example, Lambe and Whitman [60] argued that there is a logarithmic relationship between the hydraulic conductivity and the void ratio:(36)e=e0+Cklogkk0,where Ck denotes the coefficient of permeability change and e0 and k0 denote the initial void ratio and the initial hydraulic conductivity, respectively. This equation only considers the relationship between hydraulic conductivity and void ratio and cannot reflect the coupling effect of particle size and void ratio from the microscopic and macro perspective, respectively.This section takes the four-particle combination as an example and further derives the relation between hydraulic conductivity and pore ratio. As shown in Figure10, in the initial state, the upper and lower particles are not closely connected, and there is a certain size of gap designated as l0 which is determined by the void ratio e of soil. During the process of consolidation, the gap between particles is gradually compressed. As a result, hydraulic conductivity changed with the decline of l0 and e.Figure 10
The change of gap during the process of consolidation.Hence, the area of arbitrary cross section of unit tubesAi can be rewritten as(37)Ai=2a+l02−πa2−axb2.Furthermore, the cross-sectional areas of tube wallsAsi can be expressed as follows:(38)Asi=πa2−axb2.According to the assumption that the soil layer is composed of innumerable particle arrangement systems, void ratio in four-particle combination can represent the void ratio of the whole soil layer. Thus,(39)e=∫−bbAidx∫−bbAsidx=32a+l02−2πa22πa2.Then, the size of gapl0 can be expressed by void ratio e:(40)l0=2π1+e3−2a.Thus, the equivalent cross-sectional area of the microchannelAx∗ and the polar moment of inertia Ipi can be rewritten as(41)Ax∗=∫−bbAidx2b=4al+l2+236−πa2=4a2π1+e3−2a+2π1+e3−2a2+236−πa2,(42)Ipix=2a+l046−π2−329πa2−axb22−2πa2−axb22a+l2−4a2−ax/b23π2=2π1+e227−π2−329πa2−axb22−2πa2−axb22π1+e/32−4a2−ax/b23π2.Substituting equations (41) and (42) into equation (24), we can finally obtain the predictive equation of horizontal hydraulic conductivity considering the change of void ratio. According to the same principle, we can further derive the theoretical model for vertical hydraulic conductivity and other particle combinations.


## 3. Results and Analysis

### 3.1. Hydraulic Conductivity vs. Particle Sizea, b

The change of hydraulic conductivity of soil with the parameter of soil particle size is listed in Tables1 and 2. Due to the fact that the model presented in this paper is not suitable for calculating the vertical hydraulic conductivity of soil under the combination of three and five particles, this section does not analyze the change of vertical hydraulic conductivity of soil layer under the two different soil particle arrangements. From the theoretical point of view, 10 groups of theoretical particles within the range of sand soil are studied. The detailed information is shown in Tables 1 and 2.Table 1
Particle size for theoretical calculation.
Group
1
2
3
4
5
a (mm)
0.1
0.2
0.3
0.4
0.5
b (mm)
From 2.0 mm to 3.8 mm, interval 0.2 mm
The data used are divided into five groups, each of which only changes the size of the major half-axisb.Table 2
Particle size for theoretical calculation.
Group
6
7
8
9
10
b (mm)
2.0
2.5
3.0
3.5
4.0
a (mm)
From 0.1 mm to 0.28 mm, interval 0.02 mm
The data used are divided into five groups, each of which only changes the size of the major half-axisa.


#### 3.1.1. Four-Particle Arrangement

Figure11 shows the relationship between critical velocity and particle size for four-particle arrangement under laminar flow via equation (8). Taking b = 2.0 mm as an example, when the coordinates of particle size and critical velocity fall within the shaded area, the flow pattern can be considered as laminar flow.Figure 11
Relationship between particle size and critical velocity in laminar flow under four-particle arrangement.Moreover, the variations of hydraulic conductivity with the axial length of particles in two directions are shown in Figures12 and 13, respectively. Figure 12 shows that the hydraulic conductivity of soil gradually linearly decreases with the increase of particle size. For higher a, the decrease trend of hydraulic conductivity is more significant. Figure 13 shows that the relation between a and hydraulic conductivity can be expressed as a power function. Generally, minor axis size a has a greater influence on hydraulic conductivity. This is because the hydraulic conductivity predicted by this study is mainly related to the cross-sectional area of pore, which is mainly influenced by the size of minor axis.Figure 12
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the minor axis of the particles is constant, while the size of the major axis is increasing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Figure 13
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the major axis of the particles is constant, while the size of the minor axis is changing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Furthermore, comparing the variation of hydraulic conductivity in horizontal and vertical directions, it can be seen thatkz is generally about 1.54 times of kx under the same particle size. This behavior is due to fact that when fluid moves along the vertical direction, the cross-sectional area of pore is determined by the major and minor axis of particles. Accordingly, the cross-sectional area of vertical pore is larger than that of horizontal tube, resulting in a higher permeability of soil layer vertically. It should be pointed out that the overall trends for vertical and horizontal permeability are basically the same. From formulas (24) and (30), it can be seen that the hydraulic conductivity in this model is only related to the area and shape of the pore. In horizontal and vertical directions, when the size of soil particles changes, only the size of voids between soil particles changed. Hence, the changes of hydraulic conductivity are consistent without considering the slip between particles.


#### 3.1.2. Three-Particle Arrangement and Five-Particle Arrangement

Figure14 shows the relationship between critical velocity and particle size for laminar flow under three-particle arrangement and five-particle arrangement. For soils consisting of the three-particle system and five-particle system, the effect of particle size on permeability is analyzed only in horizontal direction. With the increase of particle size, the variation of hydraulic conductivity under these two types of particle arrangements is basically the same as that of four-particle arrangement. Generally, the microchannel enclosed by three particles has smaller cross-sectional area; accordingly, its hydraulic conductivity is low (Figures 15(a) and 15(b)). Under the same size, the hydraulic conductivity of the soil layer that consists of three-particle system is about 15% of that of the four-particle system. For the soil consisting of five-particle arrangement, the permeability is obviously higher because of its greater cross-sectional area. As a result, the hydraulic conductivity computed by five-particle model is 4.57 times that of four-particle model (Figures 15(c) and 15(d)). According to the above analysis, it can be concluded that different arrangements of soil particles will exert a significant impact on the permeability of the soil layer.Figure 14
Relationship between particle size and critical velocity in laminar flow under three-particle arrangement and five-particle arrangement. (a) Three-particle arrangement. (b) Five-particle arrangement.
(a)
(b)Figure 15
The changing law of hydraulic conductivity with the increasing sizes of particles under three-particle arrangement (a, b) and five-particle arrangement (c, d).
(a)
(b)
(c)
(d)


### 3.2. Effects of Void Ratioe

To analyze the relationship between void ratio and hydraulic conductivity based on current model, the sand soil with different particle arrangement conditions are analyzed and the results are shown in Figures16 and 17. For three kinds of particle combinations, the values of void ratio are smallest when the upper and lower particles closely contact (l0=0). In that case, void ratios for different combinations are(43)e3p=33π−1=0.654,e4p=6π−1=0.910,e5p=32sin54°+1tan54°+1/cos54°+25π−1=1.539.Figure 16
The relationship between hydraulic conductivity and void ratio under different particle combinations.Figure 17
Vertical hydraulic conductivity vs horizontal hydraulic conductivity under four-particle combination.Hence, we can obtain the curve of three-particle combination frome = 0.654, that of four-particle combination from e = 0.910, and that of five-particle combination from e = 1.539. It is clear that the changes of hydraulic conductivity are not same with the increase of void ratio under the three combination conditions as shown in Figure 16. Generally, the hydraulic conductivity gradually rises with the increase of void ratio. Comparing three kinds of combination conditions, it is obvious that the hydraulic conductivity increases more rapidly and the slope of curve is lager under the five-particle arrangement. Specifically, when e = 2.0, the size of kx is 0.417 cm/s. With the process of consolidation, the void ratio decreases gradually. When e = 1.5, kx is only 0.106 cm/s, which decreases by 74.6%. In contrast, the hydraulic conductivity changes slowly under the three-particle arrangement, and the overall growth rate is low.As can be seen from Figure17, with the increase of void ratio, the changes of hydraulic conductivity are both smooth curved, and the curvature radii of curves in such two directions are different. When void ratio falls into different ranges, the slopes of the curves are totally different. Because the derivative equation of hydraulic conductivity is very complex, the change rate of hydraulic conductivity will change gradually, rather than keep constant with the increase of void ratio. According to the slopes of the curves, we divide the curves into three parts: ① when 0.95 < e < 1.55, θv<θH (θ is the angularity of tangent line which can describe the slope of curve; subscript V and H mean vertical direction and horizontal direction, respectively). In that case, with the increase of void ratio, the growth rate of hydraulic conductivity in horizontal direction is greater than that in vertical direction; ② when 1.55 < e < 1.85, θv≈θH and the growth rates of hydraulic conductivity in different directions are basically the same; ③ when 1.85 < e < 2.5, θv>θH. The growth rate of hydraulic conductivity in vertical direction is greater than that in horizontal direction, and the gap between two directions rises gradually with the increase of void ratio.In addition, the vertical hydraulic conductivity is always greater than the horizontal hydraulic conductivity with the increase of void ratio. This is because, from the perspective of the geometry of the soil layer, when the coordinate values of the two directions are equal (X = Y), the pore cross-sectional area in the vertical direction is always larger than that in the horizontal direction, so the vertical permeability is better than that in the horizontal direction; from the perspective of the model equation, the values of hydraulic conductivity in equations (24) and (30) are mainly affected by the higher order terms Ai5 and An5. Due to that An5>Ai5, the vertical hydraulic conductivity is larger than the horizontal hydraulic conductivity. Furthermore, when l0 remains unchanged, the increment of vertical pore section area ΔAn is larger than that of horizontal direction ΔAi, which also affects the permeability of soil layer.


## 3.1. Hydraulic Conductivity vs. Particle Sizea, b

The change of hydraulic conductivity of soil with the parameter of soil particle size is listed in Tables1 and 2. Due to the fact that the model presented in this paper is not suitable for calculating the vertical hydraulic conductivity of soil under the combination of three and five particles, this section does not analyze the change of vertical hydraulic conductivity of soil layer under the two different soil particle arrangements. From the theoretical point of view, 10 groups of theoretical particles within the range of sand soil are studied. The detailed information is shown in Tables 1 and 2.Table 1
Particle size for theoretical calculation.
Group
1
2
3
4
5
a (mm)
0.1
0.2
0.3
0.4
0.5
b (mm)
From 2.0 mm to 3.8 mm, interval 0.2 mm
The data used are divided into five groups, each of which only changes the size of the major half-axisb.Table 2
Particle size for theoretical calculation.
Group
6
7
8
9
10
b (mm)
2.0
2.5
3.0
3.5
4.0
a (mm)
From 0.1 mm to 0.28 mm, interval 0.02 mm
The data used are divided into five groups, each of which only changes the size of the major half-axisa.


### 3.1.1. Four-Particle Arrangement

Figure11 shows the relationship between critical velocity and particle size for four-particle arrangement under laminar flow via equation (8). Taking b = 2.0 mm as an example, when the coordinates of particle size and critical velocity fall within the shaded area, the flow pattern can be considered as laminar flow.Figure 11
Relationship between particle size and critical velocity in laminar flow under four-particle arrangement.Moreover, the variations of hydraulic conductivity with the axial length of particles in two directions are shown in Figures12 and 13, respectively. Figure 12 shows that the hydraulic conductivity of soil gradually linearly decreases with the increase of particle size. For higher a, the decrease trend of hydraulic conductivity is more significant. Figure 13 shows that the relation between a and hydraulic conductivity can be expressed as a power function. Generally, minor axis size a has a greater influence on hydraulic conductivity. This is because the hydraulic conductivity predicted by this study is mainly related to the cross-sectional area of pore, which is mainly influenced by the size of minor axis.Figure 12
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the minor axis of the particles is constant, while the size of the major axis is increasing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Figure 13
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the major axis of the particles is constant, while the size of the minor axis is changing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Furthermore, comparing the variation of hydraulic conductivity in horizontal and vertical directions, it can be seen thatkz is generally about 1.54 times of kx under the same particle size. This behavior is due to fact that when fluid moves along the vertical direction, the cross-sectional area of pore is determined by the major and minor axis of particles. Accordingly, the cross-sectional area of vertical pore is larger than that of horizontal tube, resulting in a higher permeability of soil layer vertically. It should be pointed out that the overall trends for vertical and horizontal permeability are basically the same. From formulas (24) and (30), it can be seen that the hydraulic conductivity in this model is only related to the area and shape of the pore. In horizontal and vertical directions, when the size of soil particles changes, only the size of voids between soil particles changed. Hence, the changes of hydraulic conductivity are consistent without considering the slip between particles.


### 3.1.2. Three-Particle Arrangement and Five-Particle Arrangement

Figure14 shows the relationship between critical velocity and particle size for laminar flow under three-particle arrangement and five-particle arrangement. For soils consisting of the three-particle system and five-particle system, the effect of particle size on permeability is analyzed only in horizontal direction. With the increase of particle size, the variation of hydraulic conductivity under these two types of particle arrangements is basically the same as that of four-particle arrangement. Generally, the microchannel enclosed by three particles has smaller cross-sectional area; accordingly, its hydraulic conductivity is low (Figures 15(a) and 15(b)). Under the same size, the hydraulic conductivity of the soil layer that consists of three-particle system is about 15% of that of the four-particle system. For the soil consisting of five-particle arrangement, the permeability is obviously higher because of its greater cross-sectional area. As a result, the hydraulic conductivity computed by five-particle model is 4.57 times that of four-particle model (Figures 15(c) and 15(d)). According to the above analysis, it can be concluded that different arrangements of soil particles will exert a significant impact on the permeability of the soil layer.Figure 14
Relationship between particle size and critical velocity in laminar flow under three-particle arrangement and five-particle arrangement. (a) Three-particle arrangement. (b) Five-particle arrangement.
(a)
(b)Figure 15
The changing law of hydraulic conductivity with the increasing sizes of particles under three-particle arrangement (a, b) and five-particle arrangement (c, d).
(a)
(b)
(c)
(d)


## 3.1.1. Four-Particle Arrangement

Figure11 shows the relationship between critical velocity and particle size for four-particle arrangement under laminar flow via equation (8). Taking b = 2.0 mm as an example, when the coordinates of particle size and critical velocity fall within the shaded area, the flow pattern can be considered as laminar flow.Figure 11
Relationship between particle size and critical velocity in laminar flow under four-particle arrangement.Moreover, the variations of hydraulic conductivity with the axial length of particles in two directions are shown in Figures12 and 13, respectively. Figure 12 shows that the hydraulic conductivity of soil gradually linearly decreases with the increase of particle size. For higher a, the decrease trend of hydraulic conductivity is more significant. Figure 13 shows that the relation between a and hydraulic conductivity can be expressed as a power function. Generally, minor axis size a has a greater influence on hydraulic conductivity. This is because the hydraulic conductivity predicted by this study is mainly related to the cross-sectional area of pore, which is mainly influenced by the size of minor axis.Figure 12
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the minor axis of the particles is constant, while the size of the major axis is increasing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Figure 13
Variation of hydraulic conductivity with particle size under four-particle arrangement. In the same curve, the size of the major axis of the particles is constant, while the size of the minor axis is changing. (a) Horizontal hydraulic conductivity. (b) Vertical hydraulic conductivity.
(a)
(b)Furthermore, comparing the variation of hydraulic conductivity in horizontal and vertical directions, it can be seen thatkz is generally about 1.54 times of kx under the same particle size. This behavior is due to fact that when fluid moves along the vertical direction, the cross-sectional area of pore is determined by the major and minor axis of particles. Accordingly, the cross-sectional area of vertical pore is larger than that of horizontal tube, resulting in a higher permeability of soil layer vertically. It should be pointed out that the overall trends for vertical and horizontal permeability are basically the same. From formulas (24) and (30), it can be seen that the hydraulic conductivity in this model is only related to the area and shape of the pore. In horizontal and vertical directions, when the size of soil particles changes, only the size of voids between soil particles changed. Hence, the changes of hydraulic conductivity are consistent without considering the slip between particles.


## 3.1.2. Three-Particle Arrangement and Five-Particle Arrangement

Figure14 shows the relationship between critical velocity and particle size for laminar flow under three-particle arrangement and five-particle arrangement. For soils consisting of the three-particle system and five-particle system, the effect of particle size on permeability is analyzed only in horizontal direction. With the increase of particle size, the variation of hydraulic conductivity under these two types of particle arrangements is basically the same as that of four-particle arrangement. Generally, the microchannel enclosed by three particles has smaller cross-sectional area; accordingly, its hydraulic conductivity is low (Figures 15(a) and 15(b)). Under the same size, the hydraulic conductivity of the soil layer that consists of three-particle system is about 15% of that of the four-particle system. For the soil consisting of five-particle arrangement, the permeability is obviously higher because of its greater cross-sectional area. As a result, the hydraulic conductivity computed by five-particle model is 4.57 times that of four-particle model (Figures 15(c) and 15(d)). According to the above analysis, it can be concluded that different arrangements of soil particles will exert a significant impact on the permeability of the soil layer.Figure 14
Relationship between particle size and critical velocity in laminar flow under three-particle arrangement and five-particle arrangement. (a) Three-particle arrangement. (b) Five-particle arrangement.
(a)
(b)Figure 15
The changing law of hydraulic conductivity with the increasing sizes of particles under three-particle arrangement (a, b) and five-particle arrangement (c, d).
(a)
(b)
(c)
(d)


## 3.2. Effects of Void Ratioe

To analyze the relationship between void ratio and hydraulic conductivity based on current model, the sand soil with different particle arrangement conditions are analyzed and the results are shown in Figures16 and 17. For three kinds of particle combinations, the values of void ratio are smallest when the upper and lower particles closely contact (l0=0). In that case, void ratios for different combinations are(43)e3p=33π−1=0.654,e4p=6π−1=0.910,e5p=32sin54°+1tan54°+1/cos54°+25π−1=1.539.Figure 16
The relationship between hydraulic conductivity and void ratio under different particle combinations.Figure 17
Vertical hydraulic conductivity vs horizontal hydraulic conductivity under four-particle combination.Hence, we can obtain the curve of three-particle combination frome = 0.654, that of four-particle combination from e = 0.910, and that of five-particle combination from e = 1.539. It is clear that the changes of hydraulic conductivity are not same with the increase of void ratio under the three combination conditions as shown in Figure 16. Generally, the hydraulic conductivity gradually rises with the increase of void ratio. Comparing three kinds of combination conditions, it is obvious that the hydraulic conductivity increases more rapidly and the slope of curve is lager under the five-particle arrangement. Specifically, when e = 2.0, the size of kx is 0.417 cm/s. With the process of consolidation, the void ratio decreases gradually. When e = 1.5, kx is only 0.106 cm/s, which decreases by 74.6%. In contrast, the hydraulic conductivity changes slowly under the three-particle arrangement, and the overall growth rate is low.As can be seen from Figure17, with the increase of void ratio, the changes of hydraulic conductivity are both smooth curved, and the curvature radii of curves in such two directions are different. When void ratio falls into different ranges, the slopes of the curves are totally different. Because the derivative equation of hydraulic conductivity is very complex, the change rate of hydraulic conductivity will change gradually, rather than keep constant with the increase of void ratio. According to the slopes of the curves, we divide the curves into three parts: ① when 0.95 < e < 1.55, θv<θH (θ is the angularity of tangent line which can describe the slope of curve; subscript V and H mean vertical direction and horizontal direction, respectively). In that case, with the increase of void ratio, the growth rate of hydraulic conductivity in horizontal direction is greater than that in vertical direction; ② when 1.55 < e < 1.85, θv≈θH and the growth rates of hydraulic conductivity in different directions are basically the same; ③ when 1.85 < e < 2.5, θv>θH. The growth rate of hydraulic conductivity in vertical direction is greater than that in horizontal direction, and the gap between two directions rises gradually with the increase of void ratio.In addition, the vertical hydraulic conductivity is always greater than the horizontal hydraulic conductivity with the increase of void ratio. This is because, from the perspective of the geometry of the soil layer, when the coordinate values of the two directions are equal (X = Y), the pore cross-sectional area in the vertical direction is always larger than that in the horizontal direction, so the vertical permeability is better than that in the horizontal direction; from the perspective of the model equation, the values of hydraulic conductivity in equations (24) and (30) are mainly affected by the higher order terms Ai5 and An5. Due to that An5>Ai5, the vertical hydraulic conductivity is larger than the horizontal hydraulic conductivity. Furthermore, when l0 remains unchanged, the increment of vertical pore section area ΔAn is larger than that of horizontal direction ΔAi, which also affects the permeability of soil layer.


## 4. Verification and Applications

In previous studies, the shape of soil particles is deemed as spheres rather than ellipsoids. Hence, the relative experimental data concerning the influences ofa and b on permeability are insufficient. Some scholars [61, 62] adopted average diameter of pore dpmm and specific surface Ssm2/g to indicate the influence of soil particle size. There is a strong correlation between dp, Ss, and soil particle size. Thereby, the theoretical model proposed by this study can be, to some extent, validated by dp, Ss. First, based on the volume equivalent principle, we assume the variable cross-sectional pore as the cylinder with constant cross section. Then, the relationship between mean cross-sectional area A∗ (A∗ can be expressed by particle size a, b) and mean pore diameter dp can be expressed as follows:(44)A∗a,b=πdp24,and it can be rewritten as follows:(45)dpa,b=4A∗a,bπ.However, the specific surface areaSs is not reported in most studies. We estimate the specific surface area by using other published soil descriptions. In the case of fine-particle soils, estimates were based on liquid limits wl:(46)Ss=1.8wl−34.For sandy soils made of rotund grains, the specific surface area can be estimated from the cumulative grain size distribution:(47)Ss=3Cu+74ρwGsD50,where the coefficient of uniformity is Cu=D60/D10 and D10,D50,andD60 (mm) are the grain diameters for 10%, 50%, and 60% cumulative passing fractions; Gs indicates the specific weight of soil; and ρw is the density of soil. By equation (46) or equation (47), we can calculate the values of Ss by a series of geotechnical parameters. Furthermore, the method for computing particle sizes a and b proposed by Santamarina [63] is adopted in this study, expressed as(48)Ssa,b=31−1/p2apbp+a2p1/pa2bρm1+e,where p=1.6075.Then, by referring to the experimental data in other literatures, we can obtain the values ofdp, Ss and other parameters of different soils. Combining equations (45) and (48), dp, Ss will be transformed to the sizes of particles (a, b). By comparing the theoretical calculation results with the actual data [64–66] and other existing calculation models, such as Indraratna model, Hazen model, Kozeny-Carman model, we can verify the rationality of the model. The calculation models of the three-particle system presented in this paper have good reliability in calculating the hydraulic conductivity with different sizes of soil particles (see Tables 3 and 4). It can be concluded from Table 3 that for most of the size of soil particles, the four-particle model is more suitable for calculation; for example, when particle size a falls within 0∼0.5 mm, the error between predicted and measured hydraulic conductivity calculated by the four-particle model is less than 50%. However, for soils with larger particle size (1 mm < a < 5 mm), the calculated results using the three-particle model are closer to the measured values. Generally, although there is a certain disparity between the predicted results and the measured values, the proposed model still has high reliability.Table 3
Computed sizes of particles and corresponding values of saturated hydraulic conductivities.
Test id
d
p (mm)
a (mm)
b (mm)
Saturated hydraulic conductivitiesk(cm/s)
Test data
Current model
Indraratna model
Hazen model
Kozeny-Carman model
Horizontal
Vertical
k
3
x
k
4
x
k
5
x
k
z
Indraratna et al. [67]
1
0.0310
0.0199
0.1080
0.0119
0.0012
0.0071
0.0303
0.0108
0.0122
0.0124
0.0192
2
0.0640
0.0411
0.2200
0.0303
0.0051
0.0304
0.1289
0.0460
0.04
0.0503
0.08
3
0.1860
0.1194
0.6200
0.31
0.0435
0.2581
1.0891
0.3895
0.232
0.4229
0.6761
4
0.0280
0.0180
0.0970
0.008
0.0010
0.0058
0.0247
0.0088
0.0104
0.01
0.0152
5
0.0340
0.0218
0.1190
0.0087
0.0015
0.0086
0.0364
0.0129
0.0144
0.0089
0.0182
6
0.1620
0.1040
0.5600
0.2
0.0330
0.1946
0.8262
0.2942
0.1856
0.2015
0.4085
7
0.0330
0.0212
0.1150
0.0098
0.0014
0.0081
0.0343
0.0122
0.0138
0.01
0.0187
8
0.0360
0.0231
0.1290
0.0103
0.0016
0.0096
0.0408
0.0145
0.0154
0.01
0.02
9
0.0380
0.0244
0.1300
0.011
0.0018
0.0107
0.0455
0.0162
0.0169
0.01
0.0211
10
0.0400
0.0257
0.1400
0.012
0.0020
0.0118
0.0504
0.0179
0.0182
0.01
0.0219
11
0.0410
0.0263
0.1400
0.022
0.0021
0.0125
0.0529
0.0189
0.0193
0.01
0.0223
12
0.0430
0.0276
0.1500
0.022
0.0023
0.0137
0.0582
0.0207
0.0211
0.0113
0.0251
13
0.0410
0.0263
0.1400
0.0092
0.0021
0.0125
0.0529
0.0189
0.0197
0.0109
0.024
14
0.0880
0.0565
0.3100
0.0909
0.0097
0.0573
0.2438
0.0866
0.0682
0.0488
0.1075
15
0.1110
0.0713
0.3800
0.1214
0.0155
0.0915
0.3879
0.1383
0.0989
0.0814
0.1754
16
0.1310
0.0841
0.4500
0.143
0.0216
0.1274
0.5402
0.1925
0.1312
0.1214
0.2557
Sherard et al. [68]
1
0.0640
0.0411
0.1643
0.025
0.0051
0.0316
0.1289
0.0470
0.0404
0.023
0.0683
2
0.0420
0.0270
0.1450
0.018
0.0022
0.0131
0.0555
0.0198
0.0202
0.023
0.0325
3
0.1310
0.0841
0.4500
0.049
0.0216
0.1274
0.5402
0.1925
0.1308
0.0576
0.2383
4
0.0830
0.0533
0.2880
0.022
0.0087
0.0510
0.2169
0.0772
0.0614
0.0576
0.094
5
0.1550
0.0995
0.5400
0.051
0.0302
0.1779
0.7563
0.2690
0.1722
0.0756
0.3342
6
0.0960
0.0616
0.3300
0.039
0.0116
0.0684
0.2901
0.1034
0.0789
0.0756
0.1173
7
0.3410
0.2189
1.2200
0.56
0.1461
0.8572
3.6606
1.2977
0.6291
0.7921
2.0125
8
0.2250
0.1444
0.5778
0.25
0.0636
0.3911
1.5937
0.5814
0.3168
0.7921
0.9806
9
0.7650
0.4911
2.9000
3.8
0.7353
4.2704
18.4233
6.4790
2.3766
4.1229
10.619
10
0.5060
0.3248
1.6600
2.6
0.3217
1.9147
8.0602
2.8880
1.2049
4.1229
5.2769
11
1.2400
0.7961
3.8220
5.3
1.9320
11.6102
48.4048
17.458
5.2555
5.0176
22.769
12
0.7840
0.5033
2.4600
4.8
0.7723
4.6287
19.3498
6.9663
2.4761
5.0176
9.5476
13
0.8440
0.5418
3.4000
5.6
0.8951
5.1409
22.4248
7.8156
2.7934
3.7249
11.189
14
0.5040
0.3236
1.8943
1.8
0.3192
1.8564
7.9966
2.8156
1.1949
3.7249
3.5844
15
2.2700
1.4573
7.9993
12.2
6.4747
38.0824
162.2170
57.620
14.2012
28.09
80.832
16
1.4640
0.9399
5.8595
8.7
2.6931
15.4866
67.4724
23.539
6.9094
28.09
36.87
17
4.0570
2.6046
9.9214
15.2
20.6813
127.675
518.1477
189.04
36.8931
92.16
273.06
18
2.6590
1.7071
6.8282
13.9
8.8839
54.625
222.5776
18.424
18.4239
92.16
130.09
Tsai [39]
1
0.0400
0.0257
0.138
0.025
0.0020
0.0119
0.0504
0.0179
0.0186
0.0272
0.0247
2
0.0400
0.0257
0.138
0.022
0.0020
0.0119
0.0504
0.0179
0.0184
0.0272
0.0244
3
0.0390
0.0250
0.136
0.027
0.0019
0.0113
0.0479
0.0170
0.0182
0.0272
0.0241
4
0.0390
0.0250
0.136
0.024
0.0019
0.0113
0.0479
0.0170
0.018
0.0272
0.024
5
0.0390
0.0250
0.136
0.025
0.0019
0.0113
0.0479
0.0170
0.0178
0.0272
0.0237
6
0.0390
0.0250
0.136
0.027
0.0019
0.0113
0.0479
0.0170
0.0177
0.0272
0.0236
7
0.0400
0.0257
0.138
0.034
0.0020
0.0119
0.0504
0.0168
0.0182
0.0272
0.0244
8
0.0940
0.0603
0.328
0.168
0.0111
0.0654
0.2782
0.0989
0.0752
0.2418
0.1951
9
0.0490
0.0315
0.173
0.0311
0.0030
0.0177
0.0756
0.0268
0.0259
0.0613
0.0454
10
0.0440
0.0282
0.154
0.0254
0.0024
0.0156
0.0609
0.0217
0.0219
0.0543
0.0375
11
0.0400
0.0257
0.138
0.0257
0.0020
0.0119
0.0504
0.0179
0.0184
0.0441
0.0299
12
0.0310
0.0199
0.109
0.0129
0.0012
0.0071
0.0303
0.0107
0.0121
0.0256
0.0177
13
0.0320
0.0205
0.14
0.0124
0.0013
0.0073
0.0322
0.0111
0.013
0.0262
0.0186
14
0.0310
0.0199
0.109
0.0133
0.0012
0.0071
0.0303
0.0107
0.0124
0.0262
0.0181
15
0.0320
0.0205
0.1116
0.0135
0.0013
0.0076
0.0322
0.0115
0.0129
0.0262
0.0192
16
0.0320
0.0205
0.1116
0.0126
0.0013
0.0076
0.0322
0.0115
0.0131
0.0262
0.0196
17
0.0430
0.0276
0.151
0.0218
0.0023
0.0137
0.0582
0.0207
0.021
0.0441
0.0342
18
0.0410
0.0263
0.142
0.0171
0.0021
0.0125
0.0529
0.0188
0.0197
0.0441
0.0318
19
0.0390
0.0250
0.136
0.0233
0.0019
0.0113
0.0479
0.0170
0.0177
0.0441
0.0283
20
0.0420
0.0270
0.148
0.024
0.0022
0.0130
0.0555
0.0197
0.0203
0.0441
0.0322
21
0.0410
0.0263
0.142
0.0257
0.0021
0.0125
0.0529
0.0188
0.019
0.0441
0.0313
22
0.0380
0.0244
0.133
0.0247
0.0018
0.0107
0.0455
0.0162
0.0168
0.0441
0.0284
23
0.0370
0.0238
0.131
0.0241
0.0017
0.0101
0.0431
0.0153
0.0164
0.0441
0.028
24
0.0360
0.0231
0.126
0.0279
0.0016
0.0096
0.0408
0.0145
0.0158
0.0441
0.0271Table 4
Computational error of theoretical models vs test data.
Test id
Hydraulic conductivity test data (cm/s)
Computational error of theoretical models (%)
Current model
Other models
k
3
x
k
4
x
k
5
x
Indraratna model
Hazen model
Kozeny-Carman model
Indraratna et al. [67]
1
0.0119
−89.92
−40.34
154.62
2.52
4.20
61.34
2
0.0303
−83.17
0.33
325.41
32.01
66.01
164.03
3
0.31
−85.97
−16.74
251.32
−25.16
36.42
118.10
4
0.008
−87.50
−27.50
208.75
30.00
25.00
90.00
5
0.0087
−82.76
−1.15
318.39
65.52
2.30
109.20
6
0.2
−83.50
−2.70
313.10
−7.20
0.75
104.25
7
0.0098
−85.71
−17.35
250.00
40.82
2.04
90.82
8
0.0103
−84.47
−6.80
296.12
49.51
−2.91
94.17
9
0.011
−83.64
−2.73
313.64
53.64
−9.09
91.82
10
0.012
−83.33
−1.67
320.00
51.67
−16.67
82.50
11
0.022
−90.45
−43.18
140.45
−12.27
−54.55
1.36
12
0.022
−89.55
−37.73
164.55
−4.09
−48.64
14.09
13
0.0092
−77.17
35.87
475.00
114.13
18.48
160.87
14
0.0909
−89.33
−36.96
168.21
−24.97
−46.31
18.26
15
0.1214
−87.23
−24.63
219.52
−18.53
−32.95
44.48
16
0.143
−84.90
−10.91
277.76
−8.25
−15.10
78.81
Sherard et al. [68]
1
0.025
−79.60
26.40
415.60
61.60
−8.00
173.20
2
0.018
−87.78
−27.22
208.33
12.22
27.78
80.56
3
0.049
−55.92
160.00
1002.45
166.94
17.55
386.33
4
0.022
−60.45
131.82
885.91
179.09
161.82
327.27
5
0.051
−40.78
248.82
1382.94
237.65
48.24
555.29
6
0.039
−70.26
75.38
643.85
102.31
93.85
200.77
7
0.56
−73.91
53.07
553.68
12.34
41.45
259.38
8
0.25
−74.56
56.44
537.48
26.72
216.84
292.24
9
3.8
−80.65
12.38
384.82
−37.46
8.50
179.45
10
2.6
−87.63
−26.36
210.01
−53.66
58.57
102.96
11
5.3
−63.55
119.06
813.30
−0.84
−5.33
329.60
12
4.8
−83.91
−3.57
303.12
−48.41
4.53
98.91
13
5.6
−84.02
−8.20
300.44
−50.12
−33.48
99.80
14
1.8
−82.27
3.13
344.26
−33.62
106.94
99.13
15
12.2
−46.93
212.15
1229.65
16.40
130.25
562.56
16
8.7
−69.04
78.01
675.54
−20.58
222.87
323.79
17
15.2
36.06
739.97
3308.87
142.72
506.32
1696.45
18
13.9
−36.09
292.99
1501.28
32.55
563.02
835.90
Tsai [39]
1
0.025
−92.00
−52.40
101.60
−25.60
8.80
−1.20
2
0.022
−90.91
−45.91
129.09
−16.36
23.64
10.91
3
0.027
−92.96
−58.15
77.41
−32.59
0.74
−10.74
4
0.024
−92.08
−52.92
99.58
−25.00
13.33
0.00
5
0.025
−92.40
−54.80
91.60
−28.80
8.80
−5.20
6
0.027
−92.96
−58.15
77.41
−34.44
0.74
−12.59
7
0.034
−94.12
−65.00
48.24
−46.47
−20.00
−28.24
8
0.168
−93.39
−61.07
65.60
−55.24
43.93
16.13
9
0.0311
−90.35
−43.09
143.09
−16.72
97.11
45.98
10
0.0254
−90.55
−38.58
139.76
−13.78
113.78
47.64
11
0.0257
−92.22
−53.70
96.11
−28.40
71.60
16.34
12
0.0129
−90.70
−44.96
134.88
−6.20
98.45
37.21
13
0.0124
−89.52
−41.13
159.68
4.84
111.29
50.00
14
0.0133
−90.98
−46.62
127.82
−6.77
96.99
36.09
15
0.0135
−90.37
−43.70
138.52
−4.44
94.07
42.22
16
0.0126
−89.68
−39.68
155.56
3.97
107.94
55.56
17
0.0218
−89.45
−37.16
166.97
−3.67
102.29
56.88
18
0.0171
−87.72
−26.90
209.36
15.20
157.89
85.96
19
0.0233
−91.85
−51.50
105.58
−24.03
89.27
21.46
20
0.024
−90.83
−45.83
131.25
−15.42
83.75
34.17
21
0.0257
−91.83
−51.36
105.84
−26.07
71.60
21.79
22
0.0247
−92.71
−56.68
84.21
−31.98
78.54
14.98
23
0.0241
−92.95
−58.09
78.84
−31.95
82.99
16.18
24
0.0279
−94.27
−65.59
46.24
−43.37
58.06
−2.87In order to validate the rationality of the model, we use the three-particle model and four-particle model to describe the relationship between horizontal hydraulic conductivity and void ratio for sandy soil and smectite soil, respectively, as shown in Figure18. It is clear that the predicted results have high accuracy and coincide with the test data. From Figure 18(a), we can conclude that the three-particle model has high accuracy for predicting hydraulic conductivity of sand soil. When the values of void ratio are high (e > 0.85), the predicted results of K-C model obviously deviate from the test data, while the predicted values still keep good agreement with the test data. This result also indirectly validates the previous conclusion: three-particle model is more suitable for soils with large size particles. Figure 18(b) shows that for smectite soil, predicted results by four-particle model are in good agreement with the measured values when the void ratio is small (e = 1∼4). In contrast, when the void ratio falls within the range of 4 to 10, the calculated results deviated from the measured value. Especially, when the void ratio is large enough (e > 6), the calculated results are obviously lower than the measured values and the accuracy is not as good as equations (1) and (37). Therefore, the model proposed by this paper is not suitable for the soil with high void ratio.Figure 18
Comparison of different (k ∼ e) models. The experimental data of (a) sand soil were measured by Taylor and Donald [42] and the (b) smectite by Mesri [45].
(a)
(b)To further examine the suitability of the proposed equation, referring to the data of geotechnical parameters (such as void ratio, specific surface area, and equivalent particle size) in many literatures, this study calculates abundant predicted values of hydraulic conductivity and compares them with the measured values in corresponding literatures, as shown in Figure19 [9, 61, 69–77].Figure 19
The measured hydraulic conductivityk versus the predicted k by equation (23) for three different types of soil: (a) three-particle system; (b) four-particle system; (c) five-particle system. Data source: [9, 61, 69–77].
(a)
(b)
(c)On the basis of a large number of calculations, when the measured results of a certain kind of soil are close to the values predicted by one of the three models, we think that this kind of soil is mainly composed of the corresponding particle combination of the model and assume that this model is suitable for calculating the hydraulic conductivity of this kind of soil. Based on the comparison between computed result and experimental data (as shown in Figure19), it can be conducted that the three prediction models (three-particle model, four-particle model, and five-particle model) proposed in this study are suitable for sandy soils, clay-silt mixed soils, and bentonite-silt mixed soils, respectively. There are three findings that can be drawn from Figure 19:(1)
For the sandy soil, the theoretical values computing by three-particle arrangement model is closer to the actual measured value. Most predicted values fall within0.3kmeas<kpred<3kmeas. However, for clayey-silty mixed soils and bentonite-silt mixed soils, the predicted values are obviously lower than the measured values. This observation also verifies the above conclusion that three-particle model is more suitable to predict for large granular soil. In addition, it can be speculated that for sandy soil, the arrangement of three particles may be dominant in quantity.(2)
The theoretical model of four-particle arrangement can calculate the hydraulic conductivity for the soil layer with clayey and silty mixtures more accurately. Despite certain predicted results with a large error, most calculated values still distribute in reasonable interval.(3)
Comparing with other two kinds of soil, the five-particle model can reliably predict the hydraulic conductivity of bentonite-silty mixed soil that when the size of particles is smallest, the results are of high accuracy. However,kpred is far greater than kmeas when calculating sandy soil with larger particle size.


## 5. Discussion and Limitations

Although this paper verifies the reliability of the model based on the substantial experimental data, there are still many shortcomings and limitations due to certain idealizations and assumptions. This section will discuss the limitations of the model and propose corresponding remedies. More appropriate optimization needs to be further studied.


### 5.1. Inhomogeneous Particle Geometry

Practically, the soil layer is composed of soil particles with different sizes, and the shape of soil particles is not limited to ellipsoid. The assumption in this paper is only applicable to the ideal soil layer. However, there are still some shortcomings when it is used in practice. Therefore, it is suggested that the parameters such as uniformity coefficientCu should be used to link the actual particle size with the even particle size in the computation model so that it can be applied to the actual condition circumstance.


### 5.2. Not considering Unavailable Pore

In this paper, it is assumed that the pore channels are all effective pore which means free water can flow in any pore and there is no blockage. However, because of the existence of “relatively small pores” and “dead-end pores” as well as the effect of clay surfaces and interlayers, not all of water is stored in fine-grained materials that participate in flow. Carman [65] defined water that does not circulate as immobile water and proposed the concept of effective void ratio ee. In the saturated soil, effective void ratio also represents the ratio of mobile water volume to the volume of solid [22, 78–80]. This paper ignores the effect of effective void ratio ee on hydraulic conductivity and assumes no existences of invalid pore which is different from the actual situation. In the next study, we can refer to the research methods of effective pores proposed by Ren et al. [36], assuming that there is a bound water film around the soil particles. This bound water film can affect the cross-sectional shape and size of the pore. We can establish the relationship between effective void ratio and hydraulic conductivity by observing this phenomenon.


### 5.3. Uncompleted Equation for Vertical Hydraulic Conductivity

In this paper, the formula for calculating the vertical hydraulic conductivity is proposed only for the case of four-particle arrangement, which has certain reliability. However, for the arrangement of three and five particles, because of the special arrangement form, there are dislocations among the layers of soil particles and the vertical pores are not completely connected. As a result, the proposed model cannot be applied to these two types of soil layers. In the future studies, we can introduce the concept of tortuosity. By establishing the relationship between hydraulic conductivity and tortuosity under the arrangement of three particles and five particles, we can discover the geometric characteristics of actual pore in these two cases and then obtain a revised equation for vertical hydraulic conductivity of soil layer.


## 5.1. Inhomogeneous Particle Geometry

Practically, the soil layer is composed of soil particles with different sizes, and the shape of soil particles is not limited to ellipsoid. The assumption in this paper is only applicable to the ideal soil layer. However, there are still some shortcomings when it is used in practice. Therefore, it is suggested that the parameters such as uniformity coefficientCu should be used to link the actual particle size with the even particle size in the computation model so that it can be applied to the actual condition circumstance.


## 5.2. Not considering Unavailable Pore

In this paper, it is assumed that the pore channels are all effective pore which means free water can flow in any pore and there is no blockage. However, because of the existence of “relatively small pores” and “dead-end pores” as well as the effect of clay surfaces and interlayers, not all of water is stored in fine-grained materials that participate in flow. Carman [65] defined water that does not circulate as immobile water and proposed the concept of effective void ratio ee. In the saturated soil, effective void ratio also represents the ratio of mobile water volume to the volume of solid [22, 78–80]. This paper ignores the effect of effective void ratio ee on hydraulic conductivity and assumes no existences of invalid pore which is different from the actual situation. In the next study, we can refer to the research methods of effective pores proposed by Ren et al. [36], assuming that there is a bound water film around the soil particles. This bound water film can affect the cross-sectional shape and size of the pore. We can establish the relationship between effective void ratio and hydraulic conductivity by observing this phenomenon.


## 5.3. Uncompleted Equation for Vertical Hydraulic Conductivity

In this paper, the formula for calculating the vertical hydraulic conductivity is proposed only for the case of four-particle arrangement, which has certain reliability. However, for the arrangement of three and five particles, because of the special arrangement form, there are dislocations among the layers of soil particles and the vertical pores are not completely connected. As a result, the proposed model cannot be applied to these two types of soil layers. In the future studies, we can introduce the concept of tortuosity. By establishing the relationship between hydraulic conductivity and tortuosity under the arrangement of three particles and five particles, we can discover the geometric characteristics of actual pore in these two cases and then obtain a revised equation for vertical hydraulic conductivity of soil layer.


## 6. Conclusion

(i)
Based on Akbari’s formula for flow resistance, the paper proposes a new method for calculating soil hydraulic conductivity. Compared with the traditional calculation method, the proposed model considers the influence of ellipsoidal soil particles on the hydraulic conductivity, which is closer to the actual soil layer. At the same time, the proposed method can more accurately calculate the hydraulic conductivity in different directions, reflecting the anisotropy of the soil layer.
(ii)
According to the range of void ratioe (0.6–1.5), this paper assumes that the soil layer consists of three different particle systems (three-particle system, four-particle system, and five-particle system), in which the pore cross-sectional area of five particles is the largest. As a result, the corresponding hydraulic conductivity is the highest.
(iii)
The proposed model can reflect the influence of particle size in different directions on permeability of soil layer. When the size of minor axis of soil particles increases, the hydraulic conductivity of soil layer increases significantly. In addition, the hydraulic conductivity linearly decreases with the increase of major axis size. Compared with the major axis, the particle minor axis size has greater influence on the permeability of the soil layer.
(iv)
The models can be used to calculate the varied hydraulic conductivity with different void ratios. The hydraulic conductivity increases more rapidly, and the slope of curve is larger for the five-particle arrangement. In contrast, the hydraulic conductivity changes slowly under the three-particle arrangement, and the overall growth rate is low. For the four-particle combination, the vertical hydraulic conductivity is always greater than the horizontal hydraulic conductivity with the increase of void ratio and the variation trends of hydraulic conductivity in the two directions are basically the same.
(v)
The three models deduced in this paper are of high applicability when applied to different types of soil layers. The five-particle model is more accurate to calculate hydraulic conductivity for bentonite-silt mixed soil with relatively small particle size and most predicted values fall within0.3kmeas<kpred<3kmeas. Four-particle model is more suitable for calculating silt-clay mixed soil layer with middle-size particle. When predicting the soil layer consisting of small particle, kpred is often less than kmeas. The three-particle model has high reliability in calculating the hydraulic conductivity of large-sized particles and is suitable for predicting for sandy soil and some large-sized silt.



---

*Source: 7213094-2019-07-31.xml*