Next Article in Journal
Dissolved Phosphorus Concentrations in Surface Runoff from Agricultural Land Based on Calcium–Acetate–Lactate Soluble Phosphorus Soil Contents
Previous Article in Journal
Efficient Use of Water in Tailings Management: New Technologies and Environmental Strategies for the Future of Mining
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of One-Dimensional Consolidation Considering Non-Darcian Flow Described by Non-Newtonian Index Incorporating Impeded Drainage Boundaries

1
School of Civil Engineering, Hunan University, Changsha 410012, China
2
College of Civil Engineering, Hunan City University, Yiyang 413000, China
3
Hunan Engineering Research Center of Structural Safety and Disaster Prevention for Urban Underground Infrastructure, Yiyang 413000, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(11), 1740; https://doi.org/10.3390/w14111740
Submission received: 13 April 2022 / Revised: 23 May 2022 / Accepted: 25 May 2022 / Published: 28 May 2022
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
The nonlinear flow law and soil boundaries greatly affect the dissipation process of soil consolidation. Thus, to study the impact of nonlinear flow under impeded drainage boundaries, the classical non-Darcian flow model described by non-Newtonian index was introduced. The numerical solutions are derived in detail by the finite difference method (FDM) for one-dimensional (1-D) consolidation incorporating the impeded boundaries, and the computer program is compiled. Then, comparing two analytical solutions based on Darcy’s law and a numerical case of Forchheimeer’s flow, the validity of the present method was verified. The numerical results indicate that there is a critical depth phenomenon for the non-Darcian flow incorporating impeded drainage boundaries. The excess pore water pressure of the soil below the critical depth dissipates more slowly than that of Darcy’s law, whereas the pore pressure of the soil above the critical depth dissipates more quickly than that of Darcy’s law. Moreover, considering that the non-Darcian flow with the non-Newtonian index will still delay the overall consolidation rate of the soft ground, the greater the nondimensional parameter I0 is, the more obvious the lagging phenomenon of the overall dissipation of pore pressure is.

1. Introduction

As one of the common issues in geotechnical engineering, soil consolidation has captured a lot of attention from geotechnical engineers since the 1920s [1]. To obtain the analytical solutions or numerical solutions more conveniently, the boundary of the soft soil layer is usually simplified as completely permeable (the first kind of drainage boundary) or absolutely impermeable (the second kind of drainage boundary) [2,3,4,5]. In fact, the drainage at the boundaries of the soil layer is partially permeable, i.e., is impeded [6]. For example, the sand cushion covering the soil layer and the underlying layer at the bottom of the soft soil will become partially permeable when its drainage path is clogged by the fine particles released from the soil layer. Therefore, it is of practical engineering importance to study the consolidation characteristics considering non-Darcian flow in foundations under semi-permeable boundaries.
In view of this, Gray [7] first proposed the definition of impeded drainage boundary in Darcy’s law, which can be described as:
R H u γ w = 1 γ w u z
where H is the height of the soil layer, u and γ w are the excess pore water pressure and the weight of water in order, and R is a dimensionless parameter that can reflect the drainage capacity of the impeded drainage boundary. When R tends to infinity or is equal to zero, Equation (1) can be simplified into the first kind of drainage boundary (i.e., complete drainage) or the second kind of drainage boundary (i.e., absolutely undrained). Thus, it is a more general boundary than the first and second kinds of boundary conditions. Using impeded drainage boundaries, Schiffman and Stein [8] studied the consolidation of soils with variable compressibility and permeability coefficients. Xie [9] conducted a fully explicit analytical solution of 1-D consolidation of double-layered soil under constant loads. Cai [10] used the Laplace transform to obtain the analytical solution of the viscoelastic soil layer under cyclic loadings. Recently, researchers [11,12] have also introduced the impede drainage boundaries into the consolidation theory of unsaturated soils and discussed the consolidation characteristics in 1-D and 2-D plane strain states.
So far, most of the studies, including those mentioned above, have been conducted under the assumption that the flow of pore water conforms to Darcy’s law. However, many penetration tests have shown [13,14,15,16,17,18] that for saturated soft clay or fine-grained soil with a low permeability coefficient under the low hydraulic slope, there are large deviations between actual flow and Darcy’s law at low hydraulic gradients. These nonlinear flow characteristics were first called non-Darcian flow by Hansbo [13], and he introduced a piecewise function to describe it. Based on this flow model, scholars [19,20] have successively analyzed the non-Darcian factors in 1-D or 2-D consolidation, and some reasonable explanations have been given to problems that cannot be explained by linear Darcy’s law. However, the above nonlinear flow model is described in the form of a piecewise function, which makes it complicated to determine the piecewise interval. Given this, Swartzendruber [21] proposed a non-Darcian flow model described by a continuous function, which is expressed as:
v = K [ i i 0 ( 1 e i / i 0 ) ]
where v is penetration velocity, K represents the permeability coefficient, i0 is the slope and intercept of the asymptotics of the v~i curve, and i0 is the non-Newtonian index. When i0 = 0, the upper formula degenerates into Darcian flow.
The expression of this flow model described by the non-Newtonian index is relatively simple, and it also fits well with the experimental data. Scholars [22,23] have analyzed the influence of the flow parameters on 1-D consolidation or viscoelastic consolidation. However, to the best of the knowledge of this study’s authors, these non-Darcy flow studies above ignored the impeded drainage boundary conditions, which is this paper’s starting point.
Based on the classical non-Darcian flow with non-Newtonian index, this paper deduces the finite difference expression of the 1-D consolidation equation of non-Darcian flow under impeded boundary drainage conditions. Then, the program of 1-D consolidation incorporating impeded boundaries is compiled. The method in this paper and the existing analytical solutions are used to verify the solution method’s rationality and the programming’s correctness. Finally, the program developed in this paper is used to analyze the 1-D consolidation characteristics with non-Darcian flow by incorporating a semi-permeable boundary.

