Characterization of the in vivo wall shear stress environment of human fetus umbilical arteries and veins

The endothelial cells of the umbilical vessels are frequently used in mechanobiology experiments. They are known to respond to wall shear stress (WSS) of blood flow, which influences vascular growth and remodeling. The in vivo environment of umbilical vascular WSS, however, is not well characterized. In this study, we performed detailed characterization of the umbilical vascular WSS environments using clinical ultrasound scans combined with computational simulations. Doppler ultrasound scans of 28 normal human fetuses from 32nd to 33rd gestational weeks were investigated. Vascular cross-sectional areas were quantified through 3D reconstruction of the vascular geometry from 3D B-mode ultrasound images, and flow velocities were quantified through pulse wave Doppler. WSS in umbilical vein was computed with Poiseuille’s equation, whereas WSS in umbilical artery was obtained via computational fluid dynamics simulations of the helical arterial geometry. Results showed that blood flow velocity for umbilical artery and vein did not correlate with vascular sizes, suggesting that velocity had a very weak trend with or remained constant over vascular sizes. Average WSS for umbilical arteries and vein was 2.81 and 0.52 Pa, respectively. Umbilical vein WSS showed a significant negative correlation with the vessel diameter, but umbilical artery did not show any correlation. We hypothesize that this may be due to differential regulation of vascular sizes based on WSS sensing. Due to the helical geometry of umbilical arteries, bending of the umbilical cord did not significantly alter the vascular resistance or WSS, unlike that in the umbilical veins. We hypothesize that the helical shape of umbilical arteries may be an adaptation feature to render a higher constancy of WSS and flow in the arteries despite umbilical cord bending.     

Wall shear stress and near-wall flows in the stenosed femoral artery

Stenotic artery hemodynamics are often characertised by metrics including oscillatory shear index (OSI) and residence time (RT). This analysis was conducted to clarify the link between the near-wall flow behaviour and these resultant flow metrics. A computational simulation was conducted of a stenosed femoral artery, with an idealised representative geometry and a physiologically realistic inlet profile. The overall flow behaviour was characterised through consideration of the axial flow, which was non-dimensionalised against mean flow velocity. The OSI and RT metrics, which are a useful indicator of likely atherosclerotic sites, were explained through a discussion of the WSS values at different time points, the velocity behaviour and velocity profiles, with a particular focus on the near-wall behaviour which influences wall shear stress and the transient evolution of the wall shear stress. While, the stenosis throat experiences high values of wall shear stress, the smooth flow through this contracted region results in low variation in wall shear stress vectors and limited opportunity for any particle stasis. However, regions were noted distal and proximal (though to a lesser extent), where the change in recirculation zones over the cycle created highly elevated regions of both OSI and RT.     

A study on the compliance of a right coronary artery and its impact on wall shear stress

A computational model incorporating physiological motion and uniform transient wall deformation of a branchless right coronary artery (RCA) was developed to assess the influence of artery compliance on wall shear stress (WSS). Arterial geometry and deformation were derived from modern medical imaging techniques, whereas the blood flow was solved numerically employing a moving-grid approach using a well-validated in-house finite element code. The simulation results indicate that artery compliance affects the WSS in the RCA heterogeneously, with the distal region mostly experiencing these effects. Under physiological inflow conditions, coronary compliance contributed to phase changes in the WSS time history, without affecting the temporal gradient of the local WSS nor the bounds of the WSS magnitude. Compliance does not cause considerable changes to the topology of WSS vector patterns nor to the localization of WSS minima along the RCA. We conclude that compliance is not an important factor affecting local hemodynamics in the proximal region of the RCA while the influence of compliance in the distal region needs to be evaluated in conjunction with the outflow to the myocardium through the major branches of the RCA.     

Requirements for mesh resolution in 3D computational hemodynamics

