Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid discrete-continuum approach to model hydro-mechanical behaviour of soil during desiccation
69
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Desiccation cracking in clayey soils occurs when they lose moisture, leading to an increase in their compressibility and hydraulic conductivity and hence significant reduction of soil strength. The prediction of desiccation cracking in soils is challenging due to the lack of insights into the complex coupled hydro-mechanical process at the grain scale. In this paper, a new hybrid discrete-continuum numerical framework, capable of capturing hydro-mechanical behaviour of soil at both grain and macro scales, is proposed for predicting desiccation cracking in clayey soil. In this framework, a soil layer is represented by an assembly of DEM particles, each occupies an equivalent continuum space and carries physical properties governing unsaturated flow. These particles move freely in the computational space following the discrete element method (DEM), while their contact network and the continuum mixture theory are used to model the unsaturated flow. The dependence of particle-to-particle contact behaviour on water content is represented by a cohesive-frictional contact model, whose material properties are governed by the water content. In parallel with the theoretical development is a series of experiments on 3D soil desiccation cracking to determine essential properties and provide data for the validation of mechanical and physical behaviour. Very good agreement in both physical behaviour (e.g. evolution of water content) and mechanical behaviour (e.g. occurrence and development of cracks, and distribution of compressive and tensile strains) demonstrates that the proposed framework is capable of capturing the hydro-mechanical behaviour of soil during desiccation. The capability of the proposed framework facilitates numerical experiments for insights into the hydro-mechanical behaviour of unsaturated soils that have not been possible before.
rate research
Read More
This article focuses on liquefaction of saturated granular soils, triggered by earthquakes. Liquefaction is definedhere as the transition from a rigid state, in which the granular soil layer supports structures placed on its surface, toa fluidlike state, in which structures placed initially on the surface sink to their isostatic depth within the granularlayer.We suggest a simple theoretical model for soil liquefaction and show that buoyancy caused by the presence ofwater inside a granular medium has a dramatic influence on the stability of an intruder resting at the surface of themedium.We confirm this hypothesis by comparison with laboratory experiments and discrete-element numericalsimulations. The external excitation representing ground motion during earthquakes is simulated via horizontalsinusoidal oscillations of controlled frequency and amplitude. In the experiments, we use particles only slightlydenser than water, which as predicted theoretically increases the effect of liquefaction and allows clear depth-of-sinkingmeasurements. In the simulations, a micromechanical model simulates grains using molecular dynamicswith friction between neighbors. The effect of the fluid is captured by taking into account buoyancy effects onthe grains when they are immersed. We show that the motion of an intruder inside a granular medium is mainlydependent on the peak acceleration of the ground motion and establish a phase diagram for the conditions underwhich liquefaction happens, depending on the soil bulk density, friction properties, presence of water, and peak acceleration of the imposed large-scale soil vibrations.We establish that in liquefaction conditions, most cases relaxtoward an equilibrium position following an exponential in time.We also show that the equilibrium position itself,for most liquefaction regimes, corresponds to the isostatic equilibrium of the intruder inside a medium of effectivedensity. The characteristic time to relaxation is shown to be essentially a function of the peak ground velocity.
Large scale modelling of fluid flow coupled with solid failure in geothermal reservoirs or hydrocarbon extraction from reservoir rocks usually involves behaviours at two scales: lower scale of the inelastic localization zone, and larger scale of the bulk continuum where elastic behaviour can be reasonably assumed. The hydraulic conductivities corresponding to the mechanical properties at these two scales are different. In the bulk elastic host rock, the hydraulic conductivity does not vary much with the deformation, while it significantly changes in the lower scale of the localization zone due to inelastic deformation. Increase of permeability due to fracture and/or dilation, or reduction of permeability due to material compaction can take place inside this zone. The challenge is to predict the evolution of hydraulic conductivities coupled with the mechanical behaviour of the material in all stages of the deformation process. In the early stage of diffuse deformation, the permeability of the material can be reasonably assumed to be homogenous over the whole Representative Volume Element (RVE) However, localized failure results in distinctly different conductivities in different parts of the RVE. This paper establishes a general framework and corresponding field equations to describe the hydro-mechanical coupling in both diffuse and localized stages of deformation in rocks. In particular, embedding the lower scale hydro-mechanical behaviour of the localization zone inside an elastic bulk, together with their corresponding effective sizes, helps effectively deal with scaling issues in large-scale modelling. Preliminary results are presented which demonstrate the promising features of this new approach.
Soil has been recognized as an indirect driver of global warming by regulating atmospheric greenhouse gases. However, in view of the higher heat capacity and CO2 concentration in soil than those in atmosphere, the direct contributions of soil to greenhouse effect may be non-ignorable. Through field manipulation of CO2 concentration both in soil and atmosphere, we demonstrated that the soil-retained heat and its slow transmission process within soil may cause slower heat leaking from the earth. Furthermore, soil air temperature was non-linearly affected by soil CO2 concentration with the highest value under 7500 ppm CO2. This study indicates that the soil and soil CO2, together with atmospheric CO2, play indispensable roles in fueling the greenhouse effect. We proposed that anthropogenic changes in soils should be focused in understanding drivers of the globe warming.
This work concerns the continuum basis and numerical formulation for deformable materials with viscous dissipative mechanisms. We derive a viscohyperelastic modeling framework based on fundamental thermomechanical principles. Since most large deformation problems exhibit the isochoric property, our modeling work is constructed based on the Gibbs free energy in order to develop a continuum theory using the pressure-primitive variables, which is known to be well-behaved in the incompressible limit. With a general theory presented, we focus on a family of free energies that leads to the so-called finite deformation linear model. Our derivation elucidates the origin of the evolution equations of that model, which was originally proposed heuristically. In our derivation, the thermodynamic inconsistency is clarified and rectified. We then discuss the relaxation property of the non-equilibrium stress in the thermodynamic equilibrium limit and its implication on the form of free energy. A modified version of the identical polymer chain model is then proposed, with a special case being the model proposed by G. Holzapfel and J. Simo. Based on the consistent modeling framework, a provably energy stable numerical scheme is constructed for incompressible viscohyperelasticity using inf-sup stable elements. In particular, we adopt a suite of smooth generalization of the Taylor-Hood element based on Non-Uniform Rational B-Splines (NURBS) for spatial discretization. The temporal discretization is performed via the generalized-alpha scheme. We present a suite of numerical results to corroborate the proposed numerical properties, including the nonlinear stability, robustness under large deformation, and the stress accuracy resolved by the higher-order elements.
A minimalist stochastic model of primary soil salinity is proposed, in which the rate of soil salinization is determined by the balance between dry and wet salt deposition and the intermittent leaching events caused by rainfall events. The long term probability density functions of salt mass and concentration are found by reducing the coupled soil moisture and salt mass balance equation to a single stochastic differential equation driven by multiplicative Poisson noise. The novel analytical solutions provide insight on the interplay of the main soil, plant and climate parameters responsible for long-term soil salinization. In particular, they show the existence of two distinct regimes, one where the mean salt mass remains nearly constant (or decreases) with increasing rainfall frequency, and another where mean salt content increases markedly with increasing rainfall frequency. As a result, relatively small reductions of rainfall in drier climates may entail dramatic shifts in long-term soil salinization trends, with significant consequences e.g. for climate change impacts on rain-fed agriculture.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
