Haemodynamic assessment of human coronary arteries is affected by degree of freedom of artery movement

Abnormal haemodynamic parameters are associated with atheroma plaque progression and instability in coronary arteries. Flow recirculation, shear stress and pressure gradient are understood to be important pathogenic mediators in coronary disease. The effect of freedom of coronary artery movement on these parameters is still unknown. Fluid–structure interaction (FSI) simulations were carried out in 25 coronary artery models derived from authentic human coronaries in order to investigate the effect of degree of freedom of movement of the coronary arteries on flow recirculation, wall shear stress (WSS) and wall pressure gradient (WPG). Each FSI model had distinctive supports placed upon it. The quantitative and qualitative differences in flow recirculation, maximum wall shear stress (MWSS), areas of low wall shear stress (ALWSS) and maximum wall pressure gradient (MWPG) for each model were determined. The results showed that greater freedom of movement was associated with lower MWSS, smaller ALWSS, smaller flow recirculation zones and lower MWPG. With increasing percentage diameter stenosis (%DS), the effect of degree of freedom on flow recirculation and WSS diminished. Freedom of movement is an important variable to be considered for computational modelling of human coronary arteries, especially in the setting of mild to moderate stenosis.

Abbreviations: 3D: Three-dimensional; 3DR: Three-dimensional Reconstruction; 3D-QCA: Three-dimensional quantitative coronary angiography; ALWSS: Areas of low wall shear stress; CAD: Coronary artery disease; CFD: Computational fluid dynamics; %DS: Diameter stenosis percentage; EPCS: End point of counter-rotating streamlines; FSI: Fluid–structure interaction; IVUS: Intravascular ultrasound; LAD: Left anterior descending; MWSS: Maximum wall shear stress; SST: Shear stress transport; TAWSS: Time-averaged wall shear stress; WSS: wall shear stress; WPG: Wall pressure gradient; MWPG: Maximum wall pressure gradient; FFR: Fractional flow reserve; iFR: Instantaneous wave-free ratio     

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

Computational fluid dynamics study of common stent models inside idealised curved coronary arteries

The haemodynamic behaviour of blood inside a coronary artery after stenting is greatly affected by individual stent features as well as complex geometrical properties of the artery including tortuosity and curvature. Regions at higher risk of restenosis, as measured by low wall shear stress (WSS < 0.5 Pa), have not yet been studied in detail in curved stented arteries. In this study, three-dimensional computational modelling and computational fluid dynamics methodologies were used to analyse the haemodynamic characteristics in curved stented arteries using several common stent models. Results in this study showed that stent strut thickness was one major factor influencing the distribution of WSS in curved arteries. Regions of low WSS were found behind struts, particularly those oriented at a large angle relative to the streamwise flow direction. These findings were similar to those obtained in studies of straight arteries. An uneven distribution of WSS at the inner and outer bends of curved arteries was observed where the WSS was lower at the inner bend. In this study, it was also shown that stents with a helical configuration generated an extra swirling component of the flow based on the helical direction; however, this extra swirl in the flow field did not cause significant changes on the distribution of WSS under the current setup.     

Impact of plaques in the left coronary artery on wall shear stress and pressure gradient in coronary side branches

In this study, we investigate plaques located at the left coronary bifurcation. We focus on the effect that the resulting changes in wall shear stress (WSS) and wall pressure stress gradient (WPSG) have on atherosclerotic progress in coronary artery disease. Coronary plaques were simulated and placed at the left main stem and the left anterior descending to produce >50% narrowing of the coronary lumen. Computational fluid dynamics analysis was carried out, simulating realistic physiological conditions that show the in vivo cardiac haemodynamic. WSS and WPSG in the left coronary artery were calculated and compared in the left coronary models, with and without the presence of plaques during cardiac cycles. Our results showed that WSS decreased while WPSG was increased in coronary side branches due to the presence of plaques. There is a direct correlation between coronary plaques and subsequent WSS and WPSG variations based on the bifurcation plaques simulated in the realistic coronary models.     

Understanding the fluid mechanics behind transverse wall shear stress