Computational techniques are widely used for studying large artery hemodynamics. Current trends favor analyzing flow in more anatomically realistic arteries. A significant obstacle to such analyses is generation of computational meshes that accurately resolve both the complex geometry and the physiologically relevant flow features. Here we examine, for a single arterial geometry, how velocity and wall shear stress patterns depend on mesh characteristics. A well-validated Navier-Stokes solver was used to simulate flow in an anatomically realistic human right coronary artery (RCA) using unstructured high-order tetrahedral finite element meshes. Velocities, wall shear stresses (WSS), and wall shear stress gradients were computed on a conventional "high-resolution" mesh series (60,000 to 160,000 velocity nodes) generated with a commercial meshing package. Similar calculations were then performed in a series of meshes generated through an adaptive mesh refinement (AMR) methodology. Mesh-independent velocity fields were not very difficult to obtain for both the conventional and adaptive mesh series. However, wall shear stress fields, and, in particular, wall shear stress gradient fields, were much more difficult to accurately resolve. The conventional (nonadaptive) mesh series did not show a consistent trend towards mesh-independence of WSS results. For the adaptive series, it required approximately 190,000 velocity nodes to reach an r.m.s. error in normalized WSS of less than 10 percent. Achieving mesh-independence in computed WSS fields requires a surprisingly large number of nodes, and is best approached through a systematic solution-adaptive mesh refinement technique. Calculations of WSS, and particularly WSS gradients, show appreciable errors even on meshes that appear to produce mesh-independent velocity fields.     

Numerical simulations of pulsatile flow in an end-to-side anastomosis model

A potential interaction between the local hemodynamics and the artery wall response has been suggested for vascular graft failure by intimal hyperplasia (IH). Among the various hemodynamic factors, wall shear stress has been implicated as the primary factor responsible for the development of IH. In order to explore the role of hemodynamics in the formation of IH in end-to-side anastomosis, computational fluid dynamics is employed. To validate the numerical simulations, comparisons with existing experimental data are performed for both steady and pulsatile flows. Generally, good agreement is observed with the velocity profiles whereas some discrepancies are found in wall shear stress (WSS) distributions. Using the same end-to-side anastomosis geometry, numerical simulations are extended using a femoral artery waveform to identify the possible role of unsteady hemodynamics. In the current simulations, Carreau-Yasuda model is used to account for the non-Newtonian nature of blood. Computations indicated a disturbed flow field at the artery-graft junction leading to locally elevated shear stresses on the vascular wall. Furthermore, the shear stress distribution followed the same behavior with oscillating magnitude over the entire flow cycle. Thus, distal IH observed in end-to-side artery-graft models may be caused by the fluctuations in WSS's along the wall.     

Stent design properties and deployment ratio influence indexes of wall shear stress: a three-dimensional computational fluid dynamics investigation within a normal artery.

Restenosis limits the effectiveness of stents, but the mechanisms responsible for this phenomenon remain incompletely described. Stent geometry and expansion during deployment produce alterations in vascular anatomy that may adversely affect wall shear stress (WSS) and correlate with neointimal hyperplasia. These considerations have been neglected in previous computational fluid dynamics models of stent hemodynamics. Thus we tested the hypothesis that deployment diameter and stent strut properties (e.g., number, width, and thickness) influence indexes of WSS predicted with three-dimensional computational fluid dynamics. Simulations were based on canine coronary artery diameter measurements. Stent-to-artery ratios of 1.1 or 1.2:1 were modeled, and computational vessels containing four or eight struts of two widths (0.197 or 0.329 mm) and two thicknesses (0.096 or 0.056 mm) subjected to an inlet velocity of 0.105 m/s were examined. WSS and spatial WSS gradients were calculated and expressed as a percentage of the stent and vessel area. Reducing strut thickness caused regions subjected to low WSS (<5 dyn/cm(2)) to decrease by approximately 87%. Increasing the number of struts produced a 2.75-fold increase in exposure to low WSS. Reducing strut width also caused a modest increase in the area of the vessel experiencing low WSS. Use of a 1.2:1 deployment ratio increased exposure to low WSS by 12-fold compared with stents implanted in a 1.1:1 stent-to-vessel ratio. Thinner struts caused a modest reduction in the area of the vessel subjected to elevated WSS gradients, but values were similar for the other simulations. The results suggest that stent designs that reduce strut number and thickness are less likely to subject the vessel to distributions of WSS associated with neointimal hyperplasia.     

