Flow rates in perfusion bioreactors to maximise mineralisation in bone tissue engineering in vitro

In bone tissue engineering experiments, fluid-induced shear stress is able to stimulate cells to produce mineralised extracellular matrix (ECM). The application of shear stress on seeded cells can for example be achieved through bioreactors that perfuse medium through porous scaffolds. The generated mechanical environment (i.e. wall shear stress: WSS) within the scaffolds is complex due to the complexity of scaffold geometry. This complexity has so far prevented setting an optimal loading (i.e. flow rate) of the bioreactor to achieve an optimal distribution of WSS for stimulating cells to produce mineralised ECM. In this study, we demonstrate an approach combining computational fluid dynamics (CFD) and mechano-regulation theory to optimise flow rates of a perfusion bioreactor and various scaffold geometries (i.e. pore shape, porosity and pore diameter) in order to maximise shear stress induced mineralisation. The optimal flow rates, under which the highest fraction of scaffold surface area is subjected to a wall shear stress that induces mineralisation, are mainly dependent on the scaffold geometries. Nevertheless, the variation range of such optimal flow rates are within 0.5–5 mL/min (or in terms of fluid velocity: 0.166–1.66 mm/s), among different scaffolds. This approach can facilitate the determination of scaffold-dependent flow rates for bone tissue engineering experiments in vitro, avoiding performing a series of trial and error experiments.     

Enhancing cell seeding of scaffolds in tissue engineering through manipulation of hydrodynamic parameters

The seeding of cells onto biocompatible scaffolds is a determinant step in the attainment of functional properties of engineered tissues. Efficient, fast and spatially uniform cell seeding can improve the clinical potential of engineered tissue templates. One way to approach these cell seeding requirements is through bioreactor design. In the present study, bovine chondrocytes were seeded (2.5, 5.0 or 10.0 million cells per scaffold) onto polyglycolic acid scaffolds within the hydrodynamic environments of wavy-walled and spinner flask bioreactors. Previous characterizations of the hydrodynamic environment in the vicinity of constructs cultivated in these bioreactors suggested decreased flow-induced shear stress as well as increased recirculation and magnitude of the axial fluid velocities in the wavy-walled bioreactor. Here we report more efficient and spatially uniform cell seeding in the wavy-walled bioreactor, and at intermediate initial cell densities (5 million cells per scaffold). This study constitutes an important step towards the achievement of functional tissue-engineered implants by (i) increasing our understanding of the influence of hydrodynamic parameters on the efficiency and spatial distribution of cell attachment to scaffolds and the production of extracellular matrix and (ii) introducing a comprehensive approach to the investigation of the effects of bioprocessing conditions on tissue morphology and composition.     

Dynamics of diffusivity and pressure drop in flow-through and parallel-flow bioreactors during tissue regeneration

In this study, transport characteristics in flow-through and parallel-flow bioreactors used in tissue engineering were simulated using computational fluid dynamics. To study nutrient distribution and consumption by smooth muscle cells colonizing the 100 mm diameter and 2-mm thick scaffold, effective diffusivity of glucose was experimentally determined using a two-chambered setup. Three different concentrations of chitosan-gelatin scaffolds were prepared by freezing at -80°C followed by lyophilization. Experiments were performed in both bioreactors to measure pressure drop at different flow rates. At low flow rates, experimental results were in agreement with the simulation results for both bioreactors. However, increase in flow rate beyond 5 mL/min in flow-through bioreactor showed channeling at the circumference resulting in lower pressure drop relative to simulation results. The Peclet number inside the scaffold indicated nutrient distribution within the flow-through bioreactor to be convection-dependent, whereas the parallel-flow bioreactor was diffusion-dependent. Three alternative design modifications to the parallel-flow were made by (i) introducing an additional inlet and an outlet, (ii) changing channel position, and (iii) changing the hold-up volume. Simulation studies were performed to assess the effect of scaffold thickness, cell densities, and permeability. These new designs improved nutrient distribution for 2 mm scaffolds; however, parallel-flow configuration was found to be unsuitable for scaffolds more than 4-mm thick, especially at low porosities as tissues regenerate. Furthermore, operable flow rate in flow-through bioreactors is constrained by the mechanical strength of the scaffold. In summary, this study showed limitations and differences between flow-through and parallel-flow bioreactors used in tissue engineering.     

A prediction of cell differentiation and proliferation within a collagen–glycosaminoglycan scaffold subjected to mechanical strain and perfusive fluid flow

