Translator Disclaimer
1 May 2009 Doppler optical coherence tomography imaging of local fluid flow and shear stress within microporous scaffolds
Author Affiliations +
Establishing a relationship between perfusion rate and fluid shear stress in a 3D cell culture environment is an ongoing and challenging task faced by tissue engineers. We explore Doppler optical coherence tomography (DOCT) as a potential imaging tool for in situ monitoring of local fluid flow profiles inside porous chitosan scaffolds. From the measured fluid flow profiles, the fluid shear stresses are evaluated. We examine the localized fluid flow and shear stress within low- and high-porosity chitosan scaffolds, which are subjected to a constant input flow rate of 0.5 ml·min-1. The DOCT results show that the behavior of the fluid flow and shear stress in micropores is strongly dependent on the micropore interconnectivity, porosity, and size of pores within the scaffold. For low-porosity and high-porosity chitosan scaffolds examined, the measured local fluid flow and shear stress varied from micropore to micropore, with a mean shear stress of 0.49±0.3 dyn·cm-2 and 0.38±0.2 dyn·cm-2, respectively. In addition, we show that the scaffold's porosity and interconnectivity can be quantified by combining analyses of the 3D structural and flow images obtained from DOCT.



The main objectives of tissue engineering are to provide autologous bioengineered implants developed from the patient’s own cells and to support the body’s ability to regenerate.1, 2 Most tissue engineering strategies involve the use of 3-D porous constructs (scaffolds) to host cellular proliferation activities and to promote extracellular matrix (ECM) production. During the early stages of development, scaffolds are of critical importance in terms of their capability to support cell growth and to deliver sufficient amount of nutrients throughout the volume of the scaffold. Porosity, pore size, and interconnectivity of scaffold are the relevant key parameters that influence tissue development.3

Contrary to standard 2-D cell culture, diffusion alone can not provide adequate nutrient supply and waste disposal for structures larger than approximately 100μm in thickness.4, 5, 6 Consequently, engineered tissues are often cultured in bioreactors that provide dynamic fluid exchange,7, 8 primarily with perfusion systems. Fluid perfusion through the interconnected pores within microporous scaffolds ensures adequate nutrient exchange to promote survivability and growth of the developing tissue.

An important consequence of the dynamic fluid flow inside the scaffold is to mechanically stimulate cells via fluid shear stress.9 In musculoskeletal tissues, mechanical stimuli generally increase ECM production by the cells.10, 11, 12, 13, 14 This increase in ECM alters the local fluid flow, providing a possible feedback mechanism to alter ECM production, which may lead to a complicated relationship between perfusion flow rate and tissue growth.11 This relationship has been well characterized for 2-D cultures where both simulation and experimental data are available.15, 16 However, the relationship from 2-D cultures cannot be directly translated to the studies of 3-D engineered tissue constructs, where complex 3-D porous scaffolds are used and the construct itself constantly evolves during the tissue development process. Thus, the determination of the relation between input perfusion rate and localized fluid flow inside scaffolds is an ongoing and challenging task faced by the tissue engineering community.

Currently, approximations and numerical simulations are often employed to establish the relationship between input flow rate and fluid shear stress in 3-D scaffolds. Wang and Tarbell’s formula17 has been used to assess mean shear stress within a porous construct by relying on the use of Darcy’s law—a constitutive equation that is phenomenologically derived to describe the flow of a fluid through a porous medium. Similarly, a cylindrical pore model has been used to assess the shear stress according to the ratio of the input flow rate to the scaffold surface area.18 Finite element methods have been implemented for estimating localized fluid flow and shear stress in 3-D porous scaffolds.19 These results, however, apply only to a specific and predesigned scaffold that is used in the simulations. Despite these efforts, there is still lack of experimental data to confirm the actual localized fluid flows and shear stresses within the complex 3-D structures used in tissue engineering.

Optical coherence tomography (OCT) is a noninvasive optical imaging technique20, 21, 22 that allows volumetric imaging of the internal architecture of 3-D specimens with high spatial resolution. OCT is analogous to ultrasound, except that near-infrared (NIR) light is used as opposed to sound, and its imaging contrast is based on optical properties of tissue rather than acoustic properties. Because of the use of optical wavelengths, OCT provides an imaging resolution that is of an order of magnitude higher than that of ultrasound, albeit with a much reduced imaging depth. Notwithstanding that OCT imaging depth is limited to within 3mm for highly scattering biological tissue,22 such an imaging depth is highly suitable for the monitoring of tissue growth in developing engineered tissues, where the growing tissue constructs are often only a few millimeters in thickness. A number of studies have investigated the potential of OCT technology as a noninvasive imaging tool to visualize cell dynamics in tissue models for tissue engineering.23, 24, 25, 26

By combining OCT with Doppler velocimetry, a technique called Doppler optical coherence tomography (DOCT) is possible.27, 28, 29 DOCT is a functional extension of OCT that can be used to measure localized fluid flows in highly scattering medium. This is based on the principle that the Doppler frequency shift of the light that is backscattered from the moving objects within a sample is either added to or subtracted from the probing optical frequency, depending on the flow direction. Standard OCT utilizes only the amplitude of backscattered light as a function of depth within the tissue, whereas DOCT uses the additional phase information of the available, complex OCT signals to monitor the velocities of moving particles in the backscattered spectrum. Thus, DOCT is capable of producing simultaneous images of both the tissue’s architecture and its localized fluid flow.