Alterations in wall shear stress predict sites of neointimal hyperplasia after stent implantation in rabbit iliac arteries

Restenosis resulting from neointimal hyperplasia (NH) limits the effectiveness of intravascular stents. Rates of restenosis vary with stent geometry, but whether stents affect spatial and temporal distributions of wall shear stress (WSS) in vivo is unknown. We tested the hypothesis that alterations in spatial WSS after stent implantation predict sites of NH in rabbit iliac arteries. Antegrade iliac artery stent implantation was performed under angiography, and blood flow was measured before casting 14 or 21 days after implantation. Iliac artery blood flow domains were obtained from three-dimensional microfocal X-ray computed tomography imaging and reconstruction of the arterial casts. Indexes of WSS were determined using three-dimensional computational fluid dynamics. Vascular histology was unchanged proximal and distal to the stent. Time-dependent NH was localized within the stented region and was greatest in regions exposed to low WSS and acute elevations in spatial WSS gradients. The lowest values of WSS spatially localized to the stented area of a theoretical artery progressively increased after 14 and 21 days as NH occurred within these regions. This NH abolished spatial disparity in distributions of WSS. The results suggest that stents may introduce spatial alterations in WSS that modulate NH in vivo.     

Wall shear stresses remain elevated in mature arteriovenous fistulas: a case study

Maintaining vascular access (VA) patency continues to be the greatest challenge for dialysis patients. VA dysfunction, primarily due to venous neointimal hyperplasia development and stenotic lesion formation, is mainly attributed to complex hemodynamics within the arteriovenous fistula (AVF). The effect of VA creation and the subsequent geometrical remodeling on the hemodynamics and shear forces within a mature patient-specific AVF is investigated. A 3D reconstructed geometry of a healthy vein and a fully mature patient-specific AVF was developed from a series of 2D magnetic resonance image scans. A previously validated thresholding technique for region segmentation and lumen cross section contour creation was conducted in MIMICS 10.01, allowing for the creation of a 3D reconstructed geometry. The healthy vein and AVF computational models were built, subdivided, and meshed in GAMBIT 2.3. The computational fluid dynamic (CFD) code FLUENT 6.3.2 (Fluent Inc., Lebanon, NH) was employed as the finite volume solver to determine the hemodynamics and shear forces within the healthy vein and patient-specific AVF. Geometrical alterations were evaluated and a CFD analysis was conducted. Substantial geometrical remodeling was observed, following VA creation with an increase in cross-sectional area, out of plane curvature (maximum angle of curvature in AVF=30?deg), and angle of blood flow entry. The mean flow velocity entering the vein of the AVF is dramatically increased. These factors result in complex three-dimensional hemodynamics within VA junction (VAJ) and efferent vein of the AVF. Complex flow patterns were observed and the maximum and mean wall shear stress (WSS) magnitudes are significantly elevated. Flow reversal was found within the VAJ and efferent vein. Extensive geometrical remodeling during AVF maturation does not restore physiological hemodynamics to the VAJ and venous conduit of the AVF, and high WSS and WSS gradients, and flow reversal persist. It is theorized that the vessel remodelling and the continued non-physiological hemodynamics within the AVF compound to result in stenotic lesion development.     

MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.     

Image segmentation methods for intracranial aneurysm haemodynamic research

