Three-Dimensional Trabecular Alignment Model

The Shear-slip Mesh Update Method (SSMUM) is being used in flow simulations involving large but regular displacements of one or more boundaries of the computational domain. We follow up the earlier discussion of the method with notes on practical implementation aspects. In order to establish a benchmark problem for this class of flow problems, we define and report results from a two-dimensional viscous flow around a rotating stirrer in a square chamber. The application potential of the method is demonstrated in the context of biomedical design problem, as we perform an analysis of blood flow in a centrifugal left ventricular assist device, or blood pump, which involves a rotating impeller in a non-axisymmetric housing.     

Evaluation of laser-Doppler flowmetry as a measure of tissue blood flow

In this study the technique of laser-Doppler flowmetry was evaluated for the measurement of tissue blood flow by comparing laser-Doppler flow (LDF) signal in the renal cortex, gracilis muscle, and cremaster muscle of anesthetized rats to whole-organ blood flow measured with an electromagnetic flowmeter or radioactive microspheres. In vitro, LDF signal was closely correlated (r = 0.99) to changes in erythrocyte velocity generated with a rotating wheel. Although individual LDF readings varied in situ, mean LDF signal calculated from multiple readings on the tissue surface were significantly correlated (r = 0.74-0.95) with tissue blood flows measured at various perfusion pressures. However, significant differences in the slope of the LDF signal vs. blood flow relationship were observed in different tissues and with different methods of measurement in the same tissue. This study indicates that mean laser-Doppler flow signal provides a good estimate of tissue blood flow, provided a sufficient number of points is scanned. However, there appears to be no universal calibration factor for the method.     

Development of a platelet adhesion transport equation for a computational thrombosis model

Thrombosis is a significant issue for cardiovascular device development and use. While thrombosis models are available, very few are device-related and none have been thoroughly validated experimentally. Here, we introduce a surface adherent platelet transport equation into a continuum model to account for the biomaterial interface/blood interaction. Using a rotating disc system and polyurethane-urea material, we characterize steady and pulsatile flow fields using laser Doppler velocimetry. In vitro measurements of platelet adhesion are used in combination with the LDV data to provide further experimental validation. The rotating disc system is computationally studied using the device-induced thrombosis model with the surface platelet adherent transport equation. The results indicate that the flow field is in excellent agreement to the experimental LDV data and that the platelet adhesion simulations are in good agreement with the in vitro platelet data. These results provide good evidence that this transport equation can be used to express the relationship between blood and a biomaterial if the correct platelet adhesion characteristics are known for the biomaterial. Further validation is necessary with other materials.     

No need for particle tracing: From accumulating fluid properties to novel blood coagulation model in the lattice Boltzmann method

The accumulation of a fluid property from the standpoint of a particle moving with non-steady fluid flow (i.e., platelet/blood-cell damage index in pulsating blood flow) is a challenged computational problem due to the current need for particle-tracing methods. The method we developed (dubbed VPI) enables the approximation of the Lagrangian integral in real-time for any point in space and time for the entire domain and which is easily integrated into the the lattice Boltzmann method. As an illustrative numerical example we applied our method to a blood coagulation model which was shown to accurately capture the coagulation characteristics observed in experiments, and therefore opening a door for more detailed study of systems which are currently hard to study using particle tracing methods.     

A novel technology for large vessel recanalization

Large vessels recanalization is a challenge for mechanical atherectomy devices, where the lumen debulked is close in diameter to the device crossing profile, which may be only 25% to 30% of the original lumen size; so, the procedure can restore only a fraction of the original blood flow. Moreover, small diameter lumens are prone to be repeatedly occluded after a relatively short period of time. In this article, we present a novel technology of recanalization, using a catheter-based microjet system to deliver a flux of biocompatible abrasive particles to the lesion site, resulting in microchipping of the plaque, while minimizing trauma to the vessel wall. Plaque debris is removed from the blood flow, and blood flow is restored. In contrast to rotating mechanical devices, plaque debulking can be performed up to diameters that are substantially larger than the device crossing profile, supporting superior long-term patency. As a case study, we evaluated the technology for use in the superficial femoral artery where the lesions tend to be very long and heavily calcified with high restenosis rates.     

In situ collagen gelation: a new method for constructing large tissue in rotary culture vessels

