Numerical analysis of hemodynamics in spastic middle cerebral arteries

Cerebral vasospasm (CVS) is the most common serious complication of subarachnoid hemorrhage. Among the many factors that are associated with the pathogenesis of CVS, hemodynamics plays an important role in the initiation and development of CVS. Numerical simulation was carried out to obtain the flow patterns and wall shear stress (WSS) distribution in spastic middle cerebral arteries. The blood was assumed to be incompressible, laminar, homogenous, Newtonian, and steady. Our simulations reveal that flow velocity and WSS level increase at the stenosis segment of the spastic vessels, but further downstream of stenosis, the WSS significantly decreases along the inner wall, and flow circulation and stagnation are observed. The hydrodynamic resistance increases with the increase of vessel spasm. Moreover, the change of flow field and hydrodynamic forces are not linearly proportional to the spasm level, and the rapid change of hemodynamic parameters is observed as the spasm is more than 50%. Accordingly, in the view of hemodynamic physiology, vessels with less than 30% stenosis are capable of self-restoration towards normal conditions. However, vessels with more than 50% stenosis may eventually lose their capacity to adapt to differing physiologic conditions due to the extreme non-physilogic hemodynamic environment, and the immediate expansion of the vessel lumen might be needed to minimize serious and non-reversible effects.     

The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

Numerical identification of the rupture locations in patient-specific abdominal aortic aneurysmsusing hemodynamic parameters

The rupture of an abdominal aortic aneurysm (AAA) is generally an unexpected event. Up to now, there is no agreement on an accurate criteria to predict the rupture risk of AAAs. This paper aims to numerically investigate the hemodynamics of three ruptured and one non-ruptured patient-specific AAA models to correlate local hemodynamic parameters with the rupture sites, and for the first time, this study introduced helicity as a potential index for the rupture potential of AAAs.3D reconstructions from CT scans were done. The simulation revealed that all the rupture sites were in regions of stagnation with near zero wall shear stress (WSS) but large WSS gradient (WSSG), which may explain the observation by the former researchers that the rupture site in the ruptured AAA has the lowest recorded wall thickness compared to other non-ruptured regions. Moreover, all the ruptures occurred at regions of zero helicity which represents a purely axial or circumferential flow. In addition, this study revealed that the double low region for the non-ruptured AAA was present with a thick layer of plaques, it suggests that the AAA rupture and the formation of atherosclerotic plaques may share a lot common physiological features. However, the fact that there are no plaques present in the walls of three RAAAs also indicates that AAA is not always a result of atherosclerosis. The current computational study may complement the maximum diameter, peak wall stress and other clinically relevant factors in AAA ruptures to identify the rupture sites of AAAs.     

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.     

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.     

The effects of stenosis severity on the hemodynamic parameters-assessment of the correlation between stress phase angle and wall shear stress

To study the effects of increase in the degree of stenosis severity and subsequent complexity of hemodynamic patterns on hemodynamic parameters, experimental investigations and numerical simulations were performed. The correlations between the large negative Stress Phase Angle (SPA), the low mean Wall Shear Stress (WSS) and high Oscillatory Shear Index (OSI) were investigated at the distal shoulder and post-stenotic regions as the outcomes of elevated stenosis severity. Models included non-Newtonian fluid flow in stenotic arteries with 30-80% symmetrical stenoses. To study the interactions between pulsatile WSS and pulsatile wall circumferential stress (WCS) acting on endothelial cells, SPA as the phase difference between WSS and WCS waves was used. Moreover, the distribution of SPA on the lumen axis was compared to the distributions of the mean WSS and OSI that have been regarded until now as the determinants of atherosclerosis-prone regions. Results indicate that an increase in stenosis severity, not only affects the mean WSS, mean WCS and pulse amplitudes, but also influences the phase difference between them. The SPA is large negative on the distal shoulder and post-stenotic areas where atherosclerotic plaque develops. The increasing stenosis severity and the subsequent increasing complexity of hemodynamic patterns affect the correlation between any of the low mean WSS and high OSI with large negative SPA, such that it not only leads to create and develop some regions where the correlation between any of the low mean WSS and high OSI with large negative SPA is well but also leads to create and develop other regions where such correlations fail.     

