			Home About us Contact

Unsaturated Soils (unsaturated + soil)

Distribution by Scientific Domains

Engineering	91%

Selected Abstracts

Dynamics of unsaturated soils using various finite element formulations

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 5 2009
Nadarajah Ravichandran
Abstract Unsaturated soils are three-phase porous media consisting of a solid skeleton, pore liquid, and pore gas. The coupled mathematical equations representing the dynamics of unsaturated soils can be derived based on the theory of mixtures. Solution of these fully coupled governing equations for unsaturated soils requires tremendous computational resources because three individual phases and interactions between them have to be taken into account. The fully coupled equations governing the dynamics of unsaturated soils are first presented and then two finite element formulations of the governing equations are presented and implemented within a finite element framework. The finite element implementation of all the terms in the governing equations results in the complete formulation and is solved for the first time in this paper. A computationally efficient reduced formulation is obtained by neglecting the relative accelerations and velocities of liquid and gas in the governing equations to investigate the effects of fluid flow in the overall behavior. These two formulations are used to simulate the behavior of an unsaturated silty soil embankment subjected to base shaking and compared with the results from another commonly used partially reduced formulation that neglects the relative accelerations, but takes into account the relative velocities. The stress,strain response of the solid skeleton is modeled as both elastic and elastoplastic in all three analyses. In the elastic analyses no permanent deformations are predicted and the displacements of the partially reduced formulation are in between those of the reduced and complete formulations. The frequency of vibration of the complete formulation in the elastic analysis is closer to the predominant frequency of the base motion and smaller than the frequencies of vibration of the other two analyses. Proper consideration of damping due to fluid flows in the complete formulation is the likely reason for this difference. Permanent deformations are predicted by all three formulations for the elastoplastic analyses. The complete formulation, however, predicts reductions in pore fluid pressures following strong shaking resulting in somewhat smaller displacements than the reduced formulation. The results from complete and reduced formulations are otherwise comparable for elastoplastic analyses. For the elastoplastic analysis, the partially reduced formulation leads to stiffer response than the other two formulations. The likely reason for this stiffer response in the elastoplastic analysis is the interpolation scheme (linear displacement and linear pore fluid pressures) used in the finite element implementation of the partially reduced formulation. Copyright © 2008 John Wiley & Sons, Ltd. [source]

Transformations of runoff chemistry in the Arctic tundra, Northwest Territories, Canada

HYDROLOGICAL PROCESSES, Issue 14 2006
W. L. Quinton
Abstract The transformation of snowmelt water chemical composition during melt, elution and runoff in an Arctic tundra basin is investigated. The chemistry of the water flowing along pathways from the surface of melting snow to the 95·5 ha basin outlet is related to relevant hydrological processes. In so doing, this paper offers physically based explanations for the transformation of major ion concentrations and loads of runoff water associated with snowmelt and rainfall along hydrological pathways to the stream outlet. Late-lying snowdrifts were found to influence the ion chemistry in adjacent reaches of the stream channel greatly. As the initial pulse of ion-rich melt water drained from the snowdrift and was conveyed through hillslope flowpaths, the concentrations of most ions increased, and the duration of the peak ionic pulse lengthened. Over the first 3 m of overland flow, the concentrations of all ions except for NO increased by one to two orders of magnitude, with the largest increase for K+, Ca2+ and Mg2+. This was roughly equivalent to the concentration increase that resulted from percolation of relatively dilute water through 0·25 m of unsaturated soil. The Na+ and Cl, were the dominant ions in snowmelt water, whereas Ca2+ and Mg2+ dominated the hillslope runoff. On slopes below a large melting snowdrift, ion concentrations of melt water flowing in the saturated layer of the soil were very similar to the relatively dilute concentrations found in surface runoff. However, once the snowdrift ablated, ion concentrations of subsurface flow increased above parent melt-water concentrations. Three seasonally characteristic hydrochemical regimes were identified in a stream reach adjacent to late-lying snowdrifts. In the first two stages, the water chemistry in the stream channel strongly resembled the hillslope drainage water. In the third stage, in-stream geochemical processes, including the weathering/ion exchange of Ca2+ and Mg2+, were the main control of streamwater chemistry. Copyright © 2006 John Wiley & Sons, Ltd. [source]