Mesenchymal stem cell (MSC) differentiation can be influenced by biophysical stimuli imparted by the host scaffold. Yet, causal relationships linking scaffold strain magnitudes and inlet fluid velocities to specific cell responses are thus far underdeveloped. This investigation attempted to simulate cell responses in a collagen–glycosaminoglycan (CG) scaffold within a bioreactor. CG scaffold deformation was simulated using μ-computed tomography (CT) and an in-house finite element solver (FEEBE/linear). Similarly, the internal fluid velocities were simulated using the afore-mentioned μCT dataset with a computational fluid dynamics solver (ANSYS/CFX). From the ensuing cell-level mechanics, albeit octahedral shear strain or fluid velocity, the proliferation and differentiation of the representative cells were predicted from deterministic functions. Cell proliferation patterns concurred with previous experiments. MSC differentiation was dependent on the level of CG scaffold strain and the inlet fluid velocity. Furthermore, MSC differentiation patterns indicated that specific combinations of scaffold strains and inlet fluid flows cause phenotype assemblies dominated by single cell types. Further to typical laboratory procedures, this predictive methodology demonstrated loading-specific differentiation lineages and proliferation patterns. It is hoped these results will enhance in-vitro tissue engineering procedures by providing a platform from which the scaffold loading applications can be tailored to suit the desired tissue.     

3-D computational modeling of media flow through scaffolds in a perfusion bioreactor

Media perfusion bioreactor systems have been developed to improve mass transport throughout three-dimensional (3-D) tissue-engineered constructs cultured in vitro. In addition to enhancing the exchange of nutrients and wastes, these systems simultaneously deliver flow-mediated shear stresses to cells seeded within the constructs. Local shear stresses are a function of media flow rate and dynamic viscosity, bioreactor configuration, and porous scaffold microarchitecture. We have used the Lattice-Boltzmann method to simulate the flow conditions within perfused cell-seeded cylindrical scaffolds. Microcomputed tomography imaging was used to define the scaffold microarchitecture for the simulations, which produce a 3-D fluid velocity field throughout the scaffold porosity. Shear stresses were estimated at various media flow rates by multiplying the symmetric part of the gradient of the velocity field by the dynamic viscosity of the cell culture media. The shear stress algorithm was validated by modeling flow between infinite parallel plates and comparing the calculated shear stress distribution to the analytical solution. Relating the simulation results to perfusion experiments, an average surface shear stress of 5x10(-5)Pa was found to correspond to increased cell proliferation, while higher shear stresses were associated with upregulation of bone marker genes. This modeling approach can be used to compare results obtained for different perfusion bioreactor systems or different scaffold microarchitectures and may allow specific shear stresses to be determined that optimize the amount, type, or distribution of in vitro tissue growth.     

In silico multi‐scale model of transport and dynamic seeding in a bone tissue engineering perfusion bioreactor

Computer simulations can potentially be used to design, predict, and inform properties for tissue engineering perfusion bioreactors. In this work, we investigate the flow properties that result from a particular poly‐L ‐lactide porous scaffold and a particular choice of perfusion bioreactor vessel design used in bone tissue engineering. We also propose a model to investigate the dynamic seeding properties such as the homogeneity (or lack of) of the cellular distribution within the scaffold of the perfusion bioreactor: a pre‐requisite for the subsequent successful uniform growth of a viable bone tissue engineered construct. Flows inside geometrically complex scaffolds have been investigated previously and results shown at these pore scales. Here, it is our aim to show accurately that through the use of modern high performance computers that the bioreactor device scale that encloses a scaffold can affect the flows and stresses within the pores throughout the scaffold which has implications for bioreactor design, control, and use. Central to this work is that the boundary conditions are derived from micro computed tomography scans of both a device chamber and scaffold in order to avoid generalizations and uncertainties. Dynamic seeding methods have also been shown to provide certain advantages over static seeding methods. We propose here a novel coupled model for dynamic seeding accounting for flow, species mass transport and cell advection‐diffusion‐attachment tuned for bone tissue engineering. The model highlights the timescale differences between different species suggesting that traditional homogeneous porous flow models of transport must be applied with caution to perfusion bioreactors. Our in silico data illustrate the extent to which these experiments have the potential to contribute to future design and development of large‐scale bioreactors. Biotechnol. Bioeng. 2013; 110: 1221–1230. © 2012 Wiley Periodicals, Inc.     

Influence of flow rate and scaffold pore size on cell behavior during mechanical stimulation in a flow perfusion bioreactor