2. Material Model and Methods

2.1. Problem Description

A simplified model of the 1-D consolidation incorporating impeded boundaries is shown in Figure 1. Assuming that the soil layer thickness of H has been fully consolidated under the stress of self-weight, the flow of pore water during consolidation is described in Equation (2). The thicknesses of the sand cushion at the top and the underlying layer at the bottom base are Lt and Lb, respectively, and the corresponding permeability coefficients are kb and kt, respectively. Other assumptions used in this study are parallel to those used in the consolidation theory of saturated soil [18,19,20]. The main assumptions are summarized as follows: (i) The soil layer is saturated and homogeneous; (ii) the load p(t) is an infinite uniform load applied instantaneously; (iii) the compression and flow are assumed to occur along the vertical direction only; and (iv) the coefficients of permeability is constant during the consolidation process. The consolidation deformation of the soil layer is small, and the influence of the secondary consolidation is omitted.
For the variable loading condition, the corresponding consolidation equation for 1-D consolidation of soft soil layer can be written as [20,22]:
K γ w 2 u z 2 ( 1 e 1 i 0 γ w u z ) = a 1 + e 0 ( u t p t )
where e0 is the initial void ratio and a is the compression coefficient of soil.
The initial conditions for this subject is:
u ( z , 0 ) = u 0 = p 0 ,   0 z H
The impeded drainage boundaries considering non-Darcian flow with non-Newtonian index can be rewritten as follows according to Equation (1):
When z = 0,
R t H u γ w = 1 γ w u z | z = 0 i 0 ( 1 e 1 i 0 γ w u z | z = 0 )
and when z = H,
R b H u γ w = 1 γ w u z | z = H i 0 ( 1 e 1 i 0 γ w u z | z = h )
where R t = k t H K L t , R b = k b H K L b , which can respectively reflect the drainage capacity of the top and bottom of the soft soil layer. It can be seen from Formula (5)–(6) that when i0 = 0, the above formula degenerates to Formula (1), which is the impeded drainage boundaries corresponding to Darcy’s law.
For the convenience of analysis, the following nondimensional variables are introduced:
Z = z H ,   U = u p u ,   T = C v t H 2 ,   I 0 = i 0 γ w H p u ,   P ( T ) = p ( t ) p u ,   C v = K ( 1 + e 0 ) a γ w
where Z is the nondimensional depth; Pu is the maximum applied load; and H is a reference length of soil thickness, which is nonconsistent with the drainage length in Terzaghi’s 1-D solution. U is the nondimensional pore pressure; T represents nondimensional time; I0 is nondimensional flow parameter; P(T) denotes nondimensional variable load; and Cv is the vertical consolidation coefficient of soil.
In this way, Equations (3)–(6) can be successively transformed into:
U j 1 , i + 1 2 U j , i + 1 + U j + 1 , i + 1 Δ Z 2 ( 1 e 1 I 0 U j + 1 , i + 1 U j 1 , i + 1 2 Δ Z ) = U j , i + 1 U j , i Δ T p i + 1 p i Δ T
U ( Z , 0 ) = U 0 = p 0 / p u ,   0 < Z < 1
R t U = U Z | Z = 0 I 0 ( 1 e 1 I 0 U Z | Z = 0 )
R b U = U Z | Z = H I 0 ( 1 e 1 I 0 U Z | Z = H )
where j is the discrete points along with the depth and i’ is time the axis discrete points.