Steady infiltration from buried point source into heterogeneous cross-anisotropic unsaturated soil

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 10 2004
G. J. Chen
Abstract The paper presents the analytical solution for the steady-state infiltration from a buried point source into two types of heterogeneous cross-anisotropic unsaturated half-spaces. In the first case, the heterogeneity of the soil is modelled by an exponential relationship between the hydraulic conductivity and the soil depth. In the second case, the heterogeneous soil is represented by a multilayered half-space where each layer is homogeneous. The hydraulic conductivity varies exponentially with moisture potential and this leads to the linearization of the Richards equation governing unsaturated flow. The analytical solution is obtained by using the Hankel integral transform. For the multilayered case, the combination of a special forward and backward transfer matrix techniques makes the numerical evaluation of the solution very accurate and efficient. The correctness of both formulations is validated by comparison with alternative solutions for two different cases. The results from typical cases are presented to illustrate the influence on the flow field of the cross-anisotropic hydraulic conductivity, the soil heterogeneity and the depth of the source. Copyright © 2004 John Wiley & Sons, Ltd. [source]

Numerical modelling of hydro-mechanical behaviour of collapsible soils

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 12 2008
L. M. Costa
Abstract This paper presents a numerical simulation of the response of unsaturated soil under different loads, with special attention devoted to the phenomenon of pore collapse during wetting. A coupled hydro-mechanical formulation, implemented in the computational code CODE_BRIGHT, is first presented. Then the model is applied to assess material parameter from laboratory tests and to simulate, as boundary value problems, new and appealing in situ tests designed in the Federal University of Pernambuco. The soil considered is a collapsible soil located in the semi-arid region of the northeast of Brazil. Two situations are analysed: (i) the soil is initially at its natural water content and then flooded at some determined stress level, (ii) the soil is previously flooded and then loaded. Comparison between output results of numerical simulations and experimental data shows a very good agreement, which validate to some extent the proposed experimental procedure and the model formulation. Copyright © 2007 John Wiley & Sons, Ltd. [source]

Analysis of heat transfer characteristics of an unsaturated soil bed: a simplified numerical method

INTERNATIONAL JOURNAL OF ENERGY RESEARCH, Issue 15 2001
Gopal B. Reddy
Abstract This paper is a continuation of a study reported in this Journal in February 1999. The paper presents a summary of the two-dimensional macroscopic continuity, momentum and energy equations in a cylindrical co-ordinate system that describe heat and mass transfer through unsaturated soil. The hydrodynamic laws governing flow of water through unsaturated soil are also presented. The explicit numerical procedure and the method to solve the equations are described. Characteristics of the corresponding computer program are also discussed. The results obtained with the current cylindrical governing equations are compared with the previously reported results based upon the Cartesian system of equations. It is observed that the results obtained with cylindrical formulations are in closer agreement with the experimental results. The effects of various heat transfer processes as well as the motion of fluids on heat transfer in a clay bed coupled to a heat pump are discussed. Heat diffusion into the soil by conduction is shown to be predominant through the early stage of heating, while the liquid water motion contributes to heat transfer during the intermediate times and the gas motion is shown to become significant during the last stages of drying. The contribution of the convective transport increases with the temperature and becomes equal to the contribution by conduction at moderately high temperatures. Copyright © 2001 John Wiley & Sons, Ltd. [source]

Reply to comment on ,A novel hysteresis model in unsaturated soil' by A. D. Werner and D. A. Lockington

HYDROLOGICAL PROCESSES, Issue 3 2007
Han-Chen Huang
No abstract is available for this article. [source]

Measurement of the size distribution of water-filled pores at different matric potentials by stray field nuclear magnetic resonance

