Thirty minutes of handgrip exercise enhances brachial artery flow-mediated dilation in response to two different shear stress profiles.

The walls of blood vessels are lined with a single-cell layer of endothelial cells. As blood flows through the arteries, a frictional force known as shear stress is sensed by mechanosensitive structures on the endothelium. Short and long term changes in shear stress can have a significant influence on the regulation of endothelial function. Acutely, shear stress triggers a pathway that culminates in the release of vasodilatory molecules from the endothelium and subsequent vasodilation of the artery. This endothelial response is known as flow mediated dilation (FMD). FMD is used as an index of endothelial function and is commonly assessed using reactive hyperemia (RH)-FMD, a method which elicits a large, short lived increase in shear stress following the release of a brief (5 min) forearm occlusion. A recent study found that a short term exposure (30 min) to a sustained elevation in shear stress potentiates subsequent RH-FMD. FMD can also result from a more prolonged, sustained increase in shear stress elicited by handgrip exercise (HGEX-FMD). There is evidence to suggest that interventions and conditions impact FMD resulting from sustained and transient shear stress stimuli differently, indicating that HGEX-FMD and RH-FMD provide different information about endothelial function. It is unknown whether HGEX-FMD is improved by short term exposure to shear stress. Understanding how exercise induced FMD is regulated is important because it contributes to blood flow responses during exercise. The study purpose was therefore to assess the impact of a handgrip exercise (intervention) induced sustained elevation in shear stress on subsequent brachial artery (BA) HGEX-FMD. Twenty healthy male participants (22±3yrs) preformed a 30-minute HGEX intervention on two experimental days. BA-FMD was assessed using either an RH or HGEX shear stress stimulus at 3 time points: pre-intervention, 10 min post and 60 min post. FMD and shear stress magnitude were determined via ultrasound. Both HGEX and RH-FMD increased significantly from pre-intervention to 10 min-post (p<0.01). These findings indicate that FMD stimulated by exercise induced increases in shear stress is potentiated by short term shear stress exposure. These findings advance our understanding regarding the regulation of endothelial function by shear stress.

Shear, shear stress and shearing

New simpler formulae are derived for the shear of a pair of material elements within the context of infinitesimal strain and finite strain. Also, new formulae are derived for shear stress based on the (symmetric) Cauchy stress and for the rate of shear of a pair of material elements within the rate of strain theory. These formulae are exploited to obtain results and to derive new simpler proofs of familiar classical results. In particular, a very simple short derivation is presented of the classical result of Coulomb and Hopkins on the maximum orthogonal shear stress. © 1992.

In situ characterization of sediment oxygen demand across the spectrum of non-resuspension to full resuspension with varying bed shear stress

Sediment oxygen demand (SOD) can be a significant oxygen sink in various types of water bodies, particularly slow-moving waters with substantial organic sediment accumulation. In most settings where SOD is a concern, the prevailing hydraulic conditions are such that the impact of sediment resuspension on SOD is not considered. However, in the case of Bubbly Creek in Chicago, Illinois, the prevailing slack water conditions are interrupted by infrequent intervals of very high flow rates associated with pumped combined sewer overflow (CSO) during intense hydrologic events. These events can cause resuspension of the highly organic, nutrient-rich bottom sediments, resulting in precipitous drawdown of dissolved oxygen (DO) in the water column. While many past studies have addressed the dependence of SOD on near-bed velocity and bed shear stress prior to the point of sediment resuspension, there has been limited research that has attempted to characterize the complex and dynamic phenomenon of resuspended-sediment oxygen demand. To address this issue, a new in situ experimental apparatus referred to as the U of I Hydrodynamic SOD Sampler was designed to achieve a broad range of velocities and associated bed shear stresses. This allowed SOD to be analyzed across the spectrum of no sediment resuspension associated with low velocity/ bed shear stress through full sediment resuspension associated with high velocity / bed shear stress. The current study split SOD into two separate components: (1) SODNR is the sediment oxygen demand associated with non-resuspension conditions and is a surface sink calculated using traditional methods to yield a value with units (g/m2/day); and (2) SODR is the oxygen demand associated with resuspension conditions, which is a volumetric sink most accurately characterized using non-traditional methods and units that reflect suspension in the water column (mg/L/day). In the case of resuspension, the suspended sediment concentration was analyzed as a function of bed shear stress, and a formulation was developed to characterize SODR as a function of suspended sediment concentration in a form similar to first-order biochemical oxygen demand (BOD) kinetics with Monod DO term. The results obtained are intended to be implemented into a numerical model containing hydrodynamic, sediment transport, and water quality components to yield oxygen demand varying in both space and time for specific flow events. Such implementation will allow evaluation of proposed Bubbly Creek water quality improvement alternatives which take into account the impact of SOD under various flow conditions. Although the findings were based on experiments specific to the conditions in Bubbly Creek, the techniques and formulations developed in this study should be applicable to similar sites.