2.2. Differential Iteration Format of The Governing Equation

The soil layer is evenly divided into the n − 1 layer from the top to the bottom base, with nodes numbered 1 to n. The time step is taken as Δt and numbered from 1. This section adopts an implicit format for linear partial differential equations, and the differential equation corresponding to the governing equation can be expressed as:
U j 1 , i + 1 2 U j , i + 1 + U j + 1 , i + 1 Δ Z 2 ( 1 e 1 I 0 U j + 1 , i + 1 U j 1 , i + 1 2 Δ Z ) = U j , i + 1 U j , i Δ T
The corresponding initial conditions and boundary conditions can be represented by discrete points in turn as
U j , 1 = P 0
( 1 Δ Z + R 1 ) U 1 , i + 1 + U 2 , i + 1 Δ Z = I 0 ( 1 e 1 I 0 U 2 , i + 1 U 1 , i + 1 Δ Z )
U k 1 , i + 1 Δ Z + ( R 2 + 1 Δ Z ) U k , i + 1 = I 0 ( 1 e 1 I 0 U k , i + 1 U k 1 , i + 1 Δ Z )
Thus, Equations (12)–(15) can be formed into a tridiagonal matrix equation set in the form of AX = B, and the elements of matrix A are expressions related to pore pressure, and it can be calculated and solved by an iterative method.

2.3. Degrees of Consolidation

For homogeneous elastic soils, the degree of consolidation defined by deformation and the degree of consolidation defined by pore pressure are equal. For the convenience of comparison and discussion, the consolidation degree Up, which reflects the dissipation of soil pore pressure, is introduced here:
U s = U p = p ( t ) p u = P i + 1 Δ z j = 1 n 1 U j , i + U j 1 , i 2

3. Verification of the Solution

