Explicit approximate analytical solutions of seepage‐deformation in unsaturated soils

An approximate analytical solution is presented for the coupled seepage and deformation problem of unsaturated soils. Because of the matric suction dependence of both saturation and permeability coefficient, the coupled governing equations are strongly nonlinear. To obtain an analytical solution, these coupled governing equations are linearized and analytically solved for a specified saturation using the eigenfunction method. Then, the obtained analytical solutions are extended to the entire saturation range. Comparison between the current solution and the previous theoretical solution indicates that the proposed solution yields excellent results. Due to its analytical nature, the proposed procedure can be effectively used to obtain the solution of the coupled seepage and deformation of unsaturated soils.     

A simple analytical solution to one‐dimensional consolidation for unsaturated soils

This paper presents a simple analytical solution to Fredlund and Hasan's one‐dimensional (1‐D) consolidation theory for unsaturated soils. The coefficients of permeability and volume change for unsaturated soils are assumed to remain constant throughout the consolidation process. The mathematical expression of the present solution is much simpler compared with the previous available solutions in the literature. Two new variables are introduced to transform the two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved with standard mathematical formulas. It is shown that the present analytical solution can be degenerated into that of Terzaghi consolidation for fully saturated condition. The analytical solutions to 1‐D consolidation of an unsaturated soil subjected to instantaneous loading, ramp loading, and exponential loading, for different drainage conditions and initial pore pressure conditions, are summarized in tables for ease of use by practical engineers. In the case studies, the analytical results show good agreement with the available analytical solution in the literature. The consolidation behaviors of unsaturated soils are investigated. The average degree of consolidation at different loading patterns and drainage conditions is presented. The pore‐water pressure isochrones for two different drainage conditions and three initial pore pressure distributions are presented and discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Analytical solution to 1D coupled water infiltration and deformation in two‐layer unsaturated soils

An analytical solution to 1D coupled water infiltration and deformation in layered soils is derived using a Laplace transformation. Coupling between seepage and deformation, and initial conditions defined by arbitrary continuous pore‐water pressure distributions are considered. The analytical solutions describe the transient pore‐water pressure distributions during 1D, vertical infiltration toward the water table through two‐layer unsaturated soils. The nonlinear coupled formulations are first linearized and transformed into a form that is solvable using a Laplace transformation. The solutions provide a reliable means of comparing the accuracy of various numerical methods. Parameters considered in the coupled analysis include the saturated permeability (k_s), desaturation coefficient (α), and saturated volumetric water content (θ_s) of each soil layer, and antecedent and subsequent rainfall infiltration rates. The analytical solution demonstrates that the coupling of seepage and deformation plays an important role in water infiltration in layered unsaturated soils. A smaller value of α or a smaller absolute value of the elastic modulus of the soil with respect to a change in soil suction (H) for layered unsaturated soils means more marked coupling effect. A smaller absolute value of H of the upper layer soil also tends to cause more marked coupling effect. A large difference between the saturated coefficients of permeability for the top and bottom soil layers leads to reduced rainfall infiltration into the deep soil layer. The initial conditions also play a significant role in the pore‐water pressure redistribution and coupling effect. Copyright © 2011 John Wiley & Sons, Ltd.     

Three-dimensional heat,moisture and air transfer in unsaturated soils

A new three-dimensional numerical model of coupled heat, moisture and air transfer in unsaturated soil is presented. In particular, the model accommodates moisture transfer in the form of liquid and vapour flow and heat transfer arising from conduction, convection and latent heat of vaporization. The bulk flow of dry air and the movement of air in a dissolved state are also included. The theoretical basis of the model, the finite element solution of the spatial terms and finite difference solution of the temporal terms are briefly presented. Attention is focused on the verification of the new numerical solution. This is achieved via comparisons with independent solutions of heat, moisture and air transfer in an unsaturated soil. The physical problem considered includes the highly non-linear hydraulic properties of sand. Thermal conductivity is also included as a function of soil moisture content. Excellent correlation of results is shown thus providing confidence in the new model. The new model is also applied to a number of test cases which illustrate the need for the development of a model which can fully include three-dimensional behaviour. In particular, three applications are presented each increasing in complexity. The first application illustrates three-dimensional heat transfer. This particular application is verified against existing commercial finite element software. Subsequent applications serve to illustrate how the coupled processes of heat moisture and air transfer combine to yield three-dimensional problems even within a simple geometric domain. Visualization of three-dimensional results is also addressed. © 1998 by John Wiley & Sons, Ltd.     

A closed form analytical solution for two‐dimensional plane strain consolidation of unsaturated soil stratum

