Three-dimensional open porous scaffolds are commonly used in tissue engineering (TE) applications to provide an initial template for cell attachment and subsequent cell growth and construct development. The macroscopic geometry of the scaffold is key in determining the kinetics of cell growth and thus in vitro 'tissue' formation. In this study, we developed a computational framework based on the level set methodology to predict curvature-dependent growth of the cell/extracellular matrix domain within TE constructs. Scaffolds with various geometries (hexagonal, square, triangular) and pore sizes (500 and 1,000 [Formula: see text]m) were produced in-house by additive manufacturing, seeded with human periosteum-derived cells and cultured under static conditions for 14 days. Using the projected tissue area as an output measure, the comparison between the experimental and the numerical results demonstrated a good qualitative and quantitative behavior of the framework. The model in its current form is able to provide important spatio-temporal information on final shape and speed of pore-filling of tissue-engineered constructs by cells and extracellular matrix during static culture.

Bone tissue engineering strategies use flow through perfusion bioreactors to apply mechanical stimuli to cells seeded on porous scaffolds. Cells grow on the scaffold surface but also by bridging the scaffold pores leading a fully filled scaffold following the scaffold's geometric characteristics. Current computational fluid dynamic approaches for tissue engineering bioreactor systems have been mostly carried out for empty scaffolds. The effect of 3D cell growth and extracellular matrix formation (termed in this study as neotissue growth), on its surrounding fluid flow field is a challenge yet to be tackled. In this work a combined approach was followed linking curvature driven cell growth to fluid dynamics modeling. The level-set method (LSM) was employed to capture neotissue growth driven by curvature, while the Stokes and Darcy equations, combined in the Brinkman equation, provided information regarding the distribution of the shear stress profile at the neotissue/medium interface and within the neotissue itself during growth. The neotissue was assumed to be micro-porous allowing flow through its structure while at the same time allowing the simulation of complete scaffold filling without numerical convergence issues.The results show a significant difference in the amplitude of shear stress for cells located within the micro-porous neo-tissue or at the neotissue/medium interface, demonstrating the importance of taking along the neotissue in the calculation of the mechanical stimulation of cells during culture.The presented computational framework is used on different scaffold pore geometries demonstrating its potential to be used a design as tool for scaffold architecture taking into account the growing neotissue. This article is protected by copyright. All rights reserved

No abstract