When Rt = ∞, Rb = 0, the subject of this paper degenerates to a 1-D consolidation problem considering the non-Darcian flow described by the non-Newtonian index under the condition of a single-sided full-drainage boundary. The effect of single-stage linear loading on this consolidation problem has been studied by Li [22] using the C-N finite-difference format, where the dimensionless parameter n0 = 1 corresponds to the condition where the vertical additional stress is uniformly distributed along with the depth. To verify the validity of the algorithm in this paper for nonlinear flow, the linear loading duration Tvc = 0.1 and I0 = 3 are calculated according to the method in this paper, and the results are shown in Figure 2. For the convenience of comparison, the data corresponding to n0 = 1 in the literature [22] are also presented in Figure 2. Then, it can be found that the results of both are quite consistent
Xie [9] gave an analytical solution for the 1-D consolidation of double-layered soil under the condition of Darcy’s law considering impeded drainage boundaries. If the parameters of the double-layered soil are the same, it can degenerate into a single-layer soil problem. As mentioned above, when I0 = 0 and the external load p(t) is constant, the solution in the paper is simplified to the impeded drainage boundaries of Darcy’s law under constant load. To further verify the validity of the solution in this paper and the correctness of the programming, two calculation examples in the literature [9] are analyzed, using the analytical solution and the numerical method in this paper, and the corresponding results are shown in Figure 3 and Figure 4. In the two calculation examples, the thickness of the soil layer is taken as 5 m, the permeability coefficient K = 2 × 109 m/s, the modulus of compressibility is Es = 2 MPa, and the instantaneous load is 100 kPa. A total of 100 equal-length units are divided along with the depth, the time calculation step is taken as 1 day (i.e., T = 1.3824 × 103), and the iteration error is set as 10−6.
Figure 3 and Figure 4 respectively show the dissipation of pore pressure along with the depth under the single-sided and double-sided impeded drainage boundaries at 500 days. It can be seen from the two figures that the numerical solutions in this paper are quite consistent with the literature solutions under either the single-sided or double-sided impeded drainage boundaries, which fully verifies the rationality of the solution in this paper and the correctness of the programming.

4. Consolidation Characteristics Analysis

For the convenience of discussion, this paper only discusses the case where p(t) is a constant.
Figure 5 and Figure 6 show the comparisons of pore pressure of Darcy flow (I0 = 0) and non-Darcian flow rule (I0 = 0.5) under single-sided and double-sided impeded drainage boundaries, respectively. It can be seen from Figure 5 that when Rt = 100, the pore pressure in Darcy’s law is faster than that in non-Darcian flow described by the non-Newtonian index along with the full depth. In contrast, when Rt =0.5, 1, 5, a different law is exhibited: There is a critical depth phenomenon, and below the critical depth, the pore pressure of Darcy’s law dissipates faster than it does under non-Darcian flow. In contrast, above the critical depth, the opposite phenomenon occurs, that is, when considering the non-Darcian flow, the excess pore water pressure dissipates faster than it does under Darcy’s flow rule, and the critical depth decreases gradually with the increase of Rt.
In fact, it is not difficult to understand from the definition of Equation (5) that when I0 > 0, the initial pore pressure at the drainage boundary is less than I0 = 0. That is to say, when considering non-Darcian flow, the pore pressure near the impeded drainage boundaries dissipates faster. However, in essence, this is due to the impeded drainage boundary condition because for a larger Rt, such as Rt = 100, there does not exist critical depth. For a smaller Rt, such as Rt = 1 and Rt = 5, the dissipation value of pore pressure becomes small. Then, along the depth direction of the soil layer, the dissipation rate of the pore pressure is gradually delayed when considering non-Darcian flow, i.e., it is gradually slower than Darcy’s law. It can be further seen from Figure 5 that when Rt = 100, the initial pore pressures at Z = 0 under the conditions of Darcy flow and non-Darcian flow are both equal to 0, which means it already can be treated as a completely permeable boundary. That is, in this scenario, whereas previous scholars [22,23] concluded that the dissipation of pore pressure under the condition of non-Darcian flow is always slower than that of Darcy’s law along with the depth of the soil layer, the present authors in fact believe that this phenomenon is due to the special void structure and mineral composition of cohesive soil on the one hand and the viscous resistance of bound water in cohesive soil on the other hand, which leads to the macro phenomenon of considering nonlinear flow delays the dissipation process of pore pressure.
Figure 7 and Figure 8 show the comparison of the pore pressure along with the depth for different I0 values at the single-sided and the double-sided impeded drainage boundaries, respectively. It can be seen from the figures that when Rt is 1, 5, and 100, with the increase of flow parameter I0, the dissipation of pore pressure along the full depth becomes slower, i.e., the greater the value of I0 is, the more obvious the delay phenomenon of the pore pressure along with the depth is, which also indicates that the dissipation rate of pore pressure will be overestimated if the non-Darcy flow effect of flow is ignored. As also shown in Figure 7 and Figure 8, the larger the boundary parameter R is, the stronger the boundary drainage capacity of impeded drainage boundaries and the faster the consolidation rate of the soft clay ground under the same value of percolation parameter I0.
In order to investigate the influence of impeded drainage boundary parameters on the consolidation process, the curves of the dimensionless parameter I0 on the average consolidation degree under the single-sided and double-sided impeded drainage boundaries are depicted in Figure 9 and Figure 10, respectively. It can be seen from the two figures that both the impeded drainage boundary parameters and flow parameters have a significant influence on the consolidation process. Under the same boundary drainage parameters, the consolidation rate slows with the increase of I0. As mentioned above, whether it is a single-sided or a double-sided impeded drainage boundary, there is a critical depth phenomenon in the dissipation of pore pressure, i.e., above the critical depth, the dissipation of pore pressure along the depth is faster than that of Darcy’s law when the non-Darcian flow is considered. However, it can be seen from Figure 9 and Figure 10 that the overall dissipation rate of the foundation soil under the above parameters still shows that non-Darcian flow slows down the consolidation process, and the larger the value of I0, the more obvious the delay in the overall consolidation process. The authors suggest that the pore water pressure considering non-Darcian flow above the critical depth differs little from Darcy’s law. Then, the overall influence on the foundation soil is not enough to affect the overall consolidation rate of the soil layer after the pore pressure is averaged along with the depth.