This paper discusses the excess pore‐air and pore‐water pressure dissipations and the average degree of consolidation in the 2D plane strain consolidation of an unsaturated soil stratum using eigenfunction expansion and Laplace transformation techniques. In this study, the application of a constant external loading on a soil surface is assumed to immediately generate uniformly or linearly distributed initial excess pore pressures. The general solutions consisting of eigenfunctions and eigenvalues are first proposed. The Laplace transform is then applied to convert the time variable t in partial differential equations into the Laplace complex argument s. Once the domain is obtained, a simplified set of equations with variable s can be achieved. The final analytical solutions can be computed by taking a Laplace inverse. The proposed equations predict the two‐dimensional consolidation behaviour of an unsaturated soil stratum capturing the uniformly and linearly distributed initial excess pore pressures. This study investigates the effects of isotropic and anisotropic permeability conditions on variations of excess pore pressures and the average degree of consolidation. Additionally, isochrones of excess pore pressures along vertical and horizontal directions are presented. It is found that the initial distribution of pore pressures, varying with depth, results in considerable effects on the pore‐water pressure dissipation rate whilst it has insignificant effects on the excess pore‐air pressure dissipation rate. Copyright © 2015 John Wiley & Sons, Ltd.     

A semi‐analytical approach for modelling heat transfer during salt‐cavern leaching process

The use of gas‐storage caverns in salt formations is a growing industry that continues to gain momentum because it allows gas to be injected and withdrawn at high rates compared with other underground gas‐storage systems such as porous rock systems. In order to predict cavern production performances, cavern thermodynamics behaviour must be studied by higher accuracy approaches. This behaviour is extremely related to the temperature distribution in the surrounding formations. During the leaching process, the thermal equilibrium of the rock salt surrounding the cavern is extensively disrupted. The purpose of this paper is to study the heat transfer problem during the leaching process and to develop a thermal model that can be easily used in field applications. The results of this work will be the input data for the prediction of the gas temperature and pressure during cavern gas‐storage operation phase. Moreover, the developed model can find its use in the design of salt caverns since it can be coupled with geometrical modelling of salt dissolution codes. Copyright © 2008 John Wiley & Sons, Ltd.     

Analytical solution for one‐dimensional consolidation of unsaturated soils using eigenfunction expansion method

This paper introduces an exact analytical solution for governing flow equations for one‐dimensional consolidation in unsaturated soil stratum using the techniques of eigenfunction expansion and Laplace transformation. The homogeneous boundary conditions adopted in this study are as follows: (i) a one‐way drainage system of homogenous soils, in which the top surface is considered as permeable to air and water, whereas the base is an impervious bedrock; and (ii) a two‐way drainage system where both soil ends allow free dissipation of pore‐air and pore‐water pressures. In addition, the analytical development adopts initial conditions capturing both uniform and linear distributions of the initial excess pore pressures within the soil stratum. Eigenfunctions and eigenvalues are parts of the general solution and can be obtained based on the proposed boundary conditions. Besides, the Laplace transform method is adopted to solve the first‐order differential equations. Once equations with transformed domain are all obtained, the final solutions, which are proposed to be functions of time and depth, can be achieved by taking an inverse Laplace transform. To verify the proposed solution, two worked examples are provided to present the consolidation characteristics of unsaturated soils based on the proposed method. The validation of the recent results against other existing analytical solutions is graphically demonstrated. Copyright © 2013 John Wiley & Sons, Ltd.     

Fast planar slides. A closed‐form thermo‐hydro‐mechanical solution

Heat‐induced excess pore pressures on the failure surface of a planar slide have been calculated by solving the mass and heat balance equations on the shear band. The set of differential equations and the equation of motion of the slide have been solved in closed form for the case of incompressible fluid and incompressible soil skeleton. The solution describes the accelerated motion of the slide. It has been compared with the numerical solution when soil and water stiffness terms are not disregarded. A case study, based on a well‐known translational slide (Cortes slide) has been solved. Numerical and analytical solutions are compared. Results of a sensitivity analysis indicate that the permeability of the shear band is the key parameter to control the onset of a rapid motion. For a band permeability above a threshold value, in the vicinity of 10^?15m² (10^?8m/s), fast accelerated motions are very unlikely. Copyright © 2009 John Wiley & Sons, Ltd.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Analytical solution to 1D coupled water infiltration and deformation in unsaturated soils