Patient-specific haemodynamic technology is being increasingly utilised in clinical applications. Under normal circumstances, computational haemodynamic simulation is performed using geometric results obtained via medical image segmentation. However, even when employed upon the same set of medical imaging data, both the geometry and volume of intracranial aneurysm models are highly dependent upon varying insufficiently validated vascular segmentation methods. In this study, we compared three segmentation methods to segment the geometry of the aneurysm. These include: the Region Growing Threshold (RGT), Chan-Vese model (CV) and Threshold-Based Level Set (TLS). The results obtained were evaluated via measurement of arterial volume differences (VD), local geometric shapes, and haemodynamic simulation results. In total, 45 patient-specific aneurysm cases with three different anatomy locations were assessed in this study. From this, we discovered that the average VD of all three segmentation methods lay in the vicinity of 9.3% (SD=±4.6%). The computational haemodynamic simulation was performed via the use of the vessel geometries. Analyses produced an average of 23.2% (SD=±8.7%) difference in energy loss (EL) between the varying segmentation methods, with the difference in Wall Shear Stress (WSS) averaging 24.0% (SD=±8.5%) and 126.4% (SD=±124.4%) for the highest and lowest volumes of WSS respectively. The results of the lowest WSS, was seen to be significantly dependent upon the geometry of the aneurysm surface. It is therefore essential, in order to confirm the quality of segmentation processes in the application of patient-specific analyses of cerebrovascular haemodynamics – to validate these individual segmentation methods.     

The effect of proximal artery flow on the hemodynamics at the distal anastomosis of a vascular bypass graft: computational study

The formation of distal anastomotic intimal hyperplasia (IH), one common mode of bypass graft failure, has been shown to occur in the areas of disturbed flow particular to this site. The nature of theflow in the segment of artery proximal to the distal anastomosis varies from case to case depending on the clinical situation presented. A partial stenosis of a bypassed arterial segment may allow residual prograde flow through the proximal artery entering the distal anastomosis of the graft. A complete stenosis may allow for zero flow in the proximal artery segment or retrograde flow due to the presence of small collateral vessels upstream. Although a number of investigations on the hemodynamics at the distal anastomosis of an end-to-side bypass graft have been conducted, there has not been a uniform treatment of the proximal artery flow condition. As a result, direct comparison of results from study to study may not be appropriate. The purpose of this work was to perform a three-dimensional computational investigation to study the effect of the proximal artery flow condition (i.e., prograde, zero, and retrograde flow) on the hemodynamics at the distal end-to-side anastomosis. We used the finite volume method to solve the full Navier-Stokes equations for steady flow through an idealized geometry of the distal anastomosis. We calculated the flow field and local wall shear stress (WSS) and WSS gradient (WSSG) everywhere in the domain. We also calculated the severity parameter (SP), a quantification of hemodynamic variation, at the anastomosis. Our model showed a marked difference in both the magnitude and spatial distribution of WSS and WSSG. For example, the maximum WSS magnitude on the floor of the artery proximal to the anastomosis for the prograde and zero flow cases is 1.8 and 3.9 dynes/cm2, respectively, while it is increased to 10.3 dynes/cm2 in the retrograde flow case. Similarly, the maximum value of WSSG magnitude on thefloor of the artery proximal to the anastomosis for the prograde flow case is 4.9 dynes/cm3, while it is increased to 13.6 and 24.2 dynes/cm3, respectively, in the zero and retrograde flow cases. The value of SP is highest for the retrograde flow case (13.7 dynes/cm3) and 8.1 and 12.1 percent lower than this for the prograde (12.6 dynes/cm3) and zero (12.0 dynes/cm3) flow cases, respectively. Our model results suggest that the flow condition in the proximal artery is an important determinant of the hemodynamics at the distal anastomosis of end-to-side vascular bypass grafts. Because hemodynamic forces affect the response of vascular endothelial cells, the flow situation in the proximal artery may affect IH formation and, therefore, long-term graft patency. Since surgeons have some control over the flow condition in the proximal artery, results from this study could help determine which flow condition is clinically optimal.     