The patchy distribution of atherosclerosis within arteries is widely attributed to local variation in haemodynamic wall shear stress (WSS). A recently-introduced metric, the transverse wall shear stress (transWSS), which is the average over the cardiac cycle of WSS components perpendicular to the temporal mean WSS vector, correlates particularly well with the pattern of lesions around aortic branch ostia. Here we use numerical methods to investigate the nature of the arterial flows captured by transWSS and the sensitivity of transWSS to inflow waveform and aortic geometry. TransWSS developed chiefly in the acceleration, peak systolic and deceleration phases of the cardiac cycle; the reverse flow phase was too short, and WSS in diastole was too low, for these periods to have a significant influence. Most of the spatial variation in transWSS arose from variation in the angle by which instantaneous WSS vectors deviated from the mean WSS vector rather than from variation in the magnitude of the vectors. The pattern of transWSS was insensitive to inflow waveform; only unphysiologically high Womersley numbers produced substantial changes. However, transWSS was sensitive to changes in geometry. The curvature of the arch and proximal descending aorta were responsible for the principal features, the non-planar nature of the aorta produced asymmetries in the location and position of streaks of high transWSS, and taper determined the persistence of the streaks down the aorta. These results reflect the importance of the fluctuating strength of Dean vortices in generating transWSS.     

Influence of overlapping pattern of multiple overlapping uncovered stents on the local mechanical environment: A patient-specific parameter study

BackgroundMultiple overlapping uncovered stents (MOUS) system has shown potentials in managing complex aortic aneurysms with side branches involvement. It promotes the development of thrombus by modulating local flow pattern that reduces the wall tension, while maintaining patency of side branches. However the modulation of local hemodynamic parameters depends on various factors that have not been assessed comprehensively.MethodsAneurysm 3D geometry was reconstructed based on CT images. One-way fluid-structure interaction analysis was performed to quantify structural stress concentration in the wall, and changes of blood velocity, wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT) and pressure in the sac due to the stent deployment.ResultsHigh structural stress concentration due to stent deployment was found in the landing zone and it increased linearly with the number of stents deployed. The wall tension in the sac was unaffected by the stent deployment. Stress within the wall was insensitive to the different overlapping pattern. After one stent was deployed, the mean flow velocity in the sac reduced by 36.4%. The deployment of the 2nd stent further reduced the mean sac velocity by 10%. WSS decreased while both OSI and RRT increased after stent deployment, however pressure in the sac remained nearly unchanged. Except for the cases with complete stents struts alignment, different overlapping pattern had little effect on flow parameters.ConclusionsMechanical parameters modulated by the MOUS are insensitive to different overlapping pattern suggesting that endovascular procedure can be performed with less attention to the overlapping pattern.     

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.     

Mis-sizing of stent promotes intimal hyperplasia: impact of endothelial shear and intramural stress

Stent can cause flow disturbances on the endothelium and compliance mismatch and increased stress on the vessel wall. These effects can cause low wall shear stress (WSS), high wall shear stress gradient (WSSG), oscillatory shear index (OSI), and circumferential wall stress (CWS), which may promote neointimal hyperplasia (IH). The hypothesis is that stent-induced abnormal fluid and solid mechanics contribute to IH. To vary the range of WSS, WSSG, OSI, and CWS, we intentionally mismatched the size of stents to that of the vessel lumen. Stents were implanted in coronary arteries of 10 swine. Intravascular ultrasound (IVUS) was used to size the coronary arteries and stents. After 4 wk of stent implantation, IVUS was performed again to determine the extent of IH. In conjunction, computational models of actual stents, the artery, and non-Newtonian blood were created in a computer simulation to yield the distribution of WSS, WSSG, OSI, and CWS in the stented vessel wall. An inverse relation (R(2) = 0.59, P < 0.005) between WSS and IH was found based on a linear regression analysis. Linear relations between WSSG, OSI, and IH were observed (R(2) = 0.48 and 0.50, respectively, P < 0.005). A linear relation (R(2) = 0.58, P < 0.005) between CWS and IH was also found. More statistically significant linear relations between the ratio of CWS to WSS (CWS/WSS), the products CWS × WSSG and CWS × OSI, and IH were observed (R(2) = 0.67, 0.54, and 0.56, respectively, P < 0.005), suggesting that both fluid and solid mechanics influence the extent of IH. Stents create endothelial flow disturbances and intramural wall stress concentrations, which correlate with the extent of IH formation, and these effects were exaggerated with mismatch of stent/vessel size. These findings reveal the importance of reliable vessel and stent sizing to improve the mechanics on the vessel wall and minimize IH.     

Non-Newtonian Blood Flow in Left Coronary Arteries with Varying Stenosis: A Comparative Study

