No Arabic abstract
To understand the process of pattern formation in a low-density granular flow, we propose a simple particle model. This model considers spherical particles moving over an inclined flat surface based on three forces: gravity as the driving force, repulsive force due to particle collision, and drag force as the particle-- interaction through the ambient fluid. Numerical simulations of this model are conducted in two different types of two-dimensional planes, i.e., the monolayer was treated. In the horizontal plane parallel to the slope, particles aggregate at the moving front of the granular flow; and subsequently, flow instability occurs as a wavy pattern. This flow pattern is caused by the interparticle interaction arising from the drag force. Additionally, a vortex convection of particles is formed inside the aggregations. Meanwhile, in the vertical plane on the slope, particle aggregation is also found at the moving front of the granular flow. The aggregation resembles a head--tail structure, where the frontal angle against the slope approaches 60 degree from a larger angle as time progresses. Comparing the numerical result by varying the particle size, the qualitative dynamics of the granular flow are independent of size. To elucidate this reason, we perform a nondimensionalization for this model. The result indicates that our model can be simplified to dimensionless equations with one dimensionless parameter that represents the ratio of the gravity term to the excluded volume effect.
We perform three-dimensional simulations of a granular jet impact for both frictional and frictionless grains. Small shear stress observed in the experiment[X. Cheng et al., Phys. Rev. Lett. 99, 188001 (2007) ] is reproduced through our simulation. However, the fluid state after the impact is far from a perfect fluid, and thus, similarity between granular jets and quark gluon plasma is superficial, because the observed viscosity is finite and its value is consistent with the prediction of the kinetic theory.
We consider dense rapid shear flow of inelastically colliding hard disks. Navier-Stokes granular hydrodynamics is applied accounting for the recent finding cite{Luding,Khain} that shear viscosity diverges at a lower density than the rest of constitutive relations. New interpolation formulas for constitutive relations between dilute and dense cases are proposed and justified in molecular dynamics (MD) simulations. A linear stability analysis of the uniform shear flow is performed and the full phase diagram is presented. It is shown that when the inelasticity of particle collision becomes large enough, the uniform sheared flow gives way to a two-phase flow, where a dense solid-like striped cluster is surrounded by two fluid layers. The results of the analysis are verified in event-driven MD simulations, and a good agreement is observed.
We show that a freshly sedimented granular bed settles and creeps forward over extended periods of time under an applied hydrodynamic shear stress, which is below the critical value for bedload transport. The rearrangements are found to last over a time scale which is millions of times the sedimentation time scale of a grain in the fluid. Compaction occurs uniformly throughout the bed, but creep is observed to decay exponentially with depth, and decreases over time. The granular volume fraction in the bed is found to increase logarithmically, saturating at the random close packing value $phi_{rcp} approx 0.64$, while the surface roughness is observed to remain essentially unchanged. We demonstrate that an increasingly higher shear stress is required to erode the bed after a sub-critical shear is applied which results in an increase in its volume fraction. Thus, we find that bed armoring occurs due to a deep shear-induced relaxation of the bed towards the volume fraction associated with the glass transition.
Flexible barriers are increasingly used for the protection from debris flow in mountainous terrain due to their low cost and environmental impact. However, a numerical tool for rational design of such structures is still missing. In this work, a hybrid computational framework is presented, using a total Lagrangian formulation of the Finite Element Method (FEM) to represent a flexible barrier. The actions exerted on the structure by a debris flow are obtained from simultaneous simulations of the flow of a fluid-grain mixture, using two conveniently coupled solvers: the Discrete Element Method (DEM) governs the motion of the grains, while the free-surface non-Newtonian fluid phase is solved using the Lattice-Boltzmann Method (LBM). Simulations on realistic geometries show the dependence of the momentum transfer on the barrier on the composition of the debris flow, challenging typical assumptions made during the design process today. In particular, we demonstrate that both grains and fluid contribute in a non-negligible way to the momentum transfer. Moreover, we show how the flexibility of the barrier reduces its vulnerability to structural collapse, and how the stress is distributed on its fabric, highlighting potential weak points.
Experimental measurements of particle dynamics on the lower surface of a 3D Couette cell containing monodisperse spheres are reported. The average radial density and velocity profiles are similar to those previously measured within the bulk and on the lower surface of the 3D cell filled with mustard seeds. Observations of the evolution of particle velocities over time reveal distinct motion events, intervals where previously stationary particles move for a short duration before jamming again. The cross-correlation between the velocities of two particles at a given distance $r$ from the moving wall reveals a characteristic lengthscale over which the particles are correlated. The autocorrelation of a single particles velocity reveals a characteristic timescale $tau$ which decreases with distance from the inner moving wall. This may be attributed to the increasing rarity at which the discrete motion events occur and the reduced duration of those events at large $r$. The relationship between the RMS azimuthal velocity fluctuations, $delta v_theta(r)$, and average shear rate, $dotgamma(r)$, was found to be $delta v_theta propto dotgamma^alpha$ with $alpha = 0.52 pm 0.04$. These observations are compared with other recent experiments and with the modified hydrodynamic model recently introduced by Bocquet et al.