Mohammad Saadatfar

Geometrical structure of disordered sphere packings

The three-dimensional structure of large packings of monosized spheres with volume fractions ranging between 0.58 and 0.64 has been studied with x-ray computed tomography. We search for signatures of organization, classifying local arrangements and exploring the effects of local geometrical constrains on the global packing. This study is the largest and the most accurate empirical analysis of disordered packings at the grain-scale to date, mapping over 380,000 sphere coordinates with precision within 0.1% of the sphere diameters. We discuss topological and geometrical methods to characterize and classify these systems emphasizing the implications that local geometry can have on the mechanisms of formation of these amorphous structures. Some of the main results are (1) the observation that the average number of contacts increases with the volume fraction; (2) the discovery that these systems have a very compact contact network; (3) the finding that disordered packing can be locally more efficient than crystalline packings; (4) the observation that the peaks of the radial distribution function follow power law divergences; (5) the discovery that geometrical frustration plays no role in the formation of such amorphous packings.

Investigating the geometrical structure of disordered sphere packings

Bead packs of up to 150,000 mono-sized spheres with packing densities ranging from 0.58 to 0.64 have been studied by means of X-ray computed tomography. These studies represent the largest and the most accurate description of the structure of disordered packings at the grain-scale ever attempted. We investigate the geometrical structure of such packings looking for signatures of disorder. We discuss ways to characterize and classify these systems and the implications that local geometry can have on densification dynamics.

Onset of mechanical stability in random packings of frictional spheres

Using sedimentation to obtain precisely controlled packings of noncohesive spheres, we find that the volume fraction phiRLP of the loosest mechanically stable packing is in an operational sense well defined by a limit process. This random loose packing volume fraction decreases with decreasing pressure p and increasing interparticle friction coefficient mu. Using x-ray tomography to correct for a container boundary effect that depends on particle size, we find for rough particles in the limit p-->0 a new lower bound, phiRLP=0.550+/-0.001.

Elastic and transport properties of cellular solids derived from three-dimensional tomographic images

We describe a three-dimensional imaging and analysis study of eight industrial cellular foam morphologies. The foam morphologies were generated by differing industrial processing methods. Tomograms are acquired on an X-ray micro-computed tomography facility at scales of approximately equal to at resolutions down to 7 μm. The image quality is sufficient in all cases to measure local structure and connectivity of the foamed material, and the field of view large enough to calculate a range of material properties. Phase separation into solid and porous components is straightforward. Three-dimensional structural characteristics are measured directly on the porous and solid phases of the images. A number of morphological parameters are obtained, including pore volume-to-surface-area ratio, connectivity, the pore and solid phase size distributions defined by maximal sphere openings and chord length measurements. We further calculate the pore size distribution associated with capillary pressure via simulating of mercury drainage on the digital images. The binarized microstructures are used as a basis for calculations of transport properties (fluid permeability, diffusivity and thermal conductivity) and elastic moduli. From the data, we generate property–porosity relationships for the range of foam morphologies imaged and quantitatively analyse the effects of porosity and microstructure on the resultant properties of the foams. We compare our numerical data to commonly used theoretical and empirical property–porosity relationships. For thermal conductivity, we find that the numerical results agree extremely well with an empirical expression based on experimental data of various foams. The upper Hashin–Shtrikman bound also provides an excellent prediction of the data across all densities. From simulation of the diffusivity, we can define the tortuosity of the pore space within the cellular solid. We find that different processing methods lead to strong variations in the tortuosity of the pore space of the foams. For elastic properties, our results show that for the Young modulus, E, both the differential effective medium theory and the classical correlation give a good correlation. Assuming a constant Poissons ratio leads to reasonable agreement. The best correlation for is given by assuming a slight variation in as a linear function of porosity. The permeability of the foams varies over three orders of magnitude. Correlations for permeability based on the classical Kozeny–Carman equation lead to reasonable agreement, except at the lowest porosities. Permeability estimations based on mercury porosimetry give excellent agreement for all foams.

Local and global relations between the number of contacts and density in monodisperse sphere packs

