Selection of material for X-ray tomography analysis and DEM simulations: comparison between granular materials of biological and non-biological origins

Physical properties and X-ray tomography images of ﬁve different granular materials: glass bead, glass grit, short grain white rice, sorghum and clay granules (Seramis) are investigated to select the most promising materials for numerical simulations and time-lapse X-ray tomography imaging. The examined materials represent granular materials of non-biological origin (glass bead, glass grit and Seramis) and granular plant materials (rice and sorghum). The choice of materials was dictated by their mechanical and packing properties, on one hand, and 3D imaging qualities of granular materials and quantitative particle shape analysis, on the other hand. It was found that materials of biological origin, i.e. rice and sorghum, showed promising characteristics and may be considered as the most appropriate materials for comparison between time-lapse X-ray tomography experiments and numerical simulation.


Introduction
Investigation of granular materials has been driving research and development interest for decades, so as to better understand mechanics an interparticle interactions during different technological processes. There are numerous types of processes involving granular materials (e.g. compression, shearing, sintering, silo filling and silo discharge), and bulk solids appear in various forms (e.g. dry grains and powders, slurries, suspensions or pastes), which make the spectrum of investigation of granular materials very large. A review of the literature reveals that various granular materials, used to different industrial processes, have been investigated over last few decades. The studies were conducted for grain beddings composed of particles of different geometrical and mechanical properties. The studies included grain beddings composed of small particles with submillimeter size [1][2][3][4][5] and the ones comprising larger elements [6][7][8][9][10]. Properties of granular assemblies were found to be different due to their origination. Granular materials of biological origin (e.g. seeds and grains), commonly used in agriculture and food engineering, reveal mechanical characteristics different than the ones obtained for mineral or plastic granular materials. Granular plant materials are usually hydroscopic and change their mechanical properties through the absorption of moisture [9,11]. Their peculiar behavior requires both, specific approach in determining the parameters and specific setup of simulations. A review of the literature shows that numerous studies have been performed so far to measure material properties of particulate materials; however, values of parameters for materials of biological origin are uncertain and dispersed across various reports [11].
For many years, extensive efforts have been made to examine properties of granular packings and analyze effects occurring in grain bedding by using experimental methods; however, in many cases, these methods appeared to be insufficient. Because experimental methods do not provide essential insight into the properties of particulate assemblies and interactions between grains, non destructive imaging methods such as X-ray tomography and computational approaches are increasingly employed to examine granular media.
In mechanics and physics, two different approaches can be used to describe and model granular materials [22]. The first approach is based on continuum theory and relies on empirical assumptions about the macroscopic material behavior. The second approach, in which the macroscopic analysis is complemented by a microscopic description of the material, enables modeling individual particles and their interactions. Application of computational techniques for modeling mechanical processes involving particulate solids requires knowledge of their structural and mechanical properties. That requirement frequently results in problems with selection of research material. A selection of material is not easy matter, since it can be governed by constraints, which are sometimes not compliant. For instance, application of numerical method such as discrete element method (DEM) is mainly constrained by number and shape of modelled particles. Bharadwaj [23] reported that current models were restricted to modelling 1 million spherical particles which could be considered to be a limit number of particles in DEM simulations. The degree of complexity increases when taking into account secondary forces (e.g. cohesion), nonspherical shapes of particles, mobile and flexible boundaries or complex contact force models. Among the many parameters which determine effectiveness of the discrete element method, material density and size of particle are important parameters also for an X-ray tomography experiments. It is commonly assumed that the scanned object of interest should be defined by at least 10 voxels along its minimum principal direction in order to be easily detected and processed. This is also directly correlated with the X-ray tomography spatial resolution in use. Last but not least, mechanical properties such as Young modulus or Poisson's ratio, static friction and internal friction angle are important parameters to consider for flow process control, which, combined with particle shape considerations, can lead to sliding or rolling motion during silo discharging [24,25].
The present project focuses on a multi-parameter approach to select materials that fulfill requirements regarding mechanical, physical and geometrical properties of grain beddings for X-ray tomography experiments and numerical simulations based on the Discrete Element Method. The study was performed for glass beads, glass grits, clay granules, and rice, and sorghum seeds, which represent granular materials of non-biological and biological origin, respectively. While glass beads, glass grits and clay granules are commonly used in the building industry, rice and sorghum are important materials for agriculture and food industry. A detailed analysis of properties of selected particulate materials can provide some hint for researchers to select research material for numerical simulations and X-ray tomography experiments.
The paper is organized as follows: the next section presents the different experimental set-ups used for investigation of material properties. Comparative results are analyzed and discussed in Sect. 3 before concluding remarks.

