Mechanical characterization of partially crystallized sphere packings

Compression and fracture of ordered and disordered droplet rafts

We simulate a two-dimensional array of droplets being compressed between two walls. The droplets are adhesive due to an attractive depletion force. As one wall moves toward the other, the droplet array is compressed and eventually induced to rearrange. The rearrangement occurs via a fracture, where depletion bonds are quickly broken between a subset of droplets. For monodisperse, hexagonally ordered droplet arrays, this fracture is preceded by a maximum force exerted on the walls, which drops rapidly after the fracture occurs. In small droplet arrays a fracture is a single well-defined event, but for larger droplet arrays, competing fractures can be observed. These are fractures nucleated nearly simultaneously in different locations. Finally, we also study the compression of bidisperse droplet arrays. The addition of a second droplet size further disrupts fracture events, showing differences between ideal crystalline arrays, crystalline arrays with a small number of defects, and fully amorphous arrays. These results are in good agreement with previously published experiments.

Is there a granular potential?

Granular materials, such as sand or grain, exhibit many structural and dynamic characteristics similar to those observed in molecular systems, despite temperature playing no role in their properties. This has led to an effort to develop a statistical mechanics for granular materials that has focused on establishing an equivalent to the microcanonical ensemble and a temperature-like thermodynamic variable. Here, we expand on these ideas by introducing a granular potential into the Edwards ensemble, as an analogue to the chemical potential, and explore its properties using a simple model of a granular system. A simple kinetic Monte Carlo simulation of the model shows the effect of mass transport leading to equilibrium and how this is connected to the redistribution of volume in the system. An exact analytical treatment of the model shows that the compactivity and the ratio of the granular potential to the compactivity determine the equilibrium between two open systems that are able to exchange volume and particles, and that mass moves from high to low values of this ratio. Analysis of the granular potential shows that adding a particle to the system increases the entropy at high compactivity, but decreases the entropy at low compactivity. Finally, we demonstrate the use of a small system thermodynamics method for the measurement of granular potential differences.

Mixing properties of bi-disperse ellipsoid assemblies: mean-field behaviour in a granular matter experiment

The structure and spatial statistical properties of amorphous ellipsoid assemblies have profound scientific and industrial significance in many systems, from cell assays to granular materials. This paper uses a fundamental theoretical relationship for mixture distributions to explain the observations of an extensive X-ray computed tomography study of granular ellipsoidal packings. We study a size-bi-disperse mixture of two types of ellipsoids of revolutions that have the same aspect ratio of α ≈ 0.57 and differ in size, by about 10% in linear dimension, and compare these to mono-disperse systems of ellipsoids with the same aspect ratio. Jammed configurations with a range of packing densities are achieved by employing different tapping protocols. We numerically interrogate the final packing configurations by analyses of the local packing fraction distributions calculated from the Voronoi diagrams. Our main finding is that the bi-disperse ellipsoidal packings studied here can be interpreted as a mixture of two uncorrelated mono-disperse packings, insensitive to the compaction protocol. Our results are consolidated by showing that the local packing fraction shows no correlation beyond their first shell of neighbours in the binary mixtures. We propose a model of uncorrelated binary mixture distribution that describes the observed experimental data with high accuracy. This analysis framework will enable future studies to test whether the observed mean-field behaviour is specific to the particular granular system or the specific parameter values studied here or if it is observed more broadly in other bi-disperse non-spherical particle systems.

A versatile microtomography system to study in situ the failure and fragmentation in geomaterials

This article describes a microtomography experimental platform enabling in situ micro-mechanical study of failure and fragmentation in geomaterials. The system is based on an original high-pressure triaxial flow cell, which is fully integrated into a custom built microtomography scanner equipped with a laboratory x-ray source. The design of the high-precision mechanical apparatus was informed by the concurrent development of advanced tomographic reconstruction methods based on helical scanning and of algorithms correcting for hardware inaccuracies. This experimental system produces very high-quality 3D images of microstructural changes occurring in rocks undergoing mechanical failure and substantial fragmentation. We present the results of two experiments as case studies to demonstrate the capabilities and versatility of this instrumental platform. These experiments tackle various questions related to the onset of rock failure, the hydromechanical coupling and relaxation mechanisms in fractured rocks, or the fragmentation process in geomaterials such as copper ores.

