Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.

Effect of soil conditions on uplift of underground structures in liquefied soil

In an earthquake, underground structures located in liquefiable soil deposits are susceptible to floatation following an earthquake event due to their lower unit weight relative to the surrounding saturated soil. Such uplift response of the buoyant structure is influenced by the soil it is buried in. In the case of a liquefiable soil deposit, the soil can lose its shear strength significantly in the event of an earthquake. If the soil liquefies fully, the buoyant structure can float towards the soil surface. However, a partly liquefied soil deposit retains some of its initial shear strength and resists the uplift. This paper discusses the different soil conditions and their influence on the uplift response of buoyant structures. © 2012 World Scientific Publishing Company.

An efficient model for the dynamic behaviour of a single pile in viscoelastic soil

This paper presents a novel, three-dimensional, single-pile model, formulated in the wavenumber domain and adapted to account for boundary conditions using the superposition of loading cases. The pile is modelled as a column in axial vibration, and a Euler-Bernoulli beam in lateral vibration. The surrounding soil is treated as a viscoelastic continuum. The response of the pile is presented in terms of the stiffness and damping coefficients, and also the magnitude and phase of the pile-head frequency-response function. Comparison with existing models shows that excellent agreement is observed between this model, a boundary-element formulation, and an elastic-continuum-type formulation. This three-dimensional model has an accuracy equivalent to a 3D boundary-element model, and a runtime similar to a 2D plane-strain analytical model. Analysis of the response of the single pile illustrates a difference in axial and lateral vibration behaviour; the displacement along the pile is relatively invariant under axial loads, but in lateral vibration the pile exhibits localised deformations. This implies that a plane-strain assumption is valid for axial loadings and only at higher frequencies for lateral loadings. © 2013 Elsevier Ltd.

Estimation of soil properties and deformations in staged constructions based on mcmc method

This paper presents a Bayesian probabilistic framework to assess soil properties and model uncertainty to better predict excavation-induced deformations using field deformation data. The potential correlations between deformations at different depths are accounted for in the likelihood function needed in the Bayesian approach. The proposed approach also accounts for inclinometer measurement errors. The posterior statistics of the unknown soil properties and the model parameters are computed using the Delayed Rejection (DR) method and the Adaptive Metropolis (AM) method. As an application, the proposed framework is used to assess the unknown soil properties of multiple soil layers using deformation data at different locations and for incremental excavation stages. The developed approach can be used for the design of optimal revisions for supported excavation systems. © 2010 ASCE.

Bayesian updating of a unified soil compression model

Some amount of differential settlement occurs even in the most uniform soil deposit, but it is extremely difficult to estimate because of the natural heterogeneity of the soil. The compression response of the soil and its variability must be characterised in order to estimate the probability of the differential settlement exceeding a certain threshold value. The work presented in this paper introduces a probabilistic framework to address this issue in a rigorous manner, while preserving the format of a typical geotechnical settlement analysis. In order to avoid dealing with different approaches for each category of soil, a simplified unified compression model is used to characterise the nonlinear compression behavior of soils of varying gradation through a single constitutive law. The Bayesian updating rule is used to incorporate information from three different laboratory datasets in the computation of the statistics (estimates of the means and covariance matrix) of the compression model parameters, as well as of the uncertainty inherent in the model.

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.

Probabilistic soil identification based on cone penetration tests

In geotechnical engineering, soil classification is an essential component in the design process. Field methods such as the cone penetration test (CPT) can be used as less expensive and faster alternatives to sample retrieval and testing. Unfortunately, current soil classification charts based on CPT data and laboratory measurements are too generic, and may not provide an accurate prediction of the soil type. A probabilistic approach is proposed here to update and modify soil identification charts based on site-specific CPT data. The probability that a soil is correctly classified is also estimated. The updated identification chart can be used for a more accurate prediction of the classification of the soil, and can account for prior information available before conducting the tests, site-specific data, and measurement errors. As an illustration, the proposed approach is implemented using CPT data from the Treporti Test Site (TTS) near Venice (Italy) and the National Geotechnical Experimentation Sites (NGES) at Texas A&M University. The applicability of the site-specific chart for other sites in Venice Lagoon is assessed using data from the Malamocco test site, approximately 20 km from TTS.

