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.

Cement-fly ash stabilisation/solidification of contaminated soil: Performance properties and initiation of operating envelopes

This study was aimed at evaluating the mechanical and pH-dependent leaching performance of a mixed contaminated soil treated with a mixture of Portland cement (CEMI) and pulverised fuel ash (PFA). It also sought to develop operating envelopes, which define the range(s) of operating variables that result in acceptable performance. A real site soil with low contaminant concentrations, spiked with 3000mg/kg each of Cd, Cu, Pb, Ni and Zn, and 10,000mg/kg of diesel, was treated with one part CEMI and four parts PFA (CEMI:PFA=1:4) using different binder and water contents. The performance was assessed over time using unconfined compressive strength (UCS), hydraulic conductivity, acid neutralisation capacity (ANC) and pH-dependent leachability of contaminants. With binder dosages ranging from 5% to 20% and water contents ranging from 14% to 21% dry weight, the 28-day UCS was up to 500kPa and hydraulic conductivity was around 10-8m/s. With leachant pH extremes of 7.2 and 0.85, leachability of the contaminants was in the range: 0.02-3500mg/kg for Cd, 0.35-1550mg/kg for Cu, 0.03-92mg/kg for Pb, 0.01-3300mg/kg for Ni, 0.02-4010mg/kg for Zn, and 7-4884mg/kg for total petroleum hydrocarbons (TPHs), over time. Design charts were produced from the results of the study, which show the water and/or binder proportions that could be used to achieve relevant performance criteria. The charts would be useful for the scale-up and design of stabilisation/solidification (S/S) treatment of similar soil types impacted with the same types of contaminants. © 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.

Numerical modeling of deep soil mixing excavation support

The finite element method (FEM) is growing in popularity over the pressure diagram/hand calculation method for analysis of excavation systems in general and deep soil mixing excavations in particular. In this paper, a finite element analysis is used to study the behavior of a deep mixed excavation. Through the use of Plaxis (a FEM software program), the construction sequence is simulated by following the various construction phases allowing for deflections due to strut or anchor installation to be predicted. The numerical model used in this study simulates the soil cement columns as a continuous wall matching the bending stiffness of the actual wall. Input parameters based on laboratory tests and modeling assumptions are discussed. An example of the approach is illustrated using the Islais Creek Transport/Storage Project in San Francisco, California. Copyright ASCE 2006.

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.

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.

Relating monolithic and granular leaching from contaminated soil treated with different cementitious binders.

This work employed a clayey, silty, sandy gravel contaminated with a mixture of metals (Cd, Cu, Pb, Ni and Zn) and diesel. The contaminated soil was treated with 5 and 10% dosages of different cementitious binders. The binders include Portland cement, cement-fly ash, cement-slag and lime-slag mixtures. Monolithic leaching from the treated soils was evaluated over a 64-day period alongside granular leachability of 49- and 84-day old samples. Surface wash-off was the predominant leaching mechanism for monolithic samples. In this condition, with data from different binders and curing ages combined, granular leachability as a function of monolithic leaching generally followed degrees 4 and 6 polynomial functions. The only exception was for Cu, which followed the multistage dose-response model. The relationship between both leaching tests varied with the type of metal, curing age/residence time of monolithic samples in the leachant, and binder formulation. The results provide useful design information on the relationship between leachability of metals from monolithic forms of S/S treated soils and the ultimate leachability in the eventual breakdown of the stabilized/solidified soil.