Low WSS Induces Intimal Thickening,while Large WSS Variation and Inflammation Induce Medial Thinning,in an Animal Model of Atherosclerosis

Objective

Atherosclerotic plaque development in the arterial wall is the result of complex interaction between the wall’s endothelial layer and blood hemodynamics. However, the interaction between hemodynamic parameters and inflammation in plaque evolution is not yet fully understood. The aim of the present study was to investigate the relation between wall shear stress (WSS) and vessel wall inflammation during atherosclerotic plaque development in a minipig model of carotid stenosis.

Methods

A surgical procedure was performed to create left common carotid artery stenosis by placement of a perivascular cuff in minipigs under atherogenic diet. Animals were followed up on 3T MRI, 1 week after surgery and 3, 6, and 8 months after initiation of the diet. Computational fluid dynamics simulation estimated WSS distribution for the first imaging point. Vascular geometries were co-registered for direct comparison of plaque development and features (Gadolinium- and USPIO-Contrast Enhanced MRI, for permeability and inflammation respectively) with the initial WSS. Histological analysis was performed and sections were matched to MR images, based on spatial landmarks.

Results

Vessel wall thickening, permeability and inflammation were observed distally from the stenosis. They were eccentric and facing regions of normal wall thickness. Histological analysis confirmed eccentric plaque formation with lipid infiltration, intimal thickening and medial degradation. High phagocytic activity in the stenosis region was co-localized with high WSS, corresponding to intense medial degradation observed on histology samples.

Conclusion

Lower WSS promotes atherosclerotic plaque development distal to an induced stenosis. Vascular and perivascular inflammation locations were predominant in the high WSS stenosis segment, where medial thinning was the major consequence.     

Allometric scaling of wall shear stress from mice to humans: quantification using cine phase-contrast MRI and computational fluid dynamics

Allometric scaling laws relate structure or function between species of vastly different sizes. They have rarely been derived for hemodynamic parameters known to affect the cardiovascular system, e.g., wall shear stress (WSS). This work describes noninvasive methods to quantify and determine a scaling law for WSS. Geometry and blood flow velocities in the infrarenal aorta of mice and rats under isoflurane anesthesia were quantified using two-dimensional magnetic resonance angiography and phase-contrast magnetic resonance imaging at 4.7 tesla. Three-dimensional models constructed from anatomic data were discretized and used for computational fluid dynamic simulations using phase-contrast velocity imaging data as inlet boundary conditions. WSS was calculated along the infrarenal aorta and compared between species to formulate an allometric equation for WSS. Mean WSS along the infrarenal aorta was significantly greater in mice and rats compared with humans (87.6, 70.5, and 4.8 dyn/cm(2), P < 0.01), and a scaling exponent of -0.38 (R(2) = 0.92) was determined. Manipulation of the murine genome has made small animal models standard surrogates for better understanding the healthy and diseased human cardiovascular system. It has therefore become increasingly important to understand how results scale from mouse to human. This noninvasive methodology provides the opportunity to serially quantify changes in WSS during disease progression and/or therapeutic intervention.     

Pulsatile flow in an end-to-side vascular graft model: comparison of computations with experimental data

