Arthur Sakellariou

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.

Digital Core Laboratory: Properties of reservoir core derived from 3D images

A facility for digital imaging, visualizing and calculation of reservoir rock properties in three dimensions (3D) is described. The facility includes a high resolution X-ray micro-computed tomography system capable of acquiring 3D images made up of 2000 voxels on core plugs up to 5 cm diameter with resolutions down to 2 μm. Subsets of four sandstone reservoir core plugs (5 mm in diameter) from a single well of a producing gas field are imaged in this study. The four cores exhibit a broad range of pore and grain sizes, porosity, permeability and mineralogy. Computational results made directly on the digitized tomographic images are presented for the pore size distribution, permeability, formation factor, NMR response and drainage capillary pressure. We show that data across a range of porosity can be computed from the suite of 5 mm plugs. Computations of permeability, formation factor and drainage capillary pressure are compared to data from a comprehensive SCAL laboratory study on 70 cores from the same well. The results are in good agreement. Empirical correlations between permeability and other petrophysical parameters are made and compared to common correlations. The results demonstrate the potential to predict petrophysical properties from core material not suited for laboratory testing (e.g., drill cuttings, sidewall core or damaged core) and the feasibility of combining digitized images with numerical calculations to predict properties and derive correlations for individual reservoir rock lithologies.

Developing a virtual materials laboratory

Tomographic imaging can now be routinely performed over three orders of magnitude in length scale with correspondingly high data fidelity. This capability, coupled with the development of advanced computational algorithms for image interpretation, three-dimensional visualization, and structural characterization and computation of physical properties on image data, allows for a new numerical laboratory approach to the study of real complex materials: the Virtual Materials Laboratory. Numerical measurements performed directly on images can, in many cases, be performed with similar accuracy to equivalent laboratory measurements, but also on traditionally intractable materials. These emerging capabilities and their impact on a range of scientific disciplines and industry are explored here.

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.

Imaging honey bee brain anatomy with micro-X-ray-computed tomography

Technologies for imaging in three dimensions are greatly desired by researchers in many biological disciplines. However, when imaging small animals such as invertebrates, the achievement of satisfactory spatial resolution and adequate contrast between tissues often requires the use of expensive and time-consuming procedures. Micro-X-ray-computed tomography (muCT) is a convenient technique which is finding greater use alongside conventional microscopies. Staining with heavy metal salts, such as osmium tetroxide improves imaging in muCT, and allows visualization of the 3D structure of the honey bee brain undistorted within the intact head capsule. We obtained detailed information about the morphology of the different brain compartments and were able to show their orientations, relative to each other, within the head capsule. This technique offers a significant improvement in resolution, time, and expense for the quantitative, three-dimensional analysis of developing bee brain centers. In this article, we introduce a rapid, high-resolution, and inexpensive technique for the three-dimensional visualization of different compartments of the honey bee brain. A detailed discussion of the honey bee brain anatomy is provided, demonstrating that muCT, with osmium staining, can indeed visualise these structures. Hence, our results show that muCT is ideally suited for researchers who are interested in the 3D visualization of small invertebrate brains.

Reliable automatic alignment of tomographic projection data by passive auto-focus.

PURPOSE The authors present a robust algorithm that removes the blurring and double-edge artifacts in high-resolution computed tomography (CT) images that are caused by misaligned scanner components. This alleviates the time-consuming process of physically aligning hardware, which is of particular benefit if components are moved or swapped frequently. METHODS The proposed method uses the experimental data itself for calibration. A parameterized model of the scanner geometry is constructed and the parameters are varied until the sharpest 3D reconstruction is found. The concept is similar to passive auto-focus algorithms of digital optical instruments. The parameters are used to remap the projection data from the physical detector to a virtual aligned detector. This is followed by a standard reconstruction algorithm, namely the Feldkamp algorithm. Feldkamp et al. [J. Opt. Soc. Am. A 1, 612-619 (1984)]. RESULTS An example implementation is given for a rabbit liver specimen that was collected with a circular trajectory. The optimal parameters were determined in less computation time than that for a full reconstruction. The example serves to demonstrate that (a) sharpness is an appropriate measure for projection alignment, (b) our parameterization is sufficient to characterize misalignments for cone-beam CT, and (c) the procedure determines parameter values with sufficient precision to remove the associated artifacts. CONCLUSIONS The algorithm is fully tested and implemented for regular use at The Australian National University micro-CT facility for both circular and helical trajectories. It can in principle be applied to more general imaging geometries and modalities. It is as robust as manual alignment but more precise since we have quantified the effect of misalignment.