In this article, we report, for the first time, the use of DOCT as a noninvasive imaging modality to explore 3-D architectural morphology, fluid flow, and shear stress in microporous chitosan scaffold—a biodegradable material that is the partially deacetylated derivative of chitin and is used in tissue engineering for effective cell attachment and growth. This optical imaging technique was developed and utilized in order to gain new insights into how the localized flow, and therefore the local shear stresses, vary between the micropores within 3-D scaffolds under a constant input flow rate. After demonstrating the method on a simple phantom model, we apply DOCT to quantify the porosity, interconnectivity, local fluid flow, and fluid shear stress distribution in two types of chitosan-based scaffolds used in tissue engineering. Last, microchannel scaffolds developed specifically for tendon engineering26, 30 are examined.


Materials and Methods


DOCT System

The schematic diagram of the DOCT system used in this investigation is shown in Fig. 1 . Spectral domain DOCT31, 32 was chosen in the setup because this setup has a system sensitivity of more than 20dB as compared to the time-domain OCT systems,33, 34, 35 thus allowing rapid imaging, which is critical for 3-D volumetric optical imaging. In addition, the phase information of OCT signal is directly accessible in this configuration, which is a sine qua non for DOCT imaging of fluid flow. A detailed explanation of the system setup has been given in Refs. 36, 37. In brief, the system used a superluminescent diode with a central wavelength of 842nm and FWHM bandwidth of 45nm that yielded a measured axial resolution of 8μm in air. The sample light was coupled into a probe, consisting of a pair of X-Y galvanometer scanners and optics to both deliver the probing light onto the sample and collect the light backscattered from the sample. The lateral imaging resolution was approximately 20μm as determined by the objective lens that focused the light onto the sample. The irradiance on the sample surface was 1mW . The detection system was a custom-built high-speed spectrometer that had a designed spectral resolution of 0.055nm , resulting in an imaging depth of 3.2mm in air ( z axis shown in Fig. 1). For imaging of the fluid flow within scaffolds, the camera integration time was set at 100μs for imaging, which meant that the imaging rate was at 10,000 axial scans per second.

Fig. 1

Schematic of the OCT system used to provide the 3-D microstructural and flow maps of the chitosan scaffolds. CCD—charge coupled device; PC—polarization controller. The inset shows the DOCT measurement geometry.


The current system used a phase-resolved DOCT approach to determine the velocity of fluid flow. Details of this approach can be found in Ref. 38. Briefly, the OCT signal in spectral domain OCT is a complex function, of which the magnitude is used to provide the imaging contrast for morphological structures of the sample. However, the phase term ϕ(z) of the OCT signal is generally random along the depth z , but fixed in position for the (relatively) static scatterers. A translation of the scatterers at position z by a distance Δd(z) during the time interval Δt between two successive axial scans will induce a change in the measured phase of the OCT signal given by

Eq. 1

where n0 is the average refractive index of the sample, and k¯ is the average wave number of the light source used. This phase change is related to the Doppler frequency by fD(z)=Δϕ(z)2πΔt . Evaluating Δϕ(z) at each depth z yields depth-resolved measurements of both the magnitude and direction of the flow velocity. Assuming that the probe beam intersects the flow velocity vector at an angle α (Doppler angle), the flow velocity is given by

Eq. 2

The use of phase changes between successive OCT scans to assess the fluid flow requires the OCT signals in interferograms to be correlated.38 To meet this requirement, the interval between successive OCT scans was set to 2.5μm , while the probing beam spot was 20μm (transverse resolution).


Scaffold Fabrication

Two-percent chitosan solution was prepared by dissolving chitosan flakes (95% deacytalated, Sigma-Aldrich) in 1% v/v glacial acetic acid. The resulting solution was cast in cylindrical molds ( 1.5mm inner diameter) that were frozen overnight and subsequently freeze-dried. The degrees of porosity in the scaffolds were controlled by adjusting the freezing rate.39 The scaffolds were sequentially rehydrated in a gradient of ethanol (100%, 70%, 50%) and stored in phosphate buffered saline (PBS) before imaging under DOCT. Three sets of scaffolds were prepared for the current study: (1) low-porosity chitosan scaffolds (LPCS) with a pore size ranging from 30to100μm ; (2) high-porosity scaffolds (HPCS) with a pore size between 100 and 200μm ; and (3) microchannel chitosan scaffolds (MCS) within which a 250-μm -diam microchannel was located in the center surrounded with micropores of size ranging from 30to100μm . Microchannel scaffolds were fabricated according to the protocol published in Ref. 26 by casting the chitosan solution in a mold in which a stainless steel needle was located at its central position.


Perfusion System

The perfusion system consisted of a 1.5-mm -inner-diameter tubing connected to a precision DC pump. The sample chamber, made of a portion of the tubing, was fixed on a goniometer under the DOCT probe at an angle of 8deg from the horizontal to achieve a Doppler angle of 82deg . The system was set up to deliver a constant input flow rate of 0.5ml-min1 . A 0.5% latex microsphere ( 0.3μm in diameter) suspension was used as the light scattering medium to monitor the fluid flow through the scaffolds situated within the tube. Preliminary perfusion tests were performed to verify the flow measurements by use of a scattering flow phantom that was flowing in a 700-μm -inner-diameter transparent plastic tube.


Shear Stress

In a fluid, the tangential stress, τ , can be related to the shear rate by:40

Eq. 3