5. Discussion of the Applicability of Darcy’s Law

Li [22] discussed that the common variation range of flow parameter i0 is concentrated in the range of 1–12. Considering its degradation form, it can be concluded that the general discussion range of i0 is 0–12. Additionally, combined with the expression of I0, it is believed that its variation range is generally 0–20 [22] and that its common range is 0–5 [20,24,25]. The thickness of foundation soil is generally 0–15 m, and the overburden load usually varies in the range of 100–200 kPa, and the common variation range of I0 is summarized here. In this part of the error analysis, the variation range of I0 is 0–15.
From Figure 9 and Figure 10, it can be concluded that when I0 is small enough, the error caused by Darcy’s flow is acceptable in engineering practice; the deviation errors of Darcy’s law and non-Darcian flow with non-Newtonian index are shown in Figure 11 and Figure 12 (lower I0) and Figure 13 and Figure 14 (larger I0) under a single-sided and a double-sided impeded drainage boundary, respectively. As shown in Figure 11, under the condition of a one-sided impeded drainage boundary, when I0 = 0.1, the maximum deviation of the non-Darcian flow described by the non-Newtonian index with Darcy’s flow is less than 4%, and when I0 = 0.2, the maximum deviation of the two cases is less than 8%. Under the condition of a double-sided impeded drainage boundary, when I0 = 0.1, the maximum deviation of Darcy’s flow is about 2%, and when I0 = 0.5, the deviation is about 8.5%. As shown in Figure 13 and Figure 14, under the condition of one-sided and double-sided impeded drainage boundaries, when I0 = 12, the maximum deviations of the non-Darcian flow with Darcy’s flow are more than 56.8% and 44.9%, respectively, and when I0 = 15, the maximum deviation of the two cases exceeded 59.5% and 48.0% in sequence. Therefore, when the soil layer is very thin and the load acting on the soil layer is large (at this time, the I0 value is very small), it is acceptable to apply Darcy’s law to calculate consolidation deformation in engineering. On the contrary, when the soft soil layer is very thick, and the load acting on it is small, the consolidation deformation is calculated by non-Darcian flow to avoid causing great calculation error, especially when the time factor T changes from 0.1 to 1.5.

6. Conclusions