The Effect of Architecture and Shear Stress on Endothelialization of 3D Printed Vascular Networks

Despite significant progress in the field of tissue engineering within the last decade, a number of unsolved problems still remain. One of the most relevant issues is the lack of proper vascularization that limits the size of engineered tissues to smaller than clinically relevant dimensions. In particular, the growth of engineered tissue in vitro within bioreactors is plagued with this challenge. Specifically, the tubular perfusion system bioreactor has been used for large scale bone constructs; however these engineered constructs lack inherent vasculature and quickly develop a hypoxic core, where no nutrient exchange can occur, thus leading to cell death. Through the use of 3D printed vascular templates in conjunction with a tubular perfusion system bioreactor, we attempt to create an endothelial cell monolayer on 3D scaffolds that could potentially serve as the foundation of inherent vasculature within these engineered bone grafts.

Molecular Origins of Higher Harmonics in Large-Amplitude Oscillatory Shear Flow: Shear Stress Response

Recent work has focused on deepening our understanding of the molecular origins of the higher harmonics that arise in the shear stress response of polymeric liquids in large-amplitude oscillatory shear flow. For instance, these higher harmonics have been explained by just considering the orientation distribution of rigid dumbbells suspended in a Newtonian solvent. These dumbbells, when in dilute suspension, form the simplest relevant molecular model of polymer viscoelasticity, and this model specifically neglects interactions between the polymer molecules [R.B. Bird et al., J Chem Phys, 140, 074904 (2014)]. In this paper, we explore these interactions by examining the Curtiss-Bird model, a kinetic molecular theory designed specifically to account for the restricted motions that arise when polymer chains are concentrated, thus interacting and specifically, entangled. We begin our comparison using a heretofore ignored explicit analytical solution [Fan and Bird, JNNFM, 15, 341 (1984)]. For concentrated systems, the chain motion transverse to the chain axis is more restricted than along the axis. This anisotropy is described by the link tension coefficient, ε, for which several special cases arise: ε = 0 corresponds to reptation, ε > 1/8 to rod-climbing, 1/2 ≥ ε ≥ 3/4 to reasonable predictions for shear-thinning in steady simple shear flow, and ε = 1 to the dilute solution without hydrodynamic interaction. In this paper, we examine the shapes of the shear stress versus shear rate loops for the special cases ε = (0,1/8, 3/8,1) , and we compare these with those of rigid dumbbell and reptation model predictions.

Specific shear-dependent viscoelastic third-grade fluid model

A specific modified constitutive equation for a third-grade fluid is proposed so that the model be suitable for applications where shear-thinning or shear-thickening may occur. For that, we use the Cosserat theory approach reducing the exact three-dimensional equations to a system depending only on time and on a single spatial variable. This one-dimensional system is obtained by integrating the linear momentum equation over the cross-section of the tube, taking a velocity field approximation provided by the Cosserat theory. From this reduced system, we obtain the unsteady equations for the wall shear stress and mean pressure gradient depending on the volume flow rate, Womersley number, viscoelastic coefficient and flow index over a finite section of the tube geometry with constant circular cross-section.

Drag reduction by polyethylene glycol in the tail arterial bed of normotensive and hypertensive rats

This study was designed to evaluate the effect of drag reducer polymers (DRP) on arteries from normotensive (Wistar) and spontaneously hypertensive rats (SHR). Polyethylene glycol (PEG 4000 at 5000 ppm) was perfused in the tail arterial bed with (E+) and without endothelium (E-) from male, adult Wistar (N = 14) and SHR (N = 13) animals under basal conditions (constant flow at 2.5 mL/min). In these preparations, flow-pressure curves (1.5 to 10 mL/min) were constructed before and 1 h after PEG 4000 perfusion. Afterwards, the tail arterial bed was fixed and the internal diameters of the arteries were then measured by microscopy and drag reduction was assessed based on the values of wall shear stress (WSS) by computational simulation. In Wistar and SHR groups, perfusion of PEG 4000 significantly reduced pulsatile pressure (Wistar/E+: 17.5 ± 2.8; SHR/E+: 16.3 ± 2.7%), WSS (Wistar/E+: 36; SHR/E+: 40%) and the flow-pressure response. The E- reduced the effects of PEG 4000 on arteries from both groups, suggesting that endothelial damage decreased the effect of PEG 4000 as a DRP. Moreover, the effects of PEG 4000 were more pronounced in the tail arterial bed from SHR compared to Wistar rats. In conclusion, these data demonstrated for the first time that PEG 4000 was more effective in reducing the pressure-flow response as well as WSS in the tail arterial bed of hypertensive than of normotensive rats and these effects were amplified by, but not dependent on, endothelial integrity. Thus, these results show an additional mechanism of action of this polymer besides its mechanical effect through the release and/or bioavailability of endothelial factors.

