skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Contact network changes in ordered and disordered disk packings
We investigate the mechanical response of packings of purely repulsive, frictionless disks to quasistatic deformations. The deformations include simple shear strain at constant packing fraction and at constant pressure, “polydispersity” strain (in which we change the particle size distribution) at constant packing fraction and at constant pressure, and isotropic compression. For each deformation, we show that there are two classes of changes in the interparticle contact networks: jump changes and point changes. Jump changes occur when a contact network becomes mechanically unstable, particles “rearrange”, and the potential energy (when the strain is applied at constant packing fraction) or enthalpy (when the strain is applied at constant pressure) and all derivatives are discontinuous. During point changes, a single contact is either added to or removed from the contact network. For repulsive linear spring interactions, second- and higher-order derivatives of the potential energy/enthalpy are discontinuous at a point change, while for Hertzian interactions, third- and higher-order derivatives of the potential energy/enthalpy are discontinuous. We illustrate the importance of point changes by studying the transition from a hexagonal crystal to a disordered crystal induced by applying polydispersity strain. During this transition, the system only undergoes point changes, with no jump changes. We emphasize that one must understand point changes, as well as jump changes, to predict the mechanical properties of jammed packings.  more » « less
Award ID(s):
2002782 1755494 2002797
PAR ID:
10248611
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Soft Matter
Volume:
16
Issue:
41
ISSN:
1744-683X
Page Range / eLocation ID:
9443 to 9455
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. none (Ed.)
    The mechanical behavior of crystal-rich mush controls the dynamics and evolution of magma bodies but is poorly understood, . The presence of a semi-rigid contact network in the crystal phase greatly affects the rheology of mush but the contributions of crystal shape to the force, contact and shape fabric remains poorly characterized. This in turn influences the transmission of stress in the mush, the packing stiffness, and the volume fraction at jamming. It is also unclear whether the total amount of deformation of a mush can be quantitatively determined from the shape preferred orientation of the crystals. We performed 3D numerical simulations using a coupled computational fluid dynamics and discrete element method to illuminate the dynamic states of non-spherical crystals in a viscous melt. Simulations consisted of the simple shear of a mush under a constant pressure upper boundary. Crystals are represented by elongated cuboids having an aspect ratio of four. Our results differ from those obtained with smooth spheres and shed light on the influence of the crystal shape and orientation fabric on the mechanical properties of a mush. We found two distinct behaviors associated with the transient and steady-state, however at all times strain is nonaffine. During the transient, the strain is accommodated by the emergence of multiple shear bands and tends to concentrate on a single one. The shear bands emerge because of steric blocking and space limitations preventing the rotation of elongated particles, generating the local and temporary jamming of the crystal network. On the contrary, in the residual and steady-state, the strain is accommodated by one main shear band. The analysis of the orientation of the crystals shows that the deformation of the mush tends to increase the foliation of the crystals more than their alignment. 
    more » « less
  2. We perform computational studies of jammed particle packings in two dimensions undergoing isotropic compression using the well-characterized soft particle (SP) model and deformable particle (DP) model that we developed for bubbles and emulsions. In the SP model, circular particles are allowed to overlap, generating purely repulsive forces. In the DP model, particles minimize their perimeter, while deforming at the fixed area to avoid overlap during compression. We compare the structural and mechanical properties of jammed packings generated using the SP and DP models as a function of the packing fraction ρ, instead of the reduced number density φ. We show that near jamming onset the excess contact number Δz=z-z J and shear modulus G scale as Δρ 0.5 in the large system limit for both models, where Δρ=ρ-ρ J and z J ≈4 and ρ J ≈0.842 are the values at jamming onset. Δz and G for the SP and DP models begin to differ for ρ≥0.88. In this regime, Δz∼G can be described by a sum of two power-laws in Δρ, i.e. Δz∼G∼C 0 Δρ 0.5 +C 1 Δρ 1.0 to lowest order. We show that the ratio C 1 /C 0 is much larger for the DP model compared to that for the SP model. We also characterize the void space in jammed packings as a function of ρ. We find that the DP model can describe the formation of Plateau borders as ρ→1.0. We further show that the results for z and the shape factor A versus ρ for the DP model agree with recent experimental studies of foams and emulsions. 
    more » « less
  3. We study packings of bidispersed spherical particles on a spherical surface. The presence of curvature necessitates defects even for monodispersed particles; bidispersity either leads to a more disordered packing for nearly equal radii, or a higher fill fraction when the smaller particles are accommodated in the interstices of the larger spheres. Variation in the packing fraction is explained by a percolation transition, as chains of defects or scars previously discovered in the monodispersed case grow and eventually disconnect the neighbor graph. 
    more » « less
  4. Abstract Static granular packings play a central role in numerous industrial applications and natural settings. In these situations, fluid or fine particle flow through a bed of static particles is heavily influenced by the narrowest passage connecting the pores of the packing, commonly referred to as pore throats, or constrictions. Existing studies predominantly assume monodisperse rigid particles, but this is an oversimplification of the problem. In this work, we illustrate the connection between pore throat size, polydispersity, and particle deformation in a packed bed of spherical particles. Simple analytical expressions are provided to link these properties of the packing, followed by examples from Discrete Element Method (DEM) simulations of fine particle percolation demonstrating the impact of polydispersity and particle deformation. Our intent is to emphasize the substantial impact of polydispersity and particle deformation on constriction size, underscoring the importance of accounting for these effects in particle transport in granular packings. 
    more » « less
  5. All-nitrogen solids, if successfully synthesized, are ideal high-energy-density materials because they store a great amount of energy and produce only harmless N 2 gas upon decomposition. Currently, the only method to obtain all-nitrogen solids is to apply high pressure to N 2 crystals. However, products such as cg-N tend to decompose upon releasing the pressure. Compared to covalent solids, molecular crystals are more likely to remain stable during decompression because they can relax the strain by increasing the intermolecular distances. The challenge of such a route is to find a molecular crystal that can attain a favorable phase under elevated pressure. In this work, we show, by designing a novel N 16 molecule (tripentazolylamine) and examining its crystal structures under a series of pressures, that the aromatic units and high molecular symmetry are the key factors to achieving an all-nitrogen molecular crystal. Density functional calculations and structural studies reveal that this new all-nitrogen molecular crystal exhibits a particularly slow enthalpy increase with pressure due to the highly efficient crystal packing of its highly symmetric molecules. Vibration mode calculations and molecular dynamics (MD) simulations show that N 16 crystals are metastable at ambient pressure and could remain inactive up to 400 K. The initial reaction steps of the decomposition are calculated by following the pathway of the concerted excision of N 2 from the N 5 group as revealed by the MD simulations. 
    more » « less