Wave propagation due to sinusoidal excitation in nonlinear one-dimensional soil column

The objective of the author's on-going research is to explore the feasibility of determining reliable in situ curves of shear modulus as a function of strain using the dynamic test. The purpose of this paper is limited to investigating what material stiffness is measured from a dynamic test, focusing on the harmonic excitation test. A one-dimensional discrete model with nonlinear material properties is used for this purpose. When a sinusoidal load is applied, the cross-correlation of signals from different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loops under steady-state conditions. The variables that contributed to changing the average slope of the stress-strain loop also influence the estimate of the wave velocity from cross-correlation. Copyright ASCE 2007.

Seismic triggering of submarine slides in soft cohesive soil deposits

The geological profile of many submerged slopes on the continental shelf consists of normally to lightly overconsolidated clays with depths ranging from a few meters to hundreds of meters. For these soils, earthquake loading can generate significant excess pore water pressures at depth, which can bring the slope to a state of instability during the event or at a later time as a result of pore pressure redistribution within the soil profile. Seismic triggering mechanisms of landslide initiation for these soils are analyzed with the use of a new simplified model for clays which predicts realistic variations of the stress-strain-strength relationships as well as pore pressure generation during dynamic loading in simple shear. The proposed model is implemented in a finite element program to analyze the seismic response of submarine slopes. These analyses provide an assessment of the critical depth and estimated displacements of the mobilized materials and thus are important components for the estimation of submarine landslide-induced tsunamis. © 2003 Elsevier B.V. All rights reserved.

Robust stabilization of a nonlinear cement mill model

Plugging is well known to be a major cause of instability in industrial cement mills. A simple nonlinear model able to simulate the plugging phenomenon is presented. It is shown how a nonlinear robust controller can be designed in order to fully prevent the mill from plugging.

Preliminary laboratory-scale model auger installation and testing of carbonated soil-MgO columns

This paper presents details of the installation and performance of carbonated soil-MgO columns using a laboratory-scale model auger setup. MgO grout was mixed with the soil using the auger and the columns were then carbonated with gaseous CO2 introduced in two different ways: one using auger mixing and the other through a perforated plastic tube system inserted into the treated column. The performance of the columns in terms of unconfined compressive strength (UCS), stiffness, strain at failure and microstructure (using X-ray diffraction and scanning electron microscopy) showed that the soil-MgO columns were carbonated very quickly (in under 1 h) and yielded relatively high strength values, of 2.4-9.4 MPa, which on average were five times that of corresponding 28-day ambient cured uncarbonated columns. This confirmed, together with observations of dense microstructure and hydrated magnesium carbonates, that a good degree of carbonation had taken place. The results also showed that the carbonation method and period have a significant effect on the resulting performance, with the carbonation through the perforated pipe producing the best results. Copyright © 2013 by ASTM International.

Cement bentonite cutoff walls for polluted sites

Cement-bentonite (CB) cutoff walls have long been used to control ground water flow and contaminant migration at polluted sites. Hydraulic conductivity and unconfined compressive strength are two short-term properties often used by industry and owners in CB specification and are important parameters discussed in this paper. For polluted sites, long-term compatibility is also an important issue. These properties are coupled to a number of external factors including the mix design, construction sequence, presence/absence of contaminants at the site. Additional short-term properties for engineering assessment include the stressstrain characteristics in both drained and undrained shear in both with and without confinement as well as one-dimensional consolidation properties. Long term CB properties are affected by aging, reaction chemistry, drying, in situ stress state, and interaction with the polluted environment. © 2013 Taylor & Francis Group.