(1)
Whether it is a single-sided or a double-sided impeded drainage boundary, there will be a critical depth phenomenon when considering non-Darcian flow incorporating impeded drainage boundaries.
(2)
Below the critical depth, the distribution of pore pressure along the depth considering non-Darcian flow with the non-Newtonian index will be slower than that of Darcy’s law, while above the critical depth, the pore pressure will dissipate faster when the non-Darcian flow rule is considered. When the impeded drainage parameter Rt or Rb increases, the critical depth gradually decreases, but the effect of this above-critical-depth phenomenon on the overall dissipation of pore pressure is not enough to affect the overall consolidation rate of the soil layer, that is, the overall dissipation rate of the soft soil layer still reflects that the overall consolidation process is delayed when considering non-Darcian flow.
(3)
Both the impeded drainage boundary parameters and flow parameters have obvious effects on the 1-D consolidation of soft soils. The larger the impeded drainage boundary parameters are, the faster the overall consolidation rate of the soft clay ground is, and when the drainage is greater than 100, it no longer has a significant effect on the consolidation process; that is, the boundary can be regarded as full drainage. Considering non-Darcian flow will delay the consolidation rate of the soft soils, and the larger the nondimensional flow parameter I0 is, the more obvious this delay phenomenon will be.
(4)
When the value of I0 is small (less than 0.1 under the condition of a one-sided impeded drainage boundary and less than 0.2 under the condition of a double-sided impeded drainage boundary), the maximum deviation caused by Darcy’s law for consolidation calculation is less than 5%.

Author Contributions

