Prediction of permeability of regular scaffolds for skeletal tissue engineering: a combined computational and...

Scaffold permeability is a key parameter combining geometrical features such as pore shape, size and interconnectivity, porosity and specific surface area. It can influence the success of bone tissue engineering scaffolds, by affecting oxygen and nutrient transport, cell seeding efficiency, in vitro three-dimensional (3D) cell culture and, ultimately, the amount of bone formation. An accurate and efficient prediction of scaffold permeability would be highly useful as part of a scaffold design process. The aim of this study was (i) to determine the accuracy of computational fluid dynamics (CFD) models for prediction of the permeability coefficient of three different regular Ti6Al4V scaffolds (each having a different porosity) by comparison with experimentally measured values and (ii) to verify the validity of the semi-empirical Kozeny equation to calculate the permeability analytically. To do so, five CFD geometrical models per scaffold porosity were built, based on different geometrical inputs: either based on the scaffold's computer-aided design (CAD) or derived from 3D microfocus X-ray computed tomography (micro-CT) data of the additive manufactured (AM) scaffolds. For the latter the influence of the spatial image resolution and the image analysis algorithm used to determine the scaffold's architectural features on the predicted permeability was analysed. CFD models based on high-resolution micro-CT images could predict the permeability coefficients of the studied scaffolds: depending on scaffold porosity and image analysis algorithm, relative differences between measured and predicted permeability values were between 2% and 27%. Finally, the analytical Kozeny equation was found to be valid. A linear correlation between the ratio Φ(3)/S(s)(2) and the permeability coefficient k was found for the predicted (by means of CFD) as well as measured values (relative difference of 16.4% between respective Kozeny coefficients), thus resulting in accurate and efficient calculation of the permeability of regular AM scaffolds.