In situ collagen gelation is a method that combines a static three-dimensional culture technique with rotating bioreactors. This method was designed for large dense tissue engineering ex vivo. To challenge the current limitations on size, we combined the static collagen gel embedding method with high-aspect ratio rotating bioreactors. Rat calvarial cells in gelated collagens were cultured in rotating vessels with 5 mM beta-glycerophosphate-containing medium for 1, 2, or 3 wk and then analyzed for cell morphology, cell distribution, and viability, as well as for contraction of the collagen gel. The size of collagen gels with rat calvarial cells averaged 2.8 cm in diameter x 0.25 cm in thickness at the end of 3 wk. Scanning electron microscopy and laser scanning confocal microscopy of collagen gels revealed a homogeneous distribution of living cells. Despite the barrier effects from induced calcification, in collagen gels, cell metabolic activity (alkaline phosphatase assay and 3-[4,5-dimethyl-2-thiazolyl]-2,5-diphenyl-2H-tetrazolium bromide assay) increased over the 3 wk, and cell viability (trypan blue exclusion and flow cytometry analysis) remained at about 90% at the end of 3 wk. Based on our results, we determined that in situ collagen gelation provides a feasible method for engineering large dense tissue ex vivo.     

Analysis of non-linear pulsatile blood flow in arteries

Some aspects of the problem of computation of non-linear pulsatile blood flow in large arteries are investigated, in the context of the computational method developed by Ling and Atabek (1972). As examples, the following aspects are considered: stability of the computations; representation of higher-frequency components of the flow; effects of keeping or omitting non-linear terms in the equations; effects of varying the dimensionless parameters of the problem. The computational method is extended to include effects of viscoelasticity of arterial walls.     

Fluorescent microsphere technique to measure cerebral blood flow in the rat

The present protocol describes a method for parallel measurement of cerebral blood flow (CBF) using fluorescent microspheres and structural assessment of the same material. The method is based on the standard microsphere technique, embolizing capillaries proportional to the blood flow, but requires dissolution of the tissue to retrieve the microspheres. To link the blood flow to the tissue morphology we modified the technique to fluorescent microspheres, which are quantified in cryo- or vibratome sections, allowing structural analysis by, for example, immunohistochemistry or standard histology. The protocol takes 8 h 50 min, without pauses, to complete, but additional flow measurements or specific protocols can increase the time needed.     

Measurement of blood viscosity using mass-detecting sensor

A newly designed mass-detecting capillary viscometer is extended to measure the viscosity of whole blood over a range of shear rates without the use of anticoagulants in a clinical setting. In the present study as proof of principle, a single measurement of liquid-mass variation with time replaces the flow rate and pressure drop measurements that are usually required for the operation of a capillary tube viscometer. Using a load cell and capillary, we measured the change of mass flowing through capillary tube with respect to the time, m(t), from which viscosity and shear rate were mathematically calculated. For water and adulterated bloods, excellent agreement was found between the results from the mass-detecting capillary viscometer and those from a commercially available rotating viscometer. Also, the mass-detecting capillary viscometer measured the viscosity of unadulterated whole blood without heparin or EDTA. This new method overcomes the drawbacks of conventional viscometers in the measurement of the whole blood viscosity. First, the mass-detecting capillary viscometer can accurately and consistently measure the unadulterated blood viscosity over a range of shear rates in less than 2 min without any anticoagulants. Second, this design provides simplicity (i.e. ease of operation, no moving parts, and disposable) and low cost.     

Trypanosome Motion Represents an Adaptation to the Crowded Environment of the Vertebrate Bloodstream

Blood is a remarkable habitat: it is highly viscous, contains a dense packaging of cells and perpetually flows at velocities varying over three orders of magnitude. Only few pathogens endure the harsh physical conditions within the vertebrate bloodstream and prosper despite being constantly attacked by host antibodies. African trypanosomes are strictly extracellular blood parasites, which evade the immune response through a system of antigenic variation and incessant motility. How the flagellates actually swim in blood remains to be elucidated. Here, we show that the mode and dynamics of trypanosome locomotion are a trait of life within a crowded environment. Using high-speed fluorescence microscopy and ordered micro-pillar arrays we show that the parasites mode of motility is adapted to the density of cells in blood. Trypanosomes are pulled forward by the planar beat of the single flagellum. Hydrodynamic flow across the asymmetrically shaped cell body translates into its rotational movement. Importantly, the presence of particles with the shape, size and spacing of blood cells is required and sufficient for trypanosomes to reach maximum forward velocity. If the density of obstacles, however, is further increased to resemble collagen networks or tissue spaces, the parasites reverse their flagellar beat and consequently swim backwards, in this way avoiding getting trapped. In the absence of obstacles, this flagellar beat reversal occurs randomly resulting in irregular waveforms and apparent cell tumbling. Thus, the swimming behavior of trypanosomes is a surprising example of micro-adaptation to life at low Reynolds numbers. For a precise physical interpretation, we compare our high-resolution microscopic data to results from a simulation technique that combines the method of multi-particle collision dynamics with a triangulated surface model. The simulation produces a rotating cell body and a helical swimming path, providing a functioning simulation method for a microorganism with a complex swimming strategy.     

