Analyzing the pore-size distribution of porous materials, made up of an aggregation of interconnected pores, is a demanding task. Mercury intrusion porosimetry (MIP) is a physical method that intrudes mercury into a sample at increasing pressures to obtain a pore-size histogram. This method has been simulated in-silice with several approaches requiring prior computation of a skeleton. We present a new approach to simulate MIP that does not require skeleton computation. Our method is an iterative process that considers the diameters corresponding to pressures. At each iteration, geometric tests detect throats for the corresponding diameter and a CCL process collects the region invaded by the mercury. Additionally, a new decomposition model called CUDB, is used. This is suitable for computing the throats and performs better with the CCL algorithm than a voxel model. Our approach obtains the pore-size distribution of the porous medium, and the corresponding pore graph.
Bibliographical noteFunding Information:
This paper was partially supported by the National Project TIN2008-02903 of the Spanish Government and by Project 02-1/07 of the Institute of Bioengineering of Catalonia (IBEC) . J. Rodríguez was supported by a grant from the University of Carabobo, Venezuela. I. Cruz was supported by a MAE-AECID grant from the Spanish government. The authors thank the anonymous reviewers whose suggestions and questions gave them challenge and opportunity to greatly improve the paper.
Copyright 2011 Elsevier B.V., All rights reserved.
All Science Journal Classification (ASJC) codes
- Modelling and Simulation
- Geometry and Topology
- Computer Graphics and Computer-Aided Design