where μ is the dynamic viscosity of the fluid. For the Dulbecco’s Modified Eagle’s Medium (DMEM) that we used in this study, we assumed that μ=8.2×104Kgm1s1 Ref. 41 Shear rate γ̇ is a measure of the shear deformation rate, which is the divergence of the displacement rate vector at any point (x,y,z) in 3-D space, γ̇=v̇ , where v is the displacement and the dot over the variable explicitly indicates the time derivative. Therefore, the magnitude of shear stress can be directly related to the fluid velocity measured by DOCT:

Eq. 4

where v is the flow velocity vector, vz is one component of velocity in the z direction (probe beam direction) measured by DOCT, n is the normal direction of flow, which is perpendicular to the wall and v . It is clear that to accurately calculate the shear stress, the Doppler angle must be known a priori. This condition can be met in simple geometries, such as pipe flow; however, for the fluid that flows in the complex micropores within the scaffold, it is not a trivial matter because α is not generally known. To deal with this matter, we approximate α as the Doppler angle as if the fluid flows in the bioreactor tube without the scaffold. This approximation, although crude, facilitates comparison of the OCT measurements with the previous published approximation approaches using Darcy’s law17 and the cylindrical pore model.18


Porosity and Interconnectivity Assessment

DOCT provides both the 3-D microstructural morphology and 3-D flow map for each measured scaffold. The combination of the structural and flow images permitted the direct measurement of porosity and interconnectivity for the chitosan scaffolds. The porosity was defined as the ratio of total pore volume to total chitosan volume as determined by the OCT structural images; it was determined by a block image processing.42 Briefly, in block processing for each B-scan, a threshold was set to discriminate between true empty pores and the pore walls. A binarized image was then produced and split into several blocks of equal area. The local porosity in each block, represented by the density of pixels, was mapped, ranging from 100% to 0%. As fluid flow existed only in the interconnected pores, it was used to assess interconnectivity. We defined the interconnectivity as the ratio of the total pore volume, from the OCT structural image (Fig. 1), to the interconnected volume, from the DOCT flow image. To improve the estimation accuracy, the OCT structural images were obtained from the scaffolds situated in the tube with perfusion of the PBS solution, while the OCT flow images were acquired from the same scaffold at the same location with perfusion of the light scattering medium.



Following scanning by DOCT, the chitosan scaffolds were fixed for 10min in 4% glutaraldehyde, embedded in paraffin, and cross-sectioned ( 5μm thick) through the center of the scaffolds. The sections were stained with hematoxylin and eosin (H&E) to show the microstructures of chitosan scaffolds.




Measurements of Scattering Fluid Flow in Phantom

Figure 2a shows a representative structural OCT image of the tube phantom in which scattering fluid was flowing at 0.15mlmin1 controlled by a precision DC pump. The relatively transparent wall of the tubing is seen as the almost black region in the image, as the associated scattering is low compared to the latex particle solution flowing in the center of the tube. The flow velocity image obtained from DOCT and the associated shear stress, respectively, are shown in Figs. 2b and 2c. The flow velocity is at its maximum in the center and decreases with radial distance to the wall. After the Doppler angle correction, the measured maximum velocity in the center of lumen was 42.21±0.02 ( mean±std. dev., n=5 ) mm/s. This result is consistent with that calculated from hydrodynamic theory, the value of which is 44.44mms under the assumption of laminar flow in the tube.43 Since the Reynolds number in this case was significantly less than the critical value of 1800, the assumption of laminar flow in the tube is valid.44 The distribution of the velocity of the fluid flow inside the cylindrical channel has a parabolic profile, as is clearly revealed in Fig. 2d, which corresponds to the region marked as the dashed box in Fig. 2b. The corresponding shear stress was calculated to be approximately zero near the center. The shear stress was a maximum at the wall with a calculated value of 26.88±0.30dyncm2 , while the theoretical prediction was 27.81dyncm2 , using the assumption of laminar flow in the tube. Note that in order to minimize the noise artifacts present in the OCT velocity and shear stress images, a 3×3 median filter was first applied to the velocity image to obtain the subsequent shear stress image.

Fig. 2

Typical cross-sectional DOCT images of the flow characteristics in a standard circular plastic tube: (a) the structural image, (b) the flow map, and (c) the corresponding shear stress map. (d) The 3-D flow velocity plot for the area marked by the dashed box in (b); the inset shows the velocity plot across the channel of the tube, i.e., along a line passing through the center, where it is seen that the flow is parabolic with the maximum in the center and minimum at the walls. The units for color bars shown are (a) dB, (b) mm/s, (c) dyncm2 , and (d) mm/s. The white bars indicate 700μm .



Low-porosity Chitosan Scaffolds

Figure 3a shows a typical cross-sectional image of LPCS acquired by the OCT system. The bidirectional flow velocity map obtained by DOCT is presented in Fig. 3b. Both the distribution and the magnitudes of the flow can be seen. To show the variations of velocities within pores, the magnitudes of bidirectional flow velocities were calculated and plotted in Fig. 3c, where the heterogeneous distribution of flow velocity in the porous structures is clearly revealed. Although the input flow rate was constant, the local fluid flow in this complex construct varied greatly in both magnitude and direction. Furthermore, the flow in the micropores did not show parabolic distributions. Consequently, the fluid shear stress, as shown in Fig. 3d, differed from pore to pore, with values ranging from 0to1.65dyncm2 . Nevertheless, the maximum fluid shear stress was generally located at the pore walls.

Fig. 3

