Drag law for an intruder in granular sediments

Soft matter physics of the ground beneath our feet

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

Drag on a circular intruder traversing a shape-heterogeneous granular mixture

The main aim of our work is to explore the effect of particle shape heterogeneity on the dynamics of an intruder moving through a two-dimensional mixture of dumbbells and disks. In spite of similar physical conditions (the mass of the dumbbell is the same as that of the disk) and the same area fraction, we noticed a significant difference in the drag experienced by the intruder as the mixture concentration varies. The propagation of stress from the intruder to the granular grains manifests in the form of force chains, and interestingly these force chains can vary significantly depending on the shape of the grains. These differences, however, appear to be suppressed in the frictionless case where the force chains cannot extend very far from the initial point of contact. Apart from particle shape, the effect of the area fraction of the system and the size of the intruder have also been explored. As the area fraction increases, the drag force on the intruder increases owing to the increase in the contact forces. Finally, we present the velocity and stress fields at different dumbbell fractions and for various intruder diameters to show the effect of the moving intruder on its surrounding particles.

The concept of the mobilized domain: how it can explain and predict the forces exerted by a cohesive granular avalanche on an obstacle

ABSTRACT

The calculation of the impact pressure on obstacles in granular flows is a fundamental issue of practical relevance, e.g. for snow avalanches impacting obstacles. Previous research shows that the load on the obstacle builds up, due to the formation of force chains originating from the obstacle and extending into the granular material. This leads to the formation of a mobilized domain, wherein the flow is influenced by the presence of the obstacle. To identify the link between the physical mobilized domain properties and the pressure exerted on obstacles, we simulate subcritical cohesionless and cohesive avalanches of soft particles past obstacles with circular, rectangular or triangular cross-section using the Discrete Element Method. Our results show that the impact pressure decreases non-linearly with increasing obstacle width, regardless of the obstacle's cross-section. While the mobilized domain size is proportional to the obstacle width, the pressure decrease with increasing width originates from the jammed material inside the mobilized domain. We provide evidence that the compression inside the mobilized domain governs the pressure build-up for cohesionless subcritical granular flows. In the cohesive case, the stress transmission in the compressed mobilized domain is further enhanced, causing a pressure increase compared with the cohesionless case. Considering a kinetic and a gravitational contribution, we are able to calculate the impact pressure based on the properties of the mobilized domain. The equations used for the pressure calculation in this article may be useful in future predictive pressure calculations based on mobilized domain properties.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10035-021-01196-1.

Mobility in immersed granular materials upon cyclic loading

We study the mobility of objects embedded in an immersed granular packing and subjected to cyclic loadings. With this aim, we conducted uplift experiments whereby a horizontal plate is embedded in the packing and subjected to a vertical cyclic force oscillating between zero and a maximum amplitude. Tests performed at different cyclic force frequencies and amplitudes evidence the development of three mobility regimes whereby the plate stays virtually immobile, moves up steadily, or slowly creeps upwards. Results show that steady plate uplift can occur at lower force magnitudes when the frequency is increased. We propose an interpretation of this frequency-weakening behavior based on force relaxation experiments and on the analysis of the mobility response of theoretical viscoelastoplastic mechanical analog. These results and analysis point out inherent differences in mobility response between steady and cyclic loadings in immersed granular materials.

Simulation of dense non-Brownian suspensions with the lattice Boltzmann method: shear jammed and fragile states

Dense non-Brownian suspensions, including both hydrodynamic interactions and frictional contacts between particles, are numerically studied under simple and oscillatory shears in terms of the lattice Boltzmann method. We successfully reproduce the discontinuous shear thickening (DST) under a simple shear for bulk three-dimensional systems. For our simulation of an oscillatory shear in a quasi-two-dimensional system, we measure the mechanical response after the reduction of the strain amplitude from the initial oscillations. Here, we find the existence of a shear-jammed state under this protocol in which the storage modulus G' is only finite for high initial strain amplitude γI0. We also find the existence of a fragile state in which both fluid-like and solid-like responses can be detected for an identical area fraction and an initial strain amplitude γI0 depending on the initial phase Θ (or the asymmetricity of the applied strain) of the oscillatory shear. We also observe a DST-like behavior under the oscillatory shear in the fragile state. Moreover, we find that the stress anisotropy becomes large in the fragile state. Finally, we confirm that a stress formula based on the angular distribution of the contact force recovers the contact contributions to the stress tensors for both simple and oscillatory shears with large strains.

Burrowing dynamics of aquatic worms in soft sediments

We investigate the dynamics of Lumbriculus variegatus in water-saturated sediment beds to understand limbless locomotion in the benthic zone found at the bottom of lakes and oceans. These slender aquatic worms are observed to perform elongation-contraction and transverse undulatory strokes in both water-saturated sediments and water. Greater drag anisotropy in the sediment medium is observed to boost the burrowing speed of the worm compared to swimming in water with the same stroke using drag-assisted propulsion. We capture the observed speeds by combining the calculated forms based on resistive-force theory of undulatory motion in viscous fluids and a dynamic anchor model of peristaltic motion in the sediments. Peristalsis is found to be effective for burrowing in noncohesive sediments which fill in rapidly behind the moving body inside the sediment bed. Whereas the undulatory stroke is found to be effective in water and in shallow sediment layers where anchoring is not possible to achieve peristaltic motion. We show that such dual strokes occur as well in the earthworm Eisenia fetida which inhabits moist sediments that are prone to flooding. Our analysis in terms of the rheology of the medium shows that the dual strokes are exploited by organisms to negotiate sediment beds that may be packed heterogeneously and can be used by active intruders to move effectively from a fluid through the loose bed surface layer which fluidizes easily to the well-consolidated bed below.

Effective drag of a rod in fluid-saturated granular beds

We measure the drag encountered by a vertically oriented rod moving across a sedimented granular bed immersed in a fluid under steady-state conditions. At low rod speeds, the presence of the fluid leads to a lower drag because of buoyancy, whereas a significantly higher drag is observed with increasing speeds. The drag as a function of the depth is observed to decrease from being quadratic at low speeds to appearing more linear at higher speeds. By scaling the drag with the average weight of the grains acting on the rod, we obtain the effective friction μ_{e} encountered over six orders of magnitude of speeds. While a constant μ_{e} is found when the grain size, rod depth, and fluid viscosity are varied at low speeds, a systematic increase is observed as the speed is increased. We analyze μ_{e} in terms of the inertial number I and viscous number J to understand the relative importance of inertia and viscous forces, respectively. For sufficiently high fluid viscosities, we find that the effect of varying the speed, depth, and viscosity can be described by the empirical function μ_{e}=μ_{o}+kJ^{n}, where μ_{o} is the effective friction measured in the quasistatic limit, and k and n are material constants. The drag is then analyzed in terms of the effective viscosity η_{e} and found to decrease systematically as a function of J. We further show that η_{e} as a function of J is directly proportional to the fluid viscosity and the μ_{e} encountered by the rod.