An analytical solution to 1D coupled water infiltration and deformation is derived using a Fourier integral transform. Exponential functional forms are used to represent the hydraulic conductivity–pore‐water pressure relationship and the soil‐water characteristic curve. Fredlund's incremental‐linear constitutive model for unsaturated soils is adopted. The analytical solution considers arbitrary initial pore‐water pressure distributions and flux and pressure boundary conditions. The corresponding analytical solutions to coupled steady‐state problems are also obtained. The analytical solutions demonstrate that the coupling of seepage and deformation plays an important role in water infiltration in unsaturated soils. In the early stages of infiltration, the difference between uncoupled and coupled conditions becomes marked over time, and in late stages, the difference caused by the coupling effects diminishes toward the steady state. The difference between the uncoupled and coupled conditions increases with decreasing desaturation coefficient (α). Pore‐water pressure or deformation changes caused by the coupling effects are mainly controlled by the degree of soil volume change due to a change in soil suction (H). The smaller the absolute value of H, the greater the effect of coupling on the infiltration and deformation. The ratio of rainfall intensity to saturated permeability (q/k_s) also has a strong influence on the coupled seepage and deformation. Copyright © 2008 John Wiley & Sons, Ltd.     

An analytical solution for rainfall infiltration into an unsaturated infinite slope and its application to slope stability analysis

Surficial slope failures in residual soils are common in tropical and subtropical regions as a result of rainfall infiltration. This study develops an analytical solution for simulating rainfall infiltration into an infinite unsaturated soil slope. The analytical solution is based on the general partial differential equation for water flow through unsaturated soils. It can accept soil–water characteristic curve and unsaturated permeability function of the exponential form into account. Numerical simulations are conducted to verify the assumptions of the analytical solution and demonstrate that the proposed analytical solution is acceptable for the coarse soils with low air entry values. The pore‐water pressure (pwp) distributions obtained from the analytical solution can be incorporated into a limit equilibrium method to do infinite slope stability analysis for a rain‐induced shallow slip. The analysis method takes into account the influence of the water content change on unit weight and hence on factor of safety. A series of analytical parametric analyses have been performed using the developed model. The analyses indicate that when the residual soil slope, consisting of a completely decomposed granite layer underlain by a less permeable layer, is subjected to a continuous heavy rainfall, the loss of negative pwp and the reduction in factor of safety were found to be most significant for the shallow soil layer and during the first 12 h. The antecedent and subsequent rainfall intensity, depth of a less permeable layer and slope angle all have a significant influence on the pwp response and hence the slope stability. Copyright © 2012 John Wiley & Sons, Ltd.     

非饱和介质中热能传输及水分迁移的数值积分解

在给出非饱和介质热能-水分传输的耦合质量控制方程和基于Fourier热传导定律的热能平衡方程的基础上,对热能传输及水分迁移的基本特征和机理进行了分析。其中,考虑了温度势、吸力势和重力势的耦合作用影响。给出有热源时控制方程的简化形式,并对半无限体自由表面作用平面热源条件下介质内非稳态温度场、体积含水率分布场进行数值积分求解。利用这些解答给出常热源强度和变热源强度两种情况下,温度场随时间的变化特征以及水分迁移的演化过程,并分析了重力项对计算结果的影响。     

Coupled solution of heat and moisture flow in unsaturated clay barriers in a repository geometry

An analytical solution in the Laplace transform domain is obtained for the transient heat and moisture transport in an unsaturated clay buffer with a geometry simulating repository conditions. A numerical inversion scheme based on Crump's method is used to obtain the time‐domain solution. The coupled effect of thermally driven moisture transport is especially investigated because of its importance to alter the flow field in low‐permeability buffers. The practical background is based on the case of an engineering bentonite barrier placed in a drift excavated in rock in the context of underground disposal of high‐level radioactive waste. Parametric study has been performed to assess the effects of dimensionless geometry and material parameters on flow field. Despite the simplified assumptions required in order to obtain analytical expressions, the results incorporate the main mechanisms involved in the coupled thermo‐hydraulic (T–H) problem, and they may be eventually used for validation purposes. Copyright © 2006 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

均质非饱和土一维稳态流水头垂直分布的解析解

渗流场水头分布计算是进行渗流量和渗流水力坡降计算的基础,准确、有效地求取渗流场水头分布是渗流计算的关键环节。对均质非饱和土体一维稳态流的流动方程进行分析,考虑到渗透系数是与基质吸力相关的函数,通过数学变换,给出了稳定渗流场的解析通式,并基于渗透性函数中的Gardner模型,给出了非饱和土一维稳态流水头垂直分布的解析解。该解析通式表明,均质非饱和土一维稳态流水头垂直分布主要受地表水头、深度和流动率3个因素控制。分别计算了一维稳态蒸发条件下粉土和黏土两种典型土类水头沿垂直方向的分布。计算结果表明：稳态蒸发条件下粉土层和黏土层内的水头分布表现出相似的变化规律,即自地表至地下水位处随着土层深度的增加,水头分布呈现出加速递减的趋势;在相同的蒸发条件下,对于相同深度处的黏土和粉土而言,黏土层内水头更高些;对同一种土类而言,在较大的蒸发状态下同一深度处土层内水头更高。反之,则较低。     