Topological signatures and stability of hexagonal close packing and Barlow stackings

Two common representations of close packings of identical spheres consisting of hexagonal layers, called Barlow stackings, appear abundantly in minerals and metals. These motifs, however, occupy an identical portion of space and bear identical first-order topological signatures as measured by persistent homology. Here we present a novel method based on k-fold covers that unambiguously distinguishes between these patterns. Moreover, our approach provides topological evidence that the FCC motif is the more stable of the two in the context of evolving experimental sphere packings during the transition from disordered to an ordered state. We conclude that our approach can be generalised to distinguish between various Barlow stackings manifested in minerals and metals.

Mechanical properties of 2D aggregates of oil droplets as model mono-crystals

We investigate the elastic and yielding properties of two dimensional defect-free mono-crystals made of highly monodisperse droplets. Crystals are compressed between two parallel boundaries of which one acts as a force sensor. As the available space between boundaries is reduced, the crystal goes through successive row-reduction transitions. For small compression forces, the crystal responds elastically until a critical force is reached and the assembly fractures in a single catastrophic global event. Correspondingly there is a peak in the force measurement associated with each row-reduction. The elastic properties of ideal mono-crystal samples are fully captured by a simple analytical model consisting of an assembly of individual capillary springs. The yielding properties of the crystal are captured with a minimal bond breaking model.

A neutron tomography study: probing the spontaneous crystallization of randomly packed granular assemblies

We study the spontaneous crystallization of an assembly of highly monodisperse steel spheres under shaking, as it evolves from localized icosahedral ordering towards a packing reaching crystalline ordering. Towards this end, real space neutron tomography measurements on the granular assembly are carried out, as it is systematically subjected to a variation of frequency and amplitude. As expected, we see a presence of localized icosahedral ordering in the disordered initial state (packing fraction ≈ 0.62). As the frequency is increased for both the shaking amplitudes (0.2 and 0.6 mm) studied here, there is a rise in packing fraction, accompanied by an evolution to crystallinity. The extent of crystallinity is found to depend on both the amplitude and frequency of shaking. We find that the icosahedral ordering remains localized and its extent does not grow significantly, while the crystalline ordering grows rapidly as an ordering transition point is approached. In the ordered state, crystalline clusters of both face centered cubic (FCC) and hexagonal close packed (HCP) types are identified, the latter of which grows from stacking faults. Our study shows that an earlier domination of FCC gives way to HCP ordering at higher shaking frequencies, suggesting that despite their coexistence, there is a subtle dynamical competition at play. This competition depends on both shaking amplitude and frequency, as our results as well as those of earlier theoretical simulations demonstrate. It is likely that this involves the very small free energy difference between the two structures.

Dynamic modelling on the confined crystallization of mono-sized cubic particles under mechanical vibration

The dynamic crystallization of cubic granular particles under three-dimensional mechanical vibration is numerically investigated by the discrete element method. The effects of operational conditions (vibration, container shape and system size) and particle properties (gravity and friction) on the formation of crystals and defects are discussed. The results show that the formation and growth of clusters with face-to-face aligned cubic particles can be easily realized under vibrations. Especially, a single crystal with both translational and orientational ordering can be reproduced in a rectangular container under appropriate vibrations. It is also found that the gravitational effect is beneficial for the ordering of a packing; the ordering of frictional particles can be improved significantly with an enlarged gravitational acceleration. The flat walls of a rectangular container facilitate the formation of orderly layered structures. The curved walls of a cylindrical container contribute to the formation of ring-like structures, whereas they also cause distortions and defects in the packing centers. Finally, it is shown that the crystallization of inelastic particles is basically accomplished by the pursuit of a better mechanical stability of the system, with decreasing kinetic and potential energies.

Nucleation in Sheared Granular Matter