DOCT monitored fluid flow in LPCS in situ. Shown are, respectively, (a) the microstructural image, (b) the bidirectional local fluid flow map, (c) the 3-D plot of magnitudes of the flow velocities, and (d) the corresponding local shear stress distribution in a typical X-Z section. The units for color bars shown are (a) dB, (b) mm/s, (c) mm/s, and (d) dyncm2 . The white bars indicate 200μm .


Figure 4 gives the 3-D volume-rendered image of the fluid flow through the isosurface reconstruction (i.e., reconstructed by enclosing all parts of a volume that are brighter than some user-defined threshold). As no flow is present in the isolated pores, the visualization of interconnected structures within chitosan scaffold can be readily appreciated. From these results, it appears that the pores of LPCS are nearly spherical and highly interconnected.

Fig. 4

3-D volume rendering of flow in the interconnected pores of LPCS, revealing the inner connection of this type of scaffold. The physical image size was 1.0×1.5×1.0(x-y-z)mm3 .



High-porosity Chitosan Scaffolds

With the same imaging and perfusion conditions as for LPCS, Figs. 5a, 5b, 5c, 5d give, respectively, the structural image, the flow velocity map, the flow magnitude plot, and the fluid shear stress of a typical OCT cross section of HPCS. From these figures, it seems that the fluid flow exists mainly in the pores of relatively large size. The magnitude of the flow is inhomogeneously distributed within the pores, resulting in shear stresses ranging from 0to1.07dyncm2 . The 3-D reconstruction of the flow (Fig. 6 ) shows that the shape of pores tends to be ellipsoidal and elongated, and the pores are interconnected in a unidirectional way. The fluid flows preferentially along sequentially aligned pores.

Fig. 5

DOCT monitored fluid flow in HPCS in situ. Shown are, respectively, (a) the microstructural image, (b) the bidirectional local fluid flow map, (c) the 3-D plot of magnitudes of the flow velocities, and (d) the corresponding local shear stress distribution in a typical X-Z section. The units for color bars shown are (a) dB, (b) mm/s, (c) dyncm2 , and (d) mm/s. The white bars indicate 200μm .


Fig. 6

3-D volume rendering of flow in the interconnected pores of HPCS, revealing the inner connection of this type of scaffold. The physical image size was 1.0×1.5×1.0(x-y-z)mm3 .



Shear Stress Distribution within Pores

The distributions of shear stresses at the pore walls of both LPCS and HPCS in 12 random cross-sectional scans are shown in Figs. 7a and 7b, respectively. The shear stresses within LPCS were distributed more broadly than those within HPCS. This indicates that the constant input flow rate imparts a wide range of shear stress within the low-porosity chitosan scaffold. In Table 1, the mean shear stress of both LPCS and HPCS (n=12) at the pore walls are compared with the values from the Wang and Tarbell approach17 and from the cylindrical pore model.18 All the values from DOCT are represented by mean±std. . dev. For a constant input flow rate (0.5mlmin1) , the mean shear stress is significantly different between LPCS and HPCS (unpaired two-tailed t -test, P<0.005 , n=12 ), and with different distributions. The shear stress is higher in the pores of LPCS compared to HPCS. In both cases, the experimental measurements were approximately one order of magnitude smaller than the mean shear stresses yielded by the Wang and Tarbell equation and the cylindrical pore model approximation.

Fig. 7

Wall shear stress distribution in a representative X-Y scan (n=12) for (a) LPCS and (b) HPCS, respectively.


Table 1

Comparison of the wall shear stress obtained from DOCT (n=12) with Wang and Tarbell and cylindrical pore models.

Mean shear stress (dyn⋅cm−2) DOCTWang and Tarbell (Ref. 17)Cylindrical poremodel (Ref. 18)
LPCS 0.49±0.3 8.595.46
HPCS 0.38±0.2 7.372.73


Microchannel Chitosan Scaffold

Figure 8a illustrates an OCT structure image of the microchannel chitosan scaffold (MCS) where the microchannel is seen as a relatively dark volume inside the scattering porous chitosan core. The flow velocity map indicates the fluid preferentially flowed in the microchannel [Fig. 8b]. The maximum measured velocity was 10.92mms1 . Low-velocity flows with magnitudes less than 0.5mms existed in the surrounding pores. However, the flow velocity distribution in the microchannel was not parabolic, as emphasized by a close-up of the ROI in Fig. 8c. Consequently, the fluid shear stress was nonuniformly distributed within the channel, with a non-null value at its center [see Fig. 8d]. The mean shear stress within the microchannel was 0.37±0.15dyncm2 .

Fig. 8

The flow characteristics for a chitosan scaffold with 250-μm central microchannel measured by DOCT, where (a) the microstructural information and (b) the corresponding fluid flow can be visualized in situ. (c) The 3-D velocity plot of the region marked with the dashed box in (b), where the flow is nonparabolic within the microchannel. (d) The corresponding local shear stress distribution within the cross section. The units for color bars shown are (a) dB, (b) mm/s, (c) mm/s, and (d) dyncm2 . The white bars indicate 200μm .




H&E staining of the LPCS and HPCS sections, Figs. 9a and 9b, agree with the structural information acquired by DOCT. The pore size and shape can be observed. The pores in LPCS have sizes ranging from 30to100μm , whereas they are between 100to200μm in HPCS. The pore shapes revealed in the histological section may be inferred from the power spectra of the grayscale images of Figs. 9a and 9b through the 2-D Fourier transformation approach. The results are shown in the insets (upper-left corner) of the individual images, from which it can be seen that the shape of pores in LPCS is nearly spherical (nearly rounded power spectrum), whereas it tends to ellipsoidal in HPCS (elongated power spectrum). In both the cases, the pores are highly interconnected.