单层非饱和多孔介质一维瞬态响应半解析解

基于Zienkiewicz提出的非饱和多孔介质波动理论,考虑两相流体和固体颗粒的压缩性以及惯性、黏滞和机械耦合作用,采用半解析的方法获得了一类典型边界条件下单层非饱和多孔介质一维瞬态响应解。首先推导出无量纲化后以位移表示的控制方程,并将其写成矩阵形式;然后,将边界条件齐次化,求解控制方程所对应的特征值问题,得到了满足齐次边界条件的特征值和相对应的特征函数。根据变异系数法并利用特征函数的正交性,得到了一系列仅黏滞耦合的关于时间的二阶常微分方程及相应的初始条件。在此基础上,运用精细时程积分法给出了常微分方程组的数值解。最后,通过若干算例验证了结果的正确性并探讨了单层非饱和多孔介质一维瞬态动力响应的特点。该方法可推广应用于其他典型的边界条件。     

A semi‐analytical modeling approach for three‐dimensional heat transfer in sparsely fractured rocks with water flow and distributed heat source

A semi‐analytical approach is developed for modeling 3D heat transfer in sparsely fractured rocks with prescribed water flow and heat source. The governing differential equations are formulated, and the corresponding integral equations over the fracture faces and the distributed heat source are established in the Laplace transformed domain using the Green function method with local systems of coordinates. The algebraic equations of the Laplace transformed temperatures of water in the fractures are formed by dividing the integrals into elemental ones; in particular, the fracture faces are discretized into rectangular elements, over which the integrations are carried out either analytically for singular integrals when the base point is involved or numerically for regular integrals when otherwise. The solutions of the algebraic equations are inverted numerically to obtain the real‐time temperatures of water in the fractures, which may be employed to calculate the temperatures at prescribed locations of the rock matrix. Three example calculations are presented to illustrate the workability of the developed approach. The calculations found that water flux in the fractures may decrease the rate of temperature rise in regions close to the distributed heat source and increase the rate of temperature rise in regions downstream away from the distributed heat source and that the temperature distribution and evolvement in a sparsely fractured rock mass may be significantly influenced by water flow exchange at intersection of fractures. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical analysis of a one‐dimensional infiltration problem in unsaturated soil by a seepage–deformation coupled method

A multiphase coupled elasto‐viscoplastic finite element analysis formulation, based on the theory of porous media, is used to describe the rainfall infiltration process into a one‐dimensional soil column. Using this framework, we have numerically analyzed the generation of pore water pressure and deformations when rainfall is applied to the soil. A parametric study, including rainfall intensity, soil–water characteristic curves, and permeability, is carried out to observe their influence on the changes in pore water pressure and volumetric strain. From the numerical results, it is shown that the generation of pore water pressure and volumetric strain is mainly controlled by material parameters α and n′ that describe the soil–water characteristic curve. A comparison with the laboratory results shows that the proposed method can describe very well the characteristics observed during the experiments of one‐dimensional water infiltration into a layered unsaturated soil column. Copyright © 2010 John Wiley & Sons, Ltd.     

Exact solutions for one‐dimensional consolidation of single‐layer unsaturated soil

Based on the Fredlund consolidation theory of unsaturated soil, exact solutions of the governing equations for one‐dimensional consolidation of single‐layer unsaturated soil are presented, in which the water permeability and air transmission are assumed to be constants. The general solution of two coupled homogeneous governing equations is first obtained. This general solution is expressed in terms of two functions psi₁ and ψ₂, where ψ₁ and ψ₂, respectively, satisfy two second‐order partial differential equations, which are in the same form. Using the method of separation of variables, the two partial differential equations are solved and exact solutions for three typical homogeneous boundary conditions are obtained. To obtain exact solutions of nonhomogeneous governing equations with three typical nonhomogeneous boundary conditions, the nonhomogeneous boundary conditions are first transformed into homogeneous boundary conditions. Then according to the method of undetermined coefficients and exact solutions of homogenous governing equations, the series form exact solutions are put forward. The validity of the proposed exact solutions is verified against other analytical solutions in the literature. Copyright © 2011 John Wiley & Sons, Ltd.