Mechanically stimulating cell-seeded scaffolds by flow-perfusion is one approach utilized for developing clinically applicable bone graft substitutes. A key challenge is determining the magnitude of stimuli to apply that enhances cell differentiation but minimizes cell detachment from the scaffold. In this study, we employed a combined computational modeling and experimental approach to examine how the scaffold mean pore size influences cell attachment morphology and subsequently impacts upon cell deformation and detachment when subjected to fluid-flow. Cell detachment from osteoblast-seeded collagen-GAG scaffolds was evaluated experimentally across a range of scaffold pore sizes subjected to different flow rates and exposure times in a perfusion bioreactor. Cell detachment was found to be proportional to flow rate and inversely proportional to pore size. Using this data, a theoretical model was derived that accurately predicted cell detachment as a function of mean shear stress, mean pore size, and time. Computational modeling of cell deformation in response to fluid flow showed the percentage of cells exceeding a critical threshold of deformation correlated with cell detachment experimentally and the majority of these cells were of a bridging morphology (cells stretched across pores). These findings will help researchers optimize the mean pore size of scaffolds and perfusion bioreactor operating conditions to manage cell detachment when mechanically simulating cells via flow perfusion.     

Flow modelling within a scaffold under the influence of uni-axial and bi-axial bioreactor rotation

The problem of donor scarcity has led to the recent development of tissue engineering technologies, which aim to create implantable tissue equivalents for clinical transplantation. These replacement tissues are being realised through the use of biodegradable polymer scaffolds; temporary/permanent substrates, which facilitate cell attachment, proliferation, retention and differentiated tissue function. To optimise gas transfer and nutrient delivery, as well as to mimic the fluid dynamic environment present within the body, a dynamic system might be chosen. Experiments have shown that dynamic systems enhance tissue growth, with the aid of scaffolds, as compared to static culture systems. Very often, tissue growth within scaffolds is only seen to occur at the periphery. The present study utilises the Computational Fluid Dynamics package FLUENT, to provide a better understanding of the flow phenomena in scaffolds, within our novel bioreactor system. The uni-axial and bi-axial rotational schemes are studied and compared, based on a vessel rotating speed of 35 rpm. The wall shear stresses within and without the constructs are also studied. Findings show that bi-axial rotation of the vessel results in manifold increases of fluid velocity within the constructs, relative to uni-axial rotation about the X- and Z-axes, respectively.     

Epithelial Cell Repopulation and Preparation of Rodent Extracellular Matrix Scaffolds for Renal Tissue Development

This protocol details the generation of acellular, yet biofunctional, renal extracellular matrix (ECM) scaffolds that are useful as small-scale model substrates for organ-scale tissue development. Sprague Dawley rat kidneys are cannulated by inserting a catheter into the renal artery and perfused with a series of low-concentration detergents (Triton X-100 and sodium dodecyl sulfate (SDS)) over 26 hr to derive intact, whole-kidney scaffolds with intact perfusable vasculature, glomeruli, and renal tubules. Following decellularization, the renal scaffold is placed inside a custom-designed perfusion bioreactor vessel, and the catheterized renal artery is connected to a perfusion circuit consisting of: a peristaltic pump; tubing; and optional probes for pH, dissolved oxygen, and pressure. After sterilizing the scaffold with peracetic acid and ethanol, and balancing the pH (7.4), the kidney scaffold is prepared for seeding via perfusion of culture medium within a large-capacity incubator maintained at 37 °C and 5% CO₂. Forty million renal cortical tubular epithelial (RCTE) cells are injected through the renal artery, and rapidly perfused through the scaffold under high flow (25 ml/min) and pressure (~230 mmHg) for 15 min before reducing the flow to a physiological rate (4 ml/min). RCTE cells primarily populate the tubular ECM niche within the renal cortex, proliferate, and form tubular epithelial structures over seven days of perfusion culture. A 44 µM resazurin solution in culture medium is perfused through the kidney for 1 hr during medium exchanges to provide a fluorometric, redox-based metabolic assessment of cell viability and proliferation during tubulogenesis. The kidney perfusion bioreactor permits non-invasive sampling of medium for biochemical assessment, and multiple inlet ports allow alternative retrograde seeding through the renal vein or ureter. These protocols can be used to recellularize kidney scaffolds with a variety of cell types, including vascular endothelial, tubular epithelial, and stromal fibroblasts, for rapid evaluation within this system.     