Does Lower Limb Exercise Worsen Renal Artery Hemodynamics in Patients with Abdominal Aortic Aneurysm?

Renal artery stenosis (RAS) and renal complications emerge in some patients after endovascular aneurysm repair (EVAR) to treat abdominal aorta aneurysm (AAA). The mechanisms for the causes of these problems are not clear. We hypothesized that for EVAR patients, lower limb exercise could negatively influence the physiology of the renal artery and the renal function, by decreasing the blood flow velocity and changing the hemodynamics in the renal arteries. To evaluate this hypothesis, pre- and post-operative models of the abdominal aorta were reconstructed based on CT images. The hemodynamic environment was numerically simulated under rest and lower limb exercise conditions. The results revealed that in the renal arteries, lower limb exercise decreased the wall shear stress (WSS), increased the oscillatory shear index (OSI) and increased the relative residence time (RRT). EVAR further enhanced these effects. Because these parameters are related to artery stenosis and atherosclerosis, this preliminary study concluded that lower limb exercise may increase the potential risk of inducing renal artery stenosis and renal complications for AAA patients. This finding could help elucidate the mechanism of renal artery stenosis and renal complications after EVAR and warn us to reconsider the management and nursing care of AAA patients.     

Patient-specific analysis of displacement forces acting on fenestrated stent grafts for endovascular aneurysm repair

Treatment options for abdominal aortic aneurysm (AAA) include highly invasive open surgical repair or minimally invasive endovascular aneurysm repair (EVAR). Despite being minimally invasive, some patients are not suitable for EVAR due to hostile AAA morphology. Fenestrated-EVAR (F-EVAR) was introduced to address these limitations of standard EVAR, where AAA is treated using a Fenestrated Stent Graft (FSG). In order to assess durability of F-EVAR, displacement forces acting on FSGs were analysed in this study, based on patient-specific geometries reconstructed from computed tomography (CT) scans. The magnitude and direction of the resultant displacement forces acting on the FSG were numerically computed using computational fluid dynamics (CFD) with a rigid wall assumption. Although displacement force arises from blood pressure and friction due to blood flow, numerical simulations elucidated that net blood pressure is the dominant contributor to the overall displacement force; as a result, time dependence of the resultant displacement force followed pressure waveform very closely. The magnitude of peak displacement force varied from 1.9 N to 14.3 N with a median of 7.0 N. A strong positive correlation was found between inlet cross-sectional area (CSA), anterior/posterior (A/P) angle and the peak displacement force i.e. as inlet CSA or A/P angle increases, the magnitude of resultant displacement increases. This study manifests that while loads exerted by the pulsatile flow dictates the cyclic variation of the displacement force, its magnitude depends not only on blood pressure but also the FSG morphology, with the latter determining the direction of the displacement force.     

Hemodynamic investigation of a patient-specific abdominal aortic aneurysm with iliac artery tortuosity

Abstract

This paper describes a systematic investigation on the hemodynamic environment in a patient-specific AAA with tortuous common iliac artery(CIA) and external iliac artery (EIA). 3D reconstructions from CT scans and subsequent computational simulation are carried out. It is found out that the Newtonian and non-Newtonian models have very similar flow field and WSS distribution. More importantly, it is revealed that the torturous CIA maintained its helical flow. It is concluded that the assumption of Newtonian blood is adequate in capturing the intra-aneurysmal hemodynamics. Moreover, it is speculated that the physiological spiral flow protects the twisted CIA from the thrombosis formation.     

Impedance-based outflow boundary conditions for human carotid haemodynamics

In this study, we develop structured tree outflow boundary conditions for modelling the human carotid haemodynamics. The model geometry was reconstructed through computerised tomography scan. Unsteady-state computational fluid dynamic analyses were performed under different conditions using a commercial software package ADINA R&D, Inc., (Watertown, MA, USA) in order to assess the impact of the boundary conditions on the flow variables. In particular, the results showed that the peripheral vessels massively impact the pressure while the flow is relatively unaffected. As an example of application of these outflow conditions, an unsteady fluid-structure interaction (FSI) simulation was carried out and the dependence of the wall shear stress (WSS) on the arterial wall compliance in the carotid bifurcation was studied. In particular, a comparison between FSI and rigid-wall models was conducted. Results showed that the WSS distributions were substantially affected by the diameter variation of the arterial wall. In particular, even similar WSS distributions were found for both cases, and differences in the computed WSS values were also found.     