Fig. 9

H&E staining of a cross section of (a) LPCS and (b) HPCS. Chitosan is stained red. The inset indicates the power spectrum of the corresponding images performed through the Fourier transformation approach. (Color online only.)



Porosity and Interconnectivity

The porosities of the LPCS and the HPCS calculated from a number of 3-D volumetric DOCT scans (n=6) were 0.85±0.04 and 0.90±0.02 , respectively. These two values were significantly different ( t -test, P<0.005 , n=6 ). Corresponding interconnectivity values were 0.97±0.02 and 0.79±0.02 , respectively, which were also significantly different ( t -test, P<0.005 , n=6 ).



Optical coherence tomography (OCT) has been demonstrated as a potential tool for in situ monitoring of volumetric architecture in tissue engineering constructs, and the addition of Doppler capability provides experimental fluid flow data within porous scaffolds. This information is valuable, as there is a consensus on the advantages of culturing tissue engineered constructs inside perfusion bioreactors.7, 8 Nutrient transport is greatly improved by dynamic fluid flow as compared to diffusion alone for scaffolds thicker than approximately 100μm .4, 5, 6 In this study, we have shown that the interconnectivity of a scaffold may be assessed by DOCT. The flow data can be correlated with the structural image to identify closed pores or regions. We have proposed a measure of the degree of interconnection based on the ratio of the volume where the fluid actually flows to the total volume of the pores. Such measurements may be useful in the optimization of the manufacturing of scaffolds. Our study shows that the interconnectivity varies between LPCS and HPCS scaffolds, with a lower interconnectivity for the latter.

OCT provided 3-D microstructural and flow maps for an imaging depth up to 1.0mm within the highly scattering chitosan scaffolds. The imaging depth is limited by the scattering properties of the material and by the probing wavelength used. It was reported that a probing light with a wavelength longer than 840nm —for example, 1300nm —would enhance the imaging depth when OCT images highly scattering tissue.45, 46 Therefore, increased imaging depth into the chitosan scaffolds might be obtained by using a probing light of 1300-nm wavelength as opposed to 840nm used in this study. By relying on the Doppler effect of the moving particles, theoretically DOCT could measure the fluid flow as small as tens of microns per second, largely limited by the system noise floor.21, 22 However, because of the optical heterogeneity of the sample, the sensitivity of the DOCT system to fluid flow is practically limited by the heterogeneous texture of the sample that appeared as the velocity noise floor in the DOCT flow images.47 For the current study with a 2.5-μm separation between adjacent axial scans, the velocity noise floor was 0.3radians ; thus, the minimum measurable velocity is 0.178mms . Therefore, fluid flows with velocities less than 0.178mms were not detected by the current OCT system. Despite this problem, the measured shear stress using DOCT (26.88dyncm2) from the well-controlled phantom experiment exhibited only a 3% deviation from the theoretical expectation (27.81dyncm2) . This implies that the current system may be sufficient to describe the shear stress within the simple geometries.

However, the fact that the minimal resolvable flow velocity was at 0.178mms would likely lead to an underestimation of the calculated interconnectivity values for the scaffolds studied. This is the case for HPCS, where the flow is predominant in larger pores. We note that if the interconnectivity is the only parameter of interest, there are two possible solutions to correct this underestimation: (1) the flow rate is physically increased, and (2) the minimal resolvable flow velocity is improved through a careful system design or reduction of separation distance between adjacent axial scans during OCT imaging. We believe that a measure of the interconnectivity based on the pores where the nutritive fluid flows at appropriate values is a suitable measure for the optimization of 3-D cell culture.

A constant perfusion rate induces fluid shear stress within a scaffold and therefore on the cells within the scaffold. Studies have shown that this mechanical stimulus could enhance ECM production for bone cells within values ranging from 1.1to24dyncm2 (Refs. 48, 49, 50, 51, 52). Shear stresses higher than the physiological window can decrease cell attachment and viability.30 Therefore, the determination of fluid shear stresses inside porous scaffolds is critical for tissue engineering. So far, several models have been explored in an attempt to establish a relationship between the input flow rate and the velocity within the porous constructs. Darcy’s law has been applied to 3-D collagen models with an assumption of homogeneous pore distribution, but it yields only a mean wall shear stress.17 The cylindrical pore approximation model18 does not reflect the complexity of the porous structures used in tissue engineering. Both models give values of the mean wall shear stress that are one order of magnitude greater than that measured by DOCT in this study. Numerical simulations by finite element methods compute the local fluid and the resulting shear stress flow in each pore according to an input flow rate.19 Although these results have not been yet validated by experimental data, they have shown the importance of the porosity and pore size on the distribution of the fluid shear stresses. This underlines the fact that input flow rate is not an independent parameter governing mechanical stimulation. Unfortunately, optimization of the input flow rate in terms of tissue generation in 3-D scaffolds by use of this model is valid only for its specific scaffold architecture. Furthermore, the current study revealed that the measured shear stress is different between LPCS and HPCS for the same input flow rate. This difference implies that the input flow rate alone might not be used to describe the mechanical stimuli applied to the cells.