Research material
The present study has considered a comparison of five granular materials, which differ in origin, shape and chemical structure, but exhibit some similarities in size. Glass bead and grit, representing materials of non-biological origin, are made of the same silica glass type, i.e. soda-potassium glass with barium additive (SiO 2 : 0.694-Na 2 O: 0.104-K 2 O: 0.063-BaO: 0.031-B 2 O 3 : 0.108). Monodispersed spherical glass marbles with average diameter of 3 mm has been chosen, while the glass grits had irregular shapes and size range similar to glass spheres. These materials are often considered for friction/shear studies [26][27][28], packing studies [29,30], gravitational flow experiments [31][32][33] and numerical mod-eling [6,28,31,33,34]. Among the materials of non-biological origin, Seramis was also examined, which is composed of highly porous, processed clay granules. Seramis is composed of kaolinite (Al 2 Si 2 O 5 (OH) 4 ), illite (K,H 3 O) and quartz, which has very high water retention qualities.
The second type of research material included granular plant materials: rice and sorghum seeds. These granular media has also been considered in many research projects, since numerous studies of silo design and discharge are related to agricultural applications [8,[35][36][37][38][39]. Short grain white rice and sorghum has been selected because of a size range corresponding to size of glass bead.

Physical property measurements
Different physical properties, such as density, Young Modulus, Poisson's ratio and coefficients of friction have been measured for single grain and bulk material, which determine behavior of granular materials. A method for determination of densities of glass beads, glass grit and seeds was based on the Archimedes' principle, which is one of the most commonly recommended methods for measuring densities of materials of irregular size and shape. A measurement of density of Seramis was conducted for samples placed in cylindrical chamber of volume of146 cm 3 , by using a Helium pycnometer. Determination of density of Seramis was performed according to ISO recommendations. The elastic moduli and Poisson's ratios of single grains were measured with a material testing machine LLOYD Instruments, as recommended by ASAE Standards [40]. The densities of consolidated materials, and mechanical and frictional parameters of bulk materials were measured during uniaxial compression testing, according to Eurocode 1 recommendations [41]. Measurements were made for samples placed in cylindrical chamber of diameter of 6 cm and height of 2.55 cm. Additionally, frictional properties of selected materials were determined in direct shear test. A Jenike shear tester was used to measure particle-particle and particle-wall friction coefficients for the granular materials and Plexiglas sheet. Plexiglas material is frequently used for X-ray tomography in situ tests, because of its low X-ray attenuation coefficient, high mechanical strength to sustain load or normal pressure and abrasion resistance.

X-ray computed tomography
The task of X-ray computed tomography (XCT) is to reconstruct an object cross section that is related to a so-called attenuation map. The X-ray beam passing through the layer of matter is weakened by dissipation and absorption. It involves irradiating a sample with a small X-ray source and acquiring a series of radiographs while the object is rotated through 360 • in equally spaced steps. A radiograph corresponds to a map of the transmitted intensities. This is detected by an imaging ADC device with a linear response, and is related to the integral attenuation of the different phases composing the object. For a specific energy and for a given pixel sensitivity of the detector, the intensity attenuation is given by the Beer-Lambert law as follows: where I is the beam intensity after passing through the medium; I 0 is the intensity of the incident beam; μ(l) stands for the attenuation coefficient of the examined material as a function of the X-ray path l, i.e. from the source to the considered pixel. After the acquisition of the series of radiographs, a standard filtered back projection algorithm is usually used to reconstruct the 3D map of the attenuation coefficient [42,43]. The value of this coefficient is related to the nature and the density of the absorbing material, its atomic number as well as the selected X-ray energy [44]. As a result, for single phase material and a given energy, the linear attenuation coefficient is proportional to the local density of the material.
In this particular study, the granular material samples were scanned using the Phoenix v|tome|x s240 X-ray computed scanner, which is a versatile high-resolution system for 2D X-ray inspection, 3D failure and structure analyzes. A view of the equipment is shown in Fig. 1. Table 1 summarizes the different experimental conditions used to scan the five tested granular materials. The X-ray energy for glass materials was larger than the one applied for the other materials, because of larger density of glass. The X-ray energy level was also chosen for an X-ray transmission value of about 20% along the longest dimension of the cross-sectional area for all scanned granular materials.
The data were collected and processed using Phoenix datos|x 2 acq software, but for the reconstruction, Phoenix datos|x 2 rec has been used. The reconstruction strategy is based on the well-known filtered-back projection algo-  [43]. This reconstruction platform also includes a wide variety of modules for optimizing CT results for superior precision and quality, especially a beam hardening module. In the present case, a beam hardening correction factor of 7.2 was used for each tested materials. Reconstructed data were further processed and analyzed using Matlab R2013a and Amira 5.2. Note that the radiographs and reconstructed volumes of all tested materials presented in this paper are made available under ODC Attribution Licence [45].

Electron microscope observations
Secondary Electron Microscopy (SEM) images of various granular materials at different magnifications are shown in Fig. 2 to reveal surface characteristics. While the surface of the glass bead is perfectly smooth (despite observations of small cavities), glass grit has a sharp edges and smooth surfaces in the length-width plane (Fig. 2b). Highly magnified images of rice and sorghum reveal a presence of an oily layer at their surface. White rice grain is characterized by compact structure, without any clear visible capillaries and dents. In the case of sorghum, longitudinal recesses centered radially at the pole of spheroid grains are visible (Fig. 2d). The width of the recesses is 10 μm with a maximum depth of 5 μm. Seramis has an irregular shape with structure similar to pumice (Fig. 2e). A high number of pore spaces with diameters ranging from 5 to 20 μm can be observed on the surface of this material.

Physical properties
The following analysis concerns physical properties of both, single grain and bulk material. Glass bead is taken as the reference material for the comparisons.
Results for single grains of the five tested granular materials are shown in Table 2. The results show that glass bead and glass grit have similar densities; however, density of the latter one is smaller by about 3%. Seramis grain has the highest density, which is about 20% higher than reference material. This is associated with its chemical composition, which contains potassium-based illite that is denser than the other silicatebased chemical components. Rice and sorghum have similar densities, which are about 50% smaller than the one measured for reference material. As far as particle-wall friction coefficient μ s,w is concerned, among all tested materials, the smallest value of parameter was obtained for glass grit. In turn, the largest particle-wall friction coefficient was obtained for Seramis (i.e. above 30% larger than for glass), which was due to its rough and porous structure. It can be noticed that rice and sorghum have similar friction coefficients and elastic properties, which are, respectively, around 20 and 25% larger than the ones determined for reference material. The μ s,w value for glass grit is smaller than one measured for glass marble and for natural materials, despite its highly irregular shape and sharp edges. This can be explained by its flake-like shape and opposite planes with smooth surface (see Fig. 2b), which slide down easier on Plexiglas plane during friction tests. An analysis of elastic properties of examined materials has shown that Seramis had the smallest Young modulus, which was due to its porous nature. It is worth highlighting that the Young modulus of glass grit was about 74% smaller than the one determined for the reference material, which can be attributed to its highly non-isotropic shape and heterogeneous structure of glass grit resulting from cracks (see Fig. 2b). Figure 3 and Table 3 present physical/mechanical properties for the bulk materials. It can be observed in Fig. 3 that all bulk density profiles increase with increasing consolidation stress, but at different rates between 0 and 100 kPa of consolidation stress: it corresponds to about 0.55 for glass materials, around 0.3 for biological materials (white rice and sorghum) and about 0.1 for Seramis. In the case of glass grit, the bulk density is smaller w.r.t. glass bead by about 10%. White rice and sorghum have similar bulk density values, which is more than two times larger than the one measured for Seramis. The different values of compressibility ratio s b , which is the relative difference between tapped and poured densities, were also obtained for various materials. In the case of bulk solids composed of irregularly shaped grains, such as glass grit and Seramis, s b equals to 5% and is almost twofold higher than value of the compressibility of the glass bead. In turn, the organic materials (i.e. white rice and sorghum) are characterized by a much greater compressibility ratio, equaled to  Standard deviation values are shown in brackets () *Density measured with helium pycnometer 8.6%. A smaller porosity was obtained for isotropic or nearly isotropic materials like glass beads and sorghum, as compared to materials composed of irregularly shaped grains, such as glass grit or white rice. These results corroborate to findings of Cho et al. [46] who reported that irregularity hinders particle mobility and their ability to attain dense packing configuration. In this study, the porosity of granular packing comprising Seramis was not determined, due to difficulties resulting from porous nature of Seramis grains. A comparison of frictional properties of examined materials has also provided interesting results, which confirmed an effect of cohesiveness and shape of particles on macroscopic friction coefficient in granular assemblies [47][48][49]. Glass grit, which has the largest natural angle of repose Φ b , is characterized by an internal friction angle Φ b,i only 10% larger than the one determined for reference material. This is probably linked with the morphology of the grits that have polished sides. Interestingly, the glass grit also exhibits an increase  Standard deviation values are shown in brackets () also observed that sorghum has the second smallest angle of repose after glass bead and the smallest Φ b,i and Φ b,w values, which are, respectively, 18 and 14% lower than the one determined for glass bead. These small values of angle of repose and friction angles are due to oily film that has been observed in SEM images, combined with the spheroid shape of the grains. Despite its fragility, Seramis has the largest friction angles and the second largest angle of repose, which is due to its porous rough surface and irregular shape. It also has the largest cohesion value (> 2.5 kPa), while glass bead, rice and sorghum can be considered as cohesionless materials. Regarding mechanical properties, differences between results obtained for single grain and bulk materials can be also observed. For instance, a Young modulus of Seramis grain is about 95% lower than the one determined for a single glass marble. However, at the level of the bulk material, this difference is reduced to about 12%. Interestingly, the Young modulus E b of grain assembly composed of glass grit is larger by 12% than the one calculated for glass bead, despite being smaller by 75% at the level of grain. Regarding the natural materials, elastic modulus of grain bedding composed of white rice is similar to one determined for the reference material (< 10%), while E b value for sorghum is smaller by 34%. Simultaneously, single grains of these materials have similar elastic properties.

Image contrast evaluation
The first part of the X-ray tomography analysis concerns the suitability of the selected granular materials for further quantitative studies related to concentration changes, individual particle and packing characterizations. Besides a particle size that is in agreement with the spatial resolution of the X-ray device, a good contrast between the granular material and the surrounding air is a required condition to facilitate material segmentation. Therefore, the first basic analysis concerns intensity histogram of the reconstructed raw volumes, before any image enhancement such as equalization or filtering are applied to increase image contrast. Figure 4a shows that all materials present two strongly separated peaks that correspond to air and granular materials. Seramis presents the case where the peaks are the closest, while glass grit corresponds to the opposite case. This is related to values of bulk density shown in Table 3, but also to the use of stronger Xray energies for the two glass materials as compared to the other granular materials. Sorghum and rice have very similar histograms due to their similar chemical composition. However, in the case of glass materials, different histograms were obtained, since glass grit has its intensity peak shifted to the right, indicating a more pronounced X-ray absorption as compared to glass bead. This goes somehow against the fact that both materials have theoretically the same chemical composition and similar density, as presented in Table 2. Table 2 shows that density of glass grit is smaller by 3% than density of glass marble. It is expected that a further analysis, involving X-ray photoelectron spectroscopy, could bring some clarifications. However, this difference in X-ray absorption is of importance to quantitative analysis, as it is discussed hereafter in Sect. 3.4. Usually, images are quantified after contrast has been enhanced. As mentioned previously, this can be simply done by applying histogram equalization, as it is illustrated in Fig. 4b. The contrast level in the image can be easily calculated from the histogram, using the following expression: where I min and I max are the modal values for the population of the lowest and the largest intensities, respectively, as Fig. 4 Intensity histogram of the five tested materials a from raw data, b after 8-bit conversion and intensity equalization (color figure online) illustrated in Fig. 4b for glass grit. The comparison of histograms obtained for five different granular materials shows that the largest contrast is found for white rice and sorghum (C = 0.82) while the lowest one is observed for Seramis (C = 0.58). Moreover, a larger contrast value is found for glass bead (C = 0.76) than for glass grit (C = 0.65). Cross-sections of the different tested materials are shown in Fig. 5, after volume intensity ranges have been equalized and rescaled to 8-bit (i.e. intensity value between 0 and 255). One can notice the difference in shape of the granular materials, that is spheroid for glass bead and sorghum, ellipsoid for white rice and irregular shape for glass grit and Seramis.
A first qualitative analysis indicates that all materials have similar cross-sectional size range (i.e. ∼ 3 mm), except Seramis that has larger grains. Also, Fig. 5f shows an example of granular mixture of glass bead and glass grit (4% in mass) to illustrate the difference of contrast between the two materials discussed above (Fig. 4b). The images also show contours of grains in red after manual thresholding have been applied for segmentation. No particular advance segmentation method was needed due to the good contrast revealed by the intensity histograms. However, the particle packing induces strong contact between grains in the X-ray tomography images and their individual study needs a digital separation approach to be applied. The watershed algorithm has been applied numerous times for X-ray tomography images [29,[50][51][52][53][54] of granular materials and proved to be an efficient method to separate convex monodispersed objects, but limited for more irregular particulate systems and usually associated to other post-processing method. For instance, the combination with level-set method [51] and contour probability map (stochastic watershed) [54] are good alternatives to refine shape and make extraction more reliable. However, both proposed method are time consuming and would show performance limitation for large volume data.
In the present case, authors have chosen classical markerbased watershed method [55], but with a pre-processing emphasis on selecting good island markers. For that, the following strategy has been chosen: • local region minima, which are related to catchment basin revealed by inversed Euclidean distance map (see Fig. 6a) for each voxel of the segmented granular objects to the "air" background, are extracted and a histogram of their distance is computed, as shown in Fig. 6b. • The largest bin distance corresponding to the frequency that is 75% of the histogram mode is used as a threshold distance value, as shown in Fig. 6b, below which corresponding voxels are considered as markers and set with the same inversed distance value, i.e. the one of the histogram mode.
This procedure creates large patches or island markers, as shown in Fig. 6c, which facilitate the region growing stage of the watershed approach. One can judge on Fig. 6d-e how the method improves the separation of the sorghum grains when compared with the original approach. Typical 3D rendering of the different granular materials are also shown in Fig. 7 to support the overall applicability of the method. Inspection of the data shows that around 10% of the grains are oversegmented, if one excludes particles that touch the region of interest (ROI) borders. Work is in progress to improve the algorithm by including feature classification in order to isolate over-segmented particles, which need further processing steps.

Particle shape analysis
After separation of all grains in the image processing step, a shape analysis can be made in next stage. The main purpose of this particle shape analysis is the estimation of the goodness of the particles of the examined materials to be modelled by ellipsoidal shape (referred hereafter as the ellipsoidal shape goodness, γ ). This is evaluated using the following expression: where V p is the volume of the particle and V c is its so-called coffin box. We understand the notion of the coffin box as the tightest parallelepiped box that surrounds the particle. In eq. (3), a perfect ellipsoid particle will have γ = 1 since V c = 6V p /π . The coffin box is obtained by calculating the matrix of inertia of a particle and getting its eigenvectors that correspond to its principal axes. Along each principal direction the voxel boundaries are searched, which have the largest negative and positive vector rejections [m i ; M i ] i={x,y,z} with respect to the particle centroid C p , as schematized in Fig. 8a. Exemplar results for the five different tested granular materials are shown in Fig. 8b, where a good alignment of the coffin boxes along the main axis of the corresponding particles can be observed. The example of the coffin box of the glass bead is obviously not a matter of discussion, because of the isotropic nature of the particle. The ellipsoidal shape goodness distribution for the five different tested granular materials is depicted in Fig. 9a. It is perfectly noticeable that the three materials which present the largest ellipsoidal shape goodness values are glass grit ( γ = 0.9), sorghum ( γ = 0.86) and white rice ( γ = 0.81), while the lowest goodness value was obtained for glass grit ( γ = 0.6). Valuable, complementary parameters allowing for the comparison of different granular materials are size and aspect ratio. Regarding the size, the smallest dimension of the coffin box, being its thickness, is considered, since it is the one that needs to be compared with voxel size for the sake of image quantification. The histogram shown in Fig. 9b indicates that sorghum presents the average size, which is the closest to monodispersed glass bead. Besides, its size scatter is small and it does not exceed 0.3 mm. Rice and glass grit have the particularity to present the smallest average thickness (35-50% smaller than glass bead diameter), having very large aspect ratio. Figure 9 shows a bivariate histogram of aspect ratio F ab = b/a and F ac = c/a, where a, b and c are the thickness, width and length of particle coffin box, respectively. A strong anisotropy of rice grain shape can be  Fig. 9c, where F ac ≈ 2 F ab = 3.1 with a relatively narrow scatter for F ab ( F ab = 0.14). The values are more scattered for glass grit, with the presence of three different shape populations, indicated by the red arrows on the histogram projection of F ac superimposed in Fig. 9b. The selected elongated glass grit particle, shown in Fig. 8b, is a part of the smallest third population. Interestingly, we can see that glass bead and sorghum present the least scattered shape, with the latter presenting an oblate spheroid shape (i.e. F ac ≈ F ab = 1.6) (Fig. 10).

Discussion: granular material selection
The topic of selection of granular material may look controversial from the point of view of material handling and conveying, since it is rather the container or silo design that is usually determined based on the given material properties and selected flow type. Jenike charts or DEM simulations are usually the principal tools applied for determination of mechanical characteristics of granular materials. However, when the numerical technic is considered, simulations are based not only on experimental data when they can be calculated (for instance macroscopic/microscopic friction properties μ s,w , Φ b,i , Φ b,w ), but also mechanical contact model that includes different forces reflecting sliding fric-tion or rolling friction. Recent models have considered the cases where particle rotation is allowed or not allowed to reflect the cases of spherical or non-spherical particle, the latter having the tendency to only rely on sliding to relieve stress under shear [47]. In other words, the lack of particle rotation would approximate more closely the movement of non-spherical particles, which have less rotational freedom due to mechanical interlocking. Knowing the particle-wall or particle-particle friction coefficients, one can estimate the macroscopic wall or internal friction angle, assuming that particles will or will not undergo rotation during gravitational movement. This is illustrated in Fig. 11 which shows how "no rotation" and "with rotation" relations can be used to predict wall friction angle Φ b,w when particle-wall friction coefficient μ s,w and particle shape are known. However, one can also confront this relationship graph with experimental measurements, as the ones presented in Sect. 3.2. While the results for glass bead, sorghum (oblate shape, nearlyspherical) or glass grit and Seramis (angular shape grains) are in agreement with the predictions, white rice, despite its anisotropic shape, perfectly lies on the "with rotation" curve. Therefore, the predicted ability of sorghum and white rice to roll during gravitational flow becomes an interesting topic to be investigated using quantitative X-ray tomography due to their non-spherical shape. Moreover, not only their X-ray tomography images present very good contrast for boundary detection, but these two materials, as analysed in Sect. 3.3.2, also present high γ values, which make them adapted for DEM modelling with ellipsoids, even if the control of the contact model and its computation during silo discharging is  10 Bivariate histogram of aspect ratio parameters F ab and F ac for the tested granular materials. a Glass bead, b glass grit, c white rice and d sorghum. In the case of (b), the histogram projection for F ac is also displayed (color figure online) more complex and more time consuming than for spherical bead [56].
This multi-parameter approach also suggests that, because of its low Young modulus E s (fragile, tendency to erode during friction), large cohesion c (cohesive material due to rough surface) and lowest average ellipsoidal shape goodness γ (irregular shape that is difficult to model), Seramis is the granular material that is the least adapted to the future comparison between time-lapse X-ray tomography data and results of DEM analysis during silo discharging process of cohensionless material. The cohesiveness of the material during gravitational flow in silo has been recently observed, among other things, using the ultra-fast X-ray tomography system ROFEX [57]. The spherical shape of glass bead is also a drawback to investigate their rolling properties during flow. However, the fact that glass grit and glass bead present Fig. 11 Comparison of experimental wall friction angle and particlewall coefficients of the five tested granular materials with DEM prediction of particle-wall interaction considering both freely rotating and non-rotating spheres [47] (color figure online) a different attenuation in X-ray images makes the former a potential candidate for tracking purposes. Indeed, by injecting a small amount of glass grit with low aspect ratio, one can anticipate that the movements of these tracer particles will be representative of the surrounding media, i.e. glass bead. This is obviously a hypothesis that remains to be verified in future.

Conclusion
The following conclusions can be drawn from the present study: • Granular materials of different origin: glass bead, glass grit, clay granules and short grain white rice, and sorghum, have been examined taking into account physical property measurements and X-ray tomographic imaging. • The granular materials of biological origin, i.e. short grain white rice and sorghum, presented similar material and mechanical properties, e.g. density, compressibility, friction properties and elastic properties, which compete with properties of classical glass materials. The macroscopic and microscopic frictional properties of plant materials were found to fit well with DEM prediction of gravitational flow including rolling friction component. • All tested materials are good candidates for X-ray tomography investigations, based on their relative low density and large grain size, which are adapted for most of the laboratory set-ups available on the market. However, the materials of biological origin are the ones which present the best contrast ratio. This is an important parameter to consider for image segmentation strategies, such as the watershed method, which aims at separating each grain for quantitative analysis. • The quantitative shape analysis, which has mainly been based on the matrix of inertia of the extracted particles, has focused on studying particle aspect ratios and ellipsoidal shape modelling. This quantification has confirmed that, besides glass bead, rice and sorghum can also be satisfyingly modelled by ellipsoid. • All the above-mentioned consideration put rice and sorghum as good research material for DEM simulations and time-lapse X-ray tomography imaging during gravitational flow studies such as silo discharge.