Flow Visualization in Arteriovenous Fistula and Aneurysm Using Computational Fluid Dynamics

The arteriovenous fistula (AVF) is characterized by enhanced blood flow and is the most widely used vascular access for chronic haemodialysis (Sivanesan et al., 1998). A large proportion of the AVF late failures are related to local haemodynamics (Sivanesan et al., 1999a). As in AVF, blood flow dynamics plays an important role in growth, rupture, and surgical treatment of aneurysm. Several techniques have been used to study the flow patterns in simplified models of vascular anastomose and aneurysm. In the present investigation, Computational Fluid Dynamics (CFD) is used to analyze the flow patterns in AVF and aneurysm through the velocity waveform obtained from experimental surgeries in dogs (Galego et al., 2000), as well as intra-operative blood flow recordings of patients with radiocephalic AVF ( Sivanesan et al., 1999b) and physiological pulses (Aires, 1991), respectively. The flow patterns in AVF for dog and patient surgeries data are qualitatively similar. Perturbation, recirculation and separation zones appeared during cardiac cycle, and these were intensified in the diastole phase for the AVF and aneurysm models. The values of wall shear stress presented in this investigation of AVF and aneurysm models oscillated in the range that can both cause damage to endothelial cells and develop atherosclerosis.

Computational fluid dynamics investigation of a centrifugal blood pump

In the development of a ventricular assist device, computational fluid dynamics (CFD) analysis is an efficient tool to obtain the best design before making the final prototype. In this study, different designs of a centrifugal blood pump were developed to investigate flow characteristics and performance. This study assumed the blood flow as being an incompressible homogeneous Newtonian fluid. A constant velocity was applied at the inlet; no slip boundary conditions were applied at device wall; and pressure boundary conditions were applied at the outlet. The CFD code used in this work was based on the finite volume method. In the future, the results of CFD analysis can be compared with flow visualization and hemolysis tests.

Turbulent blood flow plays an essential localizing role in the development of atherosclerotic lesions in experimentally induced hypercholesterolaemia in rats

Taking into account that atherosclerosis is a focal disease and high levels of plasma cholesterol are closely correlated with its pathogenesis, it is a challenge to explain how equal concentrations of cholesterol bathing the endothelium can produce local, rather than global, effects on arteries. The focal distribution of atherosclerotic lesions has been considered to be dependent, at least in part, on hydrodynamic factors. The present study was carried out to further test the hypothesis that these forces are an important localizing factor in rats feeding a hypercholesterolaemic diet and submitted to infra-diaphragmatic aortic constriction. These animals develop a normotensive prestenotic region with laminar blood flow that serves as control for a normotensive poststenotic region with turbulent blood flow. Our findings clearly demonstrated that the combination of turbulent blood flow and low wall shear stress (WSS) in the presence of hypercholesterolaemia and oxidative stress creates conditions to the formation of focally distributed incipient atherosclerotic lesions observed in the poststenotic segment. In contrast, only diffuse fatty streaks could be observed in the normotensive prestenotic segment with laminar blood flow and normal WSS in the presence of hypercholesterolaemia and oxidative stress. Although haemodynamic forces are not by themselves responsible for the pathogenesis of atherosclerosis, they prime the local vascular wall in which the lesion develop. Further studies are required to establish how haemodynamic forces are detected and transduced into chemical signalling by the cells of the artery wall and then converted into pathophysiologically relevant phenotypic changes.

Numerical investigation of turbulent flow around a rotating stepped cylinder for corrosion study

Direct numerical simulation has been carried out for turbulent flow set up by a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction. This flow geometry creates a qualitatively similar flow pattern as observed near, a sudden, pipe expansion or a plane backward-facing step, characterized by flow separation and reattachment. A region of intense turbulence intensity and high wall-shear-stress fluctuations is formed in,the recirculating I region downstream of the step, where high mass-transfer capacity was also experimentally observed. Since, corrosion is frequently mass-transfer., controlled, our findings, put forward this apparatus as useful tool for future corrosion research.

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

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.