DOCT allows real-time in situ measurements of the fluid shear stress at localized positions within a porous scaffold, and the results shown in this study emphasize the need to experimentally measure the fluid velocity within the scaffolds. We found that within the porous scaffolds, the local fluid flow is not uniformly distributed among the pores and there is a large variation in the magnitude of the flow velocities. This results in a variation of the fluid shear stresses from pore to pore, which may lead to a nonhomogeneous cell attachment, cell proliferation, and ECM production. The measured mean shear stress (around 0.5dyncm2 ) for a typical input flow rate of 0.5mlmin1 is an order of magnitude below the values previously reported.17, 18 This result may indicate that cells could be stimulated by even lower levels of shear stress than that stated in previous studies. Once the porosity has decreased due to tissue growth, the fluid shear stress could reach physiological values.

In a previous study,26 it has been postulated that microporous scaffolds with an inner microchannel could be used for tendon tissue engineering. It was shown that this architecture supports cell growth and ECM deposition within the microchannel when cultured in a perfusion bioreactor, leading eventually to collagen bundle formation.26 In the present paper, DOCT results support this postulation, as the flow predominantly exists in the microchannel that would encourage the tissue growth at this location. However the impact of the nonhomogeneous distribution of the fluid shear stress inside the microchannel requires further investigation.

The challenge in optimizing the parameters of a perfusion bioreactor is that there is a feedback loop of tissue growth on the distribution and magnitude of the fluid flow velocity. For a constant input flow rate, the mean shear stress may increase continuously until a nonviable value is reached, resulting in localized cell and ECM detachment and tissue resorption. This, in turn, would decrease the local fluid flow velocity, until equilibrium is reached. Additionally, tissue growth could result in a loss of interconnectivity, by isolating large portions of the scaffolds from nutritive flow. Ultimately, only noninvasive, in situ measurements of both structural information and actual fluid shear stress during tissue development will allow the optimization of perfusion bioreactors for tissue engineering.



We have demonstrated that the localized fluid flow and shear stress within porous scaffolds can be experimentally evaluated by Doppler optical coherence tomograph (DOCT). Three types of porous chitosan scaffolds have been explored by DOCT in order to gain new insight into how the fluid flow, and therefore the shear stress, varies within the micropores when the scaffolds were subjected to a constant input flow rate. With the input flow rate kept at 0.5mlmin1 , we have shown that the behavior of the fluid flow and shear stress in micropores is strongly dependent on the micropore interconnectivity, porosity, and size of pores within the scaffolds. The fluid flows were found to have greater variation within the low-porosity scaffold than in the high-porosity scaffold, leading to a broader distribution of the shear stress for the former, with a mean value of 0.49±0.3 versus 0.38±0.2dyncm2 at the pore walls (Table 1). For the microporous scaffold in which there was a microchannel 250μm in diameter, the fluid flow predominantly existed in the microchannel. In addition, we have also shown that with the DOCT images, a number of key parameters relevant to the scaffold—for example, porosity and interconnectivity—can be quantified.


The work was supported in part by research grants from the National Heart, Lung, and Blood Institute (1R01 HL093140-01), the American Heart Association Grant-in-Aid (0855733G), and the Biotechnological and Biological Sciences Research Council, UK (BBS/B/04242). The content is solely the responsibility of the authors and does not necessarily represent the official views of grant-giving bodies.



R. Langer and J. P. Vacanti, “Tissue engineering,” Science, 260 920 –926 (1993). 0036-8075 Google Scholar


L. G. Griffith and G. Naughton, “Tissue engineering—current challenges and expanding opportunities,” Science, 295 1009 –1014 (2002). 0036-8075 Google Scholar


T. S. Karande, J. L. Ong, and C. M. Agrawal, “Diffusion in musculoskeletal tissue engineering scaffolds: design issues related to porosity, permeability, architecture, and nutrient mixing,” Ann. Biomed. Eng., 32 1728 –1743 (2004). 0090-6964 Google Scholar


L. G. Cima and R. Langer, “Engineering human tissue,” Chem. Eng. Prog., 89 46 –50 (1993). 0360-7275 Google Scholar


C. K. Colton, “Implantable biohybrid artificial organs,” Cell Transplant, 4 415 –436 (1995). 0963-6897 Google Scholar


R. L. Carrier, “Perfusion improves tissue architecture of engineered cardiac muscle,” Tissue Eng., 8 175 –188 (2002). 1076-3279 Google Scholar


I. Martin, D. Wendt, and M. Heberer, “The role of bioreactors in tissue engineering,” Trends Biotechnol., 22 80 –86 (2004). 0167-7799 Google Scholar


L. E. Freed and G. Vunjak-Novakovic, “Tissue engineering bioreactors,” Principles of Tissue Engineering, 143 –156 Academic Press, San Diego (2000). Google Scholar


J. A. Frangos, L. V. Mclntire, and S. G. Eskin, “Shear stress induced stimulation of mammalian cell metabolism,” Biotechnol. Bioeng., 32 1053 –1060 (1988). 0006-3592 Google Scholar


T. Davisson, R. L. Sah, and A. Ratcliffe, “Perfusion increases cell content and matrix synthesis in chondrocyte three-dimensional cultures,” Tissue Eng., 8 807 –816 (2002). 1076-3279 Google Scholar


G. N. Bancroft, V. I. Sikavitsas, J. van den Dolder, T. L. Sheffield, C. G. Ambrose, J. A. Jansen, and A. G. Mikos, “Fluid flow increases mineralized matrix deposition in 3-D perfusion culture of marrow stromal osteoblasts in a dose-dependent manner,” Proc. Natl. Acad. Sci. U.S.A., 99 12600 –12605 (2002). 0027-8424 Google Scholar