Optimization of cardiovascular stent design using computational fluid dynamics

Coronary stent design affects the spatial distribution of wall shear stress (WSS), which can influence the progression of endothelialization, neointimal hyperplasia, and restenosis. Previous computational fluid dynamics (CFD) studies have only examined a small number of possible geometries to identify stent designs that reduce alterations in near-wall hemodynamics. Based on a previously described framework for optimizing cardiovascular geometries, we developed a methodology that couples CFD and three-dimensional shape-optimization for use in stent design. The optimization procedure was fully-automated, such that solid model construction, anisotropic mesh generation, CFD simulation, and WSS quantification did not require user intervention. We applied the method to determine the optimal number of circumferentially repeating stent cells (N(C)) for slotted-tube stents with various diameters and intrastrut areas. Optimal stent designs were defined as those minimizing the area of low intrastrut time-averaged WSS. Interestingly, we determined that the optimal value of N(C) was dependent on the intrastrut angle with respect to the primary flow direction. Further investigation indicated that stent designs with an intrastrut angle of approximately 40 deg minimized the area of low time-averaged WSS regardless of vessel size or intrastrut area. Future application of this optimization method to commercially available stent designs may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes.     

Characterization of the hemodynamic wall shear stresses in human umbilical vessels from normal and intrauterine growth restricted pregnancies

Significant reductions in blood flow and umbilical diameters were reported in pregnancies affected by intrauterine growth restriction (IUGR) from placental insufficiency. However, it is not known if IUGR umbilical blood vessels experience different hemodynamic wall shear stresses (WSS) compared to normal umbilical vessels. As WSS is known to influence vasoactivity and vascular growth and remodeling, which can regulate flow rates, it is important to study this parameter. In this study, we aim to characterize umbilical vascular WSS environment in normal and IUGR pregnancies, and evaluate correlation between WSS and vascular diameter, and gestational age. Twenty-two normal and 21 IUGR pregnancies were assessed via ultrasound between the 27th and 39th gestational week. IUGR was defined as estimated fetal weight and/or abdominal circumference below the 10th centile, with no improvement during the remainder of the pregnancy. Vascular diameter was determined by 3D ultrasound scans and image segmentation. Umbilical artery (UA) WSS was computed via computational flow simulations, while umbilical vein (UV) WSS was computed via the Poiseuille equation. Univariate multiple regression analysis was used to test for the differences between normal and IUGR cohort. UV volumetric flow rate, UA and UV diameters were significantly lower in IUGR fetuses, but flow velocities and WSS trends in UA and UV were very similar between normal and IUGR groups. In both groups, UV WSS showed a significant negative correlation with diameter, but UA WSS had no correlation with diameter, suggesting a constancy of WSS environment and the existence of WSS homeostasis in UA, but not in UV. Despite having reduced flow rate and vascular sizes, IUGR UAs had hemodynamic mechanical stress environments and trends that were similar to those in normal pregnancies. This suggested that endothelial dysfunction or abnormal mechanosensing was unlikely to be the cause of small vessels in IUGR umbilical cords.     

Flow-induced wall shear stress in abdominal aortic aneurysms: Part I--steady flow hemodynamics

Numerical predictions of blood flow patterns and hemodynamic stresses in Abdominal Aortic Aneurysms (AAAs) are performed in a two-aneurysm, axisymmetric, rigid wall model using the spectral element method. Homogeneous, Newtonian blood flow is simulated under steady conditions for the range of Reynolds numbers 10 < or =Re < or =2265. Flow hemodynamics are quantified by calculating the distributions of wall pressure (p(w)), wall shear stress (tau(w)), Wall Shear Stress Gradient (WSSG). A correlation between maximum values of hemodynamic stresses and Reynolds number is established, and the spatial distribution of WSSG is considered as a hemodynamic force that may cause damage to the arterial wall at an intermediate stage of AAA growth. The temporal distribution of hemodynamic stresses in pulsatile flow and their physical implications in AAA rupture are discussed in Part II of this paper.     

Can temporal fluctuation in spatial wall shear stress gradient initiate a cerebral aneurysm? A proposed novel hemodynamic index,the gradient oscillatory number (GON)