Various hemodynamic factors have been implicated in vascular graft intimal hyperplasia, the major mechanism contributing to chronic failure of small-diameter grafts. However, a thorough knowledge of the graft flow field is needed in order to determine the role of hemodynamics and how these factors affect the underlying biological processes. Computational fluid dynamics offers much more versatility and resolution than in vitro or in vivo methods, yet computations must be validated by careful comparison with experimental data. Whereas numerous numerical and in vitro simulations of arterial geometries have been reported, direct point-by-point comparisons of the two techniques are rare in the literature. We have conducted finite element computational analyses for a model of an end-to-side vascular graft and compared the results with experimental data obtained using laser-Doppler velocimetry. Agreement for velocity profiles is found to be good, with some clear differences near the recirculation zones during the deceleration and reverse-flow segments of the flow waveform. Wall shear stresses are determined from velocity gradients, whether by computational or experimental methods, and hence the agreement for this quantity, while still good, is less consistent than for velocity itself from the wall shear stress numerical results, we computed four variables that have been cited in the development of intiimal hyperplasia-the time-averaged wall shear stress, an oscillating shear index, and spatial and temporal wall shear stress gradients in order to illustrate the versatility of numerical methods. We conclude that the computational approach is a valid alternative to the experimental approach for quantitative hemodynamic studies. Where differences in velocity were found by the two methods, it was generally attributed to the inability of the numerical method to model the fluid dynamics when flow conditions are destabilizing. Differences in wall shear, in the absence of destabilizing phenomena, were more likely to be caused by difficulties in calculating wall shear from relatively low resolution in vitro data.     

Towards the improved quantification of in vivo abnormal wall shear stresses in BAV-affected patients from 4D-flow imaging: Benchmarking and application to real data

Bicuspid aortic valve (BAV), i.e. the fusion of two aortic valve cusps, is the most frequent congenital cardiac malformation. Its progression is often characterized by accelerated leaflet calcification and aortic wall dilation. These processes are likely enhanced by altered biomechanical stimuli, including fluid-dynamic wall shear stresses (WSS) acting on both the aortic wall and the aortic valve. Several studies have proposed the exploitation of 4D-flow magnetic resonance imaging sequences to characterize abnormal in vivo WSS in BAV-affected patients, to support prognosis and timing of intervention. However, current methods fail to quantify WSS peak values.On this basis, we developed two new methods for the improved quantification of in vivo WSS acting on the aortic wall based on 4D-flow data.We tested both methods separately and in combination on synthetic datasets obtained by two computational fluid-dynamics (CFD) models of the aorta with healthy and bicuspid aortic valve. Tests highlighted the need for data spatial resolution at least comparable to current clinical guidelines, the low sensitivity of the methods to data noise, and their capability, when used jointly, to compute more realistic peak WSS values as compared to state-of-the-art methods.The integrated application of the two methods on the real 4D-flow data from a preliminary cohort of three healthy volunteers and three BAV-affected patients confirmed these indications. In particular, quantified WSS peak values were one order of magnitude higher than those reported in previous 4D-flow studies, and much closer to those computed by highly time- and space-resolved CFD simulations.     

Impact of hemodynamics on lumen boundary displacements in abdominal aortic aneurysms by means of dynamic computed tomography and computational fluid dynamics

The aim of the present work is to quantitatively assess the three-dimensional distributions of the displacements experienced during the cardiac cycle by the luminal boundary of abdominal aortic aneurysm (AAA) and to correlate them with the local bulk hemodynamics. Ten patients were acquired by means of time resolved computed tomography, and each patient-specific vascular morphology was reconstructed for all available time frames. The AAA lumen boundary motion was tracked, and the lumen boundary displacements (LBD) computed for each time frame. The intra-aneurysm hemodynamic quantities, specifically wall shear stress (WSS), were evaluated with computational fluid dynamics simulations. Co-localization of LBD and WSS distributions was evaluated by means of Pearson correlation coefficient. A clear anisotropic distribution of LBD was evidenced in both space and time; a combination of AAA lumen boundary inward- and outward-directed motions was assessed. A co-localization between largest outward LBD and high WSS was demonstrated supporting the hypothesis of a mechanistic relationship between anisotropic displacement and hemodynamic forces related to the impingement of the blood on the lumen boundary. The presence of anisotropic displacement of the AAA lumen boundary and their link to hemodynamic forces have been assessed, highlighting a new possible role for hemodynamics in the study of AAA progression.     

Effects of different stent designs on local hemodynamics in stented arteries