Investigation of microstructural features in regenerating bone using micro computed tomography.

We illustrate some of the uses of micro-computed tomography (micro-CT) to study tissue-engineered bone using a micro-CT facility for imaging and visualizing biomaterials in three dimensions (3-D). The micro-CT is capable of acquiring 3D X-ray CT images made up of 20003 voxels on specimens up to 5 cm in extent with resolutions down to 2 μm. This allows the 3-D structure of tissue-engineered materials to be imaged across orders of magnitude in resolution. This capability is used to examine an explanted, tissue-engineered bone material based on a polycaprolactone scaffold and autologous bone marrow cells. Imaging of the tissue-engineered bone at a scale of 1 cm and resolutions of 10 μm allows one to visualize the complex ingrowth of bone into the polymer scaffold. From a theoretical viewpoint the voxel data may also be used to calculate expected mechanical properties of the tissue-engineered implant. These observations illustrate the benefits of tomography over traditional techniques for the characterization of bone morphology and interconnectivity. As the method is nondestructive it can perform a complimentary role to current histomorphometric techniques.

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.

The architecture of cartilage: Elemental maps and scanning transmission ion microscopy/tomography

Articular cartilage is not just a jelly-like cover of the bone within the joints but a highly sophisticated architecture of hydrated macromolecules, collagen fibrils and cartilage cells. Influences on the physiological balance due to age-related or pathological changes can lead to malfunction and subsequently to degradation of the cartilage. Many activities in cartilage research are dealing with the architecture of joint cartilage but have limited access to elemental distributions. Nuclear microscopy is able to yield spatially resolved elemental concentrations, provides density information and can visualise the arrangement of the collagen fibres. The distribution of the cartilage matrix can be deduced from the elemental and density maps. The findings showed a varying content of collagen and proteoglycan between zones of different cell maturation. Zones of higher collagen content are characterised by aligned collagen fibres that can form tubular structures. Recently we focused on STIM tomography to investigate the three dimensional arrangement of the collagen structures.

Micro-CT Facility for Imaging Reservoir Rocks at Pore Scales

A micro-CT facility for imaging, visualizing and calculating sedimentary rock properties in three dimensions (3D) is described. The facility is capable of acquiring 3D Xray CT images made up of 2000 voxels on core plugs up to 5 cm diameter with resolutions down to 2 μm. This allows the 3D pore-space of a rock to be imaged across several orders of magnitude. In parallel with standard microscopic techniques, the spatial distributions of different mineralogies can be identified. We demonstrate the capabilities by imaging a reservoir carbonate core at different resolutions. First, an image of a 4 cm diameter plug is analysed at a resolution of 42 μm. This allows one to deduce the size, shape and spatial distribution of the disconnected vug porosity. Within the imaged volume over 30000 separate vugs are identified and a broad vug size distribution is measured. From higher resolution images (2.5-20 μm) on a 5 mm diameter subset of the core one can measure characteristic (intergranular) pore sizes. The apparent porosity of the core increases with enhanced image resolution. This behaviour implies a continuum of pore sizes exist within the core at these resolutions. Carbonate sediments have been conventionally described by a discrete bior tri-modal pore size distribution; in contrast our analysis exhibits no distinct pore sizes but a broad distribution of pore size spanning over more than two orders of magnitude.