We present an experiment on crystallization of packings of macroscopic granular spheres. This system is often considered to be a model for thermally driven atomic or colloidal systems. Cyclically shearing a packing of frictional spheres, we observe a first order phase transition from a disordered to an ordered state. The ordered state consists of crystallites of mixed fcc and hcp symmetry that coexist with the amorphous bulk. The transition, initiated by homogeneous nucleation, overcomes a barrier at 64.5% volume fraction. Nucleation consists predominantly of the dissolving of small nuclei and the growth of nuclei that have reached a critical size of about ten spheres.

Pore configuration landscape of granular crystallization

Uncovering grain-scale mechanisms that underlie the disorder–order transition in assemblies of dissipative, athermal particles is a fundamental problem with technological relevance. To date, the study of granular crystallization has mainly focussed on the symmetry of crystalline patterns while their emergence and growth from irregular clusters of grains remains largely unexplored. Here crystallization of three-dimensional packings of frictional spheres is studied at the grain-scale using X-ray tomography and persistent homology. The latter produces a map of the topological configurations of grains within static partially crystallized packings. Using numerical simulations, we show that similar maps are measured dynamically during the melting of a perfect crystal. This map encodes new information on the formation process of tetrahedral and octahedral pores, the building blocks of perfect crystals. Four key formation mechanisms of these pores reproduce the main changes of the map during crystallization and provide continuous deformation pathways representative of the crystallization dynamics.

Emergence and growth of crystalline domains in granular media remains under-explored. Here, the authors analyse tomographic snapshots from partially recrystallized packings of spheres using persistent homology and find agreement with proposed transitions based on continuous deformation of octahedral and tetrahedral voids.

Analyzing X-ray tomographies of granular packings

Starting from three-dimensional volume data of a granular packing, as, e.g., obtained by X-ray Computed Tomography, we discuss methods to first detect the individual particles in the sample and then analyze their properties. This analysis includes the pair correlation function, the volume and shape of the Voronoi cells, and the number and type of contacts formed between individual particles. We mainly focus on packings of monodisperse spheres, but we will also comment on other monoschematic particles such as ellipsoids and tetrahedra. This paper is accompanied by a package of free software containing all programs (including source code) and an example three-dimensional dataset which allows the reader to reproduce and modify all examples given.

From Crystals to Disordered Crystals: A Hidden Order-Disorder Transition

To distinguish between order and disorder is of fundamental importance to understanding solids. It becomes more significant with recent observations that solids with high structural order can behave like disordered solids, while properties of disordered solids can approach crystals under certain circumstance. It is then imperative to understand when and how disorder takes effect to deviate the properties of a solid from crystals and what the correct factors are to control the behaviours of solids. Here we answer these questions by reporting the finding of a hidden order-disorder transition from crystals to disordered crystals for static packings of frictionless spheres. While the geometric indicators are mostly blind to the transition, disordered crystals already exhibit properties apart from crystals. The transition approaches the close packing of hard spheres, giving rise to the singularity of the close packing point. We evidence that both the transition and properties of disordered crystals are jointly determined by the structural order and density. Near the transition, the elastic moduli and coordination number of disordered crystals show particular pressure dependence distinct from known behaviours of both crystals and jammed solids. The discovery of the transition therefore reveals some unknown aspects of solids.

Structural and mechanical features of the order-disorder transition in experimental hard-sphere packings

Here we present an experimental and numerical investigation on the grain-scale geometrical and mechanical properties of partially crystallized structures made of macroscopic frictional grains. Crystallization is inevitable in arrangements of monosized hard spheres with packing densities exceeding Bernal's limiting density ϕ(Bernal)≈0.64. We study packings of monosized hard spheres whose density spans over a wide range (0.59<ϕ<0.72). These experiments harness x-ray computed tomography, three-dimensional image analysis, and numerical simulations to access precisely the geometry and the 3D structure of internal forces within the sphere packings. We show that clear geometrical transitions coincide with modifications of the mechanical backbone of the packing both at the grain and global scale. Notably, two transitions are identified at ϕ(Bernal)≈0.64 and ϕ(c)≈0.68. These results provide insights on how geometrical and mechanical features at the grain scale conspire to yield partially crystallized structures that are mechanically stable.