ised in terms of porosity, coordination number and radial distribution function, and the particle–particle or pore–pore connectivity is analysed with the Voronoi or Delaunay tessellation. The force structures and their link with jamming states are also discussed. Structural properties are shown to be correlated with porosity which links directly to the cohesive force acting on a particle.

We present a numerical study of the packing of uniform spheres under three-dimensional vibration using the discrete element method (DEM), focusing on the effects of vibration condition (amplitude and frequency) and inter-particle frictions (sliding and rolling frictions). The results are analysed in terms of packing density, coordination number (CN), radial distribution function (RDF) and pore structure. It is shown that increasing either the vibration amplitude or frequency causes packing density to increase initially to a maximum and then decrease. Both vibration frequency and amplitude should be considered to characterize the effect of vibration process on packing structure. The sliding and rolling frictions between particles can decrease packing density since they dissipate energy, although the effect of rolling friction is less significant. In line with the change of packing density, microstructural properties such as CN, RDF and pore distribution also change: a looser packing often corresponds to smaller CN, less peaked RDF and larger but more widely distributed pores.

We present a numerical method capable of reproducing the densification process from the so-called random loose to dense packing of uniform spheres under vertical vibration. The effects of vibration amplitude and frequency are quantified, and the random close packing is shown to be achieved only if both parameters are properly controlled. Two densification mechanisms are identified: pushing filling by which the contact between spheres is maintained and jumping filling by which the contact between particles is periodically broken. In general, pushing filling occurs when the vibration intensity is low and jumping filling becomes dominant when the vibration intensity is high.

It is shown that by properly controlling vibrational and charging conditions, the transition from disordered to ordered, densest packing of particles can be obtained consistently. The method applies to both spherical and nonspherical particles. For spheres, face centered cubic packing with different orientations can be achieved by monitoring the vibration amplitude and frequency, and the structure of the bottom layer, in particular. The resultant force structures are ordered but do not necessarily correspond to the packing structures fully. The implications of the findings are also discussed.

Micro structures of equal sphere packing (ranging from loose to dense packing) generated numerically by discrete element method under different vibration conditions are characterized using Voronoi/Delaunay tessellation, which is applied on a wide range of packing densities. The analysis on micro properties such as the total perimeter, surface area, and the face number distribution of each Voronoi polyhedron, and the pore size distribution in each Voronoi/Delaunay subunit is systematically carried out. The results show that with the increasing density of sphere packing, the Voronoi/Delaunay pore size distribution is narrowed. That indicates large pores to be gradually substituted by small uniformed ones during densification. Meanwhile, the distributions of face number, total perimeter, and surface area of Voronoi polyhedra at high packing densities tend to be narrower and higher, which is in good agreement with those in random loose packing.

This paper presents a study of the structural transition of the packing of uniform spheres with a wide range of packing fractions. Several critical states are identified including the onset of local jamming, onset of global jamming and maximally random jammed (MRJ) states. With the increase of packing fraction , the packing structure transforms from one-dimensional chains to two-dimensional triangles and finally three-dimensional tetrahedra, correspondingly undergoing phase changes from non-jamming to local jamming, global jamming, MRJ and finally crystal structure. There is a competition between FCC (face-centred cubic) and HCP (hexagonal closed packed) in the transition from disordered to ordered structure with the increase of . HCP can transform to FCC after reaching its maximum at =0.7

The crystallization, corresponding to the fcc structure (with packing density p≈ 0.74), of smooth equal hard spheres under batch-wised feeding and three-dimensional interval vibration is numerically obtained by using the discrete element method. The numerical experiment shows that the ordered packing can be realized by proper control of the dynamic parameters such as batch of each feeding ξ and vibration amplitude A. The radial distribution function and force network are used to characterize the ordered structure. The defect formed during vibrated packing is characterized as well. The results in our work fill the gap of getting packing density between random close packing and fcc packing in phase diagram which provides an effective way of theoretically investigating the complex process and mechanism of hard sphere crystallization and its dynamics.