V. I. Sikavitsas, G. N. Bancroft, H. L. Holtorf, J. A. Jansen, and A. G. Mikos, “Mineralized matrix deposition by marrow stromal osteoblasts in 3-D perfusion culture increases with increasing fluid shear forces,” Proc. Natl. Acad. Sci. U.S.A., 100 14683 –14688 (2003). 0027-8424 Google Scholar


A. M. Freyria, Y. Yang, H. Chajra, C. F. Rousseau, M. C. Ronzière, D. Herbage, and A. J. El Haj, “Optimization of dynamic culture conditions: effects on biosynthetic activities of chondrocytes grown in collagen sponges,” Tissue Eng., 11 674 –684 (2005). 1076-3279 Google Scholar


Y. Yang and A. J. El Haj, “Bone tissue engineering,” SCOPE, 13 33 –35 (2004). 0271-972X Google Scholar


C. R. Jacobs, C. E. Yellowley, B. R. Davis, Z. Zhou, J. M. Cimbala, and H. J. Donahue, “Differential effect of steady versus oscillating flow on bone cells,” J. Biomech., 31 969 –976 (1998). 0021-9290 Google Scholar


M. Edlich, C. E. Yellowley, C. R. Jacobs, and H. J. Donahue, “Oscillating fluid flow regulates cytosolic calcium concentration in bovine articular chondrocytes,” J. Biomech., 34 59 –65 (2001). 0021-9290 Google Scholar


S. Wang and J. M. Tarbell, “Effect of fluid flow on smooth muscle cells in a 3-dimensional collagen gel model,” Arterioscler., Thromb., Vasc. Biol., 20 2220 –2225 (2000). 1079-5642 Google Scholar


A. S. Goldstein, T. M. Juarez, C. D. Helmke, M. C. Gustin, and A. G. Mikos, “Effect of convection on osteoblastic cell growth and function in biodegradable polymer foam scaffolds,” Biomaterials, 22 1279 –1288 (2001). 0142-9612 Google Scholar


F. Boschetti, M. T. Raimondi, F. Migliavacca, and G. Dubini, “Prediction of the micro-fluid dynamic environment imposed to three-dimensional engineered cell systems in bioreactors,” J. Biomech., 39 418 –425 (2006). 0021-9290 Google Scholar


D. Huang, E. Swanson, C. Lin, J. Schuman, W. Stinson, W. Chang, M. Hee, T. Flotte, K. Gregory, C. Puliafito, and J. Fujimoto, “Optical coherence tomography,” Science, 254 1178 –1181 (1991). 0036-8075 Google Scholar


A. F. Fercher, W. Drexler, C. K. Hitzenberger, and T. Lasser, “Optical coherence tomography—principles and applications,” Rep. Prog. Phys., 66 239 –303 (2003). 0034-4885 Google Scholar


P. H. Tomlins and R. K. Wang, “Theory, development and applications of optical coherence tomography,” J. Phys. D, 38 2519 –2535 (2005). 0022-3727 Google Scholar


X. Xu, R. K. Wang, and A. J. El Haj, “Investigation of changes in optical attenuation of bone and neuronal cells in organ culture or three-dimensional constructs in vitro with optical coherence tomography: relevance to cytochrome oxidase monitoring,” Eur. Biophys. J., 32 355 –362 (2003). 0175-7571 Google Scholar


W. Tan, A. Sendemir-Urkmez, L. J. Fahrner, R. Jamison, D. Leckband, and S. A. Boppart, “Structural and functional optical imaging of three-dimensional engineered tissue development,” Tissue Eng., 10 1747 –1756 (2004). 1076-3279 Google Scholar


Y. Yang, A. Dubois, X. P. Qin, J. Li, A. J. El Haj, and R. K. Wang, “Investigation of optical coherence tomography as an imaging modality in tissue engineering,” Phys. Med. Biol., 51 1649 –1659 (2006). 0031-9155 Google Scholar


P. O. Banganinchi, Y. Yang, N. Maffulli, R. K. Wang, and A. J. El Haj, “Chitosan microchannel scaffolds for tendon tissue engineering characterized by optical coherence tomography,” Tissue Eng., 13 313 –322 (2007). 1076-3279 Google Scholar


Z. Chen, T. Milner, D. Dave, and J. Nelson, “Optical Doppler tomographic imaging of fluid flow velocity in highly scattering media,” Opt. Lett., 22 64 –66 (1997). 0146-9592 Google Scholar


Z. Chen, Y. Zhao, S. Srinivas, J. Nelson, N. Prakash, and R. Frostig, “Optical Doppler tomography,” IEEE J. Sel. Top. Quantum Electron., 5 1134 –1142 (1999). 1077-260X Google Scholar


R. K. Wang, “High-resolution visualization of fluid dynamics with Doppler optical coherence tomography,” Meas. Sci. Technol., 15 725 –733 (2004). 0957-0233 Google Scholar


P. O. Bagnaninchi, Y. Yang, A. J. El Haj, and N. Maffuli, “Tissue engineering for tendon repair,” Br. J. Sports Med., 41 e10 (2007). 0306-3674 Google Scholar


A. F. Fercher, C. K. Hitzenberger, G. Kamp, and S. Y. El-Zaiat, “Measurement of intraocular distances by backscattering spectral interferometry,” Opt. Commun., 117 43 –48 (1995). 0030-4018 Google Scholar