The topological structure resulting from the network of contacts between grains (contact network) is studied for large samples of monosized spheres with densities (fraction of volume occupied by the spheres) ranging from 0.59 to 0.64. We retrieve the coordinates of each bead in the pack and we calculate the average coordination number by using three different methods. We show that, in the range of density investigated, the coordination number is larger than four and it increases with the packing fraction. At local level we also observe a positive correlation between local packing fraction and number of neighbours. We discover a dependence between the local densities of configurations with few neighbours in contact and the global sample densities. This might indicate that local configurations with small numbers of neighbours are able to deform plastically when the sample is compactifying.

An x-ray tomography facility for quantitative prediction of mechanical and transport properties in geological, biological, and synthetic systems

A fully integrated X-ray tomography facility with the ability to generate tomograms with 20483 voxels at 2 micron spatial resolution was built to satisfy the requirements of a virtual materials testing laboratory. The instrument comprises of a continuously pumped micro-focus X-ray gun, a milli-degree rotation stage and a high resolution and large field X-ray camera, configured in a cone beam geometry with a circular trajectory. The purpose of this facility is to routinely analyse and investigate real world biological, geological and synthetic materials at a scale in which the traditional domains of physics, chemistry, biology and geology merge. During the first 2 years of operation, approximately 4 Terabytes of data have been collected, processed and analysed, both as static and in some cases as composite dynamic data sets. This incorporates over 300 tomograms with 10243 voxels and 50 tomograms with 20483 voxels for a wide range of research fields. Specimens analysed include sedimentary rocks, soils, bone, soft tissue, ceramics, fibre-reinforced composites, foams, wood, paper, fossils, sphere packs, bio-morphs and small animals. In this paper, the flexibility of the facility is highlighted with some prime examples.

Local Origin of Global Contact Numbers in Frictional Ellipsoid Packings

In particulate soft matter systems the average number of contacts Z of a particle is an important predictor of the mechanical properties of the system. Using x-ray tomography, we analyze packings of frictional, oblate ellipsoids of various aspect ratios α, prepared at different global volume fractions ϕg. We find that Z is a monotonically increasing function of ϕg for all α. We demonstrate that this functional dependence can be explained by a local analysis where each particle is described by its local volume fraction ϕl computed from a Voronoi tessellation. Z can be expressed as an integral over all values of ϕl: Z(ϕg,α,X)=∫Zl(ϕl,α,X)P(ϕl|ϕg)dϕl. The local contact number function Zl(ϕl,α,X) describes the relevant physics in term of locally defined variables only, including possible higher order terms X. The conditional probability P(ϕl|ϕg) to find a specific value of ϕl given a global packing fraction ϕg is found to be independent of α and X. Our results demonstrate that for frictional particles a local approach is not only a theoretical requirement but also feasible.

Non-universal Voronoi cell shapes in amorphous ellipsoid packs

In particulate systems with short-range interactions, such as granular matter or simple fluids, local structure determines the macroscopic physical properties. We analyse local structure metrics derived from the Voronoi diagram of oblate ellipsoids, for various aspect ratios and global packing fractions φg. We focus on jammed static configurations of frictional ellipsoids, obtained by tomographic imaging and by discrete element method simulations. The rescaled distribution of local packing fractions φl, defined as the ratio of particle volume and its Voronoi cell volume, is found to be independent of the particle aspect ratio, and coincide with results for sphere packs. By contrast, the typical Voronoi cell shape, quantified by the Minkowski tensor anisotropy index β = β02,0, points towards a difference between random packings of spheres and those of oblate ellipsoids. While the average cell shape β of all cells with a given value of is similar in dense and loose jammed sphere packings, the structure of dense and loose ellipsoid packings differs substantially such that this does not hold true.

Percolating length scales from topological persistence analysis of micro-CT images of porous materials

This work was funded in part by ARCgrant DP110102888 (to A.P.S. and V.R.),and ARC Future FellowshipsFT100100470 (A.P.S.) and FT140100604(V.R.). Financial support was alsoprovided by the member companiesof the ANU/UNSW DigicoreConsortium, who have alsocontributed samples for this study.

Tomographic Analysis of Jammed Ellipsoid Packings

Disordered packings of ellipsoidal particles are an important model for disordered granular matter. Here we report a way to determine the average contact number of ellipsoid packings from tomographic analysis. Tomographic images of jammed ellipsoid packings prepared by vertical shaking of loose configurations are recorded and the positions and orientations of the ellipsoids are reconstructed. The average contact number can be extracted from a contact number scaling (CNS) function. The size of the particles, that may vary due to production inaccuracies, can also be determined by this method.