Three-dimensional rotational plasma flows near solid surfaces in an axial magnetic field

A rotational flow of a conducting viscous medium near an extended dielectric disk in a uniform axial magnetic field is analyzed in the magnetohydrodynamic (MHD) approach. An analytical solution to the system of nonlinear differential MHD equations of motion in the boundary layer for the general case of different rotation velocities of the disk and medium is obtained using a modified Slezkin–Targ method. A particular case of a medium rotating near a stationary disk imitating the end surface of a laboratory device is considered. The characteristics of a hydrodynamic flow near the disk surface are calculated within the model of a finite-thickness boundary layer. The influence of the magnetic field on the intensity of the secondary flow is studied. Calculations are performed for a weakly ionized dense plasma flow without allowance for the Hall effect and plasma compressibility. An MHD flow in a rotating cylinder bounded from above by a retarding cap is considered. The results obtained can be used to estimate the influence of the end surfaces on the main azimuthal flow, as well as the intensities of circulating flows in various devices with rotating plasmas, in particular, in plasma centrifuges and laboratory devices designed to study instabilities of rotating plasmas.     

Analysis of a high‐throughput cone‐and‐plate apparatus for the application of defined spatiotemporal flow to cultured cells

The shear stresses derived from blood flow regulate many aspects of vascular and immunobiology. In vitro studies on the shear stress‐mediated mechanobiology of endothelial cells have been carried out using systems analogous to the cone‐and‐plate viscometer in which a rotating, low‐angle cone applies fluid shear stress to cells grown on an underlying, flat culture surface. We recently developed a device that could perform high‐throughput studies on shear‐mediated mechanobiology through the rotation of cone‐tipped shafts in a standard 96‐well culture plate. Here, we present a model of the three‐dimensional flow within the culture wells with a rotating, cone‐tipped shaft. Using this model we examined the effects of modifying the design parameters of the system to allow the device to create a variety of flow profiles. We first examined the case of steady‐state flow with the shaft rotating at constant angular velocity. By varying the angular velocity and distance of the cone from the underlying plate we were able to create flow profiles with controlled shear stress gradients in the radial direction within the plate. These findings indicate that both linear and non‐linear spatial distributions in shear stress can be created across the bottom of the culture plate. In the transition and “parallel shaft” regions of the system, the angular velocities needed to provide high levels of physiological shear stress (5 Pa) created intermediate Reynolds number Taylor‐Couette flow. In some cases, this led to the development of a flow regime in which stable helical vortices were created within the well. We also examined the system under oscillatory and pulsatile motion of the shaft and demonstrated minimal time lag between the rotation of the cone and the shear stress on the cell culture surface. Biotechnol. Bioeng. 2013; 110: 1782–1793. © 2013 Wiley Periodicals, Inc.     

An approach to automatic blood vessel image registration of microcirculation for blood flow analysis on nude mice

Image registration is often a required and a time-consuming step in blood flow analysis of large microscopic video sequences in vivo. In order to obtain stable images for blood flow analysis, frame-to-frame image matching as a preprocessing step is a solution to the problem of movement during image acquisition. In this paper, microscopic system analysis without fluorescent labelling is performed to provide precise and continuous quantitative data of blood flow rate in individual microvessels of nude mice. The performance properties of several matching metrics are evaluated through simulated image registrations. An automatic image registration programme based on Powell's optimisation search method with low calculation redundancy was implemented. The matching method by variance of ratio is computationally efficient and improves the registration robustness and accuracy in practical application of microcirculation registration. The presented registration method shows acceptable results in close requisition to analyse red blood cell velocities, confirming the scientific potential of the system in blood flow analysis.     

Finite element analysis of the lift on a slightly deformable and freely rotating and translating cylinder in two-dimensional channel flow

Motivated by the lateral migration phenomena of fresh and glutaraldehyde-fixed red blood cells in a field flow fractionation (FFF) separation system, we studied the transverse hydrodynamic lift on a slightly flexible cylinder in a two-dimensional channel flow. The finite element method was used to analyze the flow field with the cylinder at different transverse locations in the channel. The shape of the cylinder was determined by the pressure on the surface of the cylinder from the flow field solution and by the internal elastic stress. The cylinder deformation and the flow field were solved simultaneously. The transverse lift exerted on the cylinder was then calculated. The axial and angular speed of the cylinder were iterated such that the drag and torque on the cylinder were nulled to represent a freely translating and rotating state. The results showed that the transverse lift on a deformable cylinder increased greatly and the equilibrium position moved closer to the center of the channel compared to a rigid cylinder. Also, with the same elastic modulus but a higher flow rate, a larger deformation and higher equilibrium location were found. The maximum deformation of the cylinder occurred when the cylinder was closest to the wall where a larger shear rate existed. The numerical results and experimental studies are discussed.     

