Computational evaluation of oxygen and shear stress distributions in 3D perfusion culture systems: Macro-scale and micro-structured models

We present a combined macro-scale/micro-scale computational approach to quantify oxygen transport and flow-mediated shear stress to human chondrocytes cultured in three-dimensional scaffolds in a perfusion bioreactor system. A macro-scale model was developed to assess the influence of the bioreactor design and to identify the proper boundary conditions for the micro-scale model. The micro-scale model based on a micro-computed tomography (μCT) reconstruction of a poly(ethylene glycol terephthalate)/poly(butylene terephthalate) (PEGT/PBT) foam scaffold, was developed to assess the influence of the scaffold micro-architecture on local shear stress and oxygen levels within the scaffold pores. Experiments were performed to derive specific oxygen consumption rates for constructs perfused under flow rates of 0.3 and 0.03 ml min−1. While macro-scale and micro-scale models predicted similar average oxygen levels at different depths within the scaffold, μCT models revealed small local oxygen variations within the scaffold micro-architecture. The combined macro-scale/micro-scale approach indicated that 0.3 ml min−1, which subjected 95% of the cells to less than 6.3 mPa shear, would maintain the oxygen supply throughout the scaffold above anoxic levels (>1%), with 99.5% of the scaffold supplied with 8–2% O2. Alternatively, at 0.03 ml min−1, the macro-scale model predicted 6% of the cells would be supplied with 0.5–1% O2, although this region of cells was confined to the periphery of the scaffold. Together with local variations predicted by the micro-scale model, the simulations underline that in the current model system, reducing the flow below 0.03 ml min−1 would likely have dire consequences on cell viability to pronounced regions within the engineered construct. The presented approach provides a sensitive tool to aid efficient bioreactor optimization and scaffold design.