We propose a new hemodynamic index for the initiation of a cerebral aneurysm, defined by the temporal fluctuations of tension/compression forces acting on endothelial cells. We employed a patient-specific geometry of a human internal carotid artery (ICA) with an aneurysm, and reconstructed the geometry of the ICA before aneurysm formation by artificially removing the aneurysm. We calculated the proposed hemodynamic index and five other hemodynamic indices (wall shear stress (WSS) at peak systole, time-averaged WSS, time-averaged spatial WSS gradient, oscillatory shear index (OSI), and potential aneurysm formation indicator (AFI)) for the geometry before aneurysm formation using a computational fluid dynamics technique. By comparing the distribution of each index at the location of aneurysm formation, we discussed the validity of each. The results showed that only the proposed hemodynamic index had a significant correlation with the location of aneurysm formation. Our findings suggest that the proposed index may be useful as a hemodynamic index for the initiation of cerebral aneurysms.     

Simulations of morphotype-dependent hemodynamics in non-dilated bicuspid aortic valve aortas

Bicuspid aortic valves (BAVs) generate flow abnormalities that may promote aortopathy. While positive helix fraction (PHF) index, flow angle (θ), flow displacement (d) and wall shear stress (WSS) exhibit abnormalities in dilated BAV aortas, it is unclear whether those anomalies stem from the abnormal valve anatomy or the dilated aorta. Therefore, the objective of this study was to quantify the early impact of different BAV morphotypes on aorta hemodynamics prior to dilation. Fluid-structure interaction models were designed to quantify standard peak-systolic flow metrics and temporal WSS characteristics in a realistic non-dilated aorta connected to functional tricuspid aortic valve (TAV) and type-I BAVs. While BAVs generated increased helicity (PHF>0.68) in the middle ascending aorta (AA), larger systolic flow skewness (θ>11.2°) and displacement (d>6.8 mm) relative to the TAV (PHF=0.51; θ<5.5°; d<3.3 mm), no distinct pattern was observed between morphotypes. In contrast, WSS magnitude and directionality abnormalities were BAV morphotype- and site-dependent. Type-I BAVs subjected the AA convexity to peak-systolic WSS overloads (up to 1014% difference vs. TAV). While all BAVs increased WSS unidirectionality on the proximal AA relative to the TAV, the most significant abnormality was achieved by the BAV with left-right-coronary cusp fusion on the wall convexity (up to 0.26 decrease in oscillatory shear index vs. TAV). The results indicate the existence of strong hemodynamic abnormalities in non-dilated type-I BAV AAs, their colocalization with sites vulnerable to dilation and the superior specificity of WSS metrics over global hemodynamic metrics to the valve anatomy.     

Developed pulsatile flow in a deployed coronary stent

The patho-physiologic process of restenosis and tissue growth may not be completely eliminated and is the primary concern of clinicians performing angioplasty and stent implantation procedures. Recent evidence suggests that the restenosis process is influenced by several factors: (1) geometry and size of vessel; (2) stent design; and (3) it's location that alter hemodynamic parameters, including local wall shear stress (WSS) distributions. The present three-dimensional (3D) analysis of pulsatile flow in a deployed coronary stent: (1) shows complex 3D variation of hemodynamic parameters; and (2) quantifies the changes in local WSS distributions for developed flow and compares with recently published WSS data for developing flow. Higher order of magnitude of WSS of 290 dyn/cm(2) is observed on the surface of cross-link intersections at the entrance of the stent for developed flow, which is about half of that for developing flow. Low WSS of 0.8 dyn/cm(2) and negative WSS of -8 dyn/cm(2) are seen at the immediate upstream and downstream regions of strut intersections. Persistent recirculation is observed at the downstream region of each strut cross-link and the regions of low and negative WSS may lead to patho-physiologic conditions near the stented region. The key finding of this study is that the location of stent in the coronary artery determines the developing or developed nature of the flow, which in turn, results in varied level of WSS.     

Comparison of the hemodynamics in 6mm and 4-7 mm hemodialysis grafts by means of CFD

