Microscopic origins of shear stress in dense fluid-grain mixtures

A numerical model is used to simulate rheometer experiments at constant normal stress on dense suspensions of spheres. The complete model includes sphere-sphere contacts using a soft contact approach, short range hydrodynamic interactions defined by frame-invariant expressions of forces and torques in the lubrication approximation, and drag forces resulting from the poromechanical coupling computed with the DEM-PFV technique. Series of simulations in which some of the coupling terms are neglected highlight the role of the poromechanical coupling in the transient regimes. They also reveal that the shear component of the lubrication forces, though frequently neglected in the literature, has a dominant effect in the volume changes. On the other hand, the effects of lubrication torques are much less significant. The bulk shear stress is decomposed into contact stress and hydrodynamic stress terms whose dependency on a dimensionless shear rate - the so called viscous number $I_v$ - are examined. Both contributions are increasing functions of $I_v$, contacts contribution dominates at low viscous number ($I_v$ 0.15, consistently with a phenomenological law infered by other authors. Statistics of microstructural variables highlight a complex interplay between solid contacts and hydrodynamic interactions. In contrast with a popular idea, the results suggest that lubrication may not necessarily reduce the contribution of contact forces to the bulk shear stress. The proposed model is general and applies directly to sheared immersed granular media in which pore pressure feedback plays a key role (triggering of avalanches, liquefaction).Comment: to appear in Granular Matte

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

A micro-hydromechanical model for granular materials is presented. It combines the discrete element method (DEM) for the modeling of the solid phase and a pore-scale finite volume (PFV) 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 DEM-PFV coupling and Biot's 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.Comment: accepted in Int. Journal for Numerical and Analytical Methods in Geomechanic

Discrete modelling of capillary mechanisms in multi-phase granular media

A numerical study of multi-phase granular materials based upon micro-mechanical modelling is proposed. Discrete element simulations are used to investigate capillary induced effects on the friction properties of a granular assembly in the pendular regime. Capillary forces are described at the local scale through the Young-Laplace equation and are superimposed to the standard dry particle interaction usually well simulated through an elastic-plastic relationship. Both effects of the pressure difference between liquid and gas phases and of the surface tension at the interface are integrated into the interaction model. Hydraulic hysteresis is accounted for based on the possible mechanism of formation and breakage of capillary menisci at contacts. In order to upscale the interparticular model, triaxial loading paths are simulated on a granular assembly and the results interpreted through the Mohr-Coulomb criterion. The micro-mechanical approach is validated with a capillary cohesion induced at the macroscopic scale. It is shown that interparticular menisci contribute to the soil resistance by increasing normal forces at contacts. In addition, more than the capillary pressure level or the degree of saturation, our findings highlight the importance of the density number of liquid bonds on the overall behaviour of the material

Localized fluidization in granular materials: Theoretical and numerical study

We present analytical and numerical results on localized fluidization within a granular layer subjected to a local injection of fluid. As the injection rate increases the three different regimes previously reported in the literature are recovered: homogeneous expansion of the bed, fluidized cavity in which fluidization starts developing above the injection area, and finally the chimney of fluidized grains when the fluidization zone reaches the free surface. The analytical approach is at the continuum scale, based on Darcy's law and Therzaghi's effective stress principle. It provides a good description of the phenomenon as long as the porosity of the granular assembly remains relatively homogeneous, i.e. for small injection rates. The numerical approach is at the particle scale based on the coupled DEM-PFV method. It tackles the more heterogeneous situations which occur at larger injection rates. The results from both methods are in qualitative agreement with data published independently. A more quantitative agreement is achieved by the numerical model. A direct link is evidenced between the occurrence of the different regimes of fluidization and the injection aperture. While narrow apertures let the three different regimes be distinguished clearly, larger apertures tend to produce a single homogeneous fluidization regime. In the former case, it is found that the transition between the cavity regime and the chimney regime for an increasing injection rate coincides with a peak in the evolution of inlet pressure. Finally, the occurrence of the different regimes is defined in terms of the normalized flux and aperture

Saint-Thomé – Le Bouzil

Le Bouzil est un des très rares grands sites de plein air de la région Rhône-Alpes attribuable au Périgordien. L’extension de la fouille, en 1991, a confirmé l’existence de trois niveaux d’occupation, marqués par des amas de débitage, un foyer (niveau 1) et par des structures de blocs probablement liées à un habitat (niveau 3). La zone décapée est maintenant protégée par un hangar métallique de 28 m2. Fig. 1 – Plan de la fouille en 1991 Dessin : équipe de fouille. Fig. 2 – Industrie lithique..

Influence of homocysteine on the physical structure and molecular mobility of elastin network in cultured arteries

The thermal and dielectric properties of the elastin network were investigated in arteries cultured with physiological and pathological concentrations of homocysteine, an aminoacid responsible of histological impairments in human arteries. The glass transition of this amorphous protein was investigated by Differential Scanning Calorimetry (DSC). To explore the molecular dynamics of the elastin network in the nanometer range, we used Thermally Stimulated Currents (TSC), a dielectric technique running at low frequency and measuring the dipolar reorientations in proteins subjected to a static electrical field. Combining TSC and DSC experiments with determination of the activation parameters of relaxation times reveals the molecular mobility of the proteins. The major differences in the relaxation behavior of elastin between arteries cultured with physiological and pathological concentrations of homocysteine are discussed

Changes in the physical structure and chain dynamics of elastin network in homocysteine-cultured arteries

The thermal and dielectric properties of the elastin network were investigated in arteries cultured with physiological and pathological concentrations of homocysteine, an aminoacid responsible of histological impairments in human arteries. The physical structure of this amorphous protein was investigated by differential scanning calorimetry (DSC). To explore the molecular dynamics of the elastin network in the nanometer range, we used thermally stimulated currents (TSC), a dielectric technique running at low frequency, and measuring the dipolar reorientations in proteins subjected to a static electrical field. Combining DSC and TSC experiments reveals the molecular mobility of the proteins, both in the glassy state and in the liquid state. Significant differences are evidenced in the physical structure and relaxation behavior of elastin network in cultured arteries (physiological and pathological concentrations of homocysteine) and discussed

Numerical simulations of dense suspensions rheology using a dem-fluid coupled model

The understanding of dense suspensions rheology is of great practical interest for both industrial and geophysical applications and has led to a large amount of publications over the past decades. This problem is especially difficult as it is a two-phase media in which particle-particle interactions as well as fluid-particle interactions are significant. In this contribution, the plane shear flow of a dense fluid-grain mixture is studied using the DEM-PFV coupled model. We further improve the original model: including the deviatoric part of the stress tensor on the basis of the lubrication theory, and extending the solver to periodic boundary conditions. Simulations of a granular media saturated by an incompressible fluid and subjected to a plane shear at imposed vertical stress are presented. The shear stress is decomposed in different contributions which can be examined separately: contact forces, lubrication forces and drag forces associated to the poromechanical couplings

From elasto-plasticity to visco-elasto-plasticity for saturated granular materials

A recent extension of the discrete element method is reported for the simulation of dense mixtures of non-colloidal particles and viscous fluids in the non-inertial regime. The numerical model includes sphere-sphere contacts using a soft contact ap- proach [2], short range hydrodynamic interactions defined by frame-invariant expressions of forces and torques in the lubrication approximation, and drag forces resulting from the poromechanical coupling computed with the DEM-PFV technique [3]. The proposed model is general and applies directly to sheared satured granular media in which pore pressure feedback plays a key role. A partitioned solver makes the algorithm trivially parallel, which enables the coupled problems to be solved with nearly the same wall-clock time as uncoupled dry materials simulations. The shear stress in a dense suspension is analyzed, and decomposed into contact stress and hydrodynamic stress. Both contributions are shown to be increasing functions of a dimensionless shear rate Iv, in agreement with experimental results [4]. In contrast with a popular idea, the results suggest that lubrication may not necessarily reduce the contribution of contact forces to the bulk shear stress