No Arabic abstract
We report a combined experimental and simulation study of deformation-induced diffusion in compacted two-dimensional amorphous granular pillars, in which thermal fluctuations play negligible role. The pillars, consisting of bidisperse cylindrical acetal plastic particles standing upright on a substrate, are deformed uniaxially and quasistatically by a rigid bar moving at a constant speed. The plastic flow and particle rearrangements in the pillars are characterized by computing the best-fit affine transformation strain and non-affine displacement associated with each particle between two stages of deformation. The non-affine displacement exhibits exponential crossover from ballistic to diffusive behavior with respect to the cumulative deviatoric strain, indicating that in athermal granular packings, the cumulative deviatoric strain plays the role of time in thermal systems and drives effective particle diffusion. We further study the size-dependent deformation of the granular pillars by simulation, and find that different-sized pillars follow self-similar shape evolution during deformation. In addition, the yield stress of the pillars increases linearly with pillar size. Formation of transient shear lines in the pillars during deformation becomes more evident as pillar size increases. The width of these elementary shear bands is about twice the diameter of a particle, and does not vary with pillar size.
The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using emph{in situ} imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale Non-Smooth Contact Dynamics (NSCD) simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.
Nonequilibrium steady states of vibrated inelastic frictionless spheres are investigated in quasi-two-dimensional confinement via molecular dynamics simulations. The phase diagram in the density-amplitude plane exhibits a fluidlike disordered and an ordered phase with threefold symmetry, as well as phase coexistence between the two. A dynamical mechanism exists that brings about metastable traveling clusters and at the same time stable clusters with anisotropic shapes at low vibration amplitude. Moreover, there is a square bilayer state which is connected to the fluid by BKTHNY-type two-step melting with an intermediate tetratic phase. The critical behavior of the two continuous transitions is studied in detail. For the fluid-tetratic transition, critical exponents of $tilde{gamma}=1.73$, $eta_4 approx 1/4$, and $z=2.05$ are obtained.
Randomly crumpled sheets have shape memory. In order to understand the basis of this form of memory, we simulate triangular lattices of springs whose lengths are altered to create a topography with multiple potential energy minima. We then deform these lattices into different shapes and investigate their ability to retain the imposed shape when the energy is relaxed. The lattices are able to retain a range of curvatures. Under moderate forcing from a state of local equilibrium, the lattices deform by several percent but return to their retained shape when the forces are removed. By increasing the forcing until an irreversible motion occurs, we find that the transitions between remembered shapes show co-operativity among several springs. For fixed lattice structures, the shape memory tends to decrease as the lattice is enlarged; we propose ways to counter this decrease by modifying the lattice geometry. We survey the energy landscape by displacing individual nodes. An extensive fraction of these nodes proves to be bistable; they retain their displaced position when the energy is relaxed. Bending the lattice to a stable curved state alters the pattern of bistable nodes. We discuss this shapeability in the context of other forms of material memory and contrast it with the shapeability of plastic deformation. We outline the prospects for making real materials based on these principles.
Contrary to the theoretical predictions that all waves in two-dimensional disordered materials are localized, Anderson localization is observed only for sufficiently high frequencies in an isotropically jammed two-dimensional disordered granular packing of photoelastic disks. More specifically, we have performed an experiment in analyzing the level statistics of normal mode vibrations. We observe delocalized modes in the low-frequency boson-peak regime and localized modes in the high frequency regime with the crossover frequency just below the Debye frequency. We find that the level-distance distribution obeys Gaussian-Orthogonal-Ensemble (GOE) statistics, i.e. Wigner-Dyson distribution, in the boson-peak regime, whereas those in the high-frequency regime Poisson statistics is observed. The scenario is found to coincide with that of harmonic vibrational excitations in three-dimensional disordered solids.
Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state flow rate is reduced, the force distribution $P(f)$ changes only slightly, while there is a large change in the impulse distribution $P(i)$. In both cases, the distributions show an increase in small forces or impulses as the systems approach jamming, the opposite of that seen in previous Lennard-Jones simulations. This occurs dynamically as well for a hopper that transitions from a flowing to a jammed state over time. The final jammed $P(f)$ is quite distinct from a poured packing $P(f)$ in the same geometry. The change in $P(i)$ is a much stronger indicator of the approach to jamming. The formation of a peak or plateau in $P(f)$ at the average force is not a general feature of the approach to jamming.