Design and validation of a dynamic flow perfusion bioreactor for use with compliant tissue engineering scaffolds

In tissue engineering, flow perfusion bioreactors can be used to enhance nutrient diffusion while mechanically stimulating cells to increase matrix production. The goal of this study was to design and validate a dynamic flow perfusion bioreactor for use with compliant scaffolds. Using a non-permanent staining technique, scaffold perfusion was verified for flow rates of 0.1-2.0 mL/min. Flow analysis revealed that steady, pulsatile and oscillatory flow profiles were effectively transferred from the pump to the scaffold. Compared to static culture, bioreactor culture of osteoblast-seeded collagen-GAG scaffolds led to a 27-34% decrease in cell number but stimulated an 800-1200% increase in the production of prostaglandin E(2), an early-stage bone formation marker. This validated flow perfusion bioreactor provides the basis for optimisation of bioreactor culture in tissue engineering applications.     

Recellularization of Decellularized Lung Scaffolds Is Enhanced by Dynamic Suspension Culture

Strategies are needed to improve repopulation of decellularized lung scaffolds with stromal and functional epithelial cells. We demonstrate that decellularized mouse lungs recellularized in a dynamic low fluid shear suspension bioreactor, termed the rotating wall vessel (RWV), contained more cells with decreased apoptosis, increased proliferation and enhanced levels of total RNA compared to static recellularization conditions. These results were observed with two relevant mouse cell types: bone marrow-derived mesenchymal stromal (stem) cells (MSCs) and alveolar type II cells (C10). In addition, MSCs cultured in decellularized lungs under static but not bioreactor conditions formed multilayered aggregates. Gene expression and immunohistochemical analyses suggested differentiation of MSCs into collagen I-producing fibroblast-like cells in the bioreactor, indicating enhanced potential for remodeling of the decellularized scaffold matrix. In conclusion, dynamic suspension culture is promising for enhancing repopulation of decellularized lungs, and could contribute to remodeling the extracellular matrix of the scaffolds with subsequent effects on differentiation and functionality of inoculated cells.     

Numerical simulation of fluid field and in vitro three-dimensional fabrication of tissue-engineered bones in a rotating bioreactor and in vivo implantation for repairing segmental bone defects

In this paper, two-dimensional flow field simulation was conducted to determine shear stresses and velocity profiles for bone tissue engineering in a rotating wall vessel bioreactor (RWVB). In addition, in vitro three-dimensional fabrication of tissue-engineered bones was carried out in optimized bioreactor conditions, and in vivo implantation using fabricated bones was performed for segmental bone defects of Zelanian rabbits. The distribution of dynamic pressure, total pressure, shear stress, and velocity within the culture chamber was calculated for different scaffold locations. According to the simulation results, the dynamic pressure, velocity, and shear stress around the surface of cell-scaffold construction periodically changed at different locations of the RWVB, which could result in periodical stress stimulation for fabricated tissue constructs. However, overall shear stresses were relatively low, and the fluid velocities were uniform in the bioreactor. Our in vitro experiments showed that the number of cells cultured in the RWVB was five times higher than those cultured in a T-flask. The tissue-engineered bones grew very well in the RWVB. This study demonstrates that stress stimulation in an RWVB can be beneficial for cell/bio-derived bone constructs fabricated in an RWVB, with an application for repairing segmental bone defects.     

Biomaterial Scaffolds with Biomimetic Fluidic Channels for Hepatocyte Culture

Biomaterial scaffolds play an important role in maintaining the viability and biological functions of highly metabolic hepatocytes in liver tissue engineering. One of the major challenges involves building a complex microchannel network inside three-dimensional (3D) scaffolds for efficient mass transportation. Here we presented a biomimetic strategy to generate a microchannel network within porous biomaterial scaffolds by mimicking the vascular tree of rat liver. The typical parameters of the blood vessels were incorporated into the biomimetic design of the microchannel network such as branching angle and diameter. Silk fibroin-gelatin scaffolds with biomimetic vascular tree were fabricated by combining micromolding, freeze drying and 3D rolling techniques. The relationship between the micro-channeled design and flow pattern was revealed by a flow experiment, which indicated that the scaffolds with biomimetic vascular tree exhibited unique capability in improving mass transportation inside the 3D scaffold. The 3D scaffolds, preseeded with primary hepatocytes, were dynamically cultured in a bioreactor system. The results confirmed that the pre-designed biomimetic microchannel network facilitated the generation and expansion of hepatocytes.     