In Vitro Model of Physiological and Pathological Blood Flow with Application to Investigations of Vascular Cell Remodeling

Vascular disease is a common cause of death within the United States. Herein, we present a method to examine the contribution of flow dynamics towards vascular disease pathologies. Unhealthy arteries often present with wall stiffening, scarring, or partial stenosis which may all affect fluid flow rates, and the magnitude of pulsatile flow, or pulsatility index. Replication of various flow conditions is the result of tuning a flow pressure damping chamber downstream of a blood pump. Introduction of air within a closed flow system allows for a compressible medium to absorb pulsatile pressure from the pump, and therefore vary the pulsatility index. The method described herein is simply reproduced, with highly controllable input, and easily measurable results. Some limitations are recreation of the complex physiological pulse waveform, which is only approximated by the system. Endothelial cells, smooth muscle cells, and fibroblasts are affected by the blood flow through the artery. The dynamic component of blood flow is determined by the cardiac output and arterial wall compliance. Vascular cell mechano-transduction of flow dynamics may trigger cytokine release and cross-talk between cell types within the artery. Co-culture of vascular cells is a more accurate picture reflecting cell-cell interaction on the blood vessel wall and vascular response to mechanical signaling. Contribution of flow dynamics, including the cell response to the dynamic and mean (or steady) components of flow, is therefore an important metric in determining disease pathology and treatment efficacy. Through introducing an in vitro co-culture model and pressure damping downstream of blood pump which produces simulated cardiac output, various arterial disease pathologies may be investigated.     

Variability of Microcirculation Detected by Blood Pulsation Imaging

The non-invasive assessment of blood flow is invaluable for the diagnostic and monitoring treatment of numerous vascular and neurological diseases. We developed a non-invasive and non-contact method of blood pulsation imaging capable of visualizing and monitoring of the two-dimensional distribution of two key parameters of peripheral blood flow: the blood pulsation amplitude and blood pulsation phase. The method is based on the photoplethysmographic imaging in the reflection mode. In contrast with previous imaging systems we use new algorithm for data processing which allows two dimensional mapping of blood pulsations in large object''s areas after every cardiac cycle. In our study we carried out the occlusion test of the arm and found (i) the extensive variability of 2D-distribution of blood pulsation amplitude from one cardiac cycle to another, and (ii) existence of the adjacent spots to which the blood is asynchronously supplied. These observations show that the method can be used for studying of the multicomponent regulation of peripheral blood circulation. The proposed technique is technologically simple and cost-effective, which makes it applicable for monitoring the peripheral microcirculation in clinical settings for example, in diagnostics or testing the efficiency of new medicines.     

Contrast enhancement of speckle patterns from blood in synchrotron X-ray imaging

Hemodynamics has been a very important factor in understanding and diagnosing various vascular diseases. Recently, the X-ray particle image velocimetry (X-ray PIV) method using speckle patterns of blood has been introduced as a new quantitative visualization method for blood flows without any seeding tracer or contrast agents. In this study, the peculiar optical characteristics of blood on the synchrotron X-ray imaging method, which were not presented in previous studies, were investigated in depth and systematically. The experimental conditions required for X-ray PIV application were found to be the distance between the sample and the scintillator (～40 cm), the thickening of the blood sample (>0.3 mm), and hematocrit (20.0–80.0%). In addition, we verified that the X-ray PIV method is reliable as an advanced flow velocimetry by comparing the flow rate evaluated from the X-ray PIV result and the input flow rate supplied from a syringe pump with an error of less than 1%. Through this study, based on the understanding of contrast enhancement mechanisms of speckle patterns from blood, we could establish a trustworthy flow visualization method that can be used effectively in hemodynamic studies.     

A method for analyzing non-Newtonian blood viscosity data in low shear rates

Viscosity measurements were made using a coaxial rotating cylinder viscometer for blood having various volume fractions of red cells. A method is described for analyzing non-Newtonian blood viscosity in low shear rates by taking account of an increase or a decrease in size of red cell aggregates induced by shear. Our results are compared with empirical formulae presented by Scott Blair, Weaver-Evans-Walder and Thurston.     

The inverse Womersley problem for pulsatile flow in straight rigid tubes

In this study a numerical solution for the problem of pulsating flow in rigid tubes is described. The method applies to the case of known flow rate waveform, as opposed to Womersley solution where the pressure gradient was the known quantity. The solution provides the pressure gradient and wall shear stress waveforms as well as the instantaneous velocity profiles. Results show that the method can be used to study the blood flow characteristics in large arteries.