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.

Stiffness and strength of tridimensional periodic lattices

This paper presents a method for the linear analysis of the stiffness and strength of open and closed cell lattices with arbitrary topology. The method hinges on a multiscale approach that separates the analysis of the lattice in two scales. At the macroscopic level, the lattice is considered as a uniform material; at the microscopic scale, on the other hand, the cell microstructure is modelled in detail by means of an in-house finite element solver. The method allows determine the macroscopic stiffness, the internal forces in the edges and walls of the lattice, as well as the global periodic buckling loads, along with their buckling modes. Four cube-based lattices and nine cell topologies derived by Archimedean polyhedra are studied. Several of them are characterized here for the first time with a particular attention on the role that the cell wall plays on the stiffness and strength properties. The method, automated in a computational routine, has been used to develop material property charts that help to gain insight into the performance of the lattices under investigation. © 2012 Elsevier B.V.

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.

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.

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.

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.

Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Carbonating magnesia for soil stabilization

This paper investigates the potential for carbonating reactive magnesia (MgO) to serve as a more sustainable soil stabilization method by providing rapid and significant strength development of the stabilized soil through absorbing substantial quantities of CO2. Gaseous CO2 was forced through laboratory-prepared reactive MgO-treated soil samples in a triaxial cell set-up, and their resulting mechanical and microstructural properties were investigated using unconfined compressive strength, X-ray diffraction, and scanning electron microscopy. The results showed that adequately carbonated MgO-treated soils could, in a few hours, reach a similar strength range to corresponding 28 day Portland cement (PC)-stabilized soils. Hydrated magnesium carbonates, namely nesquehonite and hydromagnesite-dypingite, were the main products of the carbonated MgO in the soil, and were responsible for the significant strength development.

Vibration from railways: Can we achieve better than +/-10dB prediction accuracy?

This paper examines the sources of uncertainly in models used to predict vibration from underground railways. It will become clear from this presentation that by varying parameters by a small amount, consistent with uncertainties in measured data, the predicted vibration levels vary significantly, often by more than 10dB. This error cannot be forecast. Small changes made to soil parameters (Compressive and Shear Wave velocities and density), to slab bending stiffness and mass and to the measurement position give rise to changes in vibration levels of more than lOdB. So if 10dB prediction error results from small uncertainties in soil parameters and measurement position it cannot be sensible to rely on prediction models for accuracy better than 10dB. The presentation will demonstrate in real time the use of the new - and freely-available - PiP software for calculating vibration from railway tunnels in real time.