This paper presents Computational fluid dynamic (CFD) analysis of blood flow in three different 3-D models of left coronary artery (LCA). A comparative study of flow parameters (pressure distribution, velocity distribution and wall shear stress) in each of the models is done for a non-Newtonian (Carreau) as well as the Newtonian nature of blood viscosity over a complete cardiac cycle. The difference between these two types of behavior of blood is studied for both transient and steady states of flow. Additionally, flow parameters are compared for steady and transient boundary conditions considering blood as non-Newtonian fluid. The study shows that the highest wall shear stress (WSS), velocity and pressure are found in artery having stenosis in all the three branches of LCA. The use of Newtonian blood model is a good approximation for steady as well as transient blood flow boundary conditions if shear rate is above 100 s-1. However, the assumption of steady blood flow results in underestimating the values of flow parameters such as wall shear stress, pressure and velocity.     

Comparison of physiological and post-endovascular aneurysm repair infrarenal blood flow

Endovascular aneurysm repair (EVAR) of abdominal aortic aneurysms results in redirection of blood through the deployed endograft (EG). Even though EVAR is clinically effective, the absolute flow restoration is not warranted. Our purpose was to compare the physiological with the post-EVAR infrarenal flow conditions. We developed patient-specific models based on computed tomography data of five healthy volunteers and ten patients treated with the Endurant^® stent-graft system. Wall shear stress (WSS), helicity, pressure and velocity fields were calculated using computational fluid dynamics. The results showed a decrease of peak WSS on the part of the EG that resides in the iliac arteries, compared to the physiological value (p = 0.01). At the abdominal part, the average helicity seems to increase after EVAR, while at the iliac arteries part, the intensity of helical flow seems physiological. Pressure drop and peak velocity in the iliac arteries part are lower than the physiological values (p = 0.04). The comparison revealed that most hemodynamic properties converge to normal levels at the abdominal part whereas statistically significant variations were observed in the iliac arteries part. The delineation of the differences between physiological and postoperative flow data could pave the way for the improvement of EG designs.     

Flow patterns in three-dimensional porcine epicardial coronary arterial tree

The branching pattern of epicardial coronary arteries is clearly three-dimensional, with correspondingly complex flow patterns. The objective of the present study was to perform a detailed hemodynamic analysis using a three-dimensional finite element method in a left anterior descending (LAD) epicardial arterial tree, including main trunk and primary branches, based on computed tomography scans. The inlet LAD flow velocity was measured in an anesthetized pig, and the outlet pressure boundary condition was estimated based on scaling laws. The spatial and temporal wall shear stress (WSS), gradient of WSS (WSSG), and oscillatory shear index (OSI) were calculated and used to identify regions of flow disturbances in the vicinity of primary bifurcations. We found that low WSS and high OSI coincide with disturbed flows (stagnated, secondary, and reversed flows) opposite to the flow divider and lateral to the junction orifice of the main trunk and primary branches. High time-averaged WSSG occurs in regions of bifurcations, with the flow divider having maximum values. Low WSS and high OSI were found to be related through a power law relationship. Furthermore, zones of low time-averaged WSS and high OSI amplified for larger diameter ratio and high inlet flow rate. Hence, different focal atherosclerotic-prone regions may be explained by different physical mechanism associated with certain critical levels of low WSS, high OSI, and high WSSG, which are strongly affected by the diameter ratio. The implications of the flow patterns for atherogenesis are enumerated.     

Numerical simulation of steady flow fields in a model of abdominal aorta with its peripheral branches

In the present study, a numerical calculation procedure based on a finite volume method was developed to simulate steady flow fields in a model of abdominal aorta with its peripheral branches. The study focused on the steady baseline flow fields and the wall shear stress (WSS) distribution as well as the localization of the reversed flow regions and results were compared to those obtained by other investigators. In the case of resting conditions, the existence of a region of reversed flow of about one to two diameters in size and next to the renal arteries and along the posterior wall as observed by other researchers was confirmed. However, under the exercise conditions this region could be wiped out. The flow reversal along the lateral walls proximal to the bifurcation persisted in both rest and exercise conditions. The WSS distribution and the wall shear stress gradient distribution were obtained. The lowest WSS occurred near the ostia of the renal arteries and the lateral walls of the iliac arteries. And the highest is always at the turn to the branch. The results were generally consistent with those obtained experimentally and numerically by other investigators. It was also shown that the steady flow might be used to depict the averaged behavior of pulsatile flow. The present computer code provides a platform for the future more realistic simulations.     