EUROPEAN JOURNAL OF SOIL SCIENCE, Issue 1 2005
N. R. A. Bird
Summary The water retention characteristic provides the traditional data set for the derivation of a soil's pore-size distribution. However, the technique employed to achieve this requires that assumptions be made about the way pores interconnect. We explore an alternative approach based on stray field nuclear magnetic resonance (STRAFI-NMR) to probe the water-filled pores of both saturated and unsaturated soils, which does not require information relating to pore connectivity. We report the relative size distributions of water-occupied pores in saturated and unsaturated samples of two sets of glass beads of known particle size, two sands, and three soils (a silty loam, a sandy loam and a loamy sand), using measurements of the NMR T1 proton relaxation time of water. The T1 values are linearly related to pore size and consequently measured T1 distributions provide a measure of the pore-size distribution. For both the sands and the glass beads at saturation the T1 distributions are unimodal, and the samples with small particle sizes show a shift to small T1 values indicating smaller voids relative to the samples with larger particles. Different matric potentials were used to reveal how the water-occupied pore-size distribution changes during drainage. These changes are inconsistent with, and demonstrate the inadequacies of, the commonly employed parallel-capillary tube model of a soil pore space. We find that not all pores of the same size drain at the same matric potential. Further, we observe that the T1 distribution is shifted to smaller values beyond the distribution at saturation. This shift is explained by a change in the weighted average of the relaxation rates as the proportion of water in the centre of water-filled pores decreases. This is evidence for the presence of pendular structures resulting from incomplete drainage of pores. For the soils the results are similar except that at saturation the T1 distributions are bimodal or asymmetrical, indicative of inter-aggregate and intra-aggregate pore spaces. We conclude that the NMR method provides a characterization of the water-filled pore space which complements that derived from the water retention characteristic and which can provide insight into the way pore connectivity impacts on drainage. [source]

Relationship between thermal conductivity and water content of soils using numerical modelling

EUROPEAN JOURNAL OF SOIL SCIENCE, Issue 3 2003
P. Cosenza
Summary There is no simple and general relationship between the thermal conductivity of a soil, ,, and its volumetric water content, ,, because the porosity, n, and the thermal conductivity of the solid fraction, ,s, play a major part. Experimental data including measurements of all the variables are scarce. Using a numerical modelling approach, we have shown that the microscopic arrangement of water influences the relation between , and ,. Simulated values for n ranging from 0.4 to 0.6, ,s ranging from 2 to 5 W m,1 K,1 and , from 0.1 to 0.4 can be fitted by a simple linear formula that takes into account n, ,s and ,. The results given by this formula and by the quadratic parallel (QP) model widely used in physical property studies are in satisfactory agreement with published data both for saturated rocks and for unsaturated soils. Consequently, the linear formula and the QP model can be used as practical and efficient tools to investigate the effects of water content and porosity on the thermal conductivity of the soil and hence to optimize the design of thermal in situ techniques for monitoring water content. [source]

The Lisse Effect Revisited

GROUND WATER, Issue 6 2002
Edwin P. Weeks
The Lisse effect is a rarely noted phenomenon occurring when infiltration caused by intense rain seals the surface soil layer to airflow, trapping air in the unsaturated zone. Compression of air by the advancing front results in a pressure increase that produces a water-level rise in an observation well screened below the water table that is several times as large as the distance penetrated by the wetting front. The effect is triggered by intense rains and results in a very rapid water-level rise, followed by a recession lasting a few days. The Lisse effect was first noted and explained by Thal Larsen in 1932 from water-level observations obtained in a shallow well in the village of Lisse, Holland. The original explanation does not account for the increased air pressure pushing up on the bottom of the wetting front. Analysis of the effect of this upward pressure indicates that a negative pressure head at the base of the wetting front, ,f, analogous to that postulated by Green and Ampt (1911) to explain initially rapid infiltration rates into unsaturated soils, is involved in producing the Lisse effect. Analysis of recorded observations of the Lisse effect by Larsen and others indicates that the water-level rise, which typically ranges from 0.10 to 0.55 m, should be only slightly larger than |,f| and that the depth of penetration of the wetting front is no more than several millimeters. [source]

