Bruno Chareyre

On the capillary stress tensor in wet granular materials

This paper presents a micromechanical study of unsaturated granular media in the pendular regime, based on numerical experiments using the discrete element method, compared with a microstructural elastoplastic model. Water effects are taken into account by adding capillary menisci at contacts and their consequences in terms of force and water volume are studied. Simulations of triaxial compression tests are used to investigate both macro and micro-effects of a partial saturation. The results provided by the two methods appear to be in good agreement, reproducing the major trends of a partially saturated granular assembly, such as the increase in the shear strength and the hardening with suction. Moreover, a capillary stress tensor is exhibited from capillary forces by using homogenization techniques. Both macroscopic and microscopic considerations emphasize an induced anisotropy of the capillary stress tensor in relation with the pore fluid distribution inside the material. Insofar as the tensorial nature of this fluid fabric implies shear effects on the solid phase associated with suction, a comparison has been made with the standard equivalent pore pressure assumption. It is shown that water effects induce microstructural phenomena that cannot be considered at the macro level, particularly when dealing with material history. Thus, the study points out that unsaturated soil stress definitions should include, besides the macroscopic stresses such as the total stress, the microscopic interparticle stresses such as the ones resulting from capillary forces, in order to interpret more precisely the implications of the pore fluid on the mechanical behaviour of granular materials. Copyright

Pore-scale modeling of fluid-particles interaction and emerging poromechanical effects

SUMMARY A micro-hydromechanical model for granular materials is presented. It combines the discrete element method for the modeling of the solid phase and a pore-scale finite volume formulation for the flow of an incompressible pore fluid. The coupling equations are derived and contrasted against the equations of conventional poroelasticity. An analogy is found between the discrete element method pore-scale finite volume coupling and Biots theory in the limit case of incompressible phases. The simulation of an oedometer test validates the coupling scheme and demonstrates the ability of the model to capture strong poromechanical effects. A detailed analysis of microscale strain and stress confirms the analogy with poroelasticity. An immersed deposition problem is finally simulated and shows the potential of the method to handle phase transitions. Copyright

Pore-Scale Modeling of Viscous Flow and Induced Forces in Dense Sphere Packings

We propose a method for effectively upscaling incompressible viscous flow in large random polydispersed sphere packings: the emphasis of this method is on the determination of the forces applied on the solid particles by the fluid. Pore bodies and their connections are defined locally through a regular Delaunay triangulation of the packings. Viscous flow equations are upscaled at the pore level, and approximated with a finite volume numerical scheme. We compare numerical simulations of the proposed method to detailed finite element simulations of the Stokes equations for assemblies of 8–200 spheres. A good agreement is found both in terms of forces exerted on the solid particles and effective permeability coefficients.

A minimal coupled fluid-discrete element model for bedload transport

A minimal Lagragian two-phase model to study turbulent bedload transport focusing on the granular phase is presented, and validated with experiments. The model intends to describe bedload transport of massive particles in fully rough flows at relatively low Shields numbers, for which no suspension occurs. A discrete element method for the granular phase is coupled with a one dimensional volume-averaged two-phase momentum equation for the fluid phase. The coupling between the discrete granular phase and the continuous fluid phase is discussed, and a consistent averaging formulation adapted to bedload transport is introduced. An original simple discrete random walk model is proposed to account for the fluid velocity fluctuations. The model is compared with experiments considering both classical sediment transport rate as a function of the Shields number, and depth profiles of solid velocity, volume fraction, and transport rate density, from existing bedload transport experiments in inclined flume. The results successfully reproduce the classical 3/2 power law, and more importantly describe well the depth profiles of the granular phase, showing that the model is able to reproduce the particle scale mechanisms. From a sensitivity analysis, it is shown that the fluctuation model allows to reproduce a realistic critical Shields number, and that the influence of the granular parameters on the macroscopic results are weak. Nevertheless, the analysis of the corresponding depth profiles reveals an evolution of the depth structure of the granular phase with varying restitution and friction coefficients, which denotes the non-trivial underlying physical mechanisms.

Pore-scale simulations of drainage in granular materials: Finite size effects and the representative elementary volume

Abstract A pore-scale model is introduced for two-phase flow in dense packings of polydisperse spheres. The model is developed as a component of a more general hydromechanical coupling framework based on the discrete element method, which will be elaborated in future papers and will apply to various processes of interest in soil science, in geomechanics and in oil and gas production. Here the emphasis is on the generation of a network of pores mapping the void space between spherical grains, and the definition of local criteria governing the primary drainage process. The pore space is decomposed by Regular Triangulation, from which a set of pores connected by throats are identified. A local entry capillary pressure is evaluated for each throat, based on the balance of capillary pressure and surface tension at equilibrium. The model reflects the possible entrapment of disconnected patches of the receding wetting phase. It is validated by a comparison with drainage experiments. In the last part of the paper, a series of simulations are reported to illustrate size and boundary effects, key questions when studying small samples made of spherical particles be it in simulations or experiments. Repeated tests on samples of different sizes give evolution of water content which are not only scattered but also strongly biased for small sample sizes. More than 20,000 spheres are needed to reduce the bias on saturation below 0.02. Additional statistics are generated by subsampling a large sample of 64,000 spheres. They suggest that the minimal sampling volume for evaluating saturation is one hundred times greater that the sampling volume needed for measuring porosity with the same accuracy. This requirement in terms of sample size induces a need for efficient computer codes. The method described herein has a low algorithmic complexity in order to satisfy this requirement. It will be well suited to further developments toward coupled flow-deformation problems in which evolution of the microstructure require frequent updates of the pore network.

The Effects of Swelling and Porosity Change on Capillarity: DEM Coupled with a Pore-Unit Assembly Method.

In this study, a grain-scale modelling technique has been developed to generate the capillary pressure–saturation curves for swelling granular materials. This model employs only basic granular properties such as particles size distribution, porosity, and the amount of absorbed water for swelling materials. Using this model, both drainage and imbibition curves are directly obtained by pore-scale simulations of fluid invasion. This allows us to produce capillary pressure–saturation curves for a large number of different packings of granular materials with varying porosity and/or amount of absorbed water. The algorithm is based on combining the Discrete Element Method for generating different particle packings with a pore-unit assembly approach. The pore space is extracted using a regular triangulation, with the centres of four neighbouring particles forming a tetrahedron. The pore space within each tetrahedron is referred to as a pore unit. Thus, the pore space of a particle packing is represented by an assembly of pore units for which we construct drainage and imbibition capillary pressure–saturation curves. A case study on Hostun sand is conducted to test the model against experimental data from literature and to investigate the required minimum number of particles to have a Representative Elementary Volume. Then, the capillary pressure–saturation curves are constructed for Absorbent Gelling Material particles, for different combinations of porosity values and amounts of absorbed water. Each combination yields a different configuration of pore units, and thus distinctly different capillary pressure–saturation curves. All these curves are shown to collapse into one curve for drainage and one curve for imbibition when we normalize capillary pressure and saturation values. We have developed a formula for the Van Genuchten parameter

DEM Modeling of a Flexible Barrier Impacted by a Dry Granular Flow

\alpha

Role of blockages in particle transport through homogeneous granular assemblies

α (which is related to the inverse of the entry pressure) as a function of porosity and the amount of absorbed water.