Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Proposal of a test specimen to evaluate the shear strength of vertical interfaces of running bond masonry walls

There is no normalized test to assess the shear strength of vertical interfaces of interconnected masonry walls. The approach used to evaluate this strength is normally indirect and often unreliable. The aim of this study is to propose a new test specimen to eliminate this deficiency. The main features of the proposed specimen are failure caused by shear stress on the vertical interface and a small number of units (blocks). The paper presents a numerical analysis based on the finite element method, with the purpose of showing the theoretical performance of the designed specimen, in terms of its geometry, boundary conditions, and loading scheme, and describes an experimental program using the specimen built with full- and third-scale clay blocks. The main conclusions are that the proposed specimen is easy to build and is appropriate to evaluate the sheaf strength of vertical interfaces of masonry walls.

Bending fatigue of stainless steel shear pins belonging to a hydroelectric plant

Coaracy Nunes was the first hydroelectric power plant in the Amazon region, being located in Araguari River, Amapa State, Brazil. The plant operates since 1976, presenting now a nominal capacity of 78 MW. The shear pins, which are installed in the turbine hydraulic arms to control the wicket gate and regulate the water flow into the turbine blades, suffered several breakdowns since 2004. These shear pins are made of an ASTM 410 stainless steel and were designed to break by a shear overload of 120 kN. Fractographic investigation of the pins, however, revealed two types of fracture topographies: a region of stable crack propagation area, with non-pronounced striation and secondary cracks; and a region of unstable propagation, featuring elongated dimples. These results indicated that the stable crack propagation occurred by fatigue (bidirectional bending), which was nucleated at machining marks under high nominal load. Finite element analysis was carried out using two loading conditions (pure shear and a combination of shear and bending) and the results indicated that the presence of a bending stress strongly increased the stress concentration factor (85% rise in the shear stress and 130% rise in the Von Mises stress). Misalignment during shear pins assembly associated with vibration might have promoted the premature failure of the shear by bending fatigue. (C) 2008 Elsevier Ltd. All rights reserved.

A new approach to heart valve tissue engineering: mimicking the heart ventricle with a ventricular assist device in a novel bioreactor