Conceptualization, Z.X. and P.C.; methodology, W.C.; writing—original draft preparation, Z.X. and P.C.; writing and editing, Y.W. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Foundation of China (Grant Nos. 51879104; 5207080378; 51909086), Hunan Natural Science Research Fund (Grant No. 51879104), Hunan Provincial Department of Education Fund (Grant Nos. 20C0371; 18C0505), National College Students Innovation and Entrepreneurship Fund (Grant No. 202111527033), The support is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Terzaghi, K. Theoretical Soil Mechanics; John Wiley: New York, NY, USA, 1943. [Google Scholar] [CrossRef]
  2. Schiffman, R.L. Consolidation of soil under time-dependent loading and varying permeability. In Highway Research Board Proceedings; National Academy of Sciences-National Research Council: Washington, DC, USA.
  3. Baligh, M.M.; Levadoux, J.N. Consolidation theory for cyclic loading. J. Geotech. Eng. Divis. 1978, 104, 415–431. [Google Scholar] [CrossRef]
  4. Xie, K.-H.; Xie, X.-Y.; Li, X.-B. Analytical theory for one-dimensional consolidation of clayey soils exhibiting rheological characteristics under time-dependent loading. Int. J. Numer. Anal. Methods Géoméch. 2008, 32, 1833–1855. [Google Scholar] [CrossRef]
  5. Liu, J.; Griffiths, D. A general solution for 1-D consolidation induced by depth-and time-dependent changes in stress. Geotechnique 2015, 65, 66–72. [Google Scholar] [CrossRef]
  6. Mesri, G. One-dimensional consolidation of a clay layer with impeded drainage boundaries. Water Resour. Res. 1973, 9, 1090–1093. [Google Scholar] [CrossRef]
  7. Gray, H. Simultaneous Consolidation contiguous layers of unlike compressible soil. Trans. Am. Soc. Civ. Eng. 1945, 110, 1327–1356. [Google Scholar]
  8. Schiffman, R.L.; Stein, J.R. One-dimensional consolidation of layered systems. J. Soil Mech. Found. Div. 1970, 96, 1499–1504. [Google Scholar] [CrossRef]
  9. Xie, K.-H.; Xie, X.-Y.; Gao, X. Theory of one dimensional consolidation of two-layered soil with partially drained boundaries. Comput. Geotech. 1999, 24, 265–278. [Google Scholar] [CrossRef]
  10. Cai, Y.Q.; Liang, X.; Wu, S.M. One-dimensional consolidation of layered soils with impeded boundaries under time-dependent loadings. Appl. Math. Mech. 2004, 25, 937–944. [Google Scholar] [CrossRef]
  11. Wang, L.; Sun, D.A.; Li, L.Z.; Xu, Y. Semi-analytical solutions to one-dimensional consolidation for unsaturated soils with symmetric semi-permeable drainage boundary. Comput. Geotech. 2017, 89, 71–80. [Google Scholar] [CrossRef]
  12. Wang, L.; Xu, Y.; Xia, X.; Zhou, A. Semi-analytical solutions to the two-dimensional plane strain consolidation for unsaturated soil with the lateral semi-permeable drainage boundary under time-dependent loading. Comput. Geotech. 2020, 124, 103562. [Google Scholar] [CrossRef]
  13. Hansbo, S. Consolidation of clay, with special reference to influence of vertical sand drains. Swed. Geotech. Inst. Proceeding 1960, 18, 45–50. [Google Scholar]
  14. Olsen H, W. Osmosis: A cause of apparent deviation from Darcy’s law. Can. Geotech. J. 1985, 22, 238–241. [Google Scholar] [CrossRef]
  15. Ing, T.C.; Xiaoyan, N. Coupled consolidation theory with non-Darcian flow. Comput. Geotech. 2002, 29, 169–209. [Google Scholar] [CrossRef]
  16. Hansbo, S. Consolidation equation valid for both Darcian and non-Darcian flow. Geotechnique 2001, 51, 51–54. [Google Scholar] [CrossRef]
  17. Deng, Y.-B.; Liu, G.-B.; Zheng, R.-Y.; Xie, K.-H. Finite element analysis of biot’s consolidation with a coupled nonlinear flow model. Math. Probl. Eng. 2016, 2016, 1–13. [Google Scholar] [CrossRef]
  18. Indraratna, B.; Kianfar, K.; Rujikiatkamjorn, C. Radial consolidation model incorporating the effects of vacuum preloading and non-Darcian flow. Geotechnique. 2013, 63, 1060–1073. [Google Scholar] [CrossRef] [Green Version]
  19. Zhao, X.; Gong, W. Model for large strain consolidation with non-Darcian flow described by a flow exponent and threshold gradient. Int. J. Numer. Anal. Methods Géoméch 2019, 43, 2251–2269. [Google Scholar] [CrossRef]
  20. Liu, Z.; Cui, P.; Zhang, J.; Xia, Y. Consolidation analysis of ideal sand-drained ground with fractional-derivative Merchant model and non-Darcian flow described by non-Newtonian index. Math. Probl. Eng. 2019, 2019, 1–12. [Google Scholar] [CrossRef]
  21. Swartzendruber, D. Modification of Darcyʼs law for the flow of water in soils. Soil Sci. 1962, 93, 22–29. [Google Scholar] [CrossRef]
  22. Li, C.-X.; Xie, K.-H.; Lu, M.-M.; Wang, K. Analysis of one-dimensional consolidation with non-Darcy flow described by non-Newtonian index. Rock Soil Mech. 2011, 32, 281–287. (In Chinese) [Google Scholar] [CrossRef]
  23. Xu, C.; Liu, Z.; Zhang, J.; Huang, J. Analysis of large-strain elastic viscoplastic consolidation for soft clay with vertical drains considering non-Darcian flow. Appl. Math. Model. 2021, 92, 770–784. [Google Scholar] [CrossRef]
  24. Cui, P.; Cao, W.; Xu, Z.; Wei, Y.; Li, H. Analysis of one-dimensional nonlinear rheological consolidation of saturated clay considering self-weight stress and Swartzendruber’s flow rule. Int. J. Numer. Anal. Methods Géoméch. 2022, 46, 1205–1223. [Google Scholar] [CrossRef]
  25. Cui, P.L.; Liu, Z.Y.; Zhang, J.C.; Fan, J.C. Analysis of one-dimensional rheological consolidation of double-layered soil with fractional-derivative Merchant model and non-Darcian flow described by non-Newtonian index. J. Cent. South Univ. 2021, 28, 284–296. [Google Scholar] [CrossRef]