Haemodynamic impact of stent–vessel (mal)apposition following carotid artery stenting: mind the gaps!

Carotid artery stenting (CAS) has emerged as a minimally invasive alternative to endarterectomy but its use in clinical treatment is limited due to the post-stenting complications. Haemodynamic actors, related to blood flow in the stented vessel, have been suggested to play a role in the endothelium response to stenting, including adverse reactions such as in-stent restenosis and late thrombosis. Accessing the flow-related shear forces acting on the endothelium in vivo requires space and time resolutions which are currently not achievable with non-invasive clinical imaging techniques but can be obtained from image-based computational analysis. In this study, we present a framework for accurate determination of the wall shear stress (WSS) in a mildly stenosed carotid artery after the implantation of a stent, resembling the commercially available Acculink (Abbott Laboratories, Abbott Park, Illinois, USA). Starting from angiographic CT images of the vessel lumen and a micro-CT scan of the stent, a finite element analysis is carried out in order to deploy the stent in the vessel, reproducing CAS in silico. Then, based on the post-stenting anatomy, the vessel is perfused using a set of boundary conditions: total pressure is applied at the inlet, and impedances that are assumed to be insensitive to the presence of the stent are imposed at the outlets. Evaluation of the CAS outcome from a geometrical and haemodynamic perspective shows the presence of atheroprone regions (low time-average WSS, high relative residence time) colocalised with stent malapposition and stent strut interconnections. Stent struts remain unapposed in the ostium of the external carotid artery disturbing the flow and generating abnormal shear forces, which could trigger thromboembolic events.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

An approach to quantitative assessment of hemodynamic differences between unruptured and ruptured ophthalmic artery aneurysms

Background and purpose: Hemodynamic parameters are important in the pathogenesis, evolution and rupture of intracranial aneurysm. Energy loss (EL) has been applied for the rupture risk prediction of artery aneurysms recently. We proposed a new EL and further investigate its effects on the rupture of aneurysms. Materials and methods: Sixty-four patient-specific ophthalmic aneurysm datasets were divided into ruptured and unruptured groups based on their clinical history. Based on patient-specific 3D-DSA data, realistic models were retrospectively reconstructed and then analyzed by using computational fluid dynamic method. Results: The flow field feature EL in ruptured cases was significantly higher than that in unruptured cases. The average wall shear stress (WSS) and the maximum WSS in ruptured cases were higher than those in unruptured cases. Modified pressure loss coefficient (PLC_M) in ruptured cases was slight higher than that in unruptured cases but the difference has no statistical significance. Multivariate logistic regression analysis demonstrated flow field feature EL (p < 0.05) and the maximum WSS (p < 0.05) were the only independently significant variables to predict rupture of ophthalmic aneurysm. There were no differences in PLC_M, the maximum oscillatory shear index (OSI), the average OSI and AR between the two groups. Conclusion: Flow field feature EL may be a reliable factor to predict the rupture risk of aneurysms.     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

Modelling evolution and the evolving mechanical environment of saccular cerebral aneurysms

A fluid–solid-growth (FSG) model of saccular cerebral aneurysm evolution is developed. It utilises a realistic two-layered structural model of the internal carotid artery and explicitly accounts for the degradation of the elastinous constituents and growth and remodelling (G&R) of the collagen fabric. Aneurysm inception is prescribed: a localised degradation of elastin results in a perturbation in the arterial geometry; the collagen fabric adapts, and the artery achieves a new homeostatic configuration. The perturbation to the geometry creates an altered haemodynamic environment. Subsequent degradation of elastin is explicitly linked to low wall shear stress (WSS) in a confined region of the arterial domain. A sidewall saccular aneurysm develops, the collagen fabric adapts and the aneurysm stabilises in size. A quasi-static analysis is performed to determine the geometry at diastolic pressure. This enables the cyclic stretching of the tissue to be quantified, and we propose a novel index to quantify the degree of biaxial stretching of the tissue. Whilst growth is linked to low WSS from a steady (systolic) flow analysis, a pulsatile flow analysis is performed to compare steady and pulsatile flow parameters during evolution. This model illustrates the evolving mechanical environment for an idealised saccular cerebral aneurysm developing on a cylindrical parent artery and provides the guidance to more sophisticated FSG models of aneurysm evolution which link G&R to the local mechanical stimuli of vascular cells.     

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.