The influence of geometric factors on the wall shear stress distribution in realistic human coronary arteries

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para a obtenção do grau de Mestre em Engenharia Biomédica. A presente dissertação foi desenvolvida no Erasmus Medical Center em Roterdão, Holanda

Static pressure and wall shear stress distributions in air flow in a seven wire-wrapped rod bundle

An experimental investigation is performed in a turbulent flow in a seven wire-wrapped rod bundle, mounted in an open air facility. Static pressure distributions are measured on central and peripheral rods. By using a Preston tube, the wall shear stress profiles are experimentally obtained along the perimeter of the rods. The geometric parameters of the test section are P/D=1.20 and H/D=15. The measuring section is located at L/D=40 from the air inlet. It is observed that the dimensionless static pressure and wall shear stress profiles are nearly independent of the Reynolds number and strongly dependent of the wire-spacer position, with abrupt variations of the parameters in the neighborhood of the wires.

Does the degree of coarctation of the aorta influence wall shear stress focal heterogeneity?

The development of atherosclerosis in the aorta is associated with low and oscillatory wall shear stress for normal patients. Moreover, localized differences in wall shear stress heterogeneity have been correlated with the presence of complex plaques in the descending aorta. While it is known that coarctation of the aorta can influence indices of wall shear stress, it is unclear how the degree of narrowing influences resulting patterns. We hypothesized that the degree of coarctation would have a strong influence on focal heterogeneity of wall shear stress. To test this hypothesis, we modeled the fluid dynamics in a patient-specific aorta with varied degrees of coarctation. We first validated a massively parallel computational model against experimental results for the patient geometry and then evaluated local shear stress patterns for a range of degrees of coarctation. Wall shear stress patterns at two cross sectional slices prone to develop atherosclerotic plaques were evaluated. Levels at different focal regions were compared to the conventional measure of average circumferential shear stress to enable localized quantification of coarctation-induced shear stress alteration. We find that the coarctation degree causes highly heterogeneous changes in wall shear stress.

Relation between plaque type, plaque thickness, blood shear stress, and plaque stress in coronary arteries assessed by X-ray Angiography and Intravascular Ultrasound

Purpose: Atheromatic plaque progression is affected, among others phenomena, by biomechanical, biochemical, and physiological factors. In this paper, the authors introduce a novel framework able to provide both morphological (vessel radius, plaque thickness, and type) and biomechanical (wall shear stress and Von Mises stress) indices of coronary arteries. Methods: First, the approach reconstructs the three-dimensional morphology of the vessel from intravascular ultrasound(IVUS) and Angiographic sequences, requiring minimal user interaction. Then, a computational pipeline allows to automatically assess fluid-dynamic and mechanical indices. Ten coronary arteries are analyzed illustrating the capabilities of the tool and confirming previous technical and clinical observations. Results: The relations between the arterial indices obtained by IVUS measurement and simulations have been quantitatively analyzed along the whole surface of the artery, extending the analysis of the coronary arteries shown in previous state of the art studies. Additionally, for the first time in the literature, the framework allows the computation of the membrane stresses using a simplified mechanical model of the arterial wall. Conclusions: Circumferentially (within a given frame), statistical analysis shows an inverse relation between the wall shear stress and the plaque thickness. At the global level (comparing a frame within the entire vessel), it is observed that heavy plaque accumulations are in general calcified and are located in the areas of the vessel having high wall shear stress. Finally, in their experiments the inverse proportionality between fluid and structural stresses is observed.

Relation between plaque type, plaque thickness, blood shear stress, and plaque stress in coronary arteries assessed by X-ray Angiography and Intravascular Ultrasound

Purpose: Atheromatic plaque progression is affected, among others phenomena, by biomechanical, biochemical, and physiological factors. In this paper, the authors introduce a novel framework able to provide both morphological (vessel radius, plaque thickness, and type) and biomechanical (wall shear stress and Von Mises stress) indices of coronary arteries. Methods: First, the approach reconstructs the three-dimensional morphology of the vessel from intravascular ultrasound(IVUS) and Angiographic sequences, requiring minimal user interaction. Then, a computational pipeline allows to automatically assess fluid-dynamic and mechanical indices. Ten coronary arteries are analyzed illustrating the capabilities of the tool and confirming previous technical and clinical observations. Results: The relations between the arterial indices obtained by IVUS measurement and simulations have been quantitatively analyzed along the whole surface of the artery, extending the analysis of the coronary arteries shown in previous state of the art studies. Additionally, for the first time in the literature, the framework allows the computation of the membrane stresses using a simplified mechanical model of the arterial wall. Conclusions: Circumferentially (within a given frame), statistical analysis shows an inverse relation between the wall shear stress and the plaque thickness. At the global level (comparing a frame within the entire vessel), it is observed that heavy plaque accumulations are in general calcified and are located in the areas of the vessel having high wall shear stress. Finally, in their experiments the inverse proportionality between fluid and structural stresses is observed.

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.

Can eccentric arterial plaques alone cause flow stagnation points and favour thrombus incorporation?

We have used an experimental model of aorta stenosis, with a Plexiglas plug, simulating a stable atheromatous plaque that promotes local turbulence and thrombosis. With animal survival of more than 24 h, we followed the partial fibrinolysis of the thrombus as well as its posterior organization and incorporation to the arterial wall as a neointima for up to 30 days. The mushroom plug form permitted the development of recirculation and stasis areas around it, favouring this evolution. Despite noted limitations, this study demonstrates that thrombus incorporation can contribute to plaque extension, as it can promote recirculation and stasis areas.

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.

The pulsed wire anemometer: review and further developments

The purpose of the present paper is to review work that has been done on the pulsed wire anemometer technique and also suggest further developments that could be made in its range of application. The aper discusses the three types of probes that have been used in pulsed wire anemometry: the crossed wire velocity probe, the parallel wire wall shear stress probe and the parallel wire velocity probe. The work shows that the crossed wire and the parallel wire techniques can be used to make velocity, turbulence and wall shear stress measurements in highly turbulent flows without any upper restriction on turbulence level. Comments are also made on the potential of a parallel wire probe for use in highly turbulent flows that would enable higher order velocity cross-product terms to be measured.