The `biomimetic` approach to tissue engineering usually involves the use of a bioreactor mimicking physiological parameters whilst supplying nutrients to the developing tissue. Here we present a new heart valve bioreactor, having as its centrepiece a ventricular assist device (VAD), which exposes the cell-scaffold constructs to a wider array of mechanical forces. The pump of the VAD has two chambers: a blood and a pneumatic chamber, separated by an elastic membrane. Pulsatile air-pressure is generated by a piston-type actuator and delivered to the pneumatic chamber, ejecting the fluid in the blood chamber. Subsequently, applied vacuum to the pneumatic chamber causes the blood chamber to fill. A mechanical heart valve was placed in the VAD`s inflow position. The tissue engineered (TE) valve was placed in the outflow position. The VAD was coupled in series with a Windkessel compliance chamber, variable throttle and reservoir, connected by silicone tubings. The reservoir sat on an elevated platform, allowing adjustment of ventricular preload between 0 and 11 mmHg. To allow for sterile gaseous exchange between the circuit interior and exterior, a 0.2 mu m filter was placed at the reservoir. Pressure and flow were registered downstream of the TE valve. The circuit was filled with culture medium and fitted in a standard 5% CO(2) incubator set at 37 degrees C. Pressure and flow waveforms were similar to those obtained under physiological conditions for the pulmonary circulation. The `cardiomimetic` approach presented here represents a new perspective to conventional biomimetic approaches in TE, with potential advantages. Copyright (C) 2010 John Wiley & Sons, Ltd.

EGG ALBUMIN AND GUAR GUM INFLUENCE ON FOAM THIXOTROPY

Food foams may be defined as products containing a gaseous phase stabilized in a matrix containing water and proteins. The rheology of a foam can be altered by how the air bubbles inside are organized, its size and shape, allowing a product that was initially fluid to be molded, like marshmallow. Rheology makes possible to understand how these properties are affected using dynamic oscillatory measurements to evaluate the behavior through G` and G `, and rotational to characterize the product through stress, and viscosity curves under controlled shear stress. These two methodologies combined are useful to analyze thixotropy. Nine formulations containing albumin and gelatin, and five others containing guar gum were evaluated by oscillatory rheology focusing thixotropy determination. Replacing gelatin by guar gum elasticity and thixotropy were increased, improving film stability around the air bubbles. It is possible to say that gelatin can be replaced by guar gum at concentrations of 0.4%, creating a product with a better ability of structure recovery. PRACTICAL APPLICATION This study indicated an important role of thixotropic analysis choice to understand the behavior of structure of products like foam. It was tested with three alternatives to get information about thixotropy of the products. It is very important for those who work with rheological analysis and products formulation.

Rheological behaviour of the mushy zone and its effect on the formation of casting defects during solidification

A new conceptual framework has been developed which explains the formation of shear-related casting defects such as porosity, segregation and tears. The theory relates defect formation to the mechanical behaviour of the partially solidified microstructure when shear stresses are developed during the filling of a casting and by the subsequent feeding processes during solidification. Two transition points, the dendrite coherency point and the maximum packing solid fraction, divide the mushy zone into three regions of different mechanical and feeding behaviours. The response of the mush to shear is related to the presence of these zones during solidification of a casting. The resulting defects are rationalized by considering the governing local shear stress and shear rate, local strength and time available for fluid flow. The design of the casting, the casting process used and the alloy composition all influence the relative importance of shearing on defect formation. (C) 1998 Acta Metallurgica Inc. Published by Elsevier Science Ltd. All rights reserved.

Shear deformation at 29% solid during solidification of magnesium alloy AZ91 and aluminium alloy A356

Partially solid commercial Al-Si and Mg-Al alloys have been deformed in shear during solidification using vane rheometry. The dendritic mush was deformed for a short period at 29% solid and allowed to cool naturally after deformation. Both alloys exhibited yield point behaviour and deformation was highly localised at the surface of maximum shear stress. The short period of deformation was found to have a distinct impact on the as-cast microstructure leading to fragmented dendrites in the deformation region of both alloys. In the case of the Mg-Al alloy, a concentrated region of interdendritic porosity was also observed in the deformation region. Concentrated porosity was not observed in the Al-Si alloy. (c) 2005 Elsevier B.V. All rights reserved.

Living bacteria rheology: Population growth, aggregation patterns, and collective behaviour under diferente shear flows

The activity of growing living bacteria was investigated using real-time and in situ rheology-in stationary and oscillatory shear. Two different strains of the human pathogen Staphylococcus aureus-strain COL and its isogenic cell wall autolysis mutant, RUSAL9-were considered in this work. For low bacteria density, strain COL forms small clusters, while the mutant, presenting deficient cell separation, forms irregular larger aggregates. In the early stages of growth, when subjected to a stationary shear, the viscosity of the cultures of both strains increases with the population of cells. As the bacteria reach the exponential phase of growth, the viscosity of the cultures of the two strains follows different and rich behaviors, with no counterpart in the optical density or in the population's colony-forming units measurements. While the viscosity of strain COL culture keeps increasing during the exponential phase and returns close to its initial value for the late phase of growth, where the population stabilizes, the viscosity of the mutant strain culture decreases steeply, still in the exponential phase, remains constant for some time, and increases again, reaching a constant plateau at a maximum value for the late phase of growth. These complex viscoelastic behaviors, which were observed to be shear-stress-dependent, are a consequence of two coupled effects: the cell density continuous increase and its changing interacting properties. The viscous and elastic moduli of strain COL culture, obtained with oscillatory shear, exhibit power-law behaviors whose exponents are dependent on the bacteria growth stage. The viscous and elastic moduli of the mutant culture have complex behaviors, emerging from the different relaxation times that are associated with the large molecules of the medium and the self-organized structures of bacteria. Nevertheless, these behaviors reflect the bacteria growth stage.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Velocity-jump instabilities in Hele-Shaw flow of associating polymer solutions

We study fracturelike flow instabilities that arise when water is injected into a Hele-Shaw cell filled with aqueous solutions of associating polymers. We explore various polymer architectures, molecular weights, and solution concentrations. Simultaneous measurements of the finger tip velocity and of the pressure at the injection point allow us to describe the dynamics of the finger in terms of the finger mobility, which relates the velocity to the pressure gradient. The flow discontinuities, characterized by jumps in the finger tip velocity, which are observed in experiments with some of the polymer solutions, can be modeled by using a nonmonotonic dependence between a characteristic shear stress and the shear rate at the tip of the finger. A simple model, which is based on a viscosity function containing both a Newtonian and a non-Newtonian component, and which predicts nonmonotonic regions when the non-Newtonian component of the viscosity dominates, is shown to agree with the experimental data.

The Effects of Headcut and Knickpoint Propagation on Bridges in Iowa, TR-541, 2008

The overarching goal of the proposed research was to provide a predictive tool for knickpoint propagation within the HCA (Hungry Canyon Alliance) territory. Knickpoints threaten the stability of bridge structures in Western Iowa. The study involved detailed field investigations over two years in order to monitor the upstream migration of a knickpoint on Mud Creek in Mills County, IA and identify the key mechanisms triggering knickpoint propagation. A state-of-the-art laser level system mounted on a movable truss provided continuous measurements of the knickpoint front for different flow conditions. A pressure transducer found in proximity of the truss provided simultaneous measurements of the flow depth. The laser and pressure transducer measurements led to the identification of the conditions at which the knickpoint migration commences. It was suggested that negative pressures developed by the reverse roller flow near the toe of the knickpoint face triggered undercutting of the knickpoint at this location. The pressure differential between the negative pressure and the atmospheric pressure also draws the impinging jet closer to the knickpoint face producing scour. In addition, the pressure differential may induce suction of sediment from the face. Other contributing factors include slump failure, seepage effects, and local fluvial erosion due to the exerted fluid shear. The prevailing flow conditions and soil information along with the channel cross-sectional geometry and gradient were used as inputs to a transcritical, one dimensional, hydraulic/geomorphic numerical model, which was used to map the flow characteristics and shear stress conditions near the knickpoint. Such detailed flow calculations do not exist in the published literature. The coupling of field and modeling work resulted in the development of a blueprint methodology, which can be adopted in different parts of the country for evaluating knickpoint evolution. This information will assist local government agencies in better understanding the principal factors that cause knickpoint propagation and help estimate the needed response time to control the propagation of a knickpoint after one has been identified.

Analysis of flow dynamics in right ventricular outflow tract.

BACKGROUND: The mechanism behind early graft failure after right ventricular outflow tract (RVOT) reconstruction is not fully understood. Our aim was to establish a three-dimensional computational fluid dynamics (CFD) model of RVOT to investigate the hemodynamic conditions that may trigger the development of intimal hyperplasia and arteriosclerosis. METHODS: Pressure, flow, and diameter at the RVOT, pulmonary artery (PA), bifurcation of the PA, and left and right PAs were measured in 10 normal pigs with a mean weight of 24.8 ± 0.78 kg. Data obtained from the experimental scenario were used for CFD simulation of pressure, flow, and shear stress profile from the RVOT to the left and right PAs. RESULTS: Using experimental data, a CFD model was obtained for 2.0 and 2.5-L/min pulsatile inflow profiles. In both velocity profiles, time and space averaged in the low-shear stress profile range from 0-6.0 Pa at the pulmonary trunk, its bifurcation, and at the openings of both PAs. These low-shear stress areas were accompanied to high-pressure regions 14.0-20.0 mm Hg (1866.2-2666 Pa). Flow analysis revealed a turbulent flow at the PA bifurcation and ostia of both PAs. CONCLUSIONS: Identified local low-shear stress, high pressure, and turbulent flow correspond to a well-defined trigger pattern for the development of intimal hyperplasia and arteriosclerosis. As such, this real-time three-dimensional CFD model may in the future serve as a tool for the planning of RVOT reconstruction, its analysis, and prediction of outcome.

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

Computational modeling of stented coronary arteries

The application of computational fluid dynamics (CFD) and finite element analysis (FEA) has been growing rapidly in the various fields of science and technology. One of the areas of interest is in biomedical engineering. The altered hemodynamics inside the blood vessels plays a key role in the development of the arterial disease called atherosclerosis, which is the major cause of human death worldwide. Atherosclerosis is often treated with the stenting procedure to restore the normal blood flow. A stent is a tubular, flexible structure, usually made of metals, which is driven and expanded in the blocked arteries. Despite the success rate of the stenting procedure, it is often associated with the restenosis (re-narrowing of the artery) process. The presence of non-biological device in the artery causes inflammation or re-growth of atherosclerotic lesions in the treated vessels. Several factors including the design of stents, type of stent expansion, expansion pressure, morphology and composition of vessel wall influence the restenosis process. Therefore, the role of computational studies is crucial in the investigation and optimisation of the factors that influence post-stenting complications. This thesis focuses on the stent-vessel wall interactions followed by the blood flow in the post-stenting stage of stenosed human coronary artery. Hemodynamic and mechanical stresses were analysed in three separate stent-plaque-artery models. Plaque was modeled as a multi-layer (fibrous cap (FC), necrotic core (NC), and fibrosis (F)) and the arterial wall as a single layer domain. CFD/FEA simulations were performed using commercial software packages in several models mimicking the various stages and morphologies of atherosclerosis. The tissue prolapse (TP) of stented vessel wall, the distribution of von Mises stress (VMS) inside various layers of vessel wall, and the wall shear stress (WSS) along the luminal surface of the deformed vessel wall were measured and evaluated. The results revealed the role of the stenosis size, thickness of each layer of atherosclerotic wall, thickness of stent strut, pressure applied for stenosis expansion, and the flow condition in the distribution of stresses. The thicknesses of FC, and NC and the total thickness of plaque are critical in controlling the stresses inside the tissue. A small change in morphology of artery wall can significantly affect the distribution of stresses. In particular, FC is the most sensitive layer to TP and stresses, which could determine plaque’s vulnerability to rupture. The WSS is highly influenced by the deflection of artery, which in turn is dependent on the structural composition of arterial wall layers. Together with the stenosis size, their roles could play a decisive role in controlling the low values of WSS (<0.5 Pa) prone to restenosis. Moreover, the time dependent flow altered the percentage of luminal area with WSS values less than 0.5 Pa at different time instants. The non- Newtonian viscosity model of the blood properties significantly affects the prediction of WSS magnitude. The outcomes of this investigation will help to better understand the roles of the individual layers of atherosclerotic vessels and their risk to provoke restenosis at the post-stenting stage. As a consequence, the implementation of such an approach to assess the post-stented stresses will assist the engineers and clinicians in optimizing the stenting techniques to minimize the occurrence of restenosis.

Understanding the mechanisms of lung mechanical stress

Physical forces affect both the function and phenotype of cells in the lung. Bronchial, alveolar, and other parenchymal cells, as well as fibroblasts and macrophages, are normally subjected to a variety of passive and active mechanical forces associated with lung inflation and vascular perfusion as a result of the dynamic nature of lung function. These forces include changes in stress (force per unit area) or strain (any forced change in length in relation to the initial length) and shear stress (the stress component parallel to a given surface). The responses of cells to mechanical forces are the result of the cell's ability to sense and transduce these stimuli into intracellular signaling pathways able to communicate the information to its interior. This review will focus on the modulation of intracellular pathways by lung mechanical forces and the intercellular signaling. A better understanding of the mechanisms by which lung cells transduce physical forces into biochemical and biological signals is of key importance for identifying targets for the treatment and prevention of physical force-related disorders.