3D Bio-Plotted Composite Scaffold Made of Collagen Treated Hydroxyapatite-Tricalciumphosphate for Rabbit Tibia Bone Regeneration

Biphasic calcium phosphate scaffolds with 20/80 HA/TCP ratio were fabricated using the 3D-Bioplotting system to heal critical size defects in rabbit tibia bone. Four different architectures were printed in a layer by layer fashion with lay down patterns viz. (a) 0°– 90°, (b) 0°– 45°– 90°– 135°, (c) 0°–108°– 216° and (d) 0°– 60°– 120°. After high-temperature sintering scaffolds were coated with collagen and were further characterized by (FTIR) Fourier Transform Infrared Spectroscopy, (SEM) Scanning Electron Microscopy, (XRD) X-Ray diffraction, Porosity analysis and Mechanical testing. Scaffold samples were tested for its ability to induce cytotoxicity in Balb/c 3T3 cells at in vitro condition using elution method. Skin sensitization potential of scaffolds was evaluated in male guinea pigs using guinea pig maximization test (GPMT). Further, scaffolds were implanted in eight rabbit tibia bones and biocompatibility and histological evaluations were carried out after 4 and 8 weeks implantation periods. In-vitro results include bonding, surface morphology, phases, porosity, mechanical strength and Cytotoxicity. In-vivo results include sensitization, capsule formation, inflammation, presence of polymorphonuclear cells, giant cells, plasma cells, X-Rays and degradation of the material. It was concluded that HA/TCP/Collagen scaffold with 0°– 45°– 90°– 135° architecture exhibits the most excellent properties in healing critical size bone defects in rabbits.     

Human osteoblast-like cells in three-dimensional culture with fluid flow

In a previous study, we showed that the combination of appropriately designed three-dimensional (3D) microcarrier scaffolds and fluid flow through and around the scaffolds during high aspect ratio vessel (HARV) rotation enhances the elaboration of mineralized bone matrix by osteoblast-like cells. In this study, we describe the ongoing characterization of our 3D culture system, including the investigation of interior fluid flow within the scaffolds and early stage integrin expression during hydrodynamic culture. Using theoretical and experimental methods, we have estimated that cells cultured on the interior of microcarrier scaffolds experience an interior nutrient flow velocity between 1 x 10(-3) and 1 x 10(-2) cm/s and maximum shear stress of 0.03 N/m(2). Under these conditions, osteoblast-like cells grew extensively in the interior regions of the scaffold and retained their osteoblastic phenotype as measured by alkaline phosphatase. In addition, flow cytometric analysis of the overall cell population showed that cells constitutively expressed integrin alpha3beta1 during 3D hydrodynamic culture.     

Flow mixing enhancement from balloon pulsations in an intravenous oxygenator

Computational investigations of flow mixing and oxygen transfer characteristics in an intravenous membrane oxygenator (IMO) are performed by direct numerical simulations of the conservation of mass, momentum, and species equations. Three-dimensional computational models are developed to investigate flow-mixing and oxygen-transfer characteristics for stationary and pulsating balloons, using the spectral element method. For a stationary balloon, the effect of the fiber placement within the fiber bundle and the number of fiber rings is investigated. In a pulsating balloon, the flow mixing characteristics are determined and the oxygen transfer rate is evaluated. For a stationary balloon, numerical simulations show two well-defined flow patterns that depend on the region of the IMO device. Successive increases of the Reynolds number raise the longitudinal velocity without creating secondary flow. This characteristic is not affected by staggered or non-staggered fiber placement within the fiber bundle. For a pulsating balloon, the flow mixing is enhanced by generating a three-dimensional time-dependent flow characterized by oscillatory radial, pulsatile longitudinal, and both oscillatory and random tangential velocities. This three-dimensional flow increases the flow mixing due to an active time-dependent secondary flow, particularly around the fibers. Analytical models show the fiber bundle placement effect on the pressure gradient and flow pattern. The oxygen transport from the fiber surface to the mean flow is due to a dominant radial diffusion mechanism, for the stationary balloon. The oxygen transfer rate reaches an asymptotic behavior at relatively low Reynolds numbers. For a pulsating balloon, the time-dependent oxygen-concentration field resembles the oscillatory and wavy nature of the time-dependent flow. Sherwood number evaluations demonstrate that balloon pulsations enhance the oxygen transfer rate, even for smaller flow rates.     

Modeling and design of optimal flow perfusion bioreactors for tissue engineering applications

