Validation of a fluid-structure interaction model of solute transport in pores of cyclically deformed tissue scaffolds

Convection induced by repetitive compression of porous tissue scaffolds enhances solute transport inside the scaffold. Our previous experiments have shown that pore size, shape, and orientation with respect to strain direction greatly influence loading-induced solute transport. The objective of this study was to develop a computational model of deformation-induced solute transport in porous tissue scaffolds, which included the pore geometry of the scaffold. This geometry consisted of a cubic scaffold with single channel in the middle of the scaffold, immersed in a fluid reservoir. Cylindrical pores with circular or elliptic cross section, and spheroid pores were modeled. The scaffold was cyclically compressed from one side, causing fluid motion and dispersion of solute inside the scaffold pore. Scaffold deformation was solved using the finite element method, and fluid flow and solute transport were solved using the finite volume method. The distortion of the scaffold-fluid interface was transferred as a boundary condition to the fluid flow solver. Both convection and diffusion were included in the computations. The solute transport rates in the different scaffold pore geometries agreed well with our previous experimental results obtained with X-ray microimaging. This model will be used to explore transport properties of a spectrum of novel scaffold designs.