Figure 1. Sketch of 1-D consolidation with impeded boundaries.
Figure 1. Sketch of 1-D consolidation with impeded boundaries.
Water 14 01740 g001
Figure 2. The verification of non-Darcian flow described by the non-Newtonian index [22].
Figure 2. The verification of non-Darcian flow described by the non-Newtonian index [22].
Water 14 01740 g002
Figure 3. The verification of the single-sided impeded drainage boundaries at 500 days (i.e., corresponding to T = 0.6912) [9].
Figure 3. The verification of the single-sided impeded drainage boundaries at 500 days (i.e., corresponding to T = 0.6912) [9].
Water 14 01740 g003
Figure 4. The verification of the double-sided impeded drainage boundaries at 500 days [9].
Figure 4. The verification of the double-sided impeded drainage boundaries at 500 days [9].
Water 14 01740 g004
Figure 5. Darcy’s law and non-Darcian flows under a single-sided impeded drainage boundary (Rb = 0, T = 0.5).
Figure 5. Darcy’s law and non-Darcian flows under a single-sided impeded drainage boundary (Rb = 0, T = 0.5).
Water 14 01740 g005
Figure 6. Darcy and non-Darcian flows under a double-sided impeded drainage boundary (Rb = Rt, T = 0.2).
Figure 6. Darcy and non-Darcian flows under a double-sided impeded drainage boundary (Rb = Rt, T = 0.2).
Water 14 01740 g006
Figure 7. The influence of Rt on the dissipation of pore water pressure Rb = 0, T = 0.5.
Figure 7. The influence of Rt on the dissipation of pore water pressure Rb = 0, T = 0.5.
Water 14 01740 g007
Figure 8. The influence of Rt on the dissipation of pore water pressure Rb = Rt, T = 0.2.
Figure 8. The influence of Rt on the dissipation of pore water pressure Rb = Rt, T = 0.2.
Water 14 01740 g008
Figure 9. The influence of I0 on the average degrees of consolidation under a single-side impeded drainage boundary.
Figure 9. The influence of I0 on the average degrees of consolidation under a single-side impeded drainage boundary.
Water 14 01740 g009
Figure 10. The influence of I0 on the average degrees of consolidation under a double-side impeded drainage boundaries.
Figure 10. The influence of I0 on the average degrees of consolidation under a double-side impeded drainage boundaries.
Water 14 01740 g010
Figure 11. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a single-side impeded drainage boundary with lower I0.
Figure 11. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a single-side impeded drainage boundary with lower I0.
Water 14 01740 g011
Figure 12. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a double-side impeded drainage boundary with lower I0.
Figure 12. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a double-side impeded drainage boundary with lower I0.
Water 14 01740 g012
Figure 13. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a single-side impeded drainage boundary with larger I0.
Figure 13. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a single-side impeded drainage boundary with larger I0.
Water 14 01740 g013
Figure 14. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a double-side impeded drainage boundary with larger I0.
Figure 14. Deviations in the consolidation degree between Darcy’s law and non-Darcian flow under a double-side impeded drainage boundary with larger I0.
Water 14 01740 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, Z.; Cao, W.; Cui, P.; Li, H.; Wei, Y. Analysis of One-Dimensional Consolidation Considering Non-Darcian Flow Described by Non-Newtonian Index Incorporating Impeded Drainage Boundaries. Water 2022, 14, 1740. https://doi.org/10.3390/w14111740

AMA Style

Xu Z, Cao W, Cui P, Li H, Wei Y. Analysis of One-Dimensional Consolidation Considering Non-Darcian Flow Described by Non-Newtonian Index Incorporating Impeded Drainage Boundaries. Water. 2022; 14(11):1740. https://doi.org/10.3390/w14111740

Chicago/Turabian Style

Xu, Zan, Wengui Cao, Penglu Cui, Huixin Li, and Yunbo Wei. 2022. "Analysis of One-Dimensional Consolidation Considering Non-Darcian Flow Described by Non-Newtonian Index Incorporating Impeded Drainage Boundaries" Water 14, no. 11: 1740. https://doi.org/10.3390/w14111740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop