Use of a 3D perfusion bioreactor with osteoblasts and osteoblast/endothelial cell co-cultures to improve tissue-engineered bone

The delivery of oxygen, nutrients, and the removal of waste are essential for cellular survival. Culture systems for 3D bone tissue engineering have addressed this issue by utilizing perfusion flow bioreactors that stimulate osteogenic activity through the delivery of oxygen and nutrients by low-shear fluid flow. It is also well established that bone responds to mechanical stimulation, but may desensitize under continuous loading. While perfusion flow and mechanical stimulation are used to increase cellular survival in vitro, 3D tissue-engineered constructs face additional limitations upon in vivo implantation. As it requires significant amounts of time for vascular infiltration by the host, implants are subject to an increased risk of necrosis. One solution is to introduce tissue-engineered bone that has been pre-vascularized through the co-culture of osteoblasts and endothelial cells on 3D constructs. It is unclear from previous studies: 1) how 3D bone tissue constructs will respond to partitioned mechanical stimulation, 2) how gene expression compares in 2D and in 3D, 3) how co-cultures will affect osteoblast activity, and 4) how perfusion flow will affect co-cultures of osteoblasts and endothelial cells. We have used an integrated approach to address these questions by utilizing mechanical stimulation, perfusion flow, and a co-culture technique to increase the success of 3D bone tissue engineering. We measured gene expression of several osteogenic and angiogenic genes in both 2D and 3D (static culture and mechanical stimulation), as well as in 3D cultures subjected to perfusion flow, mechanical stimulation and partitioned mechanical stimulation. Finally, we co-cultured osteoblasts and endothelial cells on 3D scaffolds and subjected them to long-term incubation in either static culture or under perfusion flow to determine changes in gene expression as well as histological measures of osteogenic and angiogenic activity. We discovered that 2D and 3D osteoblast cultures react differently to shear stress, and that partitioning mechanical stimulation does not affect gene expression in our model. Furthermore, our results suggest that perfusion flow may rescue 3D tissue-engineered constructs from hypoxic-like conditions by reducing hypoxia-specific gene expression and increasing histological indices of both osteogenic and angiogenic activity. Future research to elucidate the mechanisms behind these results may contribute to a more mature bone-like structure that integrates more quickly into host tissue, increasing the potential of bone tissue engineering.

Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Exploratory modelling of grinding pressure within a compressed particle bed

With increasing industry interest in high pressure roll grinding (HPGR) technology, there is a strong incentive for improved understanding of the nature of grinding pressure that exists in the interior of a compressed particle bed. This corresponds to the crushing region of the HPGR. The relationship between applied pressure (stress) to the particle bed and induced pressure (stress) within particles and at contact points between particles is of particular interest. A detailed parametric investigation is beyond the scope of this exploratory paper. However, this exploratory investigation does suggest some interesting behaviour. The compressed particle bed within an 80 turn diameter piston has been modelled using Particle Flow Code for three dimensions. PFC3D is a discrete element code. The total number of simulated particles was 1225 and 2450 for two beds of different thickness. Particle diameters were uniformly distributed between 4 and 4.5 mm. The results of the simulations show that stress intensity within the simulated particle beds and within the observed particles increased with increase of the applied stress. The intensity of the average vertical stress in the selected particles tended to be comparable with the intensity of the pressure applied to the surface of particle bed and was only occasionally higher. However, the stress at contact points between particles could be several times higher. In a real crusher, such high stress amplification at contacts will quickly decrease due to local crushing and a resultant increase the size of the contact area. Therefore, its significance is likely to be relatively small in an industrial context. The modelling results also suggest that failure within the particle bed will progress from the crushing surface towards the depth of the bed. (c) 2006 Elsevier Ltd. All rights reserved.

Gravity spun polycaprolactone fibers for applications in vascular tissue engineering:proliferation and function of human vascular endothelial cells

Poly(ε-caprolactone) (PCL) fibers produced by wet spinning from solutions in acetone under low-shear (gravity-flow) conditions resulted in fiber strength of 8 MPa and stiffness of 0.08 Gpa. Cold drawing to an extension of 500% resulted in an increase in fiber strength to 43 MPa and stiffness to 0.3 GPa. The growth rate of human umbilical vein endothelial cells (HUVECs) (seeded at a density of 5 × 104 cells/mL) on as-spun fibers was consistently lower than that measured on tissue culture plastic (TCP) beyond day 2. Cell proliferation was similar on gelatin-coated fibers and TCP over 7 days and higher by a factor of 1.9 on 500% cold-drawn PCL fibers relative to TCP up to 4 days. Cell growth on PCL fibers exceeded that on Dacron monofilament by at least a factor of 3.7 at 9 days. Scanning electron microscopy revealed formation of a cell layer on samples of cold-drawn and gelatin-coated fibers after 24 hours in culture. Similar levels of ICAM-1 expression by HUVECs attached to PCL fibers and TCP were measured using RT-PCR and flow cytometry, indicative of low levels of immune activation. Retention of a specific function of HUVECs attached to PCL fibers was demonstrated by measuring their immune response to lipopolysaccharide. Levels of ICAM-1 expression increased by approximately 11% in cells attached to PCL fibers and TCP. The high fiber compliance, favorable endothelial cell proliferation rates, and retention of an important immune response of attached HUVECS support the use of gravity spun PCL fibers for three-dimensional scaffold production in vascular tissue engineering. © Mary Ann Liebert, Inc.

Novel biodegradable fibres for applications in tissue engineering and drug delivery

The preparation and characterisation of novel biodegradable polymer fibres for application in tissue engineering and drug delivery are reported. Poly(e-caprolactone) (PCL) fibres were produced by wet spinning from solutions in acetone under low shear (gravity flow) conditions. The tensile strength and stiffness of as-spun fibres were highly dependent on the concentration of the spinning solution. Use of a 6% w/v solution resulted in fibres having strength and stiffness of 1.8 MPa and 0.01 GPa respectively, whereas these values increased to 9.9 MPa and 0.1 GPa when fibres were produced from 20% w/v solutions. Cold drawing to an extension of 500% resulted in further increases in fibre strength (up to 50 MPa) and stiffness (0.3 GPa). Hot drawing to 500% further increased the fibre strength (up to 81 MPa) and stiffness (0.5 GPa). The surface morphology of as-spun fibres was modified, to yield a directional grooved pattern by drying in contact with a mandrel having a machined topography characterised by a peak-peak separation of 91 mm and a peak height of 30 mm. Differential scanning calorimetery (DSC) analysis of as-spun fibres revealed the characteristic melting point of PCL at around 58°C and a % crystallinity of approximately 60%. The biocompatibility of as-spun fibres was assessed using cell culture. The number of attached 3T3 Swiss mouse fibroblasts, C2C12 mouse myoblasts and human umbilical vein endothelial cells (HUVECs) on as-spun, 500% cold drawn, and gelatin coated PCL fibres were observed. The results showed that the fibres promoted cell proliferation for 9 days in cell culture and was slightly lower than on tissue culture plastic. The morphology of all cell lines was assessed on the various PCL fibres using scanning electron microscopy. The cell function of HUVECs growing on the as-spun PCL fibres was evaluated. The ability HUVECs to induce an immune response when stimulated with lipopolysaccaride (LPS) and thereby to increase the amount of cell surface receptors was assessed by flow cytometry and reverse transcription-polymerase chain reaction (RT-PCR). The results showed that PCL fibres did not inhibit this function compared to TCP. As-spun PCL fibres were loaded with 1 % ovine albumin (OVA) powder, 1% OVA nanoparticles and 5% OVA nanoparticles by weight and the protein release was assessed in vitro. PCL fibres loaded with 1 % OVA powder released 70%, 1% OVA nanoparticle released 60% and the 5% OVA nanoparticle released 25% of their protein content over 28 days. These release figures did not alter when the fibres were subjected to lipase enzymatic degradation. The OVA released was examined for structural integrity by SDS-PAGE. This showed that the protein molecular weight was not altered after incorporation into the fibres. The bioactivity of progesterone was assessed following incorporation into PCL fibres. Results showed that the progesterone released had a pronounced effect on MCF-7 breast epithelial cells, inhibiting their proliferation. The PCL fibres display high fibre compliance, a potential for controlling the fibre surface architecture to promote contact guidance effects, favorable proliferation rate of fibroblasts, myoblasts and HUVECs and the ability to release pharmaceuticals. These properties recommended their use for 3-D scaffold production in soft tissue engineering and the fibres could also be exploited for controlled presentation and release of biopharmaceuticals such as growth factors.

Gravity spun polycaprolactone fibres for soft tissue engineering:interaction with fibroblasts and myoblasts in cell culture

Poly(ε-caprolactone) (PCL) fibres were produced by wet spinning from solutions in acetone under low shear (gravity flow) conditions. As-spun PCL fibres exhibited a mean strength and stiffness of 7.9 MPa and 0.1 GPa, respectively and a rough, porous surface morphology. Cold drawing to an extension of 500% resulted in increases in fibre strength (43 MPa) and stiffness (0.3 GPa) and development of an oriented, fibrillar surface texture. The proliferation rate of Swiss 3T3 mouse fibroblasts and C2C12 mouse myoblasts on as-spun, 500% cold-drawn and gelatin-modified PCL fibres was determined in cell culture to provide a basic measure of the biocompatibility of the fibres. Proliferation of both cell types was consistently higher on gelatin-coated fibres relative to as-spun fibres at time points below 7 days. Fibroblast growth rates on cold-drawn PCL fibres exceeded those on as-spun fibres but myoblast proliferation was similar on both substrates. After 1 day in culture, both cell types had spread and coalesced on the fibres to form a cell layer, which conformed closely to the underlying topography. The high fibre compliance combined with a potential for modifying the fibre surface chemistry with cell adhesion molecules and the surface architecture by cold drawing to enhance proliferation of fibroblasts and myoblasts, recommends further investigation of gravity-spun PCL fibres for 3-D scaffold production in soft tissue engineering. © 2005 Elsevier Ltd. All rights reserved.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Rheo-optical characterization of liquid crystalline acetoxypropylcellulose melt undergoing large shear flow and relaxation after flow cessation

The rheological and structural characteristics of acetoxypropylcellulose (APC) nematic melt are studied at shear rates ranging from 10 s(-1) to 1000 s(-1) which are relevant to extrusion based processes. APC shows a monotonic shear thinning behavior over the range of shear rates tested. The negative extrudate-swell shows a minimum when a critical shear rate (gamma) over dot(c) is reached. For shear rates smaller than (gamma) over dot(c), the flow-induced texture consists of two set of bands aligned parallel and normal to the flow direction. At shear rates larger than (gamma) over dot(c), the flow induced texture is reminiscent of a 2 fluids structure. Close to the shearing walls, domains elongated along the flow direction and stacked along the vorticity are imaged with POM, whereas SALS patterns indicate that the bulk of the sheared APC is made of elliptical domains oriented along the vorticity. No full nematic alignment is achieved at the largest shear rate tested. Below (gamma) over dot(c), the stress relaxation is described by a stretched exponential. Above (gamma) over dot(c), the stress relaxation is described by a fast and a slow process. The latter coincides with the growth of normal bands thicknesses, as the APC texture after flow cessation consists of two types of bands with parallel and normal orientations relative to the flow direction. Both bands thicknesses do not depend on the applied shear rate, in contrast to their orientation. (C) 2015 Elsevier Ltd. All rights reserved.

Stochastic resonance in a suspension of magnetic dipoles under shear flow

We show that a magnetic dipole in a shear flow under the action of an oscillating magnetic field displays stochastic resonance in the linear response regime. To this end, we compute the classical quantifiers of stochastic resonance, i.e., the signal to noise ratio, the escape time distribution, and the mean first passage time. We also discuss the limitations and role of the linear response theory in its applications to the theory of stochastic resonance.

Development of fluid conduits in the auriferous shear zones of the Hutti Gold Mine, India: evidence for spatially and temporally heterogeneous fluid flow

The gold mineralization of the Hutti Mine is hosted by nine parallel, N - S trending, steeply dipping, 2 - 10 m wide shear zones, that transect Archaean amphibolites. The shear zones were formed after peak metamorphism during retrograde ductile D, shearing in the lower amphibolite facies. They were reactivated in the lower to mid greenschist facies by brittle-ductile D-3 shearing and intense quartz veining. The development of a S-2-S-3 crenulation cleavage facilitates the discrimination between the two deformation events and contemporaneous alteration and gold mineralization. Ductile D, shearing is associated with a pervasively developed distal chlorite - sed cite alteration assemblage in the outer parts of the shear zones and the proximal biotite-plagioclase alteration in the center of the shear zones. D3 is characterized by development of the inner chlorite-K-feldspar alteration, which forms a centimeter-scale alteration halo surrounding the laminated quartz veins and replaces earlier biotite along S-3. The average size of the laminated vein systems is 30-50 m along strike as well as down-dip and 2-6 m in width. Mass balance calculations suggest strong metasomatic changes for the proximal biotite-plagioclase alteration yielding mass and volume increase of ca. 16% and 12%, respectively. The calculated mass and volume changes of the distal chlorite-sericite alteration (ca. 11%, ca. 8%) are lower. The decrease in 6180 values of the whole rock from around 7.5 parts per thousand for the host rocks to 6-7 parts per thousand for the distal chlorite-sericite and the proximal biotite-plagioclase alteration and around 5 parts per thousand for the inner chlorite-K-feldspar alteration suggests hydrothermal alteration during two-stage deformation and fluid flow. The ductile D-2 deformation in the lower amphibolite facies has provided grain scale porosities by microfracturing. The pervasive, steady-state fluid flow resulted in a disseminated style of gold-sulfide mineralization and a penetrative alteration of the host rocks. Alternating ductile and brittle D3 deformation during lower to mid greenschist facies conditions followed the fault-valve process. Ductile creep in the shear zones resulted in a low permeability environment leading to fluid pressure build-up. Strongly episodic fluid advection and mass transfer was controlled by repeated seismic fracturing during the formation of laminated quartz(-gold) veins. The limitation of quartz veins to the extent of earlier shear zones indicate the importance of preexisting anisotropies for fault-valve action and economic gold mineralization. (C) 2003 Elsevier B.V. All rights reserved.

Pulsatile Blood Flow Simulations in Computed Tomography(CT) Scan-Based and Idealized Geometries of Human Aorta

Blood flow in human aorta is an unsteady and complex phenomenon. The complex patterns are related to the geometrical features like curvature, bends, and branching and pulsatile nature of flow from left ventricle of heart. The aim of this work was to understand the effect of aorta geometry on the flow dynamics. To achieve this, 3D realistic and idealized models of descending aorta were reconstructed from Computed Tomography (CT) images of a female patient. The geometries were reconstructed using medical image processing code. The blood flow in aorta was assumed to be laminar and incompressible and the blood was assumed to be Newtonian fluid. A time dependent pulsatile and parabolic boundary condition was deployed at inlet. Steady and unsteady blood flow simulations were performed in real and idealized geometries of descending aorta using a Finite Volume Method (FVM) code. Analysis of Wall Shear Stress (WSS) distribution, pressure distribution, and axial velocity profiles were carried out in both geometries at steady and unsteady state conditions. The results obtained in thesis work reveal that the idealization of geometry underestimates the values of WSS especially near the region with sudden change of diameter. However, the resultant pressure and velocity in idealized geometry are close to those in real geometry

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.