Following the deployment of a coronary stent and disruption of an atheromatous plaque, the deformation of the arterial wall and the presence of the stent struts create a new fluid dynamic field, which can cause an abnormal biological response. In this study 3D computational models were used to analyze the fluid dynamic disturbances induced by the placement of a stent inside a coronary artery. Stents models were first expanded against a simplified arterial plaque, with a solid mechanics analysis, and then subjected to a fluid flow simulation under pulsatile physiological conditions. Spatial and temporal distribution of arterial wall shear stress (WSS) was investigated after the expansion of stents of different designs and different strut thicknesses. Common oscillatory WSS behavior was detected in all stent models. Comparing stent and vessel wall surfaces, maximum WSS values (in the order of 1Pa) were located on the stent surface area. WSS spatial distribution on the vascular wall surface showed decreasing values from the center of the vessel wall portion delimited by the stent struts to the wall regions close to the struts. The hemodynamic effects induced by two different thickness values for the same stent design were investigated, too, and a reduced extension of low WSS region (<0.5Pa) was observed for the model with a thicker strut.     

Multiphase hemodynamic simulation of pulsatile flow in a coronary artery

A multiphase transient non-Newtonian three-dimensional (3-D) computational fluid dynamics (CFD) simulation has been performed for pulsatile hemodynamics in an idealized curved section of a human coronary artery. We present the first prediction, to the authors' knowledge, of particulate buildup on the inside curvature using the multiphase theory of dense suspension hemodynamics. In this study, the particulates are red blood cells (RBCs). The location of RBC buildup on the inside curvature correlates with lower wall shear stress (WSS) relative to the outside curvature. These predictions provide insight into how blood-borne particulates interact with artery walls and hence, have relevance for understanding atherogenesis since clinical observations show that atherosclerotic plaques generally form on the inside curvatures of arteries. The buildup of RBCs on the inside curvature is driven by the secondary flow and higher residence times. The higher viscosity in the central portion of the curved vessel tends to block their flow, causing them to migrate preferentially through the boundary layer. The reason for this is the nearly neutrally buoyant nature of the dense two-phase hemodynamic flow. The two-phase non-Newtonian viscosity model predicts greater shear thinning than the single-phase non-Newtonian model. Consequently, the secondary flow induced in the curvature is weaker. The waveforms for computed hemodynamic parameters, such as hematocrit, WSS, and viscosity, follow the prescribed inlet velocity waveforms. The lower oscillatory WSS produced on the inside curvature has implications for understanding thickening of the intimal layer.     

Circumferential vascular deformation after stent implantation alters wall shear stress evaluated with time-dependent 3D computational fluid dynamics models.

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that stent geometry may cause local alterations in wall shear stress (WSS) that have been associated with neointimal hyperplasia and subsequent restenosis. However, previous CFD studies have ignored histological evidence of vascular straightening between circumferential stent struts. We tested the hypothesis that consideration of stent-induced vascular deformation may more accurately predict alterations in indexes of WSS that may subsequently account for histological findings after stenting. We further tested the hypothesis that the severity of these alterations in WSS varies with the degree of vascular deformation after implantation. Steady-state and time-dependent simulations of three-dimensional CFD arteries based on canine coronary artery measurements of diameter and blood flow were conducted, and WSS and WSS gradients were calculated. Circumferential straightening introduced areas of high WSS between stent struts that were absent in stented vessels of circular cross section. The area of vessel exposed to low WSS was dependent on the degree of circumferential vascular deformation and axial location within the stent. Stents with four vs. eight struts increased the intrastrut area of low WSS in vessels, regardless of cross-sectional geometry. Elevated WSS gradients were also observed between struts in vessels with polygonal cross sections. The results obtained using three-dimensional CFD models suggest that changes in vascular geometry after stent implantation are important determinants of WSS distributions that may be associated with subsequent neointimal hyperplasia.