The aim of our study is to investigate with computational fluid dynamics (CFD) whether different arterial anastomotic geometries result in a different hemodynamics at the arterial (AA) and venous anastomosis (VA) of hemodialysis vascular access grafts. We have studied a 6mm graft (CD) and a 4-7 mm graft (TG). A validated three-dimensional CFD model is developed to simulate flow in the two graft types. Only the arterial anastomosis (AA) geometry differs. The boundary conditions applied are a periodic velocity signal at the arterial inlet and a periodic pressure wave at the venous outlet. Flow rate is set to 1,000 ml/min. The time dependent Navier-Stokes equations are solved. Wall shear stress (WSS), wall shear stress gradient (WSSG) and pressure gradient (PG) are calculated. Anastomotic flow is asymmetric although the anastomosis geometry is symmetric. The hemodynamic parameters, WSS, WSSG and PG, values at the suture line of the arterial anastomosis of the TG are at least twice as much as in the CD. Comparing the parameters at the two AA indicate that little flow rate increase introduces the risk of hemolysis in the TG whereas the CD is completely free of hemolysis. The hemodynamic parameter values at the venous anastomosis of the CD are 24 till 35% higher compared to the values of the TG. WSS values (> 3 Pa) in the VA are in the critical range for stenosis development in both graft geometries. The zones where the parameters reach extreme values correspond to the locations where intimal hyperplasia formation is reported in literature. In all anastomoses, the hemodynamic parameter levels are in the range where leucocytes and platelets get activated. Our simulations confirm clinical results where TG did not show a better outcome when compared to the CD.     

Curvature effect on hemodynamic conditions at the inner bend of the carotid siphon and its relation to aneurysm formation

Although high-impact hemodynamic forces are thought to lead to cerebral aneurysmal change, little is known about the aneurysm formation on the inner aspect of vascular bends such as the internal carotid artery (ICA) siphon where wall shear stress (WSS) is expected to be low. This study evaluates the effect of vessel curvature and hemodynamics on aneurysm formation along the inner carotid siphon. Catheter 3D-rotational angiographic volumes of 35 ICA (10 aneurysms, 25 controls) were evaluated in 3D for radius of curvature and peak curvature of the siphon bend, followed by univariate statistical analysis. Computational fluid dynamic (CFD) simulations were performed on patient-derived models after aneurysm removal and on synthetic variants of increasing curvature. Peak focal siphon curvature was significantly higher in aneurysm bearing ICAs (0.36±0.045 vs. 0.30±0.048 mm⁻¹, p=0.003), with no difference in global radius of curvature (p=0.36). In CFD simulations, increasing parametric curvature tightness (from 5 to 3 mm radius) resulted in dramatic increase of WSS and WSS gradient magnitude (WSSG) on the inner wall of the bend. In patient-derived data, the location of aneurysms coincided with regions of low WSS (<4 Pa) flanked by high WSS and WSSG peaks. WSS peaks correlated with the aneurysm neck. In contrast, control siphon bends displayed low, almost constant, WSS and WSSG profiles with little spatial variation. High bend curvature induces dynamically fluctuating high proximal WSS and WSSG followed by regions of flow stasis and recirculation, leading to local conditions known to induce destructive vessel wall remodeling and aneurysmal initiation.     

Effects of vessel compliance on flow pattern in porcine epicardial right coronary arterial tree

The compliance of the vessel wall affects hemodynamic parameters which may alter the permeability of the vessel wall. Based on experimental measurements, the present study established a finite element (FE) model in the proximal elastic vessel segments of epicardial right coronary arterial (RCA) tree obtained from computed tomography. The motion of elastic vessel wall was measured by an impedance catheter and the inlet boundary condition was measured by an ultrasound flow probe. The Galerkin FE method was used to solve the Navier–Stokes and Continuity equations, where the convective term in the Navier–Stokes equation was changed in the arbitrary Lagrangian–Eulerian (ALE) framework to incorporate the motion due to vessel compliance. Various hemodynamic parameters (e.g., wall shear stress—WSS, WSS spatial gradient—WSSG, oscillatory shear index—OSI) were analyzed in the model. The motion due to vessel compliance affects the time-averaged WSSG more strongly than WSS at bifurcations. The decrease of WSSG at flow divider in elastic bifurcations, as compared to rigid bifurcations, implies that the vessel compliance decreases the permeability of vessel wall and may be atheroprotective. The model can be used to predict coronary flow pattern in subject-specific anatomy as determined by noninvasive imaging.