Storage dynamics and streamflow in a catchment with a variable contributing area

HYDROLOGICAL PROCESSES, Issue 16 2010
C. Spence
Abstract Storage heterogeneity effects on runoff generation have been well documented at the hillslope or plot scale. However, diversity across catchments can increase the range of storage conditions. Upscaling the influence of small-scale storage on streamflow across the usually more heterogeneous environment of the catchment has been difficult. The objective of this study was to observe the distribution of storage in a heterogeneous catchment and evaluate its significance and potential influence on streamflow. The study was conducted in the subarctic Canadian Shield: a region with extensive bedrock outcrops, shallow predominantly organic soils, discontinuous permafrost and numerous water bodies. Even when summer runoff was generated from bedrock hillslopes with small storage capacities, intermediary locations with large storage capacities, particularly headwater lakes, prevented water from transmitting to higher order streams. The topographic bounds of the basin thus constituted the maximum potential contributing area to streamflow and rarely the actual area. Topographic basin storage had little relation to basin streamflow, but hydrologically connected storage exhibited a strong hysteretic relationship with streamflow. This relationship defines the form of catchment function such that the basin can be defined by a series of storing and contributing curves comparable with the wetting and drying curves used in relating tension and hydraulic conductivity to water content in unsaturated soils. These curves may prove useful for catchment classification and elucidating predominant hydrological processes. Copyright © 2009 John Wiley & Sons, Ltd and Her Majesty the Queen in right of Canada. [source]

Dynamics of unsaturated soils using various finite element formulations

ACMEG-TS: A constitutive model for unsaturated soils under non-isothermal conditions

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 16 2008
Bertrand François
Abstract This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo-mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto-plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo-hydro-mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non-isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi-mechanism plasticity, whereas water retention characteristics are described using elasto-plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 John Wiley & Sons, Ltd. [source]

A hypoplastic model for mechanical response of unsaturated soils

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 15 2008
David Ma
Abstract A new constitutive model is developed for the mechanical behaviour of unsaturated soils based on the theory of hypoplasticity and the effective stress principle. The governing constitutive relations are presented and their application is demonstrated using several experimental data from the literature. Attention is given to the stiffening effect of suction on the mechanical response of unsaturated soils and the phenomenon of wetting-induced collapse. All model parameters have direct physical interpretation, procedures for their quantification from test data are highlighted. Quantitative predictions of the model are presented for wetting, drying and constant suction tests. Copyright © 2008 John Wiley & Sons, Ltd. [source]

Effective stress concept in unsaturated soils: Clarification and validation of a unified framework

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 7 2008
Mathieu Nuth
Abstract The effective stress principle, conventionally applied in saturated soils, is reviewed for constitutive modelling purposes. The assumptions for the applicability of Terzaghi's single effective stress are recalled and its advantages are inventoried. The possible stress frameworks applicable to unsaturated soil modelling are reassessed in a comparative manner, specifically the Bishop's single effective stress, the independent stress variables approach and the generalized stress framework. The latter considerations lead to the definition of a unified stress context, suitable for modelling soils under different saturation states. In order to qualify the implications brought by the proposed stress framework, several experimental data sets are re-examined in the light of the generalized effective stress. The critical state lines (CSLs) at different saturation states tend to converge remarkably towards a unique saturated line in the deviatoric stress versus mean effective stress plane. The effective stress interpretation is also applied to isotropic paths and compared with conventional net stress conception. The accent is finally laid on a second key feature for constitutive frameworks based on a unified stress, namely the sufficiency of a unique mechanical yield surface besides the unique CSL. Copyright © 2007 John Wiley & Sons, Ltd. [source]