Perfusion bioreactors have been used in different tissue engineering applications because of their consistent distribution of nutrients and flow-induced shear stress within the tissue-engineering scaffold. A widely used configuration uses a scaffold with a circular cross-section enclosed within a cylindrical chamber and inlet and outlet pipes which are connected to the chamber on either side through which media is continuously circulated. However, fluid-flow experiments and simulations have shown that the majority of the flow perfuses through the center. This pattern creates stagnant zones in the peripheral regions as well as in those of high flow rate near the inlet and outlet. This non-uniformity of flow and shear stress, owing to a circular design, results in limited cell proliferation and differentiation in these areas. The focus of this communication is to design an optimized perfusion system using computational fluid dynamics as a mathematical tool to overcome the time-consuming trial and error experimental method. We compared the flow within a circular and a rectangular bioreactor system. Flow simulations within the rectangular bioreactor are shown to overcome the limitations in the circular design. This communication challenges the circular cross-section bioreactor configuration paradigm and provides proof of the advantages of the new design over the existing one.     

可控多孔结构支架制备及灌注式体外培养

利用CAD和快速成形技术设计制造具有可控多孔结构的支架。构建灌注式生物反应器系统,实现氧气和营养物质的大量输送,同时产生一定流体剪应力,调节细胞功能的发挥。根据支架负型结构制造出相应的树脂原型,用磷酸钙骨水泥进行填充烧结,得到与设计相符的多孔支架。接种兔成骨细胞,分别采用静态和灌注式三维动态培养方法,观察不同培养条件下细胞在支架表面以及所构造微管道内的生长情况。试验结果表明,灌注式体外培养方法更有利于细胞在支架微管道内的存活和功能的发挥,此灌注式系统能够改善支架微管道内细胞生存的微环境,增强黏附在支架微管道内细胞的活性,促进细胞进一步的增殖和矿化基质的产生。     

Modeling of flow‐induced shear stress applied on 3D cellular scaffolds: Implications for vascular tissue engineering

Novel tissue‐culture bioreactors employ flow‐induced shear stress as a means of mechanical stimulation of cells. We developed a computational fluid dynamics model of the complex three‐dimensional (3D) microstructure of a porous scaffold incubated in a direct perfusion bioreactor. Our model was designed to predict high shear‐stress values within the physiological range of those naturally sensed by vascular cells (1–10 dyne/cm²), and will thereby provide suitable conditions for vascular tissue‐engineering experiments. The model also accounts for cellular growth, which was designed as an added cell layer grown on all scaffold walls. Five model variants were designed, with geometric differences corresponding to cell‐layer thicknesses of 0, 50, 75, 100, and 125 µm. Four inlet velocities (0.5, 1, 1.5, and 2 cm/s) were applied to each model. Wall shear‐stress distribution and overall pressure drop calculations were then used to characterize the relation between flow rate, shear stress, cell‐layer thickness, and pressure drop. The simulations showed that cellular growth within 3D scaffolds exposes cells to elevated shear stress, with considerably increasing average values in correlation to cell growth and inflow velocity. Our results provide in‐depth analysis of the microdynamic environment of cells cultured within 3D environments, and thus provide advanced control over tissue development in vitro. Biotechnol. Bioeng. 2010; 105: 645–654. © 2009 Wiley Periodicals, Inc.     

Modeling,simulations, and optimization of smooth muscle cell tissue engineering for the production of vascular grafts

The paper presents a transient, continuum, two-phase model of the tissue engineering in fibrous scaffolds, including transport equations for the flowing culture medium, nutrient and cell concentration with transverse and in-plane diffusion and cell migration, a novel feature of local in-plane transport across a phenomenological pore and innovative layer-by-layer cell filling approach. The model is successfully validated for the smooth muscle cell tissue engineering of a vascular graft using crosslinked, electrospun gelatin fiber scaffolds for both static and dynamic cell culture, the latter in a dynamic bioreactor with a rotating shaft on which the tubular scaffold is attached. Parametric studies evaluate the impact of the scaffold microstructure, cell dynamics, oxygen transport, and static or dynamic conditions on the rate and extent of cell proliferation and depth of oxygen accessibility. An optimized scaffold of 75% dry porosity is proposed that can be tissue engineered into a viable and still fully oxygenated graft of the tunica media of the coronary artery within 2 days in the dynamic bioreactor. Such scaffold also matches the mechanical properties of the tunica media of the human coronary artery and the suture retention strength of a saphenous vein, often used as a coronary artery graft.