G. Hausler and M. W. Linduer, “Coherence radar and spectral radar—new tools for dermatological diagnosis,” J. Biomed. Opt., 3 21 –31 (1998). 1083-3668 Google Scholar


J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett., 28 2067 –2069 (2003). 0146-9592 Google Scholar


M. A. Choma, M. V. Sarunic, C. Yang, and J. A. Izatt, “Sensitivity advantage of swept source and Fourier domain optical coherence tomography,” Opt. Express, 11 2183 –2189 (2003). 1094-4087 Google Scholar


R. Leitgeb, C. K. Hitzenberger, and A. F. Fercher, “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express, 11 889 –894 (2003). 1094-4087 Google Scholar


R. K. Wang and Z. H. Ma, “A practical approach to eliminate autocorrelation artifacts for volume-rate spectral domain optical coherence tomography,” Phys. Med. Biol., 51 3231 –3239 (2006). 0031-9155 Google Scholar


R. K. Wang, S. Jacques, Z. Ma, S. Hurst, S. Hanson, and A. Gruber, “Three-dimensional optical angiography,” Opt. Express, 15 4083 –4097 (2007). 1094-4087 Google Scholar


Y. Zhao, Z. P. Chen, C. Saxer, S. Xiang, J. F. de Boer, and J. S. Nelson, “Phase-resolved optical coherence tomography and optical Doppler tomography for imaging blood flow in human skin with fast scanning speed and high velocity sensitivity,” Opt. Lett., 25 114 –116 (2000). 0146-9592 Google Scholar


S. V. Madihally and H. W. Matthew, “Porous chitosan scaffolds for tissue engineering,” Biomaterials, 20 1133 –1142 (1999). 0142-9612 Google Scholar


P. M. Gerhart, Fundamentals of Fluid Mechanics, Addison-Wesley, Boston (1992). Google Scholar


M. T. Raimondi, F. Boschetti, L. Falcone, G. B. Fiore, A. Remuzzi, M. Marazzi, E. Marinoni, and R. Pietrabissa, “Mechanobiology of engineered cartilage cultured under a quantified fluid dynamic environment,” Biomech. Model. Mechanobiol., 1 69 –82 (2002). 1617-7959 Google Scholar


Y. Yang, P. Bagnaninchi, M. Wood, A. El Haj, E. Guyot, A. Dubois, and R. K. Wang, “Monitoring cell profile in tissue engineered constructs by OCT,” Proc. SPIE, 5695 51 –57 (2005). 0277-786X Google Scholar


R. K. Wang, “High resolution visualization of fluid dynamics with Doppler optical coherence tomography,” Meas. Sci. Technol., 15 725 –733 (2004). 0957-0233 Google Scholar


L. D. Landau and E. M. Lifshitz, Fluid Mechanics, 106 2nd ed.Pergamon Press, Oxford, UK (1987). Google Scholar


Y. Pan and D. L. Farkas, “Noninvasive imaging of living human skin with dual–wavelength optical coherence tomography in two and three dimensions,” J. Biomed. Opt., 3 446 –465 (1998). 1083-3668 Google Scholar


R. K. Wang and S. Hurst, “Mapping of cerebrovascular blood perfusion in mice with skin and cranium intact by optical micro-angiography at 1300-nm wavelength,” Opt. Express, 15 11402 –11412 (2007). 1094-4087 Google Scholar


R. K. Wang and Z. Ma, “Real-time flow imaging by removing texture pattern artifacts in ultrafast spectral domain optical Doppler tomography,” Opt. Lett., 31 3001 –3003 (2006). 0146-9592 Google Scholar


T. L. Donahue, T. R. Haut, C. E. Yellowley, H. J. Donahue, and C. R. Jacobs, “Mechanosensitivity of bone cells to oscillating fluid flow induced shear stress may be modulated by chemotransport,” J. Biomech., 36 1363 –1371 (2003). 0021-9290 Google Scholar


K. M. Reich and J. A. Frangos, “Effect of flow on prostaglandin E2 and inositol trisphosphate levels in osteoblasts,” Am. J. Physiol., 261 428 –432 (1991). 0002-9513 Google Scholar


J. Klein-Nulend, A. van der Plas, C. M. Semeins, N. E. Ajubi, J. A. Frangos, P. J. Nijweide, and E. H. Burger, “Sensitivity of osteocytes to biomechanical stress in vitro,” FASEB J., 9 441 –445 (1995). 0892-6638 Google Scholar


C. H. Turner and F. M. Pavalko, “Mechanotransduction and function response of the skeleton to physical stress: the mechanisms and mechanics of bone adaptation,” J. Orthop. Sci., 3 346 –355 (1998). 0949-2658 Google Scholar


S. Weinbaum, S. C. Cowin, and Y. A. Zeng, “A model for the excitation of osteocytes by mechanical loading induced bone fluid shear stresses,” J. Biomech., 27 339 –360 (1994). 0021-9290 Google Scholar
©(2009) Society of Photo-Optical Instrumentation Engineers (SPIE)
Yali Jia, Pierre-Olivier Bagnaninchi, Ying Yang, Alicia J. El Haj, Monica T. Hinds, Sean J. Kirkpatrick, and Ruikang Wang "Doppler optical coherence tomography imaging of local fluid flow and shear stress within microporous scaffolds," Journal of Biomedical Optics 14(3), 034014 (1 May 2009).
Published: 1 May 2009

Back to Top