A density-dependent elastoplastic hydro-mechanical model for unsaturated compacted soils

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 11 2007
D. A. Sun
Abstract This paper presents a three-dimensional elastoplastic constitutive model for predicting the hydraulic and mechanical behaviour of unsaturated soils. It is based on experimental results obtained from a series of controlled-suction triaxial tests on unsaturated compacted clay with different initial densities. Hydraulic hysteresis in the water-retention behaviour is modelled as an elastoplastic process, with the elastic part modelled by a series of scanning curves and the elastoplastic part modelled by the main drying and wetting curves. The effect of void ratio on the water-retention behaviour is studied using data obtained from controlled-suction wetting,drying cyclic tests on unsaturated compacted clay with different initial densities. The effect of the degree of saturation on the stress,strain-strength behaviour and the effect of void ratio on the water-retention behaviour are considered in the model, as is the effect of suction on the hydraulic and mechanical behaviour. The initial density dependency of the compacted soil behaviour is modelled by experimental relationships between the initial density and the corresponding yield stress and, thereby, between the initial density and the normal compression line. The model is generalized to three-dimensional stress states by assuming that the shapes of the failure and yield surfaces in the deviatoric stress plane are given by the Matsuoka,Nakai criterion. Model predictions of the stress,strain and water-retention behaviour are compared with those obtained from triaxial tests with different initial densities under isotropic compression, triaxial compression and triaxial extension, with or without variation in suction. The comparisons indicate that the model accurately predicts the hydraulic and mechanical behaviour of unsaturated compacted soils with different initial densities using the same material constant. Copyright © 2006 John Wiley & Sons, Ltd. [source]

A unified bounding surface plasticity model for unsaturated soils

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2006
A.R. Russell
Abstract A unified constitutive model for unsaturated soils is presented in a critical state framework using the concepts of effective stress and bounding surface plasticity theory. Consideration is given to the effects of unsaturation and particle crushing in the definition of the critical state. A simple isotropic elastic rule is adopted. A loading surface and a bounding surface of the same shape are defined using simple and versatile functions. The bounding surface and elastic rules lead to the existence of a limiting isotropic compression line, towards which the stress trajectories of all isotropic compression load paths approach. A non-associated flow rule of the same general form is assumed for all soil types. Isotropic hardening/softening occurs due to changes in plastic volumetric strains as well as suction for some unsaturated soils, enabling the phenomenon of volumetric collapse upon wetting to be accounted for. The model is used to simulate the stress,strain behaviour observed in unsaturated speswhite kaolin subjected to three triaxial test load paths. The fit between simulation and experiment is improved compared to that of other constitutive models developed using conventional Cam-Clay-based plasticity theory and calibrated using the same set of data. Also, the model is used to simulate to a high degree of accuracy the stress,strain behaviour observed in unsaturated Kurnell sand subjected to two triaxial test load paths and the oedometric compression load path. For oedometric compression theoretical simulations indicate that the suction was not sufficiently large to cause samples to separate from the confining ring. Copyright © 2005 John Wiley & Sons, Ltd. [source]

Finite element formulation and algorithms for unsaturated soils.

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 9 2003
Part I: Theory
Abstract This paper presents a complete finite-element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress,strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time-stepping scheme which automatically adjusts the time-step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non-linearity caused by suction-dependent plastic yielding, suction-dependent degree of saturation, and saturation-dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress,strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd. [source]

Solid,liquid,air coupling in multiphase porous media

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2003
Lyesse Laloui
Abstract This paper addresses various issues concerning the modelling of solid,liquid,air coupling in multiphase porous media with an application to unsaturated soils. General considerations based on thermodynamics permit the derivation and discussion of the general form of field equations; two cases are considered: a three phase porous material with solid, liquid and gas, and a two phase porous material with solid, liquid and empty space. Emphasis is placed on the presentation of differences in the formulation and on the role of the gas phase. The finite element method is used for the discrete approximation of the partial differential equations governing the problem. The two formulations are then analysed with respect to a documented drainage experiment carried out by the authors. The merits and shortcomings of the two approaches are shown. Copyright © 2003 John Wiley & Sons, Ltd. [source]