The script CN_Contactsurfacearea.m computes coordination number (CN) and contact surface area of a grain from a 3D micro-CT image of a granular pack.
The image should be in ASCII format that has a single column of all the voxel values, each of which represents either a grain label or a void space as 0. Here in this case, it is imported from AVIZO, a commercial software for micro-CT image processing.
A detailed description of this code and how it works are available in Ahmed and Lebedev (2018). Elastic properties of sands, Part-1: Micro-CT image analysis of grain shapes and their relationship with microstructure, Geophysical Prospecting, https://doi.org/10.1111/1365-2478.12652.
Results from this code are used to characterise elastic properties of unconsolidated granular packs using effective medium models based on Hertzian contact theory (Z. Ahmed, and M. Lebedev, 2018): Elastic properties of sands, Part 2: Implementation of contact based model to determine the elasticity of the grains from ultrasonic measurement, Geophysical Prospecting, under revision).
The file Codes_for_CN.mp4 is a video explaining the algorithm.
Note to Visitors
Please note that only current Industry Partners are able to download available script packages, publications and other content on this site.